Methods and Systems for Calculating Electromagnetic Wave Velocity in Tunnel Lining Radar Data

By using the FK offset algorithm and entropy calculation, the problem of inaccurate calculation of electromagnetic wave velocity in tunnel lining radar data is solved, achieving high resolution of radar images. This method is applicable to radar data offset processing of steel bars and steel arches.

CN117949895BActive Publication Date: 2026-05-26BEIJING JIAOTONG UNIV

Patent Information

Authority / Receiving Office
CN · China
Patent Type
Patents(China)
Current Assignee / Owner
BEIJING JIAOTONG UNIV
Filing Date
2023-12-12
Publication Date
2026-05-26

AI Technical Summary

Technical Problem

In existing technologies, the electromagnetic wave velocity in tunnel lining radar data is difficult to determine accurately, resulting in insufficient radar image resolution and limiting its widespread application.

Method used

The FK offset algorithm is used for radar image offset processing. The radar image entropy value corresponding to the electromagnetic wave velocity is calculated, and a velocity-entropy piecewise line graph is established. The accurate electromagnetic wave velocity value is obtained by iteratively searching for the minimum image entropy value.

Benefits of technology

It enables accurate calculation of electromagnetic wave velocity in radar data of tunnel lining, improves the resolution of radar images, and is applicable to radar data offset processing of steel bars and steel arches.

✦ Generated by Eureka AI based on patent content.

Smart Images

  • Figure CN117949895B_ABST
    Figure CN117949895B_ABST
Patent Text Reader

Abstract

This invention provides a method and system for calculating electromagnetic wave velocity in tunnel lining radar data, belonging to the field of tunnel lining quality inspection technology. It utilizes a frequency-wavenumber domain migration algorithm to migrate radar data, and then calculates the image entropy value of the migrated radar data to plot a velocity-entropy curve of the radar image. Finally, by iteratively searching for the electromagnetic wave velocity value corresponding to the minimum image entropy value, the precise electromagnetic wave velocity value in the radar data can be obtained. This invention, through its research on a precise calculation method for electromagnetic wave velocity in tunnel lining radar data, overcomes the deficiency of ground-penetrating radar in accurately calculating electromagnetic wave velocity in tunnel lining radar data, providing a theoretical foundation and data support for subsequent research on electromagnetic wave velocity in tunnel lining radar data.
Need to check novelty before this filing date? Find Prior Art

Description

Technical Field

[0001] This invention relates to the field of tunnel lining quality inspection technology, specifically to a method and system for calculating electromagnetic wave velocity in tunnel lining radar data. Background Technology

[0002] Tunnel engineering is a common type of project in road and bridge infrastructure, municipal facilities, and industrial and mining enterprises. Because tunnels involve the safety of road and bridge transportation and industrial and mining production, real-time monitoring and safety early warning of tunnel conditions are necessary. Among tunnel conditions, the tunnel lining is a crucial parameter reflecting the tunnel's normal operation, and its inspection is an important aspect of tunnel safety. Currently, ground penetrating radar (GPR) is commonly used for real-time monitoring of tunnel conditions.

[0003] Ground-penetrating radar (GPR), a non-destructive testing technology, primarily relies on processing radar echo signals to identify buried targets. When used for tunnel quality inspection, GPR can obtain high-resolution images without damaging the tunnel structure, and this technology has relatively low requirements for the on-site working environment. However, the electromagnetic wave velocity values ​​in radar images are often difficult to determine, leading to insufficient image resolution during data acquisition and image interpretation, thus hindering its widespread application. Summary of the Invention

[0004] The purpose of this invention is to provide a method and system for calculating electromagnetic wave velocity in radar data of tunnel lining, so as to solve at least one of the technical problems existing in the background art.

[0005] To achieve the above objectives, the present invention adopts the following technical solution:

[0006] On one hand, the present invention provides a method for calculating electromagnetic wave velocity in radar data of tunnel lining, comprising:

[0007] Given a velocity value, perform radar image offset processing based on the FK offset algorithm;

[0008] Calculate the entropy value of the radar image after offset corresponding to a given electromagnetic wave velocity value;

