Over the last decades, an important number of regional and global Digital Elevation Models (DEMs) have been released publicly. As a consequence, researchers need to choose between several of these models to perform their studies and to use these DEMs as third-party data to compute derived products (e.g., for orthorectification). However, the comparison of DEMs is not trivial. For most quantitative comparisons, DEMs require to be expressed in the same Coordinate Reference System (CRS), sampled over the same grid (i.e., be at the same ground sampling distance with the same pixel-is-area or pixel-is-point convention), and with heights relative to the same Vertical Reference System (VRS). Thankfully, many open tools allow us to perform these transformations precisely and easily. Despite these rigorous transformations, local or global planimetric displacements may still be observed from one DEM to another. These displacements, or disparities, may lead to significant biases in comparisons of DEM elevations or derived products such as slope, aspect or curvature. Therefore, the control of DEM planimetric accuracy is certainly the most important task to perform at first before any comparison. This paper presents the disparity analysis method enhanced to achieve a sub-pixel accuracy, by interpolating the linear regression coefficients computed within an exploration window. This new method is significantly faster than oversampling the input data because it uses the correlation coefficients that have been already computed in the disparity analysis. To demonstrate the robustness of this algorithm, artificial displacements have been performed in a 11x11 grid with a 0.1-pixel step in both directions. The Copernicus DEM GLO-30 has been interpolated by a “classical” bicubic to achieve these 11x11 displacements. This validation method has been applied in four 10 km x 10 km DEMIX tiles showing different height distributions. Globally, this new sub-pixel accuracy method is robust. Artificial displacements have been retrieved with with typical errors (eb) ranging from 12 to 20 % of the pixel size (worst case in Croatia).These errors in displacement retrievals are not constant in the 11x11 grid and the overall error Eb depends on the roughness (altimetry variations) encountered in the different tiles. The second idea of this paper is to assess the impact of the bicubic parameter (slope of the weight function at a distance d=1 of the interpolated point) on the accuracy of the displacement retrieval. By considering Eb as a quality indicator, tests have been performed in the four DEMIX tiles making the bicubic parameter vary between -1.5 and 0.0 by step of 0.1. For each DEMIX tile, the best bicubic (BBC) parameter b* is interpolated from the four Eb minimal values. This BBC parameter b* is low for flat areas (around -0.95) and higher in mountainous areas (around -0.75). The roughness indicator is the standard deviation of the slope norms computed from all the pixels of a tile. A logarithmic regression analysis performed between the roughness indicator and the BBC parameter b* computed in 67 DEMIX tiles shows a high correlation (r=0.717). The logarithmic regression formula estimating the BBC parameter from the roughness indicator is generic and may be applied to estimate the displacements between two different DEMs. This formula may also be used to set up a future Adaptative Best BiCubic (ABBC) that will estimate the local roughness in a sliding window to compute a local BBC.
