Description Usage Arguments Details Value Author(s) Examples
Currently there are only two GIS functions: one for producing spatial concentration maps (GIS.Concentration.matrix
) and one for using spatial
population data and concentration maps to calculate exposure (GIS.Exposure
).
1 2 3 | GIS.Concentration.matrix(Emission, LO, LA, distx = 10.5, disty = 10.5,
resolution = 1, N = 1000, dbug = FALSE, ...)
GIS.Exposure(Concentration.matrix, dbug = FALSE, ...)
|
Emission |
|
Concentration.matrix |
|
LO |
|
LA |
|
distx |
|
disty |
|
resolution |
|
N |
|
dbug |
use |
... |
excess arguments are ignored or passed to |
The concentration matrix is computed using PILTTI source-receiver-matrices (http://en.opasnet.org/w/Piltti_source-receptor_matrix). They are originally for modeling PM2.5 distributions in a few Finnish cities between the years 2000 and 2003. To produce a rudimentary probability distribution these matrices are randomized between iterations.
Exposure is calculated by matching a concentration matrix to Finnish population data (http://en.opasnet.org/w/Special:Opasnet_Base?id=op_en2949.2012).
LA and LO are not required arguments for exposure, but speed the computation significantly.
See also: http://en.opasnet.org/
GIS.Concentration.matrix
returns an ovariable
whose output is a grid defined as bins for coordinates.
GIS.Exposure
returns an ovariable
whose output is concentration * population. All population data that matches cells defined by the
concentration matrix is included in the output.
T. Rintala teemu.rintala.a@gmail.com
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 | ## Not run:
# Excerpt from http://en.opasnet.org/w/Health_impacts_of_fine_particles_in_Rauma
# (not evaluated)
# Paasto Emissions
Paasto <- new(
"ovariable",
name = "Paasto",
dependencies = data.frame(Name = "tieliikennepaastot", Key = "0194s0uuucjxq8Wi"),
formula = function(dependencies, ...) {
ComputeDependencies(dependencies, ...)
# Muutetaan paivapaasto vuosipaastoksi ja grammat tonneiksi
out <- tieliikennepaastot * 365 * 1E-6
return(out)
}
)
# Muita tarpeellisia arvoja Other important values
bg.mort <- 45182 / 5203826 # same values as used in PILTTI
## J. T. Tuomisto, A. Wilson, et al. Uncertainty in mortality response to
## airborne fine particulate matter... 2008
erf <- 0.0097
# unit: m^3 /ug
# Ovariablet
## Pitoisuudet Concentrations
Pitoisuus <- new(
"ovariable",
name = "Pitoisuus",
dependencies = data.frame(
Name = c("Paasto", "LO", "LA")
),
formula = function(dependencies, ...) {
ComputeDependencies(dependencies, ...)
temp <- GIS.Concentration.matrix(Paasto, LO, LA, ...)
return(temp)
}
)
## Altistuminen Exposure
Altistuminen <- new(
"ovariable",
name = "Altistuminen",
dependencies = data.frame(
Name = c("Pitoisuus", "LO", "LA")
),
formula = function(dependencies, ...) {
ComputeDependencies(dependencies, ...)
out <- GIS.Exposure(Pitoisuus, LO, LA, ...)
return(out)
}
)
## End(Not run)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.