Mitigation of atmospheric perturbations and solid Earth movements in a TerraSAR-X time-series

The TerraSAR-X (TSX) synthetic aperture radar (SAR) marks the recent emergence of a new generation of spaceborne radar sensors that can for the first time lay claim to localization accuracies in the sub-meter range. The TSX platform's extremely high orbital stability and the sensor's hardware timing accuracy combine to enable direct measurements of atmospheric refraction and solid Earth movements. By modeling these effects for individual TSX acquisitions, absolute pixel geolocation accuracy on the order of several centimeters can be achieved without need for even a single tiepoint. A 16-month time series of images was obtained over a fixed test site, making it possible to validate both an atmospheric refraction and a solid Earth tide model, while at the same time establishing the instrument's long-term stability. These related goals were achieved by placing trihedral corner reflectors (CRs) at the test site and estimating their phase centers with centimeter-level accuracy using differential GPS (DGPS). Oriented in pairs toward a given satellite track, the CRs could be seen as bright "points” in the images, providing a geometric reference set. SAR images from the high-resolution spotlight (HS) mode were obtained in alternating ascending and descending orbit configurations. The highest-resolution products were selected for their small sample dimensions, as positions can be more precisely determined. Based on the delivered product annotations, the CR image positions were predicted, and these predictions were compared with their measured image positions both before and after compensation for atmospheric refraction and systematic solid Earth deviations. It was possible to show that when the atmospheric distortion and Earth tides are taken into account, the TSX HS products have geolocation accuracies far exceeding the specified requirements. Furthermore, this accuracy was maintained for the duration of the 16-month test period. It could be demonstrated that with a correctly calibrated sensor, and after accounting for atmospheric and tidal effects, tiepoint-free geolocation is possible with TSX with an absolute product accuracy of about 5c

Mitigation of atmospheric perturbations and solid Earth movements in a TerraSAR-X time-series

The TerraSAR-X (TSX) synthetic aperture radar (SAR) marks the recent emergence of a new generation of spaceborne radar sensors that can for the first time lay claim to localization accuracies in the sub-meter range. The TSX platform’s extremely high orbital stability and the sensor’s hardware timing accuracy combine to enable direct measurements of atmospheric refraction and solid Earth movements. By modeling these effects for individual TSX acquisitions, absolute pixel geolocation accuracy on the order of several centimeters can be achieved without need for even a single tiepoint. A 16-month time series of images was obtained over a fixed test site, making it possible to validate both an atmospheric refraction and a solid Earth tide model, while at the same time establishing the instrument’s long-term stability. These related goals were achieved by placing trihedral corner reflectors (CRs) at the test site and estimating their phase centers with centimeter-level accuracy using differential GPS (DGPS). Oriented in pairs toward a given satellite track, the CRs could be seen as bright “points” in the images, providing a geometric reference set. SAR images from the high-resolution spotlight (HS) mode were obtained in alternating ascending and descending orbit configurations. The highest-resolution products were selected for their small sample dimensions, as positions can be more precisely determined. Based on the delivered product annotations, the CR image positions were predicted, and these predictions were compared with their measured image positions both before and after compensation for atmospheric refraction and systematic solid Earth deviations. It was possible to show that when the atmospheric distortion and Earth tides are taken into account, the TSX HS products have geolocation accuracies far exceeding the specified requirements. Furthermore, this accuracy was maintained for the duration of the 16-month test period. It could be demonstrated that with a correctly calibrated sensor, and after accounting for atmospheric and tidal effects, tiepoint-free geolocation is possible with TSX with an absolute product accuracy of about 5 cm

The Extended Timing Annotation Dataset for Sentinel-1 - Product Description and First Evaluation Results

This article introduces the extended timing annotation dataset (ETAD) product for Sentinel-1 (S-1) which was developed in a joint effort of German Aerospace Center (DLR) and European Space Agency (ESA). It allows to correct range and azimuth timing of S-1 images for geophysical effects and for inaccuracies in synthetic aperture radar (SAR) image focusing. In combination with the precise orbit solution, these effects determine the absolute geolocation accuracy of S-1 SAR images and the relative collocation accuracy of repeat pass image stacks. ETAD contains the gridded timing corrections for the tropospheric and ionospheric path delays, the tidal-based surface displacements, and the SAR processing effects, all of which are computed for each data taken using standard models from geodesy and auxiliary atmospheric data. The ETAD product helps S-1 users to significantly improve the geolocation accuracy of the S-1 SAR products to better than 0.2 m and offers a potential solution for correcting large-scale interferometric phase variations. The product layout and product generation are described schematically. This article also reports first the results for different SAR techniques: first, the improvement in geolocation accuracy down to a few centimeters by verification of accurately surveyed corner reflector positions in the range–azimuth plane; second, the well-established offset-tracking technique, which is used for systematic ice velocity monitoring of ice sheets and glaciers, where ETAD can reduce velocity biases down to subcentimetric values; and third, the correction of atmospheric phase contributions in wide-area interferograms used for national and European ground motion services. These early results proof the added value of the ETAD corrections and that the product design is well-suited to be integrated into the processing flows of established SAR applications such as absolute ranging of targets, speckle/feature tracking, and interferometry

