US20090295617A1 - System, Method, and Computer Program Product Providing Three-Dimensional Visualization of Ground Penetrating Radar Data - Google Patents
System, Method, and Computer Program Product Providing Three-Dimensional Visualization of Ground Penetrating Radar Data Download PDFInfo
- Publication number
- US20090295617A1 US20090295617A1 US11/852,030 US85203007A US2009295617A1 US 20090295617 A1 US20090295617 A1 US 20090295617A1 US 85203007 A US85203007 A US 85203007A US 2009295617 A1 US2009295617 A1 US 2009295617A1
- Authority
- US
- United States
- Prior art keywords
- scan
- depth
- data
- connections
- scans
- 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
- 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/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
Definitions
- the present description relates, in general, to techniques for analyzing and displaying radar information.
- Modern Ground Penetrating Radar (GPR) systems typically analyze and display radar data for a user to interpret.
- the data can be displayed to the user in a variety of ways, usually in “waterfall” plots that include a Two-Dimensional (2D) rendering of a most recent scan.
- prior art plots usually take one of three forms—as A-Scans, B-Scans or C-Scans.
- An A-Scan represents a single channel of a scan, and it is a line trace similar to the display on an oscilloscope.
- An A-Scan is a one dimensional data set with its progression being in the depth or time dimension.
- FIG. 9 is an illustration of prior art display 900 .
- Display 900 includes down-track display 901 and cross-track display 902 .
- Display 900 shows the data from a radar system that takes a series of 2D scans over a length.
- Down-track B-Scan 901 is a single channel of the radar array at various down-track distances as the array moves forward in time.
- Cross-track B-Scan 920 represents the plane across the width of the linear array at a single down-track distance.
- a method to analyze and display radar data employs a mathematical connection algorithm to further eliminate depth bins and also to derive more useful information from the remaining depth bins. Such further elimination of depth bins may help to limit the amount of data being processed, thereby freeing up some amount of processing power of a computer system.
- embodiments employ a display wherein a down-track B-scan and a cross-track B-scan are arranged orthogonally to each other to create two sides of a rectangular prism. Within the rectangular prism is a 3D plot of anomaly data.
- Some embodiments include a more intuitive down-track B-scan that selects channels of interest over time to create a 2D view of what is generally expected to be the most relevant data from a survey volume.
- the creation of the down-track B-scan can be enhanced by using an algorithm that smoothes the movement from one scan to the next, thereby helping to focus the B-scan on items of interest while minimizing abrupt shifts due to high frequency noise.
- Such displays tend to be more intuitive to human users than the displays of the prior art, at the very least, due to the rendering of more relevant data in a 3D plot, often in real-time.
- FIG. 1 is an illustration of an exemplary system adapted according to one embodiment of the invention
- FIG. 2 is an illustration of an exemplary scan, which is common to conventional GPR systems
- FIG. 3 is an illustration of an exemplary data set, which is common to conventional GPR systems
- FIG. 4 is an illustration of an exemplary display, adapted according to one embodiment of the present invention.
- FIG. 5 is an illustration of an exemplary display adapted according to one embodiment of the invention.
- FIG. 6 is an illustration of an exemplary flow adapted according to one embodiment of the invention for analyzing and displaying radar data
- FIG. 7 is an illustration of an exemplary flow adapted according to one embodiment of the invention for displaying radar data
- FIG. 8 illustrates an example computer system 800 adapted according to embodiments of the present invention.
- FIG. 9 is an illustration of prior art display.
- FIG. 1 is an illustration of exemplary system 100 adapted according to one embodiment of the invention.
- System 100 includes transmitter/receiver (Tx/Rx) unit 101 , processing unit 102 , and display unit 103 and is operable to produce a real-time Three-Dimensional (3D) radar display.
- Tx/Rx transmitter/receiver
- processing unit 102 processing unit 102
- display unit 103 is operable to produce a real-time Three-Dimensional (3D) radar display.
- 3D Three-Dimensional
- Tx/Rx unit 101 includes a plurality of elements that radiate radar signals and receive returned signals.
- transmitter and receiver elements may be separate or may be integrated, e.g., transceiver elements.
- system 100 is a Ground Penetrating Radar (GPR) system, and its radiation is focused on ground survey area 104 . The radiation penetrates the ground and provides data with regard to survey volume 105 .
- GPR Ground Penetrating Radar
- the transmitter/receiver elements are arranged in one or more linear arrays parallel to the ground, with each of the elements corresponding to a portion of the width of survey area 104 .
- the time between a transmit pulse and a receive pulse gives the depth of the radar sample.
- Each of the transmit/receive elements captures a series of pulses over a time period in order to provide the depth of survey volume 105 .
- Each of the captured pulses may be referred to as a “sample.”
- a series of pulses captured by the same transmit/receive element over a depth range may be referred to as a “channel.”
- the collection of channels for each of the transmit/receive elements may be referred to as a “scan.”
- a scan is shaped like a curtain—essentially Two-Dimensional (2D) and covering the width and depth of survey volume 105 .
- a scan can be rendered to show a 2D cross section of volume 105 , as shown in FIGS. 4 , 5 , and 9 .
- System 100 can be placed on a moveable platform, such as on a vehicle, so that additional scans can be taken to traverse the length dimension of survey volume 105 . Additionally or alternatively, unit 101 may be tilted up or down with respect to the ground in order to cover the length dimension of survey volume 105 .
- Processing unit 102 receives the returned signals from Tx/Rx unit 101 and performs various algorithms on the data in order to derive useful information therefrom, as well as to create image data to be displayed on display 103 .
- Processing unit 102 decreases the number of samples in the returned signals so that the information therein can be processed with speed, while at the same time retaining enough data to provide an adequate degree of resolution so that smaller items of interest can be identified. Processes for decreasing the number of samples are discussed further below.
- FIG. 2 is an illustration of exemplary scan 200 , which is common to conventional GPR systems.
- scan 200 the individual samples are represented as rectangular prisms, with twelve channels across and ten samples down in each channel.
- the dimensions of scan 200 are for convenience of illustration only, as a scan taken by a real-world GPR system is likely to be at least one or two orders of magnitude larger than that shown in FIG. 2 .
- a single scan (e.g., scan 200 ) can be thought of as being similar to a playing card. If the card is divided into a grid of equal sized blocks, each block can each be imagined as a sample point. If many cards are stacked into a deck, an observer can only see the first card in the front and the tops of the rest of the cards. In a deck of playing cards, it is generally known that some portion of the top will be white and, therefore, clear of useful information, so that it is safe to cut a small portion of the tops off from the cards and discard them. In GPR systems, this is analogous to the ground surface. Then, the top edge of each of the playing cards contains a little bit of useful information, such as red or some black, depending on the suit of the card.
- the next step is to discard depth bin data blocks that are still relatively very weak, since they are generally not expected to contain useful information.
- an example system filters out depth bin values that are close to zero. To do this the depth bin values are scored using a normalization based on depth that takes into account that deeper samples will generally have a lower value return due to attenuation. Depth bins with normalized scores less than a certain threshold are discarded. In general, this step eliminates much clutter, but still leaves High Frequency (HF) noise as well as small underground features and layers.
- HF High Frequency
- FIG. 3 is an illustration of exemplary data set 300 , which is common to conventional GPR systems.
- Data set 300 is a simplified example of two scans after depth binning and eliminating samples with normalized scores less than a given threshold.
- Data set 300 includes 12 channels across and 10 depth bins down. It should be noted, however, that most real-world data sets generally contain at least 24 channels across and 30 depth bins down. In the blocks shown, the lightest and darkest squares have the highest scores, while the middle colors have lower scores.
- the blanks are included in FIG. 3 to illustrate depth bins that have been eliminated and will not be present in future figures.
- the depth bins in FIG. 3 are common to conventional GPR systems.
- Various embodiments of the present invention are operable to eliminate even more depth bins from the data in data set 300 and to derive more useful information from the data.
- connection uses a “connection” algorithm to further analyze the data in data set 300 . For instance, for every depth bin that received a score in the above-described process, connections are made to surrounding depth bins based on that score. A depth bin with a score of five, for instance, will be mathematically connected to every depth bin in a sphere of radial distance five or less from itself. The each bin's connections are counted For an individual depth bin, any given connection counts once. However, it does not matter if that connection was initiated to or from that depth bin.
- a count is maintained of all the internal connections (i.e., connections identified and established by that depth bin) and external connections (i.e., connections identified and established by a different depth bin) as well as a total number of connected bins.
- most of the connections for the low-scored depth bin are external (i.e., the connections indicate an algorithmic efficiency rather than an actual quality of the depth bin).
- this example embodiment To connect to depth bins in future scans, this example embodiment generates a delay list that creates a list of depth bins and scores in a given scan so that connections from those depth bins can be applied to depth bins of future scans. As future scans come in, connections from the delay list are applied, and connections from future scans to past scans are applied.
- the mathematical connections can be used to filter out noise and subsurface layers in the data display, which can interfere with the detection of items of interest (e.g., mines).
- the total number of connections for a given depth bin can be compared to a connection threshold.
- a lower threshold e.g., eight connections
- HF noise which shows up as a vertical column of high value samples in the data
- basic clutter e.g., subsurface layers and other large features
- subsurface layers and other large features will generally show up, causing difficulty in distinguishing items of interest from background clutter.
- a bush or other object on the surface can also cause clutter in the display because an algorithm that eliminates everything above ground level is rarely perfect and can be expected to fail to find the ground correctly, and the strong return value of the ground will cause it to connect with a great many nearby samples.
- a higher threshold e.g., twenty-five
- the majority of strong returns e.g., mines
- Weaker and smaller targets will generally not show up.
- the system uses the F1 algorithm to return deviation values for the various depth bins in the scans. If the F1 algorithm begins to show high deviation for a subsection of the survey volume, the system reduces the connectivity threshold. So if, for example, a current threshold value is fifty, the system can reduce the threshold to twenty-five or fifteen, depending on the deviation value. On the other hand, lower deviation for a subsection may result in an increase of the connectivity threshold. Sharper increases/decreases or high, sustained increases/decreases in deviation may result in greater adjustments of the connectivity threshold.
- the system may also adjust the connectivity threshold in order to optimize detection of items of interest. Depth bins that do not have connectivity scores higher than the threshold are not displayed (e.g., discarded, or just simply not displayed). Depth bins that meet or exceed the threshold are provided to a display unit (e.g., unit 103 of FIG. 1 ) so that they can be rendered in a humanly-perceptible manner.
- a display unit e.g., unit 103 of FIG. 1
- various embodiments of the invention may use two or more algorithms to determine which depth bins to render to the user.
- An advantage of some embodiments of the invention is that they can help to eliminate lag time between the processing of successive scans by employing the delay list, discussed above.
- the delay list allows systems to continue to create mathematical connections to and from a given scan even after the system has progressed to another scan.
- various embodiments of the invention are not limited to processing all or a subset of scans at one time.
- such feature can be advantageous in that they may allow for the use of less-advanced or less powerful processors. This can be important in, e.g., military applications wherein space is at a premium and/or more complex systems may not be able to withstand the mechanical forces inherent in vehicle-based radar.
- many embodiments of the present invention can display 3D radar data in real-time. By contrast F1 displays usually lag real-time by at least twenty-four scans.
- FIG. 4 is an illustration of exemplary display 400 , adapted according to one embodiment of the present invention.
- Display 400 includes samples 401 , which correspond to data processed according to the technique described above.
- samples 401 are referred to as “mine samples.”
- Display 400 also includes cross-track B-scan 410 and down-track B-scan 420 projected orthogonally with respect to each other to create two sides of a rectangular prism.
- Samples 401 are rendered within the rectangular prism so that they appear in three-dimensional relation to where they correspond in the survey volume. For instance, note that samples 401 appear to be horizontally offset (along the width dimension) with regard to the portion of B-scan 420 that indicates the same mine, yet it has the same depth and distance placement with regard to the mine return in B-scan 410 .
- Cross-track B-scan 410 is a rendering of the data that is generated by the current scan.
- Cross-track B-scan 410 can be thought of as a single cross section of the survey volume at a given distance. It is refreshed with each new scan.
- Down-track B-scan 420 is a collection of channels, each one taken from a successive scan. However, techniques according to various embodiments of the invention generate down-track B-scan 420 by selecting a channel of interest in each scan, such that down-track B-scan 420 is generally not a true cross-section of the survey volume.
- Down-track B-scan 420 reflects the fact that it may not be possible to show every channel at once, such that a choice should be made as to which channel in a given scan to display. Instead of picking a single channel (e.g., the middle channel), channels are selected in various embodiments based on an algorithm so that each scan results in a different channel's data being written to display 400 . According to one example, if the current channel is X n , then on a subsequent scan, the channel with the highest score for interest (e.g., after a check for largest single depth bin value) is set as the target channel T. However rather than simply changing to channel T, the value of X n is added to or subtracted from by the value F, wherein F is given by Equation 1.
- ⁇ 1 (1)
- X n+1 does not “jump” the full distance to T. This may help to avoid gross discontinuities in down-track B-scan 420 .
- Down-track B-scan 420 then becomes an amalgamation of different channel scans one after another.
- the factor F is set arbitrarily as a value that minimizes discontinuities while at the same time traversing quickly enough to ensure that at least part of a given signal of interest is always displayed.
- the factor F can be changed automatically or manually, as desired.
- F is chosen based on an expected mean deviation and it is set at a level so that the channel selection “slides” rather than “jumps.”
- a scan is already on a mine, but HF noise appears in a subsequent scan. The channel selector will move to a channel that is closer to the HF noise but will come back toward the mine when the HF noise abates. This may help to keep focus near the mine in the face of HF noise.
- channel indicator 402 denotes the channel being selected from the current scan.
- FIG. 5 is an illustration of exemplary display 500 adapted according to one embodiment of the invention.
- one or more prior art algorithms can be used in addition to the mathematical connection algorithm described above in order to eliminate or enhance the data.
- the additional use of such algorithms can also be applied to display techniques to create more meaningful displays.
- display 500 corresponds to display 400 ( FIG. 4 ) slightly advanced in time.
- mine data in FIG. 4 was monochrome
- mine data in FIG. 5 can be rendered so that it indicates an amount of confidence that can be had in each returned depth bin.
- a display unit uses the returned deviation information (discussed above) to indicate confidence values. Generally, a higher deviation from the background return corresponds to a higher confidence.
- the rectangles that make up samples 501 can then be colored according to the confidence associated with each.
- the color coding is made intuitive, such that higher confidence depth bins have higher color contrast with regard to scans 410 and 420 . Additionally or alternatively, depth bins with higher confidence levels may be increased in size relative to other depth bins.
- Color coding and increasing sizes of depth bins represent two reliable methods to make some depth bins stand out, thereby indicating high relevance.
- boxes can be drawn around a group of depth bins, such as around items 401 or 501 of FIGS. 4 and 5 .
- Such boxes can also be selectively generated based on deviations. For example, if a certain number of depth bins have a deviation that passes another threshold, a box may be rendered, thereby indicating very high relevance. Further, color coding of boxes can be used to indicate higher or lower relevance of the depth bins therein.
- An advantage of displays 400 and 500 over various prior art displays, is that they are generally more intuitive to human operators. For instance, the 3D placement of samples 401 helps a human user to quickly judge the placement of an item of interest within the survey volume. Further, the positioning of mine samples with respect to down-track B-scan 420 can allow a human user to contrast samples 401 with relevant channels of past scans. This can assist in determining whether an indicated item is an item of interest or is a false positive (e.g., noise, utility pipes, and the like). Nevertheless, various embodiments may add, omit, or change some features of the displays shown in FIGS. 4 and 5 . For example, it is possible to render down-track B-scan 420 and cross-track B-scan 410 from different vantage points. It is also possible to add a C-scan as a “floor” to the rectangular prism shown in those figures. In other words, various changes may be made to the illustrated embodiments without departing from the scope of the invention.
- FIG. 6 is an illustration of exemplary flow 600 adapted according to one embodiment of the invention for analyzing and displaying radar data.
- the process represented by flow 600 may be performed, for example, by a computer or other processor-based device that is part of a radar system.
- a first scan is received that includes a plurality of channels of data over a time period, wherein each of the channels includes a plurality of samples corresponding to a portion of a surveyed volume.
- the number of samples in the scan are reduced by aggregating every N number of samples in a channel into a depth bin, thereby creating a plurality of depth bins.
- N is a positive integer, such as five, ten, or twenty, and may be higher than twenty or lower than five, depending upon the application.
- the depth bins are scored according to return strength. For example, a higher score may indicate a higher return strength. A higher return strength generally indicates a higher relevance of a piece of data, especially when searching for items that are expected to reflect electromagnetic waves.
- the score may be normalized to account for expected attenuation for deeper samples. Further, the scoring can be set to any index, such as in relation to mean background return strength, such that a score of one corresponds to an average return.
- step 604 ones of the depth bins that fail to reach a return strength score threshold are eliminated. This helps to eliminate items that are generally not expected to be relevant.
- step 605 mathematical connections are created among the remaining depth bins based on the return strength. Since the return strength of a depth bin determines the spherical radial distance allowed for identifying connections, a higher value will generally allow connections to a much larger number of other depth bins. Connections are not only made within a 2D scan, but are also made among nearby scans, so that the connections span three dimensions. In some embodiments, a delay list is used to allow connection to future scans without requiring a variable lag in displaying the data.
- step 606 for each of the remaining depth bins, connections both to and from other depth bins are totaled.
- step 607 those depth bins that do not meet a threshold number of connections are suppressed, thereby creating a reduced scan.
- the threshold can be adjusted dynamically, as described above.
- steps 601 - 607 are repeated for subsequent scans, including dynamically changing the connection threshold based upon information from another prescan algorithm (e.g., F1), thereby creating a reduced data set for the volume of interest.
- F1 prescan algorithm
- step 609 data are displayed from the reduced scans.
- Example displays include those shown in FIGS. 4 and 5 . In many embodiments, the displays are provided in real-time.
- FIG. 7 is an illustration of exemplary flow 700 adapted according to one embodiment of the invention for displaying radar data.
- the process represented by flow 700 may be performed, for example, by a computer or other processor-based device that is part of a radar system.
- a real-time display of a current scan as a cross-track B-scan is generated.
- a down-track B-scan is generated, and it includes selected channels from past scans.
- Step 702 may include selecting a channel from each of said scans as the scans are generated. The selecting can include identifying a channel with a relatively strong depth bin in a current scan and using a sliding algorithm to select a channel that lies between the identified channel and a channel identified in a previous scan.
- step 703 the cross-track B-scan and the down-track B-scan are projected so that they appear in planes orthogonal to each other, such an arrangement can be seen in FIGS. 4 and 5 .
- a display of the data from the reduced scans is generated in 3D positions in relation to corresponding areas of the down-track B-scan.
- the 3D placement of the displayed data also corresponds to the location of the data within the survey volume.
- a visual indicator is provided to the display of the data from the reduced scans.
- the visual indicator is produced based upon operation of an independent algorithm (e.g., F1) that identifies items of interest.
- Methods 600 and 700 are shown as series of discrete steps. However, other embodiments of the invention may add, delete, repeat, modify and/or rearrange various portions of methods 600 and 700 .
- steps 701 - 705 may be performed at or near the same time.
- the various components of the display may be rendered in any order.
- readable media can include any medium that can store or transfer information.
- FIG. 8 illustrates an example computer system 800 adapted according to embodiments of the present invention. That is, computer system 800 comprises an example system on which embodiments of the present invention may be implemented (such as a computer in a vehicle-based GPR system).
- Central processing unit (CPU) 801 is coupled to system bus 802 .
- CPU 801 may be any general purpose CPU. However, the present invention is not restricted by the architecture of CPU 801 as long as CPU 801 supports the inventive operations as described herein.
- CPU 801 may execute the various logical instructions according to embodiments of the present invention. For example, CPU 801 may execute machine-level instructions according to the exemplary operational flows described above in conjunction with FIGS. 6 and 7 .
- Computer system 800 also preferably includes random access memory (RAM) 803 , which may be SRAM, DRAM, SDRAM, or the like.
- Computer system 800 preferably includes read-only memory (ROM) 804 which may be PROM, EPROM, EEPROM, or the like.
- RAM 803 and ROM 804 hold user and system data and programs, as is well known in the art.
- Computer system 800 also preferably includes input/output (I/O) adapter 805 , communications adapter 811 , user interface adapter 808 , and display adapter 809 .
- I/O adapter 805 , user interface adapter 808 , and/or communications adapter 811 may, in certain embodiments, enable a user to interact with computer system 800 in order to input information, such as instructions to change a vantage point, to apply color coding to data, or the like.
- I/O adapter 805 preferably connects to storage device(s) 806 , such as one or more of hard drive, compact disc (CD) drive, floppy disk drive, tape drive, etc. to computer system 800 .
- the storage devices may be utilized when RAM 803 is insufficient for the memory requirements associated with storing radar and program data.
- Communications adapter 811 is preferably adapted to couple computer system 800 to network 812 (e.g., the Internet, a wireless network, or other network).
- User interface adapter 808 couples user input devices, such as keyboard 813 , pointing device 807 , and microphone 814 and/or output devices, such as speaker(s) 815 to computer system 800 .
- Display adapter 809 is driven by CPU 801 to control the display on display device 810 to, for example, display the user interface (such as that of FIGS. 4 and 5 ) of embodiments of the present invention.
- the present invention is not limited to the architecture of system 800 .
- any suitable processor-based device may be utilized, including without limitation personal computers, laptop computers, handheld computing devices, computer workstations, and multi-processor servers.
- embodiments of the present invention may be implemented on application specific integrated circuits (ASICs) or very large scale integrated (VLSI) circuits.
- ASICs application specific integrated circuits
- VLSI very large scale integrated circuits
- persons of ordinary skill in the art may utilize any number of suitable structures capable of executing logical operations according to the embodiments of the present invention.
- various embodiments of the invention are not limited to military applications or vehicle based applications and may find use in, e.g., archaeology, construction, and/or the like.
Abstract
Description
- The present description relates, in general, to techniques for analyzing and displaying radar information.
- Modern Ground Penetrating Radar (GPR) systems typically analyze and display radar data for a user to interpret. The data can be displayed to the user in a variety of ways, usually in “waterfall” plots that include a Two-Dimensional (2D) rendering of a most recent scan. However, prior art plots usually take one of three forms—as A-Scans, B-Scans or C-Scans. An A-Scan represents a single channel of a scan, and it is a line trace similar to the display on an oscilloscope. An A-Scan is a one dimensional data set with its progression being in the depth or time dimension. A B-Scan is a 2D slice of either a down-track or a cross-track in the ground (down-track and cross-track are explained in more detail with regard to
FIG. 9 ). B-scans are often textured or color mapped. A C-Scan is similar to what many people call a “plan view” or a “top down view”, and it is produced by an algorithm that collapses the received energy up into a 2D representation of energy under the ground as seen from a bird's eye view. -
FIG. 9 is an illustration ofprior art display 900.Display 900 includes down-track display 901 andcross-track display 902.Display 900 shows the data from a radar system that takes a series of 2D scans over a length. Down-track B-Scan 901 is a single channel of the radar array at various down-track distances as the array moves forward in time. Cross-track B-Scan 920 represents the plane across the width of the linear array at a single down-track distance. - A disadvantage of prior art solutions is that they are not intuitive. For instance, a user often is unable to discern the Three-Dimensional (3D) shape of an anomaly (e.g., a land mine) from prior art displays, such as C-scans and
scans - An additional disadvantage of prior art systems is that they tend to process an enormous amount of data, making it difficult to produce a 3D display in real-time. For example, a commonly used algorithm, known as the F1 algorithm, reduces received data by dropping, e.g., eighty percent of the data and aggregating the remaining data into units called “depth bins.” The F1 algorithm then drops a subset of the depth bins that show a weak return and creates a C-scan therefrom. However, this still leaves enough data to strain the processing power of most military computing systems. As a result, systems that provide real-time radar data rarely are able to generate a 3D display. In fact, 3D GPR displays are currently limited to systems that are not real-time systems, such as archaeological systems that take hours to give final results.
- The present invention is directed to systems, methods, and computer program products which analyze data and display data using techniques more advanced than those found in the prior art. For example, according to one embodiment, a method to analyze and display radar data employs a mathematical connection algorithm to further eliminate depth bins and also to derive more useful information from the remaining depth bins. Such further elimination of depth bins may help to limit the amount of data being processed, thereby freeing up some amount of processing power of a computer system.
- Additionally or alternatively, other embodiments employ a display wherein a down-track B-scan and a cross-track B-scan are arranged orthogonally to each other to create two sides of a rectangular prism. Within the rectangular prism is a 3D plot of anomaly data. Some embodiments include a more intuitive down-track B-scan that selects channels of interest over time to create a 2D view of what is generally expected to be the most relevant data from a survey volume. The creation of the down-track B-scan can be enhanced by using an algorithm that smoothes the movement from one scan to the next, thereby helping to focus the B-scan on items of interest while minimizing abrupt shifts due to high frequency noise. Such displays tend to be more intuitive to human users than the displays of the prior art, at the very least, due to the rendering of more relevant data in a 3D plot, often in real-time.
- The foregoing has outlined rather broadly the features and technical advantages of the present invention in order that the detailed description of the invention that follows may be better understood. Additional features and advantages of the invention will be described hereinafter which form the subject of the claims of the invention. It should be appreciated by those skilled in the art that the conception and specific embodiment disclosed may be readily utilized as a basis for modifying or designing other structures for carrying out the same purposes of the present invention. It should also be realized by those skilled in the art that such equivalent constructions do not depart from the spirit and scope of the invention as set forth in the appended claims. The novel features which are believed to be characteristic of the invention, both as to its organization and method of operation, together with further objects and advantages will be better understood from the following description when considered in connection with the accompanying figures. It is to be expressly understood, however, that each of the figures is provided for the purpose of illustration and description only and is not intended as a definition of the limits of the present invention.
- For a more complete understanding of the present invention, reference is now made to the following descriptions taken in conjunction with the accompanying drawings, in which:
-
FIG. 1 is an illustration of an exemplary system adapted according to one embodiment of the invention; -
FIG. 2 is an illustration of an exemplary scan, which is common to conventional GPR systems; -
FIG. 3 is an illustration of an exemplary data set, which is common to conventional GPR systems; -
FIG. 4 is an illustration of an exemplary display, adapted according to one embodiment of the present invention; -
FIG. 5 is an illustration of an exemplary display adapted according to one embodiment of the invention; -
FIG. 6 is an illustration of an exemplary flow adapted according to one embodiment of the invention for analyzing and displaying radar data; -
FIG. 7 is an illustration of an exemplary flow adapted according to one embodiment of the invention for displaying radar data; -
FIG. 8 illustrates anexample computer system 800 adapted according to embodiments of the present invention; and -
FIG. 9 is an illustration of prior art display. -
FIG. 1 is an illustration ofexemplary system 100 adapted according to one embodiment of the invention.System 100 includes transmitter/receiver (Tx/Rx)unit 101,processing unit 102, anddisplay unit 103 and is operable to produce a real-time Three-Dimensional (3D) radar display. - Tx/
Rx unit 101 includes a plurality of elements that radiate radar signals and receive returned signals. In some embodiments, transmitter and receiver elements may be separate or may be integrated, e.g., transceiver elements. In this example,system 100 is a Ground Penetrating Radar (GPR) system, and its radiation is focused onground survey area 104. The radiation penetrates the ground and provides data with regard tosurvey volume 105. - In this example, the transmitter/receiver elements are arranged in one or more linear arrays parallel to the ground, with each of the elements corresponding to a portion of the width of
survey area 104. Further, in GPR systems, the time between a transmit pulse and a receive pulse gives the depth of the radar sample. Each of the transmit/receive elements captures a series of pulses over a time period in order to provide the depth ofsurvey volume 105. - Each of the captured pulses may be referred to as a “sample.” A series of pulses captured by the same transmit/receive element over a depth range may be referred to as a “channel.” The collection of channels for each of the transmit/receive elements may be referred to as a “scan.” A scan is shaped like a curtain—essentially Two-Dimensional (2D) and covering the width and depth of
survey volume 105. A scan can be rendered to show a 2D cross section ofvolume 105, as shown inFIGS. 4 , 5, and 9.System 100 can be placed on a moveable platform, such as on a vehicle, so that additional scans can be taken to traverse the length dimension ofsurvey volume 105. Additionally or alternatively,unit 101 may be tilted up or down with respect to the ground in order to cover the length dimension ofsurvey volume 105. -
Processing unit 102 receives the returned signals from Tx/Rx unit 101 and performs various algorithms on the data in order to derive useful information therefrom, as well as to create image data to be displayed ondisplay 103.Processing unit 102 decreases the number of samples in the returned signals so that the information therein can be processed with speed, while at the same time retaining enough data to provide an adequate degree of resolution so that smaller items of interest can be identified. Processes for decreasing the number of samples are discussed further below. -
FIG. 2 is an illustration ofexemplary scan 200, which is common to conventional GPR systems. Inscan 200, the individual samples are represented as rectangular prisms, with twelve channels across and ten samples down in each channel. However, the dimensions ofscan 200 are for convenience of illustration only, as a scan taken by a real-world GPR system is likely to be at least one or two orders of magnitude larger than that shown inFIG. 2 . - In
scan 200, some samples are stronger than others. In general, a stronger return indicates a more reflective object in a portion of the survey volume that corresponds to the sample; however, no scan is ideal, as most scans will include some amount of noise that manifests itself as strong samples. - To aid in understanding the radar data, an analogy may be useful. A single scan (e.g., scan 200) can be thought of as being similar to a playing card. If the card is divided into a grid of equal sized blocks, each block can each be imagined as a sample point. If many cards are stacked into a deck, an observer can only see the first card in the front and the tops of the rest of the cards. In a deck of playing cards, it is generally known that some portion of the top will be white and, therefore, clear of useful information, so that it is safe to cut a small portion of the tops off from the cards and discard them. In GPR systems, this is analogous to the ground surface. Then, the top edge of each of the playing cards contains a little bit of useful information, such as red or some black, depending on the suit of the card.
- Next for each card, a copy is made that is 1/10th of the size vertically (i.e., in the depth dimension) such that for every ten grid blocks, the two brightest/strongest colors are detected and added together and drawn in the corresponding block in the copy. The other portions are discarded, as it is assumed that they contain less relevant information. This is equivalent to taking “depth bins” of the data set.
- The next step is to discard depth bin data blocks that are still relatively very weak, since they are generally not expected to contain useful information. Thus, an example system filters out depth bin values that are close to zero. To do this the depth bin values are scored using a normalization based on depth that takes into account that deeper samples will generally have a lower value return due to attenuation. Depth bins with normalized scores less than a certain threshold are discarded. In general, this step eliminates much clutter, but still leaves High Frequency (HF) noise as well as small underground features and layers. The above-described steps eliminate much data, but in real-world systems, there is usually still a large amount of data, some of which contains little useful information.
-
FIG. 3 is an illustration ofexemplary data set 300, which is common to conventional GPR systems.Data set 300 is a simplified example of two scans after depth binning and eliminating samples with normalized scores less than a given threshold.Data set 300 includes 12 channels across and 10 depth bins down. It should be noted, however, that most real-world data sets generally contain at least 24 channels across and 30 depth bins down. In the blocks shown, the lightest and darkest squares have the highest scores, while the middle colors have lower scores. The blanks are included inFIG. 3 to illustrate depth bins that have been eliminated and will not be present in future figures. - As mentioned above, the depth bins in
FIG. 3 are common to conventional GPR systems. Various embodiments of the present invention are operable to eliminate even more depth bins from the data indata set 300 and to derive more useful information from the data. - One example embodiment uses a “connection” algorithm to further analyze the data in
data set 300. For instance, for every depth bin that received a score in the above-described process, connections are made to surrounding depth bins based on that score. A depth bin with a score of five, for instance, will be mathematically connected to every depth bin in a sphere of radial distance five or less from itself. The each bin's connections are counted For an individual depth bin, any given connection counts once. However, it does not matter if that connection was initiated to or from that depth bin. For each bin a count is maintained of all the internal connections (i.e., connections identified and established by that depth bin) and external connections (i.e., connections identified and established by a different depth bin) as well as a total number of connected bins. This means a depth bin with a normalized score of as low as one might end up highly connected if nearby depth bins have high scores. In such a scenario most of the connections for the low-scored depth bin are external (i.e., the connections indicate an algorithmic efficiency rather than an actual quality of the depth bin). - It can be seen from this example that even small increases in score can cause large increases in connectivity, thereby helping to show relationships between samples. To connect to depth bins in future scans, this example embodiment generates a delay list that creates a list of depth bins and scores in a given scan so that connections from those depth bins can be applied to depth bins of future scans. As future scans come in, connections from the delay list are applied, and connections from future scans to past scans are applied.
- The mathematical connections can be used to filter out noise and subsurface layers in the data display, which can interfere with the detection of items of interest (e.g., mines). The total number of connections for a given depth bin can be compared to a connection threshold. A lower threshold (e.g., eight connections) will tend to filter out most HF noise (which shows up as a vertical column of high value samples in the data) as well as basic clutter. However, subsurface layers and other large features will generally show up, causing difficulty in distinguishing items of interest from background clutter. A bush or other object on the surface can also cause clutter in the display because an algorithm that eliminates everything above ground level is rarely perfect and can be expected to fail to find the ground correctly, and the strong return value of the ground will cause it to connect with a great many nearby samples. By contrast, at a higher threshold (e.g., twenty-five) the majority of strong returns (e.g., mines) will be visible while all but the biggest layers and surface features will be thrown away. Weaker and smaller targets will generally not show up.
- Various embodiments of the present invention employ a process that adjusts the connection threshold to optimize detection of items of interest while minimizing false positives. In one example, a radar system according to one embodiment of the invention also runs another pre-screener algorithm. An example of such an algorithm is the F1 algorithm, developed by NIITEK™ (based on work by Professor Leslie Collins at Duke University), which is a least-means-squared pre-screener algorithm that produces a C-scan of the same survey volume. Additionally or alternatively, some embodiments may apply other pre-screener algorithms, such as Constant False Alarm Rate (CFAR) algorithms, signal conditioning algorithms, and the like.
- In this example, the system uses the F1 algorithm to return deviation values for the various depth bins in the scans. If the F1 algorithm begins to show high deviation for a subsection of the survey volume, the system reduces the connectivity threshold. So if, for example, a current threshold value is fifty, the system can reduce the threshold to twenty-five or fifteen, depending on the deviation value. On the other hand, lower deviation for a subsection may result in an increase of the connectivity threshold. Sharper increases/decreases or high, sustained increases/decreases in deviation may result in greater adjustments of the connectivity threshold.
- Accordingly, in this example, as the system progresses from one scan to the next, the system may also adjust the connectivity threshold in order to optimize detection of items of interest. Depth bins that do not have connectivity scores higher than the threshold are not displayed (e.g., discarded, or just simply not displayed). Depth bins that meet or exceed the threshold are provided to a display unit (e.g.,
unit 103 ofFIG. 1 ) so that they can be rendered in a humanly-perceptible manner. Thus, various embodiments of the invention may use two or more algorithms to determine which depth bins to render to the user. - An advantage of some embodiments of the invention is that they can help to eliminate lag time between the processing of successive scans by employing the delay list, discussed above. The delay list allows systems to continue to create mathematical connections to and from a given scan even after the system has progressed to another scan. In other words, various embodiments of the invention are not limited to processing all or a subset of scans at one time. When taken in combination with the fact that various embodiments also further reduce the data when compared to prior art systems, such feature can be advantageous in that they may allow for the use of less-advanced or less powerful processors. This can be important in, e.g., military applications wherein space is at a premium and/or more complex systems may not be able to withstand the mechanical forces inherent in vehicle-based radar. Further, many embodiments of the present invention can display 3D radar data in real-time. By contrast F1 displays usually lag real-time by at least twenty-four scans.
- Various embodiments described above are operable to provide data which indicates the presence of items of interest in a survey volume. While such data can be used with a variety of display techniques, various displays are shown herein that may be especially adapted for the use of such data.
FIG. 4 is an illustration ofexemplary display 400, adapted according to one embodiment of the present invention.Display 400 includessamples 401, which correspond to data processed according to the technique described above. One particular application of GPR systems is in detecting underground mines, and the example data indisplay 400 corresponds to a survey volume containing a known mine. Thus inFIG. 4 ,samples 401 are referred to as “mine samples.” -
Display 400 also includes cross-track B-scan 410 and down-track B-scan 420 projected orthogonally with respect to each other to create two sides of a rectangular prism.Samples 401 are rendered within the rectangular prism so that they appear in three-dimensional relation to where they correspond in the survey volume. For instance, note thatsamples 401 appear to be horizontally offset (along the width dimension) with regard to the portion of B-scan 420 that indicates the same mine, yet it has the same depth and distance placement with regard to the mine return in B-scan 410. - Cross-track B-
scan 410 is a rendering of the data that is generated by the current scan. Cross-track B-scan 410 can be thought of as a single cross section of the survey volume at a given distance. It is refreshed with each new scan. - Down-track B-
scan 420 is a collection of channels, each one taken from a successive scan. However, techniques according to various embodiments of the invention generate down-track B-scan 420 by selecting a channel of interest in each scan, such that down-track B-scan 420 is generally not a true cross-section of the survey volume. - Down-track B-
scan 420 reflects the fact that it may not be possible to show every channel at once, such that a choice should be made as to which channel in a given scan to display. Instead of picking a single channel (e.g., the middle channel), channels are selected in various embodiments based on an algorithm so that each scan results in a different channel's data being written todisplay 400. According to one example, if the current channel is Xn, then on a subsequent scan, the channel with the highest score for interest (e.g., after a check for largest single depth bin value) is set as the target channel T. However rather than simply changing to channel T, the value of Xn is added to or subtracted from by the value F, wherein F is given byEquation 1. -
F=(the number of channels)/24 or F=|X n −T| if |X n −T|<=1 (1) - Thus, Xn+1 does not “jump” the full distance to T. This may help to avoid gross discontinuities in down-track B-
scan 420. Down-track B-scan 420 then becomes an amalgamation of different channel scans one after another. The factor F is set arbitrarily as a value that minimizes discontinuities while at the same time traversing quickly enough to ensure that at least part of a given signal of interest is always displayed. The factor F can be changed automatically or manually, as desired. - If abrupt discontinuities are allowed, then HF noise (which could be caused by some periodic external interference) and other anomalies can threaten to pull focus away from items of interest. In some embodiments, F is chosen based on an expected mean deviation and it is set at a level so that the channel selection “slides” rather than “jumps.” In one example, according to an embodiment, a scan is already on a mine, but HF noise appears in a subsequent scan. The channel selector will move to a channel that is closer to the HF noise but will come back toward the mine when the HF noise abates. This may help to keep focus near the mine in the face of HF noise.
- In
display 400, there are very few abrupt discontinuities between channel selections, and the majority of the mine information is clearly visible insamples 401 and in down-track B-scan 420. Further, inFIG. 4 ,channel indicator 402 denotes the channel being selected from the current scan. -
FIG. 5 is an illustration ofexemplary display 500 adapted according to one embodiment of the invention. As mentioned earlier, one or more prior art algorithms can be used in addition to the mathematical connection algorithm described above in order to eliminate or enhance the data. The additional use of such algorithms can also be applied to display techniques to create more meaningful displays. For example,display 500 corresponds to display 400 (FIG. 4 ) slightly advanced in time. Whereas the mine data inFIG. 4 was monochrome, mine data inFIG. 5 can be rendered so that it indicates an amount of confidence that can be had in each returned depth bin. A display unit uses the returned deviation information (discussed above) to indicate confidence values. Generally, a higher deviation from the background return corresponds to a higher confidence. The rectangles that make upsamples 501 can then be colored according to the confidence associated with each. In one example, the color coding is made intuitive, such that higher confidence depth bins have higher color contrast with regard toscans - Color coding and increasing sizes of depth bins represent two reliable methods to make some depth bins stand out, thereby indicating high relevance. Additionally or alternatively, boxes (not shown) can be drawn around a group of depth bins, such as around
items FIGS. 4 and 5 . Such boxes can also be selectively generated based on deviations. For example, if a certain number of depth bins have a deviation that passes another threshold, a box may be rendered, thereby indicating very high relevance. Further, color coding of boxes can be used to indicate higher or lower relevance of the depth bins therein. - An advantage of
displays samples 401 helps a human user to quickly judge the placement of an item of interest within the survey volume. Further, the positioning of mine samples with respect to down-track B-scan 420 can allow a human user to contrastsamples 401 with relevant channels of past scans. This can assist in determining whether an indicated item is an item of interest or is a false positive (e.g., noise, utility pipes, and the like). Nevertheless, various embodiments may add, omit, or change some features of the displays shown inFIGS. 4 and 5 . For example, it is possible to render down-track B-scan 420 and cross-track B-scan 410 from different vantage points. It is also possible to add a C-scan as a “floor” to the rectangular prism shown in those figures. In other words, various changes may be made to the illustrated embodiments without departing from the scope of the invention. - Various embodiments of the present invention can be thought of as processes of analyzing and displaying radar data. For instance,
FIG. 6 is an illustration ofexemplary flow 600 adapted according to one embodiment of the invention for analyzing and displaying radar data. The process represented byflow 600 may be performed, for example, by a computer or other processor-based device that is part of a radar system. - In
step 601, a first scan is received that includes a plurality of channels of data over a time period, wherein each of the channels includes a plurality of samples corresponding to a portion of a surveyed volume. Instep 602, the number of samples in the scan are reduced by aggregating every N number of samples in a channel into a depth bin, thereby creating a plurality of depth bins. N is a positive integer, such as five, ten, or twenty, and may be higher than twenty or lower than five, depending upon the application. - In
step 603, the depth bins are scored according to return strength. For example, a higher score may indicate a higher return strength. A higher return strength generally indicates a higher relevance of a piece of data, especially when searching for items that are expected to reflect electromagnetic waves. In some examples, the score may be normalized to account for expected attenuation for deeper samples. Further, the scoring can be set to any index, such as in relation to mean background return strength, such that a score of one corresponds to an average return. - In
step 604, ones of the depth bins that fail to reach a return strength score threshold are eliminated. This helps to eliminate items that are generally not expected to be relevant. - In
step 605, mathematical connections are created among the remaining depth bins based on the return strength. Since the return strength of a depth bin determines the spherical radial distance allowed for identifying connections, a higher value will generally allow connections to a much larger number of other depth bins. Connections are not only made within a 2D scan, but are also made among nearby scans, so that the connections span three dimensions. In some embodiments, a delay list is used to allow connection to future scans without requiring a variable lag in displaying the data. - In
step 606, for each of the remaining depth bins, connections both to and from other depth bins are totaled. Instep 607, those depth bins that do not meet a threshold number of connections are suppressed, thereby creating a reduced scan. The threshold can be adjusted dynamically, as described above. - In
step 608, steps 601-607 are repeated for subsequent scans, including dynamically changing the connection threshold based upon information from another prescan algorithm (e.g., F1), thereby creating a reduced data set for the volume of interest. - In
step 609, data are displayed from the reduced scans. Example displays include those shown inFIGS. 4 and 5 . In many embodiments, the displays are provided in real-time. -
FIG. 7 is an illustration ofexemplary flow 700 adapted according to one embodiment of the invention for displaying radar data. The process represented byflow 700 may be performed, for example, by a computer or other processor-based device that is part of a radar system. Instep 701, a real-time display of a current scan as a cross-track B-scan is generated. Instep 702, a down-track B-scan is generated, and it includes selected channels from past scans. Step 702 may include selecting a channel from each of said scans as the scans are generated. The selecting can include identifying a channel with a relatively strong depth bin in a current scan and using a sliding algorithm to select a channel that lies between the identified channel and a channel identified in a previous scan. - In
step 703, the cross-track B-scan and the down-track B-scan are projected so that they appear in planes orthogonal to each other, such an arrangement can be seen inFIGS. 4 and 5 . - In
step 704, using the data from reduced scans, a display of the data from the reduced scans is generated in 3D positions in relation to corresponding areas of the down-track B-scan. In some embodiments, the 3D placement of the displayed data also corresponds to the location of the data within the survey volume. - In
step 705, a visual indicator is provided to the display of the data from the reduced scans. The visual indicator is produced based upon operation of an independent algorithm (e.g., F1) that identifies items of interest. -
Methods methods - When implemented via computer-executable instructions, various elements of embodiments of the present invention are in essence the software code defining the operations of such various elements. The executable instructions or software code may be obtained from a readable medium (e.g., a hard drive media, optical media, EPROM, EEPROM, tape media, cartridge media, flash memory, ROM, memory stick, and/or the like). In fact, readable media can include any medium that can store or transfer information.
-
FIG. 8 illustrates anexample computer system 800 adapted according to embodiments of the present invention. That is,computer system 800 comprises an example system on which embodiments of the present invention may be implemented (such as a computer in a vehicle-based GPR system). Central processing unit (CPU) 801 is coupled tosystem bus 802.CPU 801 may be any general purpose CPU. However, the present invention is not restricted by the architecture ofCPU 801 as long asCPU 801 supports the inventive operations as described herein.CPU 801 may execute the various logical instructions according to embodiments of the present invention. For example,CPU 801 may execute machine-level instructions according to the exemplary operational flows described above in conjunction withFIGS. 6 and 7 . -
Computer system 800 also preferably includes random access memory (RAM) 803, which may be SRAM, DRAM, SDRAM, or the like.Computer system 800 preferably includes read-only memory (ROM) 804 which may be PROM, EPROM, EEPROM, or the like.RAM 803 andROM 804 hold user and system data and programs, as is well known in the art. -
Computer system 800 also preferably includes input/output (I/O)adapter 805,communications adapter 811,user interface adapter 808, anddisplay adapter 809. I/O adapter 805,user interface adapter 808, and/orcommunications adapter 811 may, in certain embodiments, enable a user to interact withcomputer system 800 in order to input information, such as instructions to change a vantage point, to apply color coding to data, or the like. - I/
O adapter 805 preferably connects to storage device(s) 806, such as one or more of hard drive, compact disc (CD) drive, floppy disk drive, tape drive, etc. tocomputer system 800. The storage devices may be utilized whenRAM 803 is insufficient for the memory requirements associated with storing radar and program data.Communications adapter 811 is preferably adapted to couplecomputer system 800 to network 812 (e.g., the Internet, a wireless network, or other network).User interface adapter 808 couples user input devices, such askeyboard 813, pointingdevice 807, andmicrophone 814 and/or output devices, such as speaker(s) 815 tocomputer system 800.Display adapter 809 is driven byCPU 801 to control the display ondisplay device 810 to, for example, display the user interface (such as that ofFIGS. 4 and 5 ) of embodiments of the present invention. - It shall be appreciated that the present invention is not limited to the architecture of
system 800. For example, any suitable processor-based device may be utilized, including without limitation personal computers, laptop computers, handheld computing devices, computer workstations, and multi-processor servers. Moreover, embodiments of the present invention may be implemented on application specific integrated circuits (ASICs) or very large scale integrated (VLSI) circuits. In fact, persons of ordinary skill in the art may utilize any number of suitable structures capable of executing logical operations according to the embodiments of the present invention. In fact, various embodiments of the invention are not limited to military applications or vehicle based applications and may find use in, e.g., archaeology, construction, and/or the like. - Although the present invention and its advantages have been described in detail, it should be understood that various changes, substitutions and alterations can be made herein without departing from the spirit and scope of the invention as defined by the appended claims. Moreover, the scope of the present application is not intended to be limited to the particular embodiments of the process, machine, manufacture, composition of matter, means, methods and steps described in the specification. As one of ordinary skill in the art will readily appreciate from the disclosure of the present invention, processes, machines, manufacture, compositions of matter, means, methods, or steps, presently existing or later to be developed that perform substantially the same function or achieve substantially the same result as the corresponding embodiments described herein may be utilized according to the present invention. Accordingly, the appended claims are intended to include within their scope such processes, machines, manufacture, compositions of matter, means, methods, or steps.
Claims (22)
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US11/852,030 US7675454B2 (en) | 2007-09-07 | 2007-09-07 | System, method, and computer program product providing three-dimensional visualization of ground penetrating radar data |
PCT/US2008/072303 WO2009032467A1 (en) | 2007-09-07 | 2008-08-06 | System, method, and computer program product providing three-dimensional visualization of ground penetrating radar data |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US11/852,030 US7675454B2 (en) | 2007-09-07 | 2007-09-07 | System, method, and computer program product providing three-dimensional visualization of ground penetrating radar data |
Publications (2)
Publication Number | Publication Date |
---|---|
US20090295617A1 true US20090295617A1 (en) | 2009-12-03 |
US7675454B2 US7675454B2 (en) | 2010-03-09 |
Family
ID=40429273
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US11/852,030 Expired - Fee Related US7675454B2 (en) | 2007-09-07 | 2007-09-07 | System, method, and computer program product providing three-dimensional visualization of ground penetrating radar data |
Country Status (2)
Country | Link |
---|---|
US (1) | US7675454B2 (en) |
WO (1) | WO2009032467A1 (en) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100329463A1 (en) * | 2009-06-24 | 2010-12-30 | Cisco Technology, Inc. | Group key management for mobile ad-hoc networks |
US20140002292A1 (en) * | 2012-06-28 | 2014-01-02 | Tektronix, Inc. | Almost real-time sampling for ground-penetrating radar impulses |
CN104755956A (en) * | 2012-11-02 | 2015-07-01 | 艾普斯科有限公司 | Method and apparatus for inspection of cooling towers |
CN109283523A (en) * | 2018-08-01 | 2019-01-29 | 西安交通大学 | A kind of geological radar B-scan data processing method |
CN112232392A (en) * | 2020-09-29 | 2021-01-15 | 深圳安德空间技术有限公司 | Data interpretation and identification method for three-dimensional ground penetrating radar |
CN116539643A (en) * | 2023-03-16 | 2023-08-04 | 南京京烁雷达科技有限公司 | Method and system for identifying coal rock data by using radar |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8957809B2 (en) | 2012-07-18 | 2015-02-17 | Geophysical Survey Systems, Inc. | Merged ground penetrating radar display for multiple antennas |
WO2014055732A2 (en) | 2012-10-04 | 2014-04-10 | Whaley Brian A | Shieldings for metal detector heads and manufacturing methods thereof |
US10901079B2 (en) | 2018-01-24 | 2021-01-26 | David Ira Wilens | Portable penetrating radar |
JP2019196995A (en) | 2018-05-10 | 2019-11-14 | ルネサスエレクトロニクス株式会社 | Radar device and method of controlling the same |
Citations (90)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3639784A (en) * | 1971-01-04 | 1972-02-01 | Raytheon Co | Pulse generator with storage means to maintain output transistor in saturation after removal of trigger pulse |
US4070673A (en) * | 1976-09-03 | 1978-01-24 | Sperry Rand Corporation | Radar video digital processor |
US4438404A (en) * | 1982-01-04 | 1984-03-20 | Tektronix, Inc. | Signal sampling system |
US4439765A (en) * | 1980-11-19 | 1984-03-27 | Hughes Aircraft Company | Radar video processor |
US4581715A (en) * | 1970-12-28 | 1986-04-08 | Hyatt Gilbert P | Fourier transform processor |
US4590614A (en) * | 1983-01-28 | 1986-05-20 | Robert Bosch Gmbh | Dipole antenna for portable radio |
US4591858A (en) * | 1982-12-03 | 1986-05-27 | The Mitre Corporation | Beacon/radar video generator system for air traffic simulation |
US4678345A (en) * | 1986-05-01 | 1987-07-07 | Tektronix, Inc. | Equivalent time pseudorandom sampling system |
US4686655A (en) * | 1970-12-28 | 1987-08-11 | Hyatt Gilbert P | Filtering system for processing signature signals |
US4715000A (en) * | 1985-08-06 | 1987-12-22 | General Electric Company | Digital phase-locked loop and frequency measuring device |
US4760525A (en) * | 1986-06-10 | 1988-07-26 | The United States Of America As Represented By The Secretary Of The Air Force | Complex arithmetic vector processor for performing control function, scalar operation, and set-up of vector signal processing instruction |
US4833475A (en) * | 1986-01-27 | 1989-05-23 | Raytheon Company | Raster scan radar with true motion memory |
US5003562A (en) * | 1989-10-27 | 1991-03-26 | Ncr Corporation | Digital phase lock loop decoder |
US5053983A (en) * | 1971-04-19 | 1991-10-01 | Hyatt Gilbert P | Filter system having an adaptive control for updating filter samples |
US5115245A (en) * | 1990-09-04 | 1992-05-19 | Hughes Aircraft Company | Single substrate microwave radar transceiver including flip-chip integrated circuits |
US5192886A (en) * | 1990-03-15 | 1993-03-09 | Hewlett-Packard Company | Sub-nanosecond calibrated delay line structure |
US5243343A (en) * | 1990-12-03 | 1993-09-07 | Zeelan Technology, Inc. | Signal acquisition system utilizing ultra-wide time range time base |
US5315627A (en) * | 1993-02-22 | 1994-05-24 | Hewlett-Packard Company | Pseudo-random repetitive sampling of a signal |
US5351055A (en) * | 1989-11-21 | 1994-09-27 | Furuno Electric Co., Ltd. | Radar apparatus |
US5386215A (en) * | 1992-11-20 | 1995-01-31 | Massachusetts Institute Of Technology | Highly efficient planar antenna on a periodic dielectric structure |
US5420589A (en) * | 1993-06-07 | 1995-05-30 | Wells; C. T. | System for evaluating the inner medium characteristics of non-metallic materials |
US5424735A (en) * | 1992-09-15 | 1995-06-13 | Advanced Intelligence Inc. | Analogue to digital and digital to analogue converters |
US5451894A (en) * | 1993-02-24 | 1995-09-19 | Advanced Micro Devices, Inc. | Digital full range rotating phase shifter |
US5495260A (en) * | 1993-08-09 | 1996-02-27 | Motorola, Inc. | Printed circuit dipole antenna |
US5510800A (en) * | 1993-04-12 | 1996-04-23 | The Regents Of The University Of California | Time-of-flight radio location system |
US5523760A (en) * | 1993-04-12 | 1996-06-04 | The Regents Of The University Of California | Ultra-wideband receiver |
US5552793A (en) * | 1994-12-02 | 1996-09-03 | Hughes Missile Systems Company | Self calibrated act pulse compression system |
US5661490A (en) * | 1993-04-12 | 1997-08-26 | The Regents Of The University Of California | Time-of-flight radio location system |
US5748153A (en) * | 1994-11-08 | 1998-05-05 | Northrop Grumman Corporation | Flared conductor-backed coplanar waveguide traveling wave antenna |
US5805110A (en) * | 1994-12-19 | 1998-09-08 | The Regents Of The University Of California | Impulse radar with swept range gate |
US5900761A (en) * | 1995-01-24 | 1999-05-04 | Advantest Corporation | Timing generating circuit and method |
US5900833A (en) * | 1996-04-16 | 1999-05-04 | Zircon Corporation | Imaging radar suitable for material penetration |
US5939912A (en) * | 1997-06-18 | 1999-08-17 | Lsi Logic Corporation | Recovery circuit having long hold time and phase range |
US5969667A (en) * | 1997-10-16 | 1999-10-19 | Automotive Systems Laboratory, Inc. | Radar system |
US5986600A (en) * | 1998-01-22 | 1999-11-16 | Mcewan; Thomas E. | Pulsed RF oscillator and radar motion sensor |
US6002723A (en) * | 1997-08-15 | 1999-12-14 | Lockheed Martin Corporation | Digital receiver having programmable complex threshold detectors and threshold interpreters |
US6055287A (en) * | 1998-05-26 | 2000-04-25 | Mcewan; Thomas E. | Phase-comparator-less delay locked loop |
US6137433A (en) * | 1999-03-18 | 2000-10-24 | The United States Of America As Represented By The Secretary Of Commerce | Scatterometer with adaptable spatial resolution |
US6150863A (en) * | 1998-04-01 | 2000-11-21 | Xilinx, Inc. | User-controlled delay circuit for a programmable logic device |
US6211814B1 (en) * | 1999-09-28 | 2001-04-03 | Lockhead Martin Corporation | Method for software scan conversion of target identification data in a general purpose computer |
US6239764B1 (en) * | 1998-06-09 | 2001-05-29 | Samsung Electronics Co., Ltd. | Wideband microstrip dipole antenna array and method for forming such array |
US6249242B1 (en) * | 1998-08-07 | 2001-06-19 | Hitachi, Ltd. | High-frequency transmitter-receiver apparatus for such an application as vehicle-onboard radar system |
US6281833B1 (en) * | 1999-11-10 | 2001-08-28 | Science And Applied Technology, Inc. | Range rate aiding in a pulsed radar system |
US6329929B1 (en) * | 1998-12-21 | 2001-12-11 | Medtronic Inc. | Telemetry system with phase-locking noise suppressing receiver |
US20020000946A1 (en) * | 2000-06-30 | 2002-01-03 | Jyrki Portin | Antenna feed structure |
US6342866B1 (en) * | 2000-03-17 | 2002-01-29 | The United States Of America As Represented By The Secretary Of The Navy | Wideband antenna system |
US6345099B1 (en) * | 1998-05-22 | 2002-02-05 | S3 Incorporated | System and method for copy protecting computer graphics |
US6433720B1 (en) * | 2001-03-06 | 2002-08-13 | Furaxa, Inc. | Methods, apparatuses, and systems for sampling or pulse generation |
US6501413B2 (en) * | 1999-03-22 | 2002-12-31 | Sensors & Software Inc. | Timing and control and data acquisition for a multi transducer ground penetrating radar system |
US20030043078A1 (en) * | 2001-08-23 | 2003-03-06 | Ten-Long Deng | Antenna module |
US6538614B2 (en) * | 2001-04-17 | 2003-03-25 | Lucent Technologies Inc. | Broadband antenna structure |
US6580304B1 (en) * | 2002-03-28 | 2003-06-17 | M/A-Com, Inc. | Apparatus and method for introducing signal delay |
US20030179025A1 (en) * | 2002-03-19 | 2003-09-25 | Infineon Technologies North America Corp. | Delay lock loop having an edge detector and fixed delay |
US6650661B1 (en) * | 1998-06-15 | 2003-11-18 | International Business Machines Corporation | System that compensates for variances due to process and temperature changes |
US6657577B1 (en) * | 1997-07-02 | 2003-12-02 | Malaa Geoscience Forvaltning Ab | Radar plant and measurement technique for determination of the orientation and the depth of buried objects |
US6680634B1 (en) * | 2002-12-03 | 2004-01-20 | Nokia Corporation | Self calibrating digital delay-locked loop |
US6690741B1 (en) * | 1997-05-16 | 2004-02-10 | Multispectral Solutions, Inc. | Ultra wideband data transmission system and method |
US6694273B2 (en) * | 2001-03-02 | 2004-02-17 | Fujitsu Limited | Receiving apparatus, method of compensating for waveform degradation of received signal, apparatus and method for detecting waveform degradation, and apparatus and method for measuring waveform |
US20040036655A1 (en) * | 2002-08-22 | 2004-02-26 | Robert Sainati | Multi-layer antenna structure |
US6726146B2 (en) * | 2002-04-26 | 2004-04-27 | Singapore Technologies Aerospace Limited | Multiple model target tracking with variable sampling rate |
US20040090373A1 (en) * | 2002-11-08 | 2004-05-13 | Antonio Faraone | Multi-band antennas |
US20040111650A1 (en) * | 2002-12-10 | 2004-06-10 | Chien-Ming Chen | Apparatus for delay calibration of a forward clock |
US6778000B2 (en) * | 2002-04-12 | 2004-08-17 | Samsung Electronics Co., Ltd. | Integrated circuit devices that provide constant time delays irrespective of temperature variation |
US20040178838A1 (en) * | 2003-03-13 | 2004-09-16 | International Business Machines Corporation | Variable pulse width and pulse separation clock generator |
US6836239B2 (en) * | 2003-05-15 | 2004-12-28 | Aviation Communication & Surveillance Systems, Llc | Method and system for displaying text on a radial display |
US6845458B2 (en) * | 2000-08-31 | 2005-01-18 | Micron Technology, Inc. | System and method of operation of DLL and PLL to provide tight locking with large range, and dynamic tracking of PVT variations using interleaved delay lines |
US6853338B2 (en) * | 2000-06-30 | 2005-02-08 | Sirf Technology, Inc. | Wireless GPS apparatus with integral antenna device |
US6853227B2 (en) * | 2001-04-17 | 2005-02-08 | K-Tek Corporation | Controller for generating a periodic signal with an adjustable duty cycle |
US6864833B2 (en) * | 2003-04-11 | 2005-03-08 | Siemens Milltronics Process Instruments, Inc. | Time-base generator with self-compensating control loop |
US6885343B2 (en) * | 2002-09-26 | 2005-04-26 | Andrew Corporation | Stripline parallel-series-fed proximity-coupled cavity backed patch antenna array |
US6914468B2 (en) * | 2001-06-07 | 2005-07-05 | Koninklijke Philips Electronics N.V. | Controllable delay circuit for delaying an electrical signal |
US6930528B2 (en) * | 2002-02-21 | 2005-08-16 | Broadcom Corporation | Delay circuit and method with delay relatively independent of process, voltage, and temperature variations |
US20050200549A1 (en) * | 2004-03-15 | 2005-09-15 | Realtronics Corporation | Optimal Tapered Band Positioning to Mitigate Flare-End Ringing of Broadband Antennas |
US6956422B2 (en) * | 2003-03-17 | 2005-10-18 | Indiana University Research And Technology Corporation | Generation and measurement of timing delays by digital phase error compensation |
US20050237260A1 (en) * | 2004-04-23 | 2005-10-27 | Centurion Wireless Technologies, Inc. | Microstrip Antenna |
US20060038599A1 (en) * | 2004-08-20 | 2006-02-23 | Avants Bradley S | Systems and methods for implementing delay line circuitry |
US7026979B2 (en) * | 2003-07-03 | 2006-04-11 | Hrl Labortories, Llc | Method and apparatus for joint kinematic and feature tracking using probabilistic argumentation |
US7026850B2 (en) * | 2001-05-21 | 2006-04-11 | Acuid Corporation Limited | Programmable self-calibrating vernier and method |
US20060087471A1 (en) * | 2004-10-12 | 2006-04-27 | Hintz Kenneth J | Syntactic landmine detector |
US7037266B2 (en) * | 2002-04-25 | 2006-05-02 | Medtronic, Inc. | Ultrasound methods and implantable medical devices using same |
US7042385B1 (en) * | 2003-09-16 | 2006-05-09 | Niitek, Inc. | Non-intrusive inspection impulse radar antenna |
US7053814B2 (en) * | 2004-01-15 | 2006-05-30 | Hrl Laboratories, Llc | RF lightwave coding system for radar pulse compression |
US20060119407A1 (en) * | 2000-10-06 | 2006-06-08 | Abrosimov Igor A | Timing delay generator and method using temperature stabilisation |
US20060132210A1 (en) * | 2004-12-22 | 2006-06-22 | Lsi Logic Corporation | Wide-range programmable delay line |
US20060203613A1 (en) * | 2005-02-18 | 2006-09-14 | Leon Thomsen | System and method for using time-distance characteristics in acquisition, processing, and imaging of t-CSEM data |
US20060256025A1 (en) * | 2005-05-13 | 2006-11-16 | Realtronics Corporation | Machine Producible Directive Closed-Loop Impulse Antenna |
US7161531B1 (en) * | 2003-11-04 | 2007-01-09 | Northrop Grumman Corporation | High performance radar display |
US7203600B2 (en) * | 2003-12-19 | 2007-04-10 | Westerngeco, L.L.C. | Geophysical data processing |
US20070080864A1 (en) * | 2005-10-11 | 2007-04-12 | M/A-Com, Inc. | Broadband proximity-coupled cavity backed patch antenna |
US20080001808A1 (en) * | 2004-12-30 | 2008-01-03 | Passarelli Richard E Jr | System and method for processing data in weather radar |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5341141A (en) * | 1993-03-09 | 1994-08-23 | Hughes Missile Systems Company | Three dimensional imaging radar |
CA2266222C (en) | 1998-03-20 | 2009-01-06 | Sensors & Software Inc. | Timing and control and data acquisition for a multitransducer ground penetrating radar system |
JP4416580B2 (en) | 2004-06-28 | 2010-02-17 | 株式会社リコー | Delay control device |
-
2007
- 2007-09-07 US US11/852,030 patent/US7675454B2/en not_active Expired - Fee Related
-
2008
- 2008-08-06 WO PCT/US2008/072303 patent/WO2009032467A1/en active Application Filing
Patent Citations (100)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4686655A (en) * | 1970-12-28 | 1987-08-11 | Hyatt Gilbert P | Filtering system for processing signature signals |
US4581715A (en) * | 1970-12-28 | 1986-04-08 | Hyatt Gilbert P | Fourier transform processor |
US3639784A (en) * | 1971-01-04 | 1972-02-01 | Raytheon Co | Pulse generator with storage means to maintain output transistor in saturation after removal of trigger pulse |
US5053983A (en) * | 1971-04-19 | 1991-10-01 | Hyatt Gilbert P | Filter system having an adaptive control for updating filter samples |
US4070673A (en) * | 1976-09-03 | 1978-01-24 | Sperry Rand Corporation | Radar video digital processor |
US4439765A (en) * | 1980-11-19 | 1984-03-27 | Hughes Aircraft Company | Radar video processor |
US4438404A (en) * | 1982-01-04 | 1984-03-20 | Tektronix, Inc. | Signal sampling system |
US4591858A (en) * | 1982-12-03 | 1986-05-27 | The Mitre Corporation | Beacon/radar video generator system for air traffic simulation |
US4590614A (en) * | 1983-01-28 | 1986-05-20 | Robert Bosch Gmbh | Dipole antenna for portable radio |
US4715000A (en) * | 1985-08-06 | 1987-12-22 | General Electric Company | Digital phase-locked loop and frequency measuring device |
US4833475A (en) * | 1986-01-27 | 1989-05-23 | Raytheon Company | Raster scan radar with true motion memory |
US4678345A (en) * | 1986-05-01 | 1987-07-07 | Tektronix, Inc. | Equivalent time pseudorandom sampling system |
US4760525A (en) * | 1986-06-10 | 1988-07-26 | The United States Of America As Represented By The Secretary Of The Air Force | Complex arithmetic vector processor for performing control function, scalar operation, and set-up of vector signal processing instruction |
US5003562A (en) * | 1989-10-27 | 1991-03-26 | Ncr Corporation | Digital phase lock loop decoder |
US5351055A (en) * | 1989-11-21 | 1994-09-27 | Furuno Electric Co., Ltd. | Radar apparatus |
US5351055B1 (en) * | 1989-11-21 | 1997-02-25 | Furuno Electric Co | Radar apparatus |
US5192886A (en) * | 1990-03-15 | 1993-03-09 | Hewlett-Packard Company | Sub-nanosecond calibrated delay line structure |
US5420531A (en) * | 1990-03-15 | 1995-05-30 | Hewlett-Packard Company | Sub-nanosecond calibrated delay line structure |
US5115245A (en) * | 1990-09-04 | 1992-05-19 | Hughes Aircraft Company | Single substrate microwave radar transceiver including flip-chip integrated circuits |
US5243343A (en) * | 1990-12-03 | 1993-09-07 | Zeelan Technology, Inc. | Signal acquisition system utilizing ultra-wide time range time base |
US5444459A (en) * | 1990-12-03 | 1995-08-22 | Zeelan Technology, Inc. | Signal acquisition system utilizing ultra-wide time range time base |
US5424735A (en) * | 1992-09-15 | 1995-06-13 | Advanced Intelligence Inc. | Analogue to digital and digital to analogue converters |
US5386215A (en) * | 1992-11-20 | 1995-01-31 | Massachusetts Institute Of Technology | Highly efficient planar antenna on a periodic dielectric structure |
US5315627A (en) * | 1993-02-22 | 1994-05-24 | Hewlett-Packard Company | Pseudo-random repetitive sampling of a signal |
US5451894A (en) * | 1993-02-24 | 1995-09-19 | Advanced Micro Devices, Inc. | Digital full range rotating phase shifter |
US5510800A (en) * | 1993-04-12 | 1996-04-23 | The Regents Of The University Of California | Time-of-flight radio location system |
US5523760A (en) * | 1993-04-12 | 1996-06-04 | The Regents Of The University Of California | Ultra-wideband receiver |
US5661490A (en) * | 1993-04-12 | 1997-08-26 | The Regents Of The University Of California | Time-of-flight radio location system |
US5420589A (en) * | 1993-06-07 | 1995-05-30 | Wells; C. T. | System for evaluating the inner medium characteristics of non-metallic materials |
US5495260A (en) * | 1993-08-09 | 1996-02-27 | Motorola, Inc. | Printed circuit dipole antenna |
US5748153A (en) * | 1994-11-08 | 1998-05-05 | Northrop Grumman Corporation | Flared conductor-backed coplanar waveguide traveling wave antenna |
US5552793A (en) * | 1994-12-02 | 1996-09-03 | Hughes Missile Systems Company | Self calibrated act pulse compression system |
US5805110A (en) * | 1994-12-19 | 1998-09-08 | The Regents Of The University Of California | Impulse radar with swept range gate |
US5900761A (en) * | 1995-01-24 | 1999-05-04 | Advantest Corporation | Timing generating circuit and method |
US5900833A (en) * | 1996-04-16 | 1999-05-04 | Zircon Corporation | Imaging radar suitable for material penetration |
US6690741B1 (en) * | 1997-05-16 | 2004-02-10 | Multispectral Solutions, Inc. | Ultra wideband data transmission system and method |
US5939912A (en) * | 1997-06-18 | 1999-08-17 | Lsi Logic Corporation | Recovery circuit having long hold time and phase range |
US6657577B1 (en) * | 1997-07-02 | 2003-12-02 | Malaa Geoscience Forvaltning Ab | Radar plant and measurement technique for determination of the orientation and the depth of buried objects |
US6002723A (en) * | 1997-08-15 | 1999-12-14 | Lockheed Martin Corporation | Digital receiver having programmable complex threshold detectors and threshold interpreters |
US5969667A (en) * | 1997-10-16 | 1999-10-19 | Automotive Systems Laboratory, Inc. | Radar system |
US5986600A (en) * | 1998-01-22 | 1999-11-16 | Mcewan; Thomas E. | Pulsed RF oscillator and radar motion sensor |
US6150863A (en) * | 1998-04-01 | 2000-11-21 | Xilinx, Inc. | User-controlled delay circuit for a programmable logic device |
US6345099B1 (en) * | 1998-05-22 | 2002-02-05 | S3 Incorporated | System and method for copy protecting computer graphics |
US6055287A (en) * | 1998-05-26 | 2000-04-25 | Mcewan; Thomas E. | Phase-comparator-less delay locked loop |
US6239764B1 (en) * | 1998-06-09 | 2001-05-29 | Samsung Electronics Co., Ltd. | Wideband microstrip dipole antenna array and method for forming such array |
US6650661B1 (en) * | 1998-06-15 | 2003-11-18 | International Business Machines Corporation | System that compensates for variances due to process and temperature changes |
US6249242B1 (en) * | 1998-08-07 | 2001-06-19 | Hitachi, Ltd. | High-frequency transmitter-receiver apparatus for such an application as vehicle-onboard radar system |
US6329929B1 (en) * | 1998-12-21 | 2001-12-11 | Medtronic Inc. | Telemetry system with phase-locking noise suppressing receiver |
US6137433A (en) * | 1999-03-18 | 2000-10-24 | The United States Of America As Represented By The Secretary Of Commerce | Scatterometer with adaptable spatial resolution |
US6501413B2 (en) * | 1999-03-22 | 2002-12-31 | Sensors & Software Inc. | Timing and control and data acquisition for a multi transducer ground penetrating radar system |
US6211814B1 (en) * | 1999-09-28 | 2001-04-03 | Lockhead Martin Corporation | Method for software scan conversion of target identification data in a general purpose computer |
US6281833B1 (en) * | 1999-11-10 | 2001-08-28 | Science And Applied Technology, Inc. | Range rate aiding in a pulsed radar system |
US6342866B1 (en) * | 2000-03-17 | 2002-01-29 | The United States Of America As Represented By The Secretary Of The Navy | Wideband antenna system |
US6853338B2 (en) * | 2000-06-30 | 2005-02-08 | Sirf Technology, Inc. | Wireless GPS apparatus with integral antenna device |
US20020000946A1 (en) * | 2000-06-30 | 2002-01-03 | Jyrki Portin | Antenna feed structure |
US6845458B2 (en) * | 2000-08-31 | 2005-01-18 | Micron Technology, Inc. | System and method of operation of DLL and PLL to provide tight locking with large range, and dynamic tracking of PVT variations using interleaved delay lines |
US7020794B2 (en) * | 2000-08-31 | 2006-03-28 | Micron Technology, Inc. | Interleaved delay line for phase locked and delay locked loops |
US6912666B2 (en) * | 2000-08-31 | 2005-06-28 | Micron Technology, Inc. | Interleaved delay line for phase locked and delay locked loops |
US6868504B1 (en) * | 2000-08-31 | 2005-03-15 | Micron Technology, Inc. | Interleaved delay line for phase locked and delay locked loops |
US6845459B2 (en) * | 2000-08-31 | 2005-01-18 | Micron Technology, Inc. | System and method to provide tight locking for DLL and PLL with large range, and dynamic tracking of PVT variations using interleaved delay lines |
US20060119407A1 (en) * | 2000-10-06 | 2006-06-08 | Abrosimov Igor A | Timing delay generator and method using temperature stabilisation |
US6694273B2 (en) * | 2001-03-02 | 2004-02-17 | Fujitsu Limited | Receiving apparatus, method of compensating for waveform degradation of received signal, apparatus and method for detecting waveform degradation, and apparatus and method for measuring waveform |
US6433720B1 (en) * | 2001-03-06 | 2002-08-13 | Furaxa, Inc. | Methods, apparatuses, and systems for sampling or pulse generation |
US6853227B2 (en) * | 2001-04-17 | 2005-02-08 | K-Tek Corporation | Controller for generating a periodic signal with an adjustable duty cycle |
US6538614B2 (en) * | 2001-04-17 | 2003-03-25 | Lucent Technologies Inc. | Broadband antenna structure |
US7026850B2 (en) * | 2001-05-21 | 2006-04-11 | Acuid Corporation Limited | Programmable self-calibrating vernier and method |
US6914468B2 (en) * | 2001-06-07 | 2005-07-05 | Koninklijke Philips Electronics N.V. | Controllable delay circuit for delaying an electrical signal |
US20030043078A1 (en) * | 2001-08-23 | 2003-03-06 | Ten-Long Deng | Antenna module |
US6930528B2 (en) * | 2002-02-21 | 2005-08-16 | Broadcom Corporation | Delay circuit and method with delay relatively independent of process, voltage, and temperature variations |
US20030179025A1 (en) * | 2002-03-19 | 2003-09-25 | Infineon Technologies North America Corp. | Delay lock loop having an edge detector and fixed delay |
US6798258B2 (en) * | 2002-03-28 | 2004-09-28 | M/A Com, Inc. | Apparatus and method for introducing signal delay |
US6580304B1 (en) * | 2002-03-28 | 2003-06-17 | M/A-Com, Inc. | Apparatus and method for introducing signal delay |
US6778000B2 (en) * | 2002-04-12 | 2004-08-17 | Samsung Electronics Co., Ltd. | Integrated circuit devices that provide constant time delays irrespective of temperature variation |
US7037266B2 (en) * | 2002-04-25 | 2006-05-02 | Medtronic, Inc. | Ultrasound methods and implantable medical devices using same |
US6726146B2 (en) * | 2002-04-26 | 2004-04-27 | Singapore Technologies Aerospace Limited | Multiple model target tracking with variable sampling rate |
US20040036655A1 (en) * | 2002-08-22 | 2004-02-26 | Robert Sainati | Multi-layer antenna structure |
US6885343B2 (en) * | 2002-09-26 | 2005-04-26 | Andrew Corporation | Stripline parallel-series-fed proximity-coupled cavity backed patch antenna array |
US20040090373A1 (en) * | 2002-11-08 | 2004-05-13 | Antonio Faraone | Multi-band antennas |
US6680634B1 (en) * | 2002-12-03 | 2004-01-20 | Nokia Corporation | Self calibrating digital delay-locked loop |
US20040111650A1 (en) * | 2002-12-10 | 2004-06-10 | Chien-Ming Chen | Apparatus for delay calibration of a forward clock |
US20040178838A1 (en) * | 2003-03-13 | 2004-09-16 | International Business Machines Corporation | Variable pulse width and pulse separation clock generator |
US6956422B2 (en) * | 2003-03-17 | 2005-10-18 | Indiana University Research And Technology Corporation | Generation and measurement of timing delays by digital phase error compensation |
US20060038598A1 (en) * | 2003-03-17 | 2006-02-23 | Reilly James P | Generation and measurement of timing delays by digital phase error compensation |
US6864833B2 (en) * | 2003-04-11 | 2005-03-08 | Siemens Milltronics Process Instruments, Inc. | Time-base generator with self-compensating control loop |
US6836239B2 (en) * | 2003-05-15 | 2004-12-28 | Aviation Communication & Surveillance Systems, Llc | Method and system for displaying text on a radial display |
US7026979B2 (en) * | 2003-07-03 | 2006-04-11 | Hrl Labortories, Llc | Method and apparatus for joint kinematic and feature tracking using probabilistic argumentation |
US7042385B1 (en) * | 2003-09-16 | 2006-05-09 | Niitek, Inc. | Non-intrusive inspection impulse radar antenna |
US7161531B1 (en) * | 2003-11-04 | 2007-01-09 | Northrop Grumman Corporation | High performance radar display |
US7203600B2 (en) * | 2003-12-19 | 2007-04-10 | Westerngeco, L.L.C. | Geophysical data processing |
US7053814B2 (en) * | 2004-01-15 | 2006-05-30 | Hrl Laboratories, Llc | RF lightwave coding system for radar pulse compression |
US20050200549A1 (en) * | 2004-03-15 | 2005-09-15 | Realtronics Corporation | Optimal Tapered Band Positioning to Mitigate Flare-End Ringing of Broadband Antennas |
US20050237260A1 (en) * | 2004-04-23 | 2005-10-27 | Centurion Wireless Technologies, Inc. | Microstrip Antenna |
US20060038599A1 (en) * | 2004-08-20 | 2006-02-23 | Avants Bradley S | Systems and methods for implementing delay line circuitry |
US7157952B2 (en) * | 2004-08-20 | 2007-01-02 | L-3 Integrated Systems Company | Systems and methods for implementing delay line circuitry |
US20060087471A1 (en) * | 2004-10-12 | 2006-04-27 | Hintz Kenneth J | Syntactic landmine detector |
US20060132210A1 (en) * | 2004-12-22 | 2006-06-22 | Lsi Logic Corporation | Wide-range programmable delay line |
US20080001808A1 (en) * | 2004-12-30 | 2008-01-03 | Passarelli Richard E Jr | System and method for processing data in weather radar |
US20060203613A1 (en) * | 2005-02-18 | 2006-09-14 | Leon Thomsen | System and method for using time-distance characteristics in acquisition, processing, and imaging of t-CSEM data |
US20060256025A1 (en) * | 2005-05-13 | 2006-11-16 | Realtronics Corporation | Machine Producible Directive Closed-Loop Impulse Antenna |
US20070080864A1 (en) * | 2005-10-11 | 2007-04-12 | M/A-Com, Inc. | Broadband proximity-coupled cavity backed patch antenna |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100329463A1 (en) * | 2009-06-24 | 2010-12-30 | Cisco Technology, Inc. | Group key management for mobile ad-hoc networks |
US8452014B2 (en) * | 2009-06-24 | 2013-05-28 | Cisco Technology, Inc. | Group key management for mobile ad-hoc networks |
US20140002292A1 (en) * | 2012-06-28 | 2014-01-02 | Tektronix, Inc. | Almost real-time sampling for ground-penetrating radar impulses |
US9581690B2 (en) * | 2012-06-28 | 2017-02-28 | Tektronix, Inc. | Almost real-time sampling for ground-penetrating radar impulses |
CN104755956A (en) * | 2012-11-02 | 2015-07-01 | 艾普斯科有限公司 | Method and apparatus for inspection of cooling towers |
CN109283523A (en) * | 2018-08-01 | 2019-01-29 | 西安交通大学 | A kind of geological radar B-scan data processing method |
CN112232392A (en) * | 2020-09-29 | 2021-01-15 | 深圳安德空间技术有限公司 | Data interpretation and identification method for three-dimensional ground penetrating radar |
CN116539643A (en) * | 2023-03-16 | 2023-08-04 | 南京京烁雷达科技有限公司 | Method and system for identifying coal rock data by using radar |
Also Published As
Publication number | Publication date |
---|---|
WO2009032467A1 (en) | 2009-03-12 |
US7675454B2 (en) | 2010-03-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US7675454B2 (en) | System, method, and computer program product providing three-dimensional visualization of ground penetrating radar data | |
CN105488770B (en) | A kind of airborne laser radar point cloud filtering method of object-oriented | |
EP0670049B1 (en) | Multi-dimensional signal processing and display | |
CN102508219B (en) | Turbulent current target detection method of wind profiler radar | |
US6943724B1 (en) | Identification and tracking of moving objects in detected synthetic aperture imagery | |
CN106772261A (en) | Radar detection signal multidimensional characteristic clusters visual display method | |
Gonzalez-Huici et al. | A combined strategy for landmine detection and identification using synthetic GPR responses | |
CN102201125A (en) | Method for visualizing three-dimensional imaging sonar data | |
CN104730148A (en) | Metal material inner inclusion three-dimensional reconstruction method based on ultrasonic testing technology | |
LeFeuvre et al. | Acoustic species identification in the Northwest Atlantic using digital image processing | |
CN113570005A (en) | Long-distance ship type identification method based on airborne photon radar | |
JP2009236720A (en) | Moving target detector | |
CN103954962B (en) | A kind of ISAR imaging pulse algorithm for estimating based on compressed sensing | |
Fakiris et al. | Quantification of regions of interest in swath sonar backscatter images using grey-level and shape geometry descriptors: The TargAn software | |
He et al. | High-accuracy acoustic sediment classification using sub-bottom profile data | |
Chen et al. | A ground-penetrating radar object detection method based on deep learning | |
KR101813357B1 (en) | Clutter elimination Method and Apparatus for applying variable data map and variable statistic characteristics information for identifying underwater fixed target and clutter | |
RU2736188C9 (en) | Hydroacoustic information displaying method | |
CN107657616A (en) | A kind of high score Remote Sensing Image Segmentation towards geographic object | |
US20090254302A1 (en) | Generating Images According To Points Of Intersection For Integer Multiples Of A Sample-Time Distance | |
Balabanian et al. | Sonar Explorer: A New Tool for Visualization of Fish Schools from 3D Sonar Data. | |
Nguyen | Seafloor classification with a multi-swath multi-beam echo sounder | |
Reinhardt et al. | Measuring human assessed complexity in synthetic aperture sonar imagery using the Elo rating system | |
Li et al. | HFSWR Clutter Recognition Based on Attention DCNN | |
Beech et al. | Quantifying the Impacts of Acoustic Target Detections Using a Range-Independent Model Versus Range-Dependent Model |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
AS | Assignment |
Owner name: NIITEK, INC.,VIRGINIA Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNORS:LAVEDAS, STEVEN;PATRIKAR, AJAY;HIBBARD, MARK;SIGNING DATES FROM 20071023 TO 20071102;REEL/FRAME:020152/0596 |
|
STCF | Information on status: patent grant |
Free format text: PATENTED CASE |
|
FPAY | Fee payment |
Year of fee payment: 4 |
|
SULP | Surcharge for late payment | ||
MAFP | Maintenance fee payment |
Free format text: PAYMENT OF MAINTENANCE FEE, 8TH YEAR, LARGE ENTITY (ORIGINAL EVENT CODE: M1552) Year of fee payment: 8 |
|
AS | Assignment |
Owner name: CHEMRING SENSORS AND ELECTRONIC SYSTEMS, INC., VIRGINIA Free format text: CHANGE OF NAME;ASSIGNOR:NON-INTRUSIVE INSPECTION TECHNOLOGY, INC.;REEL/FRAME:051610/0367 Effective date: 20161028 |
|
FEPP | Fee payment procedure |
Free format text: MAINTENANCE FEE REMINDER MAILED (ORIGINAL EVENT CODE: REM.); ENTITY STATUS OF PATENT OWNER: LARGE ENTITY |
|
LAPS | Lapse for failure to pay maintenance fees |
Free format text: PATENT EXPIRED FOR FAILURE TO PAY MAINTENANCE FEES (ORIGINAL EVENT CODE: EXP.); ENTITY STATUS OF PATENT OWNER: LARGE ENTITY |
|
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: 20220309 |