rGEDI (version 0.1.11)

plotWFMetrics: GEDI full waveform plot with metrics

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

# NOT RUN {
# 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)
# }