In SevillaR/aire: Grab data on air pollution in Andalusia (Spain)

library(knitr)
opts_knit$set(root.dir=normalizePath('../'))
opts_chunk$set(warning = FALSE, message = FALSE)
opts_chunk$set(fig.path = "nubepolvo_figuras/", fig.show="hide")

Entre los días 20 y 22 de Febrero de 2016 ha habido una nube de polvo sahariano sobre la Península, que se ha dejado notar también en Sevilla (http://www.abc.es/sociedad/abci-medio-ambiente-nube-polvo-sahariano-dispara-contaminacion-aire-particulas-toda-peninsula-201602231715_noticia.html).

Vamos a utilizar nuestro paquete aire para analizar la calidad del aire (más concretamente la concentración de partículas) en la última semana.

Primero cargamos el paquete:

library(aire)

Ahora descargamos los datos de la última semana para la provincia de Sevilla:

datos <- get_daily_data(province = "se", from = "2016-02-15", to = "2016-02-23")

Y los analizamos con el paquete openair:

library(openair)

Con la función summaryPlot podemos inspeccionar los datos:

summaryPlot(datos, pollutant = "PART", period = "months")

La función timePlot nos dibuja la serie temporal para cada estación más claramente:

timePlot(datos, pollutant = "PART", type = "site")

Vemos que el pico de contaminación por partículas se produjo entre el 20 y 22 de Febrero, y que se alcanzaron valores de hasta 500 microgramos por metro cúbico de aire, cuando los valores normales están por debajo de 25 microgramos (http://bit.ly/airedata).

También podemos hacer un calendarPlot:

calendarPlot(datos, pollutant = "PART")

Finalmente hacemos un mapa de la concentración de partículas por estaciones para el día de mayor contaminación (21 de Febrero).

Primero descargamos y preparamos los datos:

datos.21 <- get_daily_data("2016-02-21", "se")
data(stations)
datos.21.sp <- merge(datos.21, stations[, "site"], by = "site", all.x = TRUE)

Y luego hacemos el mapa:

GoogleMapsPlot(datos.21.sp, pollutant = "PART")  

Se puede observar que aunque la densidad de partículas fue muy alta en toda la zona, hubo estaciones (como la de la Sierra Norte) con valores algo mejores.

SevillaR/aire documentation built on May 9, 2019, 1:23 p.m.