Get GEDI Pulse Full Waveform (GEDI Level1B)
This function extracts the full waveform of a given pulse from GEDI Level1B data.
getLevel1BWF(level1b, shot_number)
level1b |
A GEDI Level1B object (output of |
shot_number |
Shot number. A scalar representing the shot number of a giving pulse. |
Shot numbers can be extracted using readLevel1B function.
Returns an S4 object of class "gedi.fullwaveform".
https://lpdaac.usgs.gov/products/gedi01_bv001/
# 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)Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.