# Investigating the Impact of Digital Elevation Models on Sentinel-1 Backscatter and Coherence Observations

Ignacio Borlaf-Mena, Maurizio Santoro, Ludovic Villard, Ovidiu Badea, Mihai Tanase

► To cite this version:

Ignacio Borlaf-Mena, Maurizio Santoro, Ludovic Villard, Ovidiu Badea, Mihai Tanase. Investigating the Impact of Digital Elevation Models on Sentinel-1 Backscatter and Coherence Observations. Remote Sensing, 2020, 12 (18), pp.3016. 10.3390/rs12183016. hal-04928749

HAL Id: hal-04928749

https://hal.science/hal-04928749v1

Submitted on 4 Feb 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Article

# Investigating the Impact of Digital Elevation Models on Sentinel-1 Backscatter and Coherence Observations

Ignacio Borlaf-Mena <fnref n="1" />,<fnref n="2" />,* ¹D, Maurizio Santoro <fnref n="3" /> ¹D, Ludovic Villard <fnref n="4" />, Ovidiu Badea <fnref n="1" />,<fnref n="5" /> and Mihai Andrei Tanase <fnref n="1" />,<fnref n="2" />

Received: 27 July 2020; Accepted: 10 September 2020; Published: 16 September 2020

<figure><img src="image_3.png" /></figure>

check for
updates

**Abstract:** Spaceborne remote sensing can track ecosystems changes thanks to continuous and systematic coverage at short revisit intervals. Active remote sensing from synthetic aperture radar (SAR) sensors allows day and night imaging as they are not affected by cloud cover and solar illumination and can capture unique information about its targets. However, SAR observations are affected by the coupled effect of viewing geometry and terrain topography. The study aims to assess the impact of global digital elevation models (DEMs) on the normalization of Sentinel-1 backscattered intensity and interferometric coherence. For each DEM, we analyzed the difference between orbit tracks, the difference with results obtained with a high-resolution local DEM, and the impact on land cover classification. Tests were carried out at two sites located in mountainous regions in Romania and Spain using the SRTM (Shuttle Radar Topography Mission, 30 m), AW3D (ALOS (Advanced Land Observation Satellite) World 3D, 30 m), TanDEM-X (12.5, 30, 90 m), and Spain national ALS (aerial laser scanning) based DEM (5 m resolution). The TanDEM-X DEM was the global DEM most suitable for topographic normalization, since it provided the smallest differences between orbital tracks, up to 3.5 dB smaller than with other DEMs for peak landform, and 1.4–1.9 dB for pit and valley landforms.

**Keywords:** synthetic aperture radar (SAR); radiometric terrain normalization; digital elevation model (DEM); coherence; backscatter; Sentinel-1; LiDAR; land cover classification

## 1. Introduction

Synthetic aperture radar (SAR) is an active imaging system with several advantages over optic sensors, such as Landsat OLI (Operational Land Imager) or Sentinel-2 MSI (Multi-Spectral Imager). SARs are independent of solar illumination and use wavelengths that can penetrate cloud cover and have unique interactions with ground targets. Furthermore, the capability to transmit and receive signals enables the use of both phase and polarization information, to monitor, among others, landslides, avalanches, snowmelt, and forests [1]. Terrain orientation affects the intensity of the backscattered signal based on lambert cosine law. The signal is further affected by the terrain scattering area. The SAR technique uses the return time to convert a table of recorded echoes into an image

<footnote n="1">¹ Romanian National Institute for Research and Development in Forestry, INCDS "Marin Drăcea", Department of Forest Monitoring, Bulevardul Eroilor 128, 077190 Voluntari, Romania; obadea@icas.ro (O.B.); mihai.tanase@uah.es (M.A.T.)</footnote>

<footnote n="2">² Department of Geology, Universidad de Alcalá, Geography and Environment, Calle Colegios 2, 28801 Alcalá de Henares, Spain</footnote>

<footnote n="3">³ Gamma Remote Sensing, Worbstrasse 225, 3073 Gümligen, Switzerland; santoro@gamma-rs.ch</footnote>

<footnote n="4">⁴ Centre d'Etudes Spatiales de la Biosphère, 31400 Toulouse, France; ludovic.villard@cesbio.cnes.fr</footnote>

<footnote n="5">5 Faculty of Silviculture and Forest Engineering, Department of Forest Engineering, "Transilvania" University, Forest Management Planning and Terrestrial Measurements, Ludwig van Beethoven Str. 1, 500123 Braşov, Romania</footnote>

<footnote>* Correspondence: ignacio.borlaf@edu.uah.es</footnote>
(focusing). These times are shortened in sensor-facing steep slopes, causing echoes to overlap and slopes to appear shortened in the focused image (a pixel represents more area). Both effects can be compensated for by using the acquisition geometry parameters and a digital elevation model (DEM). We refer to this compensation as topographic normalization. The backscattered intensity is normalized by accounting for the scattering area and the local incidence angle [2,3]. The coherence normalization is based on removing the topographic phase component from the interferogram before estimation [4]. Normalization results are heavily dependent on the DEM characteristics and quality [2,3] as they can be generated using different data sources, including remote sensing (optic, SAR, airborne laser scanning—ALS), and processing techniques (ALS point cloud, stereography, interferometry, and radargrammetry)

ALS uses the delay between emission and reception of light pulses to determine the 3D position of objects. When the pulses return, the instrument registers an intensity profile, which can be completely (full-waveform) or partially recorded (discrete return, i.e., the position of the leading edge before the peak) [5]. Airborne discrete return systems have become the source of national elevation datasets for many countries [6,7]. Spaceborne full-waveform Light Detection and Ranging (space LiDAR, such as the ICESat (Ice, Cloud, and land Elevation Satellite)) data have also been employed in the context of topographic mapping as a primary source for calibration or validation of global elevation datasets derived from other sensors [8,9]. However, global topographic mapping from space based on LiDAR is difficult, due to the sparse sampling, and the sensitivity to cloud cover.

Stereoscopic techniques are based on differences in the line of sight to objects (parallax) for common points (tie-points) in an overlapping set of images. Results are dependent on tie-point quantity, image contrast, noise, and features (such as shadows and homogeneous surfaces), which may pose problems [10–12]. These techniques have been applied over aerial [6,11] and satellite optical images [9,13], as well as SAR imagery (i.e., radargrammetry) [10,14]. Stereoscopic processing of optical images was used to generate the ASTER global DEM (ASTER GDEM) and the ALOS (Advanced Land Observation Satellite) World 3D Digital Surface Model (ALOS AW3D DSM). The ASTER GDEM was created using imagery from the Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) onboard the Terra satellite. ASTER stereo pairs were formed from two near-infrared images (nadir, backward) with 15m resolution. The AW3D DEM was based on data from the Panchromatic Remote-sensing Instrument for Stereo Mapping (PRISM) onboard the ALOS. PRISM stereo acquisitions were formed with three panchromatic images (forward, nadir, backward) with 2.5 m resolution. In both cases, the most challenging task was masking clouds, snow, ice, or water on every acquisition [9] to avoid the introduction of outliers. Afterwards, images from each individual acquisition were matched, and elevation was calculated. All height estimates from individual acquisitions were stacked to ensure continuity and reduce noise [9,13]. The AW3D was corrected for biases using preexisting data, such as ICESat shots and the preexisting Shuttle Radar Topography Mission (SRTM) DEM (see next paragraphs) [9].

Even though the processing chains for ASTER GDEM and AW3D were similar, their accuracies are different. Studies comparing global DEMs based on ground control points (GCP) report that the ASTER GDEM has larger uncertainties and is affected by striping, hummock-like artifacts, and outliers [12,15–17]. These artifacts may stem from the tie point generation, the choice of band (NIR), its relatively low spatial resolution [16,17], or unremoved cloud patterns [15]. AW3D performed better, although hillslope and step-like artifacts (scene mismatch) have been found [12,16–19].

SAR interferometry uses two co-registered SAR images acquired from close orbits. The interferogram (i.e., the phase difference between the two SAR images) relates to the 3D position of each target on the ground. Thereof, an interferogram reproduces the topographic information, which appears in the form of fringes as phase is measured between 0 and $2\pi$. To obtain absolute phase values, from which elevation can be estimated, the interferogram is unwrapped. Unwrapping may be hindered in areas of steep topography or areas affected by the lack of coherence between images in
consequence of changes between acquisitions (wind-induced motion, precipitation, etc.) [1,14]. SAR interferometry was used to generate two global DEMs, the SRTM DEM, and the TanDEM-X DEM.

The SRTM acquired data over 80% of the Earth's land surface (60° N–56° S) on an 11-day orbital flight in February 2000. SRTM operated two antennas physically separated in space by 60 m at C-band, as well as X-band. At C-band, a gap-free coverage was obtained with single-pass interferograms. The interferometric height was reconstructed (unwrapping) and re-gridded into map coordinates with variable-resolution smoothing. Data takes were combined using coincident tie points [20]. The main artifacts of the SRTM dataset were related to striping from uncompensated movements of the mast, voids in correspondence of steep slopes or for low coherence areas, or coarser than nominal spatial detail from the re-gridding step [17,20–22]. The global TanDEM-X DEM was generated using SAR data acquired during 2010–2015 by the TanDEM-X and TerraSAR-X satellites flying in formation. Individual scenes were focused, multi-looked (sample averaging) to 10–12 m pixel spacing and unwrapped. The DEM was generated in an iterative process with the first global coverage using data acquisition parameters (baseline) adequate for moderate terrain. The second global coverage was shifted half a swath, and its unwrapping was aided by the data from the first coverage. Over some areas, further coverages were acquired from a different viewing geometry to avoid errors caused by topographic distortions or volumetric scattering (e.g., forest, desert) [23,24]. Areas with height ambiguities were infilled using radargrammetric processing of the scenes. TanDEM-X DEM data are distributed with 12.5 m (0.4 arcsec, original), 30 m (1 arcsec), and 90 m resolution (3 arcsec). The latter was generated by the unweighted average of the overlapping 12.5 m pixels [17].

All DEMs are affected by contributions from elements covering the terrain, such as cities or vegetation, thus reporting elevations higher than those recorded for the ground control points [15,18]. For this reason, elevation refers to surface elevation rather than terrain elevation. For the specific case of vegetation, the main reason for this is the different penetration of each wavelength. On the one hand, the nanometric-scale wavelengths employed for generating photogrammetric DEMs have limited penetration, and thus the tie points and the generated surface tend to reflect canopy surface height. On the other hand, the centimetric wavelengths employed by SAR sensors are able to penetrate further, albeit the scattering center height depends on the frequency employed and the vegetation structure [20,25–27]. Furthermore, the quality of InSAR-based DEM also depends on the spatial and temporal baselines, or unaccounted sensor movements.

Despite the rather large range of DEMs available, their effect on terrain normalization is poorly understood. Hoekman [28] suggested that the SRTM 90 m spacing may not be adequate for the radiometric normalization of the SAR backscatter in complex terrains. Recently Truckenbrodt [29] compared and tested the SRTM (30 and 90 m), AW3D (30 m), and TanDEM-X (90 m) DEMs in the context of radiometric terrain normalization of Sentinel-1 data. They analyzed the deviation of each DEM from the pixel-wise median of all DEMs and performed a regression analysis between terrain flattened $\gamma^0$ and the local incidence. The deviation analysis showed that the SRTM DEMs have the smallest difference from the median values, but the 30 m SRTM version had high deviation artifacts at one test site. The same errors were found for SRTM and AW3D, as the latter has been infilled with data from the former due persistent cloud cover. Both AW3D and TanDEM-X DEM contained outliers or noise over water areas. The 90 m TanDEM-X DEM was found to contain several large artifacts in mountainous areas. Regression analysis showed that all the analyzed DEMs largely removed the terrain influence, with complete removal (i.e., slope of 0) being observed in some experiments using higher resolution DEMs (SRTM 1 arcsec, AW3D) [29].

