InSAR Digital Surface Model Refinement by Block Adjustment with Horizontal Constraints

Abstract. In synthetic aperture radar (SAR) technology, urban mapping and modelling have become possible with revolutionary missions TerraSAR-X (TSX) and Cosmo-SkyMed (CSK) since 2007. These satellites offer 1m spatial resolution in high-resolution spotlight imaging mode and capable for high quality digital surface model (DSM) acquisition for urban areas utilizing interferometric SAR (InSAR) technology. With the advantage of independent generation from seasonal weather conditions, TSX and CSK DSMs are much in demand by scientific users. The performance of SAR DSMs is influenced by the distortions such as layover, foreshortening, shadow and double-bounce depend up on imaging geometry. In this study, the potential of DSMs derived from convenient 1m high-resolution spotlight (HS) InSAR pairs of CSK and TSX is validated by model-to-model absolute and relative accuracy estimations in an urban area. For the verification, an airborne laser scanning (ALS) DSM of the study area was used as the reference model. Results demonstrated that TSX and CSK urban DSMs are compatible in open, built-up and forest land forms with the absolute accuracy of 8–10 m. The relative accuracies based on the coherence of neighbouring pixels are superior to absolute accuracies both for CSK and TSX.

The development of satellite technology is rapidly increasing the evolution of remote sensing. Satellite images give extensive useful information about the land structure that is easily manageable in the process of generating true, high-speed information which allows the forecasting of future environmental and urban planning. Remote sensing comprises active and passive systems. Passive sensors detect natural radiation that is emitted or reflected by the object or surrounding area being observed. Active systems which produce their own electromagnetic energy and their main properties are their ability of collecting data in nearly all atmospheric conditions, day or night. These systems are frequently used to generate a digital elevation model (DEM) because they cover large areas. DEM supplies essential data for applications that are concerned with the Earth’s surface and DEMs derived from survey data are accurate but very expensive and time consuming to create. However, the use of satellite remote sensing to provide images to generate a DEM is considered to be an efficient method of obtaining data. Interferometric Synthetic Aperture Radar (InSAR) is a new geodetic technique for determining earth topography. InSAR measurements are highly dense and they only give information in Line of Sight of Radar. In the study, interferograms were produced from the InSAR images taken by ERS satellites in 1992 and 2007 and we developed the methods to generate a DEM using the InSAR technique and present the results relating to Kayseri Province in Turkey. The accuracy of the DEM derived from the InSAR technique is evaluated in comparison with a reference DEM generated from contours in a topographical map.

Voids caused by shadow, layover, and decorrelation usually occur in digital elevation models (DEMs) of mountainous areas that are derived from interferometric synthetic aperture radar (InSAR) datasets. The presence of voids degrades the quality and usability of the DEMs. Thus, void removal is considered as an integral part of the DEM production using InSAR data. The fusion of multiple DEMs has been widely recognized as a promising way for the void removal. Because the vertical accuracy of multiple DEMs can be different, the selection of optimum weights becomes a key problem in the fusion and is studied in this article. As a showcase, two high-resolution InSAR DEMs near Mt. Qilian in northwest China are created and then merged. The two pairs of InSAR data were acquired by TerraSAR-X from an ascending orbit and COSMO-SkyMed from a descending orbit. A maximum likelihood fusion scheme with the weights optimally determined by the height of ambiguity and the variance of phase noise is adopted to syncretize the two DEMs in our study. The fused DEM has a fine spatial resolution of 10 m and depicts the landform of the study area well. The percentage of void cells in the fused DEM is only 0.13 %, while 6.9 and 5.7 % of the cells in the COSMO-SkyMed DEM and the TerraSAR-X DEM are originally voids. Using the ICESat/GLAS elevation data and the Chinese national DEM of scale 1:50,000 as references, we evaluate vertical accuracy levels of the fused DEM as well as the original InSAR DEMs. The results show that substantial improvements could be achieved by DEM fusion after atmospheric phase screen removal. The quality of fused DEM can even meet the high-resolution terrain information (HRTI) standard.

