Learn R Programming

rGEDI (version 0.1.11)

getLevel1BWF: Get GEDI Pulse Full Waveform (GEDI Level1B)

Description

This function extracts the full waveform of a given pulse from GEDI Level1B data.

Usage

getLevel1BWF(level1b, shot_number)

Arguments

level1b

A GEDI Level1B object (output of readLevel1B function). A S4 object of class "gedi.level1b".

shot_number

Shot number. A scalar representing the shot number of a giving pulse.

Value

Returns an S4 object of class "gedi.fullwaveform".

Details

Shot numbers can be extracted using readLevel1B function.

See Also

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

Examples

Run this code
# NOT RUN {
# Specifying the path to GEDI level1B data (zip file)
outdir = tempdir()
level1B_fp_zip <- system.file("extdata",
                  "GEDI01_B_2019108080338_O01964_T05337_02_003_01_sub.zip",
                  package="rGEDI")

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

# Reading GEDI level1B data (h5 file)
level1b<-readLevel1B(level1Bpath=level1Bpath)

# Extracting GEDI full waveform for a giving shotnumber
wf <- getLevel1BWF(level1b, shot_number="19640521100108408")

# Plotting GEDI Full waveform
oldpar<-par()
par(mfrow = c(1,2), cex.axis = 1.5)
plot(wf, relative=FALSE, polygon=TRUE, type="l", lwd=2, col="forestgreen",
xlab="Waveform Amplitude", ylab="Elevation (m)")

plot(wf, relative=TRUE, polygon=TRUE, type="l", lwd=2, col="forestgreen",
xlab="Waveform Amplitude (%)", ylab="Elevation (m)")

par(oldpar)
close(level1b)
# }

Run the code above in your browser using DataLab