CN112882023A - Method, medium and equipment for suppressing shielding interference of steel bar mesh in ground penetrating radar data - Google Patents
Method, medium and equipment for suppressing shielding interference of steel bar mesh in ground penetrating radar data Download PDFInfo
- Publication number
- CN112882023A CN112882023A CN202110077643.6A CN202110077643A CN112882023A CN 112882023 A CN112882023 A CN 112882023A CN 202110077643 A CN202110077643 A CN 202110077643A CN 112882023 A CN112882023 A CN 112882023A
- Authority
- CN
- China
- Prior art keywords
- data
- ground penetrating
- penetrating radar
- along
- steel bar
- 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
- 230000000149 penetrating effect Effects 0.000 title claims abstract description 83
- 238000000034 method Methods 0.000 title claims abstract description 43
- 229910000831 Steel Inorganic materials 0.000 title claims abstract description 28
- 239000010959 steel Substances 0.000 title claims abstract description 28
- 230000003014 reinforcing effect Effects 0.000 claims abstract description 43
- 230000004044 response Effects 0.000 claims abstract description 39
- 230000000694 effects Effects 0.000 claims abstract description 23
- 239000011159 matrix material Substances 0.000 claims description 49
- 238000003860 storage Methods 0.000 claims description 16
- 238000000354 decomposition reaction Methods 0.000 claims description 14
- 229910001294 Reinforcing steel Inorganic materials 0.000 claims description 13
- 230000002787 reinforcement Effects 0.000 claims description 12
- 239000013598 vector Substances 0.000 claims description 7
- 239000000126 substance Substances 0.000 claims description 5
- 238000012935 Averaging Methods 0.000 claims description 4
- 238000012856 packing Methods 0.000 claims 1
- 201000010099 disease Diseases 0.000 abstract description 23
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 abstract description 23
- 238000001514 detection method Methods 0.000 abstract description 14
- 239000004567 concrete Substances 0.000 abstract description 8
- 238000003384 imaging method Methods 0.000 abstract description 6
- 238000004590 computer program Methods 0.000 description 9
- 238000012545 processing Methods 0.000 description 9
- 238000010586 diagram Methods 0.000 description 8
- 230000006870 function Effects 0.000 description 5
- 238000003825 pressing Methods 0.000 description 4
- 238000012360 testing method Methods 0.000 description 4
- 230000007547 defect Effects 0.000 description 3
- 238000000926 separation method Methods 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 238000001914 filtration Methods 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 230000002093 peripheral effect Effects 0.000 description 2
- 241001669679 Eleotris Species 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011065 in-situ storage Methods 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 229910052704 radon Inorganic materials 0.000 description 1
- SYUHGPGVQRZVTB-UHFFFAOYSA-N radon atom Chemical compound [Rn] SYUHGPGVQRZVTB-UHFFFAOYSA-N 0.000 description 1
- 230000008707 rearrangement Effects 0.000 description 1
- 239000011150 reinforced concrete Substances 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/14—Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
- G06F17/141—Discrete Fourier transforms
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Remote Sensing (AREA)
- Mathematical Physics (AREA)
- Radar, Positioning & Navigation (AREA)
- Data Mining & Analysis (AREA)
- Computational Mathematics (AREA)
- Mathematical Optimization (AREA)
- Electromagnetism (AREA)
- Mathematical Analysis (AREA)
- Pure & Applied Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- Theoretical Computer Science (AREA)
- General Engineering & Computer Science (AREA)
- Software Systems (AREA)
- Databases & Information Systems (AREA)
- Discrete Mathematics (AREA)
- Algebra (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
The invention discloses a method, medium and equipment for suppressing interference shielding of a steel bar net in ground penetrating radar data. The invention can suppress the shielding effect of the reinforcing mesh in the ground penetrating radar data by separating the reinforcing mesh diffracted waves, so that other ground penetrating radar response signals originally shielded by the reinforcing mesh diffracted waves are more accurately displayed, and the imaging of the internal structure of a detection target such as a ballastless track subgrade is more accurate, thereby accurately detecting the diseases in the railway subgrade. The method is also suitable for detecting internal diseases of various concrete precast slabs with the reinforcing mesh by adopting ground penetrating radar equipment.
Description
Technical Field
The invention belongs to the technical field of ground penetrating radar data processing, and particularly relates to a method, medium and equipment for suppressing interference shielding of a steel bar net in ground penetrating radar data.
Background
The ground penetrating radar is used as a nondestructive detection method and is widely applied to the investigation and detection of the internal diseases of the ballastless track subgrade of the high-speed railway. The ballastless track adopts a preset track bed plate and a supporting layer which are provided with complex structures and contain reinforcing steel meshes, and a track bed plate made of mixed materials such as reinforced concrete and the like to replace a sleeper and a railway ballast of a traditional track, meanwhile, a high-speed railway tunnel body generally contains a support for keeping the stability of surrounding rocks, and the main support structure comprises the reinforcing steel meshes and a steel frame. Therefore, when the ground penetrating radar is used for detecting the ballastless track of the high-speed rail and the high-speed rail tunnel, the obtained data often has strong diffracted waves generated by the steel frame and the reinforcing mesh which are arranged in the ballastless track. The diffracted wave can seriously interfere and even shield ground penetrating radar response signals of internal diseases such as medium and small-scale cracks, cavities and the like of peripheral structures of ballastless track subgrades and high-speed rail tunnel bodies. The data of the ground penetrating radar is difficult to accurately reflect the real underground structure of the ballastless track, so that the problem is brought to the identification of the diseases in the roadbed structure of the ballastless track or the peripheral structure of the tunnel body, and the same problem can be caused when the ground penetrating radar equipment is adopted to detect the internal diseases of other various concrete precast slabs with reinforcing mesh. Therefore, a certain method needs to be adopted to effectively separate the strong diffracted waves of the reinforcing mesh in the ground penetrating radar data. By separating the ground penetrating radar data after the shielding effect of the strong diffraction wave pressing reinforcing mesh, shielded internal diseases and ground penetrating radar response signals generated by other detection targets are more clearly displayed. The data can be used for more accurate imaging of the internal structure of the detection target, so that the detection target, such as the internal diseases of the roadbed, including possible holes, cracks and other diseases in structures such as prefabricated roadbed slabs, reinforcing meshes in supporting layers and the like, can be accurately detected. In addition, when ground penetrating radar equipment is adopted to detect internal defects of various other concrete precast slabs with reinforcing steel bar nets, in order to enable ground penetrating radar data to accurately reflect internal structure conditions of various concrete precast slabs, the shielding effect of the reinforcing steel bar nets needs to be suppressed.
The traditional method for suppressing the shielding effect of the ground penetrating radar data by separating the strong diffracted waves of the steel bar mesh in the ground penetrating radar data comprises a Radon transform method, a Singular Value Decomposition (SVD) filtering method and an offset, wherein the methods all achieve a certain separation effect, but when the strong diffracted waves of the steel bar mesh and the response signals of the ground penetrating radar target are mixed in space, the strong diffracted waves of the steel bar mesh and the response signals of the ground penetrating radar target are removed from input data together, so that the existing method needs to be improved to achieve a better effect.
The existing SVD filtering method has the following disadvantages:
1. the number of singular values of other ground penetrating radar response signals except the strong diffraction wave of the reinforcing steel bar mesh needs to be repeatedly tested or calculated, and the strong diffraction wave of the reinforcing steel bar mesh is easy to remain or damage other ground penetrating radar response signals.
2. When the strong diffracted waves of the reinforcing mesh and other ground penetrating radar response signals are mixed in space, the other ground penetrating radar response signals and the diffracted waves are filtered out easily.
Disclosure of Invention
The technical problem to be solved by the present invention is to provide a method, medium and device for suppressing interference shielding of a mesh reinforcement in ground penetrating radar data, which can be used for more accurate imaging of the internal structure of a detection target, thereby accurately detecting diseases such as voids, cracks and the like which may exist in structures such as a prefabricated bed plate and a mesh reinforcement in a supporting layer, for example, the internal diseases of a roadbed; meanwhile, the method is suitable for detecting internal diseases of various concrete precast slabs with the reinforcing mesh by adopting ground penetrating radar equipment.
The invention adopts the following technical scheme:
the method for suppressing the shielding interference of the reinforcement mesh in the ground penetrating radar data comprises the following steps:
s1, reading a measured ground penetrating radar B-scan data set;
s2, dividing the B-Scan data set read in the step S1 into a plurality of windows along the space direction, and obtaining a Fourier coefficient matrix F by using one-dimensional discrete Fourier transform on each data in any sub-window X along the time direction;
s3, randomly rearranging each channel of data in the Fourier coefficient matrix F obtained in the step S2 along the spatial direction to obtain F ', arranging the F' into a Hankel matrix, performing singular value decomposition, taking first singular value reconstruction data, rearranging the reconstruction data, and performing one-dimensional inverse discrete Fourier transform on the rearranged reconstruction data to obtain data after the interference shielding of the separated reinforcing mesh;
and S4, repeating the step S2 and the step S3 until the data in all sub-windows in the B-scan data set are processed, and replacing the data at the corresponding position in the original data with the processed sub-window data to obtain the data after the steel bar net diffracted waves are separated and used as the data result after the shielding effect of the pressed steel bar net.
Specifically, in step S2, the fourier coefficient matrix F is expressed as follows:
F＝[f1,…,fm]
wherein f is1,…,fmFor each track of data in F in the spatial direction.
Furthermore, any windowing data X in the plurality of windows comprises two parts, wherein X is XN+(XS1+XS2+…)，XS1,XS2… denotes other ground penetrating radar response signals outside the strong diffracted waves generated by the mesh reinforcement electromagnetic shield, XNIndicating a strong diffraction wave of the steel mesh.
Specifically, step S3 specifically includes:
s301, randomly arranging each channel of data in the Fourier coefficient matrix F along the spatial direction to ensure that all elements are not in the original positions to obtain F';
s302, using data corresponding to the ith frequency point in the F' along the spatial direction to pile up a Hankel matrix Hi；
S303, pair matrix HiUsing singular value decomposition;
s304, reconstructing the first singular value and the corresponding singular vector thereof to obtain H'iTo H'iPerforming inverse diagonal averaging on the intermediate data to obtain reconstructed data corresponding to the ith frequency point along the spatial direction, and repeating the steps S302-S304 until the data corresponding to each frequency point along the spatial direction are processed to obtain the reconstructed data
Will be provided with
The position of each data track along the time direction is rearranged to return to the corresponding position before random arrangement; for those after rearrangement
Obtaining the shielding effect data X ″ -X of the pressed reinforcing steel bar net after one-dimensional inverse discrete Fourier transform is used for each data along the frequency directionS1+XS2+…。
Further, in step S301, the fourier coefficient matrix F' after random arrangement is as follows:
wherein, f'j∈{f′1,…,f′mF is a random Fourier coefficient matrix F without repetition1,…,fmIn the selected data, and
f′n,mis the data of the nth row and the mth column in F'.
Further, in step S302, data along the spatial direction corresponding to the ith frequency point in the fourier coefficient matrix F' after random arrangement is represented as follows:
[f′i，1,…,f′i,m]
matrix HiThe following were used:
wherein q is any positive integer of more than 0 and less than m, f'i,mIs the data of the ith row and the mth column in F'.
Further, in step S303, singular value decomposition is used as follows:
wherein [ u ]1 … uq]，
λγγ is 1, …, and m-q +1 is HiLeft and right singular matrices and singular values of (u) is HiV is HiH is a Hankel matrix which is piled up by the data of the ith frequency point along the space direction in the F'.
Further, in step S304, the obtained H 'is reconstructed'iThe following were used:
wherein the content of the first and second substances,
is H'iThe data of the ith row and the mth column; the reconstructed data expression along the spatial direction corresponding to the ith frequency point is as follows:
fourier coefficient matrix obtained from reconstructed data of all frequency points along spatial direction
The following were used:
X″＝XS1+XS2+ … is as follows:
wherein the content of the first and second substances,
is composed of
Line n and column m, x ″)n,mIs the data of the nth row and the mth column in X'.
Another aspect of the invention is a computer readable storage medium storing one or more programs, the one or more programs comprising instructions, which when executed by a computing device, cause the computing device to perform any of the methods described.
Another aspect of the present invention is a computing device, including:
one or more processors, memory, and one or more programs stored in the memory and configured to be executed by the one or more processors, the one or more programs including instructions for performing any of the methods.
Compared with the prior art, the invention has at least the following beneficial effects:
the invention relates to a method for suppressing the interference shielding of a steel bar net in ground penetrating radar data, which comprises the steps of firstly reading a B-scan data set of the ground penetrating radar, then dividing the read data into a plurality of windows along the spatial direction, obtaining a Fourier coefficient matrix by using one-dimensional Discrete Fourier Transform (DFT) on each channel of data of any windowed data along the time direction, and randomly rearranging each channel of data along the time direction in the Fourier coefficient matrix to enable all elements not to be in situ to obtain new data. Arranging data corresponding to each frequency point in the new data along the spatial direction into Hankel matrixes, performing singular value decomposition on each Hankel matrix, reconstructing by using a first singular value and a singular vector corresponding to the first singular value, and enabling the rearrangement position of each channel of data along the time direction in the reconstructed data to return to the corresponding position before random arrangement. And performing one-dimensional Inverse Discrete Fourier Transform (IDFT) on the rearranged reconstruction data along the frequency direction to obtain data from which diffraction waves are separated. And repeating the steps except for reading the ground penetrating radar ballastless track data set until all the windowing data are processed, and replacing the data of the corresponding position in the original data with the processed windowing data, so as to obtain the ground penetrating radar data with the shielding effect of the pressed reinforcing steel bar net. The shielded internal diseases in the data and ground penetrating radar response signals generated by other detection targets are more clearly displayed, and further the data can be used for more accurately imaging the internal structure of the detection targets, so that the diseases such as holes and cracks which may exist in the detection targets, such as roadbed internal diseases including prefabricated roadbed slabs, reinforcing meshes in supporting layers and other structures, can be accurately detected. The method is suitable for detecting the internal diseases of various concrete precast slabs with the reinforcing mesh by adopting ground penetrating radar equipment. And repeated tests are avoided or the number of singular values of other ground penetrating radar response signals except for reconstruction of diffracted waves is calculated and selected. When other ground penetrating radar response signals and the strong diffracted waves of the reinforcing mesh are mixed in space, the shielding effect of the reinforcing mesh can be effectively suppressed, and other ground penetrating radar response signals are prevented from being damaged.
Further, in step S2, the preliminary separation of the steel bar mesh diffracted wave and other ground penetrating radar response signals is realized through fourier transform.
Furthermore, considering the change of each response in the actual data along the time space direction, windowing along the space direction is utilized to ensure that other ground penetrating radar response signals contained in the data in the window are almost kept unchanged in the time space direction and the variation of the strong diffracted wave of the reinforcing mesh is large, and at the moment, any windowing data X in a plurality of windows is taken for processing, so that the shielding effect of the reinforcing mesh can be better suppressed.
Furthermore, the difference of the spatial coherence of other ground penetrating radar response signals and the strong diffracted wave of the reinforcing mesh is further enhanced.
Further, in step S301, by random arrangement, all elements are not at the original positions to reduce the spatial coherence of the diffracted waves;
further, the coherence of other ground penetrating radar response signals in the space direction is strengthened through the step S302;
further, in steps S303 and S304, the mesh reinforcement diffracted wave and other ground penetrating radar response signals are separated in the form of singular values and singular vectors by singular value decomposition using correlation differences, and the main energy of other ground penetrating radar response signals outside the mesh reinforcement diffracted wave can be completely reconstructed by using the first singular value and the corresponding singular vector thereof, thereby avoiding the number of singular values selected for reconstruction by repeated experimental calculation.
In conclusion, when the method is used for suppressing the shielding effect of the reinforcing steel bar mesh in the ground penetrating radar data set, the number of singular values for reconstructing the response signals of the ground penetrating radar except the strong diffracted waves of the reinforcing steel bar mesh, which are selected by repeated tests, can be avoided; when the strong diffraction wave of the reinforcing mesh and other ground penetrating radar response signals are mixed in space, the invention effectively suppresses the shielding effect of the reinforcing mesh and avoids damaging other ground penetrating radar response signals.
The technical solution of the present invention is further described in detail by the accompanying drawings and embodiments.
Drawings
FIG. 1 is a two-dimensional cross-sectional view of a measured georadar data set;
fig. 2 is a two-dimensional cross-sectional view of the isolated mesh reinforcement of fig. 1 after diffraction;
fig. 3 is a two-dimensional cross-sectional view of the isolated reinforcing mesh data from fig. 1;
FIG. 4 is a flow chart of the present invention.
Detailed Description
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, 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 will be understood that the terms "comprises" and/or "comprising," when used in this specification and the appended claims, specify the presence of stated features, integers, steps, operations, elements, and/or components, but do not preclude the presence or addition of one or more other features, integers, steps, operations, elements, components, and/or groups thereof.
It is also to be understood that the terminology used in the description of the invention herein is for the purpose of describing particular embodiments only and is not intended to be limiting of the invention. As used in the specification of the present invention and the appended claims, the singular forms "a," "an," and "the" are intended to include the plural forms as well, unless the context clearly indicates otherwise.
It should be further understood that the term "and/or" as used in this specification and the appended claims refers to and includes any and all possible combinations of one or more of the associated listed items.
Various structural schematics according to the disclosed embodiments of the invention are shown in the drawings. The figures are not drawn to scale, wherein certain details are exaggerated and possibly omitted for clarity of presentation. The shapes of various regions, layers and their relative sizes and positional relationships shown in the drawings are merely exemplary, and deviations may occur in practice due to manufacturing tolerances or technical limitations, and a person skilled in the art may additionally design regions/layers having different shapes, sizes, relative positions, according to actual needs.
The invention provides a method for suppressing the shielding interference of a steel bar mesh in ground penetrating radar data, which takes any windowing data in a plurality of windows into which a ground penetrating radar B-Scan data set is divided along the space direction as two parts: the steel bar net diffracted waves and other ground penetrating radar response signals (mainly comprising subgrade diseases and ground penetrating radar response signals generated by subgrade stratums). Firstly, reading a B-Scan data set of the ground penetrating radar, dividing the B-Scan data set into a plurality of windows along the spatial direction, obtaining a Fourier coefficient matrix by using one-dimensional Discrete Fourier Transform (DFT) on each channel of data in any one of the windows along the time direction, randomly rearranging the Fourier coefficient matrix along the spatial direction to ensure that all elements are not in the original positions, arranging data corresponding to each frequency point in the rearranged matrix along the spatial direction into a Hankel matrix, performing singular value decomposition on each Hankel matrix, reconstructing by using a first singular value and a singular vector corresponding to the first singular value, and rearranging each channel of data along the time direction in reconstructed data to ensure that the rearranged position returns to the corresponding position before random arrangement. And performing one-dimensional Inverse Discrete Fourier Transform (IDFT) on each data in the rearranged reconstruction data along the frequency direction to obtain the data after the separated reinforcing mesh shields the interference. And repeating the steps except reading the ground penetrating radar data set until all the sub-windows in the data set are processed, replacing the data of the corresponding position in the original data with the processed sub-window data, and displaying the shielded internal diseases and ground penetrating radar response signals generated by other detection targets more clearly in the ground penetrating radar data with the shielding effect of the pressed reinforcing steel bar mesh. The data can be used for more accurate imaging of the internal structure of the detection target, so that the detection target, such as the internal diseases of the roadbed, including possible holes, cracks and other diseases in structures such as prefabricated roadbed slabs, reinforcing meshes in supporting layers and the like, can be accurately detected. The method is also suitable for detecting the internal diseases of various concrete precast slabs with the reinforcing mesh by adopting ground penetrating radar equipment.
Referring to fig. 4, the method for suppressing interference shielding of the reinforcement mesh in the ground penetrating radar data of the present invention includes the following steps:
s1, reading a measured ground penetrating radar B-scan data set;
s2, dividing the B-scan data of the ground penetrating radar into a plurality of windows along the space direction, wherein any windowing data X comprises two parts, and X is equal to XN+(XS1+XS2+…)，XS1,XS2… denotes other ground penetrating radar response signals outside the strong diffracted wave, XNIndicating a strong diffraction wave of the steel mesh.
Please refer to fig. 1, which is a cross-sectional view corresponding to two-dimensional cross-sectional data X. As can be seen from the figure, the conical homophase axes are distributed on the whole two-dimensional section, which is mainly the strong diffraction wave X of the reinforcing meshN. Other ground penetrating radar response signals XS1,XS2… comprising ballastless track subgrade knotThe defects in the structure include a hole reflection wave (for example, a region indicated by a white arrow), a roadbed formation reflection wave having a nearly horizontal in-phase axis, and the like. Because the reflected waves of the holes and the roadbed stratum and the strong diffracted waves of the reinforcing mesh are overlapped in space, the in-phase axis of the holes and the roadbed stratum is broken or even shielded. And obtaining a Fourier coefficient matrix F matrix by using one-dimensional discrete Fourier transform along the time direction for each channel of X data, wherein the expression is as follows:
F＝[f1,…,fm]
wherein f is1,…,fmFor each track of data in F in the spatial direction.
S3, randomly arranging each channel of data in the F along the spatial direction to enable all elements not to be in the original positions to obtain the F ', arranging the data in the F' corresponding to each frequency point along the spatial direction into a Hankel matrix, performing singular value decomposition on the Hankel matrix, and taking the first singular value to reconstruct the data. And rearranging the position of each data in the reconstructed data to enable the data to return to the position before random arrangement, and performing one-dimensional inverse discrete Fourier transform on the rearranged data along the frequency direction to obtain the data with the shielding effect of the steel bar net pressed. The method comprises the following specific steps:
s301, randomly arranging each channel of data in the Fourier coefficient matrix F along the spatial direction to enable all elements not to be in the original positions to obtain F', wherein the expression is as follows;
wherein, f'j∈{f′1,…,f′mF is a random Fourier coefficient matrix F without repetition1,…,fmIn the selected data, and
f′n,mis the data of the nth row and the mth column in F'.
S302, using data corresponding to the ith frequency point in the F' along the spatial direction to pile up a Hankel matrix HiData representation in the spatial direction corresponding to the ith frequency point in FHas the following formula
[f′i，1,…,f′i,m]
Matrix HiThe following were used:
wherein q is any positive integer of more than 0 and less than m, f'i,mData of ith row and mth column in F';
s303, pair matrix HiUsing singular value decomposition, the expression is as follows;
wherein [ u ]1 … uq]，
λγγ is 1, …, and m-q +1 is HiLeft and right singular matrices and singular values.
S304, reconstructing the first singular value and the corresponding singular vector thereof to obtain H'iTo H'iPerforming inverse diagonal averaging on the intermediate data to obtain reconstructed data corresponding to the ith frequency point along the spatial direction, and repeating the steps S302-S304 until the data corresponding to each frequency point along the spatial direction are processed to obtain the reconstructed data
Will be provided with
The position of each data track along the time direction is rearranged to return to the corresponding position before random arrangement; for those after rearrangement
Obtaining the shielding effect number of the pressed reinforcing steel bar net by using one-dimensional inverse discrete Fourier transform of each data along the frequency directionAccording to X ″ ═ XS1+XS2+…。
Wherein the content of the first and second substances,
is H'iThe data of the ith row and the mth column; to H'iPerforming inverse diagonal averaging on the medium data to obtain a reconstructed data expression corresponding to the ith frequency point along the spatial direction as follows:
repeating steps S302-S304 until each frequency slice is processed, resulting in the following:
wherein, f'n,mIs the data of the nth row and the mth column in F'.
After one-dimensional IDFT is used for each spatial position of F', obtaining
Wherein, x ″)n,mThe data of the nth row and the mth column in X' is the data of the separated diffraction wave after the separation.
After all the windowing data are processed, replacing the data at the corresponding positions in the original data with the processed windowing data, and obtaining the data for pressing the reinforcing mesh by separating the shielding effect of the diffracted waves as shown in fig. 2, wherein compared with the same area indicated by arrows in fig. 1 and fig. 2, the strong diffracted waves of the tapered reinforcing mesh in fig. 2 are obviously pressed; and the subgrade defect shielded by the diffracted waves in fig. 1 and the response signal of the ground penetrating radar of the subgrade stratum are also more clearly shown in fig. 2. The separated strong diffraction wave data are shown in fig. 3, and only the mesh reinforcement diffraction waves can be seen in fig. 3, which shows that the mesh reinforcement diffraction waves are sufficiently separated while other ground penetrating radar response signals are not damaged.
The practical data calculation example shows that the method can avoid repeated test selection of the number of singular values for reconstructing the response signals of the ground penetrating radar except the strong diffracted waves of the reinforcing mesh. When the strong diffraction wave of the reinforcing mesh and other ground penetrating radar response signals are mixed in space, the invention can still effectively suppress the shielding effect of the reinforcing mesh and avoid damaging other ground penetrating radar response signals.
In yet another embodiment of the present invention, a terminal device is provided that includes a processor and a memory for storing a computer program comprising program instructions, the processor being configured to execute the program instructions stored by the computer storage medium. The Processor may be a Central Processing Unit (CPU), or may be other general purpose Processor, a Digital Signal Processor (DSP), an Application Specific Integrated Circuit (ASIC), an off-the-shelf Programmable gate array (FPGA) or other Programmable logic device, a discrete gate or transistor logic device, a discrete hardware component, etc., which is a computing core and a control core of the terminal, and is adapted to implement one or more instructions, and is specifically adapted to load and execute one or more instructions to implement a corresponding method flow or a corresponding function; the processor provided by the embodiment of the invention can be used for the operation of the pressing method for shielding interference of the steel bar net in the ground penetrating radar data, and comprises the following steps: reading a measured ground penetrating radar B-scan data set; dividing the B-Scan data set into a plurality of windows along the space direction, and obtaining a Fourier coefficient matrix F by using one-dimensional discrete Fourier transform on each data in any sub-window X along the time direction; randomly rearranging each channel of data in the Fourier coefficient matrix F along the spatial direction to obtain F ', arranging F' into a Hankel matrix, performing singular value decomposition, taking first singular value reconstruction data, rearranging the reconstruction data, and performing one-dimensional inverse discrete Fourier transform on the rearranged reconstruction data to obtain data after the reinforcing mesh is separated from interference; and repeating the steps until the data processing in all the sub-windows in the data set is finished, and replacing the data at the corresponding position in the original data with the processed sub-window data to obtain the data after the steel bar mesh diffracted wave is separated and used as the data result after the steel bar mesh is pressed for shielding.
In still another embodiment of the present invention, the present invention further provides a storage medium, specifically a computer-readable storage medium (Memory), which is a Memory device in a terminal device and is used for storing programs and data. It is understood that the computer readable storage medium herein may include a built-in storage medium in the terminal device, and may also include an extended storage medium supported by the terminal device. The computer-readable storage medium provides a storage space storing an operating system of the terminal. Also, one or more instructions, which may be one or more computer programs (including program code), are stored in the memory space and are adapted to be loaded and executed by the processor. It should be noted that the computer-readable storage medium may be a high-speed RAM memory, or may be a non-volatile memory (non-volatile memory), such as at least one disk memory.
The processor can load and execute one or more than one instruction stored in the computer readable storage medium to realize the corresponding steps of the method for suppressing the interference shielding of the steel bar mesh in the ground penetrating radar data in the embodiment; one or more instructions in the computer-readable storage medium are loaded by the processor and perform the steps of: reading a measured ground penetrating radar B-scan data set; dividing the B-Scan data set into a plurality of windows along the space direction, and obtaining a Fourier coefficient matrix F by using one-dimensional discrete Fourier transform on each data in any sub-window X along the time direction; randomly rearranging each channel of data in the Fourier coefficient matrix F along the spatial direction to obtain F ', arranging F' into a Hankel matrix, performing singular value decomposition, taking first singular value reconstruction data, rearranging the reconstruction data, and performing one-dimensional inverse discrete Fourier transform on the rearranged reconstruction data to obtain data after the reinforcing mesh is separated from interference; and repeating the steps until the data processing in all the sub-windows in the data set is finished, and replacing the data at the corresponding position in the original data with the processed sub-window data to obtain the data after the steel bar mesh diffracted wave is separated and used as the data result after the steel bar mesh is pressed for shielding.
The invention can suppress the shielding effect of the reinforcing mesh in the ground penetrating radar data by separating the reinforcing mesh diffracted waves, so that other ground penetrating radar response signals originally shielded by the reinforcing mesh diffracted waves are more accurately displayed, and the imaging of the internal structure of a detection target such as a ballastless track subgrade is more accurate, thereby accurately detecting the diseases in the railway subgrade. Meanwhile, the method is also suitable for detecting the internal diseases of various concrete precast slabs with the reinforcing mesh by adopting ground penetrating radar equipment.
In summary, the pressing method, medium and equipment for shielding interference of the reinforcement mesh in the ground penetrating radar data have the following advantages:
1, the shielding effect of the ballastless track reinforcing mesh is suppressed by using the method, so that the number of singular values for reconstructing the response signals of the ground penetrating radar except the strong diffracted wave of the reinforcing mesh can be avoided from being selected in repeated tests.
2, when the strong diffracted wave of the reinforcing mesh and other ground penetrating radar response signals are mixed in space, the method effectively suppresses the shielding effect of the reinforcing mesh of the ballastless track and avoids damaging other ground penetrating radar response signals.
As will be appreciated by one skilled in the art, embodiments of the present application may be provided as a method, system, or computer program product. Accordingly, the present application may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present application may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.
The present application is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the application. It will be understood that each flow and/or block of the flow diagrams and/or block diagrams, and combinations of flows and/or blocks in the flow diagrams and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
The above-mentioned contents are only for illustrating the technical idea of the present invention, and the protection scope of the present invention is not limited thereby, and any modification made on the basis of the technical idea of the present invention falls within the protection scope of the claims of the present invention.
Claims (10)
1. The method for suppressing the shielding interference of the steel bar mesh in the ground penetrating radar data is characterized by comprising the following steps of:
s1, reading a measured ground penetrating radar B-scan data set;
s2, dividing the B-Scan data set read in the step S1 into a plurality of windows along the space direction, and obtaining a Fourier coefficient matrix F by using one-dimensional discrete Fourier transform on each data in any sub-window X along the time direction;
s3, randomly rearranging each channel of data in the Fourier coefficient matrix F obtained in the step S2 along the spatial direction to obtain F ', arranging the F' into a Hankel matrix, performing singular value decomposition, taking first singular value reconstruction data, rearranging the reconstruction data, and performing one-dimensional inverse discrete Fourier transform on the rearranged reconstruction data to obtain data after the interference shielding of the separated reinforcing mesh;
and S4, repeating the step S2 and the step S3 until the data in all sub-windows in the B-scan data set are processed, and replacing the data at the corresponding position in the original data with the processed sub-window data to obtain the data after the steel bar net diffracted waves are separated and used as the data result after the shielding effect of the pressed steel bar net.
2. The method according to claim 1, wherein in step S2, the fourier coefficient matrix F is represented as follows:
F＝[f1,…,fm]
wherein f is1,…,fmFor each track of data in F in the spatial direction.
3. The method of claim 2, wherein any of the windowed data X in the plurality of windows comprises two portions, X ═ XN+(XS1+XS2+…)，XS1,XS2… denotes other ground penetrating radar response signals outside the strong diffracted waves generated by the mesh reinforcement electromagnetic shield, XNIndicating a strong diffraction wave of the steel mesh.
4. The method according to claim 1, wherein step S3 is specifically:
s301, randomly arranging each channel of data in the Fourier coefficient matrix F along the spatial direction to ensure that all elements are not in the original positions to obtain F';
s302, using the frequency point along the space direction corresponding to the ith frequency point in the FData packing Hankel matrix Hi；
S303, pair matrix HiUsing singular value decomposition;
s304, reconstructing the first singular value and the corresponding singular vector thereof to obtain H'iTo H'iPerforming inverse diagonal averaging on the intermediate data to obtain reconstructed data corresponding to the ith frequency point along the spatial direction, and repeating the steps S302-S304 until the data corresponding to each frequency point along the spatial direction are processed to obtain the reconstructed data
Will be provided with
The position of each data track along the time direction is rearranged to return to the corresponding position before random arrangement; for those after rearrangement
Obtaining the shielding effect data X ″ -X of the pressed reinforcing steel bar net after one-dimensional inverse discrete Fourier transform is used for each data along the frequency directionS1+XS2+…。
5. The method according to claim 4, wherein in step S301, the matrix F' of Fourier coefficients after random arrangement is as follows:
wherein, f'j∈{f′1,…,f′mF is a random Fourier coefficient matrix F without repetition1,…,fmIn the selected data, and
f′n,mis the data of the nth row and the mth column in F'.
6. The method according to claim 4, wherein in step S302, the data along the spatial direction corresponding to the ith frequency point in the fourier coefficient matrix F' after random arrangement is represented as follows:
[f′i，1,…,f′i,m]
matrix HiThe following were used:
wherein q is any positive integer of more than 0 and less than m, f'i,mIs the data of the ith row and the mth column in F'.
7. The method of claim 4, wherein in step S303, singular value decomposition is used as follows:
wherein [ u ]1 … uq]，
λγγ is 1, …, and m-q +1 is HiLeft and right singular matrices and singular values of (u) is HiV is HiH is a Hankel matrix which is piled up by the data of the ith frequency point along the space direction in the F'.
8. The method of claim 4, wherein in step S304, the obtained H is reconstructedi' the following:
wherein the content of the first and second substances,
is H'iThe data of the ith row and the mth column; the reconstructed data expression along the spatial direction corresponding to the ith frequency point is as follows:
fourier coefficient matrix obtained from reconstructed data of all frequency points along spatial direction
The following were used:
X″＝XS1+XS2+ … is as follows:
wherein the content of the first and second substances,
is composed of
Line n and column m of data, x'n,mIs the data of the nth row and the mth column in X'.
9. A computer readable storage medium storing one or more programs, the one or more programs comprising instructions, which when executed by a computing device, cause the computing device to perform any of the methods of claims 1-8.
10. A computing device, comprising:
one or more processors, memory, and one or more programs stored in the memory and configured for execution by the one or more processors, the one or more programs including instructions for performing any of the methods of claims 1-8.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110077643.6A CN112882023A (en) | 2021-01-20 | 2021-01-20 | Method, medium and equipment for suppressing shielding interference of steel bar mesh in ground penetrating radar data |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110077643.6A CN112882023A (en) | 2021-01-20 | 2021-01-20 | Method, medium and equipment for suppressing shielding interference of steel bar mesh in ground penetrating radar data |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN112882023A true CN112882023A (en) | 2021-06-01 |
Family
ID=76051147
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110077643.6A Pending CN112882023A (en) | 2021-01-20 | 2021-01-20 | Method, medium and equipment for suppressing shielding interference of steel bar mesh in ground penetrating radar data |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112882023A (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117115049A (en) * | 2023-10-17 | 2023-11-24 | 深圳市城市交通规划设计研究中心股份有限公司 | Noise reduction method for road interior image of ground penetrating radar, electronic equipment and storage medium |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102854533A (en) * | 2011-07-01 | 2013-01-02 | 中国石油化工股份有限公司 | Wave field separation principle based denoising method for increasing signal to noise ratio of seismic data |
| CN104360393A (en) * | 2014-10-20 | 2015-02-18 | 李闯 | Seismic data reconstruction method |
| US20180259664A1 (en) * | 2017-03-08 | 2018-09-13 | Saudi Arabian Oil Company | Automated system and methods for adaptive robust denoising of large-scale seismic data sets |
| CN108983158A (en) * | 2018-04-16 | 2018-12-11 | 中国地质大学（武汉） | A kind of Ground Penetrating Radar noise suppressing method based on Hankel Singular Value Decomposition Using |
| CN110826017A (en) * | 2019-09-25 | 2020-02-21 | 中国地质大学（武汉） | Signal denoising method based on parameter optimization Hankel matrix and singular value decomposition |
| CN111290024A (en) * | 2020-03-05 | 2020-06-16 | 吉林大学 | SVD self-adaptive seismic data noise suppression method |
| CN111736224A (en) * | 2020-07-14 | 2020-10-02 | 西安交通大学 | Method, storage medium and equipment for suppressing linear interference of pre-stack seismic data |
-
2021
- 2021-01-20 CN CN202110077643.6A patent/CN112882023A/en active Pending
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102854533A (en) * | 2011-07-01 | 2013-01-02 | 中国石油化工股份有限公司 | Wave field separation principle based denoising method for increasing signal to noise ratio of seismic data |
| CN104360393A (en) * | 2014-10-20 | 2015-02-18 | 李闯 | Seismic data reconstruction method |
| US20180259664A1 (en) * | 2017-03-08 | 2018-09-13 | Saudi Arabian Oil Company | Automated system and methods for adaptive robust denoising of large-scale seismic data sets |
| CN108983158A (en) * | 2018-04-16 | 2018-12-11 | 中国地质大学（武汉） | A kind of Ground Penetrating Radar noise suppressing method based on Hankel Singular Value Decomposition Using |
| CN110826017A (en) * | 2019-09-25 | 2020-02-21 | 中国地质大学（武汉） | Signal denoising method based on parameter optimization Hankel matrix and singular value decomposition |
| CN111290024A (en) * | 2020-03-05 | 2020-06-16 | 吉林大学 | SVD self-adaptive seismic data noise suppression method |
| CN111736224A (en) * | 2020-07-14 | 2020-10-02 | 西安交通大学 | Method, storage medium and equipment for suppressing linear interference of pre-stack seismic data |
Non-Patent Citations (4)
| Title |
|---|
| 崔少华等: "地震数据去噪的研究与分析", 《廊坊师范学院学报(自然科学版)》 * |
| 曹中林等: "基于复数域混合SVD滤波方法及在随机噪声压制中的应用", 《地球物理学进展》 * |
| 蔡文芮: "基于Hankel矩阵的随机噪音衰减法在煤田勘探中的应用", 《物探化探计算技术》 * |
| 鲍伟等: "频率域奇异值分解压制随机噪声方法", 《科技导报》 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117115049A (en) * | 2023-10-17 | 2023-11-24 | 深圳市城市交通规划设计研究中心股份有限公司 | Noise reduction method for road interior image of ground penetrating radar, electronic equipment and storage medium |
| CN117115049B (en) * | 2023-10-17 | 2024-03-29 | 深圳市城市交通规划设计研究中心股份有限公司 | Noise reduction method for road interior image of ground penetrating radar, electronic equipment and storage medium |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Dinh et al. | Automated visualization of concrete bridge deck condition from GPR data | |
| Viggiani et al. | Full-field measurements, a new tool for laboratory experimental geomechanics | |
| Wang et al. | Three-dimensional observations of meso-structural changes in bimsoil using X-ray computed tomography (CT) under triaxial compression | |
| CN112882023A (en) | Method, medium and equipment for suppressing shielding interference of steel bar mesh in ground penetrating radar data | |
| Zhang et al. | An analysis of the simulated acoustic emission sources with different propagation distances, types and depths for rail defect detection | |
| Nakahata et al. | Three dimensional image-based simulation of ultrasonic wave propagation in polycrystalline metal using phase-field modeling | |
| CN104952043B (en) | Image filtering method and CT systems | |
| Jiang et al. | Acoustic emission tomography based on simultaneous algebraic reconstruction technique to visualize the damage source location in Q235B steel plate | |
| CN111736224B (en) | Method, storage medium and equipment for suppressing linear interference of pre-stack seismic data | |
| CN102854532A (en) | Three-dimensional pre-stack offset stochastic noise suppression method | |
| Asadollahi et al. | Analytical reverse time migration: an innovation in imaging of infrastructures using ultrasonic shear waves | |
| Rucka et al. | Interpolation methods in GPR tomographic imaging of linear and volume anomalies for cultural heritage diagnostics | |
| Luo et al. | A dual-scale morphological filtering method for composite damage identification using FBP | |
| CN106019188A (en) | MR saturation taking into consideration anatomical structures to be imaged | |
| Qu et al. | Multi-bandwidth wavelet transform of impact echo signals for delamination detection of thin concrete slabs | |
| CN104181481B (en) | Magnetic resonance image acquiring and reconstructing method | |
| Vermeij et al. | Martensite plasticity and damage competition in dual-phase steel: A micromechanical experimental–numerical study | |
| Wang et al. | Unsupervised learning method for rebar signal suppression and defect signal reconstruction and detection in ground penetrating radar images | |
| Jeng et al. | Application of sub-image multiresolution analysis of Ground-penetrating radar data in a study of shallow structures | |
| Lei et al. | GPR detection localization of underground structures based on deep learning and reverse time migration | |
| Albright et al. | The benefits of using time-frequency analysis with synthetic aperture focusing technique | |
| Yeh et al. | Imaging of internal cracks in concrete structures using the surface rendering technique | |
| CN104166114B (en) | FID and/or DC artifact correction method and device, MRI collection and method for reconstructing | |
| Graziani et al. | Modern style nappe stacking in the Paleoproterozoic lower crust: An example from the snowbird tectonic zone, Canadian Shield | |
| Alebrahim | Modified wave dispersion properties in 1D and 2D state-based peridynamic media |
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 | ||
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20210601 |