Distribution by OpenStreetnMap street

Distribute emissions by streets of OpenStreetMap

streetDist(
  emission = 1,
  dist = c(1, 0, 0, 0, 0),
  grid = NULL,
  osm = NULL,
  epsg = 31983,
  warnings = FALSE
)

Arguments

emission	Numeric; emissions.
dist	Numeric; vector with length 5. The order represents motorway, trunk, primary, secondary and tertiary
grid	'sf' POLYGON; grid of polygons class sf.
osm	streets of OpenStreetMaps class sf
epsg	Numeric; spatial code for projecting spatial data
warnings	Logical; to show warnings.

Value

grid of polygon

Examples

if (FALSE) {
# Do not run
library(sf)
# Download OSM streets
streets <- st_read("path")
streets <- streets[streets$highway != "residential", ]
# Grid
grid  <- gridInfo(paste(system.file("extdata", package = "EmissV"),"/wrfinput_d02",sep=""))
names(grid)
d3 <- data.frame(x = as.numeric(grid$Lon),
                 y = as.numeric(grid$Lat))
d3 <- st_as_sf(d3, coords = c("x","y"))
st_crs(d3) <- st_crs(4326)
library(vein)
g <- st_transform(st_as_sf(vein::make_grid(as(st_transform(d3, 31983),
                  "Spatial"),
               grid$DX*1000, grid$DX*1000, T)), 4326)
streets$id <- NULL
per <- c(1, 0, 0, 0, 0)
teste <- streetDist(emission = 1000000, dist = per, grid = g,
                    osm = streets, epsg = 31983)
# Another example:
library (EmissV)
library (osmdata)
library (sf)
city <- "accra"
bb <- getbb (city)
dat <- opq (bbox = city) %>%
  add_osm_feature (key = "highway") %>%
  osmdata_sf (quiet = FALSE) %>%
  osmdata::osm_poly2line () %>%
  magrittr::extract2 ("osm_lines")
#saveRDS (dat, file = "accra-hw.Rds")
utm <- 32630 # for Accra

# Get a raster grid of population density to use for the emission distribution:
url <- paste0 ("https://github.com/ATFutures/who-data/releases/download/",
               "v0.0.2-worldpop-tif-gha-npl/accra.2fpopdens.2fGHA15adj_040213.tif")
download.file (url, "accra-pop.tif", mode = "wb")
ras <- raster::raster ("accra-pop.tif") %>%
  raster::crop (raster::extent (bb)) %>%
  as ("SpatialPolygons") %>%
  st_as_sf ()

#dat <- readRDS (file = "accra-hw.Rds")
dat <- dat[dat$highway %in% c ("motorway", "trunk", "primary",
                                "secondary", "teritary"), ]

s <- streetDist (emission = 1, dist = c (1, 0, 0, 0, 0), grid = ras,
                 osm = dat, epsg = utm)
}