raster_toMonitor: Create a PWFSLSmoke ws_monitor object from a Raster\* object

In MazamaScience/AirFireModeling: Utilities for Working with Monitoring Data and Smoke Model Output

Description

Time series associated with multiple grid cells are converted into a PWFSLSmoke ws_monitor. The monitorID associated with each location is generated with:

  MazamaLocationUtils::location_createID(x, y)

Usage

raster_toMonitor(
  raster,
  longitude = NULL,
  latitude = NULL,
  radius = 5,
  count = NULL,
  rasterName = NULL,
  verbose = TRUE
)

Arguments

raster	A Raster* object.
longitude	Target longitude from which the radius will be calculated.
latitude	Target latitude from which the radius will be calculated.
radius	Distance (km) of radius from target location.
count	Number of grid cells within radius to return.
rasterName	Optional string prepended to monitorIDs.
verbose	Logical to display messages.

Value

A PWFSLSmoke ws_monitor object representing a single monitor.

Examples

library(AirFireModeling)
setModelDataDir('~/Data/BlueSky')

# Load model data
rasterList <- raster_load(
  modelName = c("PNW-1.33km", "PNW-4km"),
  modelRun = c(2020091300),
  xlim = c(-125, -115),
  ylim = c(42, 50)
)

# MazamaSpatialUtils are needed to determine timezone and state
library(MazamaSpatialUtils)
PWFSLSmoke::initializeMazamaSpatialUtils()

model_Portland <- raster_toMonitor(
  rasterList,
  longitude = -122.68,
  latitude = 45.52,
  radius = 5
)

PWFSLSmoke::monitor_timeseriesPlot(
  model_Portland[[1]],
  type = 'l', col = 'salmon'
)
PWFSLSmoke::monitor_timeseriesPlot(
  model_Portland[[2]],
  type = 'l', col = 'dodgerblue', add = TRUE
)

legend(
  x = "topright",
  legend = names(model_Portland),
  lwd = c(1, 1),
  col= c("salmon", "dodgerblue")
)

MazamaScience/AirFireModeling documentation built on March 13, 2021, 12:02 a.m.