Create a SpatialPointsDataFrame of GOES data

Create a SpatialPointsDataFrame of GOES data including data within the specified bounding box and data quality flag level.

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_createSpatialPoints(
  nc = NULL,
  bbox = c(-125, -65, 24, 50),
  dqfLevel = NULL
)

Arguments

nc

ncdf4 handle.
bbox

Bounding box for the region of interest.
dqfLevel

Data quality flag level.

Value

SpatialPointsDataFrame

Examples

# \donttest{
library(MazamaSatelliteUtils)
setSatelliteDataDir("~/Data/Satellite")

netCDF <- goesaodc_downloadAOD(
  satID = "G17",
  datetime = "2019-10-27 10",
  timezone = "America/Los_Angeles")[1]

kincade_bbox <- c(-126, -119, 35, 40)

nc <- goesaodc_openFile(netCDF)
sp <- goesaodc_createSpatialPoints(nc, dqfLevel = 3, bbox = kincade_bbox)
maps::map(database = "state", "regions" = c("california"), xlim = c(-126, -113.5))
goesaodc_plotSpatialPoints(sp, cex = 0.2, add = TRUE)
# }