Scattering statistics of rock outcrops: Model-data comparisons and Bayesian inference using mixture distributions

The probability density function of the acoustic field amplitude scattered by the seafloor was measured in a rocky environment off the coast of Norway using a synthetic aperture sonar system, and is reported here in terms of the probability of false alarm. Interpretation of the measurements focused on finding appropriate class of statistical models (single versus two-component mixture models), and on appropriate models within these two classes. It was found that two-component mixture models performed better than single models. The two mixture models that performed the best (and had a basis in the physics of scattering) were a mixture between two K distributions, and a mixture between a Rayleigh and generalized Pareto distribution. Bayes' theorem was used to estimate the probability density function of the mixture model parameters. It was found that the K-K mixture exhibits significant correlation between its parameters. The mixture between the Rayleigh and generalized Pareto distributions also had significant parameter correlation, but also contained multiple modes. We conclude that the mixture between two K distributions is the most applicable to this dataset.Comment: 15 pages, 7 figures, Accepted to the Journal of the Acoustical Society of Americ

Bayesian off-line detection of multiple change-points corrupted by multiplicative noise : application to SAR image edge detection

This paper addresses the problem of Bayesian off-line change-point detection in synthetic aperture radar images. The minimum mean square error and maximum a posteriori estimators of the changepoint positions are studied. Both estimators cannot be implemented because of optimization or integration problems. A practical implementation using Markov chain Monte Carlo methods is proposed. This implementation requires a priori knowledge of the so-called hyperparameters. A hyperparameter estimation procedure is proposed that alleviates the requirement of knowing the values of the hyperparameters. Simulation results on synthetic signals and synthetic aperture radar images are presented

Method for detecting surface motions and mapping small terrestrial or planetary surface deformations with synthetic aperture radar

A technique based on synthetic aperture radar (SAR) interferometry is used to measure very small (1 cm or less) surface deformations with good resolution (10 m) over large areas (50 km). It can be used for accurate measurements of many geophysical phenomena, including swelling and buckling in fault zones, residual, vertical and lateral displacements from seismic events, and prevolcanic swelling. Two SAR images are made of a scene by two spaced antennas and a difference interferogram of the scene is made. After unwrapping phases of pixels of the difference interferogram, surface motion or deformation changes of the surface are observed. A second interferogram of the same scene is made from a different pair of images, at least one of which is made after some elapsed time. The second interferogram is then compared with the first interferogram to detect changes in line of sight position of pixels. By resolving line of sight observations into their vector components in other sets of interferograms along at least one other direction, lateral motions may be recovered in their entirety. Since in general, the SAR images are made from flight tracks that are separated, it is not possible to distinguish surface changes from the parallax caused by topography. However, a third image may be used to remove the topography and leave only the surface changes

Synthetic aperture radar/LANDSAT MSS image registration

Algorithms and procedures necessary to merge aircraft synthetic aperture radar (SAR) and LANDSAT multispectral scanner (MSS) imagery were determined. The design of a SAR/LANDSAT data merging system was developed. Aircraft SAR images were registered to the corresponding LANDSAT MSS scenes and were the subject of experimental investigations. Results indicate that the registration of SAR imagery with LANDSAT MSS imagery is feasible from a technical viewpoint, and useful from an information-content viewpoint

Improving InSAR geodesy using global atmospheric models

Spatial and temporal variations of pressure, temperature and water vapor content in the atmosphere introduce significant confounding delays in Interferometric Synthetic Aperture Radar (InSAR) observations of ground deformation and bias estimatesof regional strain rates. Producing robust estimates of tropospheric delays remains one of the key challenges in increasing the accuracy of ground deformation measurements using InSAR. Recent studies revealed the efficiency of global atmospheric reanalysis to mitigate the impact of tropospheric delays, motivating further exploration of their potential. Here, we explore the effectiveness of these models in several geographic and tectonic settings on both single interferograms and time series analysis products. Both hydrostatic and wet contributions to the phase delay are important to account for. We validate these path delay corrections by comparing with estimates of vertically integrated atmospheric water vapor content derived from the passive multi-spectral imager MERIS, onboard the ENVISAT satellite. Generally, the performance of the prediction depends on the vigor of atmospheric turbulence. We discuss (1) how separating atmospheric and orbital contributions allows one to better measure long wavelength deformation, (2) how atmospheric delays affect measurements of surface deformation following earthquakes and (3) we show that such a method allows us to reduce biases in multi-year strain rate estimates by reducing the influence of unevenly sampled seasonal oscillations of the tropospheric delay

Advanced InSAR atmospheric correction: MERIS/MODIS combination and stacked water vapour models

