Real-time system for imaging and object detection with a multistatic GPR array

• US 8,854,248 B2
• Filed: 08/26/2011
• Issued: 10/07/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: Paglieroni, David W., Beer, N. Reginald, Bond, Steven W., Top, Philip L., Chambers, David H., Mast, Jeffrey E., Donetti, John G., Mason, Blake C., Jones, Steven M.
• Primary Examiner(s): BYTHROW, PETER M
• Application Number:
  US13/219,410
• Publication Number:
  US 20130082856A1
• Time in Patent Office: 1,138 Days
• Field of Search: 342/22, 342/25.R, 342/25.A, 342/25.B, 342/25.C, 342/25.D, 342/25.E, 342/25.F, 342/27, 342159-164, 342/175, 342/176, 342/179, 342/191, 342/192, 342/195
• 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 detecting presence of a subsurface object in a medium, the method comprising:

for each of a plurality of down-track locations, acquiring return signals at receivers from signals emitted by transmitters into the medium, the transmitters and receivers operating in a multistatic mode;

pre-processing the return signals to suppress effects of extraneous signals;

for each down-track location, generating by the computing device a real aperture radar image from the pre-processed return signals for that down-track location;

for each down-track location, generating by the computing device a synthetic aperture radar image from the real aperture radar images of that down-track location and adjacent down-track locations;

post-processing the synthetic aperture radar images to separate objects in the foreground from clutter and medium in the background by filtering out effects of clutter followed by restoring spots within the filtered images to their pre-filtered values;

for each down-track location, generating a statistic based on energy of a restored spot within the filtered image for that down-track location; and

identifying down-track locations with peaks in the generated statistics; and

upon identifying a down-track location with a peak, indicating the presence of a subsurface object at that down-track location.

2. The method of claim 1 wherein the medium is ground and the transmitters and receivers form a linear array of transceivers.

3. The method of claim 1 wherein the return signals are acquired and the presence of a subsurface object is indicated in real-time.

4. The method of claim 1 wherein the pre-processing of current return signals includes:

calculating an estimated background signal for a current return signal based on a combination of an estimated background signal for a prior return signal for a prior down-track location that has the same transmitter and receiver pair and a similar air gap to the current return signal and an estimated gain and bias adjusted return signal for the prior return signal;

calculating a gain and bias adjusted return signal for the current return signal; and

calculating an estimated foreground signal for the current return signal based on a difference between the gain and bias adjusted return signal for the current return signal and the estimated background signal for the current return signal.

5. The method of claim 1 wherein the generating of a real aperture radar image for a down-track location includes:

performing a plane-to-plane propagation of the return signals to generate a contribution image frame;

generating a contribution from each transmitter to the real aperture image using an aperture weighting function to generate a weighted sum of a component of that contribution from each receiver; and

generating the real aperture image for that down-track location by summing the contributions from each transmitter.

6. The method of claim 1 wherein the post-processing of a synthetic aperture radar image includes:

applying a filter to the image to generate a filtered image, the filtered image having an original spot of pixels;

designating that the pixels of the original spot are within an enlarged spot and that pixels in a region adjacent to the original spot that satisfy restoration criterion are also within the enlarged spot; and

setting values of the pixels based on whether the pixels are in the enlarged spot.

7. The method of claim 1 wherein the generating of the statistic for a down-track location includes:

calculating energy of pixels in the synthetic aperture image frame for the down-track location;

calculating a spot energy of pixels within the restored spot; and

setting the statistic based on a ratio of the spot energy to the energy of pixels.

8. The method of claim 1 wherein the surface is ground and a subsurface object is a buried explosive.

9. A computer-readable memory storing computer-executable instructions for controlling a computing device to detect presence of a subsurface object in a medium, by a method comprising:

for each of a plurality of down-track locations, acquiring return signals at receivers from signals emitted by transmitters into the medium, the transmitters and receiver forming a linear array of transceivers;

pre-processing the return signals to suppress effects of extraneous signals;

for each down-track location, generating a synthetic aperture radar image representing a subsurface image at that down-track location from the pre-processed return signals;

post-processing the synthetic aperture radar images to separate objects in the foreground from clutter and medium in the background by filtering out effects of clutter;

for each down-track location, generating a statistic based on energy of a spot within the filtered image for that down-track location; and

identifying down-track locations with peaks in the generated statistics; and

upon identifying a down-track location with a peak, indicating the presence of a subsurface object at that down-track location.

10. The computer-readable memory of claim 9 wherein the medium is ground and the method is for detecting presence of buried explosives.

11. The computer-readable memory of claim 9 wherein the return signals are acquired and the presence of a subsurface object is indicated in real-time.

12. The computer-readable memory of claim 9 wherein the pre-processing of current return signals includes:

calculating an estimated background signal for a current return signal based on a combination of an estimated background signal for a prior return signal for a prior down-track location that has a similar time delay between when the transmitter emits the signal and when a return signal that bounces off a surface is acquired by the receiver and an estimated gain and bias adjusted return signal for the prior return signal;

calculating a gain and bias adjusted return signal for the current return signal; and

calculating an estimated foreground signal for the current return signal based on a difference between the gain and bias adjusted return signal for the current return signal and the estimated background signal for the current return signal.

13. The computer-readable memory of claim 9 wherein the generating of a synthetic aperture radar image includes for each down-track location generating a real aperture radar image, the generating including:

performing a plane-to-plane propagation of the return signals to generate a contribution image frame;

