clipLevel1B: Clip GEDI Level1B data by Coordinates

Description Usage Arguments Value See Also Examples

View source: R/clipLevel1B.R

Description

This function clips GEDI Level1B data (geolocated waveforms) within a given bounding coordinates

Usage

1
clipLevel1B(level1b, xmin, xmax, ymin, ymax, output)

Arguments

level1b

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

xmin

Numeric. West longitude (x) coordinate of the bounding rectangle, in decimal degrees.

xmax

Numeric. East longitude (x) coordinate of the bounding rectangle, in decimal degrees.

ymin

Numeric. South latitude (y) coordinate of the bounding rectangle, in decimal degrees.

ymax

Numeric. North latitude (y) coordinate of the bounding rectangle, in decimal degrees.

output

Optional character path where to save the new hdf5file. The default stores a temporary file only.

Value

Returns a list of S4 objects of class "gedi.level1b" containing clipped GEDI Level1B data.

See Also

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

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
# 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)

# Bounding rectangle coordinates
xmin=-44.13
xmax=-44.12
ymin=-13.74
ymax=-13.73

# Specifying output file and path
output<-file.path(outdir,"GEDI01_B_2019108080338_O01964_T05337_02_003_01_clip")

# Clipping GEDI Level1B data by extent boundary box
level1b_clip <- clipLevel1B(level1b,xmin, xmax, ymin, ymax,output)

close(level1b)
close(level1b_clip)

rGEDI documentation built on Jan. 21, 2021, 1:06 a.m.