Minimizing the residual topography effect on interferograms to improve DInSAR results: estimating land subsidence in Port-Said City, Egypt

The accurate detection of land subsidence rates in urban areas is important to identify damage-prone areas and provide decision-makers with useful information. Meanwhile, no precise measurements of land subsidence have been undertaken within the coastal Port-Said City in Egypt to evaluate its hazard in relationship to sea-level rise. In order to address this shortcoming, this work introduces and evaluates a methodology that substantially improves small subsidence rate estimations in an urban setting. Eight ALOS/PALSAR-1 scenes were used to estimate the land subsidence rates in Port-Said City, using the Small BAse line Subset (SBAS) DInSAR technique. A stereo pair of ALOS/PRISM was used to generate an accurate DEM to minimize the residual topography effect on the generated interferograms. A total of 347 well distributed ground control points (GCP) were collected in Port-Said City using the leveling instrument to calibrate the generated DEM. Moreover, the eight PALSAR scenes were co-registered using 50 well-distributed GCPs and used to generate 22 interferogram pairs. These PALSAR interferograms were subsequently filtered and used together with the coherence data to calculate the phase unwrapping. The phase-unwrapped interferogram-pairs were then evaluated to discard four interferograms that were affected by phase jumps and phase ramps. Results confirmed that using an accurate DEM (ALOS/PRISM) was essential for accurately detecting small deformations. The vertical displacement rate during the investigated period (2007–2010) was estimated to be −28 mm. The results further indicate that the northern area of Port-Said City has been subjected to higher land subsidence rates compared to the southern area. Such land subsidence rates might induce significant environmental changes with respect to sea-level rise

Ameliorative Minimum Cost Flow Algorithm for Phase Unwrapping

AbstractAmeliorative minimum cost flow algorithm is put forward in this paper Based on MCF,. As we all know, MCF is one phase unwrapping algorithm based on network flow and it is used widely because of many merits. However, efficiency of MCF is relative low and it demands upper computer capability during phase unwrapping. The ameliorative minimum cost flow algorithm is that interferograms are divided into skit sub-areas so that it can be unwrapped respectively, after this, interferograms are fused by supper wavelet based on contourlet transform. In the end, Minimum Cost Flow and its Ameliorative algorithm are validated by practical example, results indicate Ameliorative algorithm are effective phase unwrapping algorithms

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

Advanced satellite radar interferometry for small-scale surface deformation detection

Synthetic aperture radar interferometry (InSAR) is a technique that enables generation of Digital Elevation Models (DEMs) and detection of surface motion at the centimetre level using radar signals transmitted from a satellite or an aeroplane. Deformation observations can be performed due to the fact that surface motion, caused by natural and human activities, generates a local phase shift in the resultant interferogram. The magnitude of surface deformation can be estimated directly as a fraction of the wavelength of the transmitted signal. Moreover, differential InSAR (DInSAR) eliminates the phase signal caused by relief to yield a differential interferogram in which the signature of surface deformation can be seen. Although InSAR applications are well established, the improvement of the interferometry technique and the quality of its products is highly desirable to further enhance its capabilities. The application of InSAR encounters problems due to noise in the interferometric phase measurement, caused by a number of decorrelation factors. In addition, the interferogram contains biases owing to satellite orbit errors and atmospheric heterogeneity These factors dramatically reduce the stlectiveness of radar interferometry in many applications, and, in particular, compromise detection and analysis of small-scale spatial deformations. The research presented in this thesis aim to apply radar interferometry processing to detect small-scale surface deformations, improve the quality of the interferometry products, determine the minimum and maximum detectable deformation gradient and enhance the analysis of the interferometric phase image. The quality of DEM and displacement maps can be improved by various methods at different processing levels. One of the methods is filtering of the interferometric phase.However, while filtering reduces noise in the interferogram, it does not necessarily enhance or recover the signal. Furthermore, the impact of the filter can significantly change the structure of the interferogram. A new adaptive radar interferogram filter has been developed and is presented herein. The filter is based on a modification to the Goldstein radar interferogram filter making the filter parameter dependent on coherence so that incoherent areas are filtered more than coherent areas. This modification minimises the loss of signal while still reducing the level of noise. A methodology leading to the creation of a functional model for determining minimum and maximum detectable deformation gradient, in terms of the coherence value, has been developed. The sets of representative deformation models have been simulated and the associated phase from these models has been introduced to real SAR data acquired by ERS-1/2 satellites. A number of cases of surface motion with varying magnitudes and spatial extent have been simulated. In each case, the resultant surface deformation has been compared with the 'true' surface deformation as defined by the deformation model. Based on those observations, the functional model has been developed. Finally, the extended analysis of the interferometric phase image using a wavelet approach is presented. The ability of a continuous wavelet transform to reveal the content of the wrapped phase interferogram, such as (i) discontinuities, (ii) extent of the deformation signal, and (iii) the magnitude of the deformation signal is examined. The results presented represent a preliminary study revealing the wavelet method as a promising technique for interferometric phase image analysis

