CN111007572A - Automatic identification method, device and system for road underground cavity - Google Patents
Automatic identification method, device and system for road underground cavity Download PDFInfo
- Publication number
- CN111007572A CN111007572A CN201911163987.8A CN201911163987A CN111007572A CN 111007572 A CN111007572 A CN 111007572A CN 201911163987 A CN201911163987 A CN 201911163987A CN 111007572 A CN111007572 A CN 111007572A
- Authority
- CN
- China
- Prior art keywords
- reflected wave
- area
- reflection coefficient
- amplitude
- value
- 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
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
- G01V3/12—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with electromagnetic waves
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Radar, Positioning & Navigation (AREA)
- General Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geophysics (AREA)
- Environmental & Geological Engineering (AREA)
- Computer Networks & Wireless Communication (AREA)
- Radar Systems Or Details Thereof (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention provides an automatic identification method, a device and a system for road underground cavities, wherein the method comprises the following steps: acquiring sensing data of the ground penetrating radar of each area on the road, wherein the sensing data is used for representing the amplitude of a reflected wave based on the ground penetrating radar; calculating the reflected wave energy value of each area according to the sensing data of the ground penetrating radar, and determining a first area of which the reflected wave energy value is greater than a first preset threshold value; calculating a reflection coefficient of a reflected wave in the first region and a Fourier spectrum amplitude of the reflected wave according to the sensing data of the first region, and further calculating a cavity response value of the reflected wave in the first region according to the reflection coefficient of the reflected wave in the first region and the Fourier spectrum amplitude of the reflected wave; and judging the position corresponding to the first area with the cavity response value larger than the second preset threshold value as an underground cavity. The method can improve the accuracy and reliability of the detection of the underground cavities of the road.
Description
Technical Field
The invention relates to the field of road detection, in particular to an automatic identification method, device and system for underground cavities of roads.
Background
Due to vehicle vibration, road surface water seepage, underground pipeline water seepage and other reasons, the urban road frequently has the problems of cracking, deformation, settlement, collapse and the like. By means of a geophysical prospecting method, urban road detection is carried out regularly, and hidden danger of road collapse can be early warned in advance. The ground penetrating radar method is a road detection geophysical prospecting method which can simultaneously meet the requirements of rapidness, no damage and high resolution.
At present, radar road detection data are interpreted manually, and in the face of mass data generated by rapid road detection, the manual method is low in efficiency and cannot submit an interpretation result in time, interpretation accuracy rates of different personnel are greatly different, and reliability is low.
Disclosure of Invention
The invention aims to provide a method, a device and a system for automatically identifying underground road cavities, which aim to solve the problems of low accuracy and low reliability of the manual interpretation of radar road detection data in the prior art.
According to a first aspect of the invention, a method for automatically identifying a road underground cavity comprises the following steps: acquiring sensing data of the ground penetrating radar of each area on the road; the sensing data is used for representing the amplitude of a ground penetrating radar-based reflected wave; calculating the reflected wave energy value of each area according to the sensing data of the ground penetrating radar, and determining a first area of which the reflected wave energy value is greater than a first preset threshold value; calculating a reflection coefficient of a reflected wave in the first region and a Fourier spectrum amplitude of the reflected wave according to the sensing data of the first region, and further calculating a cavity response value of the reflected wave in the first region according to the reflection coefficient of the reflected wave in the first region and the Fourier spectrum amplitude of the reflected wave; and judging the position corresponding to the first area with the cavity response value larger than the second preset threshold value as an underground cavity.
Further, the step of calculating the cavity response value of the reflected wave in the first region according to the reflection coefficient of the reflected wave in the first region and the fourier spectrum amplitude of the reflected wave includes:
and calculating the void response value of the reflected wave in the first area according to the weighted summation of the reflection coefficient of the reflected wave in the first area and the Fourier spectrum amplitude of the reflected wave.
Further, before the step of calculating the hole response value of the reflected wave in the first region according to the weighted summation of the reflection coefficient of the reflected wave in the first region and the fourier spectrum amplitude of the reflected wave, the method comprises:
calculating the reflection coefficient of the reflected wave in each first area, and calculating the average value of the reflection coefficients of the reflected wave in each first area;
determining areas with the reflection coefficients larger than the average value in each first area as second areas, wherein the first areas comprise the second areas;
calculating a cavity response value of the reflected wave in the first area according to the weighted summation of the reflection coefficient of the reflected wave in the first area and the Fourier spectrum amplitude of the reflected wave, and determining the position corresponding to the first area with the cavity response value larger than a second preset threshold value as an underground cavity, wherein the step comprises the following steps:
and calculating the cavity response value of the reflected wave in the second area according to the weighted summation of the reflection coefficient of the reflected wave in the second area and the Fourier spectrum amplitude of the reflected wave, and judging the position corresponding to the second area with the cavity response value being greater than a second preset threshold value as an underground cavity.
Further, the step of acquiring sensing data of the ground penetrating radar is followed by: and filtering the sensing data of the acquired ground penetrating radar.
Further, the reflected wave energy value E is calculated as follows:wherein N is the number of sampling points, PiThe amplitude of the reflected wave of the ith sampling point is obtained;
the reflection coefficient R is calculated as follows:wherein, P0The peak value of the direct wave; the fourier spectrum amplitude F (ω) of the reflected wave is calculated as follows:
f (t) is the time domain signal of the amplitude of the reflected wave, ω is the amplitude signal frequency of the reflected wave, and t is the amplitude signal of the reflected waveSampling time;
the hole response value K of the reflected wave is calculated in the following way that K is α R + β F (omega), α is a preset reflection coefficient weight, and β is a preset Fourier spectrum amplitude weight.
According to a second aspect of the present invention, an apparatus for automatically identifying a road underground cavity, comprises: the system comprises an acquisition module, a detection module and a processing module, wherein the acquisition module is used for acquiring sensing data of the ground penetrating radar in each area on a road, and the sensing data is used for representing the amplitude of a reflected wave based on the ground penetrating radar;
the processing module is used for calculating the reflected wave energy value of each area according to the sensing data of the ground penetrating radar and determining a first area with the reflected wave energy value being larger than a first preset threshold value; calculating a reflection coefficient of a reflected wave in the first region and a Fourier spectrum amplitude of the reflected wave according to the sensing data of the first region, and further calculating a cavity response value of the reflected wave in the first region according to the reflection coefficient of the reflected wave in the first region and the Fourier spectrum amplitude of the reflected wave;
and the identification module is used for judging the position corresponding to the first area with the cavity response value larger than a second preset threshold value as an underground cavity.
Further, the processing module is further configured to calculate a void response value of the reflected wave in the first region according to a weighted summation of a reflection coefficient of the reflected wave in the first region and a fourier spectrum amplitude of the reflected wave.
Further, the processing module is further configured to calculate a reflection coefficient of the reflected wave in each first area, and calculate an average value of the reflection coefficients of the reflected waves in each first area; determining an area with a reflection coefficient larger than the average value in each first area as a second area; calculating a hole response value of the reflected wave in the second area according to the weighted summation of the reflection coefficient of the reflected wave in the second area and the Fourier spectrum amplitude of the reflected wave, and judging the position corresponding to the second area with the hole response value larger than a second preset threshold value as an underground hole; wherein the first region comprises a second region.
Further, the automatic identification device for the underground cavity of the road further comprises: the filtering module is used for filtering the sensing data of the acquired ground penetrating radar;
the reflected wave energy value E is calculated as follows:wherein N is the number of sampling points, PiThe amplitude of the reflected wave of the ith sampling point is obtained;
the reflection coefficient R is calculated as follows:wherein, P0The peak value of the direct wave; the fourier spectrum amplitude F (ω) of the reflected wave is calculated as follows:
f (t) is a time domain signal of the amplitude of the reflected wave, omega is the amplitude signal frequency of the reflected wave, and t is the amplitude signal sampling time of the reflected wave;
the hole response value K of the reflected wave is calculated in the following way that K is α R + β F (omega), α is a preset reflection coefficient weight, and β is a preset Fourier spectrum amplitude weight.
According to a third aspect of the invention, the system for automatically identifying the road underground cavity comprises the device for automatically identifying the road underground cavity and a ground penetrating radar, wherein the device for automatically identifying the road underground cavity is in communication connection with the ground penetrating radar.
According to the method, the device and the system for automatically identifying the road underground cavity, the first area with the energy accumulated value of the ground penetrating radar sensing data larger than the first preset threshold value is determined, the cavity response value of the first area is obtained through calculation according to the reflection coefficient of the first area and the Fourier spectrum amplitude, the position corresponding to the first area with the cavity response value larger than the second preset threshold value is determined as the underground cavity, and the accuracy and the reliability of radar road detection data interpretation are improved based on the analysis of the difference of the physical properties of the cavity and the surrounding medium at the interface of the underground cavity.
Other characteristic features and advantages of the invention will become apparent from the following description of exemplary embodiments, which is to be read in connection with the accompanying drawings.
Drawings
The accompanying drawings, which are incorporated in and constitute a part of the specification, illustrate embodiments of the invention and together with the description, serve to explain the principles of the invention. In the drawings, like reference numerals are used to indicate like elements. The drawings in the following description are directed to some, but not all embodiments of the invention. For a person skilled in the art, other figures can be derived from these figures without inventive effort.
FIG. 1 is a flow chart of an embodiment of a method for automatically identifying road underground cavities according to the present invention;
FIG. 2 is a flow chart of another embodiment of the method for automatically identifying underground cavities in roads of the present invention;
FIG. 3 is a schematic diagram of features of the holes identified by the automatic identification method shown in FIG. 2 in a radar image;
FIG. 4 is a schematic diagram of a Fourier magnitude spectrum of a region shown by a dashed box at a hole position in FIG. 3, wherein a peak between an abscissa 200 and an abscissa 400 is a characteristic of the hole in the Fourier spectrum;
the dotted line in FIG. 5 marks the features of the homogeneous medium in the radar image;
FIG. 6 is a Fourier magnitude spectrum at the dashed line of the uniform medium in FIG. 5, wherein the peaks between the abscissa 200 and the abscissa 400 are the characteristic of the non-voids in the Fourier spectrum;
FIG. 7 is a block diagram of an embodiment of an automatic road underground cavity recognition device according to the present invention;
FIG. 8 is a block diagram of an embodiment of an automatic identification system for underground cavities in a road according to the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the embodiments of the present invention clearer, the technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are some, but not all, embodiments of the present invention. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention. It should be noted that the embodiments and features of the embodiments in the present application may be arbitrarily combined with each other without conflict.
The research of the invention finds that: the difference of the physical properties of the cavity and the surrounding medium at the interface of the underground cavity is huge, strong reflection is generated when electromagnetic waves are transmitted, the strong reflection is mainly based on high-frequency components, and a high-frequency peak value can be formed in a frequency spectrum. When the electromagnetic wave is transmitted in a uniform medium, strong reflection cannot be generated, and corresponding high-frequency response cannot exist in a frequency spectrum. The difference between the electromagnetic physical properties of the cavity and the surrounding medium is huge, the reflection coefficient obtained through calculation is larger than that of a normal area, and meanwhile the underground cavity of the road can be accurately identified by combining Fourier frequency spectrum. The following detailed description is made in conjunction with the embodiments illustrated in the various figures.
As shown in FIG. 1, the invention relates to an automatic identification method of road underground cavities, which comprises the following steps:
step 101: acquiring sensing data of the ground penetrating radar of each area on the road, wherein the sensing data is used for representing the amplitude of a reflected wave based on the ground penetrating radar;
step 102: calculating the reflected wave energy value of each area according to the sensing data of the ground penetrating radar, and determining a first area of which the reflected wave energy value is greater than a first preset threshold value; calculating a reflection coefficient of a reflected wave in the first region and a Fourier spectrum amplitude of the reflected wave according to the sensing data of the first region, and further calculating a cavity response value of the reflected wave in the first region according to the reflection coefficient of the reflected wave in the first region and the Fourier spectrum amplitude of the reflected wave;
step 103: and judging the position corresponding to the first area with the cavity response value larger than the second preset threshold value as an underground cavity.
According to the embodiment of the automatic identification method of the road underground cavity, the first area with the energy accumulated value of the ground penetrating radar sensing data larger than the first preset threshold value is determined, the cavity response value of the first area is obtained through calculation according to the reflection coefficient of the first area and the Fourier spectrum amplitude, the position corresponding to the first area with the cavity response value larger than the second preset threshold value is determined as the underground cavity, and the accuracy and the reliability of radar road detection data interpretation are improved based on the analysis of the difference of the physical properties of the cavity and the surrounding medium at the interface of the underground cavity.
As shown in fig. 2, the present invention provides an automatic identification method for underground road cavities, which is a preferred implementation of the embodiment of the method shown in fig. 1, and the explanation of the embodiment shown in fig. 1 can be applied to this embodiment, and the automatic identification method of this embodiment includes:
step 201: carrying out data loading to acquire sensing data of the ground penetrating radar;
step 202: filtering the ground penetrating radar data to remove interference clutter; the specific filtering manner can be various, such as mean filtering, FK filtering (Frequency Wave-number domain filtering), etc.;
step 203: calculating an energy accumulated value for the filtered data; the energy accumulation value can be calculated in various ways, such as summing the squares of the amplitude values of all points of each data track; the reflected wave energy value E is calculated as follows:wherein N is the number of sampling points, PiThe amplitude of the reflected wave of the ith sampling point is obtained;
step 204: determining a region with an energy value smaller than or equal to a first preset threshold as a non-cavity region; the size of the first preset threshold is set according to actual needs, and is not limited herein;
step 205: determining an area with an energy value larger than a first preset threshold value as a first area, and calculating a reflection coefficient; calculating the average value of the reflection coefficients of all the first areas; the reflection coefficient refers to the ratio of the amplitude value domain direct wave amplitude value of each point of each channel of data; such asThe reflection coefficient R is calculated as follows:wherein, P0The peak value of the direct wave;
step 206: determining a first area with a reflection coefficient less than or equal to the average value of the reflection coefficients as a non-void area;
step 207: calculating a fourier spectrum of a region (as a second region, the first region includes the second region) in which the reflection coefficient is larger than the average value; for the area (namely a second area) with the reflection coefficient larger than the average value, carrying out weighted summation on the reflection coefficient and the Fourier spectrum amplitude to obtain a void response value; such as: the fourier spectrum magnitude F (ω) of the reflected wave is calculated as follows:
f (t) is a time domain signal of the amplitude of the reflected wave, omega is the amplitude signal frequency of the reflected wave, and t is the amplitude signal sampling time of the reflected wave;
the calculation method of the hole response value K of the reflected wave is as follows, wherein K is α R + β F (omega), α is a preset reflection coefficient weight, β is a preset Fourier spectrum amplitude weight, the specific values of α and β can be set according to actual needs, including but not limited to the fact that the weights of the two are both 1, and unequal weights can be selected according to actual experimental conditions.
Step 208: judging the position of the void response value less than or equal to a second preset threshold value as a non-underground void; the size of the second preset threshold is set according to actual needs, and is not limited herein;
step 209: and judging the position with the cavity response value larger than the second preset threshold value as an underground cavity.
The experimental results of the automatic identification method for the road underground cavity in the embodiment can be seen in fig. 3-6, and it can be seen that the accuracy and reliability of interpretation of radar road detection data can be effectively improved by the automatic identification method for the road underground cavity in the embodiment.
According to the embodiment of the automatic identification method of the road underground cavity, after the first area is determined, the reflection coefficient of each data is calculated for the whole radar data, the size of the reflection coefficient of each data is compared with the average value of the reflection coefficients, the second area with the reflection coefficient larger than the average value of the reflection coefficients is defined, FFT operation is conducted on the data of the second area, the Fourier frequency spectrum of each data is calculated, weighted summation calculation is conducted on the reflection coefficient and the Fourier frequency spectrum, the cavity response value is obtained, the position with the cavity response value larger than the threshold value can be determined as the underground cavity area, and the accuracy and reliability of radar road detection data interpretation can be effectively improved.
As shown in fig. 7, the present invention further provides an automatic identification device for underground road cavities, which is a corresponding device embodiment to the method embodiments shown in fig. 1 and 2, and the explanation of the embodiments shown in fig. 1 to 6 can be applied to this embodiment, and the automatic identification device includes:
an obtaining module 701, configured to obtain sensing data of a ground penetrating radar in each area on a road, where the sensing data is used to characterize an amplitude of a reflected wave based on the ground penetrating radar;
the processing module 702 is configured to calculate a reflected wave energy value of each region according to the sensing data of the ground penetrating radar, and determine a first region in which the reflected wave energy value is greater than a first preset threshold; calculating a reflection coefficient of a reflected wave in the first region and a Fourier spectrum amplitude of the reflected wave according to the sensing data of the first region, and further calculating a cavity response value of the reflected wave in the first region according to the reflection coefficient of the reflected wave in the first region and the Fourier spectrum amplitude of the reflected wave;
the identifying module 703 is configured to determine a position corresponding to the first area where the cavity response value is greater than the second preset threshold as an underground cavity.
Further, the processing module 702 is further configured to calculate a void response value of the first region according to a weighted summation of the reflection coefficient of the first region and the fourier spectrum amplitude.
Further, the processing module 702 is further configured to calculate a reflection coefficient of the reflected wave in each first area, and calculate an average value of the reflection coefficients of the reflected waves in each first area; determining an area with a reflection coefficient larger than the average value in each first area as a second area; calculating a hole response value of the reflected wave in the second area according to the weighted summation of the reflection coefficient of the reflected wave in the second area and the Fourier spectrum amplitude of the reflected wave, and judging the position corresponding to the second area with the hole response value larger than a second preset threshold value as an underground hole; wherein the first region comprises a second region.
Further, the automatic identification device for the underground road cavity may further include: and a filtering module 704, configured to perform filtering processing on the sensing data of the ground penetrating radar.
As shown in fig. 8, the present invention further provides an automatic identification system for underground road cavities, which comprises the automatic identification device shown in fig. 7. The explanation of the embodiment shown in fig. 1-7 can be applied to this embodiment, and the automatic identification system of this embodiment includes the automatic identification device for the underground cavity of the road and the ground penetrating radar, and the automatic identification device is connected to the ground penetrating radar in communication.
The automatic identification system for the road underground cavity has the corresponding technical effects of the automatic identification device for the road underground cavity, and the details are not repeated herein.
The above-described aspects may be implemented individually or in various combinations, and such variations are within the scope of the present invention.
Finally, it should be noted that: the above examples are only for illustrating the technical solutions of the present invention, and are not limited thereto. Although the present invention has been described in detail with reference to the foregoing embodiments, it will be understood by those of ordinary skill in the art that: the technical solutions described in the foregoing embodiments may still be modified, or some technical features may be equivalently replaced; and such modifications or substitutions do not depart from the spirit and scope of the corresponding technical solutions of the embodiments of the present invention.
Claims (10)
1. An automatic identification method for underground road cavities is characterized by comprising the following steps:
acquiring sensing data of the ground penetrating radar of each area on the road, wherein the sensing data is used for representing the amplitude of a reflected wave based on the ground penetrating radar;
calculating the reflected wave energy value of each area according to the sensing data of the ground penetrating radar, and determining a first area of which the reflected wave energy value is greater than a first preset threshold value;
calculating a reflection coefficient of a reflected wave in the first region and a Fourier spectrum amplitude of the reflected wave according to the sensing data of the first region, and further calculating a cavity response value of the reflected wave in the first region according to the reflection coefficient of the reflected wave in the first region and the Fourier spectrum amplitude of the reflected wave;
and judging the position corresponding to the first area with the cavity response value larger than the second preset threshold value as an underground cavity.
2. The method for automatically identifying an underground cavity of a road according to claim 1, wherein: calculating a hole response value of the reflected wave in the first area according to the reflection coefficient of the reflected wave in the first area and the Fourier spectrum amplitude of the reflected wave, wherein the step comprises the following steps:
and calculating the void response value of the reflected wave in the first area according to the weighted summation of the reflection coefficient of the reflected wave in the first area and the Fourier spectrum amplitude of the reflected wave.
3. The method for automatically identifying an underground cavity of a road according to claim 2, wherein: the method comprises the following steps of calculating the void response value of the reflected wave in the first area according to the weighted summation of the reflection coefficient of the reflected wave in the first area and the Fourier spectrum amplitude of the reflected wave, wherein the method comprises the following steps:
calculating the reflection coefficient of the reflected wave in each first area, and calculating the average value of the reflection coefficients of the reflected wave in each first area;
determining areas with the reflection coefficients larger than the average value in each first area as second areas, wherein the first areas comprise the second areas;
calculating a cavity response value of the reflected wave in the first area according to the weighted summation of the reflection coefficient of the reflected wave in the first area and the Fourier spectrum amplitude of the reflected wave, and determining the position corresponding to the first area with the cavity response value larger than a second preset threshold value as an underground cavity, wherein the step comprises the following steps:
and calculating the cavity response value of the reflected wave in the second area according to the weighted summation of the reflection coefficient of the reflected wave in the second area and the Fourier spectrum amplitude of the reflected wave, and judging the position corresponding to the second area with the cavity response value being greater than a second preset threshold value as an underground cavity.
4. The method of claim 3, wherein the step of obtaining sensed data of the ground penetrating radar of each area on the road is followed by:
and filtering the sensing data of the acquired ground penetrating radar.
5. The method for automatically identifying an underground cavity of a road according to any one of claims 1 to 4, wherein:
the reflected wave energy value E is calculated as follows:wherein N is the number of sampling points, PiThe amplitude of the reflected wave of the ith sampling point is obtained;
the reflection coefficient R is calculated as follows:wherein, P0The peak value of the direct wave; the fourier spectrum amplitude F (ω) of the reflected wave is calculated as follows:
f (t) is the time domain signal of the amplitude of the reflected wave, ω is the amplitude of the reflected waveThe signal frequency t is the amplitude signal sampling time of the reflected wave;
the hole response value K of the reflected wave is calculated in the following way that K is α R + β F (omega), α is a preset reflection coefficient weight, and β is a preset Fourier spectrum amplitude weight.
6. An automatic identification device of road underground cavity, characterized by, includes:
the system comprises an acquisition module, a detection module and a processing module, wherein the acquisition module is used for acquiring sensing data of the ground penetrating radar in each area on a road, and the sensing data is used for representing the amplitude of a reflected wave based on the ground penetrating radar;
the processing module is used for calculating the reflected wave energy value of each area according to the sensing data of the ground penetrating radar and determining a first area with the reflected wave energy value being larger than a first preset threshold value; calculating a reflection coefficient of a reflected wave in the first region and a Fourier spectrum amplitude of the reflected wave according to the sensing data of the first region, and further calculating a cavity response value of the reflected wave in the first region according to the reflection coefficient of the reflected wave in the first region and the Fourier spectrum amplitude of the reflected wave;
and the identification module is used for judging the position corresponding to the first area with the cavity response value larger than a second preset threshold value as an underground cavity.
7. The apparatus for automatically identifying an underground cavity of a road according to claim 6, wherein: the processing module is further configured to calculate a void response value of the reflected wave in the first region according to a weighted summation of a reflection coefficient of the reflected wave in the first region and a fourier spectrum amplitude of the reflected wave.
8. The apparatus according to claim 7, wherein the processing module is further configured to calculate a reflection coefficient of the reflected wave in each of the first areas, and calculate an average value of the reflection coefficients of the reflected waves in each of the first areas; determining an area with a reflection coefficient larger than the average value in each first area as a second area; calculating a hole response value of the reflected wave in the second area according to the weighted summation of the reflection coefficient of the reflected wave in the second area and the Fourier spectrum amplitude of the reflected wave, and judging the position corresponding to the second area with the hole response value larger than a second preset threshold value as an underground hole; wherein the first region comprises a second region.
9. The apparatus for automatically identifying an underground cavity of a road according to any one of claims 6 to 8, further comprising:
the filtering module is used for filtering the sensing data of the acquired ground penetrating radar;
the reflected wave energy value E is calculated as follows:wherein N is the number of sampling points, PiThe amplitude of the reflected wave of the ith sampling point is obtained;
the reflection coefficient R is calculated as follows:wherein, P0The peak value of the direct wave; the fourier spectrum amplitude F (ω) of the reflected wave is calculated as follows:
f (t) is a time domain signal of the amplitude of the reflected wave, omega is the amplitude signal frequency of the reflected wave, and t is the amplitude signal sampling time of the reflected wave;
the hole response value K of the reflected wave is calculated in the following way that K is α R + β F (omega), α is a preset reflection coefficient weight, and β is a preset Fourier spectrum amplitude weight.
10. An automatic identification system for underground road cavities, characterized in that it comprises an automatic identification device for underground road cavities according to any one of claims 6 to 9 and a ground penetrating radar, said automatic identification device being connected in communication with the ground penetrating radar.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911163987.8A CN111007572B (en) | 2019-11-22 | 2019-11-22 | Automatic identification method, device and system for road underground cavity |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911163987.8A CN111007572B (en) | 2019-11-22 | 2019-11-22 | Automatic identification method, device and system for road underground cavity |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111007572A true CN111007572A (en) | 2020-04-14 |
CN111007572B CN111007572B (en) | 2020-11-10 |
Family
ID=70113248
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911163987.8A Active CN111007572B (en) | 2019-11-22 | 2019-11-22 | Automatic identification method, device and system for road underground cavity |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111007572B (en) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114076943A (en) * | 2021-11-17 | 2022-02-22 | 广州市市政工程试验检测有限公司 | Rapid automatic comprehensive detection method for urban road underground cavity |
CN114137517A (en) * | 2022-02-07 | 2022-03-04 | 北京中科蓝图科技有限公司 | Ground penetrating detection method and device for road and ground penetrating radar device |
CN115267143A (en) * | 2022-09-28 | 2022-11-01 | 江苏筑升土木工程科技有限公司 | Road cavity defect detection system and detection method |
CN115793086A (en) * | 2023-02-07 | 2023-03-14 | 武汉新楚光电科技发展有限公司 | Optical cable laying environment underground cavity judgment method and system based on optical fiber sensing |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS6022681A (en) * | 1983-07-19 | 1985-02-05 | Japan Radio Co Ltd | Method and device for wave height radar observation |
JPH0587945A (en) * | 1991-09-27 | 1993-04-09 | Doro Hozen Gijutsu Center | Cavity inspection method for paved road |
JP5719075B1 (en) * | 2014-10-06 | 2015-05-13 | ジオ・サーチ株式会社 | Cavity thickness exploration method |
US20160202346A1 (en) * | 2013-06-15 | 2016-07-14 | Howard University | Using An MM-Principle to Enforce a Sparsity Constraint on Fast Image Data Estimation From Large Image Data Sets |
CN105891896A (en) * | 2016-04-25 | 2016-08-24 | 湖南科技大学 | Feature information recognition and analysis method for underground mined area |
CN109001728A (en) * | 2018-06-26 | 2018-12-14 | 石家庄铁道大学 | Method and device based on disease inside Ground Penetrating Radar detection armored concrete |
CN109521421A (en) * | 2018-01-27 | 2019-03-26 | 河南工业大学 | A kind of Ground Penetrating Radar thin layer object recognition and detection method |
CN109615692A (en) * | 2018-12-12 | 2019-04-12 | 中国矿业大学(北京) | A kind of underground piping surrounding objects interpolation modeling algorithm based on Ground Penetrating Radar |
CN110319916A (en) * | 2019-06-06 | 2019-10-11 | 浙江大学 | A kind of limited pond low frequency expanding method based on the measurement of waters interference of reflected wave |
-
2019
- 2019-11-22 CN CN201911163987.8A patent/CN111007572B/en active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS6022681A (en) * | 1983-07-19 | 1985-02-05 | Japan Radio Co Ltd | Method and device for wave height radar observation |
JPH0587945A (en) * | 1991-09-27 | 1993-04-09 | Doro Hozen Gijutsu Center | Cavity inspection method for paved road |
US20160202346A1 (en) * | 2013-06-15 | 2016-07-14 | Howard University | Using An MM-Principle to Enforce a Sparsity Constraint on Fast Image Data Estimation From Large Image Data Sets |
JP5719075B1 (en) * | 2014-10-06 | 2015-05-13 | ジオ・サーチ株式会社 | Cavity thickness exploration method |
CN105891896A (en) * | 2016-04-25 | 2016-08-24 | 湖南科技大学 | Feature information recognition and analysis method for underground mined area |
CN109521421A (en) * | 2018-01-27 | 2019-03-26 | 河南工业大学 | A kind of Ground Penetrating Radar thin layer object recognition and detection method |
CN109001728A (en) * | 2018-06-26 | 2018-12-14 | 石家庄铁道大学 | Method and device based on disease inside Ground Penetrating Radar detection armored concrete |
CN109615692A (en) * | 2018-12-12 | 2019-04-12 | 中国矿业大学(北京) | A kind of underground piping surrounding objects interpolation modeling algorithm based on Ground Penetrating Radar |
CN110319916A (en) * | 2019-06-06 | 2019-10-11 | 浙江大学 | A kind of limited pond low frequency expanding method based on the measurement of waters interference of reflected wave |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114076943A (en) * | 2021-11-17 | 2022-02-22 | 广州市市政工程试验检测有限公司 | Rapid automatic comprehensive detection method for urban road underground cavity |
CN114137517A (en) * | 2022-02-07 | 2022-03-04 | 北京中科蓝图科技有限公司 | Ground penetrating detection method and device for road and ground penetrating radar device |
CN115267143A (en) * | 2022-09-28 | 2022-11-01 | 江苏筑升土木工程科技有限公司 | Road cavity defect detection system and detection method |
CN115793086A (en) * | 2023-02-07 | 2023-03-14 | 武汉新楚光电科技发展有限公司 | Optical cable laying environment underground cavity judgment method and system based on optical fiber sensing |
CN115793086B (en) * | 2023-02-07 | 2023-06-06 | 武汉新楚光电科技发展有限公司 | Optical cable laying environment underground cavity judging method and system based on optical fiber sensing |
Also Published As
Publication number | Publication date |
---|---|
CN111007572B (en) | 2020-11-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111007572B (en) | Automatic identification method, device and system for road underground cavity | |
CN111476088B (en) | Asphalt pavement water damage identification model construction method, identification method and system | |
CN104020495B (en) | Automatic underground pipeline parameter recognizing method on basis of ground penetrating radar | |
CN110376584B (en) | Water supply pipeline leakage detection method based on ground penetrating radar image characteristic signal identification | |
KR101917374B1 (en) | Apparatus and method for processing 3d ground penetrating radar signal | |
EP0943101A1 (en) | Method of detecting and classifying objects by means of radar | |
KR101865016B1 (en) | Cavity detection from ground penetrating radar data and cavity detection apparatus | |
CN112666554A (en) | Method for identifying radar amplitude characteristic crack width of asphalt pavement | |
CN102621531A (en) | Rainfall interference suppression method based on X-band radar images | |
CN108710888A (en) | A kind of Coherent Noise in GPR Record method for registering | |
Wang et al. | Automatic asphalt layer interface detection and thickness determination from ground-penetrating radar data | |
CN111679275A (en) | Underground pipeline identification method based on ground penetrating radar | |
Economou et al. | Attenuation analysis of real GPR wavelets: The equivalent amplitude spectrum (EAS) | |
CN113514833B (en) | Sea surface arbitrary point wave direction inversion method based on sea wave image | |
CN110991507A (en) | Road underground cavity identification method, device and system based on classifier | |
CN111007464B (en) | Road underground cavity identification method, device and system based on optimal weighting | |
CN107656270A (en) | The measurement apparatus and measuring method of a kind of contactless buried pipe track forces cun | |
CN106443674A (en) | Ground penetrating radar wave velocity estimation method based on diffraction, imaging and minimum entropy technology | |
CN102518422A (en) | Method for detecting and identifying current stress of oil field downhole casing | |
CN108303745A (en) | A kind of inversion method of the buried cable detection based on electromagnetic wave saturating ground technology | |
CN111691876B (en) | Method, device and storage medium for imaging adjacent well by using acoustic logging | |
Shen et al. | An alternative method for surface current extraction from X-band marine radar images | |
RU2803396C1 (en) | Method for detecting objects and determining their location in real time using distributed optic fibre interferometric vibration sensors | |
CN112464777B (en) | Intelligent estimation method for vertical distance of optical fiber vibration source | |
JP2551952B2 (en) | Invisible object detection method |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
TA01 | Transfer of patent application right |
Effective date of registration: 20200827 Address after: No.86-a2112, Wanxing Road, Changyang, Fangshan District, Beijing 102400 Applicant after: Beijing Zhongke blueprints Technology Co.,Ltd. Applicant after: Zhongke yuntu Technology Co., Ltd Address before: No.86-a2112, Wanxing Road, Changyang, Fangshan District, Beijing 102400 Applicant before: Beijing Zhongke blueprints Technology Co.,Ltd. |
|
TA01 | Transfer of patent application right | ||
GR01 | Patent grant | ||
GR01 | Patent grant |