A tunnel measuring point fitting method, system, device and storage medium

Abstract

The application discloses a tunnel measuring point fitting method, system, device and storage medium, considers the diversity of the tunnel design contour and the complexity of the tunnel construction scene, and therefore has higher universality for tunnel actual measurement fitting calculation. An intelligent interpolation method is adopted, the tunnel design contour is comprehensively used, trend information of surrounding measuring points is combined, and abnormal values are more reasonably estimated. Elimination of the abnormal values does not affect the order of non-abnormal values, interpolation only affects the abnormal value measuring points which are eliminated, does not change the positions and the order of other measuring points, the algorithm can process various abnormal value arrangement conditions, and has strong stability.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of tunnel engineering detection, more particularly, to a tunnel measuring point fitting method, system, device and storage medium. BACKGROUND

[0002] The development of measurement robots has greatly promoted the field of measurement, especially in complex environments and large-scale projects. In tunnel engineering, the development of measurement robots marks the transition of traditional measurement methods to automation and intelligence. Measurement robots can autonomously measure cross sections inside the tunnel, greatly improving measurement efficiency and data acquisition accuracy. However, due to the increase in data acquisition and the complexity of the tunnel measurement environment, how to effectively process and optimize these cross section measurement data has become an urgent problem.

[0003] Currently, some tunnel measuring point data processing algorithms have gradually appeared on the market. The existing technology is aimed at three-dimensional laser scanned point cloud data, and only noise point processing is performed on the data. The existing technology only removes noise points from three-dimensional laser scanned point cloud data, and limits the tunnel contour to be an ellipse or a circle when processing data. However, in actual tunnel design, not only independent ellipses and circles are designed, but most cases are combinations of arcs and straight lines. SUMMARY

[0004] To overcome the defects of low measurement efficiency and low data acquisition accuracy of the prior art, the present application provides a tunnel measuring point fitting method, system, device and storage medium.

[0005] To solve the above technical problems, the technical solutions of the present application are as follows:

[0006] The present application provides a tunnel measuring point fitting method, comprising:

[0007] S1: obtaining tunnel measuring point coordinate data set and tunnel contour line data;

[0008] S2: calculating overbreak value and underbreak value according to the tunnel measuring point coordinate data set, and determining average overbreak value and average underbreak value according to the overbreak value and the underbreak value;

[0009] S3: determining the offset value of the tunnel contour line based on the average overbreak value and the average underbreak value;

[0010] S4: determining abnormal measuring points for all tunnel measuring points according to a preset measuring point optimization threshold, obtaining a continuous abnormal measuring point set and a discrete abnormal measuring point set;

[0011] S5: offsetting the tunnel contour line as a whole according to the offset value of the tunnel contour line;

[0012] S6: Grouping the continuous abnormal measuring point set, performing fitting preprocessing based on each group of abnormal measuring points to obtain the state of each group of abnormal measuring points, and executing step S7 until the preset state is met;

[0013] S7: Connecting the adjacent normal points before and after each group of abnormal measuring points to form a straight line, and calculating the bisecting point on the straight line according to the number of the group of abnormal measuring points;

[0014] S8: Projecting the bisecting point on the straight line corresponding to each group of abnormal measuring points onto the overall offset tunnel contour line to fit, calculating the projection point coordinates, modifying the coordinates of each group of abnormal measuring points to the corresponding projection point coordinates, and processing the first and last abnormal measuring points of the tunnel measuring point coordinate data set to obtain the fitted tunnel measuring point coordinate data set of all abnormal measuring points.

[0015] Preferably, in step S3, the offset value of the tunnel contour line is determined based on the average overbreak value and the average underbreak value, which includes three cases:

[0016] S31: When the average overbreak value is not close to 0, the average overbreak value is taken as the offset value of the tunnel contour line;

[0017] S32: When the average overbreak value is close to 0, the average underbreak value is taken as the offset value of the tunnel contour line;

[0018] S33: When the average overbreak value and the average underbreak value are both close to 0, there is no abnormal measuring point to be fitted in this case.

[0019] Preferably, in step S5, the abnormal measuring point determination is performed on all tunnel measuring points according to the preset measuring point optimization threshold, which includes:

[0020] If the absolute values of the overbreak value and the underbreak value of the tunnel measuring point are greater than the preset measuring point optimization threshold, the tunnel measuring point is determined as an abnormal measuring point; if the absolute values of the overbreak value and the underbreak value of the tunnel measuring point are less than or equal to the preset measuring point optimization threshold, the tunnel measuring point is a normal measuring point.

[0021] Preferably, in step S6, the continuous abnormal measuring point set is grouped, and fitting preprocessing is performed based on each group of abnormal measuring points to obtain the state of each group of abnormal measuring points, the fitting preprocessing includes four cases, and step S7 is executed until the preset state is met, which includes:

[0022] S61: When the first measuring point in the tunnel measuring point coordinate data set is an abnormal measuring point, and the index value of the first abnormal measuring point in the first group in the continuous abnormal measuring point set is 0, and the last abnormal measuring point in the group is not the last measuring point in the tunnel measuring point coordinate data set, the projection point coordinates of the first abnormal measuring point in the group on the tunnel profile line after the overall offset are calculated, the abnormal measuring point coordinates are recorded, and the coordinates of the abnormal measuring point are modified to the projection point coordinates; if there is only one abnormal measuring point in the group, the fitting pretreatment is continued, otherwise the coordinates of the normal measuring point adjacent to the last abnormal measuring point in the group are recorded, and step S7 is performed;

