Calculate the spatial distribution by a raster
rasterSource(r, grid, nlevels = "all", conservative = TRUE, verbose = TRUE)
Exemple data is a low resolution cutting from image of persistent lights of the Defense Meteorological Satellite Program (DMSP) https://pt.wikipedia.org/wiki/Defense_Meteorological_Satellite_Program
Data avaliable https://www.ospo.noaa.gov/Operations/DMSP/index.html
input raster object
grid object with the grid information
number of vertical levels off the emission array
TRUE (default) to conserve total mass, FALSE to conserve flux
display additional information
Returns a matrix
About the DMSP and example data https://en.wikipedia.org/wiki/Defense_Meteorological_Satellite_Program
gridInfo
and lineSource
grid <- gridInfo(paste(system.file("extdata", package = "EmissV"),"/wrfinput_d01",sep=""))
#> Grid information from: C:/Users/schuch/Documents/R/win-library/4.1/EmissV/extdata/wrfinput_d01
x <- raster::raster(paste(system.file("extdata", package = "EmissV"),"/dmsp.tiff",sep=""))
test <- rasterSource(x, grid)
#> Grid output: 99 columns 93 rows
image(test, axe = FALSE, main = "Spatial distribution by Persistent Nocturnal Lights from DMSP")