High-accuracy digital elevation model generation and ship monitoring from synthetic aperture radar images: innovative techniques and experimental results.

In this Thesis several state-of-the-art and innovative techniques for Digital Elevation Model (DEM) generation from Synthetic Aperture Radar (SAR) images are deeply analyzed, with a special focus on the methods which allow the improvement of the accuracy of the DEM product, which is directly related to the geolocation accuracy of geocoded images and is considered as an enabling factor for a large series of civilian and Defence applications. Furthermore, some of the proposed techniques, which are based both on phase and amplitude information, are experimented on real data, i.e. COSMO-SkyMed (CSK) data, assessing the achievable performances compared with the state-of-the-art, and pointing out and quantitatively highlighting the acquisition and processing strategies which would allow to maximize the quality of the results. Moreover, a critical analysis is performed about the main errors affecting the applied techniques, as well as the limitations of the orbital configurations, identifying several complementary techniques which would allow to overcome or mitigate the observed drawbacks. An innovative procedure for on-demand DEM production from CSK SAR data is elaborated and proposed, as well as an auto-validation technique which would enable the validation of the produced DEM also where vertical ground truths are not available. Based on the obtained results and on the consequent critical analysis, several interferometric specifications for new generation SAR satellites are identified. Finally, a literature review is proposed about the main state-of-the-art ship monitoring techniques, considered as one of the main fields of application which takes benefit from SAR data, based on single/multi-platform multi-channel SAR data, with a focus on TanDEM-X (TDX). In particular, in Chapter 1 the main concepts concerning SAR operating principles are introduced and the main characteristics and performances of CSK and TDX satellite systems are described; in Chapter 2 a review is proposed about the state-of-the-art SAR interferometric techniques for DEM generation, analyzing all the relevant processing steps and deepening the study of the main solutions recently proposed in the literature to increase the accuracy of the interferometric processing; in Chapter 3 complementary and innovative techniques respect to the interferometric processing are analyzed to mitigate disadvantages and to improve performances; in Chapter 4 experimental results are presented, obtained in the generation of high accuracy DEM by applying to a dataset of CSK images properly selected state-of-the-art interferometric techniques and innovative methods to improve DEM accuracy, exploring relevant limitations, and pointing out innovative acquisition and processing strategies. In Chapter 5, the basic principles of Ground Moving Target Indication (GMTI) are described, focusing on Displaced Phase Center Antenna (DPCA) and Along-Track Interferometry (ATI) techniques

Near real-time flood detection in urban and rural areas using high resolution Synthetic Aperture Radar images

A near real-time flood detection algorithm giving a synoptic overview of the extent of flooding in both urban and rural areas, and capable of working during night-time and day-time even if cloud was present, could be a useful tool for operational flood relief management. The paper describes an automatic algorithm using high resolution Synthetic Aperture Radar (SAR) satellite data that builds on existing approaches, including the use of image segmentation techniques prior to object classification to cope with the very large number of pixels in these scenes. Flood detection in urban areas is guided by the flood extent derived in adjacent rural areas. The algorithm assumes that high resolution topographic height data are available for at least the urban areas of the scene, in order that a SAR simulator may be used to estimate areas of radar shadow and layover. The algorithm proved capable of detecting flooding in rural areas using TerraSAR-X with good accuracy, classifying 89% of flooded pixels correctly, with an associated false positive rate of 6%. Of the urban water pixels visible to TerraSAR-X, 75% were correctly detected, with a false positive rate of 24%. If all urban water pixels were considered, including those in shadow and layover regions, these figures fell to 57% and 18% respectively

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

3D space intersection features extraction from Synthetic Aperture Radar images

The main purpose of this Thesis is to develop new theoretical models in order to extend the capabilities of SAR images space intersection techniques to generate three dimensional information. Furthermore, the study aims at acquiring new knowledge on SAR image interpretation through the three dimensional comprehension of the scene. The proposed methodologies allow to extend the known radargrammetric applications to vector data generation, exploiting SAR images acquired with every possible geometries. The considered geometries are points, circles, cylinders and lines. The study assesses the estimation accuracy of the features in terms of absolute and relative position and dimensions, analyzing the nowadays operational SAR sensors with a special focus on the national COSMO-SkyMed system. The proposed approach is original as it does not require the direct matching between homologous points of different images, which is a necessary step for the classical radargrammetric techniques; points belonging to the same feature, circular or linear, recognized in different images, are matched through specific models in order to estimate the dimensions and the location of the feature itself. This approach is robust with respect to the variation of the viewing angle of the input images and allows to better exploit archive data, acquired with diverse viewing geometries. The obtained results confirm the validity of the proposed theoretical approach and enable important applicative developments, especially in the Defence domain: (i) introducing original three dimensional measurement tools to support visual image interpretation; (ii) performing an advanced modelling of building counting only on SAR images; (iii) exploiting SAR images as a source for geospatial information and data; (iv) producing geospatial reference information, such as Ground Control Point, without any need for survey on the ground

