| title | date | author |
|---|---|---|
TanDEM-X, SRTM-C Correction and Differencing |
June 2018 |
Ben Purinton ([[email protected]]([email protected])) |
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)
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 ([email protected])
Below a suggested processing routine:
-
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).
-
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).
-
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