Retrieve hourly averages, daily maximums, or daily minimums of pollution data in the original units, by station, from the air quality server at Consulta de Concentraciones
Arguments
- criterion
Type of data to download.
HORARIOS - Hourly data
MAXIMOS" - Daily maximums
MINIMOS - Daily minimums
- pollutant
The type of pollutant to download.
SO2 - Sulfur Dioxide (parts per billion)
CO - Carbon Monoxide (parts per million)
NOX - Nitrogen Oxides (parts per billion)
NO2 - Nitrogen Dioxide (parts per billion)
NO - Nitric Oxide (parts per billion)
O3 - Ozone (parts per billion)
PM10 - Particulate matter 10 micrometers or less (micrograms per cubic meter)
PM25 - Particulate matter 2.5 micrometers or less (micrograms per cubic meter)
WSP - Wind velocity (meters per second)
WDR - Wind direction (degrees)
TMP - Temperature (degrees Celsius)
RH - Relative humidity (percentage)
- year
an integer indicating the year for which to download data (the earliest possible value is 1986)
- month
month number to download
Value
A data.frame with pollution data, the hours correspond to the Etc/GMT+6 timezone, with no daylight saving time
Details
Temperature (TMP) data was rounded to the nearest integer, but the
get_station_data function allows you to download data accurate
to one decimal point in some cases (i.e. for old data).
See also
stations for a data.frame with the location and names
of all pollution measuring stations
Other raw data functions:
get_station_data()
Examples
if (FALSE) { # \dontrun{
## Download daily hourly PM10 data (particulate matter 10 micrometers or
## less in diameter) from March 2016
df_pm10 <- get_station_month_data("HORARIOS", "PM10", 2016, 3)
head(df_pm10)
## Download daily hourly O3 data from October 2017
df_o3 <- get_station_month_data("HORARIOS", "O3", 2018, 1)
## Convert to local Mexico City time
df_o3$mxc_time <- format(as.POSIXct(paste0(df_o3$date,
" ",
df_o3$hour, ":00"),
tz = "Etc/GMT+6"),
tz = "America/Mexico_City")
head(df_o3)
} # }