Read a georeference image and extract values for using in griding

extract values from a image to be used on each region

possess(geoemiss, filename = NA, plots = T, verbose = T)

Arguments

geoemiss	a sf from geoemiss
filename	image filename
plots	plot individual regions
verbose	display additional information

Value

a sf with the area shapes and the total emission and a raster image for each area

Examples

states <- sf::st_read(paste(system.file("extdata", package = "inventory"),"/states.shp",sep=""),
                      quiet = TRUE)
Nox    <- geoemiss(geom = states,variable = "Nox",names = c("sp","rj"),values = c(1000,20))
#> calculating Nox for 2 areas
image  <- paste(system.file("extdata", package = "inventory"),"/tiny.tif",sep="")
ras    <- possess(Nox,image,plots = FALSE)
test   <- griding(ras,variable = "Nox", res = 0.1, type = "local", plot = TRUE)
#> processing sp area ...
#> dist is assumed to be in decimal degrees (arc_degrees).
#> although coordinates are longitude/latitude, st_intersection assumes that they are planar
#> processing rj area ...
#> dist is assumed to be in decimal degrees (arc_degrees).
#> although coordinates are longitude/latitude, st_intersection assumes that they are planar