The objective of this study was to investigate the impact of global DEMs on the normalization of SAR backscatter and coherence observations by Sentinel-1 at two sites characterized by complex topography. To build on previous literature [29], we analyzed the performance of 12.5 (resampled to 20 m, see Section 3), 30, and 90 m pixel Tandem-X DEMs, along with AW3D (30 m), SRTM (30 m). A very high-resolution (5 m) ALS DEM was used to benchmark results. The impact of terrain normalization was assessed by investigating the inter-orbit variability of the observations by land
cover and landforms. Then, we evaluated a land cover classification scheme based on the observations
normalized for topography.

## 2. Study Area and Satellite Data

The study area consisted of an N-S transect over the Romanian Carpathians (11,700 km²) and the National Park of Sierra Nevada in Spain (2360 km²) (Figure 1). We could not add additional sites, due to TanDEM-X scientific proposal area limitations [30]. The sites were selected to account for the different vegetation types and structures encountered in the temperate and Mediterranean climates. Due to the more humid climatic conditions, the vegetation in the Carpathians is characterized by denser, taller, and more diverse forest types (broadleaf, needleleaf, mixed) when compared to the sparser and shorter forests dominated by pine species encountered in the Sierra Nevada.

<figure><img src="image_4.png" /><figcaption>**Figure 1.** Extent of the study areas (a, Romania, Ro; c, Spain, Sp) and digital elevation models (DEMs) used for synthetic aperture radar (SAR) data processing. The yellow box indicates the location of the area covered by the Sentinel-1 dataset used for the analysis. The red box indicates the extent of the subset shown in the right-hand side panels (b, Romania; d, Spain).</figcaption></figure>

Histograms of natural land covers for each site have been plotted to show their frequencies relative to the slopes they occupy (Figure 2). Needleleaf and mixed forests occupy steep slopes, whereas grasslands and broadleaf forests occupy moderate slopes at the Romanian site. At the Spanish site, bare and needleleaf forests occupy moderate slopes, although for the former, a significant fraction of the pixels occupy steep or very steep slopes.
<figure><img src="image_5.png" /><figcaption>**Figure 2.** Distribution of the slopes for grassland, bare soils, and forest land. The latter has been separated according to leaf type.</figcaption></figure>

For each site, we assembled the SRTM 1-arcsecond, i.e., 30 m, DEM [31] from the United States Geological Service (USGS) Earth Explorer [32], the AW3D 30 m DEM from the Japan Aerospace Exploration Agency (JAXA) Earth Observation Research Center (EORC) [33] and the TanDEM-X DEM (©DLR 2019) with a pixel spacing of 12.5 m (original resolution), 30 m and 90 m (resampled) from the German Aerospace Agency [34]. In addition, for the Spanish site, we used an ALS-based DEM to benchmark the results obtained from three global DEMs. The ALS DEM was available through the Spanish national plan of orthophotography (PNOA) from the National Center of Geographic Information of Spain (CNIG) [35,36]. The ALS DEM was created from ALS point clouds with a density of 0.5 returns/m². The ALS scan over our study area was performed in 2014 with the LEICA ALS60 sensor. The points were translated from ellipsoidal to ortho-metric heights, assigned color (RGB and near-infrared) from PNOA orthophotographs, and classified automatically using TerraScan [37,38]. Classification eliminates returns considered noise and filters point to avoid oversampling due to flight strip overlap. Ground points classification is based on slope, rugosity, and return count. Vegetation and Buildings are classified based on height, separating both based on NDVI values. The DEM is generated by calculating the mean value of all ground returns within a 5m pixel [39]. The reported accuracies of the DEM products are presented in Table 1.

<table><caption>**Table 1.** Reported accuracies of the DEM used in this study.</caption><thead><tr><th rowspan="2">DEM Product</th><th rowspan="2">Pixel Spacing</th><th colspan="2">Accuracy</th><th rowspan="2">Relative Vertical Accuracy</th><th rowspan="2">Coverage</th><th rowspan="2">Reference</th></tr><tr><th>Horizontal</th><th>Vertical</th></tr></thead><tbody><tr><td>SRTM DEM</td><td>~30 m</td><td>≤12.6 m</td><td>≤9 m</td><td>≤9.8 m</td><td>Nearly global (60° N–56° S)</td><td>[40]</td></tr><tr><td>AW3D DSM</td><td>~30 m</td><td>-</td><td>&lt;7 m</td><td>&gt;3 m (slope ≤ 20%)<br/>&gt;5 m (slope &gt; 20%)</td><td>Global</td><td>[41]</td></tr><tr><td>TanDEM-X DEM</td><td>~12.5 m</td><td>&lt;10 m</td><td>&lt;10 m</td><td>2 m (slope ≤ 20%)<br/>4 m (slope &gt; 20%)</td><td>Global</td><td>[23]</td></tr><tr><td>PNOA</td><td>~5 m</td><td>≤0.5 m</td><td>≤0.5 m</td><td>-</td><td>Spain</td><td>[42]</td></tr><tr><td>LiDAR DEM</td><td></td><td></td><td></td><td></td><td></td><td></td></tr></tbody></table>

The SAR dataset consisted of a time series of Sentinel-1 dual-polarized (VV and VH) images acquired in the Interferometric Wide Swath (IWS) mode. The images were obtained in Single Look Complex format (SLC) with a pixel spacing of 14.1 m in azimuth and 2.3 m in range. All SAR images were resampled to a pixel size matching the different DEMs used (see Section 3.2). For the Romanian site, 21 images acquired between 30 December 2016 and 6 February 2017 from three relative orbits (7, 29, 131) were used. For the Spanish site, 18 images acquired between 21 August 2018 and 8 October 2018 from relative orbits 1 and 81 were used.

### 3. Methods

The following analyses and processes where carried away using the software GDAL/OGR [43], GAMMA software [44], GRASS (Geographic Resources Analysis Support System) [45], Python [46], Rasterio [47], Pandas [48], Geopandas [49], Numpy [50], Scipy [51], and Matplotlib [52].
## 3.1. DEM Assembly

The global DEMs were provided in equiangular geographic coordinates. The SRTM and the AW3D DEMs height reference had to be shifted from geoidal to ellipsoidal heights without resampling. The TanDEM-X (TDX) DEMs were provided as height above the ellipsoid. The Tandem-X DEM at 30 m (TDX30) was used as provided. The Tandem-X 12.5 m DEM was resampled (bilinear interpolation) to 20 m pixel spacing (TDX20) to reduce pixel size difference with respect to the multi-looked Sentinel-1 image. The 90 m resolution Tandem-X DEM (TDX90) was resampled to 30 m (bilinear interpolation). The ALS DEM, originally projected to ETRS89 UTM zone 30N coordinate system, was translated to a height above the ellipsoid and resampled to a 20 m pixel size (bilinear interpolation). As the ALS DEM has not been re-projected, all products geocoded with it (i.e., geocoded Sentinel-1 backscatter) share the same projection.

## 3.2. SAR Data Preparation

For each SAR image, the SLC sub-swathes were mosaicked, and the resulting image was multi-looked by a factor of 7 in range and 2 in azimuth. The objective was to reduce noise and obtain the SAR backscattered intensity at a pixel spacing close to the target 20 m used for the analysis. For a given orbit, the first acquired image was used as a master. All remaining SAR images from the same orbit were co-registered to the master image using an iterative process based on intensity matching and spectral diversity aided by each DEM [53]. For each orbit, the master image was used to generate a lookup table (LUT) relating map and range doppler coordinates. The LUT was used to orthorectify the master and the co-registered images (interferograms and SAR backscatter) from the same orbit.

Interferograms were generated for each consecutive image pair (a-b, b-c, c-d, etc.), and the DEM-estimated topographic phase was subtracted from each. The interferometric coherence was estimated in a two-step adaptive approach [54,55]. The first estimate of coherence was obtained with a 3-by-3 window. To reduce the estimation bias due to the small window size [56], the coherence was then recomputed using a window size inversely proportional to the initial estimate of the coherence. As a trade-off between preserving spatial resolution and reducing the bias, the largest window size was set to 9-by-9 pixels. In addition, when the estimation window included scatterers with a coherence level different than the coherence of the target in the center of the window, the estimator masked out such features to preserve the true coherence of the latter target.

The backscatter coefficient was calibrated to terrain flattened $\gamma^0$, considering the scattering area on the ellipsoid and on DEM surfaces ($A_{flat}$ and $A_{slope}$) [3] and the incidence angle on the ellipsoid and on DEM surfaces ($\theta_{ref}$ and $\theta_{loc}$), as reported in equation 1 [57]. The parameter $n$ can be employed to account for volume effects [57]. The parameter was set to 1, an adequate value for most land cover types, as dealing with volumetric effects was not the objective of this study. The backscatter intensity and coherence images were orthorectified using the LUT and an inverse distance resampling. LUT coordinates located more than two pixels (range) apart to its counterpart on the SAR image were masked as no data. To reduce speckle, the multi-temporal backscatter images were averaged, by polarization, in time. Seven images were averaged for each orbit for the Romanian site, while nine images were averaged for each orbit for the Spanish site. Similarly, the six and eight coherence images were averaged in time for each orbit and polarization for the Romanian and the Spanish site, respectively.

$$ \gamma^0 = \sigma^0 \frac{A_{flat}}{A_{slope}} \left( \frac{\cos \theta_{ref}}{\cos \theta_{loc}} \right)^n \quad (1) $$

## 3.3. Auxiliary Datasets

In support of the analysis, a land cover dataset was created for each site, based on the agreement between the ESA CCI land cover map (2015) [58], the DLR's global urban footprint (GUF) (2016) [59–61], the ALOS PALSAR forest map (ALOS FNF) (2017) [62], and either Corine land cover map 2012
(CLC) [63] (Romanian Site) or Spanish information system on soil occupation (SIOSE, 2014) [64] (Spanish site, more detailed and recent) (Table 2). ESA CCI land cover maps are generated at 300 m resolution using optical imagery time series from AVHRR, MERIS, SPOT-VGT, and POBA-V imagery, as well as GlobCover unsupervised classification chain and machine learning [58]. GUF is generated from TanDEM-X imagery at 12 m resolution based on amplitude and texture [59–61]. ALOS forest map is generated based on local thresholding of annual composites of ALOS PALSAR 1/2 amplitude images (25 m pixel side) [62]. Both CLC and SIOSE are based on photointerpretation of satellite and aerial imagery, with minimum polygon surface of 25 Ha [63] and 1 Ha [65], respectively.

