SNV的python代码

import numpy as np
import pandas as pd
import os
import matplotlib.pyplot as plt
a=os.getcwd()
print(a)
data=pd.read_excel("Q1.xlsx")   # x条光谱，y个光谱通道
data=np.array(data)
print(data)
m=data.shape[0]
n=data.shape[1]
print(m,n)#
#求标准差
data_std=np.std(data,axis=1)   # 每条光谱的标准差
print(data_std)
#求平均值
data_average=np.mean(data,axis=1)   # 每条光谱的平均值
print(data_average)
#SNV计算
data_snv=[[((data[i][j]-data_average[i])/data_std[i]) for j in range(n)] for i in range(m)]
plt.subplot(1,2,1)
plt.plot(data)
plt.subplot(1,2,2)
plt.plot(data_snv)
plt.show()