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

Seafloor depth estimation by means of interferometric synthetic aperture sonar

The topic of this thesis is relative depth estimation using interferometric sidelooking sonar. We give a thorough description of the geometry of interferometric sonar and of time delay estimation techniques. We present a novel solution for the depth estimate using sidelooking sonar, and review the cross-correlation function, the cross-uncertainty function and the phase-differencing technique. We find an elegant solution to co-registration and unwrapping by interpolating the sonar data in ground-range. Two depth estimation techniques are developed: Cross-correlation based sidescan bathymetry and synthetic aperture sonar (SAS) interferometry. We define flank length as a measure of the horizontal resolution in bathymetric maps and find that both sidescan bathymetry and SAS interferometry achieve theoretical resolutions. The vertical precision of our two methods are close to the performance predicted from the measured coherence. We study absolute phase-difference estimation using bandwidth and find a very simple split-bandwidth approach which outperforms a standard 2D phase unwrapper on complicated objects. We also examine advanced filtering of depth maps. Finally, we present pipeline surveying as an example application of interferometric SAS

UAVSAR: A new NASA airborne SAR system for science and technology research

NASA’s Jet Propulsion Laboratory is currently building a reconfigurable, polarimetric L-band synthetic aperture radar (SAR), specifically designed to acquire airborne repeat track SAR data for differential interferometric measurements. Differential interferometry can provide key deformation measurements, important for studies of earthquakes, volcanoes and other dynamically changing phenomena. Using precision real-time GPS and a sensor controlled flight management system, the system will be able to fly predefined paths with great precision. The expected performance of the flight control system will constrain the flight path to be within a 10 m diameter tube about the desired flight track. The radar will be designed to be operable on a UAV (Unpiloted Arial Vehicle) but will initially be demonstrated on a on a NASA Gulfstream III. The radar will be fully polarimetric, with a range bandwidth of 80 MHz (2 m range resolution), and will support a 16 km range swath. The antenna will be electronically steered along track to assure that the antenna beam can be directed independently, regardless of the wind direction and speed. Other features supported by the antenna include elevation monopulse and pulse-to-pulse re-steering capabilities that will enable some novel modes of operation. The system will nominally operate at 45,000 ft (13800 m). The program began as an Instrument Incubator Project (IIP) funded by NASA Earth Science and Technology Office (ESTO)

The European Space Agency BIOMASS mission: Measuring forest above-ground biomass from space

The primary objective of the European Space Agency's 7th Earth Explorer mission, BIOMASS, is to determine the worldwide distribution of forest above-ground biomass (AGB) in order to reduce the major uncertainties in calculations of carbon stocks and fluxes associated with the terrestrial biosphere, including carbon fluxes associated with Land Use Change, forest degradation and forest regrowth. To meet this objective it will carry, for the first time in space, a fully polarimetric P-band synthetic aperture radar (SAR). Three main products will be provided: global maps of both AGB and forest height, with a spatial resolution of 200 m, and maps of severe forest disturbance at 50 m resolution (where “global” is to be understood as subject to Space Object tracking radar restrictions). After launch in 2022, there will be a 3-month commissioning phase, followed by a 14-month phase during which there will be global coverage by SAR tomography. In the succeeding interferometric phase, global polarimetric interferometry Pol-InSAR coverage will be achieved every 7 months up to the end of the 5-year mission. Both Pol-InSAR and TomoSAR will be used to eliminate scattering from the ground (both direct and double bounce backscatter) in forests. In dense tropical forests AGB can then be estimated from the remaining volume scattering using non-linear inversion of a backscattering model. Airborne campaigns in the tropics also indicate that AGB is highly correlated with the backscatter from around 30 m above the ground, as measured by tomography. In contrast, double bounce scattering appears to carry important information about the AGB of boreal forests, so ground cancellation may not be appropriate and the best approach for such forests remains to be finalized. Several methods to exploit these new data in carbon cycle calculations have already been demonstrated. In addition, major mutual gains will be made by combining BIOMASS data with data from other missions that will measure forest biomass, structure, height and change, including the NASA Global Ecosystem Dynamics Investigation lidar deployed on the International Space Station after its launch in December 2018, and the NASA-ISRO NISAR L- and S-band SAR, due for launch in 2022. More generally, space-based measurements of biomass are a core component of a carbon cycle observation and modelling strategy developed by the Group on Earth Observations. Secondary objectives of the mission include imaging of sub-surface geological structures in arid environments, generation of a true Digital Terrain Model without biases caused by forest cover, and measurement of glacier and icesheet velocities. In addition, the operations needed for ionospheric correction of the data will allow very sensitive estimates of ionospheric Total Electron Content and its changes along the dawn-dusk orbit of the mission

Investigation of developments in interferometric synthetic aperture radar until 1994

Bibliography: p. 149-155.This thesis examines the topic of Synthetic Aperture Radar Interferometry in a historical perspective, tracing its development from its beginnings in the 1960s up until May 1994. Applications are listed and airborne and spaceborne implementations reviewed. The underlying theory of interferometry is explained, including a discussion of error sources, and a simulation for point targets is documented to illustrate the interferometric processing steps. The application of the SASAR VHF SAR system to interferometric operation is examined analytically

Research progress on geosynchronous synthetic aperture radar

