|
|
@ -1,11 +1,13 @@
|
|
|
|
## ----setup, include=FALSE------------------------------------------------
|
|
|
|
## ----setup, include=FALSE------------------------------------------------
|
|
|
|
knitr::opts_chunk$set(echo = TRUE)
|
|
|
|
knitr::opts_chunk$set(echo = TRUE)
|
|
|
|
|
|
|
|
|
|
|
|
## ------------------------------------------------------------------------
|
|
|
|
## ----results = F---------------------------------------------------------
|
|
|
|
library(magrittr)
|
|
|
|
library(magrittr)
|
|
|
|
library(opensensmapr)
|
|
|
|
library(opensensmapr)
|
|
|
|
|
|
|
|
|
|
|
|
all_sensors = osem_boxes()
|
|
|
|
all_sensors = osem_boxes()
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
## ------------------------------------------------------------------------
|
|
|
|
summary(all_sensors)
|
|
|
|
summary(all_sensors)
|
|
|
|
|
|
|
|
|
|
|
|
## ----message=F, warning=F------------------------------------------------
|
|
|
|
## ----message=F, warning=F------------------------------------------------
|
|
|
@ -22,13 +24,14 @@ str(phenoms)
|
|
|
|
## ------------------------------------------------------------------------
|
|
|
|
## ------------------------------------------------------------------------
|
|
|
|
phenoms[phenoms > 20]
|
|
|
|
phenoms[phenoms > 20]
|
|
|
|
|
|
|
|
|
|
|
|
## ------------------------------------------------------------------------
|
|
|
|
## ----results = F---------------------------------------------------------
|
|
|
|
pm25_sensors = osem_boxes(
|
|
|
|
pm25_sensors = osem_boxes(
|
|
|
|
exposure = 'outdoor',
|
|
|
|
exposure = 'outdoor',
|
|
|
|
date = Sys.time(), # ±4 hours
|
|
|
|
date = Sys.time(), # ±4 hours
|
|
|
|
phenomenon = 'PM2.5'
|
|
|
|
phenomenon = 'PM2.5'
|
|
|
|
)
|
|
|
|
)
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
## ------------------------------------------------------------------------
|
|
|
|
summary(pm25_sensors)
|
|
|
|
summary(pm25_sensors)
|
|
|
|
plot(pm25_sensors)
|
|
|
|
plot(pm25_sensors)
|
|
|
|
|
|
|
|
|
|
|
@ -45,6 +48,7 @@ berlin = st_point(c(13.4034, 52.5120)) %>%
|
|
|
|
st_transform(4326) %>% # the opensensemap expects WGS 84
|
|
|
|
st_transform(4326) %>% # the opensensemap expects WGS 84
|
|
|
|
st_bbox()
|
|
|
|
st_bbox()
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
## ----results = F---------------------------------------------------------
|
|
|
|
pm25 = osem_measurements(
|
|
|
|
pm25 = osem_measurements(
|
|
|
|
berlin,
|
|
|
|
berlin,
|
|
|
|
phenomenon = 'PM2.5',
|
|
|
|
phenomenon = 'PM2.5',
|
|
|
|