高（近外）光谱多元散射校正（MSC）

Posted 2023-02-26 志行千里90

尽管CSDN中能扒拉到的资源不少，但实际能用的不多，自己写的一个希望能帮助到各位。另外，在第二步求“光谱回归”系数的地方一直存疑。是直接用一元回归求得回归系数，还是用最小二乘拟合回归系数，了解的朋友留言。

输入数据示例图：

#!/usr/bin/env python3.7
# encoding: utf-8
"""
  @author: ISR
  @contact: [e-mail]
  @file: 
  @time: 
  @desc:
  对采集的光谱数据进行多元散射校正（MSC）
  输入格式为csv,数据格式首行重复测量名称，首列为波段名称
"""
import numpy as np
from pandas.core.frame import DataFrame
from sklearn.linear_model import LinearRegression
import pandas as pd

def msc(data):
    """"
    data 横轴方向表示重复测量的数据，纵轴方向表示波段数量
    """
    # 计算平均光谱，实际就是x值
    s_mean = np.mean(data, axis=1)

    # 行列数
    r, c = data.shape

    # 创建一个单位矩阵
    msc_x = np.ones((r, c))

    # 遍历各列，实际是各重复测量
    for i in range(c):

        # y值
        y = data[:, i]

        # 计算光谱回归系数Ki,Bi
        lin = LinearRegression()
        lin.fit(s_mean.reshape(-1, 1), y.reshape(-1, 1))
        k = lin.coef_
        b = lin.intercept_

        msc_x[:, i] = (y - b) / k

    return msc_x

if __name__ == '__main__':

    txt_path = r"C:\\Users\\79420\\Desktop\\test\\30比1\\30比1.csv"
    df = pd.read_csv(txt_path, header=0, index_col=0, engine="python", encoding="utf-8")

    # 行列名称
    columns_n = df.columns.values.tolist()
    row_n = df._stat_axis.values.tolist()
    print("0行，1列".format(len(row_n),len(columns_n)))

    df_arr = df.values

    msc_arr = msc(df_arr)

    txt_path2 = r"C:\\Users\\79420\\Desktop\\test\\30比1\\30比1—b.csv"
    dabai1_arr_data_you = DataFrame(msc_arr, index=row_n, columns=columns_n)
    dabai1_arr_data_you.to_csv(txt_path2)