Efficient SAR Raw Data Compression in Frequency Domain

SAR raw data compression is necessary to reduce huge amounts of SAR data for a memory on board a satellite, space shuttle or aircraft and for later downlink to a ground station. In view of interferometric and polarimetric applications for SAR data, it becomes more and more important to pay attention to phase errors caused by data compression. Herein, a detailed comparison of block adaptive quantization in time domain (BAQ) and in frequency domain (FFT-BAQ) is given. Inclusion of raw data compression in the processing chain allows an efficient use of the FFT-BAQ and makes implementation for on-board data compression feasible. The FFT-BAQ outperforms the BAQ in terms of signal-to-quantization noise ratio and phase error and allows a direct decimation of the oversampled data equivalent to FIR-filtering in time domain. Impacts on interferometric phase and coherency are also given

Non-Local Compressive Sensing Based SAR Tomography

Tomographic SAR (TomoSAR) inversion of urban areas is an inherently sparse reconstruction problem and, hence, can be solved using compressive sensing (CS) algorithms. This paper proposes solutions for two notorious problems in this field: 1) TomoSAR requires a high number of data sets, which makes the technique expensive. However, it can be shown that the number of acquisitions and the signal-to-noise ratio (SNR) can be traded off against each other, because it is asymptotically only the product of the number of acquisitions and SNR that determines the reconstruction quality. We propose to increase SNR by integrating non-local estimation into the inversion and show that a reasonable reconstruction of buildings from only seven interferograms is feasible. 2) CS-based inversion is computationally expensive and therefore barely suitable for large-scale applications. We introduce a new fast and accurate algorithm for solving the non-local L1-L2-minimization problem, central to CS-based reconstruction algorithms. The applicability of the algorithm is demonstrated using simulated data and TerraSAR-X high-resolution spotlight images over an area in Munich, Germany.Comment: 10 page

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

An accurate method to correct atmospheric phase delay for InSAR with the ERA5 global atmospheric model

Differential SAR Interferometry (DInSAR) has proven its unprecedented ability and merits of monitoring ground deformation on a large scale with centimeter to millimeter accuracy. However, atmospheric artifacts due to spatial and temporal variations of the atmospheric state often affect the reliability and accuracy of its results. The commonly-known Atmospheric Phase Screen (APS) appears in the interferograms as ghost fringes not related to either topography or deformation. Atmospheric artifact mitigation remains one of the biggest challenges to be addressed within the DInSAR community. State-of-the-art research works have revealed that atmospheric artifacts can be partially compensated with empirical models, point-wise GPS zenith path delay, and numerical weather prediction models. In this study, we implement an accurate and realistic computing strategy using atmospheric reanalysis ERA5 data to estimate atmospheric artifacts. With this approach, the Line-of-Sight (LOS) path along the satellite trajectory and the monitored points is considered, rather than estimating it from the zenith path delay. Compared with the zenith delay-based method, the key advantage is that it can avoid errors caused by any anisotropic atmospheric phenomena. The accurate method is validated with Sentinel-1 data in three different test sites: Tenerife island (Spain), Almería (Spain), and Crete island (Greece). The effectiveness and performance of the method to remove APS from interferograms is evaluated in the three test sites showing a great improvement with respect to the zenith-based approach.Peer ReviewedPostprint (published version

Tropospheric phase delay in interferometric synthetic aperture radar estimated from meteorological model and multispectral imagery

ENVISAT Medium Resolution Imaging Spectrometer Instrument (MERIS) multispectral data and the mesoscale meteorological model MM5 are used to estimate the tropospheric phase delay in synthetic aperture radar (SAR) interferograms. MERIS images acquired simultaneously with ENVISAT Advanced Synthetic Aperture Radar data provide an estimate of the total water vapor content W limited to cloud-free areas based on spectral bands ratio (accuracy 0.17 g cm^(−2) and ground resolution 300 m). Maps of atmospheric delay, 2 km in ground resolution, are simulated from MM5. A priori pertinent cumulus parameterization and planetary boundary layer options of MM5 yield near-equal phase correction efficiency. Atmospheric delay derived from MM5 is merged with available MERIS W product. Estimates of W measured from MERIS and modeled from MM5 are shown to be consistent and unbiased and differ by ~0.2 g cm^(−2) (RMS). We test the approach on data over the Lebanese ranges where active tectonics might contribute to a measurable SAR signal that is obscured by atmospheric effects. Local low-amplitude (1 rad) atmospheric oscillations with a 2.25 km wavelength on the interferograms are recovered from MERIS with an accuracy of 0.44 rad or 0.03 g cm^(−2). MERIS water product overestimates W in the clouds shadow due to mismodeling of multiple scattering and underestimates W on pixels with undetected semitransparent clouds. The proposed atmospheric filter models dynamic atmospheric signal which cannot be recovered by previous filtering techniques which are based on a static atmospheric correction. Analysis of filter efficiency with spatial wavelength shows that ~43% of the atmospheric signal is removed at all wavelengths

Image fusion techniqes for remote sensing applications

Image fusion refers to the acquisition, processing and synergistic combination of information provided by various sensors or by the same sensor in many measuring contexts. The aim of this survey paper is to describe three typical applications of data fusion in remote sensing. The first study case considers the problem of the Synthetic Aperture Radar (SAR) Interferometry, where a pair of antennas are used to obtain an elevation map of the observed scene; the second one refers to the fusion of multisensor and multitemporal (Landsat Thematic Mapper and SAR) images of the same site acquired at different times, by using neural networks; the third one presents a processor to fuse multifrequency, multipolarization and mutiresolution SAR images, based on wavelet transform and multiscale Kalman filter. Each study case presents also results achieved by the proposed techniques applied to real data

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