Interferometric synthetic aperture radar (InSAR) has become a key technology for producing high-precision digital surface models (DSMs) and digital orthophoto maps (DOMs) in full time and all weathers. Airborne millimeter-wave InSAR, with large-scale and high-resolution imaging, is characterized by high spatial resolution, flexibility, and immunity to loss-of-correlation. This paper introduces our modeling experiments with airborne dual-antenna, Ka-band InSAR regarding typical topographies of China. Ka-band SAR data were acquired in designated experimental areas in flat (Heyang area in Shaanxi) and mountainous areas (Shibing area in Guizhou and Qionglai area in Sichuan). The key processing of the experimental data for DSMs and DOMs is demonstrated in the paper, especially the proposed robust and efficient phase unwrapping (PU) method for the interferometric data and block adjustment method of strip calibration. The results show that the proposed unwrapping method can provide reliable unwrapped phase results in undulating areas, and the block adjustment can carry out consistent calibration for strips with sparse ground control points (GCPs). The accuracy assessment of the DSM shows that the coordinate root mean square error (RSME) of the obtained DSM is less than 2 m in height, and 2.5 m horizontally, which meets the 1:5000 requirement for topographic mapping in difficult areas.

Interferometric synthetic aperture radar (InSAR) processing, if lacking high quality ground control points (GCPs), may produce large errors in the final DEM. Some kind of alignment or registration can reduce these errors. This paper evaluates the accuracy of InSAR processed digital elevation models (DEM), against a high resolution DEM. InSAR DEMs were aligned and least-squares registered with USGS, 1/3 arc second NED. The accuracy evaluations after registration show InSAR DEMs alignment and registration can eliminate DEM errors caused by lack of accurate GCPs. In some cases, this process improves the accuracy of the InSAR DEM even more efficiently than including GCPs into InSAR processing.

Comparison and Validation of Different DEM Data Derived from InSAR

GaoFen-3 (GF-3) is the first fully polarimetric C-band synthetic aperture radar (SAR) satellite in China. Interferometric observation is an important capability of GF-3. However, the long temporal baseline is a great challenge for repeat-pass interferometric synthetic aperture radar (InSAR) processing, especially over mountainous areas. Phase noise and complex topography make phase unwrapping difficult. This paper proposes a step-wise approach for high-resolution InSAR data processing focusing on facilitating phase unwrapping. First, the coarse GF-3 digital elevation model (DEM) is reconstructed by InSAR processing with multilooking of large looks, which is used to suppress phase noise. Second, the coarse GF-3 DEM is used to remove the terrain phase from the InSAR phase with multilooking of small looks, which can reduce phase gradient and improve the accuracy of phase unwrapping. After several iterations, the GF-3 DEM is finally inversed with multilooking of 2 looks, which is of high resolution and high accuracy. The performance of GF-3 C-band InSAR data on the extraction of DEM is examined over a deciduous forested area with the proposed method. The root-mean-square error of the extracted GF-3 DEM is 3.78 m in comparison with the LiDAR DEM.

In this paper, the study on the digital elevation model (DEM) reconstruction using Gaofen-3 satellite repeat-pass synthetic aperture radar (SAR) interferometry is presented. The coherence of the single-look complex (SLC) images and residuals of interferogram provide a powerful guarantee for interferometric processing for Gaofen-3. Baseline estimation is the key factor to generate accurate DEM with SAR interferometry. This paper uses the coarse orbit data and the precise orbit data of Gaofen-3 to carry on the baseline estimation and uses the prior coarse DEM information to refine the baseline estimation. Then, the topographic height is reconstructed with more precise spatial baseline. An experiment is conducted to verify the conclusion, and some mountain and flatten region are chosen as the study area. The results show that the DEM errors based on the prior DEM method are rapidly converged compared with the ones based on the orbit data method. Therefore, benefiting from the precise orbit data of the Gaofen-3 satellite, the accuracy of the baseline estimation can be improved further by using the method based on prior DEM information; thus, higher accuracy of DEM can be achieved so that Gaofen-3 repeat-pass data have some good characteristics, such as high imaging quality, stable spatial baseline, and high correlation, which has a good repeat-pass SAR interferometry application prospect.

