Three-Dimensional Nepal Earthquake Displacement Using Hybrid Genetic Algorithm Phase Unwrapping from Sentinel-1A Satellite

Introduction: Geophysicists had forewarned for decades that Nepal was exposed to a deadly earthquake, exceptionally despite its geology, urbanization and architecture. Gorkha earthquake is the most horrible natural disaster to crash into Nepal since the 1934 Nepal-Bihar earthquake. Gorkha earthquake occurred on April 25, 2015, at 11:56 NST and killed more than 10,000 people and injured more than 23,000 population. Objective: The main objective of this work is to utilize hybrid genetic algorithm for three-dimensional phase unwrapping of Nepal earthquake displacement using Sentinel-1A satellite. The three-dimensional best-path avoiding singularity loops (3DBPASL) algorithm was implemented to perform 3D Sentinel-1A satellite phase unwrapping. The hybrid genetic algorithm (HGA) was used to achieve 3DBPASL phase matching. Advancely, the errors in phase decorrelation were reduced by optimization of 3DBPASL using HGA. Results: The findings indicate a few cm of ground deformation and vertical northern of Kathmandu. Approximately, an area of 12,000 km2 has been drifted also the northern of Kathmandu. Further, each fringe of colour represents about 2.5 cm of deformation. The large amount of fringes indicates a large deformation pattern with ground motions of 3 m. Conclusion: In conclusion, HGA can be used to produce accurate 3D quake deformation using Sentinel-1A satellite

High-speed, solid state, interferometric interrogator and multiplexer for fiber Bragg grating sensors

We report on the design and prototyping of a robust high-speed interferometric multiplexer and interrogator for fiber Bragg grating sensors. The scheme is based on the combination of active WDM channel switching and passive, instantaneous interferometry, allowing the resolution of virgin interferometric interrogators to be retained at MHz multiplexing rates. In this article the system design and operation are described, and a prototype scheme is characterized for three sensors and a multiplexing rate of 4 kHz, demonstrating a noise floor of 10 nε/√Hz and no cross-sensitivity. It is proposed that the system will be applicable to demanding monitoring applications requiring high speed and high resolution measurements across the sensor array

Earthquake source parameters of the 2009 Mw 7.8 Fiordland (New Zealand) earthquake from L-band InSAR observations

The 2009 MW7.8 Fiordland (New Zealand) earthquake is the largest to have occurred in New Zealand since the 1931 Mw 7.8 Hawke’s Bay earthquake, 1 000 km to the northwest. In this paper two tracks of ALOS PALSAR interferograms (one ascending and one descending) are used to determine fault geometry and slip distribution of this large earthquake. Modeling the event as dislocation in an elastic half-space suggests that the earthquake resulted from slip on a SSW-NNE orientated thrust fault that is associated with the subduction between the Pacific and Australian Plates, with oblique displacement of up to 6.3 m. This finding is consistent with the preliminary studies undertaken by the USGS using seismic data

Development of temporal phase unwrapping algorithms for depth-resolved measurements using an electronically tuned Ti:Sa laser

