precf <- system.file("map/Mekrou_precipitation.grd", package="rasterList")
prec <- stack(precf)
# \donttest{
## Sample L-moments
if (requireNamespace("lmom",quietly = TRUE)) {
library(lmom)
samlmom <- stack(rasterList(prec,FUN=samlmu))
## Fitting a Random Probability Distribution: it is a 'rasterList' Object
fitdist <- rasterList(samlmom,FUN=pelgam)
##### ZOOM IN
## set a mask
mask <-raster( extent(fitdist)/4 )
fitdist_masked <- crop ( x = fitdist,y=mask)
}
# }
Run the code above in your browser using DataLab