Trong những nam gần đây, kỹ thuật SAR giao thoa (Synthetic Aperture Radar Interferometry - InSAR) đa được xem như một trong những kỹ thuật hiệu quả hỗ trợ cho việc xây dựng va cập nhật dữ liệu của GIS. Dữ liệu độ cao số trong GIS thường lưu trữ theo cac dạng mo hinh độ cao số (DEM) va được sử dụng kha phổ biến để thể hiện địa hinh. Bai bao nhằm giới thiệu khả nang ứng dụng kỹ thuật InSAR trong xây dựng DEM va minh họa kết qủa thực nghiệm đạt được tại lưu vực song Kon – Ha Thanh thuộc tỉnh Binh Định. Độ chinh xac của mo hinh độ cao số (DEM) tạo bởi kỹ thuật InSAR được so sanh với cac phương phap tạo DEM khac (DEM nhận được từ cặp ảnh lập thể SPOT hay ASTER). In recent years, the Synthetic Aperture Radar Interferometry (InSAR) technique has been considered as an efficienty method in creating or updating GIS data. Digital Elevation Data referred to as Digital Elevation Models (DEMs) are digital representations of the sharpe of the earth’s surface. DEMs derived from survey data are accurate but they are very expensive and time-consuming. This paper introduces the ability of application of InSAR technique to generate a DEM and shows the result of experiment in Kon – Ha Thanh river in Binh Dinh province. The accuracy of DEM derived from InSAR technique is evaluated in comparison with other methods (DEM is derived from the stereoscopic images of SPOT or ASTER).

Since 2007, the revolutionary X-band synthetic aperture radar (SAR) satellites Cosmo-SkyMed (CSK) and TerraSAR-X (TSX) have been collecting high-resolution images that are convenient for digital surface model (DSM) acquisition by the interferometric SAR (InSAR) technique. In this study, the potential of DSMs derived from CSK and TSX high-resolution spotlight (HS) images is thoroughly analysed regarding the basic accuracy metrics absolute vertical accuracy (AVA) and relative vertical accuracy (RVA). Utilising convenient InSAR pairs, 2 m gridded DSMs are generated in Istanbul and validated with model-based comparisons using an actual airborne laser scanning (ALS) DSM. Results show that CSK and TSX DSMs are compatible in open, built-up and forest land forms. The AVAs are between 8 m and 10 m based on standard deviation of height discrepancies against the ALS model. The RVAs, calculated by the coherence of neighbouring pixels for each DSM, are superior to AVAs for both CSK and TSX.

In Antarctica, SAR interferometry has largely been used in coastal glacial areas, while in rare cases this method has been used on the Antarctic plateau. In this paper, the authors present a digital elevation and ice flow map based on SAR interferometry for an area encompassing Talos Dome (TD) and the Frontier Mountain (FM) meteorite site in North Victoria Land/Antarctica. A digital elevation model (DEM) was calculated using a double SAR interferometry method. The DEM of the region was calculated by extracting approximately 100 control points from the Reference Elevation Model of Antarctica (REMA). The two DEMs differ slightly in some areas, probably due to the penetration of the SAR-C band signal into the cold firn. The largest differences are found in the western area of TD, where the radar penetration is more pronounced and fits well with the layer structures calculated by the georadar and the snow accumulation observations. By differentiating a 70-day interferogram with the calculated DEM, a displacement interferogram was calculated that represents the ice dynamics. The resulting ice flow pattern clearly shows the catchment areas of the Priestley and Rennick Glaciers as well as the ice flow from the west towards Wilkes Basin. The ice velocity field was analysed in the area of FM. This area has become well known due to the search for meteorites. The velocity field in combination with the calculated DEM confirms the generally accepted theories about the accumulation of meteorites over the Antarctic Plateau.

Abstract. The Antarctic Peninsula (AP) is one of the widely studied polar regions because of its sensitivity to climate change and potential contribution of its glaciers to global sea level rise. Precise digital elevation models (DEMs) at a high spatial resolution are much demanded for investigating the complex glacier system of the AP at fine scales. However, the two most recent circum-Antarctic DEMs, the 12 m TanDEM-X DEM (TDM DEM) from bistatic interferometric synthetic aperture radar (InSAR) data acquired between 2013 and 2014 and the Reference Elevation Model of Antarctica mosaic (REMA mosaic) at an 8 m spatial resolution derived from optical data acquired between 2011 and 2017 have specific individual limitations in this area. The TDM DEM has the advantage of good data consistency and few data voids (approx. 0.85 %), but there exist residual systematic elevation errors such as phase-unwrapping errors in the non-edited DEM version. The REMA mosaic has high absolute vertical accuracy, but on the AP it suffers from extended areas with data voids (approx. 8 %). To generate a consistent, gapless and high-resolution topography product of the AP, we fill the data voids in the TDM DEM with newly processed TDM raw DEM data acquired in austral winters of 2013 and 2014 and detect and correct the residual systematic elevation errors (i.e., elevation biases) in the TDM DEM with the support of the accurately calibrated REMA mosaic. Instead of a pixelwise replacement with REMA mosaic elevations, these provide reference values to correct the TDM elevation biases over entire regions detected through a path propagation algorithm. The procedure is applied iteratively to gradually correct the errors in the TDM DEM from a large to small scale. The proposed method maintains the characteristics of an InSAR-generated DEM and is minimally influenced by temporal or penetration differences between the TDM DEM and REMA mosaic. The performance of the correction is evaluated with laser altimetry data from Operation IceBridge and ICESat-2 missions. The overall root mean square error (RMSE) of the corrected TDM DEM has been reduced from more than 30 m to about 10 m which together with the improved absolute elevation accuracy indicates comparable values to the REMA mosaic. The generated high-resolution DEM depicts the up-to-date topography of the AP in detail and can be widely used for interferometric applications as well as for glaciological studies on individual glaciers or at regional scales.

