Development of the TanDEM-X Calibration Concept: Analysis of Systematic Errors

The TanDEM-X mission, result of the partnership between the German Aerospace Center (DLR) and Astrium GmbH, opens a new era in spaceborne radar remote sensing. The first bistatic satellite synthetic aperture radar mission is formed by flying the TanDEM-X and TerraSAR-X in a closely controlled helix formation. The primary mission goal is the derivation of a high-precision global digital elevation model (DEM) according to High-Resolution Terrain Information (HRTI) level 3 accuracy. The finite precision of the baseline knowledge and uncompensated radar instrument drifts introduce errors that may compromise the height accuracy requirements. By means of a DEM calibration, which uses absolute height references, and the information provided by adjacent interferogram overlaps, these height errors can be minimized. This paper summarizes the exhaustive studies of the nature of the residual-error sources that have been carried out during the development of the DEM calibration concept. Models for these errors are set up and simulations of the resulting DEM height error for different scenarios provide the basis for the development of a successful DEM calibration strategy for the TanDEM-X mission

Fusion of Urban TanDEM-X raw DEMs using variational models

Recently, a new global Digital Elevation Model (DEM) with pixel spacing of 0.4 arcseconds and relative height accuracy finer than 2m for flat areas (slopes 20%) was created through the TanDEM-X mission. One important step of the chain of global DEM generation is to mosaic and fuse multiple raw DEM tiles to reach the target height accuracy. Currently, Weighted Averaging (WA) is applied as a fast and simple method for TanDEM-X raw DEM fusion in which the weights are computed from height error maps delivered from the Interferometric TanDEM-X Processor (ITP). However, evaluations show that WA is not the perfect DEM fusion method for urban areas especially in confrontation with edges such as building outlines. The main focus of this paper is to investigate more advanced variational approaches such as TV-L1 and Huber models. Furthermore, we also assess the performance of variational models for fusing raw DEMs produced from data takes with different baseline configurations and height of ambiguities. The results illustrate the high efficiency of variational models for TanDEM-X raw DEM fusion in comparison to WA. Using variational models could improve the DEM quality by up to 2m particularly in inner-city subsets.Comment: This is the pre-acceptance version, to read the final version, please go to IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing on IEEE Xplor

Elevation Extraction from Spaceborne SAR Tomography Using Multi-Baseline COSMO-SkyMed SAR Data

SAR tomography (TomoSAR) extends SAR interferometry (InSAR) to image a complex 3D scene with multiple scatterers within the same SAR cell. The phase calibration method and the super-resolution reconstruction method play a crucial role in 3D TomoSAR imaging from multi-baseline SAR stacks, and they both influence the accuracy of the 3D SAR tomographic imaging results. This paper presents a systematic processing method for 3D SAR tomography imaging. Moreover, with the newly released TanDEM-X 12 m DEM, this study proposes a new phase calibration method based on SAR InSAR and DEM error estimation with the super-resolution reconstruction compressive sensing (CS) method for 3D TomoSAR imaging using COSMO-SkyMed Spaceborne SAR data. The test, fieldwork, and results validation were executed at Zipingpu Dam, Dujiangyan, Sichuan, China. After processing, the 1 m resolution TomoSAR elevation extraction results were obtained. Against the terrestrial Lidar ‘truth’ data, the elevation results were shown to have an accuracy of 0.25 ± 1.04 m and a RMSE of 1.07 m in the dam area. The results and their subsequent validation demonstrate that the X band data using the CS method are not suitable for forest structure reconstruction, but are fit for purpose for the elevation extraction of manufactured facilities including buildings in the urban area

Elevation and Deformation Extraction from TomoSAR

