plotWFMetrics: GEDI full waveform plot with metrics

In rGEDI: NASA's Global Ecosystem Dynamics Investigation (GEDI) Data Visualization and Processing

Description

Plots the waveform with overlaid RH metrics

Usage

plotWFMetrics(level1b, level2a, shot_number, rh=c(25, 50, 75),...)

Arguments

level1b	A GEDI Level1B object (output of readLevel1B function). An S4 object of class "gedi.level1b".
level2a	A GEDI Level2A object (output of readLevel2A function). An S4 object of class "gedi.level2a".
shot_number	Shot number. A scalar representing the shot number of a giving pulse.
rh	Integer vector. Specify which RH metrics to plot except rh0 and rh100, default c(25, 50, 75).
...	Will be passed to the main plot.

Value

Returns a raster layer(s) of selected GEDI Canopy Cover and Vertical Profile Metric(s)

See Also

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

Examples

# specify the path to GEDI level1B and Level2A data (zip file)
outdir = tempdir()
level1B_fp_zip <- system.file("extdata",
                  "GEDI01_B_2019108080338_O01964_T05337_02_003_01_sub.zip",
                  package="rGEDI")

level2A_fp_zip <- system.file("extdata",
                  "GEDI02_A_2019108080338_O01964_T05337_02_001_01_sub.zip",
                  package="rGEDI")

# Unzipping GEDI level1B data
level1Bpath <- unzip(level1B_fp_zip,exdir = outdir)
level2Apath <- unzip(level2A_fp_zip,exdir = outdir)

# Reading GEDI level1B and Level2A data (h5 file)
level1b<-readLevel1B(level1Bpath=level1Bpath)
level2a<-readLevel2A(level2Apath=level2Apath)

shot_number = "19640521100108408"

plotWFMetrics(level1b, level2a, shot_number, rh=c(25, 50, 75, 90))

close(level1b)
close(level2a)

rGEDI documentation built on Jan. 19, 2021, 3:01 p.m.