光谱反射率及四种变换相关系数求取

import numpy as np
import pandas as pd


def msc(X):
    # +++++ 输入：X = m × p 矩阵，m个样本,p个特征。X 应为 ndarray 数据类型）
    # +++++ 输出：X_msc = m × p
    me = np.mean(X, axis=0)
    [m, p] = np.shape(X)
    X_msc = np.zeros((m, p))

    for i in range(m):
        poly = np.polyfit(me, X[i], 1)  # 每个样本做一次一元线性回归
        j = 0
        for j in range(p):
            X_msc[i, j] = (X[i, j] - poly[1]) / poly[0]

    return X_msc

if __name__ == '__main__':

    data1 = pd.read_excel('./sampleSOM.xlsx',header=None).values
    data2