3D SAR tomography (TomoSAR) and 4D SAR differential tomography (Diff-TomoSAR) exploit multi-baseline SAR data stacks to provide an essential innovation of SAR Interferometry for many applications, sensing complex scenes with multiple scatterers mapped into the same SAR pixel cell. However, these are still influenced by DEM uncertainty, temporal decorrelation, orbital, tropospheric and ionospheric phase distortion and height blurring. In this thesis, these techniques are explored. As part of this exploration, the systematic procedures for DEM generation, DEM quality assessment, DEM quality improvement and DEM applications are first studied. Besides, this thesis focuses on the whole cycle of systematic methods for 3D & 4D TomoSAR imaging for height and deformation retrieval, from the problem formation phase, through the development of methods to testing on real SAR data. After DEM generation introduction from spaceborne bistatic InSAR (TanDEM-X) and airborne photogrammetry (Bluesky), a new DEM co-registration method with line feature validation (river network line, ridgeline, valley line, crater boundary feature and so on) is developed and demonstrated to assist the study of a wide area DEM data quality. This DEM co-registration method aligns two DEMs irrespective of the linear distortion model, which improves the quality of DEM vertical comparison accuracy significantly and is suitable and helpful for DEM quality assessment. A systematic TomoSAR algorithm and method have been established, tested, analysed and demonstrated for various applications (urban buildings, bridges, dams) to achieve better 3D & 4D tomographic SAR imaging results. These include applying Cosmo-Skymed X band single-polarisation data over the Zipingpu dam, Dujiangyan, Sichuan, China, to map topography; and using ALOS L band data in the San Francisco Bay region to map urban building and bridge. A new ionospheric correction method based on the tile method employing IGS TEC data, a split-spectrum and an ionospheric model via least squares are developed to correct ionospheric distortion to improve the accuracy of 3D & 4D tomographic SAR imaging. Meanwhile, a pixel by pixel orbit baseline estimation method is developed to address the research gaps of baseline estimation for 3D & 4D spaceborne SAR tomography imaging. Moreover, a SAR tomography imaging algorithm and a differential tomography four-dimensional SAR imaging algorithm based on compressive sensing, SAR interferometry phase (InSAR) calibration reference to DEM with DEM error correction, a new phase error calibration and compensation algorithm, based on PS, SVD, PGA, weighted least squares and minimum entropy, are developed to obtain accurate 3D & 4D tomographic SAR imaging results. The new baseline estimation method and consequent TomoSAR processing results showed that an accurate baseline estimation is essential to build up the TomoSAR model. After baseline estimation, phase calibration experiments (via FFT and Capon method) indicate that a phase calibration step is indispensable for TomoSAR imaging, which eventually influences the inversion results. A super-resolution reconstruction CS based study demonstrates X band data with the CS method does not fit for forest reconstruction but works for reconstruction of large civil engineering structures such as dams and urban buildings. Meanwhile, the L band data with FFT, Capon and the CS method are shown to work for the reconstruction of large manmade structures (such as bridges) and urban buildings

Ocean Surface Observations Using the TanDEM-X Satellite Formation

The TanDEM-X SAR satellite formation permits improved ocean surface observations by means of bistatic along-track interferometry (ATI) when compared to single-satellite systems. The flexible imaging geometry of its two cooperating SAR sensors forms an interferometer that can achieve very high sensitivity to motions of objects on ground. This way, radar imaging of surface currents with unprecedented accuracy, high spatial resolution and wide coverage at the same time becomes possible. We demonstrate the capabilities of the sensor in the contexts of tidal current mapping, measurement of thermohaline and wind-driven ocean currents as well as detection of areas with surface films. We have developed a dedicated postprocessing system for TanDEM-X image products that allows extracting surface current information from the data. By this paper, we address bistatic data acquisition and processing aspects for sea surface imaging with TanDEM-X like interferometric baseline geometry, temporal decorrelation, and phase calibration. We present a variety of examples of data evaluation that clearly demonstrate the application potential of the methodology

An error prediction framework for interferometric SAR data

Three of the major error sources in interferometric synthetic aperture radar measurements of terrain elevation and displacement are baseline errors, atmospheric path length errors, and phase unwrapping errors. In many processing schemes, these errors are calibrated out by using ground control points (GCPs) (or an external digital elevation model). In this paper, a simple framework for the prediction of error standard deviation is outlined and investigated. Inputs are GCP position, a priori GCP accuracy, baseline calibration method along with a closed-form model for the covariance of atmospheric path length disturbances, and a model for phase unwrapping errors. The procedure can be implemented as a stand-alone add-on to standard interferometric processors. It is validated by using a set of single-frame interferograms acquired over Rome, Italy, and a double difference data set over Flevoland, The Netherlands

Multistatic SAR Imaging: First Results of a Four Phase Center Experiment with TerraSAR-X and TanDEM-X

Multichannel synthetic aperture radar (SAR) imaging offers the possibility to overcome the pulse repetition frequency (PRF) constraints inherent to single-channel SAR systems. The multichannel approach enables the acquisition of wide swathes with high azimuth resolution. Using a constellation or swarm of small satellites, a cost efficient, faulttolerant system can be envisaged. This paper describes the first results of a multistatic four phase center experiment conducted with TerraSAR-X and TanDEM-X. The experiment is intended to increase the know-how and demonstrate the capabilities for the design of future SAR systems. Key challenges are addressed, the experimental acquisition is described and an evaluation approach is presented. Finally, first results focusing on the azimuth ambiguity performance are shown

Multistatic SAR Imaging: Comparison of Simulation Results and Experimental Data

