goesaodc_createNativeGrid: Create a nativeGrid object

In MazamaScience/MazamaSatelliteUtils: Utilities for Working with Model and Satellite Data as Rasters

Description

Creates a nativeGrid object with a list of matrices for: longitude, latitude, AOD values, and DQF values. The arrays are defined in native i, j coordinates and are thus curvilinear as opposed to rectilinear geospatial coordinates.

The nc parameter can be either a single netCDF handle or a list of handles. If a list of handles is provided, grid cell values for AOD and DQF will be averaged across all nc handles. This "native grid" averaging provides a simple way to convert 5-minute data into an hourly average.

Usage

goesaodc_createNativeGrid(nc = NULL, bbox = bbox_CONUS)

Arguments

nc	ncdf4 handle or a list of handles.
bbox	Bounding box for the region of interest. Defaults to CONUS.

Value

List with lon, lat, AOD and DQF matrices

Examples

library(MazamaSatelliteUtils)

setSatelliteDataDir("~/Data/Satellite")

scanFile <- goesaodc_downloadScanFiles(
  satID = "G17", 
  datetime = "2019-10-27 14:00",
  timezone = "America/Los_Angeles"
)

nc <- goesaodc_openScanFile(scanFile)

# Kincade fire region
kincade_bbox <- c(-124, -120, 36, 39)

layout(matrix(seq(2)))

nativeGrid <- goesaodc_createNativeGrid(nc, kincade_bbox)
image(nativeGrid$AOD[,ncol(nativeGrid$AOD):1])
title("Single timestep")

nativeGrid <- goesaodc_createNativeGrid(nc, kincade_bbox)
image(nativeGrid$AOD[,ncol(nativeGrid$AOD):1])
title("Average of 12 timesteps")

layout(1)

MazamaScience/MazamaSatelliteUtils documentation built on Dec. 17, 2021, 3:20 a.m.