Co-registration and residual correction of digital elevation models: a comparative study

Abstract

Digital elevation models (DEMs) are currently one of the most widely used data sources in glacier thickness change research, due to the high spatial resolution and continuous coverage. However, raw DEM data are often misaligned with each other, due to georeferencing errors, and a co-registration procedure is required before DEM differencing. In this paper, we present a comparative analysis of the two classical co-registration methods proposed by Nuth and Kääb (2011) and Rosenholm and Torlegard (1988). The former is currently the most commonly used method in glacial studies, while the latter is a seminal work in the photogrammetric field that has not been extensively investigated by the cryosphere community. Furthermore, we also present a new residual correction method using a generalized additive model (GAM) to eliminate the remaining systematic errors in DEM co-registration results. The performance of the two DEM co-registration methods and three residual correction algorithms (the GAM-based method together with two parametric-model-based methods) was evaluated using multiple DEM pairs from the Greenland Ice Sheet and mountain glaciers, including Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) DEMs, ZiYuan-3 (ZY-3) DEMs, the Shuttle Radar Topography Mission (SRTM) DEM, and the Copernicus DEM. The experimental results confirm our theoretical analysis of the two co-registration methods. The method of Rosenholm and Torlegard has a greater ability to remove DEM misalignments (an average of 4.6 % and 13.7 % for the test datasets from Greenland Ice Sheet and High Mountain Asia, respectively) because it models the translation, scale, and rotation-induced biases, while the method of Nuth and Kääb considers translation only. The proposed GAM-based method performs statistically better than the two residual correction methods based on parametric regression models (high-order polynomials and the sum of the sinusoidal functions). A visual inspection reveals that the GAM-based method, as a non-parametric regression technique, can capture complex systematic errors in the DEM co-registration residuals.