[0009] Based on the correspondence between velocity and image entropy, a velocity-entropy line graph is constructed;

[0010] Based on the velocity-entropy line graph, the precise electromagnetic wave velocity value in the radar data is obtained by iteratively searching for the electromagnetic wave velocity value corresponding to the minimum image entropy value.

[0011] Preferably, given the range and step size of the electromagnetic wave velocity, the initial electromagnetic wave velocity value is set to 0.1 m / ns; the electromagnetic wave velocity value is determined by... The calculation yields V, where n represents the number value of the nth iteration. n This represents the electromagnetic wave speed value calculated in the nth iteration.

[0012] Preferably, the one-dimensional scan signal in the image entropy value is derived from... It means that x ij X represents the amplitude value of the sampling point; j=T Represents a column vector consisting of sampled points; the minimum image entropy value of the offset radar image is obtained through the formula... Calculate, where X represents a matrix vector consisting of sampling points; E(X) represents the image entropy value of the offset image.

[0013] Preferably, the smaller the image entropy value in the velocity-entropy curve of the radar image, the better the offset effect. Therefore, by iteratively calculating the entropy of the offset radar image, the electromagnetic wave velocity value corresponding to the minimum image entropy value is taken as the accurate electromagnetic wave velocity value.

[0014] Preferably, the number of cycles for the image entropy value can be determined by... Calculated; where V min V max Represents the preset minimum and maximum electromagnetic wave speeds; The value step size represents the preset electromagnetic wave speed; N represents the total number of iterations.

[0015] Preferably, the minimum image entropy corresponding to the electromagnetic wave velocity value can be obtained using the minimum value function in Matlab software; if two identical entropy values ​​are obtained, the larger wave velocity value is taken as the accurate electromagnetic wave velocity value.

[0016] Secondly, the present invention provides a system for calculating electromagnetic wave velocity in radar data of tunnel lining, comprising:

[0017] The offset module is used to perform radar image offset processing based on the FK offset algorithm, given a velocity value.

[0018] The calculation module is used to calculate the entropy value of the radar image after offset corresponding to a given electromagnetic wave velocity value;

[0019] A module is built to create a velocity-entropy line graph based on the relationship between velocity and image entropy.

[0020] The cyclic search module is used to obtain the precise electromagnetic wave velocity value in the radar data by cyclically searching for the electromagnetic wave velocity value corresponding to the minimum image entropy value based on the velocity-entropy line graph.

[0021] Thirdly, the present invention provides a non-transitory computer-readable storage medium for storing computer instructions, which, when executed by a processor, implement the method for calculating electromagnetic wave velocity in tunnel lining radar data as described above.

[0022] Fourthly, the present invention provides a computer device including a memory and a processor, wherein the processor and the memory communicate with each other, the memory stores program instructions that can be executed by the processor, and the processor calls the program instructions to execute the electromagnetic wave velocity calculation method in tunnel lining radar data as described above.

[0023] Fifthly, the present invention provides an electronic device, comprising: a processor, a memory, and a computer program; wherein the processor is connected to the memory, the computer program is stored in the memory, and when the electronic device is running, the processor executes the computer program stored in the memory to cause the electronic device to execute instructions for implementing the method for calculating electromagnetic wave velocity in tunnel lining radar data as described above.

[0024] The beneficial effects of this invention are: it is applicable to the offset processing of radar data of steel bars and steel arches, and can quickly obtain accurate electromagnetic wave velocity values.

[0025] The advantages of additional aspects of the invention will be set forth more clearly in the following description or will be learned by practice of the invention. Attached Figure Description

[0026] To more clearly illustrate the technical solutions of the embodiments of the present invention, the drawings used in the description of the embodiments will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present invention. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.

[0027] Figure 1 This is a flowchart illustrating the precise calculation of electromagnetic wave velocity according to an embodiment of the present invention.

[0028] Figure 2 This is a schematic diagram of the hyperbolic migration process of a discrete single target as described in an embodiment of the present invention.

[0029] Figure 3 This is a line graph of the measured rebar radar data wave velocity-entropy value as described in an embodiment of the present invention.

