Create a daytime RasterStack for a specified date

Create a RasterStack from GOES AOD data files for the date specified by datetime. Each RasterLayer contains data from one Advanced Baseline Imager (ABI) scan during the specified time period.

If data for the specified time period is not found in the directory specified by setSatelliteDataDir(), it will be downloaded in order to create the RasterStack.

The Z axis of the RasterStack is a character vector where each element is the time stamp of the scan and has the format YYYYMMDDHHMMSS. This can be accessed using the raster::getZ() function. Names of the RasterStack are also time stamps of the scan, of the format XHH.MM.SS.

The dqfLevel parameter can take a value of:

• 0 -- High quality retrieval flag

• 1 -- Medium quality retrieval flag

• 2 -- Low quality retrieval flag

• 3 -- No retrieval quality flag

The bbox parameter can be a vector of floats in c(lonLo, lonHi, latLo, latHi) order or the return value from sp::bbox() or raster::extent().

goesaodc_createDaytimeRasterStack(
  satID = NULL,
  datetime = NULL,
  var = "AOD",
  res = 0.1,
  bbox = NULL,
  dqfLevel = NULL,
  timezone = NULL,
  verbose = FALSE
)

Arguments

satID	ID of the source GOES satellite (G16 or G17).
datetime	Desired datetime in any Ymd H [MS] format or POSIXct.
var	variable ("AOD, "DQF" or "ID")
res	resolution of raster in degrees
bbox	Bounding box for the region of interest.
dqfLevel	Data quality flag level.
timezone	timezone in which to interpret the datetime.
verbose	show progress of raster stacking.

Value

RasterStack

Examples

if (FALSE) {

library(MazamaSatelliteUtils)
library(MazamaSpatialUtils)

setSatelliteDataDir("~/Data/Satellite")
setSpatialDataDir("~/Data/Spatial")

bbox      <- c(-124, -120, 36, 39) # Kincade fire region
satID     <- "G16"
datetime  <- "2019-10-27"
timezone  <- "America/Los_Angeles"
dqfLevel  <- 2

latitude  <- 38.245   # \
                      #  - Town of Tomales
longitude <- -122.906 # /

dayStack <- goesaodc_createDaytimeRasterStack(
  satID = satID,
  datetime = datetime,
  timezone = timezone,
  bbox = bbox,
  dqfLevel = dqfLevel,
  verbose = TRUE
)

tb <- raster_createLocationTimeseries(dayStack,
                                      longitude = longitude,
                                      latitude = latitude,
                                      bbox = bbox)

plot(x = tb$datetime, y = tb$aod,
     pch = 15, cex = 0.8, col = rgb(red = 0, green = 0, blue = 0, alpha = 0.8),
     main = "Tomales, CA - 2019-10-27", xlab = "Time (UTC)", ylab = "AOD")
}