US20020122000A1 - Ground penetrating radar system - Google Patents
Ground penetrating radar system Download PDFInfo
- Publication number
- US20020122000A1 US20020122000A1 US09/752,085 US75208500A US2002122000A1 US 20020122000 A1 US20020122000 A1 US 20020122000A1 US 75208500 A US75208500 A US 75208500A US 2002122000 A1 US2002122000 A1 US 2002122000A1
- Authority
- US
- United States
- Prior art keywords
- transmit
- receive
- signal
- antenna
- antennas
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/04—Display arrangements
- G01S7/06—Cathode-ray tube displays or other two dimensional or three-dimensional displays
- G01S7/20—Stereoscopic displays; Three-dimensional displays; Pseudo-three-dimensional displays
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F41—WEAPONS
- F41H—ARMOUR; ARMOURED TURRETS; ARMOURED OR ARMED VEHICLES; MEANS OF ATTACK OR DEFENCE, e.g. CAMOUFLAGE, IN GENERAL
- F41H11/00—Defence installations; Defence devices
- F41H11/12—Means for clearing land minefields; Systems specially adapted for detection of landmines
- F41H11/13—Systems specially adapted for detection of landmines
- F41H11/136—Magnetic, electromagnetic, acoustic or radiation systems, e.g. ground penetrating radars or metal-detectors
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/35—Details of non-pulse systems
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/885—Radar or analogous systems specially adapted for specific applications for ground probing
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/40—Means for monitoring or calibrating
- G01S7/4052—Means for monitoring or calibrating by simulation of echoes
- G01S7/406—Means for monitoring or calibrating by simulation of echoes using internally generated reference signals, e.g. via delay line, via RF or IF signal injection or via integrated reference reflector or transponder
- G01S7/4073—Means for monitoring or calibrating by simulation of echoes using internally generated reference signals, e.g. via delay line, via RF or IF signal injection or via integrated reference reflector or transponder involving an IF signal injection
Definitions
- the invention relates to the field of radar systems. More particularly, the invention relates to the field of ground penetrating radar systems.
- the invention is a ground penetrating radar system mounted on a cart to achieve the desired mobility.
- the system uses two offset banks of interleaved transmit and receive antennas to achieve the desired accuracy.
- the receive and transmit antennas are properly oriented with respect to each other to reduce cross coupling and maximize desired subsurface echoes.
- the system uses nearfield beam forming, which is accomplished through fully coherent signal processing and synthetic aperture reception and processing, to image buried objects in three dimensions.
- the system displays a plan, or top, view and a side view of the area being scanned to provide a three dimensional perspective on a two dimensional computer screen.
- the invention features a ground penetrating radar system which includes a cart configured to be movable along the ground.
- a computer is mechanically coupled to the cart.
- a radar electronics module is mechanically coupled to the cart and electrically coupled to the computer.
- a first antenna array is mechanically coupled to the cart, electrically coupled to the radar electronics module, and oriented to radiate into the ground and receive radiation from the ground.
- a second antenna array is mechanically coupled to the cart, electrically coupled to the radar electronics module, and oriented to radiate into the ground and receive radiation from the ground.
- a movement detector which is configured to detect movement of the cart, is coupled to the computer.
- the computer is configured to trigger the radar electronics module when the computer detects that the cart has moved a predefined distance.
- the radar electronics module may include a first radar electronics module electrically coupled to the first antenna array and a second radar electronics module electrically coupled to the second antenna array.
- the first antenna array may be configured to radiate and receive radiation from a first series of points along a first set of curves parallel to the direction of movement of the cart.
- the second antenna array may be configured to radiate and receive radiation from a second series of points along a second set of curves parallel to the direction of movement of the cart.
- the first set of curves may be interleaved with the second set of curves.
- the invention features a ground penetrating radar system including a first bank of receive antennas arranged along a first axis, a first bank of transmit antennas arranged along a second axis substantially parallel to the first axis and horizontally displaced from the first axis, a second bank of receive antennas arranged along a third axis substantially parallel to the first axis and horizontally displaced from the first axis, and a second bank of transmit antennas arranged along a fourth axis substantially parallel to the first axis and horizontally displaced from the first axis.
- a first radar electronics module is coupled to the first bank of transmit antennas and the first bank of receive antennas.
- a second radar electronics module is coupled to the second bank of transmit antennas and the second bank of receive antennas.
- the transmit antennas in the first bank of transmit antennas are interleaved with the receive antennas in the first bank of receive antennas and the transmit antennas in the second bank of transmit antennas are interleaved with the receive antennas in the second bank of receive antennas.
- the receive antennas in the first bank of transmit antennas are offset along the first axis from the receive antennas in the second bank of transmit antennas.
- Implementations of the invention may include one or more of the following.
- the first bank of transmit antennas may be offset along the second axis with respect to the second bank of transmit antennas.
- the banks of receive antennas may alternate with the banks of transmit antennas.
- Each transmit antenna may be adjacent to at least one receive antenna.
- Each transmit antenna may be oriented to minimize electromagnetic coupling to at least one of its adjacent receive antennas.
- Each transmit antenna may include at least one spiral arm of conductive material.
- Each receive antenna may include at least one spiral arm of conductive material.
- a tangent to the inside of the spiral arm at the edge of a transmit antenna may be substantially perpendicular to a tangent to the inside of the spiral arm at the edge of a receive antenna adjacent to the transmit antenna.
- Each transmit antenna may include two spiral arms of conductive material.
- Each receive antenna may include two spiral arms of conductive material.
- the transmit antennas and the receive antennas may have faces with centers. Two adjacent first bank receive antennas from the first bank of receive antennas and a first bank transmit antenna from the first bank of transmit antennas interleaved between the two adjacent first bank receive antennas may be positioned such that lines between the centers of the faces of the two adjacent first bank receive antennas and the interleaved first bank transmit antenna form a first triangle having sides of approximately the same length.
- Two adjacent second bank receive antennas from the second bank of receive antennas and a second bank transmit antenna from the second bank of transmit antennas interleaved between the two adjacent second bank receive antennas may be positioned such that lines between the centers of the faces of the two adjacent second bank receive antennas and the interleaved second bank transmit antenna form a second triangle having sides of approximately the same length.
- a vertex of the first triangle may be displaced in the direction of the first axis relative to a corresponding vertex of the second triangle by an amount substantially equal to one-half the distance from the center of one side of the first triangle to the center of another side of the first triangle.
- the third axis may be horizontally displaced from the first axis by an amount substantially equal to eight times the distance from the center of one side of the first triangle to the center of another side of the first triangle.
- the transmit antennas may not be required to be in contact with the ground when in operation.
- the receive antennas may not be required to be in contact with the ground when in operation.
- the invention features a ground penetrating radar system including a digital module.
- the digital module includes a direct digital synthesizer configured to generate a digital IF reference signal.
- a digital to analog converter is coupled to the direct digital synthesizer and is configured to convert the digital IF reference signal to an analog IF transmit signal.
- An analog to digital converter is configured to convert an analog IF receive signal to a digital IF receive signal.
- a digital down converter is configured to digitally mix the digital IF receive signal with the digital IF reference signal to produce an in-phase product and the digital IF reference signal shifted in phase by ninety degrees to produce a quadrature product.
- the ground penetrating radar system includes an RF module coupled to the digital module.
- the RF module includes an up-converter configured to convert the analog IF transmit signal into a transmit signal and a down-converter configured to convert a receive signal into an analog IF receive signal.
- the system includes a transmit antenna array coupled to the up-converter for radiating the transmit signal and a receive antenna array coupled to the down-converter for receiving the receive signal.
- the transmit antenna array may include a plurality of transmit antennas.
- the receive antenna array may include a plurality of receive antennas.
- the system may include a digital signal processor.
- the system may include a transmit switch for applying the transmit signal to one of the plurality of transmit antennas.
- the transmit switch may be controlled by the digital signal processor.
- the system may include a receiver switch for receiving the receive signal from one of the plurality of receive antennas.
- the receiver switch may be controlled by the digital signal processor.
- the digital signal processor may control the direct digital synthesizer, the digital down converter, the up-converter and the down-converter.
- the transmit signal may be a stepped-frequency transmit signal.
- the receive signal may be a stepped-frequency receive signal.
- the system may include a computer coupled to a processor through an extensible network.
- the processor may be configured to communicate with the digital signal processor.
- the extensible network may be an Ethernet network.
- the invention features a ground penetrating radar system including a digital module configured to generate an analog IF transmit signal and to receive an analog IF receive signal.
- the system includes an RF module, which includes a triple-heterodyne up-converter for converting an analog IF transmit signal into a stepped-frequency transmit signal.
- the RF module also includes a triple-heterodyne frequency converter for converting a stepped-frequency receive signal into an analog IF receive signal.
- the system includes a transmit antenna coupled to the up-converter for radiating the stepped-frequency transmit signal and a receive antenna coupled to the down-converter for receiving the stepped-frequency receive signal.
- Implementations of the invention may include one or more of the following.
- the triple-heterodyne up-converter may include a first up-converter configured to mix the analog IF transmit signal with a signal from a first local oscillator to produce a first intermediate signal and an aliased first intermediate signal.
- the triple-heterodyne up-converter may include a first filter coupled to the first up-converter for substantially rejecting the aliased first intermediate signal.
- the triple-heterodyne up-converter may include a second up-converter coupled to the first filter configured to mix the first intermediate signal with a signal from a second local oscillator to produce a second intermediate signal and an aliased second intermediate signal.
- the triple-heterodyne up-converter may include a second filter coupled to the second up-converter for substantially rejecting the aliased second intermediate signal.
- the triple-heterodyne up-converter may include a down-converter coupled to the second filter configured to mix the second intermediate signal with a stepped frequency signal to produce the stepped-frequency transmit signal and an aliased stepped-frequency transmit signal.
- the stepped-frequency transmit signal may have substantially no frequency components in the pass bands of the first filter or the second filter.
- the triple-heterodyne up-converter may include a third filter coupled to the down-converter for substantially rejecting the aliased stepped-frequency transmit signal.
- the triple-heterodyne up converter may include an up-converter configured to mix the stepped-frequency receive signal with a stepped-frequency signal to produce a first intermediate signal and an aliased first intermediate signal.
- the triple-heterodyne up-converter may include a first filter coupled to the first up-converter for substantially rejecting the aliased first intermediate signal.
- the triple-heterodyne up-converter may include a first down-converter coupled to the first filter configured to mix the first intermediate signal with a signal from a first local oscillator to produce a second intermediate signal and an aliased second intermediate signal.
- the triple-heterodyne up-converter may include a second filter coupled to the first down-converter for substantially rejecting the aliased second intermediate signal.
- the triple-heterodyne up-converter may include a second down-converter coupled to the second filter configured to mix the second intermediate signal with a second local oscillator to produce the analog IF receive signal and an aliased analog IF receive signal.
- the triple-heterodyne up-converter may include a third filter coupled to the second down-converter for substantially rejecting the aliased analog IF receive signal.
- the invention features a ground penetrating radar system including a transmitter, a receiver, an array of transmit antennas, an array of receive antennas interleaved with the array of transmit antennas, a transmit switch configured to selectively couple the transmitter to one of the array of transmit antennas and a receive switch configured to selectively couple the receiver to one of the array of receive antennas.
- the array of transmit antennas is arranged in one or more rows.
- the array of receive antennas is arranged in one or more rows. Each row is parallel to, adjacent to and offset from one of the rows of transmit antennas, so that each receive antenna in a row except one is adjacent to two transmit antennas, and each transmit antenna in a row except one is adjacent to two receive antennas.
- the transmit switch and the receive switch are configured to couple the transmitter and receiver, respectively, to a first transmit antenna and a first adjacent receive antenna, and subsequently to the first transmit antenna and a second adjacent receive antenna.
- the invention features a method for collecting and displaying data from a ground penetrating radar system, which includes a plurality of transmit antennas and a plurality of receive antennas. Each transmit antenna, except one, has two adjacent receive antennas.
- the system is mounted on a movable cart.
- the method includes collecting raw data. Collecting raw data includes (a) selecting a first of the plurality of transmit antennas. Collecting raw data further includes (b) selecting a first receive antenna that is adjacent to the selected transmit antenna. Collecting raw data further includes (c) collecting data using the selected transmit antenna and the selected receive antenna to produce raw data.
- the raw data collected at spatial location (x m , v n ) is denoted by ⁇ tilde over ( ⁇ ) ⁇ mnp where the indices m, n are used to denote position in a grid of spatial locations where data has been collected, and p is an index ranging from 1 to P corresponding to the frequency f p at which the data was collected.
- Collecting raw data includes (d) repeating step (c) for both receive antennas adjacent to the selected transmit antenna.
- Collecting raw data further includes (e) repeating steps b, c and d for all transmit antennas.
- Collecting raw data further includes repeating steps a, b, c, d, and e each time the cart moves to a new location.
- the method further includes preconditioning the raw data to produce preconditioned data, analyzing the preconditioned data, and displaying images of the analyzed data.
- Implementations of the invention may include one or more of the following.
- Preconditioning the raw data to produce preconditioned data may include (g) removing a constant frequency component and a system travel time delay, (h) removing a transmit-antenna to receive-antenna coupling effect, (i) prewhitening, and (h) repeating steps (g), (h) and (i) for each spatial location of the raw data.
- Removing the transmit-antenna to receive-antenna coupling effect may include applying the following equation:
- ⁇ mnp ⁇ circumflex over ( ⁇ ) ⁇ mnp ⁇ circumflex over ( ⁇ ) ⁇ m ⁇ p
- N 1 and N 2 define a region to be imaged.
- Removing the transmit-antenna to receive-antenna coupling effect may include applying the following equation:
- ⁇ mnp ⁇ circumflex over ( ⁇ ) ⁇ mnp ⁇ circumflex over ( ⁇ ) ⁇ ⁇ tilde over (m) ⁇ np
- Prewhitening may include applying the following equation:
- I mnw the complex image value at spatial location (X F,m , Y F,n , Z w );
- U is the SAR array size in the cross-track direction
- V is SAR array size in the along track direction
- (f p 1 , f p 2 ) is the frequency processing band
- ⁇ uvw is the travel time from source (u,v) in the SAR array down to a focal point at depth z w and back up to receiver (u, v) in the SAR array;
- x F , m 3 ⁇ d 4 + ( m - 1 ) ⁇ d 4 ;
- dy scan spacing
- Displaying images of the analyzed data may include computing a plan view image of the analyzed data, computing a side view image of the analyzed data, and displaying the plan view image and the side view image.
- Computing a plan view image of the analyzed data may include applying the following equation:
- PlanView min max w
- max w is the maximum value across all w (depths).
- Computing a side view image of the analyzed data may include applying the following equation:
- max w is the maximum value across all w (depths).
- FIG. 1 is a perspective view of a ground penetrating radar system according to the present invention.
- FIG. 2 is a block diagram of a ground penetrating radar system.
- FIG. 3 illustrates the grid of radar data collected by a ground penetrating radar system.
- FIG. 4 illustrates the geometry of antenna arrays incorporated in the ground penetrating radar system.
- FIG. 5 illustrates the switching among the antennas within the antenna arrays.
- FIG. 6 illustrates the configuration of the antennas.
- FIG. 7 illustrates the relative orientation of a transmit antenna and a receive antenna.
- FIG. 8 is a block diagram of the electronics in a ground penetrating radar system.
- FIG. 9 is a block diagram of a digital module.
- FIG. 10 is a block diagram of an embodiment of the signal synthesis and analysis module within a digital module.
- FIG. 11 is a block diagram of a digital down-converter.
- FIG. 12A is a block diagram of an embodiment of the signal synthesis and analysis module within a digital module.
- FIG. 12B is a block diagram of a field programmable gate array.
- FIG. 13 is a block diagram of the RF module.
- FIG. 14 is a frequency plan for a preferred embodiment of the ground penetrating radar system.
- FIG. 15 illustrates a stepped-frequency signal.
- FIG. 16 is a block diagram of the first stage in a triple-heterodyne frequency up-converter and the last stage in a triple-heterodyne frequency converter.
- FIG. 17 is a block diagram of a middle stage in a triple-heterodyne frequency up-converter and the middle stage in a triple-heterodyne down-converter.
- FIG. 18 illustrates the last stage in a triple-heterodyne frequency up-converter and the first stage in a triple-heterodyne frequency converter.
- FIG. 19 is a block diagram of amplifiers and filters in the triple-heterodyne up-converter and the triple-heterodyne frequency converter.
- FIG. 20 is a block diagram of a local oscillator.
- FIG. 21 is a block diagram of the transmit switch and the receive switch.
- FIG. 22 illustrates the grid of points collected by movement of the ground penetrating radar system along the ground.
- FIG. 23 is a block diagram of the raw data collection processing.
- FIG. 24 is a block diagram of the raw data analysis processing.
- FIG. 25 illustrates the geometry of depth focusing.
- FIG. 26 illustrates the geometry of collection of points midway between a transmit antenna and a receive antenna.
- FIG. 27 illustrates the geometry of synthetic aperture radar processing.
- FIG. 28 illustrates a volume of collected data and a top view and a side view of that data.
- FIG. 29 illustrates images of data prior to processing.
- FIG. 30 illustrates images of data after processing.
- FIG. 31 illustrates a physical block diagram of the processing performed in the ground penetrating radar system.
- a system according to the invention includes interleaved antenna arrays of properly oriented spiral transmit and receive antennas. Switching among the antennas allows samples to be efficiently taken with high spatial density.
- a digital frequency synthesizer and digital down converter allow fast and accurate measurements. Triple heterodyne up and down conversion reduces the likelihood that mixing products will interfere with the measurements.
- the system is extensible because of an Ethernet interconnection between system components. The system uses advanced processing and spatial filtering techniques to improve the quality of the images produced.
- a ground penetrating radar system 102 illustrated in FIG. 1, includes a cart 104 that can be moved along the ground.
- the cart 104 has four all terrain wheels 106 that give the cart 104 mobility over all types of terrain.
- Other techniques for providing mobility, such as tracks or skids, are within the scope of the invention.
- the cart 104 is a lightweight, compact design fabricated out of aluminum tubing. Four inch square, thin wall tubing is used for the frame. All joints are welded and the frame is painted with polyurethane paint. Mount points for all major system components are built into the frame with easy access to each mount for individual component removal without tools.
- the cart is fabricated out of other metals, such as steel, or other materials, such as light-weight composites.
- the system 102 includes a computer 108 mechanically coupled to the cart 104 .
- the computer 108 is a laptop computer, having its own display, memory, long term storage and other peripherals, that rests on or is secured to the cart 104 .
- the computer is a desktop-type computer.
- a radar electronics module 110 , 112 is mechanically coupled to the cart 104 and electrically coupled to the computer 108 .
- the radar electronics module 110 , 112 includes two 25 inch wide, 12 inch long, 2.5 inch high black anodized aluminum enclosures 110 , 112 , each containing two separate compartments for a digital circuit board and an analog (or RF) printed circuit board (both described in detail below).
- the circuit board compartments are separated by an aluminum plate, which physically and electrically isolates the boards.
- the barrier plate eliminates any potential for EMI and RFI noise between the two boards.
- the radar module enclosure cavities are accessible by removing the top cover plate.
- a first (or front) antenna array 114 is mechanically coupled to the cart 104 and electrically coupled to the radar electronics module 110 , 112 .
- the first antenna array 114 is oriented to radiate into the ground and receive radiation from the ground.
- a second (or back) antenna array 116 is mechanically coupled to the cart 104 and electrically coupled to the radar electronics module 110 , 112 .
- the second antenna array 116 is oriented to radiate into the ground and receive radiation from the ground.
- the two antenna arrays 114 and 116 are coupled to the cart 104 by an antenna framework 118 .
- the antenna arrays 114 and 116 and the antenna framework 118 can be removed from the cart 104 for diagnostics and transport.
- the height of the antenna arrays 114 and 116 above the ground may be adjusted by a machine screw jack 120 to be, in a preferred embodiment, between six inches and twenty inches.
- a drill motor 122 drives the machine screw jack 120 to raise and lower the antenna arrays.
- the drill motor 122 could be replaced by any type of motor that has sufficient power to raise and lower the antenna arrays.
- the antenna framework 118 offsets the antenna banks (in a preferred embodiment by 24 inches) from the cart 104 frame to eliminate any possibility of target shadowing by the frame.
- a movement detector 202 shown in FIG. 2, coupled to the computer detects movement of the cart 104 .
- the movement detector 202 is an optical encoder mounted on the front right wheel of the cart 104 .
- the computer 108 triggers the radar electronics module 110 , 112 when the computer 108 detects that the cart 104 has moved a predefined distance. In this way, data can be acquired at prescribed intervals measured by the motion detector 202 .
- the radar electronics module 110 , 112 includes two electronics modules.
- a first radar electronics module 110 is coupled to the first antenna array 114 .
- a second radar electronics module 112 is coupled to the second antenna array 116 .
- the first antenna array 114 radiates and receives radiation from a first series of points, e.g. 302 shown in FIG. 3, along a first set of curves 304 parallel to the direction of movement 306 of the cart 104 .
- the open circles represent the locations of the centers of the transmit and receive antennas when they take a radar sample, as described below.
- the closed, or filled, circles represent the location of the samples.
- the curves in the first set of curves 304 are approximately d/2, or 2.76 inches, apart.
- the second antenna array 116 radiates and receives radiation from a second series of points, e.g. 308 , along a second set of curves 310 parallel to the direction of movement 306 of the cart 104 .
- the second antenna array 116 is shown a large distance behind the first antenna array 114 .
- the two arrays are preferably arranged in much closer proximity to each other, as shown in FIGS. 1 and 2 and as discussed below in the discussion of FIG. 4.
- the curves in the second set of curves 310 are approximately d/2, or 2.76 inches, apart.
- the first set of curves 304 is interleaved with the second set of curves 310 .
- the result is a set of curves including the first set of curves 304 and the second set of curves 310 , with the curves being approximately d/4, or 1.38, inches apart.
- the parameter d in FIG. 3 denotes the outside diameter of an antenna.
- d-5.515 inches The individual transmit and receive antennas are of the log spiral type and were designed to operate over the band 800-4000 MHz but actually radiate well down to 500 MHz.
- the antenna cavity is nine inches tall and is filled with radar absorbing material to minimize reflections from the top of the antenna.
- the transmit and receive antennas have different windings in order to minimize cross coupling during the simultaneous transmission and reception required by a stepped-frequency radar.
- the first antenna array 114 includes a first bank of receive antennas 402 arranged along a first axis 404 and a first bank of transmit antennas 406 arranged along a second axis 408 substantially parallel to the first axis and horizontally displaced from the first axis, as shown in FIG. 4.
- the second antenna array 116 includes a second bank of receive antennas 410 arranged along a third axis 412 substantially parallel to the first axis 404 and horizontally displaced from the first axis 404 and a second bank of transmit antennas 414 arranged along a fourth axis 416 substantially parallel to the first axis 404 and horizontally displaced from the first axis 404 .
- the first radar electronics module 110 is coupled to the first bank of receive antennas 402 and the first bank of transmit antennas 406 .
- the second radar electronics module 112 is coupled to the second bank of receive antennas 410 and the second bank of transmit antennas 414 .
- the receive antennas in the first bank of receive antennas 402 are interleaved with the transmit antennas in the first bank of transmit antennas 406 .
- the receive antennas in the second bank of receive antennas 410 are interleaved with the transmit antennas in the second bank of transmit antennas 414 .
- the receive antennas in the first bank of receive antennas 402 are offset along the first axis 404 from the receive antennas in the second bank of receive antennas.
- the first bank of transmit antennas 406 is offset with respect to the second bank of transmit antennas 414 .
- the banks of transmit antennas 406 , 414 alternate with the banks of receive antennas 402 , 410 .
- Each receive antenna e.g. 418
- Each receive antenna, e.g. 418 is oriented to minimize electromagnetic coupling with at least one of its adjacent transmit antennas, e.g. 420 , 422 .
- the close spacing of the curves 304 and 310 is accomplished through the interleaving of the transmit and receive banks of antennas, as discussed above, and by sharing the transmit and receive antennas within the same bank.
- Each transmit antenna transmits to its two adjacent receive antennas, as shown in FIG. 5.
- the transmit/receive sequence for the first antenna array is TX 0 -RX 0 , TX 0 -RX 1 , TX 1 -RX 1 , TX 1 -RX 2 , . . ., TX 6 -RX 6 .
- the transmit/receive sequence for the second antenna array is similar.
- a sample of radar data is taken.
- a transmitter/receiver pair such as TX 0 -RX 0
- an array of data may be acquired. Further, samples at a number of frequencies are taken at each sample location.
- y(m,n,p) is the complex frequency response of the ground at the location (m,n). It is the quantity measured by the stepped-frequency radar.
- the Fourier transform of y(m,n,p) is the equivalent time domain response of the radar at location (m,n).
- Each transmit antenna 602 illustrated in FIG. 6, includes at least one spiral arm 604 of conductive material.
- each receive antenna e.g. 606 , includes at least one spiral arm 608 of conductive material.
- the transmit antenna 602 and one of its adjacent receive antennas 606 are oriented so that a tangent 702 to the inside of the spiral arm 604 at the edge of a transmit antenna 602 is substantially perpendicular to a tangent 704 to the inside of the spiral arm 608 at the edge of a receive antenna 606 adjacent to the transmit antenna. This orientation minimizes electromagnetic cross coupling between the transmit antenna 602 and the receive antenna 606 .
- each transmit antenna e.g. 602
- each receive antenna 606 includes two spiral arms 608 , 612 of conductive material.
- the transmit antennas, e.g. 424 and 426 , and the receive antennas, e.g. 428 have faces. Each of the faces has a center.
- Two adjacent first bank receive antennas 424 , 426 from the first bank of receive antennas 402 and a first bank transmit antenna 428 from the first bank of transmit antennas 406 interleaved between the two adjacent first bank receive antennas 424 , 426 are positioned such that lines between the centers of the faces of the two adjacent first bank receive antennas 424 , 426 and the interleaved first bank transmit antenna 428 form a first triangle 430 having sides of approximately the same length.
- Two adjacent second bank receive antennas 432 , 434 from the second bank of receive antennas 410 and a second bank transmit antenna 436 from the second bank of transmit antennas 414 interleaved between the two adjacent second bank receive antennas 432 , 434 are positioned such that lines between the centers of the faces of the two adjacent second bank receive antennas 432 , 434 and the interleaved second bank transmit antenna 436 form a second triangle 438 having sides of approximately the same length.
- a vertex 440 of the first triangle 430 is displaced in the direction of the first axis 404 relative to a corresponding vertex 442 of the second triangle 438 by an amount substantially equal to one-half the distance from the center of one side of the first triangle 432 to the center of another side of the first triangle 432 .
- the third axis 412 is horizontally displaced from the first axis 404 by an amount substantially equal to eight times the distance from the center of one side of the first triangle 432 to the center of another side of the first triangle 432 .
- neither the transmit antennas, e.g. 420 , nor the receive antennas, e.g. 422 , are required to be in contact with the ground when in operation.
- the ground penetrating radar system operates in a continuous wave mode.
- the first and second radar electronics modules 110 , 112 generate a frequency-stepped radar signal and receive and analyze the return signal.
- the first and second radar electronics modules 110 , 112 in the ground penetrating radar system each include a digital module 802 , 804 coupled to the computer 108 , and an RF module 806 , 808 coupled to the digital module 802 , 804 , as shown in FIG. 8.
- the digital module 802 , 804 generates an IF signal, which in the preferred embodiment has a frequency of 10.7 MHz.
- the RF module 806 , 808 under the control of the digital module 802 , 804 , converts the IF signal to a stepped-frequency signal, RFTX, which is provided to a transmit switch 810 .
- the transmit switch 810 , 812 provides the RFTX signal to one transmit antenna in its respective first or second antenna array 114 , 116 which radiates a signal into the ground.
- a stepped frequency return signal is received by the antenna array 114 , 116 and routed to a receiver switch 814 , 816 .
- the receiver switch 814 , 816 selects a receive antenna from which to receive the returned signal and routes the signal to the RF module 806 , 808 .
- the RF module 806 , 808 under the control of the digital module 802 , 804 , converts the stepped frequency return signal into a receive IF signal, which in the preferred embodiment has a frequency of 10.7 MHz.
- the digital module 802 , 804 demodulates the IF signal and provides the result to the computer 108 for storage and processing.
- the digital module 802 , 804 illustrated in more detail in FIG. 9, includes an extensible network interface 902 , which in the preferred embodiment is a Ethernet interface, and a point-to-point communication interface 904 , which in the preferred embodiment is an RS-232 interface. Both the extensible network interface 902 and the point-to-point communication interface 904 allow communication between the computer 108 and a processor 906 within the digital module 802 , 804 for communicating configuration commands and data from the computer 108 to the digital module 802 , 804 and for communicating status and collected and processed data from the digital module 802 , 804 to the computer 108 . Data and program code for the processor 906 are stored within the digital module on a memory 910 .
- a signal synthesis and analysis module 908 synthesizes the transmit IF signal, controls the transmitter switch 810 , 812 , controls the up- and down-conversion performed in the RF module 806 , 808 , demodulates the receive IF signal, and controls the receiver switch 814 , 816 .
- the signal synthesis and analysis module 908 includes a direct digital synthesizer 1002 , which generates a digital IF reference signal 1004 .
- a digital to analog converter 1006 converts the digital IF signal to an analog IF transmit signal 1008 .
- the analog IF transmit signal 1008 is a 10.7 MHz IF signal.
- the direct digital synthesizer 1002 and the digital to analog converter 1006 are incorporated into a single module 1010 , such as the AD7008 manufactured by Analog Devices.
- the AD7008 includes a 10-bit analog to digital converter.
- the digital IF reference signal 1004 provides a stable and accurate reference for both the transmit and receive paths of the apparatus. As will be seen, the digital IF reference signal 1004 is used in the down conversion process without the necessity of applying delays.
- the signal synthesis and analysis module 908 also includes an analog to digital converter 1012 which digitizes the analog IF receive signal 1014 to produce a digital IF receive signal 1016 .
- the digital IF receive signal 1016 is routed to a digital down converter 1018 which digitally mixes the digital IF receive signal 1016 with the digital IF reference signal 1004 from the direct digital synthesizer 1002 through a programmable logic device (PLD) 1020 .
- PLD 1020 provides control logic and interfacing for the direct digital synthesizer 1002 and the digital down converter 1018 and a digital signal processor (DSP) 1022 .
- the DSP 1022 controls the direct digital synthesizer 1002 , the digital down converter 1018 , and the other peripherals in the system.
- the PLD also receives an Out of Range signal from the analog to digital converter 1012 , which it communicates to the DSP 1022 .
- the DSP 1022 reacts by changing the gain in the receive chain, as discussed below.
- the digital down converter 1018 digitally mixes 1102 the digital IF reference signal 1004 with the digital IF receive signal 1016 to produce an in-phase product 1104 , as shown in FIG. 11.
- the in-phase product 1104 is then digitally filtered by a low pass filter 1106 to produce a filtered in phase component 1108 which is transferred to the DSP 1022 for processing and storage.
- the digital down converter 1018 also digitally shifts the phase of the digital IF reference signal by ninety degrees 1110 and digitally mixes 1112 it with the digital IF receive signal 1016 to produce a quadrature product 1114 .
- the quadrature product 1114 is then digitally filtered through a low pass filter 1116 to produce a filtered quadrature component 1118 which is passed to the DSP 1022 for processing and storage.
- the low pass filters 1106 and 1116 perform their functions by integrating a moving window of input data.
- the window is 128 signals wide.
- the signal synthesis and analysis module 908 includes a memory 1024 where the DSP 1022 and the PLD 1020 store data.
- the signal synthesis and analysis module 908 also includes a program memory 1026 where program code for the DSP 1022 is stored.
- FIG. 12A Another example embodiment of the digital module 908 is illustrated in FIG. 12A.
- the DSP 1022 , the memory 1024 and the program memory 1026 remain in the same configuration as shown in FIG. 10.
- the data and programs stored in the memory 1024 and the program memory 1026 may be different from the data and programs stored in the same elements in FIG. 10.
- a PLD 1202 provides a control logic interface between the DSP 1022 and a field programmable gate array (FPGA) 1204 .
- FPGA field programmable gate array
- a memory 1206 preferably a programmable read only memory (PROM), stores the FPGA boot program.
- PROM programmable read only memory
- a clock 1208 preferably an 85.6 MHz clock, provides a clock signal to the FPGA, a digital to analog converter (DAC) 1210 and an analog to digital converter (ADC) 1212 .
- the FPGA 1204 is controlled by the PLD 1202 , which in turn is controlled by the DSP 1022 , through control lines 1214 .
- the FPGA 1204 generates a digital IF reference signal 1216 .
- the DAC 1210 converts the digital IF reference signal 1216 to an analog signal which is isolated by transformer 1218 and provided as the analog IF transmit signal 1008 .
- the analog IF receive signal 1014 is converted by the ADC 1212 to a digital receive signal 1220 which is provided to the FPGA 1204 .
- the ADC 1212 also produces an out of range signal 1222 if the analog IF receive signal 1014 exceeds the range of the ADC 1212 .
- the FPGA 1204 produces received data 1224 which is transferred to the DSP 1022 for processing.
- the digital receive signal 1220 is provided to two digital mixers 1226 and 1228 .
- the first digital mixer 1226 mixes the digital receive signal 1220 with a digital cosine wave 1230 having a frequency, in one preferred embodiment, of 10.7 MHZ.
- the output of the first digital mixer 1226 is the real component 1232 of the digital receive signal 1220 .
- the real component 1232 is filtered 1234 , preferably by summing a 128 sample moving window of data, to produce a filtered real component signal 1236 which is provided to a multiplexer 1238 .
- the second digital mixer 1228 mixes the digital receive signal 1220 with a digital sine wave 1240 having a frequency, in one preferred embodiment, of 10.7 MHZ.
- the digital sine wave 1240 is identical to the digital cosine wave 1230 except that the digital sine wave 1240 is shifted ninety degrees in phase with respect to the digital cosine wave 1230 .
- the output of the first digital mixer 1228 is the imaginary component 1242 of the digital receive signal 1220 .
- the imaginary component 1242 is filtered 1244 , preferably by summing a 128 sample moving window of data, to produce a filtered imaginary component signal 1246 which is provided to the multiplexer 1238 .
- the digital sine wave 1240 is also provided as the digital IF reference signal 1216 .
- the multiplexer 1238 provides to the DSP 1022 either the filtered real component signal 1236 or the filtered imaginary component signal 1246 , depending on the select signal 1248 .
- the RF module 806 , 808 shown in more detail in FIG. 13, includes a triple-heterodyne up-converter 1302 for converting the analog IF transmit signal 1008 into a stepped-frequency transmit signal 1304 .
- the RF module 806 , 808 also includes a triple-heterodyne frequency converter 1306 for converting a stepped frequency receive signal 1308 into the analog IF receive signal 1014 .
- the triple-heterodyne up-converter 1302 includes a first up-converter 1310 .
- the first up-converter 1310 includes a mixer 1602 , as shown in FIG. 16, which mixes the analog IF transmit signal 1008 (after it has been isolated by transformer 1604 and filtered by low pass filter 1606 ) with the signal from a first local oscillator 1312 to produce a first intermediate signal 1314 and an aliased first intermediate signal 1402 , as shown in FIG. 14. These two signals are isolated by transformer 1608 and amplified by amplifiers 1610 , 1612 and 1613 .
- a filter 1614 having a pass band 1404 , shown in FIG. 14, substantially rejects the aliased first intermediate signal 1402 .
- the first local oscillator 1312 operates at a frequency of 122 MHZ and when mixed with the 10.7 MHZ IF produces mixing products at 111.3 MHZ and 132.7 MHZ.
- the filter 1614 has a center frequency of 139.75 MHZ and passes the 132.7 MHZ mixing product.
- a second up-converter 1316 includes a mixer 1702 , as shown in FIG. 17, which mixes the first intermediate signal 1314 with the signal produced by the second local oscillator 1318 to produce a second intermediate signal 1320 and an aliased second intermediate signal 1406 , as shown in FIG. 14.
- a filter 1704 substantially rejects the aliased second intermediate signal 1406 .
- the second local oscillator 1318 operates at a frequency of 2275 MHZ and when mixed with the first intermediate signal 1314 produces mixing products at 2142.3 MHZ and 2407.7 MHZ.
- the filter 1704 has a center frequency of 2450 MHZ and passes the 2407.7 MHZ product.
- a down-converter 1322 includes a mixer 1802 , shown in FIG. 18, that mixes the second intermediate signal 1320 (after being amplified by amplifier 1804 ) with a stepped frequency signal 1324 generated by a synthesizer 1326 to produce the stepped-frequency transmit signal 1328 and an aliased stepped-frequency transmit signal (not shown in FIG. 14).
- the stepped-frequency transmit signal has substantially no frequency components in the pass bands of the first filter 1614 or the second filter 1704 .
- a third filter 1902 shown in FIG. 19, substantially rejects the aliased stepped-frequency transmit signal (after it is amplified by amplifier 1904 ).
- Amplifier 1906 amplifies the output of the third filter 1902 to produce the stepped-frequency transmit signal 1304 .
- the third filter 1902 is a low pass filter with a 3 dB break point at 2000 MHZ.
- the synthesizer 1326 is controlled by the DSP 1022 through the LO Control and Configuration signals shown on FIG. 13 to produce a stepped-frequency signal 1324 , illustrated in FIG. 15.
- the frequency of the stepped-frequency signal 1324 steps through a range of frequencies 1502 under the control of the DSP 1022 .
- the range of frequencies 1502 is from 500 to 2000 MHZ.
- the stepped-frequency signal 1324 is mixed with the second intermediate signal 1322 the result is a preamplified stepped-frequency transmit signal 1328 ranging from 407.7 MHZ to 1907.7 MHZ.
- the triple-heterodyne up-converter 1302 includes amplifiers and filters 1330 which amplify and filter the preamplified stepped-frequency transmit signal 1328 before it is sent to the transmit switch 810 , 812 and then to the antenna array 114 , 116 as stepped-frequency transmit signal 1304 .
- the triple-heterodyne up converter 1306 includes an amplifier/filter 1332 that amplifies, using amplifier 1908 , and filters, using filter 1910 , the stepped-frequency receive signal 1308 to produce an amplified stepped-frequency received signal 1334 , as illustrated in FIG. 19.
- An up-converter 1336 uses mixer 1806 to mix the stepped-frequency receive signal 1334 (after amplification by amplifier 1808 ) with the stepped-frequency signal 1324 output of the synthesizer 1326 to produce a first intermediate signal 1338 and an aliased first intermediate signal 1410 , as shown in FIG. 14.
- a first filter 1810 substantially rejects the aliased first intermediate signal 1410 (after it has been amplified by amplifier 1812 ).
- a first down-converter 1340 uses a mixer 1706 to mix the first intermediate signal 1338 with the signal produced by the second local oscillator 1318 to produce a second intermediate signal 1342 and an aliased second intermediate signal 1412 .
- a second filter 1708 substantially rejects the aliased second intermediate signal 1412 (after it is amplified by amplifier 1710 ).
- a second down-converter 1344 uses a mixer 1616 to mix the second intermediate signal 1342 (after it is amplified by amplifier 1618 ) with the signal produced by the first local oscillator 1312 to produce the analog IF receive signal 1014 and an aliased analog IF receive signal (not shown).
- a third filter 1620 substantially rejects the aliased analog IF receive signal.
- An variable-gain amplifier 1622 between the mixer 1616 and the third filter 1620 allows the DSP to control the receive gain using gain control signals 1624 .
- a transformer 1626 isolates the mixer 1616 from the amplifier 1624 .
- An amplifier 1628 amplifies the output of the third filter 1620 and a transformer 1630 isolates the output of the amplifier 1628 from the output 1014 .
- the second local oscillator shown in detail in FIG. 17, is a phase locked loop including a phase locked loop (PLL) frequency synthesizer 1712 .
- PLL frequency synthesizer 1712 is a MC145202 manufactured by Motorola.
- the PLL frequency synthesizer 1712 has an SDI input, a SCLK input and an SLD input from the DSP 1022 .
- the SDI input is the serial data line from the DSP 1022 that provides the configuration data for the PLL frequency synthesizer 1712 .
- the DSP 1022 can specify the frequency at which the phase locked loop 1318 is to operate.
- the SDI signal is also used to transfer data into the local oscillator 1326 , as discussed below.
- the SCLK is used to clock the transfer of serial data into the PLL frequency synthesizer 1712 .
- the SCLK signal is also used to clock the transfer of data into the LO, as described below.
- the SLD signal is activated when the SDI data is intended for the PLL frequency synthesizer 1712 .
- the PLL frequency synthesizer 1712 provides a loop error signal 1716 as an input to the voltage controlled oscillator (“VCO”) 1714 .
- the frequency of the output 1718 of the VCO 1714 varies with the loop error signal 1716 .
- the VCO output is amplified by amplifier 1720 and provided as the output of local oscillator 1318 through amplifiers 1722 and 1724 .
- the output of amplifier 1720 is frequency divided by divider 1726 to produce an input to the PLL frequency synthesizer 1712 , thereby closing the phase locked loop.
- phase locked loop 1318 is replaced by other apparatus for generating a stable frequency signal, such as a crystal.
- the local oscillator 1326 shown in more detail in FIG. 18, includes an LO synthesizer 1814 , which generates the stepped frequency signal 1324 .
- Amplifiers 1816 and 1818 amplify the stepped frequency signal 1324 for use by the down converter 1322 and the up converter 1336 .
- the LO synthesizer 1814 is controlled by the DSP using control lines 1820 and produces a calibration signal PN 1822 , which is used by the DSP to calibrate the LO.
- the LO synthesizer 1814 illustrated in detail in FIG. 20, includes a digital to analog converter (DAC) 2002 which is provided a voltage reference 2004 .
- the level of the analog output of the DAC 2002 is determined by the SDI, SCLK and SLD inputs from the DSP. These signals were described above in the description of the phase locked loop 1318 .
- the SCLK signal clocks data in the SDI signal into the DAC 2002 .
- the data clocked into the DAC determines the level of its output.
- the DAC output is amplified by an amplifier 2005 and used to drive a voltage controlled oscillator (VCO) 2006 .
- VCO voltage controlled oscillator
- the DSP 1022 can control the frequency of the output of the VCO 2006 .
- the DSP 1022 controls the frequency of the output of the VCO 2006 to vary in steps, as shown in FIG. 15, from 1000 to 2000 MHZ.
- An amplifier 2008 amplifies the output of the VCO 2006 .
- the system generates the entire desired range of frequencies, which, in a preferred embodiment is from 500 to 2000 MHZ, by using a series of switches.
- the DSP unasserts the SHI signal which causes an input switch 2010 to switch to its “lo” position and an output switch 2012 to switch to its “lo” position.
- This configuration of switches causes the amplified output of the VCO 2006 to be routed to a frequency divider 2104 , which divides the frequency of the signal by two. The result is filtered by low pass filter 2016 and provided as the stepped frequency signal 1324 .
- the DSP 1022 can generate the higher range of frequencies (e.g., from 1000 to 2000 MHZ) by asserting the SHI signal, which causes the input switch 2010 and the output switch 2012 to switch to their respective “HI” positions.
- switch 2018 is also closed allowing the amplified output of the VCO 2006 to be filtered by a low pass filter and provided as the stepped frequency signal 1324 .
- the output of the VCO 2006 is also used to generate a calibration signal PN for use by the DSP 1022 .
- the output of the VCO 2006 is amplified by amplifier 2022 and its frequency is divided by 64 by using a series of dividers 2024 , 2026 , and 2028 .
- a comparator 2030 squares the edges of the divided signal.
- Two PALs 2032 and 2034 provide a logic interface between the divided signal and the DSP 1022 .
- the transmit switch 810 , 812 illustrated in FIG. 21, includes seven single-pole, single-throw switches 2102 with their common poles wired in parallel and connected to the stepped-frequency transmit signal 1304 .
- the other pole of each switch is connected to one of the transmit antennas. Therefore, when one of the switches 2102 is actuated, the stepped frequency transmit signal 1304 is connected to a respective transmit antenna.
- the selection of which of the switches 2102 to close is controlled by the DSP 1022 through the transmitter switch control lines, which include control lines TXCTL 1 , TXCTL 2 , and TXCTL 3 .
- the binary combination of these control lines determine which of the transmitter switches 2102 will close.
- the receiver switch 814 , 816 also includes seven single-pole, single-throw switches 2104 with their common poles wired in parallel and connected to the stepped frequency received signal 1308 .
- the other pole of each switch is connected to one of the received antennas in the receive array. Therefore, when one of the switches 2104 is closed, the corresponding receive antenna provides the stepped frequency received signal 1308 .
- the determination of which of the switches 2104 is to close is controlled by the DSP 1022 through the receive switch control signals, which includes three switch control signals RXCTL 1 , RXCTL 2 , and RXCTL 3 .
- the switch 2104 that is closed depends on the binary combination of those three control signals.
- the system records data at 13 locations spaced 2.76 inches apart for each antenna array 114 , 116 .
- Software in the computer 108 fuses this data into a grid of rows and columns of data, where each row includes data taken from 26 locations spaced 1.38 inches apart. The columns represent the data taken at different locations along the direction of cart movement 306 (see FIG. 3).
- dx is 1.38 inches and dy is operator selectable, and is preferably selected to be 1.38 inches, producing a rectangular array.
- m varies from 1 to 26 and n varies from 1 to N, where N is the total number of scans that the system performs.
- the cross track resolution when both antenna arrays are used is 1.38′′. Individual images can be formed using either the first antenna array 114 alone or the second antenna array alone. For individual images of these two types, the cross track resolution is 2.76′′.
- the processing that the system performs includes collecting raw data and analyzing the raw data. Collecting the raw data, illustrated in FIG. 23, includes detecting movement of the cart to the next data collection position (block 2302 ). The system accomplishes this by monitoring the movement detector 202 . The system then selects one of the transmit antennas (block 2304 ) and determines if all of the transmit antennas have been processed (block 2306 ). If they have, the process terminates (block 2308 ) until the cart moves to the next data collection position (block 2302 ).
- the system selects a receive antenna adjacent to the selected transmit antenna (block 2310 ). The system then determines if both of the receive antennas adjacent to the selected transmit antenna have been processed (block 2312 ). If they have, the system selects the next transmit antenna (block 2314 ) and returns to the beginning of the loop (to block 2306 ).
- the system collects data using the selected transmit antenna and the selected receive antenna (block 2316 ) to produce raw data 2318 , the raw data collected at spatial location (x m ,y n ,) being denoted by ⁇ tilde over ( ⁇ ) ⁇ mnp where the indices m, n are used to denote position in the grid of spatial locations where data has been collected, and p is an index ranging from 1 to P corresponding to the frequency f p at which the data was collected.
- the unambiguous range of the radar is c/(2 ⁇ df) where c denotes light speed in air.
- Each of the data points ⁇ tilde over ( ⁇ ) ⁇ mnp is a complex number. In a preferred embodiment, a user may specify f 1 and f p .
- the system selects the next receive antenna (block 2320 ) and returns to the beginning of the data collection loop (block 2312 ).
- Analyzing the raw data 2318 includes preconditioning the raw data 2318 to produce preconditioned data (block 2402 ), analyzing the preconditioned data (block 2404 ) and displaying images of the analyzed data (block 2406 ).
- Preconditioning the raw data to produce preconditioned data includes removing a constant frequency component and a system travel time delay (block 2408 ), removing a transmit-antenna to receive-antenna coupling effect (block 2410 ) and prewhitening (block 2412 ) for each spatial location of the raw data.
- Removing the transmit-antenna to receive-antenna coupling effect (block 2410 ) is performed to minimize the effects of ground bounce. In the preferred embodiment, this is accomplished by using a spatial high pass filter which acts in the in track direction. This changes the system from an absolute return system to a relative return system.
- ⁇ mnp ⁇ circumflex over ( ⁇ ) ⁇ mnp ⁇ circumflex over ( ⁇ ) ⁇ m ⁇ p
- the user selects the in track reference scan.
- N 1 and N 2 define a region to be imaged.
- ⁇ mnp ⁇ circumflex over ( ⁇ ) ⁇ mnp ⁇ circumflex over ( ⁇ ) ⁇ ⁇ tilde over (m) ⁇ np
- prewhitening includes applying the following equation:
- analyzing the preconditioned data includes depth focusing (block 2414 ) and synthetic aperture radar (SAR) processing (block 2416 ).
- Depth focusing includes analyzing the preprocessed data to resolve a target in depth.
- the geometry of depth focusing is illustrated in FIG. 25.
- the system measures the transfer function of the ground as a function of frequency. The Fourier transform of the transfer function yields the desired depth response function.
- depth focusing is performed at points, e.g. point 2202 in FIG. 22, directly beneath the points of tangency between the source and receive antennas, as illustrated in FIG. 26.
- SAR processing (block 2416 ), which can enhance the images eventually produced from the raw data, is accomplished by combining data from multiple locations using nearfield, delay and sum beamforming, as illustrated in FIG. 27.
- phase weights in this equation do not depend on absolute sensor position and, in the preferred embodiment, they are precomputed and reused which greatly reduces the time required to perform SAR processing.
- SAR processing (block 2416 ) is an optional procedure which can be selected by the user to enhance the images produced by the system.
- I mnw is the complex image value at spatial location (X F,m ,Y F,n ,Z w );
- U is the SAR array size in the cross-track direction
- V is SAR array size in the along track direction
- (f p 1 , f p 3 ) is the frequency processing band
- ⁇ uvw is the travel time from source (u,v) in the SAR array down to a focal point at depth Z w and back up to receiver (u,v) in the SAR array;
- x F , m 3 ⁇ d 4 + ( m - 1 ) ⁇ d 4 ;
- the preprocessing of the data (block 2402 ) and the analysis of the preprocessed data (block 2404 ) produces a block of data representing a three dimensional volume 2802 as illustrated in FIG. 28.
- the invention allows the display of a plan view of the analyzed data and a side view of the analyzed data.
- Displaying images of the analyzed data includes computing a plan view image of the analyzed data, computing a side view image of the analyzed data (block 2418 ) and displaying the plan view image and the side view image (block 2420 ).
- computing the plan view image of the analyzed data includes applying the following equation:
- PlanView mn max w
- max w is the maximum value across all w (depths).
- computing the side view image of the analyzed data includes applying the following equation:
- max w is the maximum value across all w (depths).
- plan view image and the side view image are computed using other rendering techniques such as averaging or displaying only image values that fall within defined ranges. Any method or technique for presenting a side view and a plan view fall within the scope of the invention.
- FIGS. 29 and 30 show the results of removing the transmit-antenna to receive-antenna coupling effects (block 2410 ) and SAR processing (block 2416 ).
- FIG. 29 shows the data without removing the transmit-antenna to receive-antenna coupling effects (block 2410 ) and SAR processing (block 2416 ).
- FIG. 30 shows the same data after difference referencing using the first data row (block 2410 ) and after SAR processing (block 2416 ).
- a mine is visible at about channel 7 and Y-Range 1.5 in the plan view and at a depth of about 6 and a Y-Range of about 1.5 in the side view.
- FIG. 31 A physical description of a preferred embodiment of the processing performed by the ground penetrating radar is illustrated in FIG. 31.
- the GPSAR program 3102 resides on the computer 108 and provides user control of the ground penetrating radar data acquisition. It accepts a file name, dy (the in track distance between grid points) and array offset (the distance the two antenna arrays are offset from each other).
- the ENCODER program 3104 resides on the computer 108 and triggers the system to take data based on a signal from the movement detector 202 . It also receives the scan data and stores it as raw binary data 3106 .
- the RADARBIN.C program 3108 resides on the computer 108 and converts the raw binary data 3106 to ASCII data sorted by scan 3110 .
- the MGPRVOL.F program 3112 resides on the computer 108 and computes a volumetric image 3114 from the ASCII data 3110 .
- the GPRIMAGE.C program 3116 renders a plan view 3118 and a side view 3120 from the volumetric image 3114 .
- the HTML User Interface 3122 is a web browser such as Netscape Navigator or Microsoft Explorer.
- the programs stored in the memory 908 for the processor 906 in the digital modules control the configuration of their respective components and accept as inputs the start frequency, stop frequency, number of frequency steps in the scan, and the dwell time at each scan step. These programs also provide an Ethernet interface to the computer 108 and serve up a web page that can be accessed from the computer 108 through the HTML user interface 3122 . These programs also collect data from the DSP 1022 .
- the programs stored in the program memory 1026 for the DSP 1022 control the digital module 802 , 804 , the RF module 806 , 808 , and the TX and RX switches 810 , 812 , 814 , and 816 .
Abstract
Description
- [0001] The invention was made with Government Support under Federal Contract Numbers DAAB07-98-C-G014 and DAAB15-00-C-1009 awarded by the United States Army. The Government has certain rights in the invention.
- In general, the invention relates to the field of radar systems. More particularly, the invention relates to the field of ground penetrating radar systems.
- There are numerous applications in which it is useful to view images of underground objects or objects embedded in the ground or in a material such as concrete, that would not otherwise be visible. For example, it is helpful to utility workers to see images of pipes, cables and conduits that are underground where they are about to dig a hole or a trench. It would be even more helpful if the workers could see a three-dimensional image, so that if an imaged object is deep in the ground below the depth they intend to work, they need not be concerned with the object.
- Similarly, a system that would allow workers to see images in three dimensions beneath the surface of the ground would be useful in the field of mine detection. This field has become increasingly important as populations have been moving back into previously war-torn areas of the earth. Frequently, when a war ends, civilians move back into a mined area before the authorities can mount a demining operation. Many civilians are injured or killed by such mines even after a war is over. During a war, demining operations help clear a mine field before foot soldiers or vehicles cross it.
- In such demining operations, it is useful to know accurately not only the surface coordinates of a buried mine, but also the depth at which the mine is buried. It is also important that the equipment used for demining is mobile and that it is capable of exploring a large amount of territory in a short period of time, while maintaining an accurate and detailed record of the objects detected. Such information allows personnel to make intelligent decisions regarding the methods and equipment that are necessary to remove the detected mines.
- The invention is a ground penetrating radar system mounted on a cart to achieve the desired mobility. The system uses two offset banks of interleaved transmit and receive antennas to achieve the desired accuracy. The receive and transmit antennas are properly oriented with respect to each other to reduce cross coupling and maximize desired subsurface echoes. The system uses nearfield beam forming, which is accomplished through fully coherent signal processing and synthetic aperture reception and processing, to image buried objects in three dimensions. The system displays a plan, or top, view and a side view of the area being scanned to provide a three dimensional perspective on a two dimensional computer screen.
- In general, in one aspect, the invention features a ground penetrating radar system which includes a cart configured to be movable along the ground. A computer is mechanically coupled to the cart. A radar electronics module is mechanically coupled to the cart and electrically coupled to the computer. A first antenna array is mechanically coupled to the cart, electrically coupled to the radar electronics module, and oriented to radiate into the ground and receive radiation from the ground. A second antenna array is mechanically coupled to the cart, electrically coupled to the radar electronics module, and oriented to radiate into the ground and receive radiation from the ground. A movement detector, which is configured to detect movement of the cart, is coupled to the computer. The computer is configured to trigger the radar electronics module when the computer detects that the cart has moved a predefined distance.
- Implementations of the invention may include one or more of the following. The radar electronics module may include a first radar electronics module electrically coupled to the first antenna array and a second radar electronics module electrically coupled to the second antenna array. The first antenna array may be configured to radiate and receive radiation from a first series of points along a first set of curves parallel to the direction of movement of the cart. The second antenna array may be configured to radiate and receive radiation from a second series of points along a second set of curves parallel to the direction of movement of the cart. The first set of curves may be interleaved with the second set of curves.
- In general, in another aspect, the invention features a ground penetrating radar system including a first bank of receive antennas arranged along a first axis, a first bank of transmit antennas arranged along a second axis substantially parallel to the first axis and horizontally displaced from the first axis, a second bank of receive antennas arranged along a third axis substantially parallel to the first axis and horizontally displaced from the first axis, and a second bank of transmit antennas arranged along a fourth axis substantially parallel to the first axis and horizontally displaced from the first axis. A first radar electronics module is coupled to the first bank of transmit antennas and the first bank of receive antennas. A second radar electronics module is coupled to the second bank of transmit antennas and the second bank of receive antennas. The transmit antennas in the first bank of transmit antennas are interleaved with the receive antennas in the first bank of receive antennas and the transmit antennas in the second bank of transmit antennas are interleaved with the receive antennas in the second bank of receive antennas. The receive antennas in the first bank of transmit antennas are offset along the first axis from the receive antennas in the second bank of transmit antennas.
- Implementations of the invention may include one or more of the following. The first bank of transmit antennas may be offset along the second axis with respect to the second bank of transmit antennas. The banks of receive antennas may alternate with the banks of transmit antennas. Each transmit antenna may be adjacent to at least one receive antenna. Each transmit antenna may be oriented to minimize electromagnetic coupling to at least one of its adjacent receive antennas. Each transmit antenna may include at least one spiral arm of conductive material. Each receive antenna may include at least one spiral arm of conductive material. A tangent to the inside of the spiral arm at the edge of a transmit antenna may be substantially perpendicular to a tangent to the inside of the spiral arm at the edge of a receive antenna adjacent to the transmit antenna. Each transmit antenna may include two spiral arms of conductive material. Each receive antenna may include two spiral arms of conductive material.
- The transmit antennas and the receive antennas may have faces with centers. Two adjacent first bank receive antennas from the first bank of receive antennas and a first bank transmit antenna from the first bank of transmit antennas interleaved between the two adjacent first bank receive antennas may be positioned such that lines between the centers of the faces of the two adjacent first bank receive antennas and the interleaved first bank transmit antenna form a first triangle having sides of approximately the same length. Two adjacent second bank receive antennas from the second bank of receive antennas and a second bank transmit antenna from the second bank of transmit antennas interleaved between the two adjacent second bank receive antennas may be positioned such that lines between the centers of the faces of the two adjacent second bank receive antennas and the interleaved second bank transmit antenna form a second triangle having sides of approximately the same length. A vertex of the first triangle may be displaced in the direction of the first axis relative to a corresponding vertex of the second triangle by an amount substantially equal to one-half the distance from the center of one side of the first triangle to the center of another side of the first triangle.
- The third axis may be horizontally displaced from the first axis by an amount substantially equal to eight times the distance from the center of one side of the first triangle to the center of another side of the first triangle. The transmit antennas may not be required to be in contact with the ground when in operation. The receive antennas may not be required to be in contact with the ground when in operation.
- In general, in another aspect, the invention features a ground penetrating radar system including a digital module. The digital module includes a direct digital synthesizer configured to generate a digital IF reference signal. A digital to analog converter is coupled to the direct digital synthesizer and is configured to convert the digital IF reference signal to an analog IF transmit signal. An analog to digital converter is configured to convert an analog IF receive signal to a digital IF receive signal. A digital down converter is configured to digitally mix the digital IF receive signal with the digital IF reference signal to produce an in-phase product and the digital IF reference signal shifted in phase by ninety degrees to produce a quadrature product. The ground penetrating radar system includes an RF module coupled to the digital module. The RF module includes an up-converter configured to convert the analog IF transmit signal into a transmit signal and a down-converter configured to convert a receive signal into an analog IF receive signal. The system includes a transmit antenna array coupled to the up-converter for radiating the transmit signal and a receive antenna array coupled to the down-converter for receiving the receive signal.
- Implementations of the invention may include one or more of the following. The transmit antenna array may include a plurality of transmit antennas. The receive antenna array may include a plurality of receive antennas. The system may include a digital signal processor. The system may include a transmit switch for applying the transmit signal to one of the plurality of transmit antennas. The transmit switch may be controlled by the digital signal processor. The system may include a receiver switch for receiving the receive signal from one of the plurality of receive antennas. The receiver switch may be controlled by the digital signal processor. The digital signal processor may control the direct digital synthesizer, the digital down converter, the up-converter and the down-converter. The transmit signal may be a stepped-frequency transmit signal. The receive signal may be a stepped-frequency receive signal. The system may include a computer coupled to a processor through an extensible network. The processor may be configured to communicate with the digital signal processor. The extensible network may be an Ethernet network.
- In general, in another aspect, the invention features a ground penetrating radar system including a digital module configured to generate an analog IF transmit signal and to receive an analog IF receive signal. The system includes an RF module, which includes a triple-heterodyne up-converter for converting an analog IF transmit signal into a stepped-frequency transmit signal. The RF module also includes a triple-heterodyne frequency converter for converting a stepped-frequency receive signal into an analog IF receive signal. The system includes a transmit antenna coupled to the up-converter for radiating the stepped-frequency transmit signal and a receive antenna coupled to the down-converter for receiving the stepped-frequency receive signal.
- Implementations of the invention may include one or more of the following. The triple-heterodyne up-converter may include a first up-converter configured to mix the analog IF transmit signal with a signal from a first local oscillator to produce a first intermediate signal and an aliased first intermediate signal. The triple-heterodyne up-converter may include a first filter coupled to the first up-converter for substantially rejecting the aliased first intermediate signal. The triple-heterodyne up-converter may include a second up-converter coupled to the first filter configured to mix the first intermediate signal with a signal from a second local oscillator to produce a second intermediate signal and an aliased second intermediate signal. The triple-heterodyne up-converter may include a second filter coupled to the second up-converter for substantially rejecting the aliased second intermediate signal. The triple-heterodyne up-converter may include a down-converter coupled to the second filter configured to mix the second intermediate signal with a stepped frequency signal to produce the stepped-frequency transmit signal and an aliased stepped-frequency transmit signal. The stepped-frequency transmit signal may have substantially no frequency components in the pass bands of the first filter or the second filter. The triple-heterodyne up-converter may include a third filter coupled to the down-converter for substantially rejecting the aliased stepped-frequency transmit signal.
- The triple-heterodyne up converter may include an up-converter configured to mix the stepped-frequency receive signal with a stepped-frequency signal to produce a first intermediate signal and an aliased first intermediate signal. The triple-heterodyne up-converter may include a first filter coupled to the first up-converter for substantially rejecting the aliased first intermediate signal. The triple-heterodyne up-converter may include a first down-converter coupled to the first filter configured to mix the first intermediate signal with a signal from a first local oscillator to produce a second intermediate signal and an aliased second intermediate signal. The triple-heterodyne up-converter may include a second filter coupled to the first down-converter for substantially rejecting the aliased second intermediate signal. The triple-heterodyne up-converter may include a second down-converter coupled to the second filter configured to mix the second intermediate signal with a second local oscillator to produce the analog IF receive signal and an aliased analog IF receive signal. The triple-heterodyne up-converter may include a third filter coupled to the second down-converter for substantially rejecting the aliased analog IF receive signal.
- In general, in another aspect, the invention features a ground penetrating radar system including a transmitter, a receiver, an array of transmit antennas, an array of receive antennas interleaved with the array of transmit antennas, a transmit switch configured to selectively couple the transmitter to one of the array of transmit antennas and a receive switch configured to selectively couple the receiver to one of the array of receive antennas. The array of transmit antennas is arranged in one or more rows. The array of receive antennas is arranged in one or more rows. Each row is parallel to, adjacent to and offset from one of the rows of transmit antennas, so that each receive antenna in a row except one is adjacent to two transmit antennas, and each transmit antenna in a row except one is adjacent to two receive antennas. The transmit switch and the receive switch are configured to couple the transmitter and receiver, respectively, to a first transmit antenna and a first adjacent receive antenna, and subsequently to the first transmit antenna and a second adjacent receive antenna.
- In general, in another aspect, the invention features a method for collecting and displaying data from a ground penetrating radar system, which includes a plurality of transmit antennas and a plurality of receive antennas. Each transmit antenna, except one, has two adjacent receive antennas. The system is mounted on a movable cart. The method includes collecting raw data. Collecting raw data includes (a) selecting a first of the plurality of transmit antennas. Collecting raw data further includes (b) selecting a first receive antenna that is adjacent to the selected transmit antenna. Collecting raw data further includes (c) collecting data using the selected transmit antenna and the selected receive antenna to produce raw data. The raw data collected at spatial location (xm, vn) is denoted by {tilde over (Ψ)}mnp where the indices m, n are used to denote position in a grid of spatial locations where data has been collected, and p is an index ranging from 1 to P corresponding to the frequency fp at which the data was collected. Collecting raw data includes (d) repeating step (c) for both receive antennas adjacent to the selected transmit antenna. Collecting raw data further includes (e) repeating steps b, c and d for all transmit antennas. Collecting raw data further includes repeating steps a, b, c, d, and e each time the cart moves to a new location. The method further includes preconditioning the raw data to produce preconditioned data, analyzing the preconditioned data, and displaying images of the analyzed data.
- Implementations of the invention may include one or more of the following. Preconditioning the raw data to produce preconditioned data may include (g) removing a constant frequency component and a system travel time delay, (h) removing a transmit-antenna to receive-antenna coupling effect, (i) prewhitening, and (h) repeating steps (g), (h) and (i) for each spatial location of the raw data.
-
- Removing the transmit-antenna to receive-antenna coupling effect may include applying the following equation:
- Ψmnp={circumflex over (Ψ)}mnp−{circumflex over (Ψ)}mñp
- where Ψmñp is an in track reference scan.
-
- where N1 and N2 define a region to be imaged.
-
- where aq are digital filter coefficients chosen to reject low frequency spatial energy.
- Removing the transmit-antenna to receive-antenna coupling effect may include applying the following equation:
- Ψmnp={circumflex over (Ψ)}mnp−{circumflex over (Ψ)}{tilde over (m)}np
- where {circumflex over (Ψ)}{tilde over (m)}np is a cross line reference scan.
- Prewhitening may include applying the following equation:
- γmnp =b p·Ψmnp
- where bp are frequency dependent weights.
-
- where
- Imnw the complex image value at spatial location (XF,m, YF,n, Zw);
- U is the SAR array size in the cross-track direction;
- V is SAR array size in the along track direction;
- (fp
1 , fp2 ) is the frequency processing band; -
- YF,m=0.933013d+(n−1)dy;
- d=5.52 inches; and
- dy=scan spacing.
-
- where (xs,u,ys,v) and (xr,u,yr,v) are the location of the transmit and receive antennas and cg is the speed of light in the ground.
- Displaying images of the analyzed data may include computing a plan view image of the analyzed data, computing a side view image of the analyzed data, and displaying the plan view image and the side view image.
- Computing a plan view image of the analyzed data may include applying the following equation:
- PlanViewmin=maxw |I mnw|2
- where
- maxw is the maximum value across all w (depths).
- Computing a side view image of the analyzed data may include applying the following equation:
- SideViewmw=maxm |I mnw|2
- where
- maxw is the maximum value across all w (depths).
- FIG. 1 is a perspective view of a ground penetrating radar system according to the present invention.
- FIG. 2 is a block diagram of a ground penetrating radar system.
- FIG. 3 illustrates the grid of radar data collected by a ground penetrating radar system.
- FIG. 4 illustrates the geometry of antenna arrays incorporated in the ground penetrating radar system.
- FIG. 5 illustrates the switching among the antennas within the antenna arrays.
- FIG. 6 illustrates the configuration of the antennas.
- FIG. 7 illustrates the relative orientation of a transmit antenna and a receive antenna.
- FIG. 8 is a block diagram of the electronics in a ground penetrating radar system.
- FIG. 9 is a block diagram of a digital module.
- FIG. 10 is a block diagram of an embodiment of the signal synthesis and analysis module within a digital module.
- FIG. 11 is a block diagram of a digital down-converter.
- FIG. 12A is a block diagram of an embodiment of the signal synthesis and analysis module within a digital module.
- FIG. 12B is a block diagram of a field programmable gate array.
- FIG. 13 is a block diagram of the RF module.
- FIG. 14 is a frequency plan for a preferred embodiment of the ground penetrating radar system.
- FIG. 15 illustrates a stepped-frequency signal.
- FIG. 16 is a block diagram of the first stage in a triple-heterodyne frequency up-converter and the last stage in a triple-heterodyne frequency converter.
- FIG. 17 is a block diagram of a middle stage in a triple-heterodyne frequency up-converter and the middle stage in a triple-heterodyne down-converter.
- FIG. 18 illustrates the last stage in a triple-heterodyne frequency up-converter and the first stage in a triple-heterodyne frequency converter.
- FIG. 19 is a block diagram of amplifiers and filters in the triple-heterodyne up-converter and the triple-heterodyne frequency converter.
- FIG. 20 is a block diagram of a local oscillator.
- FIG. 21 is a block diagram of the transmit switch and the receive switch.
- FIG. 22 illustrates the grid of points collected by movement of the ground penetrating radar system along the ground.
- FIG. 23 is a block diagram of the raw data collection processing.
- FIG. 24 is a block diagram of the raw data analysis processing.
- FIG. 25 illustrates the geometry of depth focusing.
- FIG. 26 illustrates the geometry of collection of points midway between a transmit antenna and a receive antenna.
- FIG. 27 illustrates the geometry of synthetic aperture radar processing.
- FIG. 28 illustrates a volume of collected data and a top view and a side view of that data.
- FIG. 29 illustrates images of data prior to processing.
- FIG. 30 illustrates images of data after processing.
- FIG. 31 illustrates a physical block diagram of the processing performed in the ground penetrating radar system.
- In the following detailed description and in the several figures of the drawings, like elements are identified with like reference numerals.
- As shown in the drawings for purposes of illustration, the invention is embodied in a novel ground penetrating radar system. A system according to the invention includes interleaved antenna arrays of properly oriented spiral transmit and receive antennas. Switching among the antennas allows samples to be efficiently taken with high spatial density. A digital frequency synthesizer and digital down converter allow fast and accurate measurements. Triple heterodyne up and down conversion reduces the likelihood that mixing products will interfere with the measurements. The system is extensible because of an Ethernet interconnection between system components. The system uses advanced processing and spatial filtering techniques to improve the quality of the images produced.
- A ground penetrating
radar system 102, illustrated in FIG. 1, includes acart 104 that can be moved along the ground. In the preferred embodiment, thecart 104 has four allterrain wheels 106 that give thecart 104 mobility over all types of terrain. Other techniques for providing mobility, such as tracks or skids, are within the scope of the invention. - In an preferred embodiment, the
cart 104 is a lightweight, compact design fabricated out of aluminum tubing. Four inch square, thin wall tubing is used for the frame. All joints are welded and the frame is painted with polyurethane paint. Mount points for all major system components are built into the frame with easy access to each mount for individual component removal without tools. In other embodiments, the cart is fabricated out of other metals, such as steel, or other materials, such as light-weight composites. - The
system 102 includes acomputer 108 mechanically coupled to thecart 104. In the preferred embodiment illustrated in FIG. 1, thecomputer 108 is a laptop computer, having its own display, memory, long term storage and other peripherals, that rests on or is secured to thecart 104. In other embodiments, the computer is a desktop-type computer. - A
radar electronics module cart 104 and electrically coupled to thecomputer 108. In a preferred embodiment, theradar electronics module anodized aluminum enclosures - A first (or front)
antenna array 114 is mechanically coupled to thecart 104 and electrically coupled to theradar electronics module first antenna array 114 is oriented to radiate into the ground and receive radiation from the ground. A second (or back)antenna array 116 is mechanically coupled to thecart 104 and electrically coupled to theradar electronics module second antenna array 116 is oriented to radiate into the ground and receive radiation from the ground. - The two
antenna arrays cart 104 by anantenna framework 118. Theantenna arrays antenna framework 118 can be removed from thecart 104 for diagnostics and transport. The height of theantenna arrays machine screw jack 120 to be, in a preferred embodiment, between six inches and twenty inches. Adrill motor 122 drives themachine screw jack 120 to raise and lower the antenna arrays. Thedrill motor 122 could be replaced by any type of motor that has sufficient power to raise and lower the antenna arrays. Theantenna framework 118 offsets the antenna banks (in a preferred embodiment by 24 inches) from thecart 104 frame to eliminate any possibility of target shadowing by the frame. - A
movement detector 202, shown in FIG. 2, coupled to the computer detects movement of thecart 104. In the preferred embodiment, themovement detector 202 is an optical encoder mounted on the front right wheel of thecart 104. - The
computer 108 triggers theradar electronics module computer 108 detects that thecart 104 has moved a predefined distance. In this way, data can be acquired at prescribed intervals measured by themotion detector 202. - In a preferred embodiment, the
radar electronics module radar electronics module 110 is coupled to thefirst antenna array 114. A secondradar electronics module 112 is coupled to thesecond antenna array 116. - The
first antenna array 114 radiates and receives radiation from a first series of points, e.g. 302 shown in FIG. 3, along a first set ofcurves 304 parallel to the direction ofmovement 306 of thecart 104. In FIG. 3, the open circles represent the locations of the centers of the transmit and receive antennas when they take a radar sample, as described below. The closed, or filled, circles represent the location of the samples. In a preferred embodiment, the curves in the first set ofcurves 304 are approximately d/2, or 2.76 inches, apart. - Similarly, the
second antenna array 116 radiates and receives radiation from a second series of points, e.g. 308, along a second set ofcurves 310 parallel to the direction ofmovement 306 of thecart 104. Note that in FIG. 3, thesecond antenna array 116 is shown a large distance behind thefirst antenna array 114. Actually, the two arrays are preferably arranged in much closer proximity to each other, as shown in FIGS. 1 and 2 and as discussed below in the discussion of FIG. 4. In a preferred embodiment, the curves in the second set ofcurves 310 are approximately d/2, or 2.76 inches, apart. The first set ofcurves 304 is interleaved with the second set ofcurves 310. The result is a set of curves including the first set ofcurves 304 and the second set ofcurves 310, with the curves being approximately d/4, or 1.38, inches apart. - The parameter d in FIG. 3 denotes the outside diameter of an antenna. In a preferred embodiment, d-5.515 inches. The individual transmit and receive antennas are of the log spiral type and were designed to operate over the band 800-4000 MHz but actually radiate well down to 500 MHz. The antenna cavity is nine inches tall and is filled with radar absorbing material to minimize reflections from the top of the antenna. The transmit and receive antennas have different windings in order to minimize cross coupling during the simultaneous transmission and reception required by a stepped-frequency radar.
- The
first antenna array 114 includes a first bank of receiveantennas 402 arranged along afirst axis 404 and a first bank of transmitantennas 406 arranged along asecond axis 408 substantially parallel to the first axis and horizontally displaced from the first axis, as shown in FIG. 4. Thesecond antenna array 116 includes a second bank of receiveantennas 410 arranged along athird axis 412 substantially parallel to thefirst axis 404 and horizontally displaced from thefirst axis 404 and a second bank of transmitantennas 414 arranged along afourth axis 416 substantially parallel to thefirst axis 404 and horizontally displaced from thefirst axis 404. - In a preferred embodiment, the first
radar electronics module 110 is coupled to the first bank of receiveantennas 402 and the first bank of transmitantennas 406. The secondradar electronics module 112 is coupled to the second bank of receiveantennas 410 and the second bank of transmitantennas 414. - As can be seen in FIG. 4, the receive antennas in the first bank of receive
antennas 402 are interleaved with the transmit antennas in the first bank of transmitantennas 406. Similarly, the receive antennas in the second bank of receiveantennas 410 are interleaved with the transmit antennas in the second bank of transmitantennas 414. - Further, the receive antennas in the first bank of receive
antennas 402 are offset along thefirst axis 404 from the receive antennas in the second bank of receive antennas. Similarly, the first bank of transmitantennas 406 is offset with respect to the second bank of transmitantennas 414. The banks of transmitantennas antennas - Each receive antenna, e.g.418, is adjacent to at least one transmit antenna, e.g. 420, 422. Each receive antenna, e.g. 418, is oriented to minimize electromagnetic coupling with at least one of its adjacent transmit antennas, e.g. 420, 422.
- The close spacing of the
curves cart movement direction 306 and taking samples at regular distance intervals (measured by the movement detector 202), an array of data may be acquired. Further, samples at a number of frequencies are taken at each sample location. - The data structure produced by this process is a complex array of numbers of the form y(m,n,p) where m is the cross-track index (m=1,2, . . . ,26), n is the row or scan index and p is the frequency index. The quantity y(m,n,p) is the complex frequency response of the ground at the location (m,n). It is the quantity measured by the stepped-frequency radar. The Fourier transform of y(m,n,p) is the equivalent time domain response of the radar at location (m,n).
- Each transmit
antenna 602, illustrated in FIG. 6, includes at least onespiral arm 604 of conductive material. Similarly, each receive antenna, e.g. 606, includes at least onespiral arm 608 of conductive material. - The transmit
antenna 602 and one of its adjacent receiveantennas 606 are oriented so that a tangent 702 to the inside of thespiral arm 604 at the edge of a transmitantenna 602 is substantially perpendicular to a tangent 704 to the inside of thespiral arm 608 at the edge of a receiveantenna 606 adjacent to the transmit antenna. This orientation minimizes electromagnetic cross coupling between the transmitantenna 602 and the receiveantenna 606. - In a preferred embodiment, each transmit antenna, e.g.602, includes two
spiral arms antenna 606 includes twospiral arms - As shown in FIG. 4, the transmit antennas, e.g.424 and 426, and the receive antennas, e.g. 428, have faces. Each of the faces has a center. Two adjacent first bank receive
antennas antennas 402 and a first bank transmitantenna 428 from the first bank of transmitantennas 406 interleaved between the two adjacent first bank receiveantennas antennas antenna 428 form afirst triangle 430 having sides of approximately the same length. - Two adjacent second bank receive
antennas antennas 410 and a second bank transmitantenna 436 from the second bank of transmitantennas 414 interleaved between the two adjacent second bank receiveantennas antennas antenna 436 form asecond triangle 438 having sides of approximately the same length. - A
vertex 440 of thefirst triangle 430 is displaced in the direction of thefirst axis 404 relative to acorresponding vertex 442 of thesecond triangle 438 by an amount substantially equal to one-half the distance from the center of one side of thefirst triangle 432 to the center of another side of thefirst triangle 432. - In a preferred embodiment, the
third axis 412 is horizontally displaced from thefirst axis 404 by an amount substantially equal to eight times the distance from the center of one side of thefirst triangle 432 to the center of another side of thefirst triangle 432. - In a preferred embodiment, neither the transmit antennas, e.g.420, nor the receive antennas, e.g. 422, are required to be in contact with the ground when in operation.
- In a preferred embodiment, the ground penetrating radar system operates in a continuous wave mode. The first and second
radar electronics modules - The first and second
radar electronics modules digital module computer 108, and anRF module digital module digital module RF module digital module switch 810. The transmitswitch second antenna array - A stepped frequency return signal is received by the
antenna array receiver switch receiver switch RF module RF module digital module digital module computer 108 for storage and processing. - The
digital module extensible network interface 902, which in the preferred embodiment is a Ethernet interface, and a point-to-point communication interface 904, which in the preferred embodiment is an RS-232 interface. Both theextensible network interface 902 and the point-to-point communication interface 904 allow communication between thecomputer 108 and aprocessor 906 within thedigital module computer 108 to thedigital module digital module computer 108. Data and program code for theprocessor 906 are stored within the digital module on amemory 910. A signal synthesis andanalysis module 908 synthesizes the transmit IF signal, controls thetransmitter switch RF module receiver switch - The signal synthesis and
analysis module 908, shown in more detail in FIG. 10, includes a directdigital synthesizer 1002, which generates a digital IFreference signal 1004. A digital toanalog converter 1006 converts the digital IF signal to an analog IF transmitsignal 1008. In a preferred embodiment, the analog IF transmitsignal 1008 is a 10.7 MHz IF signal. In a preferred embodiment, the directdigital synthesizer 1002 and the digital toanalog converter 1006 are incorporated into asingle module 1010, such as the AD7008 manufactured by Analog Devices. The AD7008 includes a 10-bit analog to digital converter. - The digital IF
reference signal 1004 provides a stable and accurate reference for both the transmit and receive paths of the apparatus. As will be seen, the digital IFreference signal 1004 is used in the down conversion process without the necessity of applying delays. - The signal synthesis and
analysis module 908 also includes an analog todigital converter 1012 which digitizes the analog IF receivesignal 1014 to produce a digital IF receivesignal 1016. The digital IF receivesignal 1016 is routed to adigital down converter 1018 which digitally mixes the digital IF receivesignal 1016 with the digital IFreference signal 1004 from the directdigital synthesizer 1002 through a programmable logic device (PLD) 1020. ThePLD 1020 provides control logic and interfacing for the directdigital synthesizer 1002 and thedigital down converter 1018 and a digital signal processor (DSP) 1022. TheDSP 1022 controls the directdigital synthesizer 1002, thedigital down converter 1018, and the other peripherals in the system. The PLD also receives an Out of Range signal from the analog todigital converter 1012, which it communicates to theDSP 1022. TheDSP 1022 reacts by changing the gain in the receive chain, as discussed below. - The
digital down converter 1018 digitally mixes 1102 the digital IFreference signal 1004 with the digital IF receivesignal 1016 to produce an in-phase product 1104, as shown in FIG. 11. The in-phase product 1104 is then digitally filtered by alow pass filter 1106 to produce a filtered inphase component 1108 which is transferred to theDSP 1022 for processing and storage. - The
digital down converter 1018 also digitally shifts the phase of the digital IF reference signal by ninetydegrees 1110 and digitally mixes 1112 it with the digital IF receivesignal 1016 to produce aquadrature product 1114. Thequadrature product 1114 is then digitally filtered through alow pass filter 1116 to produce a filteredquadrature component 1118 which is passed to theDSP 1022 for processing and storage. - In a preferred embodiment, the
low pass filters - Returning to FIG. 10, the signal synthesis and
analysis module 908 includes amemory 1024 where theDSP 1022 and thePLD 1020 store data. The signal synthesis andanalysis module 908 also includes aprogram memory 1026 where program code for theDSP 1022 is stored. - Another example embodiment of the
digital module 908 is illustrated in FIG. 12A. TheDSP 1022, thememory 1024 and theprogram memory 1026 remain in the same configuration as shown in FIG. 10. The data and programs stored in thememory 1024 and theprogram memory 1026 may be different from the data and programs stored in the same elements in FIG. 10. APLD 1202 provides a control logic interface between theDSP 1022 and a field programmable gate array (FPGA) 1204. Amemory 1206, preferably a programmable read only memory (PROM), stores the FPGA boot program. Aclock 1208, preferably an 85.6 MHz clock, provides a clock signal to the FPGA, a digital to analog converter (DAC) 1210 and an analog to digital converter (ADC) 1212. TheFPGA 1204 is controlled by thePLD 1202, which in turn is controlled by theDSP 1022, throughcontrol lines 1214. - The
FPGA 1204 generates a digital IFreference signal 1216. TheDAC 1210 converts the digital IFreference signal 1216 to an analog signal which is isolated by transformer 1218 and provided as the analog IF transmitsignal 1008. - The analog IF receive
signal 1014 is converted by theADC 1212 to a digital receivesignal 1220 which is provided to theFPGA 1204. TheADC 1212 also produces an out ofrange signal 1222 if the analog IF receivesignal 1014 exceeds the range of theADC 1212. - The
FPGA 1204 produces receiveddata 1224 which is transferred to theDSP 1022 for processing. - In one example embodiment of the FPGA, illustrated in FIG. 12B, the digital receive
signal 1220 is provided to twodigital mixers digital mixer 1226 mixes the digital receivesignal 1220 with adigital cosine wave 1230 having a frequency, in one preferred embodiment, of 10.7 MHZ. The output of the firstdigital mixer 1226 is thereal component 1232 of the digital receivesignal 1220. Thereal component 1232 is filtered 1234, preferably by summing a 128 sample moving window of data, to produce a filteredreal component signal 1236 which is provided to amultiplexer 1238. - The second
digital mixer 1228 mixes the digital receivesignal 1220 with adigital sine wave 1240 having a frequency, in one preferred embodiment, of 10.7 MHZ. Preferably, thedigital sine wave 1240 is identical to thedigital cosine wave 1230 except that thedigital sine wave 1240 is shifted ninety degrees in phase with respect to thedigital cosine wave 1230. The output of the firstdigital mixer 1228 is theimaginary component 1242 of the digital receivesignal 1220. Theimaginary component 1242 is filtered 1244, preferably by summing a 128 sample moving window of data, to produce a filteredimaginary component signal 1246 which is provided to themultiplexer 1238. - The
digital sine wave 1240 is also provided as the digital IFreference signal 1216. - The
multiplexer 1238 provides to theDSP 1022 either the filteredreal component signal 1236 or the filteredimaginary component signal 1246, depending on theselect signal 1248. - The
RF module converter 1302 for converting the analog IF transmitsignal 1008 into a stepped-frequency transmitsignal 1304. TheRF module heterodyne frequency converter 1306 for converting a stepped frequency receivesignal 1308 into the analog IF receivesignal 1014. - The triple-heterodyne up-
converter 1302 includes a first up-converter 1310. The first up-converter 1310 includes amixer 1602, as shown in FIG. 16, which mixes the analog IF transmit signal 1008 (after it has been isolated by transformer 1604 and filtered by low pass filter 1606) with the signal from a firstlocal oscillator 1312 to produce a firstintermediate signal 1314 and an aliased firstintermediate signal 1402, as shown in FIG. 14. These two signals are isolated by transformer 1608 and amplified byamplifiers filter 1614, having apass band 1404, shown in FIG. 14, substantially rejects the aliased firstintermediate signal 1402. In the preferred embodiment, the firstlocal oscillator 1312 operates at a frequency of 122 MHZ and when mixed with the 10.7 MHZ IF produces mixing products at 111.3 MHZ and 132.7 MHZ. In the preferred embodiment, thefilter 1614 has a center frequency of 139.75 MHZ and passes the 132.7 MHZ mixing product. - A second up-
converter 1316 includes a mixer 1702, as shown in FIG. 17, which mixes the firstintermediate signal 1314 with the signal produced by the secondlocal oscillator 1318 to produce a secondintermediate signal 1320 and an aliased secondintermediate signal 1406, as shown in FIG. 14. Afilter 1704 substantially rejects the aliased secondintermediate signal 1406. In the preferred embodiment, the secondlocal oscillator 1318 operates at a frequency of 2275 MHZ and when mixed with the firstintermediate signal 1314 produces mixing products at 2142.3 MHZ and 2407.7 MHZ. In the preferred embodiment, thefilter 1704 has a center frequency of 2450 MHZ and passes the 2407.7 MHZ product. - A down-
converter 1322 includes amixer 1802, shown in FIG. 18, that mixes the second intermediate signal 1320 (after being amplified by amplifier 1804) with a steppedfrequency signal 1324 generated by asynthesizer 1326 to produce the stepped-frequency transmitsignal 1328 and an aliased stepped-frequency transmit signal (not shown in FIG. 14). The stepped-frequency transmit signal has substantially no frequency components in the pass bands of thefirst filter 1614 or thesecond filter 1704. Athird filter 1902, shown in FIG. 19, substantially rejects the aliased stepped-frequency transmit signal (after it is amplified by amplifier 1904).Amplifier 1906 amplifies the output of thethird filter 1902 to produce the stepped-frequency transmitsignal 1304. In the preferred embodiment, thethird filter 1902 is a low pass filter with a 3 dB break point at 2000 MHZ. - The
synthesizer 1326 is controlled by theDSP 1022 through the LO Control and Configuration signals shown on FIG. 13 to produce a stepped-frequency signal 1324, illustrated in FIG. 15. The frequency of the stepped-frequency signal 1324 steps through a range offrequencies 1502 under the control of theDSP 1022. In the preferred embodiment the range offrequencies 1502 is from 500 to 2000 MHZ. In the preferred embodiment, when the stepped-frequency signal 1324 is mixed with the secondintermediate signal 1322 the result is a preamplified stepped-frequency transmitsignal 1328 ranging from 407.7 MHZ to 1907.7 MHZ. - The triple-heterodyne up-
converter 1302 includes amplifiers andfilters 1330 which amplify and filter the preamplified stepped-frequency transmitsignal 1328 before it is sent to the transmitswitch antenna array signal 1304. - The triple-heterodyne up
converter 1306 includes an amplifier/filter 1332 that amplifies, using amplifier 1908, and filters, usingfilter 1910, the stepped-frequency receivesignal 1308 to produce an amplified stepped-frequency receivedsignal 1334, as illustrated in FIG. 19. - An up-
converter 1336 usesmixer 1806 to mix the stepped-frequency receive signal 1334 (after amplification by amplifier 1808) with the stepped-frequency signal 1324 output of thesynthesizer 1326 to produce a firstintermediate signal 1338 and an aliased firstintermediate signal 1410, as shown in FIG. 14. Afirst filter 1810 substantially rejects the aliased first intermediate signal 1410 (after it has been amplified by amplifier 1812). - A first down-
converter 1340 uses amixer 1706 to mix the firstintermediate signal 1338 with the signal produced by the secondlocal oscillator 1318 to produce a secondintermediate signal 1342 and an aliased secondintermediate signal 1412. Asecond filter 1708 substantially rejects the aliased second intermediate signal 1412 (after it is amplified by amplifier 1710). - A second down-
converter 1344 uses amixer 1616 to mix the second intermediate signal 1342 (after it is amplified by amplifier 1618) with the signal produced by the firstlocal oscillator 1312 to produce the analog IF receivesignal 1014 and an aliased analog IF receive signal (not shown). Athird filter 1620 substantially rejects the aliased analog IF receive signal. An variable-gain amplifier 1622 between themixer 1616 and thethird filter 1620 allows the DSP to control the receive gain using gain control signals 1624. Atransformer 1626 isolates themixer 1616 from the amplifier 1624. Anamplifier 1628 amplifies the output of thethird filter 1620 and atransformer 1630 isolates the output of theamplifier 1628 from theoutput 1014. - In a preferred embodiment, the second local oscillator, shown in detail in FIG. 17, is a phase locked loop including a phase locked loop (PLL)
frequency synthesizer 1712. In a preferred embodiment,PLL frequency synthesizer 1712 is a MC145202 manufactured by Motorola. - The
PLL frequency synthesizer 1712 has an SDI input, a SCLK input and an SLD input from theDSP 1022. The SDI input is the serial data line from theDSP 1022 that provides the configuration data for thePLL frequency synthesizer 1712. Through the configuration data transferred through the SDI input, theDSP 1022 can specify the frequency at which the phase lockedloop 1318 is to operate. In the preferred embodiment, the SDI signal is also used to transfer data into thelocal oscillator 1326, as discussed below. - The SCLK is used to clock the transfer of serial data into the
PLL frequency synthesizer 1712. In the preferred embodiment, the SCLK signal is also used to clock the transfer of data into the LO, as described below. The SLD signal is activated when the SDI data is intended for thePLL frequency synthesizer 1712. - The
PLL frequency synthesizer 1712 provides aloop error signal 1716 as an input to the voltage controlled oscillator (“VCO”) 1714. The frequency of theoutput 1718 of theVCO 1714 varies with theloop error signal 1716. The VCO output is amplified byamplifier 1720 and provided as the output oflocal oscillator 1318 throughamplifiers amplifier 1720 is frequency divided bydivider 1726 to produce an input to thePLL frequency synthesizer 1712, thereby closing the phase locked loop. - In other embodiments, the phase locked
loop 1318 is replaced by other apparatus for generating a stable frequency signal, such as a crystal. - The
local oscillator 1326, shown in more detail in FIG. 18, includes anLO synthesizer 1814, which generates the steppedfrequency signal 1324.Amplifiers frequency signal 1324 for use by thedown converter 1322 and the upconverter 1336. - The
LO synthesizer 1814 is controlled by the DSP usingcontrol lines 1820 and produces acalibration signal PN 1822, which is used by the DSP to calibrate the LO. - The
LO synthesizer 1814, illustrated in detail in FIG. 20, includes a digital to analog converter (DAC) 2002 which is provided avoltage reference 2004. The level of the analog output of theDAC 2002 is determined by the SDI, SCLK and SLD inputs from the DSP. These signals were described above in the description of the phase lockedloop 1318. When the SLD signal enables theDAC 2002, the SCLK signal clocks data in the SDI signal into theDAC 2002. The data clocked into the DAC determines the level of its output. - The DAC output is amplified by an
amplifier 2005 and used to drive a voltage controlled oscillator (VCO) 2006. Thus, by selecting the data to store in theDAC 2002, theDSP 1022 can control the frequency of the output of theVCO 2006. In a preferred embodiment, theDSP 1022 controls the frequency of the output of theVCO 2006 to vary in steps, as shown in FIG. 15, from 1000 to 2000 MHZ. Anamplifier 2008 amplifies the output of theVCO 2006. - The system generates the entire desired range of frequencies, which, in a preferred embodiment is from 500 to 2000 MHZ, by using a series of switches. To generate the stepped frequency signal over the range from 500 to 1000 MHZ, the DSP unasserts the SHI signal which causes an
input switch 2010 to switch to its “lo” position and anoutput switch 2012 to switch to its “lo” position. This configuration of switches causes the amplified output of theVCO 2006 to be routed to afrequency divider 2104, which divides the frequency of the signal by two. The result is filtered bylow pass filter 2016 and provided as the steppedfrequency signal 1324. - The
DSP 1022 can generate the higher range of frequencies (e.g., from 1000 to 2000 MHZ) by asserting the SHI signal, which causes theinput switch 2010 and theoutput switch 2012 to switch to their respective “HI” positions. Similarly,switch 2018 is also closed allowing the amplified output of theVCO 2006 to be filtered by a low pass filter and provided as the steppedfrequency signal 1324. - The output of the
VCO 2006 is also used to generate a calibration signal PN for use by theDSP 1022. To accomplish this, the output of theVCO 2006 is amplified byamplifier 2022 and its frequency is divided by 64 by using a series ofdividers comparator 2030 squares the edges of the divided signal. TwoPALs DSP 1022. - The transmit
switch throw switches 2102 with their common poles wired in parallel and connected to the stepped-frequency transmitsignal 1304. The other pole of each switch is connected to one of the transmit antennas. Therefore, when one of theswitches 2102 is actuated, the stepped frequency transmitsignal 1304 is connected to a respective transmit antenna. - The selection of which of the
switches 2102 to close is controlled by theDSP 1022 through the transmitter switch control lines, which include control lines TXCTL1, TXCTL2, and TXCTL3. The binary combination of these control lines determine which of the transmitter switches 2102 will close. - The
receiver switch throw switches 2104 with their common poles wired in parallel and connected to the stepped frequency receivedsignal 1308. The other pole of each switch is connected to one of the received antennas in the receive array. Therefore, when one of theswitches 2104 is closed, the corresponding receive antenna provides the stepped frequency receivedsignal 1308. - The determination of which of the
switches 2104 is to close is controlled by theDSP 1022 through the receive switch control signals, which includes three switch control signals RXCTL1, RXCTL2, and RXCTL3. Theswitch 2104 that is closed depends on the binary combination of those three control signals. - As discussed above, in a preferred embodiment, the system records data at 13 locations spaced 2.76 inches apart for each
antenna array computer 108 fuses this data into a grid of rows and columns of data, where each row includes data taken from 26 locations spaced 1.38 inches apart. The columns represent the data taken at different locations along the direction of cart movement 306 (see FIG. 3). - In general, the system records stepped-frequency radar data over an n row by m channel grid, as illustrated in FIG. 22, with the location of the individual points in the grid being expressed as (Xm, yn) where xm=m·dx, yn=n·dy, where dx is the space between points in the cross track direction (see FIG. 3), m is an index in the cross track direction, dy is the sampling interval in the in track direction, and n is an index in the in track direction. In the preferred embodiment, dx is 1.38 inches and dy is operator selectable, and is preferably selected to be 1.38 inches, producing a rectangular array. In the preferred embodiment, m varies from 1 to 26 and n varies from 1 to N, where N is the total number of scans that the system performs.
- In the preferred embodiment, data channels m=1,3,5, . . . , 25 are recorded using the
first antenna array 114 and data channels m=2,4,6, . . . , 26 are recorded using thesecond antenna array 116. The cross track resolution when both antenna arrays are used is 1.38″. Individual images can be formed using either thefirst antenna array 114 alone or the second antenna array alone. For individual images of these two types, the cross track resolution is 2.76″. - The processing that the system performs includes collecting raw data and analyzing the raw data. Collecting the raw data, illustrated in FIG. 23, includes detecting movement of the cart to the next data collection position (block2302). The system accomplishes this by monitoring the
movement detector 202. The system then selects one of the transmit antennas (block 2304) and determines if all of the transmit antennas have been processed (block 2306). If they have, the process terminates (block 2308) until the cart moves to the next data collection position (block 2302). - If the all of the transmit antennas have not yet been processed, the system selects a receive antenna adjacent to the selected transmit antenna (block2310). The system then determines if both of the receive antennas adjacent to the selected transmit antenna have been processed (block 2312). If they have, the system selects the next transmit antenna (block 2314) and returns to the beginning of the loop (to block 2306).
- Otherwise, the system collects data using the selected transmit antenna and the selected receive antenna (block2316) to produce
raw data 2318, the raw data collected at spatial location (xm,yn,) being denoted by {tilde over (Ψ)}mnp where the indices m, n are used to denote position in the grid of spatial locations where data has been collected, and p is an index ranging from 1 to P corresponding to the frequency fp at which the data was collected. The frequency step of the system is df=(fp-f1)(P−1). The unambiguous range of the radar is c/(2·df) where c denotes light speed in air. Each of the data points {tilde over (Ψ)}mnp is a complex number. In a preferred embodiment, a user may specify f1 and fp. - Once the data has been collected for the selected transmit antenna and the selected receive antenna, the system selects the next receive antenna (block2320) and returns to the beginning of the data collection loop (block 2312).
- Analyzing the
raw data 2318, as illustrated in FIG. 24, includes preconditioning theraw data 2318 to produce preconditioned data (block 2402), analyzing the preconditioned data (block 2404) and displaying images of the analyzed data (block 2406). - Preconditioning the raw data to produce preconditioned data includes removing a constant frequency component and a system travel time delay (block2408), removing a transmit-antenna to receive-antenna coupling effect (block 2410) and prewhitening (block 2412) for each spatial location of the raw data.
-
- Removing the transmit-antenna to receive-antenna coupling effect (block2410) is performed to minimize the effects of ground bounce. In the preferred embodiment, this is accomplished by using a spatial high pass filter which acts in the in track direction. This changes the system from an absolute return system to a relative return system.
- There are several techniques for implementing the spatial high pass filter. They include single row differencing, averaging techniques, Fourier techniques, and digital filtering. Single row differencing, which is the simplest method, includes applying the following equation:
- Ψmnp={circumflex over (Ψ)}mnp{circumflex over (Ψ)}mñp
- where {circumflex over (Ψ)}mñpis an in track reference scan.
- In a preferred embodiment, the user selects the in track reference scan.
-
- where N1 and N2 define a region to be imaged.
-
- where aq are digital filter coefficients chosen to reject low frequency spatial energy.
- Removing the transmit-antenna to receive-antenna coupling effect may be accomplished by performing spatial filtering in the cross track direction by applying the following equation:
- Ψmnp={circumflex over (Ψ)}mnp−{circumflex over (Ψ)}{tilde over (m)}np
- where {circumflex over (Ψ)}{tilde over (m)}np is a cross line reference scan.
- In the preferred embodiment, prewhitening (block2442) includes applying the following equation:
- γmnp=bp·Ψmnp
- where bp are frequency dependent weights.
- In the preferred embodiment, analyzing the preconditioned data (block2404) includes depth focusing (block 2414) and synthetic aperture radar (SAR) processing (block 2416).
-
-
-
-
- The phase weights in this equation do not depend on absolute sensor position and, in the preferred embodiment, they are precomputed and reused which greatly reduces the time required to perform SAR processing.
- In the preferred embodiment, the size of the SAR array is (2Ms+1)(2Ns+2). For Ms=1 and Ns=1, the array will contain 9 points.
- SAR processing (block2416) is an optional procedure which can be selected by the user to enhance the images produced by the system.
-
- where
- Imnw is the complex image value at spatial location (XF,m,YF,n,Zw);
- U is the SAR array size in the cross-track direction;
- V is SAR array size in the along track direction;
- (fp
1 , fp3 ) is the frequency processing band; -
- YF,m=0.933013d+(n−1)dy, preferably;
- d=5.52 inches, preferably; and
- dy=scan spacing.
-
- where (Xs,nVs,v) and (Xrn,Yrv) are the location of the transmit and receive antennas and cg is the speed of light in the ground.
- The preprocessing of the data (block2402) and the analysis of the preprocessed data (block 2404) produces a block of data representing a three
dimensional volume 2802 as illustrated in FIG. 28. The invention allows the display of a plan view of the analyzed data and a side view of the analyzed data. - Displaying images of the analyzed data (block2406) includes computing a plan view image of the analyzed data, computing a side view image of the analyzed data (block 2418) and displaying the plan view image and the side view image (block 2420).
- In the preferred embodiment, computing the plan view image of the analyzed data includes applying the following equation:
- PlanViewmn=maxw |I mnw|2
- where
- maxw is the maximum value across all w (depths).
- In the preferred embodiment, computing the side view image of the analyzed data includes applying the following equation:
- SideViewnw=maxm|Imnw|2
- where
- maxw is the maximum value across all w (depths).
- In other embodiments, the plan view image and the side view image are computed using other rendering techniques such as averaging or displaying only image values that fall within defined ranges. Any method or technique for presenting a side view and a plan view fall within the scope of the invention.
- The results of removing the transmit-antenna to receive-antenna coupling effects (block2410) and SAR processing (block 2416) are illustrated in FIGS. 29 and 30. FIG. 29 shows the data without removing the transmit-antenna to receive-antenna coupling effects (block 2410) and SAR processing (block 2416). FIG. 30 shows the same data after difference referencing using the first data row (block 2410) and after SAR processing (block 2416). A mine is visible at about
channel 7 and Y-Range 1.5 in the plan view and at a depth of about 6 and a Y-Range of about 1.5 in the side view. - A physical description of a preferred embodiment of the processing performed by the ground penetrating radar is illustrated in FIG. 31.
- The
GPSAR program 3102 resides on thecomputer 108 and provides user control of the ground penetrating radar data acquisition. It accepts a file name, dy (the in track distance between grid points) and array offset (the distance the two antenna arrays are offset from each other). - The
ENCODER program 3104 resides on thecomputer 108 and triggers the system to take data based on a signal from themovement detector 202. It also receives the scan data and stores it as rawbinary data 3106. - The
RADARBIN.C program 3108 resides on thecomputer 108 and converts the rawbinary data 3106 to ASCII data sorted byscan 3110. - The
MGPRVOL.F program 3112 resides on thecomputer 108 and computes avolumetric image 3114 from theASCII data 3110. - The
GPRIMAGE.C program 3116 renders aplan view 3118 and aside view 3120 from thevolumetric image 3114. - The
HTML User Interface 3122 is a web browser such as Netscape Navigator or Microsoft Explorer. - The programs stored in the
memory 908 for theprocessor 906 in the digital modules control the configuration of their respective components and accept as inputs the start frequency, stop frequency, number of frequency steps in the scan, and the dwell time at each scan step. These programs also provide an Ethernet interface to thecomputer 108 and serve up a web page that can be accessed from thecomputer 108 through theHTML user interface 3122. These programs also collect data from theDSP 1022. - The programs stored in the
program memory 1026 for theDSP 1022 control thedigital module RF module - While some of the components have been described as being implemented in hardware and others in software or firmware, it will be apparent to persons of ordinary skill in the art that some of the hardware portions could be implemented in software or firmware and that some of the software or firmware portions could be implemented in hardware. The software or firmware portions of the system may be written in machine language, assembly language or a higher order language, including such languages as C++, FORTRAN, JAVA or PEARL.
- Although several specific embodiments of the invention have been described and illustrated, the invention is not to be limited to the specific forms or arrangements of the parts steps so described and illustrated. The invention is limited only by the claims.
Claims (36)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US09/752,085 US6445334B1 (en) | 2000-12-29 | 2000-12-29 | Ground penetrating radar system |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US09/752,085 US6445334B1 (en) | 2000-12-29 | 2000-12-29 | Ground penetrating radar system |
Publications (2)
Publication Number | Publication Date |
---|---|
US6445334B1 US6445334B1 (en) | 2002-09-03 |
US20020122000A1 true US20020122000A1 (en) | 2002-09-05 |
Family
ID=25024804
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US09/752,085 Expired - Fee Related US6445334B1 (en) | 2000-12-29 | 2000-12-29 | Ground penetrating radar system |
Country Status (1)
Country | Link |
---|---|
US (1) | US6445334B1 (en) |
Cited By (46)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20030026382A1 (en) * | 2001-08-06 | 2003-02-06 | Shinichi Takeda | Radiation image taking apparatus |
US20050062639A1 (en) * | 2001-09-15 | 2005-03-24 | The Secretary Of State For Defence | Sub-surface radar imaging |
GB2406990A (en) * | 2003-10-06 | 2005-04-13 | Rst Gmbh | Ground penetrating radar |
US20050128125A1 (en) * | 2003-08-28 | 2005-06-16 | Jian Li | Land mine detector |
US20070139261A1 (en) * | 2005-12-19 | 2007-06-21 | Honeywell International Inc. | Systems and methods for self-calibrating a radar altimeter |
US20070211077A1 (en) * | 2006-03-09 | 2007-09-13 | Ambercore Software Inc. | Fast gridding of irregular data |
US20100189312A1 (en) * | 2008-10-02 | 2010-07-29 | Certusview Technologies, Llc | Methods and apparatus for overlaying electronic locate information on facilities map information and/or other image information displayed on a locate device |
US20100188088A1 (en) * | 2008-10-02 | 2010-07-29 | Certusview Technologies, Llc | Methods and apparatus for displaying and processing facilities map information and/or other image information on a locate device |
US20100188407A1 (en) * | 2008-10-02 | 2010-07-29 | Certusview Technologies, Llc | Methods and apparatus for displaying and processing facilities map information and/or other image information on a marking device |
US20120112957A1 (en) * | 2010-11-09 | 2012-05-10 | U.S. Government As Represented By The Secretary Of The Army | Multidirectional target detecting system and method |
US20120154207A1 (en) * | 2010-12-21 | 2012-06-21 | Electronics And Telecommunications Research Institute | Apparatus and method for scanning image in image processing system |
US20120262324A1 (en) * | 2011-04-13 | 2012-10-18 | Raytheon Company | Subterranean Image Generating Device And Associated Method |
US20120268307A1 (en) * | 2011-04-19 | 2012-10-25 | Raytheon Company | Systems and Methods for Mapping the Crust of the Earth |
US8306747B1 (en) * | 2007-01-19 | 2012-11-06 | Starodub, Inc. | Travel way measurement system |
US8311765B2 (en) | 2009-08-11 | 2012-11-13 | Certusview Technologies, Llc | Locating equipment communicatively coupled to or equipped with a mobile/portable device |
CN102799131A (en) * | 2012-08-08 | 2012-11-28 | 中国科学院东北地理与农业生态研究所 | Ground penetrating radar lower computer control system based on field programmable gate array (FPGA) |
EP2544022A1 (en) * | 2011-07-07 | 2013-01-09 | Astrium GmbH | Système de radar avec ouverture synthétique |
US8361543B2 (en) | 2008-10-02 | 2013-01-29 | Certusview Technologies, Llc | Methods and apparatus for displaying an electronic rendering of a marking operation based on an electronic record of marking information |
US8374789B2 (en) | 2007-04-04 | 2013-02-12 | Certusview Technologies, Llc | Systems and methods for using marking information to electronically display dispensing of markers by a marking system or marking tool |
US8401791B2 (en) | 2007-03-13 | 2013-03-19 | Certusview Technologies, Llc | Methods for evaluating operation of marking apparatus |
US8400155B2 (en) | 2008-10-02 | 2013-03-19 | Certusview Technologies, Llc | Methods and apparatus for displaying an electronic rendering of a locate operation based on an electronic record of locate information |
EP1722251A3 (en) * | 2005-05-10 | 2013-04-03 | HILTI Aktiengesellschaft | Hand guided scanning underground detector |
US8416995B2 (en) | 2008-02-12 | 2013-04-09 | Certusview Technologies, Llc | Electronic manifest of underground facility locate marks |
US8442766B2 (en) | 2008-10-02 | 2013-05-14 | Certusview Technologies, Llc | Marking apparatus having enhanced features for underground facility marking operations, and associated methods and systems |
US8473209B2 (en) | 2007-03-13 | 2013-06-25 | Certusview Technologies, Llc | Marking apparatus and marking methods using marking dispenser with machine-readable ID mechanism |
US8478617B2 (en) | 2008-10-02 | 2013-07-02 | Certusview Technologies, Llc | Methods and apparatus for generating alerts on a locate device, based on comparing electronic locate information to facilities map information and/or other image information |
US8478523B2 (en) | 2007-03-13 | 2013-07-02 | Certusview Technologies, Llc | Marking apparatus and methods for creating an electronic record of marking apparatus operations |
US8510141B2 (en) | 2008-10-02 | 2013-08-13 | Certusview Technologies, Llc | Methods and apparatus for generating alerts on a marking device, based on comparing electronic marking information to facilities map information and/or other image information |
US8583264B2 (en) | 2008-10-02 | 2013-11-12 | Certusview Technologies, Llc | Marking device docking stations and methods of using same |
US8620572B2 (en) | 2009-08-20 | 2013-12-31 | Certusview Technologies, Llc | Marking device with transmitter for triangulating location during locate operations |
US8620616B2 (en) | 2009-08-20 | 2013-12-31 | Certusview Technologies, Llc | Methods and apparatus for assessing marking operations based on acceleration information |
EP2717068A1 (en) * | 2012-10-08 | 2014-04-09 | MBDA UK Limited | Improvements in and relating to radar receivers |
WO2014057234A1 (en) * | 2012-10-08 | 2014-04-17 | Mbda Uk Limited | Improvements in and relating to radar receivers |
US20140121964A1 (en) * | 2012-10-25 | 2014-05-01 | Massachusetts Institute Of Technology | Vehicle localization using surface penetrating radar |
AU2005234659B2 (en) * | 2005-11-18 | 2014-05-29 | Tr Lord & Associates Limited | Method for a quality of outcome-based approach to the location, mapping, and management of sub-surface utilities |
US8749239B2 (en) | 2008-10-02 | 2014-06-10 | Certusview Technologies, Llc | Locate apparatus having enhanced features for underground facility locate operations, and associated methods and systems |
US8965700B2 (en) | 2008-10-02 | 2015-02-24 | Certusview Technologies, Llc | Methods and apparatus for generating an electronic record of environmental landmarks based on marking device actuations |
US9004004B2 (en) | 2008-07-10 | 2015-04-14 | Certusview Technologies, Llc | Optical sensing methods and apparatus for detecting a color of a marking substance |
US9097522B2 (en) | 2009-08-20 | 2015-08-04 | Certusview Technologies, Llc | Methods and marking devices with mechanisms for indicating and/or detecting marking material color |
US9177403B2 (en) | 2008-10-02 | 2015-11-03 | Certusview Technologies, Llc | Methods and apparatus for overlaying electronic marking information on facilities map information and/or other image information displayed on a marking device |
US9208458B2 (en) | 2008-10-02 | 2015-12-08 | Certusview Technologies, Llc | Methods and apparatus for analyzing locate and marking operations with respect to facilities maps |
WO2017186951A1 (en) * | 2016-04-29 | 2017-11-02 | RESCH, Claudia | Image generation device for local resolution measurement of surfaces |
KR101917374B1 (en) | 2017-02-10 | 2018-11-09 | 세종대학교산학협력단 | Apparatus and method for processing 3d ground penetrating radar signal |
US20190242678A1 (en) * | 2016-10-12 | 2019-08-08 | Yeditepe Universitesi | Mine sweeping vehicle |
US10663579B2 (en) | 2017-01-27 | 2020-05-26 | Massachusetts Institute Of Technology | Extending footprint for localization using surface penetrating radar (SPR) |
WO2023041176A1 (en) * | 2021-09-17 | 2023-03-23 | Proceq Sa | Method for probing a subsurface structure |
Families Citing this family (51)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20030012411A1 (en) * | 2001-07-13 | 2003-01-16 | Sjostrom Keith Jerome | System and method for displaying and collecting ground penetrating radar data |
WO2004029562A2 (en) * | 2002-09-27 | 2004-04-08 | Innovatum, Inc. | Apparatus and method using continuous -wave radiation for detecting and locating targets hidden behind a surface |
US6690316B2 (en) | 2002-09-27 | 2004-02-10 | The United States Of America As Represented By The Secretary Of The Army | System and method for automated alerting to geospatial anomalies |
US8174429B2 (en) * | 2003-08-15 | 2012-05-08 | L-3 Communications Cyterra Corporation | Mine detection |
US7042385B1 (en) | 2003-09-16 | 2006-05-09 | Niitek, Inc. | Non-intrusive inspection impulse radar antenna |
US7230564B2 (en) * | 2004-06-11 | 2007-06-12 | 1M International Corporation | Microwave method and system for material inspection |
US7467810B2 (en) * | 2004-07-07 | 2008-12-23 | Sensors & Software Inc. | Apparatus for transporting a sensor |
CA2531423C (en) * | 2004-12-23 | 2014-06-10 | Sensors & Software Inc. | Data acquisition for a ground penetrating radar system |
US8253619B2 (en) * | 2005-02-15 | 2012-08-28 | Techtronic Power Tools Technology Limited | Electromagnetic scanning imager |
US7535407B2 (en) * | 2005-03-15 | 2009-05-19 | Prairielands Energy Marketing, Inc. | Apparatus using continuous-wave radiation for detecting and locating targets hidden behind a surface |
US8374754B2 (en) | 2005-12-05 | 2013-02-12 | Niitek, Inc. | Apparatus for detecting subsurface objects with a reach-in arm |
US7683821B1 (en) | 2006-10-25 | 2010-03-23 | Niitek, Inc. | Sensor sweeper for detecting surface and subsurface objects |
DE102007015561B3 (en) * | 2007-03-29 | 2008-08-21 | Deutsches Zentrum für Luft- und Raumfahrt e.V. | Method for exploring ice area or dry area by radar echo-sounder sounding, involves investigating area by radar sensor or by multiple compatible radar sensors of same working wavelengths |
US8040272B1 (en) | 2007-04-20 | 2011-10-18 | Niitek Inc. | Leading and lagging aperture for linear ground penetrating radar array |
US8140217B2 (en) * | 2007-07-31 | 2012-03-20 | Niitek, Inc. | Damage control system and method for a vehicle-based sensor |
US8249306B2 (en) | 2008-03-18 | 2012-08-21 | Certusview Technologies, Llc | Virtual white lines for delimiting planned excavation sites |
US8280117B2 (en) * | 2008-03-18 | 2012-10-02 | Certusview Technologies, Llc | Virtual white lines for indicating planned excavation sites on electronic images |
US8566737B2 (en) | 2009-02-11 | 2013-10-22 | Certusview Technologies, Llc | Virtual white lines (VWL) application for indicating an area of planned excavation |
CA2691780C (en) | 2009-02-11 | 2015-09-22 | Certusview Technologies, Llc | Management system, and associated methods and apparatus, for providing automatic assesment of a locate operation |
US8832565B2 (en) * | 2009-02-11 | 2014-09-09 | Certusview Technologies, Llc | Methods and apparatus for controlling access to a virtual white line (VWL) image for an excavation project |
AU2010221750A1 (en) * | 2009-03-03 | 2011-10-13 | Herbert Duvoisin Iii | Detection of surface and buried objects |
US8593329B2 (en) * | 2010-03-17 | 2013-11-26 | Tialinx, Inc. | Hand-held see-through-the-wall imaging and unexploded ordnance (UXO) detection system |
KR101191387B1 (en) * | 2010-03-22 | 2012-10-18 | 한국전자통신연구원 | Two dimensional array antenna and object detection device using the same |
US8842035B2 (en) | 2010-04-08 | 2014-09-23 | L-3 Communications Security And Detection Systems, Inc. | Sensor head |
US8854248B2 (en) * | 2010-08-26 | 2014-10-07 | Lawrence Livermore National Security, Llc | Real-time system for imaging and object detection with a multistatic GPR array |
CN102288625A (en) * | 2011-07-08 | 2011-12-21 | 李典基 | Channel concrete lining nondestructive testing system |
US8786485B2 (en) * | 2011-08-30 | 2014-07-22 | Masachusetts Institute Of Technology | Mobile coherent change detection ground penetrating radar |
US20130113648A1 (en) * | 2011-09-30 | 2013-05-09 | L-3 Communications Cyterra Corporation | Sensor head |
US9581690B2 (en) * | 2012-06-28 | 2017-02-28 | Tektronix, Inc. | Almost real-time sampling for ground-penetrating radar impulses |
US8957809B2 (en) | 2012-07-18 | 2015-02-17 | Geophysical Survey Systems, Inc. | Merged ground penetrating radar display for multiple antennas |
US9715008B1 (en) * | 2013-03-20 | 2017-07-25 | Bentley Systems, Incorporated | Visualization of 3-D GPR data in augmented reality |
US9395437B2 (en) | 2013-06-06 | 2016-07-19 | The United States Of America, As Represented By The Secretary Of The Army | Moving multi-polarization multi-transmitter/receiver ground penetrating radar system and signal processing for buried target detection |
CN104749559B (en) * | 2013-12-27 | 2017-02-08 | 中国科学院电子学研究所 | FPGA chip-based ice-penetrating radar control method |
US9733353B1 (en) * | 2014-01-16 | 2017-08-15 | L-3 Communications Security And Detection Systems, Inc. | Offset feed antennas |
US10917163B2 (en) | 2014-08-15 | 2021-02-09 | SEAKR Engineering, Inc. | Integrated mixed-signal RF transceiver with ADC, DAC, and DSP and high-bandwidth coherent recombination |
ES2846828T3 (en) | 2014-08-15 | 2021-07-29 | Seakr Eng Inc | Integrated mixed signal ASIC with ADC, DAC and DSP |
US10007996B2 (en) | 2015-03-02 | 2018-06-26 | Lawrence Livermore National Security, Llc | System for detecting objects in streaming 3D images formed from data acquired with a medium penetrating sensor |
WO2016166752A1 (en) * | 2015-04-12 | 2016-10-20 | Dov Zahavi | Method and system for locating underground targets |
CN105223621B (en) * | 2015-10-13 | 2017-10-24 | 北京中交桥宇科技有限公司 | A kind of device and method of utilization gps signal triggering geological radar collection |
US10782115B2 (en) * | 2016-09-02 | 2020-09-22 | Hossein Karami Porzani | Detection of radial deformations of transformers |
US11400956B2 (en) | 2017-02-23 | 2022-08-02 | Richard Anthony Bishel | Vehicle guidance system |
US10983209B2 (en) * | 2017-04-12 | 2021-04-20 | Ellumen, Inc. | Accurate signal compensations for UWB radar imaging in dispersive medium |
US10845461B2 (en) * | 2017-06-06 | 2020-11-24 | The United States Of America As Represented By The Secretary Of The Army | System and method for enabling radar operation in a spectrally congested environment |
RU2707419C1 (en) * | 2019-02-28 | 2019-11-26 | Ярослав Викторович Савватеев | Method for georadiolocation sounding and device for its implementation |
TR201913129A2 (en) * | 2019-08-29 | 2021-03-22 | Tuerkiye Bilimsel Ve Teknolojik Arastirma Kurumu Tuebitak | METAL AND MINE DETECTOR FITTED TO THE VEHICLE |
CA3123517A1 (en) | 2019-09-18 | 2021-03-25 | Thales Canada Inc. | Method and system for high-integrity vehicle localization and speed determination |
US10649079B1 (en) * | 2019-11-05 | 2020-05-12 | Altec Industries, Inc. | Ground-directed radar system |
CN111562574B (en) * | 2020-05-22 | 2022-08-16 | 中国科学院空天信息创新研究院 | MIMO ground penetrating radar three-dimensional imaging method based on backward projection |
CN112180452B (en) * | 2020-09-23 | 2023-09-29 | 中国建筑第八工程局有限公司 | Underground pipeline buried depth estimation method based on ground penetrating radar and three-dimensional velocity spectrum |
CN114137517B (en) * | 2022-02-07 | 2022-05-10 | 北京中科蓝图科技有限公司 | Two-three-dimensional integrated road detection method and device |
CN117388932B (en) * | 2023-12-13 | 2024-03-19 | 浙江艺信科技有限公司 | Underground pipeline detection visual positioning device and method for complex environment |
Family Cites Families (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US2517951A (en) * | 1947-06-20 | 1950-08-08 | Hazeltine Research Inc | Wide-band directive antenna |
US4218678A (en) | 1979-05-11 | 1980-08-19 | Ensco, Inc. | Synthetic pulse radar including a microprocessor based controller |
US5248975A (en) * | 1991-06-26 | 1993-09-28 | Geophysical Survey Systems, Inc. | Ground probing radar with multiple antenna capability |
US5227799A (en) * | 1992-03-09 | 1993-07-13 | Mitsui Engineering & Shipbuilding Co., Ltd. | Road void search radar system |
US5325095A (en) | 1992-07-14 | 1994-06-28 | The United States Of America As Represented By The United States Department Of Energy | Stepped frequency ground penetrating radar |
US5499029A (en) | 1992-07-14 | 1996-03-12 | Eg&G Energy Measurements, Inc. | Wide band stepped frequency ground penetrating radar |
US5557277A (en) * | 1994-09-06 | 1996-09-17 | Gde Systems, Inc. | Method for locating leakage of substances from subterranean structures |
US5553407A (en) * | 1995-06-19 | 1996-09-10 | Vermeer Manufacturing Company | Excavator data acquisition and control system and method of use |
GB2304483B (en) * | 1995-08-18 | 2000-03-29 | London Electricity Plc | System for and method of determining the location of an object in a medium |
US5647439A (en) * | 1995-12-14 | 1997-07-15 | Caterpillar Inc. | Implement control system for locating a surface interface and removing a layer of material |
US5720354A (en) * | 1996-01-11 | 1998-02-24 | Vermeer Manufacturing Company | Trenchless underground boring system with boring tool location |
US5835054A (en) * | 1996-03-01 | 1998-11-10 | The Regents Of The University Of California | Ultra wideband ground penetrating radar imaging of heterogeneous solids |
JP2001526771A (en) * | 1996-04-16 | 2001-12-18 | エム. スンリン,ウィリアム | Material transmission imaging radar |
-
2000
- 2000-12-29 US US09/752,085 patent/US6445334B1/en not_active Expired - Fee Related
Cited By (85)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6972410B2 (en) * | 2001-08-06 | 2005-12-06 | Canon Kabushiki Kaisha | Radiation image taking apparatus |
US20030026382A1 (en) * | 2001-08-06 | 2003-02-06 | Shinichi Takeda | Radiation image taking apparatus |
US20050062639A1 (en) * | 2001-09-15 | 2005-03-24 | The Secretary Of State For Defence | Sub-surface radar imaging |
US7190302B2 (en) * | 2001-09-15 | 2007-03-13 | The Secretary Of State For Defence In Her Britannic Majesty's Government Of The United Kingdom Of Great Britain And Northern Ireland | Sub-surface radar imaging |
US20050128125A1 (en) * | 2003-08-28 | 2005-06-16 | Jian Li | Land mine detector |
US7173560B2 (en) * | 2003-08-28 | 2007-02-06 | University Of Florida Research Foundation, Inc. | Land mine detector |
GB2406990A (en) * | 2003-10-06 | 2005-04-13 | Rst Gmbh | Ground penetrating radar |
EP1722251A3 (en) * | 2005-05-10 | 2013-04-03 | HILTI Aktiengesellschaft | Hand guided scanning underground detector |
AU2005234659B2 (en) * | 2005-11-18 | 2014-05-29 | Tr Lord & Associates Limited | Method for a quality of outcome-based approach to the location, mapping, and management of sub-surface utilities |
US20070139261A1 (en) * | 2005-12-19 | 2007-06-21 | Honeywell International Inc. | Systems and methods for self-calibrating a radar altimeter |
US7705773B2 (en) * | 2005-12-19 | 2010-04-27 | Honeywell International Inc. | Self-calibrating a radar altimeter based on a simulated return signal |
US7864175B2 (en) * | 2006-03-09 | 2011-01-04 | Ambercore Software Inc | Fast gridding of irregular data |
US20070211077A1 (en) * | 2006-03-09 | 2007-09-13 | Ambercore Software Inc. | Fast gridding of irregular data |
US8306747B1 (en) * | 2007-01-19 | 2012-11-06 | Starodub, Inc. | Travel way measurement system |
US8775077B2 (en) | 2007-03-13 | 2014-07-08 | Certusview Technologies, Llc | Systems and methods for using location data to electronically display dispensing of markers by a marking system or marking tool |
US8700325B2 (en) | 2007-03-13 | 2014-04-15 | Certusview Technologies, Llc | Marking apparatus and methods for creating an electronic record of marking operations |
US9086277B2 (en) | 2007-03-13 | 2015-07-21 | Certusview Technologies, Llc | Electronically controlled marking apparatus and methods |
US8407001B2 (en) | 2007-03-13 | 2013-03-26 | Certusview Technologies, Llc | Systems and methods for using location data to electronically display dispensing of markers by a marking system or marking tool |
US8478523B2 (en) | 2007-03-13 | 2013-07-02 | Certusview Technologies, Llc | Marking apparatus and methods for creating an electronic record of marking apparatus operations |
US8473209B2 (en) | 2007-03-13 | 2013-06-25 | Certusview Technologies, Llc | Marking apparatus and marking methods using marking dispenser with machine-readable ID mechanism |
US8903643B2 (en) | 2007-03-13 | 2014-12-02 | Certusview Technologies, Llc | Hand-held marking apparatus with location tracking system and methods for logging geographic location of same |
US8401791B2 (en) | 2007-03-13 | 2013-03-19 | Certusview Technologies, Llc | Methods for evaluating operation of marking apparatus |
US8374789B2 (en) | 2007-04-04 | 2013-02-12 | Certusview Technologies, Llc | Systems and methods for using marking information to electronically display dispensing of markers by a marking system or marking tool |
US8386178B2 (en) | 2007-04-04 | 2013-02-26 | Certusview Technologies, Llc | Marking system and method |
US8416995B2 (en) | 2008-02-12 | 2013-04-09 | Certusview Technologies, Llc | Electronic manifest of underground facility locate marks |
US9004004B2 (en) | 2008-07-10 | 2015-04-14 | Certusview Technologies, Llc | Optical sensing methods and apparatus for detecting a color of a marking substance |
US8457893B2 (en) | 2008-10-02 | 2013-06-04 | Certusview Technologies, Llc | Methods and apparatus for generating an electronic record of a marking operation including service-related information and/or ticket information |
US8600526B2 (en) | 2008-10-02 | 2013-12-03 | Certusview Technologies, Llc | Marking device docking stations having mechanical docking and methods of using same |
US8361543B2 (en) | 2008-10-02 | 2013-01-29 | Certusview Technologies, Llc | Methods and apparatus for displaying an electronic rendering of a marking operation based on an electronic record of marking information |
US20100189312A1 (en) * | 2008-10-02 | 2010-07-29 | Certusview Technologies, Llc | Methods and apparatus for overlaying electronic locate information on facilities map information and/or other image information displayed on a locate device |
US8442766B2 (en) | 2008-10-02 | 2013-05-14 | Certusview Technologies, Llc | Marking apparatus having enhanced features for underground facility marking operations, and associated methods and systems |
US8965700B2 (en) | 2008-10-02 | 2015-02-24 | Certusview Technologies, Llc | Methods and apparatus for generating an electronic record of environmental landmarks based on marking device actuations |
US8467969B2 (en) | 2008-10-02 | 2013-06-18 | Certusview Technologies, Llc | Marking apparatus having operational sensors for underground facility marking operations, and associated methods and systems |
US9177403B2 (en) | 2008-10-02 | 2015-11-03 | Certusview Technologies, Llc | Methods and apparatus for overlaying electronic marking information on facilities map information and/or other image information displayed on a marking device |
US8478525B2 (en) | 2008-10-02 | 2013-07-02 | Certusview Technologies, Llc | Methods, apparatus, and systems for analyzing use of a marking device by a technician to perform an underground facility marking operation |
US8478617B2 (en) | 2008-10-02 | 2013-07-02 | Certusview Technologies, Llc | Methods and apparatus for generating alerts on a locate device, based on comparing electronic locate information to facilities map information and/or other image information |
US8476906B2 (en) | 2008-10-02 | 2013-07-02 | Certusview Technologies, Llc | Methods and apparatus for generating electronic records of locate operations |
US8478524B2 (en) | 2008-10-02 | 2013-07-02 | Certusview Technologies, Llc | Methods and apparatus for dispensing marking material in connection with underground facility marking operations based on environmental information and/or operational information |
US20100188088A1 (en) * | 2008-10-02 | 2010-07-29 | Certusview Technologies, Llc | Methods and apparatus for displaying and processing facilities map information and/or other image information on a locate device |
US8510141B2 (en) | 2008-10-02 | 2013-08-13 | Certusview Technologies, Llc | Methods and apparatus for generating alerts on a marking device, based on comparing electronic marking information to facilities map information and/or other image information |
US8527308B2 (en) * | 2008-10-02 | 2013-09-03 | Certusview Technologies, Llc | Methods and apparatus for overlaying electronic locate information on facilities map information and/or other image information displayed on a locate device |
US8577707B2 (en) | 2008-10-02 | 2013-11-05 | Certusview Technologies, Llc | Methods and apparatus for overlaying electronic locate information on facilities map information and/or other image information displayed on a locate device |
US8583264B2 (en) | 2008-10-02 | 2013-11-12 | Certusview Technologies, Llc | Marking device docking stations and methods of using same |
US8589202B2 (en) | 2008-10-02 | 2013-11-19 | Certusview Technologies, Llc | Methods and apparatus for displaying and processing facilities map information and/or other image information on a marking device |
US8589201B2 (en) | 2008-10-02 | 2013-11-19 | Certusview Technologies, Llc | Methods and apparatus for generating alerts on a locate device, based on comparing electronic locate information to facilities map information and/or other image information |
US8400155B2 (en) | 2008-10-02 | 2013-03-19 | Certusview Technologies, Llc | Methods and apparatus for displaying an electronic rendering of a locate operation based on an electronic record of locate information |
US8612148B2 (en) | 2008-10-02 | 2013-12-17 | Certusview Technologies, Llc | Marking apparatus configured to detect out-of-tolerance conditions in connection with underground facility marking operations, and associated methods and systems |
US20100188407A1 (en) * | 2008-10-02 | 2010-07-29 | Certusview Technologies, Llc | Methods and apparatus for displaying and processing facilities map information and/or other image information on a marking device |
US8770140B2 (en) | 2008-10-02 | 2014-07-08 | Certusview Technologies, Llc | Marking apparatus having environmental sensors and operations sensors for underground facility marking operations, and associated methods and systems |
US8766638B2 (en) | 2008-10-02 | 2014-07-01 | Certusview Technologies, Llc | Locate apparatus with location tracking system for receiving environmental information regarding underground facility marking operations, and associated methods and systems |
US8644965B2 (en) | 2008-10-02 | 2014-02-04 | Certusview Technologies, Llc | Marking device docking stations having security features and methods of using same |
US9208458B2 (en) | 2008-10-02 | 2015-12-08 | Certusview Technologies, Llc | Methods and apparatus for analyzing locate and marking operations with respect to facilities maps |
US9542863B2 (en) | 2008-10-02 | 2017-01-10 | Certusview Technologies, Llc | Methods and apparatus for generating output data streams relating to underground utility marking operations |
US9069094B2 (en) | 2008-10-02 | 2015-06-30 | Certusview Technologies, Llc | Locate transmitter configured to detect out-of-tolerance conditions in connection with underground facility locate operations, and associated methods and systems |
US9046621B2 (en) | 2008-10-02 | 2015-06-02 | Certusview Technologies, Llc | Locate apparatus configured to detect out-of-tolerance conditions in connection with underground facility locate operations, and associated methods and systems |
US8731830B2 (en) | 2008-10-02 | 2014-05-20 | Certusview Technologies, Llc | Marking apparatus for receiving environmental information regarding underground facility marking operations, and associated methods and systems |
US8749239B2 (en) | 2008-10-02 | 2014-06-10 | Certusview Technologies, Llc | Locate apparatus having enhanced features for underground facility locate operations, and associated methods and systems |
US8311765B2 (en) | 2009-08-11 | 2012-11-13 | Certusview Technologies, Llc | Locating equipment communicatively coupled to or equipped with a mobile/portable device |
US8620616B2 (en) | 2009-08-20 | 2013-12-31 | Certusview Technologies, Llc | Methods and apparatus for assessing marking operations based on acceleration information |
US8620572B2 (en) | 2009-08-20 | 2013-12-31 | Certusview Technologies, Llc | Marking device with transmitter for triangulating location during locate operations |
US9097522B2 (en) | 2009-08-20 | 2015-08-04 | Certusview Technologies, Llc | Methods and marking devices with mechanisms for indicating and/or detecting marking material color |
US8624773B2 (en) * | 2010-11-09 | 2014-01-07 | The United States Of America As Represented By The Secretary Of The Army | Multidirectional target detecting system and method |
US20120112957A1 (en) * | 2010-11-09 | 2012-05-10 | U.S. Government As Represented By The Secretary Of The Army | Multidirectional target detecting system and method |
US20120154207A1 (en) * | 2010-12-21 | 2012-06-21 | Electronics And Telecommunications Research Institute | Apparatus and method for scanning image in image processing system |
US8917199B2 (en) * | 2011-04-13 | 2014-12-23 | Raytheon Company | Subterranean image generating device and associated method |
US20120262324A1 (en) * | 2011-04-13 | 2012-10-18 | Raytheon Company | Subterranean Image Generating Device And Associated Method |
US8884805B2 (en) * | 2011-04-19 | 2014-11-11 | Raytheon Company | Systems and methods for mapping the crust of the earth |
US20120268307A1 (en) * | 2011-04-19 | 2012-10-25 | Raytheon Company | Systems and Methods for Mapping the Crust of the Earth |
US8957806B2 (en) | 2011-07-07 | 2015-02-17 | Astrium Gmbh | Radar system with synthetic aperture |
EP2544022A1 (en) * | 2011-07-07 | 2013-01-09 | Astrium GmbH | Système de radar avec ouverture synthétique |
CN102799131A (en) * | 2012-08-08 | 2012-11-28 | 中国科学院东北地理与农业生态研究所 | Ground penetrating radar lower computer control system based on field programmable gate array (FPGA) |
WO2014057234A1 (en) * | 2012-10-08 | 2014-04-17 | Mbda Uk Limited | Improvements in and relating to radar receivers |
EP2717068A1 (en) * | 2012-10-08 | 2014-04-09 | MBDA UK Limited | Improvements in and relating to radar receivers |
US8949024B2 (en) * | 2012-10-25 | 2015-02-03 | Massachusetts Institute Of Technology | Vehicle localization using surface penetrating radar |
US20140121964A1 (en) * | 2012-10-25 | 2014-05-01 | Massachusetts Institute Of Technology | Vehicle localization using surface penetrating radar |
WO2017186951A1 (en) * | 2016-04-29 | 2017-11-02 | RESCH, Claudia | Image generation device for local resolution measurement of surfaces |
US20190242678A1 (en) * | 2016-10-12 | 2019-08-08 | Yeditepe Universitesi | Mine sweeping vehicle |
US10794666B2 (en) * | 2016-10-12 | 2020-10-06 | Yeditepe Universitesi | Mine sweeping vehicle |
US10663579B2 (en) | 2017-01-27 | 2020-05-26 | Massachusetts Institute Of Technology | Extending footprint for localization using surface penetrating radar (SPR) |
US10725171B2 (en) | 2017-01-27 | 2020-07-28 | Massachusetts Institute Of Technology | Method and system for localization of a vehicle using surface penetrating radar |
US10746867B2 (en) | 2017-01-27 | 2020-08-18 | Massachusetts Institute Of Technology | Determining surface characteristics |
JP2022024024A (en) * | 2017-01-27 | 2022-02-08 | マサチューセッツ インスティテュート オブ テクノロジー | Specification of surface characteristics |
US11402493B2 (en) | 2017-01-27 | 2022-08-02 | Massachusetts Institute Of Technology | Determining surface characteristics |
KR101917374B1 (en) | 2017-02-10 | 2018-11-09 | 세종대학교산학협력단 | Apparatus and method for processing 3d ground penetrating radar signal |
WO2023041176A1 (en) * | 2021-09-17 | 2023-03-23 | Proceq Sa | Method for probing a subsurface structure |
Also Published As
Publication number | Publication date |
---|---|
US6445334B1 (en) | 2002-09-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US6445334B1 (en) | Ground penetrating radar system | |
EP0288578B1 (en) | Method of detecting object and apparatus therefor | |
CA1277749C (en) | High resolution imaging doppler interferometer | |
US6664914B2 (en) | Ground penetrating radar | |
US7170440B1 (en) | Linear FM radar | |
US9207316B2 (en) | Radar for rejecting and looking past surface reflections | |
Witten | Present state of the art in ground-penetrating radars for mine detection | |
US20040257268A1 (en) | Doppler complex FFT police radar with direction sensing capability | |
JPH0344590A (en) | Radar device having overall aperture on base of rotary antenna | |
AU2002321657A1 (en) | Sub-surface radar imaging | |
WO2003025613A2 (en) | Sub-surface radar imaging | |
CN1311861A (en) | Radar ice sounder with parallel doppler processing | |
Fischer et al. | Detection of antipersonnel mines by using the factorization method on multistatic ground-penetrating radar measurements | |
US20190129001A1 (en) | Radar System For Detecting Profiles Of Objects, Particularly In A Vicinity Of A Machine Work Tool | |
EP1744177A1 (en) | Radar system and method for locating and identifying objects by their non-linear echo signals | |
EP1287378B1 (en) | Noncoherent integrator for signals with high acceleration uncertainty | |
Groenenboom et al. | Data processing for a land-mine-detection-dedicated GPR | |
Bradley et al. | Mine detection with a multichannel stepped-frequency ground-penetrating radar | |
CN106154264B (en) | Multi-band reconfigurable underground radar profiler system | |
Rodríguez-Morales et al. | A compact, reconfigurable, multi-UWB radar for snow thickness evaluation and altimetry: Development and field trials | |
CN112014793A (en) | Two-dimensional direction finding system and method based on directional diagram reconfigurable antenna array | |
Palmer et al. | High‐resolution wind profiling using combined spatial and frequency domain interferometry | |
Swart et al. | Collision avoidance radar able to differentiate objects | |
RU2303279C1 (en) | Method and device for sub-surface radiolocation probing | |
Koppenjan et al. | A lightweight ground penetrating radar |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
AS | Assignment |
Owner name: PLANNING SYSTEMS INCORPORATED, VIRGINIA Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNORS:BRADLEY, MARSHALL R.;MCCUMMINS, ROBERT J., JR.;CROWE, MICHAEL P.;AND OTHERS;REEL/FRAME:011406/0652 Effective date: 20001222 |
|
FEPP | Fee payment procedure |
Free format text: PAYOR NUMBER ASSIGNED (ORIGINAL EVENT CODE: ASPN); ENTITY STATUS OF PATENT OWNER: LARGE ENTITY |
|
FEPP | Fee payment procedure |
Free format text: PAT HOLDER NO LONGER CLAIMS SMALL ENTITY STATUS, ENTITY STATUS SET TO UNDISCOUNTED (ORIGINAL EVENT CODE: STOL); ENTITY STATUS OF PATENT OWNER: LARGE ENTITY |
|
FPAY | Fee payment |
Year of fee payment: 4 |
|
SULP | Surcharge for late payment | ||
FEPP | Fee payment procedure |
Free format text: PETITION RELATED TO MAINTENANCE FEES FILED (ORIGINAL EVENT CODE: PMFP); ENTITY STATUS OF PATENT OWNER: LARGE ENTITY |
|
FEPP | Fee payment procedure |
Free format text: PETITION RELATED TO MAINTENANCE FEES GRANTED (ORIGINAL EVENT CODE: PMFG); ENTITY STATUS OF PATENT OWNER: LARGE ENTITY |
|
PRDP | Patent reinstated due to the acceptance of a late maintenance fee |
Effective date: 20070717 |
|
FPAY | Fee payment |
Year of fee payment: 8 |
|
REMI | Maintenance fee reminder mailed | ||
AS | Assignment |
Owner name: FOSTER-MILLER, INC., MASSACHUSETTS Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNOR:PLANNING SYSTEMS INCORPORATED;REEL/FRAME:032807/0643 Effective date: 20140331 |
|
LAPS | Lapse for failure to pay maintenance fees | ||
STCH | Information on status: patent discontinuation |
Free format text: PATENT EXPIRED DUE TO NONPAYMENT OF MAINTENANCE FEES UNDER 37 CFR 1.362 |
|
FP | Lapsed due to failure to pay maintenance fee |
Effective date: 20140903 |