Create Densities — create_densities • raceland

Calculate local densities of subpopulations (race-specific local densities)

create_densities(x, y, window_size)

Arguments

x: SpatRaster with realizations
y: SpatRaster with shares of subpopulations
window_size: Size, expressed in the number of cells, of a square-shaped local window for which local densities will be calculated; it is recommended to use the small window_size, i.e., 10

Value

a SpatRaster containing n local densities. Local density layer is calculated for each realization

Examples

library(terra)
race_raster = rast(system.file("extdata/race_raster.tif", package = "raceland"))
real_rasters = create_realizations(race_raster, n = 5)
plot(real_rasters)

dens_raster = create_densities(real_rasters, race_raster, window_size = 10)
dens_raster
#> class       : SpatRaster 
#> size        : 60, 60, 5  (nrow, ncol, nlyr)
#> resolution  : 30, 30  (x, y)
#> extent      : 685545, 687345, 2103255, 2105055  (xmin, xmax, ymin, ymax)
#> coord. ref. : USA_Contiguous_Albers_Equal_Area_Conic 
#> source(s)   : memory
#> names       :        lyr.1,       lyr.1,       lyr.1,        lyr.1,       lyr.1 
#> min values  :  0.004819739, 0.004877809, 0.004763037,  0.004763037, 0.004991135 
#> max values  : 19.578132380, 8.912088643, 9.735491307, 14.666666985, 9.325799625 
plot(dens_raster)