localExamples/2019_WESTAR.R

In MazamaScience/PWFSLSmoke: Utilities for Working with Air Quality Monitoring Data

li# ----- Chetco Bar -----------------------------------------------

conus <-
  monitor_loadLatest() %>%
  monitor_subset(stateCodes = CONUS)

monitor_leaflet(conus)

Shoshone_Bannock <-
  monitor_subset(conus, monitorIDs = "160050020_01")

monitor_timeseriesPlot(Shoshone_Bannock)












# ----- Sacramento smoke -------------------------------------------------------

library(PWFSLSmoke)

camp_fire <-
  monitor_loadAnnual(2018) %>%
  monitor_subset(stateCodes = 'CA') %>%
  monitor_subset(tlim = c(20181108,20181123))

monitor_leaflet(camp_fire)

Sacramento <-
  camp_fire %>%
  monitor_subset(monitorIDs = '060670010_01')

monitor_timeseriesPlot(
  Sacramento,
  style='aqidots',
  pch=16,
  xlab="2017"
)
addAQIStackedBar()
addAQILines()
addAQILegend()
title("Sacramento Smoke")

Sacramento_area <-
  camp_fire %>%
  monitor_subsetByDistance(
    longitude = Sacramento$meta$longitude,
    latitude = Sacramento$meta$latitude,
    radius = 50
  )

monitor_leaflet(Sacramento_area)

Sacramento_area %>%
  monitor_collapse() %>%
  monitor_dailyStatistic() %>%
  monitor_extractData()


# ----- Kennewick School -------------------------------------------------------

library(PWFSLSmoke)
library(AirMonitorPlots)

# Get the appropriate data
wa_2017 <-
  monitor_loadAnnual(2017) %>%
  monitor_subset(stateCodes = 'WA') %>%
  monitor_subset(tlim = c(20170901,20170914))

# Find the monitorID for Kennewick
monitor_leaflet(wa_2017)

# Create a Kennewick-only ws_monitor object
Kennewick <-
  wa_2017 %>%
  monitor_subset(monitorIDs = '530050002_01')

# Create a plot
ggplot_pm25Timeseries(Kennewick) +
  stat_dailyAQCategory(alpha = .5) +
  geom_pm25Points(shape = "square", alpha = 1.0) +
  custom_aqiStackedBar(width = 0.01) +
  ggtitle("Kennewick, September 2017")

# Daily/Hourly barplot
ggplot_pm25Timeseries(Kennewick) +
  stat_AQCategory(color = NA) +
  stat_dailyAQCategory(alpha = .5, timezone = "America/Los_Angeles") +
  facet_grid(rows = vars(monitorID)) +
  ggtitle("Kennewick, September 2017")

Toppenish <-
  wa_2017 %>%
  monitor_subset(monitorIDs = "530770015_01")

Yakima_Basin <-
  wa_2017 %>%
  monitor_subsetByDistance(longitude = Toppenish$meta$longitude,
                           latitude = Toppenish$meta$latitude,
                           radius=50)

# Daily/Hourly barplot
ggplot_pm25Timeseries(Yakima_Basin) +
  stat_AQCategory(color = NA) +
  stat_dailyAQCategory(alpha = .5, timezone = "America/Los_Angeles") +
  facet_grid(rows = vars(monitorID)) +
  ggtitle("Kennewick, September 2017")