[0023] S62: When the last measuring point in the tunnel measuring point coordinate data set is an abnormal measuring point, and the last abnormal measuring point is the last measuring point in the tunnel measuring point coordinate data set, the projection point coordinates of the last abnormal measuring point in the group on the tunnel profile line after the overall offset are calculated, the abnormal measuring point coordinates are recorded, and the coordinates of the abnormal measuring point are modified to the projection point coordinates; if there is only one abnormal measuring point in the group, the fitting pretreatment is continued, otherwise the coordinates of the normal measuring point adjacent to the first abnormal measuring point in the group are recorded, and step S7 is performed;

[0024] S63: When the first measuring point and the last measuring point in the tunnel measuring point coordinate data set are abnormal measuring points, the tunnel measuring point coordinate data set has and only has one group of abnormal measuring points; the projection point coordinates of the first abnormal measuring point in the group on the tunnel profile line after the overall offset are calculated, the abnormal measuring point coordinates are recorded, and the coordinates of the abnormal measuring point are modified to the projection point coordinates; if there is only one abnormal measuring point in the group, the fitting pretreatment is continued, otherwise the coordinates of the normal measuring point adjacent to the last abnormal measuring point in the group are recorded, and a judgment is made as to whether there are two abnormal measuring points; if there are only two abnormal measuring points in the group, the fitting pretreatment is continued, otherwise the abnormal measuring point coordinates are recorded and the coordinates of the abnormal measuring point are modified to the projection point coordinates, the first abnormal measuring point and the last abnormal measuring point in the group are taken as normal measuring points, and step S7 is performed;

[0025] S64: When the first measuring point and the last measuring point in the tunnel measuring point coordinate data set are normal measuring points, the coordinates of the previous normal measuring point of the first abnormal measuring point and the coordinates of the next normal measuring point of the last abnormal measuring point in the group of abnormal measuring points are recorded, and S7 is performed.

[0026] Preferably, in step S8, the bisecting point on the straight line corresponding to each group of abnormal measuring points is projected onto the tunnel profile line after the overall offset, the projection point coordinates are calculated, the coordinates of each group of abnormal measuring points are modified to the corresponding projection point coordinates, and after the first and last abnormal measuring points of the tunnel measuring point coordinate data set are processed, the tunnel measuring point coordinate data set after fitting of all abnormal measuring points is obtained, including:

[0027] S81: Projecting the bisecting point on the straight line corresponding to each group of abnormal measuring points to the fitted tunnel contour line after the overall offset, calculating the coordinates of the projected points, and modifying the coordinates of each group of abnormal measuring points to the corresponding projected point coordinates;

[0028] S82: Processing the first and last abnormal measuring points in the tunnel measuring point coordinate data set, if the number of tunnel measuring points in the tunnel measuring point coordinate data set is less than or equal to three, then intelligent interpolation of tunnel measuring points is performed; otherwise, the first and last abnormal measuring points are processed by vector inner product;

[0029] Preferably, in step S82, if the number of tunnel measuring points in the tunnel measuring point coordinate data set is less than or equal to three, then intelligent interpolation of tunnel measuring points is performed, specifically as follows:

[0030] S821: If the number of tunnel measuring points in the tunnel measuring point coordinate data set is less than or equal to three, then the average interval is calculated according to the distance between adjacent measuring points in the tunnel measuring point coordinate set;

[0031] S822: Calculate the interval between the first and last abnormal measuring points and the adjacent measuring points in the tunnel measuring point coordinate set, if the interval between the first abnormal measuring point and the adjacent measuring point is greater than twice the average interval, then an interpolation point is inserted between the first abnormal measuring point and the adjacent measuring point according to the average interval, and the interpolation point is projected onto the tunnel contour line after the offset, and the coordinates of the first and last abnormal measuring points are modified to the corresponding projected point coordinates.

[0032] Preferably, in step S82, the first and last abnormal measuring points are processed by vector inner product, specifically as follows:

[0033] S823: If the number of tunnel measuring points in the tunnel measuring point coordinate data set is greater than three, then vectors are constructed, including a first vector, a second vector, a third vector, and a fourth vector; the direction of the first vector is from the first tunnel measuring point to the second tunnel measuring point; the direction of the second vector is from the second tunnel measuring point to the third tunnel measuring point; the direction of the third vector is from the last tunnel measuring point to the second last tunnel measuring point; and the direction of the fourth vector is from the second last tunnel measuring point to the third last tunnel measuring point;

[0034] S824: Calculate the vector inner product of the first vector and the second vector, if the vector inner product is less than 0, then remove the first tunnel measuring point; calculate the third vector and the fourth vector, if the vector inner product is less than 0, then remove the last tunnel measuring point.

[0035] The application also provides a tunnel measuring point fitting system for implementing the above method, comprising:

[0036] a data acquisition module for acquiring tunnel measuring point coordinate data set and tunnel contour line data;

[0037] The overbreak value calculation module calculates the overbreak value and the underbreak value according to the tunnel point coordinate data set, and determines the average overbreak value and the average underbreak value according to the overbreak value and the underbreak value;

