Creates a time series plot of PM2.5 data from a ws_monitor object (see note below). Optional arguments color code by AQI index, add shading to indicate nighttime, and adjust the time display (local vs. UTC).
When a named style
is used, some graphical parameters will be
overridden. Available styles include:
aqidots
-- hourly values are individually colored by 24-hr
AQI levels
gnats
-- semi-transparent dots like a cloud of gnats
monitor_timeseriesPlot(ws_monitor, monitorID = NULL, tlim = NULL,
localTime = TRUE, style = NULL, shadedNight = FALSE, add = FALSE,
gridPos = "", gridCol = "black", gridLwd = 1, gridLty = "solid",
dayLwd = 0, hourLwd = 0, hourInterval = 6, ...)
ws_monitor object.
Monitor ID for one or more monitor in the ws_monitor object.
Optional vector with start and end times (integer or character representing YYYYMMDD[HH]).
Logical specifying whether tlim
is in local time or
UTC.
Custom styling, one of "aqidots"
.
Add nighttime shading.
Logical specifying whether to add to the current plot.
Position of grid lines either "over", "under" ("" for no grid lines).
Grid line color.
Grid line width.
Grid line type.
Day marker line width.
Hour marker line width.
Interval for grid (max = 12).
Additional arguments to be passed to points()
.
# NOT RUN {
N_M <- Northwest_Megafires
# monitor_leaflet(N_M) # to identify Spokane monitorIDs
Spokane <- monitor_subsetBy(
N_M,
stringr::str_detect(N_M$meta$monitorID, "^53063")
)
monitor_timeseriesPlot(Spokane, style = "gnats")
title("Spokane PM2.5 values, 2015")
monitor_timeseriesPlot(
Spokane,
tlim = c(20150801, 20150831),
style = "aqidots",
pch = 16
)
addAQILegend()
title("Spokane PM2.5 values, August 2015")
monitor_timeseriesPlot(
Spokane,
tlim = c(20150821, 20150828),
shadedNight = TRUE,
style = "gnats"
)
abline(h = AQI$breaks_24, col = AQI$colors, lwd = 2)
addAQILegend()
title("Spokane PM2.5 values, August 2015")
# }
Run the code above in your browser using DataLab