EmissV (version 0.664.5)

read: Read NetCDF data from global inventaries

Description

Read data from global inventoris

Usage

read(file, version = "EDGAR 4.3.1 v2", as_raster = T, verbose = T)

Arguments

file

file name or names (variables are summed)

version

inventorie information

as_raster

return a raster (defoult) or matrix (with units)

verbose

display additional information

Value

Matrix or raster

See Also

rasterSource and gridInfo

Examples

# NOT RUN {
d1     <- gridInfo(paste(system.file("extdata", package = "EmissV"),"/wrfinput_d01",sep=""))
d2     <- gridInfo(paste(system.file("extdata", package = "EmissV"),"/wrfinput_d02",sep=""))
print("download and untar EDGAR data from:")
print("http://edgar.jrc.ec.europa.eu/gallery.php?release=v431_v2&substance=NOx§or=TRO")
nox    <- read("v431_v2_REFERENCE_NOx_2010_10_TRO.0.1x0.1.nc")
sp::spplot(nox, scales = list(draw=TRUE), xlab="Lat", ylab="Lon",main="NOx emissions from EDGAR")
nox_d1 <- rasterSource(nox,d1)
nox_d2 <- rasterSource(nox,d2)
image(nox_d1, main = "NOx emissions from transport from EDGAR 3.4.1 for d1")
image(nox_d2, main = "NOx emissions from transport from EDGAR 3.4.1 for d2")
# }