近红外光谱建模之光谱预处理python实现（二）

最新推荐文章于 2025-07-02 21:19:30 发布

原创

最新推荐文章于 2025-07-02 21:19:30 发布 · 1.3w 阅读

253 ·

CC 4.0 BY-SA版权

文章标签：

#近红外光谱 #预处理方法

4 多元散射矫正处理

def msc(sdata):
    n = sdata.shape[0]  # 样本数量
    k = np.zeros(sdata.shape[0])
    b = np.zeros(sdata.shape[0])

    M = np.mean(sdata, axis=0)

    from sklearn.linear_model import LinearRegression
    for i in range(n):
        y = sdata[i, :]
        y = y.reshape(-1, 1)
        M = M.reshape(-1, 1)
        model = LinearRegression()
        model.fit(M, y)
        k[i] = model.coef_
        b[i] = model.intercept_

    spec_msc = np.zeros_like(sdata)
    for i in range(n):
        bb = np.repeat(b[i], sdata.shape[1])
        kk = np.repeat(k[i], sdata.shape[1])
        temp = (sdata[i, :] - bb)/kk
        spec_msc[i, :] = temp
    return spec_msc


x = msc(x)
pp = PlotSpectrum(x)
pp.show()