fleet <- vehicles(example = TRUE)
EmissionFactors <- emissionFactor(example = TRUE)
TOTAL <- totalEmission(fleet,EmissionFactors,pol = c("CO"),verbose = TRUE)
grid <- gridInfo(paste0(system.file("extdata", package = "EmissV"),"/wrfinput_d01"))
shape <- raster::shapefile(paste0(system.file("extdata", package = "EmissV"),"/BR.shp"))
raster <- raster::raster(paste0(system.file("extdata", package = "EmissV"),"/dmsp.tiff"))
SP <- areaSource(shape[22,1],raster,grid,name = "SP")
RJ <- areaSource(shape[17,1],raster,grid,name = "RJ")
e_CO <- emission(total = TOTAL,
pol = "CO",
area = list(SP = SP, RJ = RJ),
grid = grid,
mm = 28)
Run the code above in your browser using DataLab