## Model spectra using PROSAIL
parameter <- data.frame(N = rep.int(c(1, 1.5),2), LAI = c(1,1,3,3))
spec <- PROSAIL(parameterList=parameter)
## Transform spectra
spec_clman <- transformSpeclib(spec, method = "sh", out = "raw")
## Fix points
spec_cp <- getcp(spec_clman, 1, c(400, 800))
spec_cp
Run the code above in your browser using DataLab