Zero source insertion technique to account for undersampling in GPR imaging

• US 8,659,467 B1
• Filed: 08/26/2011
• Issued: 02/25/2014
• Est. Priority Date: 08/26/2010
• Status: Active Grant

• Current Assignee: Lawrence Livermore National Security LLC (Government of the United States of America)
• Original Assignee: Lawrence Livermore National Security LLC (Government of the United States of America)
• Inventors: Chambers, David H., Mast, Jeffrey E., Paglieroni, David W.
• Primary Examiner(s): Sotomayor, John B
• Assistant Examiner(s): WINDRICH, MARCUS E
• Application Number:
  US13/219,504
• Time in Patent Office: 914 Days
• Field of Search: 342/22, 342 25 R- 25 F, 342/27, 342/162, 342/191
• US Class Current: 342/22
• CPC Class Codes: G01S 13/003 Bistatic radar systems; Multistatic radar systems
• CPC Class Codes: G01S 13/885 for ground probing prospecting or detecting using electromagnetic waves G01V3/12
• CPC Class Codes: G01S 13/90 using synthetic aperture techniques , e.g. synthetic aperture radar [SAR] techniques
• CPC Class Codes: G01V 3/12 operating with electromagnetic waves operating with millimetre waves G01V8/005
• CPC Class Codes: G01V 3/15 specially adapted for use during transport, e.g. by a person, vehicle or boat
• CPC Class Codes: G06N 99/00 Subject matter not provided for in other groups of this subclass

Abstract

A method and system for detecting the presence of subsurface objects within a medium is provided. In some embodiments, the imaging and detection system operates in a multistatic mode to collect radar return signals generated by an array of transceiver antenna pairs that is positioned across the surface and that travels down the surface. The imaging and detection system pre-processes the return signal to suppress certain undesirable effects. The imaging and detection system then generates synthetic aperture radar images from real aperture radar images generated from the pre-processed return signal. The imaging and detection system then post-processes the synthetic aperture radar images to improve detection of subsurface objects. The imaging and detection system identifies peaks in the energy levels of the post-processed image frame, which indicates the presence of a subsurface object.

Claims

1. A method performed by a computing device for generating a contribution image frame for a transmitter based on return signals acquired by multiple receivers from a signal emitted by the transmitter at a down-track location, the receivers being at cross-track locations, the method comprising:

providing a matrix of the return signals with each column representing a receiver and each row representing a fast time delay;

applying a fast Fourier transform separately to each column of the matrix to generate a column transformed matrix;

inserting columns of zeros into the column transformed matrix to account for under-sampling in a cross-track direction to generate a zero inserted matrix;

applying a fast Fourier transform separately to each row of the zero inserted matrix to generate a row transformed matrix;

multiplying the row transformed matrix by a factor to propagate the matrix to a specific depth to generate a propagated row transformed matrix;

applying an inverse fast Fourier transform separately to each row of the propagated row transformed matrix to generate an inverse row transformed matrix;

multiplying elements of the inverse row transformed matrix by a factor representing a field emitted by the transmitter to generate a modified inverse row transformed matrix;

applying an inverse fast Fourier transform separately to each column of the modified inverse row transformed matrix to generate an inverse column transformed matrix; and

creating the contribution image frame from elements of the inverse column transformed matrix.

2. The method of claim 1 wherein the number of columns of zeros inserted between each column of the column transformed matrix is sufficient to ensure that the total number of columns in the cross-track direction is greater than the width of the contribution image frame divided by a pixel width in an image frame.

3. The method of claim 2 wherein the number of columns of zeros is based on the sampling theorem.

4. The method of claim 1 wherein the emitted signal is a ground-penetrating radar signal that is emitted into the ground for detecting a buried explosive device.

5. The method of claim 1 including creating a contribution image frame for each of a plurality of transmitters at the down-track location and generating a real aperture image for the down-track location from the created contribution image frames.

6. The method of claim 5 wherein the transmitters and receivers operate in a multistatic mode.