[0038] The offset value acquisition module determines the offset value of the tunnel contour line based on the average overbreak value and the average underbreak value.

[0039] The abnormal point determination module determines all tunnel points as abnormal points according to a preset point optimization threshold, and obtains a continuous abnormal point set and a discrete abnormal point set.

[0040] The offset module offsets the tunnel contour line as a whole according to the offset value of the tunnel contour line.

[0041] The abnormal point state determination module groups the continuous abnormal point set, performs fitting preprocessing based on each group of abnormal points, and obtains the state of each group of abnormal points until the preset state is met, and then jumps to the bisecting point calculation module.

[0042] The bisecting point calculation module connects adjacent normal points before and after each group of abnormal points to form a straight line, and calculates a bisecting point on the straight line according to the number of the group of abnormal points.

[0043] The point fitting module projects the bisecting point on the straight line corresponding to each group of abnormal points onto the tunnel contour line after overall offsetting, fits the projection point coordinates, modifies the coordinates of each group of abnormal points to the corresponding projection point coordinates, and processes the first and last abnormal points of the tunnel point coordinate data set, to obtain the tunnel point coordinate data set after fitting of all abnormal points.

[0044] The present application also provides an electronic device comprising:

[0045] One or more processors;

[0046] Memory for storing one or more programs;

[0047] When the one or more programs are executed by the one or more processors, the processor implements the above method.

[0048] The present application also provides a computer-readable storage medium comprising a computer program, which stores a computer program, characterized in that the program is executed by one or more processors to implement the above method.

[0049] Compared with the prior art, the beneficial effects of the technical scheme of the present application are:

[0050] The application provides a tunnel measuring point fitting method, system, device and storage medium, which considers the diversity of tunnel design profiles and the complexity of tunnel construction scenes, and thus has high universality for tunnel actual measurement fitting calculation. An intelligent interpolation method is adopted, tunnel design profiles are comprehensively considered, trend information of surrounding measuring points is combined, and abnormal values are more reasonably estimated. The abnormal values are removed without affecting the order of non-abnormal values, the interpolation only affects the removed abnormal value measuring points, does not change the positions and order of other measuring points, the algorithm can process various abnormal value arrangement conditions, and has high stability. BRIEF DESCRIPTION OF DRAWINGS

[0051] Figure 1 A flowchart of a tunnel measuring point fitting method described in Embodiment 1 is provided.

[0052] Figure 2 A flowchart of a tunnel measuring point fitting method described in Embodiment 2 is provided.

[0053] Figure 3 A structural schematic diagram of a tunnel measuring point fitting system described in Embodiment 3 is provided. DETAILED DESCRIPTION

[0054] The drawings are only used for illustrative description, and cannot be understood as limitation to the patent;

[0055] In order to better illustrate the embodiments, some components in the drawings may be omitted, enlarged or reduced, and do not represent the actual product size;

[0056] For those skilled in the art, it is understandable that some known structures and their descriptions in the drawings may be omitted.

[0057] The technical solutions of the application will be further described below in combination with the drawings and embodiments.

[0058] Embodiment 1

[0059] The embodiment provides a tunnel measuring point fitting method, as shown in the figure, which comprises: Figure 1

[0060] S1: Obtain tunnel measuring point coordinate data set and tunnel profile line data;

[0061] S2: Calculate overbreak value and underbreak value according to the tunnel measuring point coordinate data set, and determine average overbreak value and average underbreak value according to the overbreak value and the underbreak value;

[0062] S3: Determine the offset value of the tunnel profile line based on the average overbreak value and the average underbreak value;

[0063] S4: Determine abnormal measuring points according to a preset measuring point optimization threshold for all tunnel measuring points, and obtain a continuous abnormal measuring point set and a discrete abnormal measuring point set; ​

[0064] S5: overall shift the tunnel contour line according to the offset value of the tunnel contour line;

[0065] S6: group the continuous abnormal measuring point set, perform fitting preprocessing based on each group of abnormal measuring points, obtain the state of each group of abnormal measuring points, and until the preset state is met, perform step S7;

[0066] S7: connect the adjacent normal points before and after each group of abnormal measuring points to form a straight line, and calculate the bisecting point on the straight line according to the number of the group of abnormal measuring points;

[0067] S8: project the bisecting point on the straight line corresponding to each group of abnormal measuring points onto the overall shifted tunnel contour line for fitting, calculate the projection point coordinates, modify the coordinates of each group of abnormal measuring points to the corresponding projection point coordinates, and after processing the first and last abnormal measuring points of the tunnel measuring point coordinate data set, obtain the tunnel measuring point coordinate data set after fitting of all abnormal measuring points.

[0068] The embodiment considers the diversity of tunnel design contours and the complexity of tunnel construction scenes, and therefore has high universality for tunnel measurement fitting calculation. The intelligent interpolation method is adopted to comprehensively consider the tunnel design contour and the trend information of the surrounding measuring points to more reasonably estimate the abnormal values. Removing abnormal values does not affect the order of non-abnormal values, interpolation only affects the removed abnormal value measuring points, does not change the positions and orders of other measuring points, the algorithm can handle various abnormal value arrangement situations, and has high stability.

[0069] Embodiment 2

[0070] The embodiment provides a tunnel measuring point fitting method, comprising:

[0071] S1: obtaining a tunnel measuring point coordinate data set and a tunnel contour line data;

[0072] S2: calculating overbreak values and underbreak values according to the tunnel measuring point coordinate data set, and determining average overbreak values and average underbreak values according to the overbreak values and the underbreak values;

[0073] S3: determining an offset value of the tunnel contour line based on the average overbreak values and the average underbreak values;

[0074] In step S3, the offset value of the tunnel contour line is determined based on the average overbreak values and the average underbreak values, which includes three cases:

[0075] S31: when the average overbreak value is not close to 0, the average overbreak value is taken as the offset value of the tunnel contour line;

[0076] S32: when the average overbreak value is close to 0, the average underbreak value is taken as the offset value of the tunnel contour line;

[0077] S33: When the average overbreak value and the average underbreak value are close to 0, there is no abnormal point to be fitted in this case.

[0078] S4: Abnormal point determination is performed on all tunnel points according to a preset point optimization threshold to obtain a continuous abnormal point set and a discrete abnormal point set.

[0079] S5: The overall offset of the tunnel contour line is performed according to the offset value of the tunnel contour line.

[0080] In step S5, abnormal point determination is performed on all tunnel points according to a preset point optimization threshold, including:

[0081] If the absolute values of the overbreak value and the underbreak value of the tunnel point are greater than the preset point optimization threshold, the tunnel point is determined as an abnormal point; if the absolute values of the overbreak value and the underbreak value of the tunnel point are less than or equal to the preset point optimization threshold, the tunnel point is a normal point.

[0082] S6: The continuous abnormal point set is grouped, and fitting preprocessing is performed based on each group of abnormal points to obtain the state of each group of abnormal points, until the preset state is met, and step S7 is performed.

[0083] In step S6, the continuous abnormal point set is grouped, and fitting preprocessing is performed based on each group of abnormal points to obtain the state of each group of abnormal points. The fitting preprocessing is divided into four cases, until the preset state is met, and step S7 is performed, including:

[0084] S61: When the first point in the tunnel point coordinate data set is an abnormal point, and the index value of the first abnormal point in the first group in the continuous abnormal point set is 0, and the last abnormal point in the group is not the last point in the tunnel point coordinate data set, the projection point coordinates of the first abnormal point in the group are calculated, the abnormal point coordinates are recorded, and the coordinates of the abnormal point are modified to the projection point coordinates. If there is only one abnormal point in the group, fitting preprocessing is continued, otherwise the coordinates of the normal point adjacent to the last abnormal point in the group are recorded, and step S7 is performed.

[0085] S62: When the last point in the tunnel point coordinate data set is an abnormal point, and the last abnormal point is the last point in the tunnel point coordinate data set, the projection point coordinates of the last abnormal point in the group are calculated, the abnormal point coordinates are recorded, and the coordinates of the abnormal point are modified to the projection point coordinates. If there is only one abnormal point in the group, fitting preprocessing is continued, otherwise the coordinates of the normal point adjacent to the first abnormal point in the group are recorded, and step S7 is performed.

[0086] S63: When the first and last points in the tunnel point coordinate data set are both abnormal points, the tunnel point coordinate data set has and only has one group of abnormal points; the first abnormal point in the group is projected onto the overall offset tunnel profile line to calculate the projection point coordinates, the abnormal point coordinates are recorded and modified to the projection point coordinates; if there is only one abnormal point in the group, the fitting pretreatment is continued, otherwise the coordinates of the normal points adjacent to the last abnormal point in the group are recorded, and whether there are two abnormal points is judged; if there are only two abnormal points in the group, the fitting pretreatment is continued, otherwise the abnormal point coordinates are recorded and modified to the projection point coordinates, the first and last abnormal points in the group are taken as normal points, and step S7 is executed;

[0087] S64: When the first and last points in the tunnel point coordinate data set are both normal points, the coordinates of the previous normal point of the first abnormal point and the next normal point of the last abnormal point in the group of abnormal points are recorded, and S7 is executed.

[0088] S7: Connecting the adjacent normal points before and after each group of abnormal points to form a straight line, and calculating the bisector point on the straight line according to the number of abnormal points in the group;

[0089] S8: Projecting the bisector point on the straight line corresponding to each group of abnormal points onto the overall offset tunnel profile line to fit, calculating the projection point coordinates, modifying the coordinates of each group of abnormal points to the corresponding projection point coordinates, and after processing the first and last abnormal points of the tunnel point coordinate data set, obtaining the tunnel point coordinate data set after fitting of all abnormal points.

[0090] In step S8, the bisector point on the straight line corresponding to each group of abnormal points is projected onto the overall offset tunnel profile line to fit, the projection point coordinates are calculated, the coordinates of each group of abnormal points are modified to the corresponding projection point coordinates, and after processing the first and last abnormal points of the tunnel point coordinate data set, the tunnel point coordinate data set after fitting of all abnormal points is obtained, including:

[0091] S81: Projecting the bisector point on the straight line corresponding to each group of abnormal points onto the overall offset tunnel profile line to fit, calculating the projection point coordinates, and modifying the coordinates of each group of abnormal points to the corresponding projection point coordinates;

[0092] S82: Processing the first and last abnormal points of the tunnel point coordinate data set, if the number of tunnel points in the tunnel point coordinate data set is less than or equal to three, intelligent interpolation of tunnel points is performed; otherwise, the first and last abnormal points are processed by vector inner product;

