Package: demcoreg
demcoregPython and shell scripts for co-registration of rasters, specifically horizontal and vertical alignment of digital elevation models (DEMs).OverviewAll DEMs have some horizontal and vertical geolocation error.  It is important to remove relative offsets when differencing DEMs for elevation change analyses.  These tools offer several options to solve this problem.  Most solve for the sub-pixel horizontal shift and vertical offset required to minimize errors over "static" control surfaces.  The ASPpc_aligntool can also solve for more complex transformations with rotations and scaling.FeaturesMultiple co-registration algorithms:Automated, iterative implementation of the algorithm outlined byNuth and Kaab (2011)Wrappers around theAmes Stereo PipelineIterative Closest Point (ICP)pc_alignutility (https://stereopipeline.readthedocs.io/en/latest/tools/pc_align.html)Normalized cross-correlation (NCC) with sub-pixel refinmentSum of absolute differences (SAD)Command-line utilities for raster differencing with necessary resampling (compute_diff.py), raster stats (robust_stats.py), and raster sampling at point locations (sample_raster_at_pts.py)Mask preparation and automatic determination of static control surfaces (i.e., exposed bedrock) for a user-specified combination of:RGI glacier polygonsLand-use/Land-cover classification:National Land Cover Database (NLCD) 30-m products for CONUSwith pre-configured filter combinations (e.g., 'not_forest+not_water')Global Bare Ground 30-m productsSnow mask:Thresholded MODSCAG fSCA from ~2-week period around DEM timestampThresholded SNODAS modelThresholded Top-of-atmosphere reflectance values from corresponding orthoimage (requires pregeneration)Some useful command-line utilities (run with-hoption for complete usage)dem_align.py- robust raster DEM co-registration (e.g., Nuth and Kaab [2011]) for surfaces with variable slope and aspect (e.g., mountains)dem_mask.py- pre-generate mask to identify "stable" surfaces to use during co-registrationpc_align_wrapper.sh- wrapper around NASA Ames Stereo Pipeline pc_align utility for iterative closest point co-registrationapply_dem_translation.py- update raster geotransform and remove vertical offsetcompute_diff.py- simple DEM difference calculation with intuitive resampling optionsrobust_stats.py- print out robust raster statistics (e.g,. for DEM difference map before/after co-registration)Sample outputdem_align.pySample command:dem_align.py ref_dem.tif src_dem.tifdem_mask.pySample command:dem_mask.py --toa --bareground --glaciers src_dem.tiffilter_glas.py outputExample applicationsHigh-mountain AsiaCo-registration of ~35000 high-resolution DEMs from multiple sensorshttps://github.com/dshean/hma_mb_paperPine Island Glacier, AntarcticaLeast-squares optimization to correct for offset and "tilt" of ~800 high-resolution DEMs with limited ground controlhttps://github.com/dshean/pig_dem_meltrateInstallationWe are hoping to clean up the code, remove unnecessary dependencies, and streamline installation using conda. For now, we recommend following the "Building from Latest Source" instructions below, to obtain latest features/bugfixes.If unfamiliar with this process, or if you are new to Python, bash, and/or git/github, start with these more detailed instructions and notes:Beginner's guide for installation and basic usageBuilding from Latest Source (recommended)Assuming you have working Python3 install with GDAL and NumPy, installpygeotoolsClone thedemcoregrepository:git clone https://github.com/dshean/demcoreg.gitPerform developer install with pip:pip install -e demcoregThe -e flag ("editable mode", setuptools "develop mode") will allow you to modify source code and immediately see changes. Useful if you need to make minor tweaks or bugfixes (please submit a Pull Request!)Optionally, append the demcoreg subdirectory containing scripts to your PATH:export PATH=${PATH}:$PWD/demcoreg/demcoreg(replacing$PWDwith the absolute path to the cloned demcoreg repository)To make this permanent, add that line to your shell config file (e.g., ~/.bashrc).Simple install with PyPIpip install demcoregDocumentationhttp://demcoreg.readthedocs.io(autogenerated from source code, may be out of date)LicenseThis project is licensed under the terms of the MIT License.CitationIf you use any of this software for research applications that result in publications, please cite:Shean, D. E., O. Alexandrov, Z. Moratto, B. E. Smith, I. R. Joughin, C. C. Porter, Morin, P. J., An automated, open-source pipeline for mass production of digital elevation models (DEMs) from very high-resolution commercial stereo satellite imagery, ISPRS J. Photogramm. Remote Sens, 116, 101-117, doi:10.1016/j.isprsjprs.2016.03.012, 2016.