spatialEco (version 1.3-7)

raster.moments: Raster moments

Description

Calculates focal statistical moments of a raster

Usage

raster.moments(x, type = "mean", s = 3, p = 0.75)

Value

raster class object of the local distributional moment

Arguments

x

raster object

type

The global statistic to represent the local deviation options are: "min", "min", "mean", "median", "var, "sd", "mad", "kurt", "skew", "quantile"

s

Size of matrix (focal window), can be single value or two values defining the [x,y] dimensions of the focal matrix

p

if type="quantile", the returned percentile.

Author

Jeffrey S. Evans <[e-mail]>

Examples

# \donttest{
  library(raster)
  r <- raster(nrows=100, ncols=100, xmn=571823, xmx=616763, 
              ymn=4423540, ymx=4453690)
  proj4string(r) <- crs("+proj=utm +zone=12 +datum=NAD83 +units=m +no_defs")
  r[] <- runif(ncell(r), 1000, 2500)
  r <- focal(r, focalWeight(r, 150, "Gauss") )

# Calculate 10th percentile for 3x3 window
r.p10 <- raster.moments(r, type="quantile", p=0.10) 
# }