Create a tibble from ncdf4 handle

Create a tibble with columns: AOD, DQF, lon and lat. This information is sufficient to plot as points or create a raster object.

goesaodc_createTibble(nc = NULL, bbox = bbox_CONUS, verbose = FALSE)

Arguments

nc	ncdf4 handle or a list of handles.
bbox	Geographic extent of area of interest; Defaults to CONUS.
verbose	Logical flag to increase messages while processing data.

Value

Tibble (dataframe) with NetCDF variables and associated locations.

Examples

# \donttest{
# Tibble based on full extent of Gridfile
library(MazamaSatelliteUtils)

setSatelliteDataDir("~/Data/Satellite")

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

nc <- goesaodc_openFile(ncFile)

tbl <- goesaodc_createTibble(nc)
head(tbl)
#> # A tibble: 6 x 4
#>     lon   lat   AOD   DQF
#>   <dbl> <dbl> <dbl> <dbl>
#> 1 -125.  50.0 0.113     2
#> 2 -125.  50.0 0.667     2
#> 3 -124.  50.0 1.45      2
#> 4 -125.  49.9 0.138     2
#> 5 -125.  49.9 0.393     2
#> 6 -125.  50.0 0.849     2

# Tibble based on BBOX filtered extent of tibble
library(MazamaSatelliteUtils)

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

nc <- goesaodc_openFile(ncFile)

# 2019 Kincade fire in Sonoma county
bbox <- c(-124, -120, 36, 39)

filtered_tbl <- goesaodc_createTibble(nc, bbox)
head(filtered_tbl)
#> # A tibble: 6 x 4
#>     lon   lat    AOD   DQF
#>   <dbl> <dbl>  <dbl> <dbl>
#> 1 -124.  39.0 0.0371     1
#> 2 -124.  39.0 0.0234     1
#> 3 -124.  39.0 0.0159     1
#> 4 -124.  39.0 0.0304     1
#> 5 -124.  39.0 0.0368     0
#> 6 -124.  39.0 0.0283     1

# }