[0030] Figure 4 This is a measured rebar radar data offset diagram according to an embodiment of the present invention. Figure 4 (a) is the under-offset plot of the measured rebar radar data V = 0.06 m / ns. Figure 4(b) is the correct offset diagram of the measured rebar radar data V = 0.12 m / ns. Figure 4 (c) is the offset map of the measured rebar radar data V = 0.20 m / ns.

[0031] Figure 5 This is a line graph showing the measured radar data wave velocity-entropy value of the steel arch frame as described in an embodiment of the present invention.

[0032] Figure 6 This is a radar data offset diagram of the measured steel arch frame described in an embodiment of the present invention. Figure 6 (a) is the under-offset plot of the measured steel arch frame radar data V=0.02m / ns. Figure 6 (b) is the correct offset diagram of the measured radar data of the steel arch frame V = 0.07 m / ns. Figure 6 (c) is the offset diagram of the measured steel arch radar data V = 0.12m / ns of the present invention. Detailed Implementation

[0033] Embodiments of the present invention are described in detail below, examples of which are shown in the accompanying drawings, wherein the same or similar reference numerals denote the same or similar elements or elements having the same or similar functions throughout. The embodiments described below with reference to the accompanying drawings are exemplary and are only used to explain the present invention, and should not be construed as limiting the present invention.

[0034] It will be understood by those skilled in the art that, unless otherwise defined, all terms used herein (including technical and scientific terms) have the same meaning as commonly understood by one of ordinary skill in the art to which this invention pertains.

[0035] It should also be understood that terms such as those defined in general dictionaries should be understood to have meanings consistent with their meanings in the context of the prior art, and should not be interpreted in an idealized or overly formal sense unless defined as here.

[0036] Those skilled in the art will understand that, unless specifically stated otherwise, the singular forms “a,” “an,” “the,” and “the” used herein may also include the plural forms. It should be further understood that the term “comprising” as used in this specification means the presence of the stated features, integers, steps, operations, elements, and / or components, but does not exclude the presence or addition of one or more other features, integers, steps, operations, elements, and / or groups thereof.

[0037] In the description of this specification, references to terms such as "one embodiment," "some embodiments," "example," "specific example," or "some examples," etc., indicate that a specific feature, structure, material, or characteristic described in connection with that embodiment or example is included in at least one embodiment or example of the present invention. Furthermore, the specific features, structures, materials, or characteristics described may be combined in any suitable manner in one or more embodiments or examples. Moreover, without contradiction, those skilled in the art can combine and integrate the different embodiments or examples described in this specification, as well as the features of those different embodiments or examples.

[0038] To facilitate understanding of the present invention, the present invention will be further explained and described below with reference to the accompanying drawings and specific embodiments. However, the specific embodiments do not constitute a limitation on the embodiments of the present invention.

[0039] Those skilled in the art should understand that the accompanying drawings are merely schematic diagrams of embodiments, and the components in the drawings are not necessarily essential for implementing the present invention.

[0040] Example 1

[0041] In this embodiment 1, an electromagnetic wave velocity calculation system for tunnel lining radar data is first provided, including: an offset module for performing radar image offset processing based on the FK offset algorithm for a given velocity value; a calculation module for calculating the entropy value of the offset radar image corresponding to the given electromagnetic wave velocity value; a construction module for establishing a velocity-entropy line graph based on the correspondence between velocity and image entropy; and a loop search module for obtaining the accurate electromagnetic wave velocity value in the radar data by iteratively searching for the electromagnetic wave velocity value corresponding to the minimum image entropy value based on the velocity-entropy line graph.

[0042] In this embodiment, the above-described system is used to implement a method for calculating electromagnetic wave velocity in tunnel lining radar data, including: using an offset module to provide a velocity value and performing radar image offset processing based on the FK offset algorithm; using a calculation module to calculate the entropy value of the offset radar image corresponding to the given electromagnetic wave velocity value; using a construction module to establish a velocity-entropy line graph based on the correspondence between velocity and image entropy; and finally, using a loop search module to obtain the accurate electromagnetic wave velocity value in the radar data by iteratively searching for the electromagnetic wave velocity value corresponding to the minimum image entropy value based on the velocity-entropy line graph.

