Access NASADEM Water Body layer (SWB)

[access-cloud-based-insar]

Access NASADEM Water Body layer (SWB)

[cmarshak]

Hey Konstantin,

Glad it's helpful.

Unfortunately, it's not in the scope of this package as is. I can share with you the rough idea of how you might approach this.

1. The urls will be the same as for NASADEM.

from dem_stitcher.datasets import get_overlapping_dem_tiles

urls = get_overlapping_dem_tiles([-121.5, 34.95, -120.2, 36.25], 'nasadem').url.tolist()

Note: if you unzip a url, you will find the *.swb layer in the unzipped directory. Here is an example url: https://e4ftl01.cr.usgs.gov/MEASURES/NASADEM_HGT.001/2000.02.11/NASADEM_HGT_n34w121.zip

2. Modify these mostly stand-alone functions (or reorganize them) to read the swb layer instead of the hgt layer. Might need a little bit of work (haven't tried it) because SRTM has its own gdal driver: https://gdal.org/drivers/raster/srtmhgt.html, but swb does not. You could infer the geometadata from the hgt file and then just generate a dataset.

Hope this helps. Good luck!

[konstantinklemmer]

Thanks for the swift reply!

[cmarshak]

Here is another option from @mgovorcin using ISCE2:

import isce
from isceobj.Alos2Proc.runDownloadDem import download_wbd
# Download SRTM-SWDB water mask
# ISCE2 requires South, North, East, West extents (or ymin, ymax, xmin, xmax)
mask_filename = download_wbd(extent[1], extent[3], extent[0], extent[2])