data(spectral_data)
## Transform speclib
bd <- transformSpeclib(spectral_data, method = "sh", out = "bd")
## Define features automatically
features <- define.features(bd)
##Example to isolate the features around 450nm, 700nm, 1200nm and 1500nm.
featureSelection <- specfeat(features, c(450,700,1200,1500))
## Plot features
plot(featureSelection, 1:4)
## Advanced plotting example
plot(featureSelection, 1:4, stylebysubset = "season")
plot(featureSelection, 1:4, stylebysubset = "season", changecol = FALSE,
changetype = TRUE)
Run the code above in your browser using DataLab