Description Usage Arguments Value References Examples
View source: R/decom.weibull.R
The function allows you to fit waveforms with the Weibull function and get parameters estimated.
1 | decom.weibull(x, smooth = TRUE, thres = 0.22, width = 3)
|
x |
is a waveform with a index at the begining and followed with intensities. |
smooth |
is tell whether you want to smooth the waveform to remove some obvious outliers. Default is TRUE. |
thres |
is to determine if the detected peak is the real peak whose intensity should be higher than threshold*maximum intensity. Default is 0.22. |
width |
width of moving window.Default is 3, must be odd integer between 1 and n.This parameter ONLY work when the smooth is TRUE. |
A list contains estimates of A (amplitude parameter), u (location parameter), sigma (scale parameter) and k (shape parameter) after decomposition. these parameters are different from the adaptive Gaussian and Gaussian decomposition.
Tan Zhou, and Sorin C. Popescu, 2017. Bayesian decomposition of full waveform LiDAR data with uncertainty analysis. Remote Sensing of Environment 200 (2017): 43-62.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 | ##import return waveform data
data(return)
lr<-nrow(return)
ind<-c(1:lr)
return<-data.frame(ind,return)
x<-return[1,] ###must be a dataset including intensity with index at the beginning.
r1<-decom.weibull(x)
r2<-decom.weibull(return[2,])
# for the whole dataset
dr3<-apply(return,1,decom.weibull)
dd<-return[10:20,]
dr4<-apply(dd,1,decom.weibull)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.