2019-12-20

1.利用R语言的plsr包，对全波段光谱数据进行最小二乘法

install.packages("plsr")
library(pls)
data<-read.csv("C:/Users/ZWPC/Desktop/pls.csv",header=F)
pls <- plsr(V1~.,data=data,30,validation="CV")#全波段最小二乘法，VI为因变量，data为剩余的全波段光谱反射率，30可根据自己需求自定义，我这里选择的是30个进行一次CV(交叉验证）
selectNcomp(pls,plot=T)#可查看从30里选出来的RMSEP的变化状态
cor(data$V1,pls$fitted.values[,,10])^2#查看R2值
summary(pls)
plot(data$V1,pls$fitted.values[,,10])#可查看因变量与个波段的一个相关性

2.利用matlab对波段的一个选择

SPA法即连续投影法