[0093] In step S82, if the number of tunnel survey points in the tunnel survey point coordinate data set is less than or equal to three, intelligent interpolation of the tunnel survey points is performed, specifically as follows:

[0094] S821: If the number of tunnel survey points in the tunnel survey point coordinate data set is less than or equal to three, the average interval is calculated according to the distance between adjacent survey points in the tunnel survey point coordinate set.

[0095] S822: The interval between the first and last abnormal survey points and the adjacent survey points in the tunnel survey point coordinate set is calculated. If the interval between the first abnormal survey point and the adjacent survey point is greater than twice the average interval, an interpolated point is interpolated between the first abnormal survey point and the adjacent survey point according to the average interval, and the first and last abnormal survey points are projected onto the offset tunnel contour line, and the coordinates of the first and last abnormal survey points are modified to the coordinates of the corresponding projection points.

[0096] In step S82, the first and last abnormal survey points are processed by vector inner product, specifically as follows:

[0097] S823: If the number of tunnel survey points in the tunnel survey point coordinate data set is greater than three, vectors are constructed, including a first vector, a second vector, a third vector, and a fourth vector. The direction of the first vector is from the first tunnel survey point to the second tunnel survey point. The direction of the second vector is from the second tunnel survey point to the third tunnel survey point. The direction of the third vector is from the last but one tunnel survey point to the last but two tunnel survey point. The direction of the fourth vector is from the last but two tunnel survey point to the last tunnel survey point.

[0098] S824: The vector inner product of the first and second vectors is calculated. If the vector inner product is less than 0, the first tunnel survey point is removed. The third and fourth vectors are calculated. If the vector inner product is less than 0, the last tunnel survey point is removed.

[0099] In a specific embodiment, before the tunnel survey point fitting process is performed, the input and output of the data are introduced. The input of the data includes a survey point data point set, tunnel contour line data, and a survey point optimization threshold. The survey point data point set is a point set converted to the tunnel coordinate system from the tunnel section point set measured by the measurement robot. The tunnel coordinate system is a plane rectangular coordinate system XOY established with the tunnel contour design surface as the plane. The tunnel contour line data is the design data of the tunnel, which is generally considered to be a contour line composed of circular curves and straight line elements. Since the tunnel is mostly symmetrical, only the line element data on the right side of the tunnel contour line needs to be input. The circular curve in the line element data needs to input the center, the radius, and the terminal point coordinate. When the radius is 0, it represents a straight line element. The survey point optimization threshold is an index for determining whether the tunnel survey point is an abnormal point and whether it participates in subsequent fitting optimization. In this method, the overbreak and underbreak values of the survey point are compared. The output of the data is the survey point data point set after optimization fitting.

[0100] As Figure 2 shown, the process of tunnel point fitting processing is divided into the following 9 steps: S1 traverses the point set to calculate overbreak and underbreak, and calculates the average overbreak value and the average underbreak value; S2 calculates the tunnel contour line offset value; S3 groups the continuous abnormal points that exceed the threshold value; S4 offsets the entire tunnel contour line according to the tunnel contour line offset value; S5 traverses the abnormal point grouping set, and pre-processes the fitting; S6 connects the adjacent normal points before and after the abnormal points, and calculates the equidivision points; S7 point fitting, the equidivision points are projected to the offset tunnel, and the projection points are calculated; S8. Process the first and last abnormal points with vector inner product; S9 intelligent interpolation of points. The following is a detailed description of the above steps.

[0101] S1 traverses the point set to calculate overbreak and underbreak, and calculates the average overbreak value and the average underbreak value. The calculated overbreak and underbreak values are positive for overbreak, negative for underbreak, and 0 for the tunnel contour line. The average overbreak value is the average of all overbreak values, and the average underbreak value is the average of all underbreak values.

[0102] S2 calculates the tunnel contour line offset value. The offset value of the tunnel contour line is the value of the subsequent tunnel to be expanded or contracted, and a positive number indicates an external opening, and a negative number represents internal contraction. Because there is a situation of "prefer overbreak to underbreak" in tunnel excavation, in most cases, overbreak is matched with the actual situation, and the average overbreak value is used as the tunnel line offset value. The calculation is divided into three cases:

[0103] (S21) The average overbreak value is not close to 0, and the average overbreak value is used as the offset value.

[0104] (S22) The average overbreak value is close to 0, and the average underbreak value is used as the offset value.

[0105] (S23) The average overbreak value and the underbreak value are both close to 0, and there is no abnormal point to be fitted in this case, and the process directly jumps to S9 intelligent interpolation of points.

[0106] S3 groups the continuous abnormal points that exceed the threshold value. The tunnel point data point set is a set of continuous points, and the points may be continuously gathered or appear alone, which forms an alternating interlaced condition with non-abnormal points. The continuously appearing abnormal points are grouped as a group, and the single abnormal points are also grouped as a group. Traverse the entire tunnel point data point set to obtain the grouped abnormal point grouping set, wherein the grouping of abnormal points only records the index of the abnormal point. When grouping the abnormal points, the set point optimization threshold (positive number) is used for judgment. If the absolute value of the overbreak and underbreak value is greater than the set point optimization threshold, it is marked as an abnormal point.