On the use of COSMO/SkyMed data and Weather Models for interferometric DEM generation

AbstractThis work experiments the potentialities of COSMO/SkyMed (CSK) data in providing interferometric Digital Elevation Model (DEM). We processed a stack of CSK data for measuring with meter accuracy the ground elevation on the available coherent targets, and used these values to check the accuracy of DEMs derived from 5 tandem-like CSK pairs. In order to suppress the atmospheric signal we experimented a classical spatial filtering of the differential phase as well as the use of numerical weather prediction (NWP) model RAMS. Tandem-like pairs with normal baselines higher than 300 m allows to derive DEMs fulfilling the HRTI Level 3 specifications on the relative vertical accuracy, while the use of NWP models still seems unfeasible especially for X-band

Long-term monitoring of geodynamic surface deformation using SAR interferometry

Thesis (Ph.D.) University of Alaska Fairbanks, 2014Synthetic Aperture Radar Interferometry (InSAR) is a powerful tool to measure surface deformation and is well suited for surveying active volcanoes using historical and existing satellites. However, the value and applicability of InSAR for geodynamic monitoring problems is limited by the influence of temporal decorrelation and electromagnetic path delay variations in the atmosphere, both of which reduce the sensitivity and accuracy of the technique. The aim of this PhD thesis research is: how to optimize the quantity and quality of deformation signals extracted from InSAR stacks that contain only a low number of images in order to facilitate volcano monitoring and the study of their geophysical signatures. In particular, the focus is on methods of mitigating atmospheric artifacts in interferograms by combining time-series InSAR techniques and external atmospheric delay maps derived by Numerical Weather Prediction (NWP) models. In the first chapter of the thesis, the potential of the NWP Weather Research & Forecasting (WRF) model for InSAR data correction has been studied extensively. Forecasted atmospheric delays derived from operational High Resolution Rapid Refresh for the Alaska region (HRRRAK) products have been compared to radiosonding measurements in the first chapter. The result suggests that the HRRR-AK operational products are a good data source for correcting atmospheric delays in spaceborne geodetic radar observations, if the geophysical signal to be observed is larger than 20 mm. In the second chapter, an advanced method for integrating NWP products into the time series InSAR workflow is developed. The efficiency of the algorithm is tested via simulated data experiments, which demonstrate the method outperforms other more conventional methods. In Chapter 3, a geophysical case study is performed by applying the developed algorithm to the active volcanoes of Unimak Island Alaska (Westdahl, Fisher and Shishaldin) for long term volcano deformation monitoring. The volcano source location at Westdahl is determined to be approx. 7 km below sea level and approx. 3.5 km north of the Westdahl peak. This study demonstrates that Fisher caldera has had continuous subsidence over more than 10 years and there is no evident deformation signal around Shishaldin peak.Chapter 1. Performance of the High Resolution Atmospheric Model HRRR-AK for Correcting Geodetic Observations from Spaceborne Radars -- Chapter 2. Robust atmospheric filtering of InSAR data based on numerical weather prediction models -- Chapter 3. Subtle motion long term monitoring of Unimak Island from 2003 to 2010 by advanced time series SAR interferometry -- Chapter 4. Conclusion and future work

Interferometric Processing of TanDEM-X Images for Forest Height Estimation

Biomass is one of the most desired parameters for applications like climate modelling, resource assessment or wood industry. By using allometry equations (82) it is possible to obtain biomass information from canopy height. Some studies have demonstrated that current interferometric techniques applied to airborne Synthetic Aperture Radar (SAR) images can provide fairly accurate estimates of tree height (45, 52, 53, 54). Space based interferometric methods can provide global estimates of canopy height but they require very accurate orbit information. In this work the ability of the recently launched SAR satellites TerraSAR-X and TanDEM-X to estimate canopy height is evaluated.To do this, a complete interferometric processing chain is created including SAR data reading into memory, complex interferogram calculation, interferogram flattening by at Earth approximation and image transformation to geographical coordinates.Finally the resulting phase height maps are compared with a digital elevation model and a canopy height model of the terrain under study as well as with X-band E-SAR data from the FINSAR campaign (52, 53, 54) of the same area