GriddedEmissionsArray returns a tranformed object with class "EmissionsArray" with 4 dimensios.

GriddedEmissionsArray(x, ..., cols, rows, times = ncol(x),
  rotate = FALSE)

# S3 method for GriddedEmissionsArray
print(x, ...)

# S3 method for GriddedEmissionsArray
summary(object, ...)

# S3 method for GriddedEmissionsArray
plot(x, ..., times = 1)

Arguments

x

Object with class "SpatialPolygonDataFrame", "sf" "data.frame" or "matrix"

...

ignored

cols

Number of columns

rows

Number of rows

times

Number of times

rotate

Logical to rotate TRUE or not FALSE the array

object

object with class "EmissionsArray'

Value

Objects of class "GriddedEmissionsArray"

Examples

if (FALSE) { data(net) data(pc_profile) data(fe2015) data(fkm) PC_G <- c(33491,22340,24818,31808,46458,28574,24856,28972,37818,49050,87923, 133833,138441,142682,171029,151048,115228,98664,126444,101027, 84771,55864,36306,21079,20138,17439, 7854,2215,656,1262,476,512, 1181, 4991, 3711, 5653, 7039, 5839, 4257,3824, 3068) veh <- data.frame(PC_G = PC_G) pc1 <- my_age(x = net$ldv, y = PC_G, name = "PC") pcw <- temp_fact(net$ldv+net$hdv, pc_profile) speed <- netspeed(pcw, net$ps, net$ffs, net$capacity, net$lkm, alpha = 1) pckm <- units::set_units(fkm[[1]](1:24), "km") pckma <- cumsum(pckm) cod1 <- emis_det(po = "CO", cc = 1000, eu = "III", km = pckma[1:11]) cod2 <- emis_det(po = "CO", cc = 1000, eu = "I", km = pckma[12:24]) #vehicles newer than pre-euro co1 <- fe2015[fe2015$Pollutant=="CO", ] #24 obs!!! cod <- c(co1$PC_G[1:24]*c(cod1,cod2),co1$PC_G[25:nrow(co1)]) lef <- ef_ldv_scaled(co1, cod, v = "PC", t = "4S", cc = "<=1400", f = "G",p = "CO", eu=co1$Euro_LDV) E_CO <- emis(veh = pc1,lkm = net$lkm, ef = lef, speed = speed, agemax = 41, profile = pc_profile, simplify = TRUE) class(E_CO) E_CO_STREETS <- emis_post(arra = E_CO, pollutant = "CO", by = "streets", net = net) g <- make_grid(net, 1/102.47/2, 1/102.47/2) #500m in degrees net@data <- net@data[,- c(1:9)] names(net) E_CO_g <- emis_grid(spobj = net, g = g, sr= 31983) head(E_CO_g) #class sf library(mapview) mapview(E_CO_g, zcol= "V1", legend = T, col.regions = cptcity::cptcity(1)) gr <- GriddedEmissionsArray(E_CO_g, rows = 19, cols = 23, times = 168, T) plot(gr) # For some cptcity color gradients: devtools::install_github("ibarraespinosa/cptcity") plot(gr, col = cptcity::cptcity(1)) }