This thesis is concerned with (a) the development of full-field, multi-axis and phase contrast wavelength scanning interferometer, using an electronically tuned CW Ti:Sa laser for the study of depth resolved measurements in composite materials such as GFRPs and (b) the development of temporal phase unwrapping algorithms for depth re-solved measurements. Item (a) was part of the ultimate goal of successfully extracting the 3-D, depth-resolved, constituent parameters (Young s modulus E, Poisson s ratio v etc.) that define the mechanical behaviour of composite materials like GFRPs. Considering the success of OCT as an imaging modality, a wavelength scanning interferometer (WSI) capable of imaging the intensity AND the phase of the interference signal was proposed as the preferred technique to provide the volumetric displacement/strain fields (Note that displacement/strain fields are analogous to phase fields and thus a phase-contrast interferometer is of particular interest in this case). These would then be passed to the VFM and yield the sought parameters provided the loading scheme is known. As a result, a number of key opto-mechanical hardware was developed. First, a multiple channel (x6) tomographic interferometer realised in a Mach-Zehnder arrangement was built. Each of the three channels would provide the necessary information to extract the three orthogonal displacement/strain components while the other three are complementary and were included in the design in order to maximize the penetration depth (sample illuminated from both sides). Second, a miniature uniaxial (tensile and/or compression) loading machine was designed and built for the introduction of controlled and low magnitude displacements. Last, a rotation stage for the experimental determination of the sensitivity vectors and the re-registration of the volumetric data from the six channels was also designed and built. Unfortunately, due to the critical failure of the Ti:Sa laser data collection using the last two items was not possible. However, preliminary results at a single wavelength suggested that the above items work as expected. Item (b) involved the development of an optical sensor for the dynamic monitoring of wavenumber changes during a full 100 nm scan. The sensor is comprised of a set of four wedges in a Fizeau interferometer setup that became part of the multi-axis interferometer (7th channel). Its development became relevant due to the large amount of mode-hops present during a full scan of the Ti:Sa source. These are associated to the physics of the laser and have the undesirable effect of randomising the signal and thus preventing successful depth reconstructions. The multi-wedge sensor was designed so that it provides simultaneously high wavenumber change resolution and immunity to the large wavenumber jumps from the Ti:Sa. The analysis algorithms for the extraction of the sought wavenumber changes were based on 2-D Fourier transform method followed by temporal phase unwrapping. At first, the performance of the sensor was tested against that of a high-end commercial wavemeter for a limited scan of 1nm. A root mean square (rms) difference in measured wavenumber shift between the two of ∼4 m-1 has been achieved, equivalent to an rms wavelength shift error of ∼0.4 pm. Second, by resampling the interference signal and the wavenumber-change axis onto a uniformly sampled k-space, depth resolutions that are close to the theoretical limits were achieved for scans of up to 37 nm. Access of the full 100 nm range that is characterised by wavelength steps down to picometers level was achieved by introducing a number of improvements to the original temporal phase unwrapping algorithm reported in ref [1] tailored to depth resolved measurements. These involved the estimation and suppression of intensity background artefacts, improvements on the 2-D Fourier transform phase detection based on a previously developed algorithm in ref [2] and finally the introduction of two modifications to the original TPU. Both approaches are adaptive and involve signal re-referencing at regular intervals throughout the scan. Their purpose is to compensate for systematic and non-systematic errors owing to a small error in the value of R (a scaling factor applied to the lower sensitivity wedge phase-change signal used to unwrap the higher sensitivity one), or small changes in R with wavelength due to the possibility of a mismatch in the refractive dispersion curves of the wedges and/or a mismatch in the wedge angles. A hybrid approach combining both methods was proposed and used to analyse the data from each of the four wedges. It was found to give the most robust results of all the techniques considered, with a clear Fourier peak at the expected frequency, with significantly reduced spectral artefacts and identical depth resolutions for all four wedges of 2.2 μm measured at FWHM. The ability of the phase unwrapping strategy in resolving the aforementioned issues was demonstrated by successfully measuring the absolute thickness of four fused silica glasses using real experimental data. The results were compared with independent micrometer measurements and showed excellent agreement. Finally, due to the lack of additional experimental data and in an attempt to justify the validity of the proposed temporal phase unwrapping strategy termed as the hybrid approach, a set of simulations that closely matched the parameters characterising the real experimental data set analysed were produced and were subsequently analysed. The results of this final test justify that the various fixes included in the hybrid approach have not evolved to solve the problems of a particular data set but are rather of general nature thereby, highlighting its importance for PC-WSI applications concerning the processing and analysis of large scans

On the integration of deformation and relief measurement using ESPI

The combination of relief and deformation measurement is investigated for improving the accuracy of Electronic Speckle-Pattern Interferometry (ESPI) data. The nature of sensitivity variations within different types of interferometers and with different shapes of objects is analysed, revealing significant variations for some common interferometers. Novel techniques are developed for real-time measurement of dynamic events by means of carrier fringes. This allows quantification of deformation and relief, where the latter is used in the correction of the sensitivity variations of the former

