Calculates the chi maps for one rasterStack or all raster all the raster pairs stored in two rasterStacks. As observed values, the first stack is used. The expected value is either set to the mean of the first stack, or if given to be the values of the second stack.

chiMaps(stack1, stack2 = NULL, verbose = FALSE)

Arguments

stack1

rasterStack

stack2

rasterStack NULL or containing the same number of rasterLayers and has euqal extent and resolution.

verbose

logical: print currently processed height band in raster stack?

Value

A rasterStack containing the chi maps.

Examples

chiMaps(raster::stack(dem))
#> class      : RasterStack 
#> dimensions : 175, 501, 87675, 1  (nrow, ncol, ncell, nlayers)
#> resolution : 251.2, 370.8  (x, y)
#> extent     : 2543252, 2669103, 1177108, 1241998  (xmin, xmax, ymin, ymax)
#> crs        : +init=epsg:2056 +proj=somerc +lat_0=46.95240555555556 +lon_0=7.439583333333333 +k_0=1 +x_0=2600000 +y_0=1200000 +ellps=bessel +towgs84=674.374,15.056,405.346,0,0,0,0 +units=m +no_defs 
#> names      : chiMap.srtm_38_03 
#> min values :         -14.30596 
#> max values :          65.88701 
#>