Calidad del aire cuarentena
Nano Castro
17 de mayo de 2020
Leer datos de los equipos de referencia
##
## Attaching package: 'lubridate'## The following object is masked from 'package:base':
##
## date## -- Attaching packages ---------------------------------------- tidyverse 1.2.1 --## v ggplot2 3.0.0 v purrr 0.2.5
## v tibble 1.4.2 v dplyr 0.7.6
## v readr 1.1.1 v stringr 1.3.1
## v ggplot2 3.0.0 v forcats 0.3.0## -- Conflicts ------------------------------------------- tidyverse_conflicts() --
## x lubridate::as.difftime() masks base::as.difftime()
## x lubridate::date() masks base::date()
## x dplyr::filter() masks stats::filter()
## x lubridate::intersect() masks base::intersect()
## x dplyr::lag() masks stats::lag()
## x lubridate::setdiff() masks base::setdiff()
## x lubridate::union() masks base::union()Lectura de datos CO - Thermo 48i
Lectura de datos NOx - Thermo 42i
## Lectura de datos O3 - Thermo 49i
Lectura datos PM - Thermo
Lectura de datos estación meteorologica
## 'data.frame': 1396 obs. of 39 variables:
## $ Date : Factor w/ 31 levels "1/5/2020","10/4/2020",..: 24 24 24 24 24 24 24 24 24 24 ...
## $ Time : Factor w/ 48 levels "00:00:00","00:30:00",..: 39 40 41 42 43 44 45 46 47 48 ...
## $ Temp : Factor w/ 204 levels "10.1","10.2",..: 77 73 70 66 62 57 52 47 43 41 ...
## $ Hi.Temp : Factor w/ 207 levels "10","10.1","10.3",..: 78 77 72 69 64 61 56 51 47 43 ...
## $ Low.Temp : Factor w/ 203 levels "10","10.1","10.2",..: 78 73 70 65 62 57 52 47 44 40 ...
## $ Hum : int 34 34 34 34 35 39 42 44 46 47 ...
## $ DewPt : Factor w/ 155 levels "-0.1","-0.2",..: 28 25 22 18 19 30 82 83 85 86 ...
## $ WindSpeed: Factor w/ 8 levels "0","0.4","0.9",..: 1 1 1 1 1 1 1 1 1 1 ...
## $ WindDir : Factor w/ 17 levels "---","E","ENE",..: 2 2 3 2 3 1 1 1 1 1 ...
## $ WindRun : Factor w/ 8 levels "0","0.8","1.61",..: 1 1 1 1 1 1 1 1 1 1 ...
## $ HiSpeed : Factor w/ 21 levels "0","0.4","0.9",..: 5 6 7 5 4 1 1 1 1 1 ...
## $ HiDir : Factor w/ 17 levels "---","E","ENE",..: 3 3 3 2 3 1 1 1 1 1 ...
## $ WindChill: Factor w/ 204 levels "10.1","10.2",..: 77 73 70 66 62 57 52 47 43 41 ...
## $ HeatInd : Factor w/ 198 levels "10","10.1","10.2",..: 59 56 52 48 45 43 39 35 32 30 ...
## $ THW : Factor w/ 198 levels "10","10.1","10.2",..: 59 56 52 48 45 43 39 35 32 30 ...
## $ THSW : Factor w/ 1 level "---": 1 1 1 1 1 1 1 1 1 1 ...
## $ Press : Factor w/ 353 levels "1000","1000.1",..: 326 329 331 333 335 337 339 342 344 343 ...
## $ Rain : Factor w/ 7 levels "0","0.25","0.51",..: 1 1 1 1 1 1 1 1 1 1 ...
## $ Rate : Factor w/ 6 levels "0","1.3","11.7",..: 1 1 1 1 1 1 1 1 1 1 ...
## $ Rad. : Factor w/ 1 level "---": 1 1 1 1 1 1 1 1 1 1 ...
## $ Energy : Factor w/ 1 level "---": 1 1 1 1 1 1 1 1 1 1 ...
## $ Rad..1 : Factor w/ 1 level "---": 1 1 1 1 1 1 1 1 1 1 ...
## $ Index : Factor w/ 1 level "---": 1 1 1 1 1 1 1 1 1 1 ...
## $ Dose : Factor w/ 1 level "---": 1 1 1 1 1 1 1 1 1 1 ...
## $ UV : Factor w/ 1 level "---": 1 1 1 1 1 1 1 1 1 1 ...
## $ D.D : Factor w/ 187 levels "0","0.001","0.002",..: 11 17 24 31 38 46 55 64 72 76 ...
## $ D.D.1 : Factor w/ 143 levels "0","0.001","0.002",..: 1 1 1 1 1 1 1 1 1 1 ...
## $ Temp.1 : Factor w/ 60 levels "20.8","21.1",..: 16 15 22 29 29 30 31 29 29 31 ...
## $ Hum.1 : int 30 32 32 32 32 33 33 33 33 34 ...
## $ Dew : Factor w/ 140 levels "1.9","10","10.1",..: 81 90 96 102 102 107 108 106 107 113 ...
## $ Heat : Factor w/ 72 levels "19.7","19.8",..: 15 15 25 33 33 35 36 34 35 37 ...
## $ X2nd : Factor w/ 252 levels "---","10.1","10.11",..: 1 1 1 1 1 1 1 1 1 1 ...
## $ X : Factor w/ 345 levels "---","1.1413",..: 1 1 1 1 1 1 1 1 1 1 ...
## $ X.1 : Factor w/ 1 level "---": 1 1 1 1 1 1 1 1 1 1 ...
## $ ET : int 0 0 0 0 0 0 0 0 0 0 ...
## $ Samp : int 225 702 702 703 702 703 702 703 702 702 ...
## $ Tx : int 1 1 1 1 1 1 1 1 1 1 ...
## $ Recept : Factor w/ 3 levels "100","32.9","63.3": 2 1 1 1 1 1 1 1 1 1 ...
## $ Int. : int 30 30 30 30 30 30 30 30 30 30 ...## Warning: Column `WindDir` joining factors with different levels, coercing
## to character vector
Graficos temporales
Grafico para todo el periodo (6/4 hasta 18/5) valores horarios
Grafico para todo el periodo (6/4 hasta 18/5) promedios diarios
Grafico de PM y CO para ver tendencia (6/4 hasta 18/5)
Graficos de NOx y O3 (6/4 hasta 18/5). HabrÃÂa que estudiar mejor la dinámica de la formación de ozono…?
Grafico para la semana del 12/4 de NOx y O3
Graficos de calendario
Promedios de 24 horas Para PM10
Para CO
Para NOx
Para O3
Grafico de la variación diaria y semanal de contaminantes para todo el periodo
Para PM10
Para CO
Para NOx
Para O3
Este grafico esta bien interesante porque da justamente al reves de lo experado… hay que analizar mejor que pasa con los precursores y tambien con la dinamica del viento, porque es como si se produjera el ozono en otro lado y volviera a la ciudad en la noche cuando cambia el viento (dinamica de la circulacion valle montaña)
Concentraciones normalizadas
