The new Global Ecosystem Dynamics Investigation GEDI launched on December 5th, 2018 is operating onboard the International Space Station ISS producing a high-resolution laser collecting 3D data around Earth on forest canopy height, canopy vertical structure, and surface elevation. As it is collecting daily data, a stable and fast platform is essential. For this reason, the library pyGEDI is developed in Python as it can utilize multiple CPUs, GPUs, and is supported by C and GDAL.
pyGEDI provides a high performance, lower cognitive load, and cleaner and more transparent code for data extraction, analysis, processing, and visualization of GEDI's products.
pyGEDI has multiple functions for visualization, processing, analysis, and data extraction all in one package. pyGEDI package contains the following functions:
- Connection to NASA's server.
- Downloading GEDI data.
- Clipping of your specific area.
- Generate multiple files like
.csv
,.shp
, and.tif
for GEDI's products. - Visualization of waveforms with canopy height and profile metrics.
- Process
.h5
files. - Generate histograms, raster’s of different GEDI metrics.
- Graph in 3D cloud points.
- Demonstrative products (correlation GEDI and Airborne Laser Scanning).
numpy
gdal
h5py
pandas
matplotlib
pip install pyGEDI
from pyGEDI import *
Connect with NASA's server Register for an Earthdata Login Profile (if you do not have an account register for a free account by providing your email).
username='write your user name'
password='write your password'
session=sessionNASA(username,password)
Define the coordinates of your specific area of interest. The next line will exhibit an example using coordinates in Colombia.
ul_lat= 2.96845
ul_lon=-73.32586
lr_lat=-1.26845
lr_lon=-70.23869
bbox=[ul_lat,ul_lon,lr_lat,lr_lon]
Now that the box is defined, download all GEDI trajectories that have flown over the area(s) of interest. Currently, there are only three products available with 25m resolution.
This includes:
- Level 1 - Geolocated Waveforms
- Level 2- Footprint Level Canopy Height
- Level 2- Profile Metrics
For more information check GEDI Finder
Note: GEDI Trajectories contain heavy amounts of data. It is recommended to have good internet for it to download in a more timely manner. Depending on your area of interest and trajectories this may take a couple hours illustrated by the download bar.
product_1B='GEDI01_B'
product_2A='GEDI02_A'
product_2B='GEDI02_B'
version='001'
outdir_1B='data/'+product_1B+'.'+version+'/'
outdir_2A='data/'+product_2A+'.'+version+'/'
outdir_2B='data/'+product_2B+'.'+version+'/'
gediDownload(outdir_1B,product_1B,version,bbox,session)
gediDownload(outdir_2A,product_2A,version,bbox,session)
gediDownload(outdir_2B,product_2B,version,bbox,session)
Once the download starts it will show a download bar for each file
Created the subdirectory
data/GEDI01_B.001/2019.08.05/
GEDI01_B_2019217151359_O03661_T02309_02_003_01.h5 | 7.405GB | 24.19%
[███████████████████......................................................]
Access the following notebooks for each respective product.
Notebook:
- GEDI 2A or Colab.research
- GEDI 2B or Colab.research
- WaveForm GEDI 1B,2A,2B or Colab.research
- Demonstrative products or Colab.research
- Time series.
- Export data to Google Earth Engine.
- Analysis for LEVEL 3 gridded canopy height metrics and variability.
- Analysis for LEVEL 4A and 4B Footprint and gridded aboveground carbon estimates.
- GEDI webpage. Accessed on April 4th 2020 https://gedi.umd.edu/
- GEDI 1B version 001. Accessed on April 4th 2020 https://lpdaac.usgs.gov/products/gedi01_bv001/
- GEDI 2A version 001. Accessed on April 4th 2020 https://lpdaac.usgs.gov/products/gedi02_av001/
- GEDI 2B version 001. Accessed on April 4th 2020 https://lpdaac.usgs.gov/products/gedi02_bv001/
- GEDI Finder. Accessed on April 4th 2020 https://lpdaacsvc.cr.usgs.gov/services/gedifinder
- G-Liht webpage. Accessed on April 4th 2020 https://gliht.gsfc.nasa.gov/
- The files
.h5
for the notebooks are temporally taken from: https://github.com/carlos-alberto-silva/rGEDI/tree/master/inst/extdata
We would like to thank The University of Maryland and NASA's Goddard Space Flight Center for developing GEDI's mission and for providing free and open data.
Please report any issue regarding the pyGEDI package to:
- 137eduin[at]gmail[dot]com
- andreshs[at]umd[dot]edu
pyGEDI blog forthcoming.
pyGEDI appreciates help from a wide range of different backgrounds. Small improvements or fixes are always appreciated. Kindly report any issues with labeling or processing. If you are considering larger contributions outside the traditional coding work, please contact us through the mailing list.
pyGEDI has several functions and is a continuously improved package. If it does not exist in pyGEDI please ask us for a new feature, and depending on the feasibility we will be glad to implement your requested feature.
Eduin H.SERNA; Andres Hernandez-Serna. pyGEDI: NASA's Global Ecosystem Dynamics Investigation (GEDI) mission data extraction, analysis, processing and visualization. version 0.2, April. 5th 2020, available at: https://pypi.org/project/pyGEDI/
pyGEDI package has not been developed by the GEDI team. It comes with no guarantee, expressed or implied, and the authors hold no responsibility for its use or reliability of its outputs.