An Unsupervised Generative Neural Approach for InSAR Phase Filtering and Coherence Estimation

Phase filtering and pixel quality (coherence) estimation is critical in producing Digital Elevation Models (DEMs) from Interferometric Synthetic Aperture Radar (InSAR) images, as it removes spatial inconsistencies (residues) and immensely improves the subsequent unwrapping. Large amount of InSAR data facilitates Wide Area Monitoring (WAM) over geographical regions. Advances in parallel computing have accelerated Convolutional Neural Networks (CNNs), giving them advantages over human performance on visual pattern recognition, which makes CNNs a good choice for WAM. Nevertheless, this research is largely unexplored. We thus propose "GenInSAR", a CNN-based generative model for joint phase filtering and coherence estimation, that directly learns the InSAR data distribution. GenInSAR's unsupervised training on satellite and simulated noisy InSAR images outperforms other five related methods in total residue reduction (over 16.5% better on average) with less over-smoothing/artefacts around branch cuts. GenInSAR's Phase, and Coherence Root-Mean-Squared-Error and Phase Cosine Error have average improvements of 0.54, 0.07, and 0.05 respectively compared to the related methods.Comment: to be published in a future issue of IEEE Geoscience and Remote Sensing Letter

Monitoring land subsidence of airport using InSAR time-series techniques with atmospheric and orbital error corrections

Land subsidence is one of the common geological hazards worldwide and mostly caused by human activities including the construction of massive infrastructures. Large infrastructure such as airport is susceptible to land subsidence due to several factors. Therefore, monitoring of the land subsidence at airport is crucial in order to prevent undesirable loss of property and life. Remote sensing technique, especially Interferometric Synthetic Aperture Radar (InSAR) has been successfully applied to measure the surface deformation over the past few decades although atmospheric artefact and orbital errors are still a concerning issue in this measurement technique. Multi-temporal InSAR, an extension of InSAR technique, uses large sets of SAR scenes to investigate the temporal evolution of surface deformation and mitigate errors found in a single interferogram. This study investigates the long-term land subsidence of the Kuala Lumpur International Airport (KLIA), Malaysia and Singapore Changi Airport (SCA), Singapore by using two multi-temporal InSAR techniques like Small Baseline Subset (SBAS) and Multiscale InSAR Time Series (MInTS). General InSAR processing was conducted to generate interferogram using ALOS PALSAR data from 2007 until 2011. Atmospheric and orbital corrections were carried out for all interferograms using weather model, namely European Centre for Medium Range Weather Forecasting (ECMWF) and Network De-Ramping technique respectively before estimating the time series land subsidence. The results show variation of subsidence with respect to corrections (atmospheric and orbital) as well as difference between multi-temporal InSAR techniques (SBAS and MInTS) used. After applying both corrections, a subsidence ranging from 2 to 17 mm/yr was found at all the selected areas at the KLIA. Meanwhile, for SCA, a subsidence of about less than 10 mm/yr was found. Furthermore, a comparison between two techniques (SBAS and MInTS) show a difference rate of subsidence of about less than 1 mm/yr for both study area. SBAS technique shows more linear result as compared to the MInTS technique which shows slightly scattering pattern but both techniques show a similar trend of surface deformation in both study sites. No drastic deformation was observed in these two study sites and slight deformation was detected which about less than 20mm/yr for both study areas probably occurred due to several reasons including conversion of the land use from agricultural land, land reclamation process and also poor construction. This study proved that InSAR time series surface deformation measurement techniques are useful as well as capable to monitor deformation of large infrastructure such as airport and as an alternative to costly conventional ground measurement for infrastructure monitoring

Non-local methods for InSAR parameters estimation

In the thesis work the nonlocal paradigm has been investigated in the framework of Multitemporal SAR Interferometry, e.g. Differential Interferometry, Tomography, etc., and single InSAR pair, e.g. DEM generation. In the former, Adaptive Multi-Looking methods have been developed for the generation of interferometric data-stacks. Following the nonlocal approach, the proposed methods rely only on similar pixels according to a suitable similarity measure that exploits the stack's temporal information. An hybrid approach that jointly uses the nonlocal paradigm and transform domain filtering has been investigated for InSAR pair phase estimation. On the track of the BM3D and SARBM3D algorithms, different approaches to the filtering in the transform domain are investigated. Furthermore, a novel approach to the similarity computation and filtering, based on a relative-topography content of the interferometric phase rather than its absolute value, is proposed