[0043] Wherein, given the range and step size of the electromagnetic wave velocity, the initial electromagnetic wave velocity value is set to 0.1 m / ns; the electromagnetic wave velocity value is determined by... The calculation yields V, where n represents the number value of the nth iteration. n This represents the electromagnetic wave speed value calculated in the nth iteration.

[0044] Among them, the one-dimensional scan signal in the image entropy value is composed of It means that x ij X represents the amplitude value of the sampling point; j=T Represents a column vector consisting of sampled points; the minimum image entropy value of the offset radar image is obtained through the formula... Calculate, where X represents a matrix vector consisting of sampling points; E(X) represents the image entropy value of the offset image.

[0045] The smaller the image entropy value in the velocity-entropy curve of the radar image, the better the offset effect. Therefore, by iteratively calculating the entropy of the offset radar image, the electromagnetic wave velocity value corresponding to the minimum image entropy value is taken as the accurate electromagnetic wave velocity value.

[0046] The number of cycles for image entropy can be determined by Calculated; where V min V max Represents the preset minimum and maximum electromagnetic wave speeds; The value step size represents the preset electromagnetic wave speed; N represents the total number of iterations.

[0047] The minimum image entropy corresponding to the electromagnetic wave speed value can be obtained in Matlab software using the minimum value function; if two identical entropy values ​​are obtained, the larger wave speed value is taken as the accurate electromagnetic wave speed value.

[0048] Example 2

[0049] like Figure 1 As shown in Embodiment 2, this invention provides a method for accurately calculating the electromagnetic wave velocity in tunnel lining radar data, thereby solving the problem of accurate calculation of electromagnetic wave velocity in tunnel lining radar data. The method includes the following steps:

[0050] Step S1: Radar image offset processing based on the FK offset algorithm. Given a velocity value, the original radar image is offset using the FK offset algorithm.

[0051] Step S2: Calculate the entropy value of the offset radar image. Calculate the entropy value of the offset radar image corresponding to a given electromagnetic wave velocity value.

[0052] Step S3: Plot the velocity-entropy curve of the radar image. Based on the correspondence between velocity and image entropy, plot a velocity-entropy line graph.

[0053] Step S4: Obtain the precise electromagnetic wave velocity value. Based on the velocity-entropy line graph, the precise electromagnetic wave velocity value in the radar data is obtained by iteratively searching for the electromagnetic wave velocity value corresponding to the minimum image entropy value.

[0054] The offset process in step S1 can be completed in Reflexw software; and when the range and step size of the electromagnetic wave velocity are given, the initial electromagnetic wave velocity value can be set to 0.1 m / ns based on engineering experience; the electromagnetic wave velocity value is determined by... The calculation yields V, where n represents the number value of the nth iteration. n This represents the electromagnetic wave speed value calculated in the nth iteration.

[0055] When the FK offset algorithm in step S1 is insufficient in offsetting the hyperbola in the radar image, an 'under-offset' phenomenon occurs. Conversely, when the electromagnetic wave velocity is too high, the FK offset algorithm causes excessive hyperbola offset in the radar image, i.e., an 'over-offset' phenomenon occurs. Based on this characteristic, the optimal electromagnetic wave velocity value for achieving the best radar image offset effect can be obtained by comparing the offset effects of radar images under different electromagnetic wave velocities.

[0056] The one-dimensional scan signal in the image entropy value of step S2 can be obtained from... It means that x ij X represents the amplitude value of the sampling point; j=T Represents a column vector consisting of sampled points; the minimum image entropy value of the offset radar image can be expressed by the formula... Calculate, where X represents a matrix vector consisting of sampling points; E(X) represents the image entropy value of the offset image.

[0057] The smaller the image entropy value in the velocity-entropy curve of the radar image in step S3, the better the offset effect. Therefore, by iteratively calculating the entropy of the offset radar image, the electromagnetic wave velocity value corresponding to the minimum image entropy value can be used as the accurate electromagnetic wave velocity value.

