# \donttest{
library(terra)
r <- rast(nrows=500, ncols=500, xmin=571823, xmax=616763,
ymin=4423540, ymax=4453690)
crs(r) <- "epsg:9001"
r[] <- runif(ncell(r), 1000, 2500)
# Calculate 10th percentile for 3x3 window
r.p10 <- raster.moments(r, type="quantile", p=0.10)
# }
Run the code above in your browser using DataLab