topohack's Introduction

Topohack

Taking ICESat-2 to the mountains: a workflow using satellite laser altimetry to resolve topography over complex terrain

General Objective

Compare and evaluate ICESat-2 data with high resolution DEMs (airborne lidar/satellite stereo) collected at lower latitudes over bare ground.

Collaborators

Vibhor Agarwal
Michelle Hu
Friedrich Knuth
HP Marshall
Justin Pflug
Mariama Dryak
Will Kochtitzky

Team Lead:

Shashank Bhushan

Data Science Leads:

Joachim Meyer
Amy Steicker

Datasets

ICESat-2 ATL06 (20 m resolution)
ICESat-2 ATL03 (geolocated points along track)
Bareground data (LULC for US, Global Bareground Data Product for areas outside US)

Validation Datasets

Tools

Python: geopandas, rasterio, numpy, scipy, pandas, pygeotools, demcoreg
NASA Ames Stereo Pipeline

High-level Goals

Learn how to download the ICESat-2 data by lat lon bounding box
Create library with some basic convenience functions
Explore OpenAltimetry and utility with multiple locations and time periods
Explore ATL03 and ATL06 products and theoretical basis
Learn how to subset ATL03 and ATL06 data based on flags
Intersect ICESat-2 tracks with RGI glacier polygons to get a sense of bare ground coverage near glaciers.
Evaluate/Compare the topography resolved by ICESat-2 profiles along steep mountains with topographic profiles returned from high-resolution DEMs.
Get a sense of snow accumulation (depth) by comparing Snow-off DEM over Grand Mesa with winter (October to February) IceSat-2 collects.
Create notebook that shows how to pull all cloud-free data for entire mission for target area bounding box.

Study Sites:

Dependent on ICESat-2 coverage.
Potential study sites: Cascades and Olympic Ranges (western WA), Rocky Mountains (CO), High Mountain Asia, Grand Mesa (CO)

Team Wiki

Further information can be found on our team wiki page

Test Drive

To launch this repository on a Pangeo Binder instance and try out topolib using the examples under notebooks/, click here:

topohack's People

Contributors

Stargazers

Watchers

topohack's Issues

TypeError: 'numpy.float64' object does not support item assignment

Hello,

Runs fine until this point:

gdf_list = [(gda_lib.ATL08_2_gdf(x,dataset_dict)) for x in ATL08_list]
gdf_colombia = gda_lib.concat_gdf(gdf_list)

TypeError Traceback (most recent call last)
d:\DRIVE\SCIENCE_Projects\ICESATDEM\IS2_DEM_comparison_WIP.ipynb Cell 40 in 1
----> 1 gdf_list = [(gda_lib.ATL08_2_gdf(x,dataset_dict)) for x in ATL08_list]
2 gdf_colombia = gda_lib.concat_gdf(gdf_list)

d:\DRIVE\SCIENCE_Projects\ICESATDEM\IS2_DEM_comparison_WIP.ipynb Cell 40 in 1
----> 1 gdf_list = [(gda_lib.ATL08_2_gdf(x,dataset_dict)) for x in ATL08_list]
2 gdf_colombia = gda_lib.concat_gdf(gdf_list)

File d:\drive\github\topohack\topolib\gda_lib.py:129, in ATL08_2_gdf(ATL06_fn, dataset_dict)
126 dataset_dict['land_segments'].append('longitude')
127 #use Ben's Scripts to convert to dict
--> 129 data_dict = ATL08_to_dict(ATL06_fn,dataset_dict)
130 #this will give us 6 tracks
131 i = 0

File d:\drive\github\topohack\topolib\gda_lib.py:105, in ATL08_to_dict(filename, dataset_dict)
103 bad = temp[dataset]==fill_value
104 temp[dataset]=np.float64(temp[dataset])
--> 105 temp[dataset][bad]=np.NaN
106 except KeyError as e:
107 pass