[0058] The number of iterations for the image entropy value in step S4 can be determined by... Calculated. Where, V min V max Represents the preset minimum and maximum electromagnetic wave speeds; The value step size represents the preset electromagnetic wave speed; N represents the total number of iterations.

[0059] The minimum image entropy corresponding to the electromagnetic wave velocity value in step S4 can be obtained using the minimum value function in Matlab software. If two identical entropy values ​​are obtained, the larger wave velocity value is taken as the accurate electromagnetic wave velocity value.

[0060] In this embodiment, Figure 2This diagram illustrates the hyperbolic migration process for a discrete single target. The true position of target B is M, which is formed by electromagnetic wave diffraction, creating a hyperbola ABC. Taking sampling point A on the left side of the diffraction hyperbola as an example, it has a horizontal offset x and a vertical offset t relative to the true position B of the target. Therefore, a migration algorithm is needed to reposition the sampling point signal at position A to position B. This process is repeated until all sampling points on the 'tail' of the diffraction hyperbola are repositioned to position B, thus obtaining the true position B of the discrete entity target.

[0061] In this embodiment, the measured radar data line length for the reinforcing bars is 3.7m, containing 14 reinforcing bars. Based on engineering experience, the electromagnetic wave velocity value can be set to 0.11m / ns, the electromagnetic wave velocity step size is set to 0.01m / ns, and the value range is [0.01m / ns, 0.20m / ns]. When the electromagnetic wave velocity values ​​are 0.1m / ns and 0.12m / ns respectively, the minimum radar image entropy is 0.9491. Figure 3 The red dots are shown in the velocity-entropy line graph.

[0062] like Figure 4 As shown in -b, when the electromagnetic wave velocity is 0.12 m / ns, the radar image offset effect is the best and close to the empirical value of 0.11 m / ns, which is suitable for accurate calculation of electromagnetic wave velocity in rebar radar data. When the electromagnetic wave velocity is less than 0.12 m / ns, the radar image shows an 'under-offset' state, such as... Figure 4 -a is shown; conversely, when the electromagnetic wave velocity is greater than 0.12 m / ns, the radar image shows an 'over-offset' state, as shown in Figure 1. Figure 4 -c is shown.

[0063] The measured radar data for the steel arches has a survey line length of 8.2m, encompassing 16 steel arches. Based on engineering experience, the electromagnetic wave velocity can be set to 0.10 m / ns, with a step size of 0.01 m / ns, and a range of [0.01 m / ns, 0.15 m / ns]. When the electromagnetic wave velocity is 0.07 m / ns, the corresponding minimum radar image entropy is 0.999936. Figure 5 The red dots are shown in the velocity-entropy line graph.

[0064] like Figure 6 As shown in -b, the radar image offset effect is best when the electromagnetic wave velocity is 0.07 m / ns. The calculated electromagnetic wave velocity value is slightly lower than the empirical value, which may be related to the strong interference of the internal steel arch frame on the electromagnetic wave signal. When the electromagnetic wave velocity is 0.02 m / ns, the radar image shows an 'under-offset' characteristic, such as... Figure 6 -a is shown; conversely, when the electromagnetic wave velocity is 0.12 m / ns, the radar image shows an 'over-offset' characteristic, as shown in Figure 1. Figure 6 -c is shown.

[0065] Example 3

[0066] This embodiment 3 provides a non-transitory computer-readable storage medium for storing computer instructions. When executed by a processor, the computer instructions implement the electromagnetic wave velocity calculation method in tunnel lining radar data as described above. The method includes:

[0067] Given a velocity value, perform radar image offset processing based on the FK offset algorithm;

[0068] Calculate the entropy value of the radar image after offset corresponding to a given electromagnetic wave velocity value;

[0069] Based on the correspondence between velocity and image entropy, a velocity-entropy line graph is constructed;

[0070] Based on the velocity-entropy line graph, the precise electromagnetic wave velocity value in the radar data is obtained by iteratively searching for the electromagnetic wave velocity value corresponding to the minimum image entropy value.