<table><caption>**Table 2.** Composition of the analyzed land covers based on preexisting datasets. When a higher level of Corine Land Cover or SIOSE has been employed (CLC (Corine land cover map) Lvl.1), the rest have been filled as "x". CCI LC forest types are further disaggregated by the fractional cover and were therefore aggregated. GUF, global urban footprint; ALOS FNF, Advanced Land Observation Satellite forest map; SIOSE, Spanish information system on soil occupation; CODIIGE, Board of directors of the geographic information infrastructure of Spain.</caption><thead><tr><th></th><th>CLC 2012 [63]</th><th>SIOSE 2014 [64] (CODIIGE)</th><th>CCI LC 2015 [66]</th><th>GUF 2016 [59–61]</th><th>ALOS FNF 2017 [62]</th></tr></thead><tbody><tr><td>Urban</td><td>1xx: Artificial surfaces</td><td>1xx: Artificial surfaces</td><td>190: Urban areas</td><td>Urban</td><td>-</td></tr><tr><td>Crops</td><td>21x: Arable land</td><td>210: Crops (herbaceous)</td><td>10, 20: Cropland, 11: Herbaceous cover</td><td>Other</td><td>Other</td></tr><tr><td>Pasture</td><td>23x: Pastures</td><td>320: Pastures</td><td>11: Herbaceous 130: Grassland</td><td>Other</td><td>Other</td></tr><tr><td>Grassland</td><td>321: Grassland</td><td>320: Pastures</td><td>11: Herbaceous 130: Grassland</td><td>Other</td><td>Other</td></tr><tr><td>Bare</td><td>33x: Open spaces</td><td>354: Bare</td><td>200: Bare areas</td><td>Other</td><td>Other</td></tr><tr><td>Broadleaf forest</td><td>311: Forest (broadleaf)</td><td>311: Forest (broadleaf)</td><td>50–62: Tree cover, broadleaved</td><td>Other</td><td>Forest</td></tr><tr><td>Needleleaf forest</td><td>312: Forest (needleleaf)</td><td>312: Forest (needleleaf)</td><td>70–82: Tree cover, needle leaved</td><td>Other</td><td>Forest</td></tr><tr><td>Mixed forest</td><td>313: Forest (mixed)</td><td>313: Forest (mixed)</td><td>90: Tree cover, mixed leaf type</td><td>Other</td><td>Forest</td></tr><tr><td>Water</td><td>5xx: Water bodies</td><td>5xx: Water bodies</td><td>210: Water</td><td>Other</td><td>Other</td></tr></tbody></table>

The ALOS FNF disagreed with the remaining data sets as part of the cities were classified as forest, and part of the water was classified as "other" (not water, nor forest). For cities, no ALOS FNF condition was applied, whereas the rest of the non-forest classes on other datasets were considered compatible with non-forest classes from ALOS FNF (non-forest, water). The polygons with the agreement were dissolved by the land cover to eliminate internal borders, and a negative buffer of 40 m was applied to avoid edge effects.

The analyses were undertaken in this study are also related to landforms (Figure 3), i.e., features of the terrain surface with a distinct and identifiable shape [67]. Landforms were labeled using the GRASS GIS add-on "r geomorphon" [68], with a search window of 25 pixels and a "flatness" threshold of 5 degrees applied to the highest spatial resolution DEMs available for each site, i.e., the TanDEM-X DEM at 30 m for the Romanian site and the PNOA DEM aggregated to 30 m for the Spanish site.
<figure><img src="image_6.png" /><figcaption>**Figure 3.** Landform classification. The general shape of the landforms (a), modified from GRASS documentation, based on [68]) and landforms over the shaded relief for a subset of the Spanish study area (b).</figcaption></figure>

### 3.4. Inter-Orbital Data Analysis

The topographic normalization (radiometric terrain normalization, topographic phase removal) and distortion masking (e.g., foreshortening, layover, shadows) of each SAR-derived variable (backscatter coefficients and coherence) were assessed using the inter-orbit range (IOR). IOR was calculated pixelwise for each SAR variable by subtracting the maximum and the minimum values available from all orbits. For example, for a pixel with data available from orbits a, b, and c, the IOR would be max(a, b, c) - min(a, b, c). The inter-orbit range was plotted by land cover class (boxplots). The analysis was repeated by landforms for the needleleaf forest, the only common forest type between both sites, and classes appearing near the mountain tops (grassland, bare). Because the only difference in image processing is the DEM employed, low IOR values reflect improved topographic effects removal.

In addition, the scattering area estimates, backscatter coefficient, and the interferometric coherence obtained using the ALS DEM (the most detailed) were employed as a reference to assess the performance of the global DEMs at the Spanish site. Said assessment was based on the root mean square deviation (RMSD, Equation (2)), relative RMSD (RMSDrel, Equation (3)), mean absolute deviation (MAD, Equation (4)), and Offset (Equation (5)) for the products obtained from the global DEMs.

$$RMSD = \sqrt{\frac{1}{P} \sum_{p=1}^{P} (v_p - r_p)^2} \qquad (2)$$

$$RMSDrel = \frac{RMSD}{\bar{r}} \qquad (3)$$

$$MAD = \frac{1}{P} \sum_{p=1}^{P} |v_p - r_p| \qquad (4)$$

$$Offset = \bar{v} - \bar{r} \qquad (5)$$
where *P* is the total number of pixels, *p* is a specific pixel and *v*<sub>n</sub> and *r*<sub>n</sub> are the variable values (i.e., backscatter) obtained for said pixel using a global DEM (*v*) and the reference DEM (*r*), whereas $\bar{v}$ and $\bar{r}$ are the mean value for said variables.

### 3.5. Land Cover Classification

A Linear support vector machine (LinearSVM) classifier was selected for its robustness and short execution times. We employed the Scikit-Learn implementation [69] with default options with a regularization parameter of 1, primal problem optimization, 0.001 tolerance for stopping criteria, and 10,000 iterations maximum. The classifier was trained per orbit/DEM pair using 96,000 samples, using VV- and VH-polarized backscatter and co-pol coherence as features.

For each orbit and land cover class, 70% of the valid sample (foreshortened and shadowed pixels were masked during SAR processing) was used to calculate the median and the median absolute deviation (MAD) of each SAR variable. Median and MAD values were then employed to calculate the z-score for each predictor (VV- and VH-polarized backscatter and coherence) by land cover class. Only samples with an absolute z-score below three were retained. Depending on the land cover class, the number of pixels retained varied from several millions (forest and low vegetation) down to tens of thousands (urban and water). For each class, 12,000 pixels (the number of pixels available for the less extended class, i.e., water) were randomly selected and used for training (n). As low vegetation and forest classes were further split into three sub-classes each (crops, pastures, grasslands, and broadleaf, needleleaf, and mixed forests), the total numbers of training samples selected were thrice as much (3n) as for urban and water land cover classes.

