ShaleVolume.py 文件源码

def GRsd (GRlog):
    em = GMM(n_components=3)
    em.fit(GRlog.reshape((-1, 1)))
    idxminmeangr = np.argmin(em.means_)
    grsd = em.means_[idxminmeangr] - em.covars_[idxminmeangr]**0.5

    return grsd