[0071] Example 4

[0072] This embodiment 4 provides a computer device, including a memory and a processor, wherein the processor and the memory communicate with each other, and the memory stores program instructions that can be executed by the processor. The processor calls the program instructions to execute a method for calculating the electromagnetic wave velocity in tunnel lining radar data, the method including:

[0073] Given a velocity value, perform radar image offset processing based on the FK offset algorithm;

[0074] Calculate the entropy value of the radar image after offset corresponding to a given electromagnetic wave velocity value;

[0075] Based on the correspondence between velocity and image entropy, a velocity-entropy line graph is constructed;

[0076] Based on the velocity-entropy line graph, the precise electromagnetic wave velocity value in the radar data is obtained by iteratively searching for the electromagnetic wave velocity value corresponding to the minimum image entropy value.

[0077] Example 5

[0078] This embodiment 5 provides an electronic device, including: a processor, a memory, and a computer program; wherein, the processor is connected to the memory, and the computer program is stored in the memory. When the electronic device is running, the processor executes the computer program stored in the memory to cause the electronic device to execute instructions for implementing the method for calculating electromagnetic wave velocity in tunnel lining radar data as described above. The method includes:

[0079] Given a velocity value, perform radar image offset processing based on the FK offset algorithm;

[0080] Calculate the entropy value of the radar image after offset corresponding to a given electromagnetic wave velocity value;

[0081] Based on the correspondence between velocity and image entropy, a velocity-entropy line graph is constructed;

[0082] Based on the velocity-entropy line graph, the precise electromagnetic wave velocity value in the radar data is obtained by iteratively searching for the electromagnetic wave velocity value corresponding to the minimum image entropy value.

[0083] In summary, the method and system for calculating electromagnetic wave velocity in tunnel lining radar data described in this invention solve the problem of accurate calculation of electromagnetic wave velocity in tunnel lining radar data. The accurate electromagnetic wave velocity calculation method designed in this invention first requires offsetting the radar data using a frequency wavenumber domain offset algorithm. Then, by calculating the image entropy value of the offset radar data, a velocity-entropy curve of the radar image is plotted. Finally, by iteratively searching for the electromagnetic wave velocity value corresponding to the minimum image entropy value, the accurate electromagnetic wave velocity value in the radar data can be obtained. This invention, through its research on the accurate calculation method of electromagnetic wave velocity in tunnel lining radar data, overcomes the deficiency of accurate calculation of electromagnetic wave velocity in tunnel lining radar data by ground-penetrating radar, providing a theoretical basis and data support for subsequent research on electromagnetic wave velocity in tunnel lining radar data.

[0084] Those skilled in the art will understand that embodiments of the present invention can be provided as methods, systems, or computer program products. Therefore, the present invention can take the form of a completely hardware embodiment, a completely software embodiment, or an embodiment combining software and hardware aspects. Furthermore, the present invention can 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, etc.) containing computer-usable program code.

[0085] This invention is described with reference to flowchart illustrations and / or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the invention. It will be understood that each block of the flowchart illustrations and / or block diagrams, and combinations of blocks in the flowchart illustrations and / or block diagrams, can be implemented by computer program instructions. These computer program instructions can 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, generate instructions for implementing the flowchart illustrations and / or block diagrams. Figure 1 One or more processes and / or boxes Figure 1 A device that provides the functions specified in one or more boxes.

[0086] These computer program instructions may also be stored in a computer-readable storage medium that can direct a computer or other programmable data processing device to function in a particular manner, such that the instructions stored in the computer-readable storage medium produce an article of manufacture including instruction means, which are implemented in a process Figure 1 One or more processes and / or boxes Figure 1 The function specified in one or more boxes.

[0087] These computer program instructions may also be loaded onto a computer or other programmable data processing equipment, whereby a series of operational steps are performed to produce a computer-implemented process, thereby providing instructions that execute on the computer or other programmable equipment for implementing the process. Figure 1 One or more processes and / or boxes Figure 1 The steps of the function specified in one or more boxes.

