CN104698503A - Radar data processing method - Google Patents
Radar data processing method Download PDFInfo
- Publication number
- CN104698503A CN104698503A CN201510155943.6A CN201510155943A CN104698503A CN 104698503 A CN104698503 A CN 104698503A CN 201510155943 A CN201510155943 A CN 201510155943A CN 104698503 A CN104698503 A CN 104698503A
- Authority
- CN
- China
- Prior art keywords
- data
- radar
- disposal route
- digital processing
- route according
- 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.)
- Pending
Links
Landscapes
- Radar Systems Or Details Thereof (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention provides a radar data processing method. The radar data processing method comprises the steps of original data collection, data editing, data processing, image processing and data displaying and interpretation. The proper processing method and process are provided, various kinds of interference noise can be suppressed, effective signals can be enhanced, the signal to noise ratio of data can be improved, a high-resolution radar image section can be obtained, and the interpretation accuracy of a detection result can be improved.
Description
Technical field
The present invention relates to field of radar, relate to a kind of radar data disposal route specifically.
Background technology
The 1980s, ground penetrating radar is introduced in inspection of railway subgrade, has carried out rapidly relevant test and study, and achieve certain progress in countries in the world.Although along with development constantly updated by ground penetrating radar hardware, its sample rate, in good time display technique, investigation depth and precision aspect have and improve to some extent.But nearly all ground penetrating radar corollary system software can only provide data acquisition and simple data processing function, this had both limited the usable range (can only be confined to professional) of ground penetrating radar on the one hand, and the performance of its end result also cannot allow on-the-spot workmen understand easily on the other hand.
According to the difference of working method, ground penetrating radar is divided into ground ground penetrating radar and borehole radar.In borehole radar passing hole, antenna in depth detects down in inspecting hole, closer descends objective body, has larger investigative range.Because underground medium is equivalent to a complicated wave filter, it absorbs to some extent electromagnetic wave and reflects, and the unevenness etc. of itself, when making radar pulse get back to receiving antenna, wave amplitude reduces, and waveform also becomes and has a greater change with original transmitted waveform.In addition, various interference in various degree and random noise, also distorted measured data.
Summary of the invention
For above-mentioned the deficiencies in the prior art, the invention provides a kind of radar data disposal route, improve the precision of radar detection result, simultaneously for final explanation provides clear and legible radar detection image.
In order to achieve the above object, the invention provides a kind of radar data disposal route, the method comprises the steps:
Step one: raw data acquisition, gathers the data of radar monitoring;
Step 2: data edition, edits according to speed of detection, geological state and the surrounding environment data to record;
Step 3: digital processing, processes appearance skew composition, weak signal in data, and processes the medium filtering of very low frequency (VLF) distortion filtering and Space-time Domain that goes of data;
Step 4: image procossing, improves picture quality, carries out enhancing process to image information;
Step 5: data display and explanation, makes an explanation in conjunction with geology, environment to the data of image display.
Described data edition comprises each survey line data consistent, data merge, insert or delete road, create equidistant road and data reset.
In described step 3, digital processing comprises zero correction, and zero correction is used for the initial point position of track datas all in radar cross-section to move to time zero.
In described step 3, digital processing comprises the gain to radar wave attenuation signal.The mode of gain can be temporal gain SEC, energy attenuation gain, linear gain and exponential gain.
Preferably, described digital processing also comprises and processing bandpass filtering, mean filter, extraction mean-trace.
Described digital processing also comprises for eliminating earth filtering deconvolution.Preferably, deconvolution can be Least square deconvolution, predictive deconvolution, minimum entropy deconvolution.
Described digital processing also comprises for making up the earth absorption coefficient difference and damnous frequency spectrum compensation of decaying.
Preferably, in described digital processing, offset correction based on ray theory and kirchhoff wave equation migration method are comprised to the mode shifting into divisional processing.
Beneficial effect: the invention provides suitable disposal route and flow process can suppress various interference noise, strengthens useful signal, improves data signal to noise ratio (S/N ratio), and obtains high-resolution radar image section, improves the Explanation Accuracy of result of detection.
Embodiment
Below by description, the specific embodiment of the present invention is as the effect of the mutual alignment between the shape of involved each component, structure, each several part and annexation, each several part and principle of work, manufacturing process and operation using method etc., be described in further detail, have more complete, accurate and deep understanding to help those skilled in the art to inventive concept of the present invention, technical scheme.
A kind of radar data disposal route, the method comprises the steps:
Step one: raw data acquisition, gathers the data of radar monitoring;
Step 2: data edition, edits according to speed of detection, geological state and the surrounding environment data to record;
Step 3: digital processing, processes appearance skew composition, weak signal in data, and processes the medium filtering of very low frequency (VLF) distortion filtering and Space-time Domain that goes of data;
Step 4: image procossing, improves picture quality, carries out enhancing process to image information;
Step 5: data display and explanation, makes an explanation in conjunction with geology, environment to the data of image display.
Described data edition comprises each survey line data consistent, data merge, insert or delete road, create equidistant road and data reset.In described step 3, digital processing comprises zero correction, and zero correction is used for the initial point position of track datas all in radar cross-section to move to time zero.In described step 3, digital processing comprises the gain to radar wave attenuation signal.The mode of gain can be temporal gain SEC, energy attenuation gain, linear gain and exponential gain.Preferably, described digital processing also comprises and processing bandpass filtering, mean filter, extraction mean-trace.Described digital processing also comprises for eliminating earth filtering deconvolution, and compression radar ripple forms burst pulse form time delay, improve radar image cut open and resolution.Preferably, deconvolution can be Least square deconvolution, predictive deconvolution, minimum entropy deconvolution.Described digital processing also comprises for making up the earth absorption coefficient difference and damnous frequency spectrum compensation of decaying.Compensate by artificial mode, widen the spectral range of signal, improve the resolution of radar cross-section.Preferably, in described digital processing, offset correction based on ray theory and kirchhoff wave equation migration method are comprised to the mode shifting into divisional processing.The effect of skew moves on to home position by gathering each reflected signal obtained, and obtains the real image of search coverage.
Concrete processing mode: the first rank of data processing are data editions.The situations such as maloperation, omission or redundant data are had unavoidably in view of in raw data, editing and processing is exactly for these problems, data are reorganized and revises (merged file, the change of data file or background information, the dipole inversion of data track, changes the location parameter in road, to delete in data some bad track or inserts some empty roads etc.), after completing above-mentioned editor, usually to carry out reset (Onset point alignment) and the topographic correction of data.If signal amplitude value changes greatly on survey line section, the normalization of signal amplitude also should be carried out.During field detection, conveniently image data, the detection direction of different survey line is just in time contrary, and in the radar cross-section causing same objective body to be formed at different survey line, position is inconsistent, for ease of the comparative analysis of radar profile, need the data consistent of each survey line.When equidistant data acquisition, uneven due to antenna translational speed, each track data collected and physical location there are differences, and when especially survey line is very long, this phenomenon is inevitable.The method of process is the trace mark according to distance during collection in worksite, and the track data to whole cuing open carries out equidistant distribution process.Bandpass filtering is combined with low-pass filtering and high-pass filtering, and the signal made way in low-frequency cut-off frequency and this frequency bandwidth range of high-frequency cut-off frequency reaches signal to be passed through.Each sampling point amplitude sampling point amplitude average each in window time this, selecting in timing window, replaces, carries out data smoothing by mean filter.This filtering method can better remove high frequency noise extremely strong in data, and effect is similar to low-pass filtering, and data processing is determined than low-pass filtering.Extract mean-trace and extraction is carried out on average to each road chosen, then replace selected time range with sought radar data by this average, the effect of this filtering is the energy that suppression level is consistent, is particularly suitable for the horizontal signal removing the initial direct wave formation arrived.
The radar cross-section image that radar record signal is formed, can calculate detection range according to the time of reception of reflected signal, therefore whether the time of reception of reflected signal accurately determines the accuracy of detection range.For determining reflection interval, the starting point of reflected signal time of reception first to be determined.The effect of zero correction is exactly that the initial point position of track datas all in radar cross-section is moved on to time zero, guarantees the accurate identification of objective body distance.The general standard that time zero is determined is, the crest that in tracer signal, first direct wave amplitude is larger or trough.
Temporal gain mode is SEC (spherical wave exponential backoff).It carries out signal amplitude compensation according to electromagnetic wave in the attenuation characteristic of underground propagation, in use, user need according to underground medium situation, select suitable attenuation coefficient and gain maximum, wherein attenuation coefficient can be surveyed or table look-up acquisition, gain maximum then depending on signal to noise ratio (S/N ratio) and needs, generally value in 50 ~ 2000.
The method of graphics process is determined on a case-by-case basis, and usually has manual gain amplification, amplitude analysis and lineups pickup etc.Pickup lineups shows with synchronous fading and match, and energy reduced data section, gives prominence to some key character.Certainly, key character mentioned here depends on the processing procedure of data and the understanding to geological condition, and the physical difference of processing target body.
Generally, the net result of data processing adopts the display translation modes such as colored charting, gray scale charting or bargraphs.Wherein, colored charting can select different colour codes to carry out the response of outstanding object as required, has and shows image effect preferably.Gray scale charting and bargraphs are black white image, also commonly use in the hard copy of data image.The geologic interpretation of data image comprises many-sided knowledge and experiences such as high frequency technique, geography and geology, engineering and humanity.When image interpretation is evaluated, except the FORWARD AND INVERSE PROBLEMS feature should grasping various model, also should carry out interpretation of anomaly in conjunction with concrete geology, ambient conditions.
Claims (10)
1. a radar data disposal route, is characterized in that, the method comprises the steps:
Step one: raw data acquisition, gathers the data of radar monitoring;
Step 2: data edition, edits according to speed of detection, geological state and the surrounding environment data to record;
Step 3: digital processing, processes appearance skew composition, weak signal in data, and processes the medium filtering of very low frequency (VLF) distortion filtering and Space-time Domain that goes of data;
Step 4: image procossing, improves picture quality, carries out enhancing process to image information;
Step 5: data display and explanation, makes an explanation in conjunction with geology, environment to the data of image display.
2. radar data disposal route according to claim 1, is characterized in that: described data edition comprises each survey line data consistent, data merge, insert or delete road, create equidistant road and data reset.
3. radar data disposal route according to claim 1, is characterized in that: in described step 3, digital processing comprises zero correction, and zero correction is used for the initial point position of track datas all in radar cross-section to move to time zero.
4. radar data disposal route according to claim 1, is characterized in that: in described step 3, digital processing comprises the gain to radar wave attenuation signal.
5. radar data disposal route according to claim 4, is characterized in that: described gain comprises temporal gain SEC, energy attenuation gain, linear gain and exponential gain.
6. radar data disposal route according to claim 1, is characterized in that: described digital processing also comprises and processing bandpass filtering, mean filter, extraction mean-trace.
7. radar data disposal route according to claim 1, is characterized in that: described digital processing also comprises for eliminating earth filtering deconvolution.
8. radar data disposal route according to claim 7, is characterized in that: described deconvolution comprises Least square deconvolution, predictive deconvolution, minimum entropy deconvolution.
9. radar data disposal route according to claim 1, is characterized in that: described digital processing also comprises for making up the earth absorption coefficient difference and damnous frequency spectrum compensation of decaying.
10. radar data disposal route according to claim 1, is characterized in that: comprise offset correction based on ray theory and kirchhoff wave equation migration method to the mode shifting into divisional processing in described digital processing.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510155943.6A CN104698503A (en) | 2015-04-02 | 2015-04-02 | Radar data processing method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510155943.6A CN104698503A (en) | 2015-04-02 | 2015-04-02 | Radar data processing method |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN104698503A true CN104698503A (en) | 2015-06-10 |
Family
ID=53345821
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510155943.6A Pending CN104698503A (en) | 2015-04-02 | 2015-04-02 | Radar data processing method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104698503A (en) |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105975990A (en) * | 2016-05-17 | 2016-09-28 | 华北水利水电大学 | Reclaimed soil profile fine expression method |
| CN109031431A (en) * | 2018-08-10 | 2018-12-18 | 中国铁道科学研究院集团有限公司铁道建筑研究所 | A kind of data processing method and system for Coherent Noise in GPR Record |
| CN112285790A (en) * | 2020-03-30 | 2021-01-29 | 中国科学院地质与地球物理研究所 | Method and device for determining attenuation coefficient of electromagnetic wave field and medium |
| CN112666554A (en) * | 2020-12-17 | 2021-04-16 | 江苏中路工程技术研究院有限公司 | Method for identifying radar amplitude characteristic crack width of asphalt pavement |
| CN112684440A (en) * | 2020-12-17 | 2021-04-20 | 江苏中路工程技术研究院有限公司 | Crack evaluation method based on three-dimensional ground penetrating radar detection |
| CN115373035A (en) * | 2022-08-15 | 2022-11-22 | 江苏南京地质工程勘察院 | Frequency extension method and system for low-frequency signal of ground penetrating radar |
| CN116360001A (en) * | 2023-04-04 | 2023-06-30 | 中国地质科学院地球物理地球化学勘查研究所 | Method for eliminating low-frequency unshielded ground penetrating radar multiple signal interference in frozen soil region |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101566687A (en) * | 2008-04-23 | 2009-10-28 | 中国矿业大学(北京) | Numerical control acquisition system for geological radar |
| JP2011185834A (en) * | 2010-03-10 | 2011-09-22 | Nikko Tankai Kk | Sar data processing method and sar data processing system |
| WO2014106755A1 (en) * | 2013-01-07 | 2014-07-10 | Bae Systems Plc | Image processing |
-
2015
- 2015-04-02 CN CN201510155943.6A patent/CN104698503A/en active Pending
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101566687A (en) * | 2008-04-23 | 2009-10-28 | 中国矿业大学(北京) | Numerical control acquisition system for geological radar |
| JP2011185834A (en) * | 2010-03-10 | 2011-09-22 | Nikko Tankai Kk | Sar data processing method and sar data processing system |
| WO2014106755A1 (en) * | 2013-01-07 | 2014-07-10 | Bae Systems Plc | Image processing |
Non-Patent Citations (2)
| Title |
|---|
| 王帮兵: "《地球物理实验教程》", 31 May 2014 * |
| 邵景柱: "《复杂条件下的冻结并筒施工成套技术》", 30 November 2010 * |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105975990A (en) * | 2016-05-17 | 2016-09-28 | 华北水利水电大学 | Reclaimed soil profile fine expression method |
| CN105975990B (en) * | 2016-05-17 | 2019-05-28 | 华北水利水电大学 | A kind of fine expression of Reclaimed Soil section |
| CN109031431A (en) * | 2018-08-10 | 2018-12-18 | 中国铁道科学研究院集团有限公司铁道建筑研究所 | A kind of data processing method and system for Coherent Noise in GPR Record |
| CN112285790A (en) * | 2020-03-30 | 2021-01-29 | 中国科学院地质与地球物理研究所 | Method and device for determining attenuation coefficient of electromagnetic wave field and medium |
| CN112666554A (en) * | 2020-12-17 | 2021-04-16 | 江苏中路工程技术研究院有限公司 | Method for identifying radar amplitude characteristic crack width of asphalt pavement |
| CN112684440A (en) * | 2020-12-17 | 2021-04-20 | 江苏中路工程技术研究院有限公司 | Crack evaluation method based on three-dimensional ground penetrating radar detection |
| CN115373035A (en) * | 2022-08-15 | 2022-11-22 | 江苏南京地质工程勘察院 | Frequency extension method and system for low-frequency signal of ground penetrating radar |
| CN115373035B (en) * | 2022-08-15 | 2024-08-23 | 江苏南京地质工程勘察院 | Frequency extension method and system for ground penetrating radar low-frequency signal |
| CN116360001A (en) * | 2023-04-04 | 2023-06-30 | 中国地质科学院地球物理地球化学勘查研究所 | Method for eliminating low-frequency unshielded ground penetrating radar multiple signal interference in frozen soil region |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104698503A (en) | Radar data processing method | |
| Cui et al. | Estimating tree-root biomass in different depths using ground-penetrating radar: Evidence from a controlled experiment | |
| Cassidy et al. | Ground penetrating radar data processing, modelling and analysis | |
| Walter et al. | Spatial prediction of topsoil salinity in the Chelif Valley, Algeria, using local ordinary kriging with local variograms versus whole-area variogram | |
| US7492664B2 (en) | Method for processing acoustic reflections in array data to image near-borehole geological structure | |
| CN102937721B (en) | Limited frequency tomography method for utilizing preliminary wave travel time | |
| CN102680575B (en) | A kind of impact mapping method of Complicate soil medium and system | |
| Leckebusch | Two‐and three‐dimensional ground‐penetrating radar surveys across a medieval choir: a case study in archaeology | |
| Wei et al. | Autofocusing techniques for GPR data from RC bridge decks | |
| CN101813644B (en) | Method for measuring biomass of plant underground roots by using ground penetrating radar | |
| CN105182416A (en) | Seismic inversion method and device based on frequency division data | |
| CN107479042B (en) | Estimation method for spatial water storage capacity of surface karst zone | |
| CN104516018A (en) | Porosity inversion method under lithological constraint in geophysical exploration | |
| CN107678064B (en) | Real-time extraction method for sound wave time difference | |
| Aboudourib et al. | A processing framework for tree-root reconstruction using ground-penetrating radar under heterogeneous soil conditions | |
| CN111708084A (en) | High-precision karst identification method based on seismic mapping data | |
| CN109765626A (en) | A kind of lunar exploration radar data processing method based on the offset of least square kirchhoff | |
| CN106019377B (en) | A kind of two-dimensional seismic survey noise remove method based on time-space domain frequency reducing model | |
| CN111142165A (en) | Method for acquiring water level information of aquifer by using ground penetrating radar | |
| Liu et al. | Detection of karst cavity beneath cast-in-place pile using the instantaneous phase difference of two receiver recordings | |
| Chen et al. | Improving the Precision of Surface Seismic Data Processing by Walkaway VSP | |
| Zhang et al. | Data processing of a wide-azimuth, broadband, high-density 3D seismic survey using a low-frequency vibroseis: a case study from Northeast China | |
| CN111880223B (en) | Attribute analysis method for well-side sound wave remote detection imaging fracture-cave carving | |
| Cowan et al. | Screening kimberlite magnetic anomalies in magnetically active areas | |
| Tassis et al. | Comparison between Sensors and Software and Malå GPR equipment based on test measurements at Bøaøyna, Stryn Municipality, Norway |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20150610 |
|
| RJ01 | Rejection of invention patent application after publication |