Clip GEDI Plant Area Volume Density profile by geometry
This function clips GEDI level2B derived Plant Area Index profile within a given geometry
clipLevel2BPAVDProfileGeometry(level2BPAVDProfile, polygon_spdf, split_by)
level2BPAVDProfile |
A GEDI Level2B object (output of |
polygon_spdf |
Polygon. An object of class |
split_by |
Polygon id. If defined, GEDI data will be clipped by each polygon using the attribute specified by |
Returns an S4 object of class data.table-class
containing the Plant Area Volume Density profile data.
https://lpdaac.usgs.gov/products/gedi02_bv001/
# Specifying the path to GEDI level2B data (zip file) outdir = tempdir() level2B_fp_zip <- system.file("extdata", "GEDI02_B_2019108080338_O01964_T05337_02_001_01_sub.zip", package="rGEDI") # Unzipping GEDI level2A data level2Bpath <- unzip(level2B_fp_zip,exdir = outdir) # Reading GEDI level2B data (h5 file) level2b<-readLevel2B(level2Bpath=level2Bpath) # Extracting GEDI Plant Area Volume Density profile level2BPAVDProfile<-getLevel2BPAVDProfile(level2b) # Specifying the path to shapefile polygon_filepath <- system.file("extdata", "stands_cerrado.shp", package="rGEDI") # Reading shapefile as SpatialPolygonsDataFrame object library(raster) polygon_spdf<-shapefile(polygon_filepath) # Clipping GEDI Plant Area Volume Density profile by geometry level2BPAVDProfile_clip <- clipLevel2BPAVDProfileGeometry( level2BPAVDProfile, polygon_spdf, split_by="id") close(level2b)
Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.