/TanDEM-SRTM-dh

Python scripts for the correction of raw SRTM-C tiles with TanDEM-X and potential vertical land-level change mapping

Primary LanguagePythonGNU General Public License v3.0GPL-3.0

title date author
TanDEM-X, SRTM-C Correction and Differencing
June 2018
Ben Purinton ([purinton@uni-potsdam.de](purinton@uni-potsdam.de))

TanDEM-X, SRTM-C Correction and Differencing

These Python codes are intended for the correction (co-registration, destriping, block shifting) of raw SRTM-C tiles to raw TanDEM-X tiles and also for the mapping of potential vertical land-level changes in conjunction with TanDEM-X auxiliary rasters. The companion paper is:

Purinton, B., and Bookhagen, B.: Measuring Decadal Vertical Land-level Changes from SRTM-C (2000) and TanDEM-X (~2015) in the South-Central Andes, Earth Surface Dynamics, 2018. (https://doi.org/10.5194/esurf-6-971-2018)

Processing

The scripts assume some knowledge of Python coding and access to the following packages: numpy, scipy, osgeo (gdal + gdalnumeric), matplotlib, scikit-image.

If you have trouble with installing or importing the scripts contact me at (purinton@uni-potsdam.de)

Below a suggested processing routine:

  1. Create a working directory with a folder for SRTM-C tiles ("/path/to/working/directory/srtms/"), TanDEM-X tiles ("/path/to/working/directory/tandems/"), and TanDEM-X auxiliary rasters ("/path/to/working/directory/tandems/auxiliary/") including the water indication mask (WAM), height error mask (HEM) and consistency mask (COM).

  2. Make sure all tiles and auxiliary rasters are in WGS84 unprojected geographic coordinates and are of the same resolution and dimensions and place them in the appropriate working directory outlined in (1).

  3. Run the scripts in order (01-05) setting the variable names appropriately before running each step. Note the instructions in comments for each variable. Variables can be found in each script at line numbers:

    • 01_coregistration.py - 43-69
    • 02_fft_destriping.py - 33-53
    • 03_blocked_shifting.py - 36-61
    • 04_change_mapping_channels.py - 42-80
    • 05_change_mapping_hillslopes.py - 33-62