ShaleVolume.py 文件源码

def GRsh (GRlog):
    em = GMM(n_components=3)
    em.fit(GRlog.reshape((-1, 1)))
    idxmaxmeangr = np.argmax(em.means_)
    grsh = em.means_[idxmaxmeangr] + em.covars_[idxmaxmeangr]**0.5

    return grsh