raster_subsetByPolygon: Subset a Raster* by polygon

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

Description

This function subsets a raster object using cookie cutter defined by polygon. All cells within the polygon are returned.

Usage

raster_subsetByPolygon(raster = NULL, polygon = NULL)

Arguments

raster	Raster\* object.
polygon	Spatial\* object.

Value

A Raster* object.

Examples

## Not run: 
library(AirFireModeling)
setModelDataDir('~/Data/BlueSky')

# Load model data
rasterList <- raster_load(
  modelName = "PNW-4km",
  modelRun = c(2020091300, 2020091400, 2020091500, 2020091600),
  xlim = c(-125, -116.5),
  ylim = c(42, 46.2)
)

library(MazamaSpatialUtils)
setSpatialDataDir("~/Data/Spatial")
loadSpatialData("USCensusStates")

OregonList <- raster_subsetByPolygon(
  rasterList,
  polygon = subset(USCensusStates, stateCode == "OR")
)

raster_ggmap(
  OregonList,
  palette = "Spectral",
  breaks = c(-1, 0, 4, 8, 16, 32, 64, 128, Inf),
  direction = -1
)

## End(Not run)

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