我正在尝试使用以下公式计算亚氯酸盐光谱成熟度(CSM)和伊利石光谱成熟度(ISM)标量。到目前为止,我已经能够使用下面的R脚本读取单个.asd数据,但我正在努力编写批处理脚本来一次读取所有数据,我不知道如何实现我的公式。拜托,我需要帮助。谢谢。
式:
CSM Scalar
= (1 − ((2310:2370)) ) + 1.7 ∗ (1 − ((2240:2270)))/(((1810:2140)) − ((1900:2020)))
ISM Scalar
= ((2136:2296) − ((2140:2270)))/(((1810:2140)) − ((1900:2020)))
R SCRIPT BELOW:
data <- asdreader::get_spectra("data.ASD")
data
asdreader::get_metadata("data.ASD")
print(data)
matplot(t(data), type = 'l')
table = as.data.frame.table(data)