7. The method of claim 5 including generating real aperture images for a plurality of down-track locations and generating a synthetic aperture image for a target down-track location from the real aperture images of adjacent down-track locations.

8. A computer-readable storage device storing computer-executable instructions for controlling a computing device to generate a contribution image frame for a transmitter based on return signals acquired by multiple receivers from a signal emitted by the transmitter, the receivers being at locations, the return signal for each receiver including a number of samples indicating amplitude of the return signal, each sample corresponding to a fast time delay at which the signal emitted by the transmitter was acquired by that receiver, by a method comprising:

transforming the return signal of each receiver from a time domain into a frequency domain to generate first transformed return signals;

augmenting the first transformed return signals with zero return signals to represent virtual receivers to account for under-sampling by the receivers to generate zero inserted return signals;

transforming the zero inserted return signals across the receivers for each fast time delay in a frequency domain to generate second transformed return signals;

inversely transforming the zero inserted return signals across the receivers for each fast time delay to generate first inverse transformed return signals;

inversely transforming the first inverse transformed return signals for each receiver from the frequency domain to a time domain to generate second inverse transformed return signals; and

creating the contribution image frame from elements of the inverse transformed return signals.

9. The computer-readable storage device of claim 8 wherein the transforming is performed using a fast Fourier transform and the inverse transforming is performed using an inverse fast Fourier transform.

10. The computer-readable storage device of claim 8 wherein the number of virtual receivers is sufficient to ensure that the total number of receivers is greater than a width of the contribution image frame divided by a pixel width of an image frame.

11. The computer-readable storage device of claim 10 wherein the number of virtual receivers is based on the sampling theorem.

12. The computer-readable storage device of claim 8 wherein the emitted signal is a ground-penetrating radar signal that is emitted into the ground for detecting a buried explosive device.

13. The computer-readable storage device of claim 8 including creating a contribution image frame for each of a plurality of transmitters and generating a real aperture image from the created contribution image frames.

14. The computer-readable storage device of claim 8 wherein the transmitters and receivers operate in a multistatic mode.

15. The computer-readable storage device of claim 8 including generating real aperture images from the contribution image frames transmitters at cross-track locations for a plurality of down-track locations and generating a synthetic aperture image for a target down-track location from the real aperture images of adjacent down-track locations.

16. A computing device for generating a contribution image frame for a transmitter based on return signals acquired by multiple receivers from a signal emitted by the transmitter at a down-track location, the receivers being at cross-track locations, comprising:

a data store that stores a matrix of the return signals with each column representing a receiver and each row representing a fast time delay;

a component that applies a fast Fourier transform to each column of the matrix to generate a column transformed matrix;

a component that inserts columns of zeros into the column transformed matrix to account for under-sampling in a cross-track direction to generate a zero inserted matrix;

a component that applies a fast Fourier transform to each row of the matrix to generate a row transformed matrix;

a component that applies an inverse fast Fourier transform to each row of the zero inserted matrix to generate a inverse row transformed matrix;

a component that modifies the matrix based on a field from the transmitter to generate a modified inverse row transformed matrix;

a component that applies an inverse fast Fourier transform to each column of the modified inverse row transformed matrix to generate an inverse column transformed matrix; and

a component that creates the contribution image frame from elements of the inverse column transformed matrix.

17. The computing device of claim 16 wherein the number of columns of zeros inserted between each column of the column transformed matrix is sufficient to ensure that the total number of columns in the cross-track direction is greater than the width of the contribution image frame divided by a pixel width in an image frame.

18. The computing device of claim 16 wherein the emitted signal is a ground-penetrating radar signal that is emitted into the ground for detecting a buried explosive device.

19. The computing device of claim 16 including a component that creates a contribution image frame for each of a plurality of transmitters at the down-track location and generating a real aperture image for the down-track location from the created contribution image frames.

20. The computing device of claim 16 wherein the transmitters and receivers operate in a multistatic mode and including a component that generates real aperture images for a plurality of down-track locations and a component that generates a synthetic aperture image for a target down-track location from the real aperture images of adjacent down-track locations.