generating a contribution from each transmitter to the real aperture image using an aperture weighting function to generate a weighted sum of a component of that contribution from each receiver; and

generating the real aperture image for that down-track location by summing the contributions from each transmitter.

14. The computer-readable memory of claim 9 wherein the post-processing of a synthetic aperture radar image includes:

applying a filter to the image to generate a filtered image, the filtered image having an original spot of pixels;

designating that the pixels of the original spot are within an enlarged spot and that pixels in a region adjacent to the original spot that satisfy restoration criterion are also within the enlarged spot; and

setting values of the pixels based on whether the pixels are in the enlarged spot.

15. The computer-readable memory of claim 9 wherein the generating of the statistic for a down-track location includes:

calculating energy of pixels in the synthetic aperture image frame for the down-track location;

calculating a spot energy of pixels within the restored spot; and

setting the statistic based on a ratio of the spot energy to the energy of pixels.

16. A computing device for detecting presence of a subsurface object in a medium from return signal acquired at a plurality of down-track locations by receivers from signals emitted by transmitters into the medium, the transmitters and receivers operating in a multistatic mode, comprising:

a component that pre-processes the return signals to suppress effects of extraneous signals;

a component that, for each down-track location, generates a real aperture radar image from the pre-processed return signals for that down-track location;

a component that, for each down-track location, generates a synthetic aperture radar image from the real aperture radar images of that down-track location and adjacent down-track locations;

a component that post-processes the synthetic aperture radar images to separate objects in the foreground from clutter and medium in the background by filtering out effects of clutter followed by restoring spots within the filtered images to their pre-filtered values;

a component that, for each down-track location, generates a statistic based on energy of a restored spot within the filtered image for that down-track location; and

a component that indicates the presence of a subsurface object at a down-track location having a statistic that represents a location peak in the statistics.

17. The computing device of claim 16 wherein the component that pre-processes the return signals:

calculates an estimated background signal for a current return signal based on a combination of an estimated background signal for a prior return signal for a prior down-track location that has a similar air gap to the current return signal and an estimated gain and bias adjusted return signal for the prior return signal;

calculates a gain and bias adjusted return signal for the current return signal; and

calculates an estimated foreground signal for the current return signal based on a difference between the gain and bias adjusted return signal for the current return signal and the estimated background signal for the current return signal.

18. The computing device of claim 17 wherein the component that generates a real aperture radar image for a down-track location:

performs a plane-to-plane propagation of the return signals to generate a contribution image frame;

generates a contribution from each transmitter to the real aperture image using an aperture weighting function to generate a weighted sum of a component of that contribution from each receiver; and

generates the real aperture image for that down-track location by summing the contributions from each transmitter.

19. The computing device of claim 18 wherein the component that post-processes a synthetic aperture radar image:

applies a filter to the image to generate a filtered image, the filtered image having an original spot of pixels;

designates that the pixels of the original spot are within an enlarged spot and that pixels in a region adjacent to the original spot that satisfy restoration criterion are also within the enlarged spot; and

sets values of the pixels based on whether the pixels are in the enlarged spot.

20. The computing device of claim 19 wherein the component that generates the statistic for a down-track location:

calculates energy of pixels in the synthetic aperture image frame for the down-track location;

calculates a spot energy of pixels within the restored spot; and

sets the statistic based on a ratio of the spot energy to the energy of pixels.

21. A computing device for detecting presence of an object in a medium from return signals acquired at a plurality of locations by receivers from signals emitted by transmitters into the medium, comprising:

a component that, for each location, generates a synthetic aperture radar image from return signals for that location and adjacent locations;

a component that processes the synthetic aperture radar images to filter out effects of clutter and restore spots within the filtered images to their pre-filtered values;

a component that, for each location, generates a statistic based on energy of a restored spot within the filtered image for that location; and

a component that indicates the presence of an object in the medium at a location having a statistic that represents a location peak in the statistics.

22. The computing device of claim 21 including a component that processes the return signals to suppress effects of extraneous signals.

23. The computing device of claim 22 wherein the component that processes the return signals:

calculates an estimated background signal for a current return signal based on a combination of an estimated background signal for a prior return signal for a prior location that has a similar air gap to the current return signal and an estimated gain and bias adjusted return signal for the prior return signal;

calculates a gain and bias adjusted return signal for the current return signal; and

calculates an estimated foreground signal for the current return signal based on a difference between the gain and bias adjusted return signal for the current return signal and the estimated background signal for the current return signal.

24. The computing device of claim 21 including a component that generates a real aperture radar image for a location by:

performing a plane-to-plane propagation of the return signals to generate a contribution image frame;

generating a contribution from each transmitter to the real aperture image using an aperture weighting function to generate a weighted sum of a component of that contribution from each receiver; and

generating the real aperture image for that location by summing the contributions from each transmitter.

25. The computing device of claim 21 wherein the component that processes a synthetic aperture radar image:

applies a filter to the image to generate a filtered image, the filtered image having an original spot of pixels;

designates that the pixels of the original spot are within an enlarged spot and that pixels in a region adjacent to the original spot that satisfy restoration criterion are also within the enlarged spot; and

sets values of the pixels based on whether the pixels are in the enlarged spot.

26. The computing device of claim 21 wherein the component that generates the statistic for a location:

calculates energy of pixels in the synthetic aperture image frame for the location;

calculates a spot energy of pixels within the restored spot; and

sets the statistic based on a ratio of the spot energy to the energy of pixels.