Slow-moving landslides can trigger severe disasters when activated by earthquakes, torrential rains, or typhoons. Early detection is crucial for mitigating loss of life and property damage. Interferometric Synthetic Aperture Radar (InSAR) technology is among the most effective techniques for detecting slow-moving landslides, though its accuracy can be further improved through integration with optical imagery and Digital Elevation Models (DEM). Current machine learning approaches that combine InSAR and optical data suffer from limited efficiency, poor transferability, and challenges in regional-scale application. To address these limitations, this study proposes a multimodal dual-path network that integrates InSAR products with textural information from optical imagery to detect slow-moving landslides. One path processes InSAR deformation rates and topographic factors, while the other incorporates texture information and auxiliary data. Together, these paths extract semantic information from high-dimensional spatial features and condense it into low-dimensional representations. A pyramid pooling module is employed to capture multi-scale features during low-level semantic extraction. For feature fusion, a rate-constrained adaptive module is introduced to enhance the contribution of deformation rates to slow-moving landslides. According to the results, the proposed method improves the F1-score for landslide detection by 6% compared to using InSAR products alone. These results provide reliable technical support for regional landslide inventory compilation and disaster management, as well as new insights for regional-scale surveys in slow-moving landslide-prone areas.

Urban flood models that use Digital Elevation Models (DEMs) to simulate extent and depth of flood inundation rely on the accuracy of DEMs for predicting flood events. Despite recent advances in developing vegetation corrected DEMs, the effect of building height and density errors in global DEMs in urban areas are still poorly understood, and their correction remains a challenge. In this research we developed a methodology for building error correction that can be applied to any other case study, where building density data and a local reference DEM data are available. This methodology was applied to Nairobi, Kenya using six global DEMs (SRTM, MERIT, ALOS, NASADEM, TanDEM-X 12 m, and TanDEM-X 90 m DEM). Our results show building error at highest building density varying between 1.25 m and 5.07 m for the DEMs used, with the MERIT DEM showing the smallest vertical height error from the reference DEM. The six DEMs were corrected by deriving a linear relationship between building density and DEM error. Our findings show that the removal of building density error resulted in the improvement of the vertical height accuracy of the global DEMs of up to 45% for MERIT and 40% for ALOS. This methodology was also applied to the Central Business District (CBD) area of Nairobi, characterized by taller buildings and high building density. The error parameters in the CBD area resulted to be between 15 to 45% higher than those of the Nairobi city wide area for the six global DEMs, thus providing further insights into the contribution of building heights to errors in global DEMs. Building height data is still unavailable on a global scale and our results show that global DEMs can be usefully corrected for building density errors in urban areas, even where specific building height data are not available.

The Digital Elevation Models (DEMs) are an important source of topographical data for many scientific and engineering applications. Where topographical data are unavailable, global coverage elevation data sets, typically DEMs from remotely sensed data, are the main sources of such information. Interferometric SAR (InSAR), is a useful method for low-cost, relatively precise and wide-coverage surface DEM generation. However, ground deformation should somehow be excluded from InSAR-based DEM generation. To identify surface deformation areas, the so-called Differential InSAR (DInSAR) is a commonly used method. In this paper, the authors propose a two-step DEM generation method: the ground deformation area detection using DInSAR technique and deformation area exclusion in InSAR DEM generation by detected mask.