Synthetic aperture radar (SAR) systems in a multistatic configuration are a promising candidate for future Earth observation and reconnaissance radar systems. They overcome the sampling constraints inherent to single-channel SAR systems. Thus, a multistatic SAR system enables the acquisition of high-resolution images while maintaining wide-swath coverage. Employing several small satellites instead of a single large one, a cost-efficient system with graceful degradation characteristics can be envisaged. Additionally, such a constellation or swarm of sensors offers interferometric and tomographic capabilities, which a single-satellite system is not able to provide. This paper shows results of multistatic experiments obtained with TerraSAR-X and TanDEM-X and compares these results with theoretical simulations. The key parameters analyzed are the Doppler spectrum and the azimuth ambiguity suppression

Retrieval of Ocean Surface Currents and Winds Using Satellite SAR backscatter and Doppler frequency shift

Ocean surface winds and currents play an important role for weather, climate, marine life, ship navigation, oil spill drift and search and rescue. In-situ observations of the ocean are sparse and costly. Satellites provide a useful complement to these observations. Synthetic aperture radar (SAR) is particularly attractive due to its high spatial resolution and its capability to extract both sea surface winds and currents day and night and almost independent of weather.The work in this thesis involves processing of along-track interferometric SAR (ATI-SAR) data, analysis of the backscatter and Doppler frequency shift, and development of wind and current retrieval algorithms. Analysis of the Doppler frequency shift showed a systematic bias. A calibration method was proposed and implemented to correct for this bias. Doppler analysis also showed that the wave contribution to the SAR Doppler centroid often dominates over the current contribution. This wave contribution is estimated using existing theoretical and empirical Doppler models. For wind and current retrieval, two methods were developed and implemented.The first method, called the direct method, consists of retrieval of the wind speed from SAR backscatter using an empirical backscatter model. In order to retrieve the radial current, the retrieved wind speed is used to correct for the wave contribution. The current retrieval was assessed using two different (theoretical and empirical) Doppler models and wind inputs (model and SAR-derived). It was found that the results obtained by combining the Doppler empirical model with the SAR-derived wind speed were more consistent with ocean models.The second method, called Bayesian method, consists of blending the SAR observables (backscatter and Doppler shift) with an atmospheric and an oceanic model to retrieve the total wind and current vector fields. It was shown that this method yields more accurate estimates, i.e. reduces the models biases against in-situ measurements. Moreover, the method introduces small scale features, e.g. fronts and meandering, which are weakly resolved by the models.The correlation between the surface wind vectors and the SAR Doppler shift was demonstrated empirically using the Doppler shift estimated from over 300 TanDEM-X interferograms and ECMWF reanalysis wind vectors. Analysis of polarimetric data showed that theoretical models such as Bragg and composite surface models over-estimate the backscatter polarization ratio and Doppler shift polarization difference. A combination of a theoretical Doppler model and an empirical modulation transfer function was proposed. It was found that this model is more consistent with the analyzed data than the pure theoretical models.The results of this thesis will be useful for integrating SAR retrievals in ocean current products and assimilating SAR observables in the atmospheric, oceanic or coupled models. The results are also relevant for preparation studies of future satellite missions

Lava flow morphology at an erupting andesitic stratovolcano: a satellite perspective on El Reventador, Ecuador

Lava flows pose a significant hazard to infrastructure and property located close to volcanoes, and understanding how flows advance is necessary to manage volcanic hazard during eruptions. Compared to low-silica basaltic flows, flows of andesite composition are infrequently erupted and so relatively few studies of their characteristics and behaviour exist. We use El Reventador, Ecuador as a target to investigate andesitic lava flow properties during a 4.5 year period of extrusive eruption between February 2012 and August 2016. We use satellite radar to map the dimensions of 43 lava flows and look at variations in their emplacement behaviour over time. We find that flows descend the north and south flanks of El Reventador, and were mostly emplaced during durations shorter than the satellite repeat interval of 24 days.Flows ranged in length from 0.3 to 1.7 km, and the length of these flows decreased over the observation period. We measure a decrease in flow volume with time that is correlated with a long-term exponential decrease in eruption rate, and propose that this behaviour is caused by temporary magma storage in the conduit acting as a melt capacitor between the magma reservoir and the surface. We use the dimensions of the flow levees and widths to estimate the flow yield strengths, which were of the order of 10-100 kPa. We observe that some flows were diverted by topographic obstacles, and compare measurements of decreased channel width and increased flow thickness at the obstacles with observations from laboratory experiments. Radar observations, such as those presented here, could be used to map and measure properties of evolving lava flow fields at other remote or difficult to monitor volcanoes