An official website of the United States government
Multistatic GPR has the advantages of reducing survey time and leverages more comprehensive data collection. Traditionally in multistatic GPR data processing, the 2D Bscan image obtained from each receive antenna are simply stacked for 3D image reconstructions. However, such approach is typically inadequate as the multistatic GPR receivers are mounted with spatial offsets, causing back-scattering signals from the same target to have differing time of arrivals. For proper fusion of multistatic GPR data, migration methods that consider the transmitters and receivers spatial offset and data variations among different receiving antennas may be employed. In this study, the back-projection algorithm (BPA) is investigated. The algorithm consists of determining the wave travel path and associated travel time, and projecting the corresponding signal value back into space domain. Furthermore, antenna radiation pattern is incorporated. The BPA enables scatter shape reconstruction and is prone to parallel computing. For validation, multistatic GPR 3D tomographic image reconstruction is successfully applied to laboratory data.
more »
« less
3D tomography for multistatic GPR subsurface sensing
Ground penetrating radar (GPR) subsurface sensing is a promising nondestructive evaluation (NDE) technique for inspecting and surveying underground utilities in complex urban environments, as well as for monitoring other key infrastructure such as bridges and railroads. A challenge of such technique lies on image formation from the recorded GPR data. In this work, a fast back projection algorithm (BPA) for three-dimensional GPR image construction is explored. The BPA is a time-domain migration method that has been effectively used in GPR image formation. However, most of the studies in the literature apply a computationally intensive BPA to a two-dimensional dataset under the assumption that an in-plane scattering occurs underneath the GPR antennas. This assumption is not precise for 3D GPR image formation as the GPR radiation scatters in multiple directions as it reaches the ground. In this study, a generalized form for an approximation to determine the scattering point in an air-coupled GPR system is developed which considerably reduces the required computations and can accurately localize the scattering point position. The algorithm is evaluated by applications on GPR data synthesized using GprMax, a finite-difference time domain (FDTD) simulator.
more »
« less
- PAR ID:
- 10074795
- Date Published:
- Journal Name:
- SPIE Defense + Security Radar Sensor Technology XXII
- Volume:
- 10633
- Page Range / eLocation ID:
- 02
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Ground penetrating radar (GPR) is a remote geophysical sensing method that has been applied in the localization of underground utilities, bridge deck survey, localization of landmines, mapping of terrain for aid in driverless cars, etc. Multistatic GPR can deliver a faster survey, wider spatial coverage, and multiple viewpoints of the subsurface. However, because of the transmit and receive antennas spatial offset, formation of 3D GPR image by simple stacking of the acquired A-scans is inaccurate. Also, averaging of different receivers data may lead to destructive interference of back-scattered waves due to different time delays implied by the spatial offset, so averaging does not lead to higher SNR in general. Furthermore, the energy back-scattered by scatter points are spread in hyperbolas in the GPR raw data. Migration or imaging algorithms are employed to increase SNR by focusing the hyperbolas. This focusing process also leads to better accuracy in target localization. In this paper, a computationally efficient synthetic aperture radar (SAR) imaging algorithm that properly integrates multistatic GPR data in both ground and air-coupled cases is presented. The algorithm is successfully applied on two synthetic datasets.more » « less
-
Abstract Image subtraction is essential for transient detection in time-domain astronomy. The point-spread function (PSF), photometric scaling, and sky background generally vary with time and across the field of view for imaging data taken with ground-based optical telescopes. Image subtraction algorithms need to match these variations for the detection of flux variability. An algorithm that can be fully parallelized is highly desirable for future time-domain surveys. Here we introduce the saccadic fast Fourier transform (SFFT) algorithm we developed for image differencing. SFFT uses aδ-function basis for kernel decomposition, and the image subtraction is performed in Fourier space. This brings about a remarkable improvement in computational performance of about an order of magnitude compared to other published image subtraction codes. SFFT can accommodate the spatial variations in wide-field imaging data, including PSF, photometric scaling, and sky background. However, the flexibility of theδ-function basis may also make it more prone to overfitting. The algorithm has been tested extensively on real astronomical data taken by a variety of telescopes. Moreover, the SFFT code allows for the spatial variations of the PSF and sky background to be fitted by spline functions.more » « less
-
The phase-field method is an attractive computational tool for simulating microstructural evolution during phase separation, including solidification and spinodal decomposition. However, the high computational cost associated with solving phase-field equations currently limits our ability to comprehend phase transformations. This article reports a novel phase-field emulator based on the tensor decomposition of the evolving microstructures and their corresponding two-point correlation functions to predict microstructural evolution at arbitrarily small time scales that are otherwise nontrivial to achieve using traditional phase-field approaches. The reported technique is based on obtaining a low-dimensional representation of the microstructures via tensor decomposition, and subsequently, predicting the microstructure evolution in the low-dimensional space using Gaussian process regression (GPR). Once we obtain the microstructure prediction in the low-dimensional space, we employ a hybrid input–output phase-retrieval algorithm to reconstruct the microstructures. As proof of concept, we present the results on microstructure prediction for spinodal decomposition, although the method itself is agnostic of the material parameters. Results show that we are able to predict microstructure evolution sequences that closely resemble the true microstructures (average normalized mean square of 6.78×10^−7) at time scales half of that employed in obtaining training data. Our data-driven microstructure emulator opens new avenues to predict the microstructural evolution by leveraging phase-field simulations and physical experimentation where the time resolution is often quite large due to limited resources and physical constraints, such as the phase coarsening experiments previously performed in microgravity.more » « less
-
Isaacs, Jason C.; Bishop, Steven S. (Ed.)Ultra-wideband (UWB) ground penetrating radar (GPR) is an effective, widely used tool for detection and mapping of buried targets. However, traditional ground penetrating radar systems struggle to resolve and identify congested target configurations and irregularly shaped targets. This is a significant limitation for many municipalities who seek to use GPR to locate and image underground utility pipes. This research investigates the implementation of orbital angular momentum (OAM) control in an UWB GPR, with the goal of addressing these limitations. Control of OAM is a novel technique which leverages an additional degree of freedom offered by spatially structured helical waveforms. This paper examines several free-space and buried target configurations to determine the ability of helical OAM waveforms to improve detectability and distinguishability of buried objects including those with symmetric, asymmetric, and chiral geometries. Microwave OAM can be generated using a uniform circular array (UCA) of antennas with phase delays applied according to azimuth angle. Here, a four-channel network analyzer transceiver is connected to a UCA to enable UWB capability. The characteristic phase delays of OAM waveforms are implemented synthetically via signal processing. The viability demonstrated with the method opens design and analysis degrees of freedom for penetrating radar that may help with discerning challenging targets, such as buried landmines and wires.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Conference Paper:
- https://doi.org/10.1117/12.2304423