Based on its ability to obtain two-dimensional (2D) high-resolution images in all-time and all-weather conditions, spaceborne synthetic aperture radar (SAR) has become an important remote sensing technique and the study of such systems has entered a period of vigorous development. Advanced imaging modes such as radar interferometry, tomography, and multi-static imaging, have been demonstrated. However, current in-orbit spaceborne SARs, which all operate in low Earth orbits, have relatively long revisit times ranging from several days to dozens of days, restricting their temporal sampling rate. Geosynchronous SAR (GEO SAR) is an active research area because it provides significant new capability, especially its much-improved temporal sampling. This paper reviews the research progress of GEO SAR technologies in detail. Two typical orbit schemes are presented, followed by the corresponding key issues, including system design, echo focusing, main disturbance factors, repeat-track interferometry, etc, inherent to these schemes. Both analysis and solution research of the above key issues are described. GEO SAR concepts involving multiple platforms are described, including the GEO SAR constellation, GEO-LEO/airborne/unmanned aerial vehicle bistatic SAR, and formation flying GEO SAR (FF-GEO SAR). Due to the high potential of FF-GEO SAR for three-dimensional (3D) deformation retrieval and coherence-based SAR tomography (TomoSAR), we have recently carried out some research related to FF-GEO SAR. This research, which is also discussed in this paper, includes developing a formation design method and an improved TomoSAR processing algorithm. It is found that GEO SAR will continue to be an active topic in the aspect of data processing and multi-platform concept in the near future

Coherent Change Detection Under a Forest Canopy

Coherent change detection (CCD) is an established technique for remotely monitoring landscapes with minimal vegetation or buildings. By evaluating the local complex correlation between a pair of synthetic aperture radar (SAR) images acquired on repeat passes of an airborne or spaceborne imaging radar system, a map of the scene coherence is obtained. Subtle disturbances of the ground are detected as areas of low coherence in the surface clutter. This thesis investigates extending CCD to monitor the ground in a forest. It is formulated as a multichannel dual-layer coherence estimation problem, where the coherence of scattering from the ground is estimated after suppressing interference from the canopy by vertically beamforming multiple image channels acquired at slightly different grazing angles on each pass. This 3D SAR beamforming must preserve the phase of the ground response. The choice of operating wavelength is considered in terms of the trade-off between foliage penetration and change sensitivity. A framework for comparing the performance of different radar designs and beamforming algorithms, as well as assessing the sensitivity to error, is built around the random-volume-over-ground (RVOG) model of forest scattering. If the ground and volume scattering contributions in the received echo are of similar strength, it is shown that an L-band array of just three channels can provide enough volume attenuation to permit reasonable estimation of the ground coherence. The proposed method is demonstrated using an RVOG clutter simulation and a modified version of the physics-based SAR image simulator PolSARproSim. Receiver operating characteristics show that whilst ordinary single-channel CCD is unusable when a canopy is present, 3D SAR CCD permits reasonable detection performance. A novel polarimetric filtering algorithm is also proposed to remove contributions from the ground-trunk double-bounce scattering mechanism, which may mask changes on the ground near trees. To enable this kind of polarimetric processing, fully polarimetric data must be acquired and calibrated. Motivated by an interim version of the Ingara airborne imaging radar, which used a pair of helical antennas to acquire circularly polarised data, techniques for the estimation of polarimetric distortion in the circular basis are investigated. It is shown that the standard approach to estimating cross-talk in the linear basis, whereby expressions for the distortion of reflection-symmetric clutter are linearised and solved, cannot be adapted to the circular basis, because the first-order effects of individual cross-talk parameters cannot be distinguished. An alternative approach is proposed that uses ordinary and gridded trihedral corner reflectors, and optionally dihedrals, to iteratively estimate the channel imbalance and cross-talk parameters. Monte Carlo simulations show that the method reliably converges to the true parameter values. Ingara data is calibrated using the method, with broadly consistent parameter estimates obtained across flights. Genuine scene changes may be masked by coherence loss that arises when the bands of spatial frequencies supported by the two passes do not match. Trimming the spatial-frequency bands to their common area of support would remove these uncorrelated contributions, but the bands, and therefore the required trim, depend on the effective collection geometry at each pixel position. The precise dependence on local slope and collection geometry is derived in this thesis. Standard methods of SAR image formation use a flat focal plane and allow only a single global trim, which leads to spatially varying coherence loss when the terrain is undulating. An image-formation algorithm is detailed that exploits the flexibility offered by back-projection not only to focus the image onto a surface matched to the scene topography but also to allow spatially adaptive trimming. Improved coherence is demonstrated in simulation and using data from two airborne radar systems.Thesis (Ph.D.) -- University of Adelaide, School of Electrical & Electronic Engineering, 202

Joint multi-baseline SAR interferometry

We propose a technique to provide interferometry by combining multiple images of the same area. This technique differs from the multi-baseline approach in literature as (a) it exploits all the images simultaneously, (b) it performs a spectral shift preprocessing to remove most of the decorrelation, and (c) it exploits distributed targets. The technique is mainly intended for DEM generation at centimetric accuracy, as well as for differential interferometry. The problem is framed in the contest of single-input multiple-output (SIMO) channel estimation via the cross-relations (CR) technique and the resulting algorithm provides significant improvements with respect to conventional approaches based either on independent analysis of single interferograms or multi-baselines phase analysis of single pixels of current literature, for those targets that are correlated in all the images, like for long-term coherent areas, or for acquisitions taken with a short revisit time (as those gathered with future satellite constellations)

Interferometric Synthetic Aperture Sonar Signal Processing for Autonomous Underwater Vehicles Operating Shallow Water

The goal of the research was to develop best practices for image signal processing method for InSAS systems for bathymetric height determination. Improvements over existing techniques comes from the fusion of Chirp-Scaling a phase preserving beamforming techniques to form a SAS image, an interferometric Vernier method to unwrap the phase; and confirming the direction of arrival with the MUltiple SIgnal Channel (MUSIC) estimation technique. The fusion of Chirp-Scaling, Vernier, and MUSIC lead to the stability in the bathymetric height measurement, and improvements in resolution. This method is computationally faster, and used less memory then existing techniques