The validation sample, formed by the remaining pixels (30%) of each class, was used to compute the confusion matrix and associated error metrics (i.e., user and producer accuracies, Cohen's Kappa, as described below). Error metrics for valleys were calculated after resampling (nearest neighbor) the landform layer to match the spatial resolution of the DEMs.

The confusion matrix C (Equation (6)) represents the occurrences of the predicted (rows) against the actual land cover class (columns) (*r*·*r* dimensions, where *r* is the number of classes). Diagonal cells (*c*<sub>ii</sub>) count pixels with the same class in the classification and the reference dataset (True Positive, TP). Cells over the diagonal count pixels of class *i* that have received other class (False Negative, FN), whereas cells under the diagonal count pixels that have been classified as *i*, when they have other class in the reference dataset (False Positive, FP). Accuracies for a specific class *i* are the count of correctly classified pixels for the class (*c*<sub>ii</sub>) divided by the number of pixels classified as *i* (count across columns, *c*<sub>i+</sub>), in the case of user accuracy (*UA*<sub>i</sub>, also called 'precision' in machine learning literature, Formula (7)), or by the number actual *i* pixels (count across columns, *c*<sub>+i</sub>), in the case of producer accuracy (*PA*<sub>i</sub>, also called 'recall' in machine learning literature, Formula (8)).

$$C = \begin{bmatrix} c_{1,1} & c_{1,j} & c_{1,r} \\ c_{i,1} & c_{i,j} & c_{i,r} \\ c_{r,1} & c_{r,j} & c_{r,r} \end{bmatrix} \qquad (6)$$

$$UA_i = precision_i = \sum_{i=1}^{r} \frac{c_{ii}}{c_{i+}} = \frac{TP}{TP + FP} \qquad (7)$$

$$PA_i = recall_i = \sum_{i=1}^{r} \frac{c_{ii}}{c_{+i}} = \frac{TP}{TP + FN} \qquad (8)$$

Cohen's Kappa [70] is a measure of agreement between the predicted cover and the one appearing on the reference dataset.

$$\hat{K} = \frac{N \sum_{i=1}^{r} c_{ii} - \sum_{i=1}^{r} (c_{i+} \cdot c_{+i})}{N^2 - \sum_{i=1}^{r} (c_{i+} \cdot c_{+i})} \qquad (9)$$
where *r* is the number of rows, *c*<sub>ii</sub> is the number of pixels where there is agreement between the classification and the reference dataset (cells on the diagonal, with row *i* and column *i*), *c*<sub>*i*+</sub> and *c*<sub>*i*+</sub> are the totals for row *i* (count of pixels classified as *i*) and column *i* (count of reference pixels with class *i*), and finally, *N* is the total number of observations [71]. A Kappa value of 1 represents the complete agreement between both, 0 represents a classifier performance similar to random guessing, and values under 0 indicate results worse than random guessing.

# 4. Results

In Sections 4.1 and 4.2, we show results of the IOR analysis by land cover and land form, respectively, Section 4.3 contains the ALS reference-based analysis, and Section 4.4 describes the results of the classification comparison. Through Sections 4.1–4.3, coherence showed very small differences based on the DEM employed, under 0.01 for the IOR based analyses, and under 0.02 for ALS reference-based analysis. For this reason, these tables have been omitted, as they carried little to no information.

## 4.1. Inter-Orbital Range by Land Cover

The inter-orbit range (IOR, Table 3) was analyzed as an indicator of the residual terrain influence on the normalized SAR metrics (higher IOR, higher influence). The mean IOR for urban cover varied very little between DEMs, except for TDX90. Crops and broadleaf forests presented little difference (up to 0.3 dB) depending on the DEM employed, but increased for mixed forests (0.5 dB, Carpathians), and classes appearing near mountain peaks, such as grassland (0.6 dB) and bare soil (1 dB). Needleleaf forests had similarly high differences at the Romanian site (0.7 dB), whereas they were lower at the Spanish site (0.4 dB).

<table><caption>**Table 3.** Backscatter Inter-orbit range (IOR) by polarization and land cover class at each study site (GL, grassland; BLF, NLF, and MLF, are broadleaf, needleleaf, and mixed forest). Cell color shows the gradient between the lowest (green) and the highest value (yellow). “M.D.” column represents the maximum difference between global DEMs for each specific land cover.</caption><thead><tr><th rowspan="3"></th><th colspan="5">Romania</th><th colspan="5">Spain</th><th rowspan="3">ALS</th><th rowspan="3">M.<br>D.</th></tr><tr><th>AW</th><th>SR</th><th>TDX</th><th>TDX</th><th>TDX</th><th rowspan="2">M.<br>D.</th><th>AW</th><th>SR</th><th>TDX</th><th>TDX</th><th>TDX</th></tr><tr><th>3D</th><th>TM</th><th>20</th><th>30</th><th>90</th><th>3D</th><th>TM</th><th>20</th><th>30</th><th>90</th></tr></thead><tbody><tr><td>Urban</td><td>5.4</td><td>5.4</td><td>5.3</td><td>5.4</td><td>6.9</td><td>1.6</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Crops</td><td>2.3</td><td>2.3</td><td>2.2</td><td>2.2</td><td>2.3</td><td>0.1</td><td>1.3</td><td>1.4</td><td>1.3</td><td>1.3</td><td>1.5</td><td>1.2</td><td>0.2</td></tr><tr><td>GL</td><td>3.4</td><td>3.4</td><td>2.9</td><td>3.1</td><td>3.5</td><td>0.6</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Bare</td><td></td><td></td><td></td><td></td><td></td><td></td><td>3.1</td><td>3.5</td><td>2.9</td><td>3.1</td><td>3.9</td><td>2.7</td><td>1.0</td></tr><tr><td>BLF</td><td>2.0</td><td>2.0</td><td>1.7</td><td>1.8</td><td>2.0</td><td>0.3</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>NLF</td><td>2.5</td><td>2.5</td><td>1.8</td><td>1.9</td><td>2.5</td><td>0.7</td><td>1.2</td><td>1.5</td><td>1.1</td><td>1.2</td><td>1.5</td><td>1.0</td><td>0.4</td></tr><tr><td>MLF</td><td>2.4</td><td>2.4</td><td>1.9</td><td>2.0</td><td>2.4</td><td>0.5</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr></tbody></table>

TDX20/30 was the global DEM with the lowest IOR, whereas the highest IOR values were observed for the TDX90 DEM, followed by SRTM DEM. Results with AW3D over natural covers (bare, grasslands, forest) varied depending on the site: At the Romanian site, its IOR was close to SRTM or TDX90 values; whereas at the Spanish site, it was closer to the values obtained with TDX20/30 DEMs. ALS-based PNOA DEM had the lowest IOR, with a large improvement over results with SRTM or TDX90 (up to 1.2 dB for bare areas, up to 0.5 for needleleaf forest), and a slight improvement over TDX20 results (0.2 dB for bare areas, and 0.1 dB for needleleaf forest). IOR values for PNOA and TDX20/30 DEMs showed less spread (Figure 4), concentrating around lower values, whereas the spread was larger for SRTM and TDX90. The main differences between sites were the lower spread of IOR at the Spanish site, and the behavior of AW3D IOR, which was closer to the TDX20/30 values at the Spanish site, but closer to SRTM/TDX90 values at the Romanian site.
<figure><img src="image_7.png" /><figcaption>**Figure 4.** Boxplot representing the VV backscatter IOR for grassland, bare areas, and needleleaf forests at both sites: Mean and median values (green triangle, orange line) and inter quantile-ranges (whiskers) for 5–95%.</figcaption></figure>

## 4.2. Inter-Orbital Ranges by Landform

Analyzing IOR values by land cover may have dampened potential differences as all pixels, regardless of landform, were averaged by class. A more detailed analysis was conducted by disaggregating mean IOR values by landforms for the grassland, bare land, and needleleaf forests (Table 4). These classes were selected as they often occur on steep slopes (Figure 4). Only landforms with more than 1000 pixels have been analyzed to limit spurious results, due to small sample size.

<table><caption>**Table 4.** IOR values disaggregated by landform for classes on mountain tops (grasslands, Romania; bare, Spain). Cell color shows the gradient between the lowest (green) and the highest value (yellow) for each row, which represents the IOR value for a specific landform when a certain was DEM employed. “M.D.” column represents the maximum difference between global DEMs for each specific landform.</caption><thead><tr><th rowspan="3"></th><th colspan="5">Romania (Grasslands)</th><th colspan="5">Spain (Bare)</th><th rowspan="2">ALS</th><th rowspan="2">M.<br>D.</th></tr><tr><th>AW</th><th>SR</th><th>TDX</th><th>TDX</th><th>TDX</th><th>M.</th><th>AW</th><th>SR</th><th>TDX</th><th>TDX</th><th>TDX</th></tr><tr><th>3D</th><th>TM</th><th>20</th><th>30</th><th>90</th><th>D.</th><th>3D</th><th>TM</th><th>20</th><th>30</th><th>90</th></tr></thead><tbody><tr><th>peak</th><td>3.1</td><td>3.0</td><td>2.5</td><td>2.6</td><td>3.5</td><td>1.0</td><td>2.9</td><td>3.3</td><td>2.6</td><td>2.7</td><td>6.1</td><td>2.2</td><td>3.5</td></tr><tr><th>ridge</th><td>3.0</td><td>2.9</td><td>2.5</td><td>2.7</td><td>3.0</td><td>0.5</td><td>2.7</td><td>3.1</td><td>2.5</td><td>2.7</td><td>3.6</td><td>2.3</td><td>1.1</td></tr><tr><th>spur</th><td>3.2</td><td>3.2</td><td>2.7</td><td>2.9</td><td>3.1</td><td>0.5</td><td>3.1</td><td>3.4</td><td>2.9</td><td>3.0</td><td>3.9</td><td>2.7</td><td>1.0</td></tr><tr><th>slope</th><td>3.6</td><td>3.6</td><td>3.0</td><td>3.3</td><td>3.6</td><td>0.6</td><td>3.3</td><td>3.7</td><td>3.1</td><td>3.3</td><td>3.8</td><td>2.9</td><td>0.7</td></tr><tr><th>hollow</th><td>4.2</td><td>4.3</td><td>3.6</td><td>3.9</td><td>4.6</td><td>1.0</td><td>3.2</td><td>3.8</td><td>3.1</td><td>3.2</td><td>3.9</td><td>2.8</td><td>0.8</td></tr><tr><th>valley</th><td>4.5</td><td>4.7</td><td>4.1</td><td>4.4</td><td>5.5</td><td>1.4</td><td>3.0</td><td>3.4</td><td>2.9</td><td>3.1</td><td>4.7</td><td>2.6</td><td>1.8</td></tr></tbody></table>

Both grasslands and bare areas presented large differences between DEMs. In the case of grasslands, the landforms valley, hollow, and peak showed the largest differences between DEMs (1.0–1.4 dB). In the case of bare areas, peak and valley showed the largest differences (3.5 dB and 1.8 dB, respectively), followed by ridge and spur (1.1 and 1.0 dB). For needleleaf forest, the largest differences between DEMs were observed for concave landforms (Table 5): Hollow (0.8–0.9 dB at either site), valleys (up to 0.8 dB at the Spanish site, and up to 1.4 dB at the Romanian site) and pits (up to 1.9 dB at the Romanian site).
<table><caption>**Table 5.** IOR values disaggregated by landform for needleleaf forests. Cell color shows the gradient between the lowest (green) and the highest value (yellow) for each row, which represents the IOR value for a specific landform when a certain was DEM employed. "M.D." column represents the maximum difference between global DEMs for each specific landform.</caption><thead><tr><th rowspan="3"></th><th colspan="6">Romania</th><th colspan="6">Spain</th></tr><tr><th colspan="2">AW</th><th colspan="2">SR</th><th colspan="2">TDX</th><th rowspan="2">M.</th><th colspan="2">AW</th><th colspan="2">SR</th><th colspan="2">TDX</th><th rowspan="2">M.</th></tr><tr><th>3D</th><th>TM</th><th>20</th><th>30</th><th>90</th><th>D.</th><th>3D</th><th>TM</th><th>20</th><th>30</th><th>90</th><th>D.</th></tr></thead><tbody><tr><td>peak</td><td>2.3</td><td>2.2</td><td>2.0</td><td>2.1</td><td>2.3</td><td>0.3</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>ridge</td><td>2.2</td><td>2.2</td><td>1.9</td><td>2.0</td><td>2.2</td><td>0.4</td><td>1.2</td><td>1.3</td><td>1.1</td><td>1.1</td><td>1.3</td><td>1.0</td><td>0.2</td></tr><tr><td>spur</td><td>2.2</td><td>2.2</td><td>1.7</td><td>1.8</td><td>2.2</td><td>0.5</td><td>1.2</td><td>1.3</td><td>1.1</td><td>1.1</td><td>1.4</td><td>1.0</td><td>0.3</td></tr><tr><td>slope</td><td>2.3</td><td>2.3</td><td>1.7</td><td>1.8</td><td>2.4</td><td>0.7</td><td>1.2</td><td>1.5</td><td>1.1</td><td>1.2</td><td>1.5</td><td>1.0</td><td>0.4</td></tr><tr><td>hollow</td><td>2.8</td><td>2.8</td><td>1.9</td><td>2.0</td><td>2.9</td><td>0.9</td><td>1.4</td><td>1.9</td><td>1.3</td><td>1.4</td><td>2.1</td><td>1.1</td><td>0.8</td></tr><tr><td>valley</td><td>3.7</td><td>3.5</td><td>2.4</td><td>2.5</td><td>3.8</td><td>1.4</td><td>1.6</td><td>2.1</td><td>1.5</td><td>1.6</td><td>2.3</td><td>1.2</td><td>0.8</td></tr><tr><td>pit</td><td>4.5</td><td>3.9</td><td>2.6</td><td>2.7</td><td>4.0</td><td>1.9</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr></tbody></table>

Using the TDX20/30 generally yielded the lowest IOR among global DEMs, whereas using TDX90 yielded the highest. The absolute minimum was always observed when using the ALS DEM (0.1–0.4 dB lower when compared to TDX20/30). IOR values obtained with the ALS DEM for bare areas were up to 3.9 dB lower when compared to the TDX90 DEM and up to 1.1 dB lower when compared with SRTM. For needleleaf forest, IOR values for valleys using ALS DEM were up to 1.1 dB lower than those obtained with TDX90.

### 4.3. Differences with and ALS-Derived DEM

The PNOA ALS-derived DEM provided the lowest IOR in all previous analyses at the Spanish site and was used as a reference for quantitative analysis of the global DEMs (Table 6). For the scattering area, the highest deviation (RMSD) was observed for the TDX90, followed by the SRTM DEMs. The relative RMSD obtained with these DEMs was at least 10% higher when compared to the remaining DEMs (TDX20, TDX30, and AW3D). The lowest RMSD were observed for AW3D (orbit 1) and TDX20 (orbit 81), followed by TDX30. The MAD for the scattering area was higher when using the SRTM or the TDX90 DEMs when compared to AW3D, TDX30, and TDX20. For both orbits, the SRTM-derived scattering area was the least biased when compared to the ALS DEM (–4 m² for orbit 1 and 0.4 m² for orbit 81), followed by TDX20 (–11 m² for orbit 1 and –12 m² for orbit 81).

<table><caption>**Table 6.** Quality assessment for needleleaf forests using PNOA as a reference. Cell color shows the gradient between the lowest (green) and the highest value (yellow).</caption><thead><tr><th rowspan="3"></th><th rowspan="3">Statistic</th><th colspan="5">O001</th><th colspan="5">O081</th></tr><tr><th colspan="2">AW</th><th colspan="2">SR</th><th colspan="2">TDX</th><th colspan="2">AW</th><th colspan="2">SR</th><th colspan="2">TDX</th></tr><tr><th>3D</th><th>TM</th><th>20</th><th>30</th><th>90</th><th>3D</th><th>TM</th><th>20</th><th>30</th><th>90</th><th>20</th><th>30</th></tr></thead><tbody><tr><td rowspan="4">Sc. area (m²)</td><td>Abs. RMSD</td><td>97</td><td>188</td><td>102</td><td>109</td><td>213</td><td>115</td><td>214</td><td>113</td><td>123</td><td>245</td><td>245</td><td>245</td></tr><tr><td>Rel. RMSD</td><td>11%</td><td>21%</td><td>11%</td><td>12%</td><td>24%</td><td>12%</td><td>22%</td><td>11%</td><td>12%</td><td>25%</td><td>25%</td><td>25%</td></tr><tr><td>MAD</td><td>38</td><td>81</td><td>32</td><td>36</td><td>87</td><td>47</td><td>96</td><td>38</td><td>43</td><td>104</td><td>104</td><td>104</td></tr><tr><td>Offset</td><td>18</td><td>-4</td><td>-11</td><td>14</td><td>-21</td><td>25</td><td>0.4</td><td>-12</td><td>20</td><td>-18</td><td>-18</td><td>-18</td></tr><tr><td rowspan="4">VV (dB)</td><td>Abs. RMSD</td><td>0.54</td><td>0.96</td><td>0.63</td><td>0.54</td><td>2.89</td><td>0.56</td><td>0.96</td><td>0.64</td><td>0.55</td><td>3.35</td><td>3.35</td><td>3.35</td></tr><tr><td>Rel. RMSD</td><td>5%</td><td>10%</td><td>6%</td><td>5%</td><td>29%</td><td>6%</td><td>10%</td><td>7%</td><td>6%</td><td>34%</td><td>34%</td><td>34%</td></tr><tr><td>MAD</td><td>0.38</td><td>0.64</td><td>0.41</td><td>0.36</td><td>0.74</td><td>0.38</td><td>0.63</td><td>0.41</td><td>0.35</td><td>0.75</td><td>0.75</td><td>0.75</td></tr><tr><td>Offset</td><td>-0.19</td><td>-0.27</td><td>-0.27</td><td>-0.17</td><td>-0.35</td><td>-0.24</td><td>-0.30</td><td>-0.35</td><td>-0.23</td><td>-0.41</td><td>-0.41</td><td>-0.41</td></tr></tbody></table>
For both orbits (1 and 81), the use of AW3D, TDX30, and TDX20 DEMs resulted in backscatter coefficient (VV) values closest to those obtained using the ALS DEM, with relative RMSD under 8%. For both orbits, the smallest offset with respect to the ALS DEM was obtained using TDX30, followed by AW3D. In all cases, backscattering coefficient was underestimated when compared with ALS DEM results.

## 4.4. Land Cover Classification

The DEM used for radiometric terrain normalization and topographic phase removal showed little effect on the overall quality of the classification, regardless of the Sentinel-1 relative orbit (Table A1). The Cohen's Kappa was between 0.94 and 0.96. Analyzing the confusion matrices and the associated error metrics (i.e., user accuracy, UE, producer accuracy, PE) showed that classes with a larger spatial extent (low vegetation or forest) had very high accuracies (>95%). However, most pixels misclassified as low vegetation (>85%) were, in fact, forest pixels according to the reference data. Water had a reasonable accuracy (>80%). Depending on the DEM, the "source class" of pixels misclassified water pixels varied. When using TDX20 for normalization, 21-27% of misclassified water pixels were recorded as forest pixels on the validation dataset, 20-33% with TDX30, 33-50% with SRTM or AW3D, and 36-58% with TDX90. The urban class had relatively small omission errors (10-15%), but the commission error was high (around 40% for orbit 7, around 60% for orbit 29, and around 50% for orbit 31), mostly due to the misclassification of forests (52-78% of the pixels misclassified as urban were forest pixels in the validation dataset).

Classification results were also analyzed by landform (Table A2), and in particular, for valleys. Valleys were selected as they showed the largest differences between DEMs in previous tests, and a reasonable number of samples were available (10 times more when compared to the number of samples available for pit landform). When only valley pixels were considered, Cohen Kappa was low (0.57-0.70). Using the TDX20 DEM allows for a marginal increase of the Cohen Kappa values over the value obtained using the rest of the DEMs (0.05 for orbits 7, 0.03 for orbit 131, and 0.02 for orbit 29). For valleys, user accuracy for low vegetation class dropped. In this context, the use of the TDX20 DEM reduced commission errors up to 11% for low vegetation and up to 21.9% for water when compared to the remaining DEMs (Figure 5). Using the TDX20 DEM also reduced the number of forest pixels misclassified as water, representing a smaller percentage of the pixels misclassified as such (17-27% less, depending on the orbit). Commission errors for urban cover increased for valley landform, especially when terrain normalization is performed with any of the TDX DEMs.
<figure><img src="image_8.png" /><figcaption>**Figure 5.** A small subset of the data around Leaota Peak. The first row shows the impact of terrain normalization on the imagery (a,b). The dotted box is the area shown for classification maps. The second row represents classified maps (c,d). White pixels indicate no data.</figcaption></figure>

The backscatter coefficient (VV) for forests located on valleys was examined to better understand the results obtained with the TDX DEMs (Figure 6). The boxplot showed that products normalized using lower resolution DEMs (AW3D, SRTM, TDX90) had an increased frequency of low values on the valley when compared to TDX20/30.
<figure><img src="image_9.png" /><figcaption>**Figure 6.** Boxplot representing the VV backscatter coefficient for forests located on valleys by Sentinel-1 relative orbit: Mean value (triangle) median value (orange) and inter quantile-ranges (whiskers) for 5–95%.</figcaption></figure>

**5. Discussion**

The influence of the DEM employed for terrain normalization of backscatter and coherence data variability was analyzed in three ways: (a) Comparing several orbital tracks (inter-orbit range, IOR); (b) using the results obtained with an ALS-derived DEM as a reference; and (c) assessing land cover classification results after a specific DEM is employed for normalization. Coherence varied very little with the DEM employed, whereas the effect was larger on the backscatter coefficient.

Terrain normalization was better served by high-resolution DEMs (i.e., TDX20, ALS DEM, AW3D at Spanish site), in agreement with prior research [17,29] (SRTM-1arcsec and AW3D outperformed SRTM-3arcsec and TanDEM-X 90m). Higher resolution DEMs have reduced vertical uncertainties (under 5 m over sloping terrain for TDX12.5 and AW3D) when compared to the SRTM DEM (Table 1), which may have contributed to reducing IOR values. Some DEMs (TDX20, TDX90, ALS) needed resampling prior, which may have impacted their performance. TDX20 and ALS DEMs, were down-sampled, which may have reduced the DEM detail with an associated increase of IOR values. However, the IOR values obtained after resampling was still smaller than those observed for the lower resolution DEMs, underlining the importance of the vertical uncertainty of the original DEM. TDX90 was resampled to a finer resolution. However, as resampling is a destructive operation, the only expected impact was the smooth interpolation of the original data to a denser grid, which does not provide additional information over the original DEM. The advantage provided by high-resolution DEMs was dependent on the specific land cover, and the landform it occupies. For instance, IOR for urban and crops showed little difference, as they occupy near-flat areas. Land cover classes occupying steeper slopes (i.e., forests, grasslands) received the largest benefits of using a more detailed DEM (minimum IOR). At the Romanian site, broadleaf forests showed smaller differences than mixed and needleleaf forest, as the latter grew on steeper slopes. Variability for needleleaf forests was smaller at the Spanish site, as it occupied milder slopes (Figure 3). Results were disaggregated by landform, as the “typical” slope of each land cover might obscure landform related effects. Peak, hollow, valley, and pit landforms showed the largest differences between the analyzed DEMs. Such differences can be attributed to the sensitivity of SAR to remote sensing artifacts, such as shadowing and foreshortening appearing with increasing slope. This affects DEM accuracy in sloped terrain, as shown in Table 1, and propagates into any analysis based on pixel neighborhood, such as slope, orientation [12], or terrain normalization. On these landforms, TDX20/30 clearly outperformed the rest of the DEMs, as it provided an improved characterization of smaller terrain forms supporting results reported by Grohmann et al. [17].
The use of an ALS-derived DEM resulted in the smallest IOR, pointing it as a suitable candidate to benchmark global DEMs at the Spanish site. Taking the ALS DEM as a reference, the lowest deviation and bias of the SAR metrics were observed for the AW3D and TDX20/30 DEMs. The results obtained by AW3D were explained by the combined effect of the resolution employed for its generation (5 m) and the low cover and height of Mediterranean forests. These factors may have eased the detection of vegetation-free pixels, "pushing" the reported data nearer to the true terrain surface once resampled to 30 m, as described by References [17,18]. TDX20/30 provided similar results thanks to the high spatial resolution of the X-band sensor employed for its creation. In addition, shorter wavelengths (X- as opposed to C-band) can capture finer spatial details [1,14,72]. In some cases, better results were observed for the TDX30 DEM when compared to the TDX20 DEM, pointing to a trade-off between the spatial detail and the effect of the DEM noise, with some improvement for slightly coarser resolutions (30 m instead of 20), but a high dispersion when resolution becomes too coarse (TDX 90). TDX90 and SRTM DEMs showed similar dispersion when compared to the ALS PNOA DEM reference values, possibly due to the variable-resolution smoothing [20] employed during the reprojection of SRTM to map coordinates, which may have decreased its detail, as shown in References [17,21,22].

DEM performance also varied across sites, with a larger inter-track variability being observed for needleleaf forests at the Romanian site. Such differences were explained by the steeper slopes this land cover class occupied in the Carpathians, as well as its characteristics (height and structure), which may have complicated DEM generation, due to volume decorrelation, as reported by References [20,24]: The Carpathians are covered by dense temperate forests, whereas sierra Nevada is populated by Mediterranean forests with lower tree height and canopy density. Among all DEMs, AW3D IOR had a distinct behavior. While at the Romanian site, the AW3D results were close to those observed to the SRTM DEM, at the Spanish, the results were closer to those observed when using the TDX20/30 DEMs. The quality mask layer showed that, at the Romanian site, the AW3D DEM has a large strip where missing data (due persistent cloud cover) have been infilled with SRTM data, a problem mentioned by Truckenbrodt [29]. For this reason, it is not possible to draw conclusions on the influence of canopy characteristics on the performance of AW3D.

The characteristics of each DEM propagated into the land cover classification results. Overall classification accuracy was similar regardless of the DEM, with reasonable accuracies for all classes except urban. The low accuracy of urban surfaces was attributed to (i) the prevalence of steep slopes, which difficult terrain normalization (not fully accounted scattering area) [29] and may introduce DEM artifacts (mischaracterizing terrain surface) [12,14,72], and (ii) the prevalence of discontinuous urban fabric at the Romanian site, which may cause confusion with ornamental and fruit tree cover present around residential areas, and genuine forested lands.

When analyzing classification results for specific landforms (i.e., valleys), the AW3D, SRTM and TDX90 showed larger commission errors (CE), due to the misclassification of the forest as low vegetation or water. Such errors can be explained by the mischaracterization of the thin crevices of the drainage network on the DEMs [12]. This is propagated to the lookup table and impacts both, the distortion masking process and the terrain normalization. Distortion masking is affected because the distance between the pixels in the range is altered. Therefore, such pixels are not marked as distorted or shadowed. Terrain normalization is affected as it uses the LUT, orientation, and slope layers to estimate the scattering area. In turn, the scattering area is overestimated, under-compensating the radiometric effects, keeping a pseudo-shadow with lowered values (as opposed to true shadow, caused by occlusion, which cannot be compensated). Therefore, forests are misclassified due to the apparently low backscatter.

Using the TDX20 reduced commission error for low vegetation (11%) and water (12–22%) on valleys, as well as the percentage of forest pixels among all pixels misclassified as water (decrease of 17–27%). This was explained by reduced pseudo-shadows when using the TDX20/30 DEMs. However, the improvement came at the cost of a 3–6% increase of commission error for urban cover located on valleys, which was caused by a slight increase in backscatter values over forests (overcompensation),
and the large size difference between the validation sample for urban and forest (<10³ vs. 7.5·10⁴ at valleys), as misclassification of a small subset of the latter would be much larger when compared with the sample size for the former. Even with this trade-off, classification results using TDX20 were marginally better (0.01–0.05 higher Cohen’s Kappa). Furthermore, increased CE for urban cover on valleys did not affect the overall CE for the urban class, which was reduced by 2–3% when using the TDX20 DEM.

## 6. Conclusions

SAR observations are heavily affected by sensor-terrain geometry, which can be corrected using a DEM. Choosing a DEM for SAR data terrain normalization is not a trivial choice, as it affects backscattering coefficient variability, and mapping products generated downstream. High-resolution TanDEM-X DEM (20 or 30 m resolution) was the global DEM providing the largest reduction of terrain induced variability, followed by AW3D in sparse vegetation areas. Natural land covers (i.e., forest, bare areas, grasslands) occupying steeper slopes and complex landforms (i.e., peaks, pits, valleys) received the largest benefits. These benefits were felt on classification, where more forest pixels were classified correctly due to a better compensation of low values (valley pseudo-shadow). An ALS-based DEM was able to provide slightly better results (i.e., marginally reduced IOR) when compared to AW3D and TDX20/30 DEMs. However, AW3D and TDX20/30 DEMs seem suitable candidates to replace ALS-based local DEMs. However, AW3D should be checked for data infilling from older datasets (i.e., SRTM) as over such areas, its performance may be degraded.

This study showed the effect of several global DEMs on terrain normalization, highlighting their advantages and shortcomings when normalizing Sentinel-1 imagery. Further research should expand this analysis by including the recent NASADEM dataset (from re-processed SRTM) [73], using a reference ALS-based DEM for temperate forests, and studying the DEM-dependent normalization effects on SAR imagery acquired at different wavelengths. Finally, the effect of terrain normalization could be tested on downstream quantitative products, such as biomass estimates or canopy cover.

**Author Contributions:** Conceptualization and Methodology, I.B.-M. and M.A.T.; software, validation, formal analysis, investigation, data curation, visualization, I.B.-M.; resources, O.B., M.A.T.; writing—original draft preparation, I.B.-M.; writing—review and editing, I.B.-M., M.A.T., M.S. and L.V.; supervision, project administration and funding acquisition, O.B. and M.A.T. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research was funded by the Romanian National Authority for Scientific Research and Innovation and the European Regional Development Fund under the project "Prototyping an Earth-Observation based monitoring and forecasting system for the Romanian forests" (EO-ROFORMON, project ID P_37_651/105058).

**Acknowledgments:** TanDEM-X DEM data were provided by the German Aerospace Center (DLR), under the proposal DEM_FOREST2614. We also would like to thank the publishers of the open datasets employed in this study: DLR for making available the TanDEM-X Global Urban Footprint dataset, European Environment Agency for providing the Corine Land Cover Dataset, European Space Agency for providing Sentinel-1 data and CCI land cover map, National Geographic Institute of Spain for providing PNOA LiDAR DEM and SIOSE land cover dataset, JAXA, for providing both AW3D DEM and ALOS forest/non-forest map, and USGS, for providing SRTM. Finally, we would like to thank INCDS Marin Dracea and University of Alcalá for hosting the research activities.

**Conflicts of Interest:** The authors declare no conflict of interest.
Appendix A

<table><caption>**Table A1.** Confusion matrices by DEM and orbit (Reference > columns; Classified > Rows) at overall level. (Ur, urban; LV, low vegetation; Fo, forest; Wa, water. UA, user accuracy, PA, producer accuracy).</caption><thead><tr><th rowspan="2"></th><th rowspan="2"></th><th colspan="6">O007</th><th colspan="6">O029</th><th colspan="6">O131</th></tr><tr><th>Ur</th><th>LV</th><th>Fo</th><th>Wa</th><th>UA</th><th>Ur</th><th>LV</th><th>Fo</th><th>Wa</th><th>UA</th><th>Ur</th><th>LV</th><th>Fo</th><th>Wa</th><th>UA</th></tr></thead><tbody><tr><td rowspan="5">AW3D</td><td>Ur</td><td>19147</td><td>6602</td><td>7310</td><td>43</td><td>57.84</td><td>18698</td><td>6465</td><td>19317</td><td>24</td><td>42.01</td><td>18749</td><td>5044</td><td>12423</td><td>44</td><td>51.71</td></tr><tr><td>LV</td><td>2206</td><td>1548500</td><td>19344</td><td>271</td><td>98.61</td><td>2063</td><td>1544334</td><td>20493</td><td>321</td><td>98.54</td><td>2366</td><td>1545097</td><td>24852</td><td>358</td><td>98.25</td></tr><tr><td>Fo</td><td>250</td><td>10075</td><td>811265</td><td>380</td><td>98.70</td><td>759</td><td>13244</td><td>787900</td><td>395</td><td>98.21</td><td>357</td><td>9861</td><td>739556</td><td>389</td><td>98.59</td></tr><tr><td>Wa</td><td>0</td><td>506</td><td>506</td><td>5054</td><td>83.32</td><td>0</td><td>611</td><td>355</td><td>5007</td><td>83.83</td><td>0</td><td>688</td><td>468</td><td>4956</td><td>81.09</td></tr><tr><td>PA</td><td>88.63</td><td>98.90</td><td>96.76</td><td>87.93</td><td></td><td>86.89</td><td>98.70</td><td>95.15</td><td>87.12</td><td></td><td>87.32</td><td>99.00</td><td>95.14</td><td>86.24</td><td></td></tr><tr><td rowspan="5">SRTM</td><td>Ur</td><td>18868</td><td>6245</td><td>7204</td><td>40</td><td>58.31</td><td>18439</td><td>6055</td><td>19377</td><td>24</td><td>42.01</td><td>18725</td><td>5048</td><td>12952</td><td>48</td><td>50.92</td></tr><tr><td>LV</td><td>2278</td><td>1540604</td><td>19090</td><td>262</td><td>98.62</td><td>2108</td><td>1537069</td><td>20272</td><td>318</td><td>98.54</td><td>2283</td><td>1538477</td><td>23950</td><td>340</td><td>98.30</td></tr><tr><td>Fo</td><td>265</td><td>9900</td><td>806632</td><td>379</td><td>98.71</td><td>774</td><td>12931</td><td>782981</td><td>386</td><td>98.23</td><td>285</td><td>9564</td><td>735257</td><td>371</td><td>98.63</td></tr><tr><td>Wa</td><td>0</td><td>546</td><td>549</td><td>4944</td><td>81.87</td><td>0</td><td>566</td><td>283</td><td>4895</td><td>85.22</td><td>0</td><td>618</td><td>358</td><td>4861</td><td>83.28</td></tr><tr><td>PA</td><td>88.12</td><td>98.93</td><td>96.78</td><td>87.89</td><td></td><td>86.48</td><td>98.74</td><td>95.15</td><td>87.05</td><td></td><td>87.94</td><td>99.02</td><td>95.18</td><td>86.49</td><td></td></tr><tr><td rowspan="5">TDX20</td><td>Ur</td><td>29990</td><td>9446</td><td>11850</td><td>73</td><td>58.39</td><td>28857</td><td>8597</td><td>30849</td><td>36</td><td>42.23</td><td>29327</td><td>7184</td><td>18985</td><td>73</td><td>52.78</td></tr><tr><td>LV</td><td>3168</td><td>2411865</td><td>22525</td><td>363</td><td>98.93</td><td>3440</td><td>2410056</td><td>31421</td><td>475</td><td>98.55</td><td>3426</td><td>2407209</td><td>39924</td><td>475</td><td>98.21</td></tr><tr><td>Fo</td><td>364</td><td>16915</td><td>1270072</td><td>619</td><td>98.61</td><td>1088</td><td>17561</td><td>1188573</td><td>552</td><td>98.41</td><td>439</td><td>12460</td><td>1060469</td><td>540</td><td>98.75</td></tr><tr><td>Wa</td><td>1</td><td>1112</td><td>419</td><td>7772</td><td>83.53</td><td>0</td><td>960</td><td>251</td><td>7535</td><td>86.15</td><td>0</td><td>1063</td><td>311</td><td>7085</td><td>83.76</td></tr><tr><td>PA</td><td>89.46</td><td>98.87</td><td>97.33</td><td>88.05</td><td></td><td>86.44</td><td>98.89</td><td>95.00</td><td>87.64</td><td></td><td>88.36</td><td>99.15</td><td>94.71</td><td>86.69</td><td></td></tr><tr><td rowspan="5">TDX30</td><td>Ur</td><td>19156</td><td>6430</td><td>8530</td><td>42</td><td>56.08</td><td>18677</td><td>6280</td><td>21653</td><td>24</td><td>40.05</td><td>18793</td><td>4847</td><td>13386</td><td>46</td><td>50.69</td></tr><tr><td>LV</td><td>2034</td><td>1540695</td><td>16591</td><td>239</td><td>98.79</td><td>2010</td><td>1538410</td><td>20705</td><td>301</td><td>98.53</td><td>2205</td><td>1539752</td><td>27089</td><td>303</td><td>98.11</td></tr><tr><td>Fo</td><td>221</td><td>9505</td><td>808033</td><td>390</td><td>98.76</td><td>642</td><td>11234</td><td>773058</td><td>350</td><td>98.44</td><td>276</td><td>8255</td><td>712706</td><td>356</td><td>98.77</td></tr><tr><td>Wa</td><td>0</td><td>665</td><td>321</td><td>4954</td><td>83.40</td><td>0</td><td>718</td><td>220</td><td>4867</td><td>83.84</td><td>0</td><td>680</td><td>171</td><td>4678</td><td>84.61</td></tr><tr><td>PA</td><td>89.47</td><td>98.93</td><td>96.95</td><td>88.07</td><td></td><td>87.57</td><td>98.83</td><td>94.78</td><td>87.82</td><td></td><td>88.34</td><td>99.11</td><td>94.60</td><td>86.90</td><td></td></tr><tr><td rowspan="5">TDX90</td><td>Ur</td><td>19027</td><td>6599</td><td>7449</td><td>46</td><td>57.45</td><td>18439</td><td>5990</td><td>18774</td><td>26</td><td>42.65</td><td>18727</td><td>5161</td><td>13647</td><td>53</td><td>49.82</td></tr><tr><td>LV</td><td>2179</td><td>1542343</td><td>18571</td><td>245</td><td>98.66</td><td>2138</td><td>1538972</td><td>19112</td><td>309</td><td>98.62</td><td>2294</td><td>1540926</td><td>23537</td><td>344</td><td>98.33</td></tr><tr><td>Fo</td><td>240</td><td>9780</td><td>808189</td><td>396</td><td>98.73</td><td>794</td><td>13421</td><td>789226</td><td>382</td><td>98.18</td><td>309</td><td>9608</td><td>743312</td><td>370</td><td>98.63</td></tr><tr><td>Wa</td><td>9</td><td>666</td><td>921</td><td>4937</td><td>75.57</td><td>0</td><td>567</td><td>343</td><td>4905</td><td>84.35</td><td>0</td><td>641</td><td>357</td><td>4854</td><td>82.95</td></tr><tr><td>PA</td><td>88.68</td><td>98.91</td><td>96.77</td><td>87.78</td><td></td><td>86.28</td><td>98.72</td><td>95.38</td><td>87.25</td><td></td><td>87.80</td><td>99.01</td><td>95.19</td><td>86.35</td><td></td></tr></tbody></table>
<table><caption>**Table A2.** Confusion matrices by DEM and orbit (Reference > columns; Classified > Rows) for valley. (Ur, urban; LV, low vegetation; Fo, forest; Wa, water. UA, user accuracy, PA, producer accuracy).</caption><thead><tr><th rowspan="2"></th><th colspan="6">O007</th><th colspan="6">O029</th><th colspan="6">O131</th></tr><tr><th>Ur</th><th>LV</th><th>Fo</th><th>Wa</th><th>UA</th><th>Ur</th><th>LV</th><th>Fo</th><th>Wa</th><th>UA</th><th>Ur</th><th>LV</th><th>Fo</th><th>Wa</th><th>UA</th></tr></thead><tbody><tr><td rowspan="5">AW3D</td><td>Ur</td><td>304</td><td>188</td><td>827</td><td>2</td><td>23.01</td><td>305</td><td>225</td><td>2478</td><td>1</td><td>10.14</td><td>298</td><td>175</td><td>1292</td><td>2</td><td>16.86</td></tr><tr><td>LV</td><td>25</td><td>5447</td><td>4117</td><td>26</td><td>56.65</td><td>20</td><td>5232</td><td>4292</td><td>30</td><td>54.65</td><td>27</td><td>5214</td><td>3984</td><td>33</td><td>56.32</td></tr><tr><td>Fo</td><td>2</td><td>537</td><td>76000</td><td>19</td><td>99.27</td><td>1</td><td>697</td><td>73596</td><td>18</td><td>99.04</td><td>3</td><td>606</td><td>70952</td><td>20</td><td>99.12</td></tr><tr><td>Wa</td><td>0</td><td>32</td><td>176</td><td>514</td><td>71.19</td><td>0</td><td>27</td><td>81</td><td>512</td><td>82.58</td><td>0</td><td>27</td><td>54</td><td>506</td><td>86.20</td></tr><tr><td>PA</td><td>91.84</td><td>87.80</td><td>93.69</td><td>91.62</td><td></td><td>93.56</td><td>84.65</td><td>91.48</td><td>91.27</td><td></td><td>90.85</td><td>86.58</td><td>93.01</td><td>90.20</td><td></td></tr><tr><td rowspan="5">SRTM</td><td>Ur</td><td>319</td><td>239</td><td>816</td><td>1</td><td>23.20</td><td>318</td><td>193</td><td>1974</td><td>2</td><td>12.79</td><td>318</td><td>177</td><td>949</td><td>1</td><td>22.01</td></tr><tr><td>LV</td><td>29</td><td>5195</td><td>3821</td><td>18</td><td>57.32</td><td>30</td><td>5137</td><td>4540</td><td>28</td><td>52.77</td><td>29</td><td>5102</td><td>4404</td><td>28</td><td>53.35</td></tr><tr><td>Fo</td><td>3</td><td>598</td><td>72313</td><td>17</td><td>99.15</td><td>0</td><td>689</td><td>70070</td><td>14</td><td>99.01</td><td>1</td><td>583</td><td>68204</td><td>15</td><td>99.13</td></tr><tr><td>Wa</td><td>0</td><td>29</td><td>128</td><td>429</td><td>73.21</td><td>0</td><td>27</td><td>84</td><td>421</td><td>79.14</td><td>0</td><td>29</td><td>58</td><td>420</td><td>82.84</td></tr><tr><td>PA</td><td>90.88</td><td>85.71</td><td>93.82</td><td>92.26</td><td></td><td>91.38</td><td>84.97</td><td>91.39</td><td>90.54</td><td></td><td>91.38</td><td>86.61</td><td>92.65</td><td>90.52</td><td></td></tr><tr><td rowspan="5">TDX20</td><td>Ur</td><td>495</td><td>410</td><td>1795</td><td>1</td><td>18.33</td><td>492</td><td>371</td><td>4221</td><td>2</td><td>9.67</td><td>485</td><td>277</td><td>1923</td><td>1</td><td>18.06</td></tr><tr><td>LV</td><td>35</td><td>8150</td><td>3874</td><td>29</td><td>67.42</td><td>35</td><td>8080</td><td>5511</td><td>47</td><td>59.09</td><td>33</td><td>7983</td><td>5496</td><td>44</td><td>58.89</td></tr><tr><td>Fo</td><td>3</td><td>994</td><td>113029</td><td>34</td><td>99.10</td><td>1</td><td>1022</td><td>104714</td><td>21</td><td>99.01</td><td>4</td><td>844</td><td>95033</td><td>26</td><td>99.09</td></tr><tr><td>Wa</td><td>0</td><td>38</td><td>108</td><td>724</td><td>83.22</td><td>0</td><td>31</td><td>36</td><td>673</td><td>90.95</td><td>0</td><td>27</td><td>48</td><td>596</td><td>88.82</td></tr><tr><td>PA</td><td>92.87</td><td>84.97</td><td>95.14</td><td>91.88</td><td></td><td>93.18</td><td>85.02</td><td>91.47</td><td>90.58</td><td></td><td>92.91</td><td>87.43</td><td>92.72</td><td>89.36</td><td></td></tr><tr><td rowspan="5">TDX30</td><td>Ur</td><td>329</td><td>274</td><td>1267</td><td>1</td><td>17.58</td><td>323</td><td>258</td><td>2910</td><td>1</td><td>9.25</td><td>324</td><td>205</td><td>1377</td><td>1</td><td>16.99</td></tr><tr><td>LV</td><td>19</td><td>5160</td><td>2846</td><td>15</td><td>64.18</td><td>25</td><td>5091</td><td>3604</td><td>24</td><td>58.22</td><td>19</td><td>5068</td><td>3732</td><td>23</td><td>57.32</td></tr><tr><td>Fo</td><td>3</td><td>594</td><td>72894</td><td>21</td><td>99.16</td><td>0</td><td>651</td><td>68921</td><td>12</td><td>99.05</td><td>1</td><td>553</td><td>64973</td><td>18</td><td>99.13</td></tr><tr><td>Wa</td><td>0</td><td>33</td><td>71</td><td>428</td><td>80.45</td><td>0</td><td>29</td><td>51</td><td>411</td><td>83.71</td><td>0</td><td>29</td><td>33</td><td>378</td><td>85.91</td></tr><tr><td>PA</td><td>93.73</td><td>85.13</td><td>94.57</td><td>92.04</td><td></td><td>92.82</td><td>84.44</td><td>91.30</td><td>91.74</td><td></td><td>94.19</td><td>86.56</td><td>92.67</td><td>90.00</td><td></td></tr><tr><td rowspan="5">TDX90</td><td>Ur</td><td>305</td><td>258</td><td>1046</td><td>0</td><td>18.96</td><td>304</td><td>243</td><td>2341</td><td>0</td><td>10.53</td><td>309</td><td>187</td><td>1256</td><td>0</td><td>17.64</td></tr><tr><td>LV</td><td>24</td><td>5299</td><td>3304</td><td>21</td><td>61.27</td><td>22</td><td>5183</td><td>3809</td><td>39</td><td>57.25</td><td>21</td><td>5216</td><td>3669</td><td>39</td><td>58.31</td></tr><tr><td>Fo</td><td>3</td><td>606</td><td>70640</td><td>19</td><td>99.12</td><td>2</td><td>725</td><td>68702</td><td>14</td><td>98.93</td><td>2</td><td>610</td><td>67360</td><td>12</td><td>99.08</td></tr><tr><td>Wa</td><td>0</td><td>26</td><td>259</td><td>451</td><td>61.28</td><td>0</td><td>28</td><td>116</td><td>437</td><td>75.22</td><td>0</td><td>25</td><td>106</td><td>439</td><td>77.02</td></tr><tr><td>PA</td><td>91.87</td><td>85.62</td><td>93.88</td><td>91.85</td><td></td><td>92.68</td><td>83.88</td><td>91.64</td><td>89.18</td><td></td><td>93.07</td><td>86.39</td><td>93.05</td><td>89.59</td><td></td></tr></tbody></table>
References

1. Ouchi, K. Recent trend and advance of synthetic aperture radar with selected topics. *Remote Sens.* **2013**, *5*, 716–807. [CrossRef]

2. Small, D. Flattening gamma: Radiometric terrain correction for SAR imagery. *IEEE Trans. Geosci. Remote Sens.* **2011**, *49*, 3081–3093. [CrossRef]

3. Frey, O.; Santoro, M.; Werner, C.L.; Wegmuller, U. DEM-based SAR pixel-area estimation for enhanced geocoding refinement and radiometric normalization. *IEEE Geosci. Remote Sens. Lett.* **2013**, *10*, 48–52. [CrossRef]

4. Askne, J.I.H.; Dammert, P.B.G.; Ulander, L.M.H.; Smith, G. C-band repeat-pass interferometric SAR observations of the forest. *IEEE Trans. Geosci. Remote Sens.* **1997**, *35*, 25–35. [CrossRef]

5. Lefsky, M.A.; Cohen, W.B.; Parker, G.G.; Harding, D.J. Lidar remote sensing for ecosystem studies. *BioScience* **2002**, *52*, 192002. [CrossRef]

6. CNIG. Centro Nacional de Información Geográfica National Plan for Aerial Orthophotography. Available online: http://pnoa.ign.es/ (accessed on 21 May 2019).

7. ODP Aerial LiDAR Scanning Datasets. Available online: http://data.opendataportal.at/dataset?tags=ALS&sort=score+desc (accessed on 29 December 2019).

8. Shuman, C.A.; Zwally, H.J.; Schutz, B.E.; Brenner, A.C.; DiMarzio, J.P.; Suchdeo, V.P.; Fricker, H.A. ICESat Antarctic elevation data: Preliminary precision and accuracy assessment. *Geophys. Res. Lett.* **2006**, *33*, L07501. [CrossRef]

9. Tadono, T.; Ishida, H.; Oda, F.; Naito, S.; Minakawa, K.; Iwamoto, H. Precise global DEM generation by ALOS PRISM. *ISPRS Ann. Photogramm. Remote Sens. Spat. Inf. Sci.* **2014**, *2*, 71–76. [CrossRef]

10. Felicísimo, A.M. *Modelos Digitales del Terreno: Introducción y Aplicaciones en las Ciencias Ambientales*; Pentalfa Ediciones: Oviedo, Spain, 1994; ISBN 978-84-7848-475-1.

11. Aber, J.S.; Marzolff, I.; Ries, J.B.; Aber, S.E.W. Principles of photogrammetry. In *Small-Format Aerial Photography and UAS Imagery*; Elsevier: Amsterdam, The Netherlands, 2010.

12. Purinton, B.; Bookhagen, B. Validation of digital elevation models (DEMs) and comparison of geomorphic metrics on the southern Central Andean Plateau. *Earth Surf. Dynam.* **2017**, *5*, 211–237. [CrossRef]

13. Tachikawa, T.; Hato, M.; Kaku, M.; Iwasaki, A. Characteristics of ASTER GDEM version 2. In *Proceedings of the IEEE International Geoscience and Remote Sensing Symposium, Vancouver, BC, Canada, 24–29 July 2011*; pp. 3657–3660.

14. Woodhouse, I.H. *Introduction to Microwave Remote Sensing*; Taylor&Francis: Boca Raton, FL, USA, 2006; ISBN 978-0-415-27123-3.

15. Hirt, C.; Filmer, M.S.; Featherstone, W.E. Comparison and validation of the recent freely available ASTER-GDEM ver1, SRTM ver4.1 and GEODATA DEM-9S ver3 digital elevation models over Australia. *Aust. J. Earth Sci.* **2010**, *57*, 337–347. [CrossRef]

16. Florinsky, I.V.; Skrypitsyna, T.N.; Luschikova, O.S. Comparative accuracy of the AW3D30 DSM, ASTER GDEM, and SRTM1 DEM: A case study on the Zaoksky testing ground, Central European Russia. *Remote Sens. Lett.* **2018**, *9*, 706–714. [CrossRef]

17. Grohmann, C.H. Evaluation of TanDEM-X DEMs on selected Brazilian sites: Comparison with SRTM, ASTER GDEM and ALOS AW3D30. *Remote Sens. Environ.* **2018**, *212*, 121–133. [CrossRef]

18. Santillan, J.R.; Makinano-Santillan, M. Vertical accuracy assessment of 30-m resolution ALOS, ASTER, AND SRTM global DEMs over northeastern Mindanao, Philippines. *Int. Arch. Photogramm. Remote Sens. Spat. Inf. Sci.* **2016**, *41*, 149–156. [CrossRef]

19. Courty, L.G.; Soriano-Monzalvo, J.C.; Pedrozo-Acuña, A. Evaluation of open-access global digital elevation models (AW3D30, SRTM, and ASTER) for flood modelling purposes. *J. Flood Risk Manag.* **2019**, *12*, e12550. [CrossRef]

20. Farr, T.G.; Rosen, P.A.; Caro, E.; Crippen, R.; Duren, R.; Hensley, S.; Kobrick, M.; Paller, M.; Rodriguez, E.; Roth, L.; et al. The shuttle radar topography mission. *Rev. Geophys.* **2007**, *45*. [CrossRef]

21. Smith, B. Accuracy and resolution of shuttle radar topography mission data. *Geophys. Res. Lett.* **2003**, *30*, 1467. [CrossRef]

22. Guth, P.L. Geomorphometry from SRTM. *Photogramm. Eng Remote Sens.* **2006**, *72*, 269–277. [CrossRef]
23. Wessel, B. *TanDEM-X Ground Segment-DEM Products Specification Document*; DLR: Oberpfaffenhofen, Germany, 2016.

24. Rizzoli, P.; Martone, M.; Gonzalez, C.; Wecklich, C.; Borla Tridon, D.; Bräutigam, B.; Bachmann, M.; Schulze, D.; Fritz, T.; Huber, M.; et al. Generation and performance assessment of the global TanDEM-X digital elevation model. *ISPRS J. Photogramm. Remote Sens.* **2017**, *132*, 119–139. [CrossRef]

25. Wallington, E.D.; Woodhouse, I.H. Forest height retrieval from commercial X-band SAR products. *IEEE Trans. Geosci. Remote Sens.* **2006**, *44*, 863–870. [CrossRef]

26. Santoro, M.; Askne, J.; Dammert, P.B.G. Tree height influence on ERS interferometric phase in boreal forest. *IEEE Trans. Geosci. Remote Sens.* **2005**, *43*, 207–217. [CrossRef]

27. Thirion-Lefevre, L.; Colin-Koeniguer, E. Investigating attenuation, scattering phase center, and total height using simulated interferometric SAR images of forested areas. *IEEE Trans. Geosci. Remote Sens.* **2007**, *45*, 3172–3179. [CrossRef]

28. Hoekman, D.H.; Reiche, J. Multi-model radiometric slope correction of SAR images of complex terrain using a two-stage semi-empirical approach. *Remote Sens. Environ.* **2015**, *156*, 1–10. [CrossRef]

29. Truckenbrodt, J.; Freemantle, T.; Williams, C.; Jones, T.; Small, D.; Dubois, C.; Thiel, C.; Rossi, C.; Syriou, A.; Giuliani, G. Towards sentinel-1 SAR analysis-ready data: A best practices assessment on preparing backscatter data for the cube. *Data* **2019**, *4*, 932019. [CrossRef]

30. DLR. TanDEM-X DEM Proposal Submission Procedure. Available online: https://tandemx-science.dlr.de/cgi-bin/wcm.pl?page=TDM-Proposal-Submission-Procedure (accessed on 8 September 2020).

31. EROS. *Shuttle Radar Topography Mission (SRTM) 1 Arc-Second Global*; EROS: Sioux Falls, SD, USA, 2017.

32. USGS EarthExplorer. Available online: https://earthexplorer.usgs.gov/ (accessed on 13 January 2020).

33. JAXA ALOS Global Digital Surface Model "ALOS World 3D-30m" (AW3D30). Available online: https://www.eorc.jaxa.jp/ALOS/en/aw3d30/index.htm (accessed on 13 January 2020).

34. DLR EOWEB GeoPortal. Available online: https://eoweb.dlr.de/egp/ (accessed on 13 January 2020).

35. CNIG. Centro Nacional de Información Geográfica Centro de Descargas del CNIG (IGN). Available online: http://centrodedescargas.cnig.es (accessed on 13 January 2020).

36. CNIG. Centro Nacional de Información Geográfica Primera Cobertura del Plan Nacional de Fotografía Aérea. Available online: https://pnoa.ign.es/estado-del-proyecto-lidar/primera-cobertura (accessed on 21 August 2020).

37. Soininen, A. *TerraScan User's Guide*; Terrasolid: Helsinki, Finland, 2004.

38. Lorite Martínez, S.; Moreno Jabato, J.; de Tejada, J.M.G.S.; Rodríguez-Cuenca, B. Automatic classification of bridges and continental water bodies from 3D point clouds (aerial lidar). *Int. Arch. Photogramm. Remote Sens. Spat. Inf. Sci.* **2019**, *42*, 1047–1051. [CrossRef]

39. Lorite Martínez, S.; Ojeda Manrique, J.C.; Rodríguez-Cuenca, B.; González Cristóbal, E.; Muñoz, P. Procesado y distribución de nubes de puntos en el proyecto PNOA-LiDAR. In *Proceedings of the XVII Congreso de la Asociación Española de Teledetección, Murcia, Spain*, 3–7 October 2017; pp. 329–332.

40. Rodríguez, E.; Morris, C.S.; Belz, J.E. A Global Assessment of the SRTM Performance. *Photogramm. Eng. Remote Sens.* **2006**, *72*, 249–260. [CrossRef]

41. Takaku, J.; Tadono, T.; Tsutsui, K.; Ichikawa, M. Validation of "AW3D" global DSM generated from Alos Prism. *ISPRS Ann. Photogramm. Remote Sens. Spat. Inf. Sci.* **2016**, *3*, 25–31. [CrossRef]

42. IGN. Plan Nacional de Ortofotografía Aérea. Especificaciones Técnicas. Available online: http://pnoa.ign.es/caracteristicas-tecnicas (accessed on 21 May 2019).

43. GDAL/OGR Contributors. *GDAL/OGR Geospatial Data Abstraction Software Library*; Open Source Geospatial Foundation: Beaverton, OR, USA, 2020.

44. Wegnüller, U.; Werner, C.; Strozzi, T.; Wiesmann, A.; Frey, O.; Santoro, M. Sentinel-1 support in the GAMMA software. *Proc. Comput. Sci.* **2016**, *100*, 1305–1312. [CrossRef]

45. GRASS Development Team. *Geographic Resources Analysis Support System (GRASS GIS) Software*; Open Source Geospatial Foundation: Beaverton, OR, USA, 2017.

46. Python Software Foundation. *Python Language Reference, Version 3*; Python Software Foundation: Fredericksburg, VA, USA, 2020.

47. Gillies, S.; Ward, B.; Petersen, A.S. *Rasterio: Geospatial raster I/O for Python Programmers*; Mapbox: Washington, DC, USA, 2013.

48. Pandas Contributors. *Pandas*; Zenodo: Geneva, Switzerland, 2020.
49. GeoPandas Contributors. *GeoPandas: Python Tools for Geographic Data*; Zenodo: Geneva, Switzerland, 2020.

50. Van der Walt, S.; Colbert, S.C.; Varoquaux, G. The NumPy array: A structure for efficient numerical computation. *Comput. Sci. Eng.* **2011**, *13*, 22–30. [CrossRef]

51. Jones, E.; Oliphant, T.; Peterson, P. *SciPy: Open Source Scientific Tools for Python*; Python Software Foundation: Fredericksburg, VA, USA, 2014.

52. Hunter, J.D. Matplotlib: A 2D graphics environment. *Comput. Sci. Eng.* **2007**, *9*, 90–95. [CrossRef]

53. Wegmüller, U.; Werner, C.; Strozzi, T.; Wiesmann, A. Automated and precise image registration procedures. In *Analysis of Multi-Temporal Remote Sensing Images*; Series in Remote Sensing; World Scientific: Singapore, 2002; Volume 2, pp. 37–49. ISBN 978-981-02-4955-7.

54. Wegmüller, U.; Werner, C. Land applications using ERS-1/2 tandem data. In Proceedings of the Fringe Workshop, ERS SAR Interferometry, European Space Agency, Zurich, Switzerland, 30 September–2 October 1996; pp. 97–112.

55. Werner, C.; Wegmüller, U.; Strozzi, T.; Wiesmann, A. Gamma SAR and interferometric processing software. In Proceedings of the ERS-Envisat Symposium, Gothenburg, Switzerland, 16–20 October 2000.

56. Touzi, R.; Lopes, A.; Bruniquel, J.; Vachon, P.W. Coherence estimation for SAR imagery. *IEEE Trans. Geosci. Remote Sens.* **1999**, *37*, 135–149. [CrossRef]

57. Castel, T.; Beaudoin, A.; Stach, N.; Stussi, N.; Le Toan, T.; Durand, P. Sensitivity of space-borne SAR data to forest parameters over sloping terrain. Theory and experiment. *Int. J. Remote Sens.* **2001**, *22*, 2351–2376. [CrossRef]

58. Kirches, G.; Brockmann, C.; Boettcher, M.; Peters, M.; Bontemps, S.; Lamarche, C.; Schlerf, M.; Santoro, M.; Defourny, P. *Land Cover CCI-Product User Guide-Version 2*; ESA: Libin, Belgium, 2017.

59. Esch, T.; Schenk, A.; Ullmann, T.; Thiel, M.; Roth, A.; Dech, S. Characterization of land cover types in TerraSAR-X images by combined analysis of speckle statistics and intensity information. *IEEE Trans. Geosci. Remote Sens.* **2011**, *49*, 1911–1925. [CrossRef]

60. Esch, T.; Heldens, W.; Hirner, A.; Keil, M.; Marconcini, M.; Roth, A.; Zeidler, J.; Dech, S.; Strano, E. Breaking new ground in mapping human settlements from space—The Global Urban Footprint. *ISPRS J. Photogramm. Remote Sens.* **2017**, *134*, 30–42. [CrossRef]

61. Esch, T.; Bachofer, F.; Heldens, W.; Hirner, A.; Marconcini, M.; Palacios-Lopez, D.; Roth, A.; Üreyen, S.; Zeidler, J.; Dech, S.; et al. Where we live—A summary of the achievements and planned evolution of the global urban footprint. *Remote Sens.* **2018**, *10*, 895. [CrossRef]

62. Shimada, M.; Itoh, T.; Motooka, T.; Watanabe, M.; Shiraishi, T.; Thapa, R.; Lucas, R. New global forest/non-forest maps from ALOS PALSAR data (2007–2010). *Remote Sens. Environ.* **2014**, *155*, 13–31. [CrossRef]

63. CRC. *European Landscape Dynamics: CORINE Land Cover Data*; Feranec, J., Soukup, T., Feranec, G., Jaffrain, G., Eds.; CRC Press: Boca Raton, FL, USA, 2016; ISBN 978-1-4822-4466-3.

64. Del Bosque González, I.; Arozarena Villar, A.; Villa Alcázar, G.; Valcárcel Sanz, N.; Porcuna Fernández Monasterio, A. Creación de un sistema de información geográfico de ocupación del suelo en España. proyecto siose. In Proceedings of the Teledetección: Avances en la Observación de la Tierra, Florence, Italy, 21–23 September 2005.

65. SIOSE. Sistema de Ocupación del Suelo de España Especificaciones Técnicas SIOSE. Available online: https://www.siose.es/especificaciones-tecnicas (accessed on 27 August 2020).

66. Li, W.; MacBean, N.; Ciais, P.; Defourny, P.; Lamarche, C.; Bontemps, S.; Houghton, R.A.; Peng, S. Gross and net land cover changes in the main plant functional types derived from the annual ESA CCI land cover maps (1992–2015). *Earth Syst. Sci. Data* **2018**, *10*, 219–234. [CrossRef]

67. MacMillan, R.A.; Shary, P.A. Landforms and landform elements in geomorphometry. In *Developments in Soil Science*; Elsevier: Amsterdam, The Netherlands, 2009; Volume 33, pp. 227–254. ISBN 978-0-12-374345-9.

68. Jasiewicz, J.; Stepinski, T.F. Geomorphons—A pattern recognition approach to classification and mapping of landforms. *Geomorphology* **2013**, *182*, 147–156. [CrossRef]

69. Buitinck, L.; Louppe, G.; Blondel, M.; Pedregosa, F.; Mueller, A.; Grisel, O.; Niculae, V.; Prettenhofer, P.; Gramfort, A.; Grobler, J.; et al. API design for machine learning software: Experiences from the scikit-learn project. *arXiv* **2013**, arXiv:1309.0238.

70. Cohen, J. A coefficient of agreement for nominal scales. *Educ. Psychol. Measur.* **1960**, *20*, 37–46. [CrossRef]
71. Congalton, R.G. A review of assessing the accuracy of classifications of remotely sensed data. *Remote Sens. Environ.* **1991**, *37*, 35–46. [CrossRef]

72. Moreira, A.; Prats-Iraola, P.; Younis, M.; Krieger, G.; Hajnsek, I.; Papathanassiou, K.P. A tutorial on synthetic aperture radar. *IEEE Geosci. Remote Sens. Mag.* **2013**, *1*, 6–43. [CrossRef]

73. NASA. USGS Release of NASADEM Data Products. Available online: https://lpdaac.usgs.gov/news/release-nasadem-data-products/ (accessed on 11 September 2020).

<figure><img src="image_10.png" /></figure>

© 2020 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).