[0107] S4 offset the tunnel profile line as a whole according to the tunnel profile line offset value. According to the interval of 1 millimeter, the points on the tunnel profile line are obtained, and the single-side buffer zone is established by connecting these points according to the tunnel line offset value calculated in S2, that is, the offset tunnel profile line is obtained. The offset tunnel profile line will be used as the reference line for subsequent fitting.

[0108] S5 traverse the abnormal point group set and perform fitting preprocessing. Determine the state of each group of abnormal points, and process abnormal points in different cases to prepare for the calculation and processing of S6. The fitting preprocessing includes the following four cases.

[0109] (S51) The starting point of the tunnel point data set is an abnormal point, the index of the first point in the first group in the abnormal point group set is 0, and the last point is not the last point of the tunnel point data set. Then this case can be determined. Calculate the projection point (the nearest point) of the first point in the abnormal group to the tunnel profile offset line, record and modify the coordinates of the measuring point. Determine whether there is only one element in the group, if yes, proceed to the next loop of S5, otherwise record the coordinates of the adjacent non-abnormal measuring points of the last element in the group, and execute S6 processing.

[0110] (S52) The ending point of the tunnel point data set is an abnormal point, the index of the last point in the last group in the abnormal point group set is recorded as the last point of the tunnel point data set. Then this case can be determined. Calculate the projection point (the nearest point) of the last point in the abnormal group to the tunnel profile offset line, record and modify the coordinates of the measuring point. Determine whether there is only one element in the group, if yes, proceed to the next loop of S5, otherwise record the coordinates of the adjacent non-abnormal measuring points of the first element in the group, and execute S6 processing.

[0111] (S53) The starting point and the ending point of the tunnel point data set are abnormal points, and there is only one group of abnormal points in this case, which contains all the measuring points. This case can be determined. Calculate the projection point (the nearest point) of the first point in the abnormal group to the tunnel profile offset line, record and modify the coordinates of the measuring point. Determine whether there is only one element in the group, if yes, proceed to the next loop of S5, otherwise calculate the projection point (the nearest point) of the last point in the abnormal group to the tunnel profile offset line; determine whether there are only two measuring points in the group, if yes, proceed to the next loop of S5, otherwise record and modify the coordinates of the measuring point, record the first element and the last element in the group as non-abnormal measuring point coordinates, and execute S6 processing.

[0112] (S54) The general case other than S51, S52, and S53, record the measuring point coordinates of the previous non-abnormal point of the first point in the abnormal group, record the measuring point coordinates of the next non-abnormal point of the last point in the abnormal group, and execute S6 processing.

[0113] S6 connecting the normal points adjacent to the abnormal point before and after the abnormal point to form a straight line, calculating the bisector point on the straight line according to the number of abnormal points, and performing S7 processing.

[0114] S7 measuring point fitting, projecting the bisector point to the offset tunnel, and calculating the projection point. Traversing the current abnormal grouping point, projecting the abnormal grouping point to the offset tunnel contour line, and modifying the point coordinates to the coordinates of the projection point.

[0115] S8. Processing the first and last abnormal points with vector inner product. If the number of tunnel measuring point data points is less than or equal to three, perform S9 processing. Otherwise, calculate the vectors, the direction of vector 1 is the first measuring point pointing to the second measuring point, the direction of vector 2 is the second measuring point pointing to the third measuring point, the direction of vector 3 is the first measuring point pointing to the second last measuring point, and the direction of vector 4 is the second last measuring point pointing to the third last measuring point. Calculate the vector inner product of vectors 1 and 2, and if the vector inner product is less than zero, remove the first measuring point. Calculate the vector inner product of vectors 3 and 4, and if the vector inner product is less than zero, remove the last measuring point.

[0116] S9 intelligent interpolation of measuring points, traversing the distance between adjacent measuring points, and calculating the average distance; traversing the distance between adjacent measuring points, and cyclically judging if the distance between adjacent measuring points is greater than twice the average distance. If so, a point is inserted between the two measuring points according to the difference between the average distance, and the inserted point is projected onto the offset tunnel contour line, and the coordinates of the projected point are saved.

[0117] Finally, the calculation is ended and the optimized measuring point data point set is returned.

[0118] Example 3

[0119] The embodiment also provides a tunnel measuring point fitting system for implementing the method of the embodiment 1 or the embodiment 2, as shown in the accompanying drawings, comprising: Figure 3

[0120] A data acquisition module acquires tunnel measuring point coordinate data set and tunnel contour line data;

[0121] An overbreak value and underbreak value calculation module calculates overbreak value and underbreak value according to the tunnel measuring point coordinate data set, and determines average overbreak value and average underbreak value according to the overbreak value and the underbreak value;

[0122] An offset value acquisition module determines the offset value of the tunnel contour line based on the average overbreak value and the average underbreak value;

[0123] An abnormal measuring point determination module determines all tunnel measuring points as abnormal measuring points according to a preset measuring point optimization threshold, and obtains continuous abnormal measuring point set and discrete abnormal measuring point set;

[0124] An offset module offsets the tunnel contour line as a whole according to the offset value of the tunnel contour line; ​

[0125] The abnormal measuring point state determination module groups the continuous abnormal measuring point set, performs fitting pretreatment based on each group of abnormal measuring points, obtains the state of each group of abnormal measuring points, and jumps to the bisecting point calculation module when the state meets a preset state;