A major source of error for repeat-pass Interferometric Synthetic Aperture Radar (InSAR) is the phase delay in radio signal propagation through the atmosphere (especially the part due to tropospheric water vapour). Based on experience with the Global Positioning System (GPS)/Moderate Resolution Imaging Spectroradiometer (MODIS) integrated model and the Medium Resolution Imaging Spectrometer (MERIS) correction model, two new advanced InSAR water vapour correction models are demonstrated using both MERIS and MODIS data: (1) the MERIS/MODIS combination correction model (MMCC); and (2) the MERIS/MODIS stacked correction model (MMSC). The applications of both the MMCC and MMSC models to ENVISAT Advanced Synthetic Aperture Radar (ASAR) data over the Southern California Integrated GPS Network (SCIGN) region showed a significant reduction in water vapour effects on ASAR interferograms, with the root mean square (RMS) differences between GPS- and InSAR-derived range changes in the line-of-sight (LOS) direction decreasing from ,10mm before correction to ,5mm after correction, which is similar to the GPS/MODIS integrated and MERIS correction models. It is expected that these two advanced water vapour correction models can expand the application of MERIS and MODIS data for InSAR atmospheric correction. A simple but effective approach has been developed to destripe Terra MODIS images contaminated by radiometric calibration errors. Another two limiting factors on the MMCC and MMSC models have also been investigated in this paper: (1) the impact of the time difference between MODIS and SAR data; and (2) the frequency of cloud-free conditions at the global scale

A robust nonlinear scale space change detection approach for SAR images

In this paper, we propose a change detection approach based on nonlinear scale space analysis of change images for robust detection of various changes incurred by natural phenomena and/or human activities in Synthetic Aperture Radar (SAR) images using Maximally Stable Extremal Regions (MSERs). To achieve this, a variant of the log-ratio image of multitemporal images is calculated which is followed by Feature Preserving Despeckling (FPD) to generate nonlinear scale space images exhibiting different trade-offs in terms of speckle reduction and shape detail preservation. MSERs of each scale space image are found and then combined through a decision level fusion strategy, namely "selective scale fusion" (SSF), where contrast and boundary curvature of each MSER are considered. The performance of the proposed method is evaluated using real multitemporal high resolution TerraSAR-X images and synthetically generated multitemporal images composed of shapes with several orientations, sizes, and backscatter amplitude levels representing a variety of possible signatures of change. One of the main outcomes of this approach is that different objects having different sizes and levels of contrast with their surroundings appear as stable regions at different scale space images thus the fusion of results from scale space images yields a good overall performance

Post-failure evolution analysis of a rainfall-triggered landslide by multi-temporal interferometry SAR approaches integrated with geotechnical analysis

Persistent Scatterers Interferometry (PSI) represents one of the most powerful techniques for Earth's surface deformation processes' monitoring, especially for long-term evolution phenomena. In this work, a dataset of 34 TerraSAR-X StripMap images (October 2013–October 2014) has been processed by two PSI techniques - Coherent Pixel Technique-Temporal Sublook Coherence (CPT-TSC) and Small Baseline Subset (SBAS) - in order to study the evolution of a slow-moving landslide which occurred on February 23, 2012 in the Papanice hamlet (Crotone municipality, southern Italy) and induced by a significant rainfall event (185 mm in three days). The mass movement caused structural damage (buildings' collapse), and destruction of utility lines (gas, water and electricity) and roads. The results showed analogous displacement rates (30–40 mm/yr along the Line of Sight – LOS-of the satellite) with respect to the pre-failure phase (2008–2010) analyzed in previous works. Both approaches allowed detect the landslide-affected area, however the higher density of targets identified by means of CPT-TSC enabled to analyze in detail the slope behavior in order to design possible mitigation interventions. For this aim, a slope stability analysis has been carried out, considering the comparison between groundwater oscillations and time-series of displacement. Hence, the crucial role of the interaction between rainfall and groundwater level has been inferred for the landslide triggering. In conclusion, we showed that the integration of geotechnical and remote sensing approaches can be seen as the best practice to support stakeholders to design remedial works.Peer ReviewedPostprint (author's final draft

Charge-coupled device data processor for an airborne imaging radar system

Processing of raw analog echo data from synthetic aperture radar receiver into images on board an airborne radar platform is discussed. Processing is made feasible by utilizing charge-coupled devices (CCD). CCD circuits are utilized to perform input sampling, presumming, range correlation and azimuth correlation in the analog domain. These radar data processing functions are implemented for single-look or multiple-look imaging radar systems

A review of RFI mitigation techniques in microwave radiometry

Radio frequency interference (RFI) is a well-known problem in microwave radiometry (MWR). Any undesired signal overlapping the MWR protected frequency bands introduces a bias in the measurements, which can corrupt the retrieved geophysical parameters. This paper presents a literature review of RFI detection and mitigation techniques for microwave radiometry from space. The reviewed techniques are divided between real aperture and aperture synthesis. A discussion and assessment of the application of RFI mitigation techniques is presented for each type of radiometer.Peer ReviewedPostprint (published version