areaSource: Distribution of emissions by area

In Schuch666/EmissV: Tools for Create Emissions for Air Quality Models

Distribution of emissions by area

Description

Calculate the spatial distribution by a raster masked by shape/model grid information.

Usage

areaSource(s, r, grid = NA, name = "", as_frac = FALSE, verbose = TRUE)

Arguments

s	input shape object
r	input raster object
grid	grid with the output format
name	area name
as_frac	return a fraction instead of the raster value
verbose	display additional data

Details

About the DMSP and example data https://en.wikipedia.org/wiki/Defense_Meteorological_Satellite_Program

Value

a raster object containing the spatial distribution of emissions

Source

Data avaliable https://www.nesdis.noaa.gov/current-satellite-missions/currently-flying/defense-meteorological-satellite-program

Examples

shape  <- raster::shapefile(paste(system.file("extdata", package = "EmissV"),
                            "/BR.shp",sep=""),verbose = FALSE)
shape  <- shape[22,1] # subset for Sao Paulo - BR
raster <- raster::raster(paste(system.file("extdata", package = "EmissV"),
                         "/dmsp.tiff",sep=""))
grid   <- gridInfo(paste(system.file("extdata", package = "EmissV"),"/wrfinput_d02",sep=""))
SP     <- areaSource(shape,raster,grid,name = "SPMA")

raster::plot(SP,ylab="Lat",xlab="Lon",
             main ="Spatial Distribution by Lights for Sao Paulo - Brazil")

Schuch666/EmissV documentation built on Dec. 24, 2024, 3:15 a.m.