[0126] The bisecting point calculation module connects the adjacent normal points before and after each group of abnormal measuring points to form a straight line, and calculates a bisecting point on the straight line according to the number of the group of abnormal measuring points.

[0127] The measuring point fitting module projects the bisecting point on the straight line corresponding to each group of abnormal measuring points onto the tunnel profile line after offset to fit, calculates the coordinates of the projection points, modifies the coordinates of each group of abnormal measuring points to the corresponding projection point coordinates, and processes the first and last abnormal measuring points of the tunnel measuring point coordinate data set to obtain the tunnel measuring point coordinate data set after fitting of all abnormal measuring points.

[0128] The embodiment also provides an electronic device, comprising:

[0129] One or more processors;

[0130] Memory for storing one or more programs;

[0131] When the one or more programs are executed by the one or more processors, the processor implements the method described above.

[0132] The embodiment also provides a computer readable storage medium containing a computer program, and the computer program is stored on the computer readable storage medium, and the program is executed by one or more processors to implement the method described above.

[0133] The same or similar reference signs correspond to the same or similar components;

[0134] The terms describing the positional relationship in the drawings are only used for exemplary illustration, and cannot be understood as a limitation on the patent;

[0135] Obviously, the above embodiments of the application are only examples for clearly illustrating the application, and are not intended to limit the implementation modes of the application. Based on the above description, other different forms of changes or variations can be made by those skilled in the art. Here, all the implementation modes are not required or can not be exhausted. Any modification, equivalent replacement and improvement made within the spirit and principle of the application should be included in the protection scope of the claims of the application.

Claims

1. A method of tunnel point fitting, the method comprising: The method comprises the following steps: S1: obtaining tunnel survey point coordinate data set and tunnel profile data; S2: calculating overbreak value and underbreak value according to the tunnel survey point coordinate data set, and determining average overbreak value and average underbreak value according to the overbreak value and the underbreak value; S3: determining the offset value of the tunnel profile based on the average overbreak value and the average underbreak value; S4: determining abnormal survey points according to a preset survey point optimization threshold, and obtaining continuous abnormal survey point set and discrete abnormal survey point set; S5: offsetting the tunnel profile based on the offset value of the tunnel profile; S6: grouping the continuous abnormal survey point set, performing fitting preprocessing based on each group of abnormal survey points, and obtaining the state of each group of abnormal survey points until a preset state is met, and then performing step S7; S7: connecting adjacent normal points before and after each group of abnormal survey points to form a straight line, and calculating the equidivision point on the straight line according to the number of the group of abnormal survey points; S8: projecting the equidivision point on the straight line corresponding to each group of abnormal survey points onto the offset tunnel profile, fitting the projection point coordinates, modifying the coordinates of each group of abnormal survey points to the corresponding projection point coordinates, and processing the first and last abnormal survey points of the tunnel survey point coordinate data set to obtain the tunnel survey point coordinate data set after fitting of all abnormal survey points.

2. The tunnel site fitting method of claim 1, wherein, In step S3, the offset value of the tunnel profile is determined based on the average overbreak value and the average underbreak value, which includes three cases: S31: when the average overbreak value is not close to 0, the average overbreak value is taken as the offset value of the tunnel profile; S32: when the average overbreak value is close to 0, the average underbreak value is taken as the offset value of the tunnel profile; S33: when the average overbreak value and the average underbreak value are both close to 0, there is no abnormal survey point to be fitted in this case.

3. The tunnel site fitting method of claim 1, wherein, In step S5, abnormal survey points are determined according to a preset survey point optimization threshold, which includes: If the absolute values of the overbreak value and the underbreak value of the tunnel survey point are greater than the preset survey point optimization threshold, the tunnel survey point is determined as an abnormal survey point; if the absolute values of the overbreak value and the underbreak value of the tunnel survey point are less than or equal to the preset survey point optimization threshold, the tunnel survey point is a normal survey point.

4. The tunnel site fitting method of claim 3, wherein, In step S6, the continuous abnormal survey point set is grouped, fitting preprocessing is performed based on each group of abnormal survey points, and the state of each group of abnormal survey points is obtained, which includes four cases until a preset state is met, and then step S7 is performed, which includes: S61: when the first survey point in the tunnel survey point coordinate data set is an abnormal survey point, the index value of the first abnormal survey point in the first group in the continuous abnormal survey point set is 0, and the last abnormal survey point in the group is not the last survey point in the tunnel survey point coordinate data set, the projection point coordinates of the first abnormal survey point in the group are calculated, the abnormal survey point coordinates are recorded, and the abnormal survey point coordinates are modified to the projection point coordinates; if there is only one abnormal survey point in the group, the fitting preprocessing is continued, otherwise the coordinates of the adjacent normal survey points of the last abnormal survey point in the group are recorded, and step S7 is performed; S62: When the last point in the tunnel point coordinate dataset is an abnormal point, and the last abnormal point is the last point in the tunnel point coordinate dataset, the projection point coordinate of the last abnormal point in the group of abnormal points is calculated, the coordinate of the abnormal point is recorded and modified to the projection point coordinate; if there is only one abnormal point in the group, the fitting preprocessing is continued, otherwise the coordinates of the normal points adjacent to the first abnormal point in the group are recorded, and step S7 is performed; S63: When the first point and the last point in the tunnel point coordinate dataset are abnormal points, the tunnel point coordinate dataset has and only has one group of abnormal points; the projection point coordinate of the first abnormal point in the group is calculated, the coordinate of the abnormal point is recorded and modified to the projection point coordinate; if there is only one abnormal point in the group, the fitting preprocessing is continued, otherwise the coordinates of the normal points adjacent to the last abnormal point in the group are recorded, and whether there are two abnormal points is judged; if there are only two abnormal points in the group, the fitting preprocessing is continued, otherwise the coordinates of the abnormal points are recorded and modified to the projection point coordinates, the first abnormal point and the last abnormal point in the group are taken as normal points, and step S7 is performed; S64: When the first point and the last point in the tunnel point coordinate dataset are normal points, the coordinates of the normal points adjacent to the first abnormal point and the last abnormal point in the group of abnormal points are recorded, and S7 is performed.

