ShaleVolume.py 文件源码

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

    return spsh