clipLevel2A: Clip GEDI Level2A data by Coordinates

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

Clip GEDI Level2A data by Coordinates

Description

This function clips GEDI Level2A data within a given bounding coordinates

Usage

clipLevel2A(level2a, xmin, xmax, ymin, ymax, output)

Arguments

level2a	A GEDI Level2A object (output of readLevel2A() function). An S4 object of class "gedi.level2a".
xmin	Numeric. West longitude (x) coordinate of the bounding rectangle, in decimal degrees.
xmax	Numeric. East longitude (x) coordinate of the bounding rectangle, in decimal degrees.
ymin	Numeric. South latitude (y) coordinate of the bounding rectangle, in decimal degrees.
ymax	Numeric. North latitude (y) coordinate of the bounding rectangle, in decimal degrees.
output	Optional character path where to save the new hdf5file. The default stores a temporary file only.

Value

Returns a list of S4 objects of class "gedi.level2a" containing clipped GEDI Level2A data.

See Also

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

Examples

outdir <- tempdir()

# Specifying the path to GEDI level2A data (zip file)
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)

# Bounding rectangle coordinates
xmin <- -44.13
xmax <- -44.12
ymin <- -13.74
ymax <- -13.73

print(level2a)

# Specifying output file and path
output <- file.path(outdir, "GEDI02_A_2019108080338_O01964_T05337_02_001_01_clip.h5")

# Clipping GEDI Level2A data by boundary box extent
level2a_clip <- clipLevel2A(level2a, xmin, xmax, ymin, ymax, output)

close(level2a)
close(level2a_clip)

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