5. The tunnel site fitting method of claim 1, wherein, In step S8, the bisecting point on the straight line corresponding to each group of abnormal points is projected onto the tunnel profile line after the overall offset is fitted, the projection point coordinate is calculated, the coordinates of each group of abnormal points are modified to the corresponding projection point coordinates, and after the first and last abnormal points in the tunnel point coordinate dataset are processed, the tunnel point coordinate dataset after fitting of all abnormal points is obtained, including: S81: The bisecting point on the straight line corresponding to each group of abnormal points is projected onto the tunnel profile line after the overall offset is fitted, the projection point coordinate is calculated, and the coordinates of each group of abnormal points are modified to the corresponding projection point coordinates; S82: The first and last abnormal points in the tunnel point coordinate dataset are processed, if the number of tunnel points in the tunnel point coordinate dataset is less than or equal to three, intelligent interpolation of tunnel points is performed; otherwise, the first and last abnormal points are processed by vector inner product.

6. The tunnel site fitting method of claim 5, wherein, In step S82, if the number of tunnel points in the tunnel point coordinate dataset is less than or equal to three, intelligent interpolation of tunnel points is performed, which is specifically: S821: If the number of tunnel points in the tunnel point coordinate dataset is less than or equal to three, the average interval is calculated according to the distance between the adjacent points in the tunnel point coordinate set; S822: Calculate the distance between the first and last abnormal points and the adjacent points in the tunnel point coordinate set. If the distance between the first abnormal point and the adjacent point is greater than twice the average distance, an interpolation point is inserted between the first abnormal point and the adjacent point according to the average distance, the interpolation point is projected onto the offset tunnel contour line, and the coordinates of the first and last abnormal points are modified to the corresponding projection point coordinates.

7. The tunnel site fitting method of claim 5, wherein, In step S82, the first and last abnormal points are processed by vector inner product, which is specifically: S823: If the number of tunnel points in the tunnel point coordinate data set is greater than three, construct vectors including a first vector, a second vector, a third vector, and a fourth vector. The direction of the first vector is from the first tunnel point to the second tunnel point. The direction of the second vector is from the second tunnel point to the third tunnel point. The direction of the third vector is from the last tunnel point to the second last tunnel point. The direction of the fourth vector is from the second last tunnel point to the third last tunnel point. S824: Calculate the vector inner product of the first and second vectors. If the vector inner product is less than 0, remove the first tunnel point. Calculate the third and fourth vectors. If the vector inner product is less than 0, remove the last tunnel point.

8. A tunnel survey point fitting system, characterized by, For implementing the method of any one of claims 1-7, comprising: a data acquisition module for acquiring the tunnel point coordinate data set and the tunnel contour line data; an overbreak and underbreak value calculation module for calculating the overbreak and underbreak values based on the tunnel point coordinate data set, and determining the average overbreak and underbreak values based on the overbreak and underbreak values; an offset value acquisition module for determining the offset value of the tunnel contour line based on the average overbreak and underbreak values; an abnormal point determination module for determining the abnormal points of all tunnel points based on a preset point optimization threshold, and obtaining the continuous abnormal point set and the discrete abnormal point set; an offset module for offsetting the tunnel contour line as a whole based on the offset value of the tunnel contour line; an abnormal point state determination module for grouping the continuous abnormal point set, fitting and preprocessing based on each group of abnormal points, and obtaining the state of each group of abnormal points until the preset state is met, and then jumping to the bisector point calculation module; a bisector point calculation module for connecting the adjacent normal points before and after each group of abnormal points to form a straight line, and calculating the bisector point on the straight line based on the number of abnormal points in the group; a point fitting module for projecting the bisector point on the straight line corresponding to each group of abnormal points onto the offset tunnel contour line for fitting, calculating the projection point coordinates, modifying the coordinates of each group of abnormal points to the corresponding projection point coordinates, and obtaining the tunnel point coordinate data set after fitting all abnormal points after processing the first and last abnormal points in the tunnel point coordinate data set.

9. An electronic device, comprising: Comprising: one or more processors; a memory for storing one or more programs; when the one or more programs are executed by the one or more processors, the processor implements the method of any one of claims 1-7.

10. A computer-readable storage medium containing a computer program, on which a computer program is stored, characterized in that, The program is executed by one or more processors to implement the method of any one of claims 1-7. The program is executed by one or more processors to implement the method of any one of claims 1-7.

Images