画光谱图——R语言graphics包matplot等函数

由于前几天一直找如何将光谱图画出，为二分类光谱数据，找了半天，，，，
最后在循环的条件下，一条条画完，，，，，，
最近在prospectr包中忽然发现matplot等函数

数据格式

原始数据处理

##设置读取数据位置
setwd("C:/Users/chengdehe/Desktop/processing")
##读取数据
data <- read.csv("H3.csv",check.names=F)
##去除表头，只剩下光谱数据
data_ref <- data[ ,-(1:2)]

光谱预处理

#####S_G处理进行光谱平滑#####
SG <- savitzkyGolay(data_ref, p = 3, w = 11, m = 0)

处理前后光谱对比——即将两种光谱图画到统一坐标系中

###画原始数据
matplot(as.numeric(colnames(data_ref)), t(data_ref[ , ]),
        type = 'l', xlab = 'Wavelength', ylab = 'Absorbance',col="red")
###画处理后数据
matlines(as.numeric(colnames(SG)),t(SG[ , ]),
              lty = 1,,col="blue")
 ###as.numeric(colnames(data_ref))将光谱列表头，即光谱波长转换为数值，作为X轴

同样也可以将转换数据换为，二分类结果，只是需要提前将数据集分开，分别赋值
原始光谱图

叠加光谱图

从图可见转换前后并无太大变化