Ocean Surface Observations Using the TanDEM-X Satellite Formation

The TanDEM-X SAR satellite formation permits improved ocean surface observations by means of bistatic along-track interferometry (ATI) when compared to single-satellite systems. The flexible imaging geometry of its two cooperating SAR sensors forms an interferometer that can achieve very high sensitivity to motions of objects on ground. This way, radar imaging of surface currents with unprecedented accuracy, high spatial resolution and wide coverage at the same time becomes possible. We demonstrate the capabilities of the sensor in the contexts of tidal current mapping, measurement of thermohaline and wind-driven ocean currents as well as detection of areas with surface films. We have developed a dedicated postprocessing system for TanDEM-X image products that allows extracting surface current information from the data. By this paper, we address bistatic data acquisition and processing aspects for sea surface imaging with TanDEM-X like interferometric baseline geometry, temporal decorrelation, and phase calibration. We present a variety of examples of data evaluation that clearly demonstrate the application potential of the methodology

An interferometric phase noise reduction method based on modified denoising convolutional neural network

Traditional interferometric synthetic aperture radar (InSAR) denoising methods normally try to estimate the phase fringes directly from the noisy interferogram. Since the statistics of phase noise are more stable than the phase corresponding to complex terrain, it could be easier to estimate the phase noise. In this paper, phase noises rather than phase fringes are estimated first, and then they are subtracted from the noisy interferometric phase for denoising. The denoising convolutional neural network (DnCNN) is introduced to estimate phase noise and then a modified network called IPDnCNN is constructed for the problem. Based on the IPDnCNN, a novel interferometric phase noise reduction algorithm is proposed, which can reduce phase noise while protecting fringe edges and avoid the use of filter windows. Experimental results using simulated and real data are provided to demonstrate the effectiveness of the proposed method

Leveraging power of deep learning for fast and efficient elite pixel selection in time series SAR interferometry

This work proposes an improved convolutional long short-term memory (ConvLSTM) based architecture for selection of elite pixels (i.e., less noisy) in time series interferometric synthetic aperture radar (TS-InSAR). Compared to previous version, the model can process InSAR stacks of variable time steps and select both persistent (PS) and distributed scatterers (DS). We trained the model on ~20,000 training images (interferograms), each of size 100 by 100 pixels, extracted from InSAR time series interferograms containing both artificial features (buildings and infrastructure) and objects of natural environment (vegetation, forests, barren or agricultural land, water bodies). Based on such categorization, we developed two deep learning models, primarily focusing on urban and coastal sites. Training labels were generated from elite pixel selection outputs generated from the wavelet-based InSAR (WabInSAR) software developed by Shirzaei (2013) and improved in Lee and Shirzaei (2023). With 4 urban and 7 coastal sites used for training and validation, the predicted elite pixel selection maps reveal that the proposed models efficiently learn from WabInSAR-generated labels, reaching a validation accuracy of 94%. The models accurately discard pixels affected by geometric and temporal decorrelation while selecting pixels corresponding to urban objects and those with stable phase history unaffected by temporal and geometric decorrelation. The density of pixels in urban areas is comparable to and higher for coastal areas compared to WabInSAR outputs. With significantly reduced time computation (order of minutes) and improved selection of elite pixels, the proposed models can efficiently process long InSAR time series stacks and generate rapid deformation maps

Efficient Phase Unwrapping Architecture for Digital Holographic Microscopy

This paper presents a novel phase unwrapping architecture for accelerating the computational speed of digital holographic microscopy (DHM). A fast Fourier transform (FFT) based phase unwrapping algorithm providing a minimum squared error solution is adopted for hardware implementation because of its simplicity and robustness to noise. The proposed architecture is realized in a pipeline fashion to maximize throughput of the computation. Moreover, the number of hardware multipliers and dividers are minimized to reduce the hardware costs. The proposed architecture is used as a custom user logic in a system on programmable chip (SOPC) for physical performance measurement. Experimental results reveal that the proposed architecture is effective for expediting the computational speed while consuming low hardware resources for designing an embedded DHM system