[0088] While the specific embodiments of the present invention have been described above in conjunction with the accompanying drawings, this is not intended to limit the scope of protection of the present invention. Those skilled in the art should understand that, based on the technical solutions disclosed in the present invention, various modifications or variations that can be made by those skilled in the art without creative effort should be included within the scope of protection of the present invention.

Claims

1. A method for calculating electromagnetic wave velocity in radar data of tunnel lining, characterized in that, include: Given a velocity value, perform radar image offset processing based on the FK offset algorithm; Calculate the entropy value of the radar image after offset corresponding to a given electromagnetic wave velocity value; Based on the correspondence between velocity and image entropy, a velocity-entropy line graph is constructed; Based on the velocity-entropy line graph, the electromagnetic wave velocity value corresponding to the minimum image entropy value is obtained by iteratively searching for the precise electromagnetic wave velocity value in the radar data. Wherein, given the range and step size of the electromagnetic wave velocity, the initial electromagnetic wave velocity value is set to 0.1 m / ns; the electromagnetic wave velocity value is determined by... The calculation yielded, where Representing the The number value calculated in the next loop. Representing the The electromagnetic wave velocity value calculated in the next iteration; This represents the minimum preset electromagnetic wave speed; The step size represents the preset electromagnetic wave velocity. The one-dimensional scan signal in the image entropy value is composed of It means that among them The amplitude value represents the sampling point; Represents a column vector consisting of sampled points; the minimum image entropy value of the offset radar image is obtained through the formula... Calculate, where, Represents a matrix vector consisting of sampling points; The image entropy value represents the offset image.

2. The method for calculating electromagnetic wave velocity in tunnel lining radar data according to claim 1, characterized in that, The smaller the image entropy value in the velocity-entropy curve of the radar image, the better the offset effect. Therefore, by iteratively calculating the entropy of the offset radar image, the electromagnetic wave velocity value corresponding to the minimum image entropy value is taken as the accurate electromagnetic wave velocity value.

3. The method for calculating electromagnetic wave velocity in tunnel lining radar data according to claim 2, characterized in that, The number of cycles for image entropy can be determined by Calculated; where, , Represents the preset minimum and maximum electromagnetic wave speeds; The step size represents the preset electromagnetic wave velocity. This represents the total number of iterations.

4. The method for calculating electromagnetic wave velocity in tunnel lining radar data according to claim 2, characterized in that, The minimum image entropy corresponding to the electromagnetic wave speed value can be obtained in Matlab software using the minimum value function; if two identical entropy values ​​are obtained, the larger wave speed value is taken as the accurate electromagnetic wave speed value.

5. A system for calculating electromagnetic wave velocity in radar data of tunnel lining, used to implement the method as described in any one of claims 1-4, characterized in that, include: The offset module is used to perform radar image offset processing based on the FK offset algorithm, given a velocity value. The calculation module is used to calculate the entropy value of the radar image after offset corresponding to a given electromagnetic wave velocity value; A module is built to create a velocity-entropy line graph based on the relationship between velocity and image entropy. The cyclic search module is used to obtain the precise electromagnetic wave velocity value in the radar data by cyclically searching for the electromagnetic wave velocity value corresponding to the minimum image entropy value based on the velocity-entropy line graph.

6. A non-transitory computer-readable storage medium, characterized in that, The non-transitory computer-readable storage medium is used to store computer instructions, which, when executed by a processor, implement the method for calculating electromagnetic wave velocity in tunnel lining radar data as described in any one of claims 1-4.

7. A computer device, characterized in that, The system includes a memory and a processor, which communicate with each other. The memory stores program instructions that can be executed by the processor, and the processor calls the program instructions to execute the electromagnetic wave velocity calculation method in tunnel lining radar data as described in any one of claims 1-4.

8. An electronic device, characterized in that, include: The device includes a processor, a memory, and a computer program; wherein the processor is connected to the memory, the computer program is stored in the memory, and when the electronic device is running, the processor executes the computer program stored in the memory to cause the electronic device to execute instructions for implementing the method for calculating electromagnetic wave velocity in tunnel lining radar data as described in any one of claims 1-4.