clipLevel2AM: Clip GEDI Elevation and Height Metrics by Coordinates

In carlos-alberto-silva/rGEDI: NASA's Global Ecosystem Dynamics Investigation (GEDI) Data Visualization and Processing

Clip GEDI Elevation and Height Metrics by Coordinates

Description

This function clips GEDI Level2A extracted Elevation and Height Metrics (Level2AM) within a given bounding coordinates

Usage

clipLevel2AM(level2AM, xmin, xmax, ymin, ymax)

Arguments

level2AM	A GEDI Level2A object (output of readLevel2A() function). An S4 object of class "gedi.level2a".
xmin	Numeric. West longitude (x) coordinate of bounding rectangle, in decimal degrees.
xmax	Numeric. East longitude (x) coordinate of bounding rectangle, in decimal degrees.
ymin	Numeric. South latitude (y) coordinate of bounding rectangle, in decimal degrees.
ymax	Numeric. North latitude (y) coordinate of bounding rectangle, in decimal degrees.

Value

Returns an S4 object of class data.table::data.table containing the clipped elevation and relative heights metrics.

See Also

https://lpdaac.usgs.gov/products/gedi02_av002/

Examples

# Specifying the path to GEDI level2A data (zip file)
outdir <- tempdir()
level2A_fp_zip <- system.file("extdata",
  "GEDI02_A_2019108080338_O01964_T05337_02_001_01_sub.zip",
  package = "rGEDI"
)

# Unzipping GEDI level2A data
level2Apath <- unzip(level2A_fp_zip, exdir = outdir)

# Reading GEDI level2A data (h5 file)
level2a <- readLevel2A(level2Apath = level2Apath)

# Extracting GEDI Elevation and Height Metrics
level2AM <- getLevel2AM(level2a)

# Bounding rectangle coordinates
xmin <- -44.15036
xmax <- -44.10066
ymin <- -13.75831
ymax <- -13.71244

# Clipping GEDI data by boundary box extent
level2AM_clip <- clipLevel2AM(level2AM, xmin, xmax, ymin, ymax)

close(level2a)

carlos-alberto-silva/rGEDI documentation built on Oct. 18, 2024, 4:46 a.m.