Automatic fitting method for toe board line of concrete face rockfill dam

Abstract

The invention discloses an automatic fitting method for a toe board line of a concrete face rockfill dam, and belongs to the technical field of water conservancy and hydropower engineering. According to the method, on the basis of a small amount of input data and original topographic data, sectioning is carried out on the basis of the topographic data to obtain a space curve, the toe slab line of the concrete faced rockfill dam is obtained through fitting of data points of the space curve, key coordinate points of the toe slab line are obtained, and fitting of the toe slab line of the concrete faced rockfill dam is completed. The time and labor cost is reduced, a three-dimensional concrete faced rockfill dam model can be quickly built based on the fitted concrete faced rockfill dam toe slab line subsequently, and then the concrete faced rockfill dam modeling design efficiency is greatly improved. The problems that in the current construction excavation process, when revealed geological conditions change, the toe board type in the space range needs to be adjusted in time, and a traditional mode is difficult to respond in time are solved.

Description

Technical Field

[0001] The invention belongs to the technical field of water conservancy and hydropower engineering, and in particular relates to an automatic fitting method for a toe plate line of a concrete face rockfill dam. Background Art

[0002] The toe plate is an important part of the anti-seepage system of the face rockfill dam. It is arranged around the face plate and is generally built on the riverbed and the bedrock on both sides. It not only ensures the impermeable connection between the face plate and the dam foundation, but also serves as a cover plate for foundation grouting. It is also the starting working surface of the face plate slipform. The design quality of the toe plate is directly related to the quality, safety and economy of the entire face rockfill dam. In the practice of face rockfill dam engineering in my country, the intersection line of the extended surface of the bottom surface of the face plate and the toe plate foundation is generally used as the toe plate line (i.e., "X line"), as the control line for the toe plate body design.

[0003] Concrete face rockfill dam is a local material dam with the characteristics of relatively simple structure, convenient construction, low cost and strong adaptability. It has been widely used in water conservancy and hydropower projects. The toe plate is one of the most important anti-seepage structures of the face dam and is the main reference structure for the design and construction of the face dam.

[0004] The panel rockfill dam involves the toe plate structure and multiple material partitions of the three-dimensional space alignment. Similarly, the toe plate line is a spatial multi-segment line. Its alignment needs to take into account the principles of foundation conditions, small engineering volume, and convenient construction, and minimize the turning angle as much as possible. Conventional toe plate alignment is carried out using CAD. First, the toe plate plane axis is determined according to the terrain and geological conditions, and then the profile is designed in two dimensions. When the dam design slope, toe plate starting position and any other structural parameters and terrain and geological conditions are adjusted, the designers need to adjust the toe plate line and re-match the horizontal, longitudinal and cross-sectional points. The entire design process is cumbersome, inefficient, and has a high error rate. It requires close cooperation with the geological profession, and the amount of repeated work between the two professions is large. At the same time, manual step-by-step modeling is complicated, the modeling steps are cumbersome, and the learning cost is high. When the geological conditions change during the construction and excavation process, it is necessary to adjust the toe plate type within the space in a timely manner. The traditional method is difficult to respond in time. Summary of the invention

[0005] The purpose of the present invention is to address the problem that when geological conditions change during the current construction and excavation process, it is necessary to timely adjust the toe plate type within the spatial range, and the traditional method is difficult to respond in time. A method for automatically fitting the toe plate line of a concrete panel rockfill dam is proposed. The present invention is based on a small amount of input data and original terrain data, and cuts based on terrain data to obtain a spatial curve. The panel dam toe plate line is obtained by fitting the spatial curve data points, and the key coordinate points of the toe plate line are obtained to complete the fitting of the panel dam toe plate line. This solution can quickly and automatically complete the drawing of the panel dam toe plate line. Subsequently, the panel dam toe plate line based on the fitting can be quickly built into a three-dimensional panel dam model, which greatly improves the efficiency of panel dam modeling and design. At the same time, combined with other automatic modeling methods of panel rockfill dams, the purpose of automatically generating a panel rockfill dam based on the existing terrain and only through a dam axis is finally achieved.

[0006] In order to solve the above technical problems, the present invention provides a method for automatically fitting the toe line of a concrete face rockfill dam, comprising the following steps:

[0007] S1. Obtain the terrain data within the target design range and the spatial curve of the plane where the toe line is located, determine the offset direction according to the proposed panel dam top axis line segment, propose the buried depth distance and extract the three-dimensional spatial coordinate points of the spatial curve, determine the target dam top elevation and filter the data points to obtain valid data points;

[0008] S2, dividing the valid data points into riverbed data points and riverbank data points, and dividing the riverbank data points into left bank data points and right bank data points according to a judgment rule;

[0009] S3, determining the leftmost point and the rightmost point of the riverbed based on the riverbed data points, classifying the riverbed according to the leftmost point and the rightmost point of the riverbed, and performing fitting to obtain a riverbed topographic line;

[0010] S4, determining left bank feature points and right bank feature points through the left bank data points and the right bank data points, and segmenting the left bank data points and the right bank data points according to the left bank feature points and the right bank feature points, thereby obtaining left bank fitting points and right bank fitting points;

[0011] S5, based on the riverbed topography line and the left bank burial depth points and the right bank burial depth points, combined with the proposed burial depth distance, the toe plate line is shifted in the proposed shift direction as a whole to obtain the toe plate line, and then the toe plate line junction point is obtained;

[0012] S6. Based on the toe board line joint point, the left bank buried depth point and the right bank buried depth point are combined to form a complete panel dam toe board line.

[0013] Preferably, the S1 comprises:

[0014] S11, obtaining topographic coordinate point cloud data and topographic line graphics within a certain range around the target dam axis, using the proposed panel dam axis and upstream surface to cut the original topographic line, and obtaining a spatial curve on the plane where the toe board line is located;

[0015] S12, according to the proposed face plate dam top axis line segment, determine the center point of the face plate dam top axis, and offset it to the left bank by a custom distance, recorded as an offset point, propose a buried depth distance, and extract the three-dimensional space coordinate point of the space curve;

[0016] S13, determine the dam crest elevation, and filter out data points that are less than the dam crest elevation through data identification to obtain valid data points.

[0017] Preferably, S2 comprises:

[0018] S21, based on the valid data points, using a scripting language to select the three-dimensional spatial coordinate point with the lowest elevation among the valid data points as the riverbed data point, and the remaining valid data points as the riverbank data points;

[0019] S22, according to the river bank data points, further dividing the river bank data points into left bank data points and right bank data points by a judgment rule;

[0020] S23, selecting the lowest point of the riverbed elevation among the riverbed data points through data identification.

[0021] Preferably, in S22, the judgment rule includes, based on the lowest point of the riverbed elevation and the center point of the dam top axis of the panel dam, connecting the two points to obtain a straight line equation, calculating the Y value, if Y>0, dividing the river bank data points into left bank data points; if Y<0, dividing the river bank data points into right bank data points.

[0022] Preferably, the calculated Y value is calculated by the following formula:

[0023] Y=(C m XB 1 X)*(PY-B 1 Y)-(PX-B 1 X)*(C m YB 1 .Y); where C m is the coordinate of the lowest point of the riverbed elevation, and its coordinate is C m (C m X, C m Y); B 1 is the coordinate of the center point of the top axis of the panel dam, and its coordinate is B 1 (B 1 X, B 1 Y); P is the coordinate of any point in the river bank data point.

[0024] Preferably, S3 includes:

[0025] S31, connecting any two points among the left bank data points to form a first straight line, selecting any point among the riverbed data points to draw a perpendicular line to the first straight line, and obtaining the perpendicular line distance;

[0026] S32, selecting a point among the riverbed data points corresponding to the smallest vertical distance as the leftmost point of the riverbed, and similarly obtaining the rightmost point of the riverbed;

[0027] S33, classifying the leftmost point of the riverbed into the left bank data point, and classifying the rightmost point of the riverbed into the right bank data point;

[0028] S34, connecting the leftmost point of the riverbed and the rightmost point of the riverbed to form a riverbed topographic line.

[0029] Preferably, in S31, the first straight line is specifically expressed as:

[0030] Q: Y = kX + b;

[0031] The vertical distance is specifically expressed as:

[0032] Where k is the slope of the first straight line Q, and b is the intercept; (X 0 , Y 0 ) is any point in the riverbed data.

[0033] Preferably, S4 includes:

[0034] S41, obtaining the highest point of the left bank elevation and the lowest point of the left bank elevation among the left bank data points through data recognition, and taking the highest point of the left bank elevation and the lowest point of the left bank elevation as starting points, connecting the highest point of the left bank elevation and the lowest point of the left bank elevation to obtain a second straight line;

[0035] S42, based on any point in the left bank data points, draw a perpendicular line to the second straight line to obtain a vertical distance, repeat the above process until all points in the left bank data points are traversed, and construct a second straight line vertical distance set;

[0036] S43, based on the second straight line perpendicular distance set, selecting a point corresponding to a perpendicular line distance with a maximum value as a first feature point, and the corresponding perpendicular line distance is a first perpendicular line distance;

[0037] S44, based on the first feature point, the second straight line is divided into two segments, and the highest point of the left bank elevation and the first feature point of the left bank are respectively connected to form two straight lines, so as to obtain a third straight line and a fourth straight line; all points of the left bank data points are respectively perpendicular to the third straight line and the fourth straight line to obtain a set of vertical distances of the third straight line;

[0038] S45, selecting the point with the largest distance value in the third straight line vertical distance set as the second feature point, and the corresponding vertical line distance is the second vertical line distance; and so on until the third feature point and its corresponding third vertical line distance, the fourth feature point and its corresponding fourth vertical line distance are obtained;

[0039] S46, based on the first vertical distance, the second vertical distance, the third vertical distance and the fourth vertical distance, an average value is calculated, and the first vertical distance, the second vertical distance, the third vertical distance and the fourth vertical distance are respectively compared with the average value, and characteristic points corresponding to the distances greater than the average value are selected as left bank characteristic points, and the left bank data points are segmented by combining the lowest point of the left bank elevation and the highest point of the left bank elevation to obtain left bank data point segments;

[0040] S47, obtaining the right bank characteristic points according to the derivation method of S41-S46, combining the right bank lowest elevation point and the right bank lowest elevation point, segmenting the right bank data points to obtain right bank data point segments;

[0041] S48, based on the left bank data point segments, fitting each data point segment to obtain a corresponding straight line;

[0042] S49, bringing the horizontal coordinates of the highest point of the left bank elevation and the left bank feature point into the corresponding straight line fitting, respectively, to obtain the left bank fitting point, and fitting between two adjacent points to obtain the left bank terrain fitting line;

[0043] S410, obtaining right bank feature points according to the derivation method of S48-S49, and then obtaining right bank fitting points and right bank terrain fitting lines.

[0044] Preferably, the S5 comprises:

[0045] S51, based on the riverbed topographic line, offset the proposed burial depth distance along the perpendicular direction of the riverbed topographic line to obtain the rightmost burial depth point of the riverbed and the leftmost burial depth point of the riverbed;

[0046] S52, based on the left bank fitting point, the proposed burial depth distance is offset along the dam crest axis direction (offset direction) of the panel dam to obtain left bank burial depth points, namely, a first left bank burial depth point, a second left bank burial depth point, a third left bank burial depth point and a fourth left bank burial depth point;

[0047] S53, similarly, the right bank burial depth points can be obtained, which are the first right bank burial depth point, the second right bank burial depth point, the third right bank burial depth point and the fourth right bank burial depth point;

[0048] S54, connect the first left bank buried depth point and the second left bank buried depth point to obtain the left bank toe line, connect the rightmost buried depth point of the riverbed and the leftmost buried depth point of the riverbed to obtain the riverbed toe line, extend the left bank toe line and the riverbed toe line to intersect them to obtain the junction point of the riverbed and the left bank toe line, and similarly obtain the junction point of the riverbed and the right bank toe line, both of which are the toe line junction points.

[0049] Preferably, in S51, the rightmost buried depth point of the riverbed is:

[0050]

[0051] The leftmost buried depth of the riverbed is:

[0052]

[0053] In the formula, C L X, C L Y is C L Coordinate, C R X, C R Y is C R Coordinate. k C 、b C For straight line C L C R coefficient.

[0054] Preferably, S6 includes connecting the first left bank buried depth point, the second left bank buried depth point, the third left bank buried depth point and the fourth left bank buried depth point, the junction point of the riverbed and the left bank toe board line, the junction point of the riverbed and the right bank toe board line, the first right bank buried depth point, the second right bank buried depth point, the third right bank buried depth point and the fourth right bank buried depth point in S5 to form a complete panel dam toe board line.

[0055] Compared with the prior art, the present invention has the following beneficial effects:

[0056] 1. This scheme uses a small amount of input data and original terrain data, and obtains a spatial curve by sectioning based on the terrain data. The panel dam toe line is obtained by fitting the spatial curve data points, and the key coordinate points of the toe line are obtained to complete the fitting of the panel dam toe line. This scheme can quickly and automatically complete the drawing of the panel dam toe line. Subsequently, the panel dam toe line based on the fitting can be quickly built into a three-dimensional panel dam model, which greatly improves the efficiency of panel dam modeling and design. Finally, a panel rockfill dam can be automatically generated based on the existing terrain through only one dam axis, which is very efficient and convenient.

[0057] 2. This solution uses scripting language to identify the existing spatial curve coordinate points to select valid data points, find the terrain line feature points among the valid data points, perform segmented fitting on the feature points, and automatically extract the panel dam toe line. Using the panel dam toe line as the basis for 3D modeling greatly improves the efficiency of modeling.

[0058] 3. This scheme automatically identifies the terrain features around the panel dam building based on the spatial curve and the attached terrain feature information of the spatial curve, constructs the toe line, assists manual design, and greatly improves the efficiency of panel dam design. BRIEF DESCRIPTION OF THE DRAWINGS

[0059] Other features, objects and advantages of the present invention will become more apparent by reading the detailed description of non-limiting embodiments made with reference to the following drawings. The drawings are only for the purpose of illustrating preferred embodiments and are not to be considered as limiting the present invention. Also, the same reference symbols are used throughout the drawings to represent the same parts.

[0060] Figure 1 It is a schematic diagram of the system flow of the present invention;

[0061] Figure 2 It is a schematic diagram of the compilation process of the present invention;

[0062] Figure 3 It is a schematic diagram of a straight line formed by the highest point of the left bank elevation and the lowest point of the left bank elevation in the present invention;

[0063] Figure 4 It is a segmented schematic diagram of characteristic points of the present invention;

[0064] Figure 5 It is a further segmented schematic diagram of the characteristic points of the present invention. DETAILED DESCRIPTION

[0065] In order to make the objectives, technical solutions and advantages of the present invention more clearly understood, the present invention is further described in detail below in conjunction with the accompanying drawings and embodiments. It should be understood that the specific implementation method described herein is only an optimal embodiment of the present invention, which is only used to explain the present invention and does not limit the scope of protection of the present invention. All other embodiments obtained by ordinary technicians in this field without making creative work are within the scope of protection of the present invention.

[0066] Example 1: Figure 1 - Figure 5 As shown, a method for automatically fitting the toe plate line of a concrete face rockfill dam comprises:

[0067] Step 1: Prepare the topographic coordinate point cloud data (i.e. topographic data points) within a certain range around the dam axis of this project. Based on the topographic line graphics within a certain range around the dam axis of this project, use the proposed panel dam axis and upstream surface to cut the original topographic line, and obtain a spatial curve LL on the plane where the toe line is located. According to the proposed panel dam top axis line segment, determine the midpoint B of the line segment 1 , and shift it to the left bank by 10m, denoted as B 2 Point. Connect B 1 B 2 As the proposed offset direction and the proposed buried depth distance. Extract the series of point coordinates (X, Y, Z) of the LL space curve to construct a data point set, which is a three-dimensional space coordinate point. After determining the dam crest elevation, the data points below the dam crest elevation are filtered out through script language data recognition, and the remaining data points are used as valid data points.

[0068] Step 2: Based on the valid data points, filter out the riverbed data points according to the specified rules, and use the remaining valid data points as riverbank data points, and further divide the riverbank data points into left bank data points and right bank data points.

[0069] Further, step 2 includes the following steps:

[0070] Step 2.1 Based on the valid data points {M} in step 1, the script language is used to automatically select the series of coordinate points with the lowest elevation among the valid data points {M} as the riverbed data points {HC}. The valid data points {M}, except the riverbed data points {HC}, are the riverbank data points {HA}. The mathematical expression is: {HA}U{HC}={M}; where, valid data points {M}; riverbed data points {HC}; riverbank data points {HA}.

[0071] In step 2.2, the riverbank data point set {HA} is further split into left bank data points {HL} and right bank data points {HR} through judgment rules.

[0072] Step 2.3. Further, among the riverbed data points {HC}, the lowest riverbed elevation point C is identified by the data m .

[0073] Specifically, in S22, the judgment rules include: based on C m The center point B of the top axis of the panel dam entered in step 1 1 , connect B 1 C mThe straight line equation is used to divide the river bank data point {HA} into the left bank data point {HL} and the right bank data point {HR}, and calculate the Y value. If Y>0, the river bank data point is divided into the left bank data point; if Y<0, the river bank data point is divided into the right bank data point. The Y value is calculated by the following formula:

[0074] Y=(C m XB 1 X)*(PY-B 1 Y)-(PX-B 1 X)*(C m YB 1 .Y); where C m is the coordinate of the lowest point of the riverbed elevation, and its coordinate is C m (C m X, C m Y); B 1 is the coordinate of the center point of the top axis of the panel dam, and its coordinate is B 1 (B 1 X, B 1 Y); P is the coordinate of any point in the river bank data point.

[0075] Step 3: Based on the riverbed data points, find its key points, that is, the leftmost point and the rightmost point of the riverbed. The key points divide the riverbed (riverbed data points) into segments, and fit the terrain lines to the terrain data points in the segments.

[0076] Step 3.1 Connect any two points in the left bank data points {HL} to form a straight line Q (first straight line). Select any point in the riverbed data points {M} and draw a perpendicular line to the straight line. The perpendicular distance is D. Q ;

[0077] Q: Y = kX + b;

[0078] In the formula, k is the slope of the straight line Q, and b is the intercept; (X 0 , Y 0 ) is any point in the riverbed data points {M}. Select the point in the riverbed data points corresponding to the smallest vertical distance as the leftmost point C of the riverbed L Similarly, we can get the rightmost point C of the riverbed R .

[0079] Step 3.2 Leftmost point C of the riverbed L Add the left bank data point {HL} and the rightmost point C of the riverbed R Add the right bank data point {HR}.

[0080] Step 3.3 Connect C L C R , since most of the riverbed data points are along CL C R Evenly distributed, so directly connect C L C R Form a riverbed fitting terrain line.

[0081] Step 4: Determine the left bank feature points and the right bank feature points through the left bank data points and the right bank data points, and segment the left bank data points and the right bank data points according to the left bank feature points and the right bank feature points, so as to obtain the left bank fitting points and the right bank fitting points.

[0082] Step 4 is further based on the riverbank data point {HA} determined in step 3. Use the algorithm to find the feature point in the left bank data point {HL}. The feature point is the point farthest from the distribution in the set. The specific steps are as follows:

[0083] Step 4.1, the data identifies the left bank data points {HL}, the highest point L on the left bank a The lowest point on the left bank is L i As the starting point of the algorithm. Connect L a L i , L a L i The equation of the straight line is:

[0084] L a L i :Ax+by+C=0 Where A, b, C are straight line L a L i The coefficient of .

[0085] Step 4.2, use any point (X 0 , Y 0 ) a L i Draw a vertical line with a vertical distance of H. Repeat the above process until all elements in the left bank data point {HL} are traversed. Where H:

[0086] Where A, B, C are straight line L a L i The coefficient of (X 0 , Y 0 ) is a point among the left bank data points {HL}.

[0087] Determine the point corresponding to the maximum vertical distance as the first feature point L on the left bank 1 , the corresponding vertical distance is H 1 .

[0088] Step 4.3, based on the feature point L found in step 4.2 1 L a L iDivided into two sections. Connect L a L 1 and L 1 L i Then draw perpendicular lines from all the left bank data points {HL} to the two straight lines, and take the point with the largest perpendicular distance as the second feature point L 2 , the corresponding vertical distance is H 2 And so on until the third feature point L is obtained. 3 , the corresponding vertical distance is H 3 , the fourth feature point L 4 , the corresponding vertical distance is H 4 .

[0089] Step 4.4, based on H 1 , H 2 , H 3 , H 4 Calculate the average value H v , determine the final feature point, H v Average distance:

[0090] Where L 1 , L 2 , L 3 , L 4 The corresponding distance H 1 , H 2 , H 3 , H 4 .

[0091] Step 4.5, select H 1 , H 2 , H 3 , H 4 The feature points corresponding to values ​​greater than the average value are taken as the left bank feature points. The left bank data points are segmented by combining the lowest point in elevation of the left bank and the highest point in elevation of the left bank to obtain the left bank data point segments.

[0092] It should be noted that “Select H 1 , H 2 , H 3 , H 4 There are many possible situations in which the feature point corresponding to the value greater than the average value is regarded as the left bank feature point. There may be only one feature point, two feature points, or three feature points.

[0093] For example, if H1>Hv (only the average value of the first vertical line distance), the corresponding L1 (first feature point) is the left bank feature point. Combined with the lowest elevation point and the highest elevation point of the left bank, the left bank data points are divided into two sections, namely, from the lowest elevation point of the left bank to the first feature point, and from the first feature point to the highest elevation point of the left bank.

[0094] For example, if it is determined that H 1> H v , H 2> H v , then the corresponding L 1 , L 2 is the characteristic point of the left bank, combined with the lowest point L on the left bank i 、The highest point on the left bank is L a .L i , L 1 , L 2 , L a The left bank data points {HL} are divided into three segments: L i L 1 , L 1 L 2 , L 2 L a .

[0095] For example, if H1>Hv, H2>Hv, H3>Hv (the first vertical line distance, the second vertical line distance, and the third vertical line distance are all greater than the average value), then the corresponding L1, L2, L3 (the first feature point, the second feature point, and the third feature point) are the left bank feature points. Combined with the lowest elevation point and the highest elevation point of the left bank, the left bank data points are divided into four sections, namely, from the lowest elevation point of the left bank to the first feature point, the first feature point to the second feature point, the second feature point to the third feature point, and the third feature point to the highest elevation point of the left bank.

[0096] Among them, the case with only one characteristic point is to choose one from four; the case with two characteristic points is to choose two from four, and the case with three characteristic points is to choose three from four. It is impossible for all four results to be greater than the average.

[0097] Similarly, the right bank feature point R can be obtained i , R 1 , R 2 , R a , combined with the lowest point of the right bank elevation and the lowest point of the right bank elevation, the right bank data points are segmented to obtain the right bank data point segments.

[0098] Step 4.6, based on the left bank data point segment, each data point segment is fitted to obtain a corresponding straight line. The horizontal coordinates of the highest point of the left bank elevation and the left bank feature point are respectively substituted into the corresponding straight line fitting to obtain the left bank fitting point, and the left bank terrain fitting line is obtained by fitting between two adjacent points.

[0099] Similarly, the right bank feature points can be obtained, and then the right bank fitting points and the right bank terrain fitting line can be obtained.

[0100] Further, in step 4.5, based on the three segments L divided into the left bank in step 4.3i L 1 , L 1 L 2 , L 2 L a The terrain data points are fitted. The three straight lines are:

[0101] L i L 1 :Y 1 =k 1 x+b 1

[0102] L 1 L 2 :Y 2 =k 2 x+b 2

[0103] L 2 L a :Y 2 =k 3 x+b 3 ;

[0104] Where k and b are the coefficients in the straight line equation.

[0105] It is understandable that step 4.6 can be specifically as follows: 1 Horizontal axis L 1 X into Y 1 =k 1 x+b 1 Fitting point N 1 (L 1 X,k 1 L 1 X+b 1 ). Similarly, we can get the fitting point N i、 N 2 (L 2 X,k 2 L 2 X+b 2 ), N a (L a X,k 3 L a X+b 3 ). Left bank terrain fitting line N i N 1 、N 1 N 2 、N 2 Na. Among them, N i The lowest point on the left bank is L i, because a continuous straight line is formed as the left bank terrain fitting line, the lowest point of the left bank elevation is fixed as the riverbed coordinate point, so no fitting is required. At the same time, because the connection between the riverbed and the riverbank is the same point, it is used as the base point of the terrain fitting line. Therefore, C l The same point as Ni and Li. Figure 2 As shown, Cr, Ri, and Ei in the figure are the same point and the principle is similar.

[0106] Similarly, the right bank fitting point E can be obtained i 、E 1 、E 2 、E a and the right bank topographic fitting line E i E 1 、E 1 E 2 、E 2 E a .

[0107] Step 5: Based on the riverbed topography line and the left bank burial depth points and the right bank burial depth points, combined with the proposed burial depth distance, the toe line is shifted as a whole in the proposed offset direction to obtain the toe line, and then the toe line junction point is obtained.

[0108] In step 5.1, the riverbed topography line C fitted in step 3 is further L C R , along C L C R The vertical line direction is offset by the burial depth d proposed in step 1, and S is obtained. 1 The deepest point on the right side of the riverbed, S 2 The deepest point on the left side of the riverbed.

[0109] C L C R :Y C =k C x+b C ;

[0110]

[0111] In the formula, C L X, C L Y is C L Coordinate, C R X, C R Y is C R Coordinate. k C 、b C For straight line C L C R coefficient.

[0112] Step 5.2 For the left bank terrain fitting line N in step 4 i N 1、N 1 N 2 、N 2 N a Middle fitting point N i 、N 1 、N 2 、N a Along the top axis B of the concrete face dam 1 B 2 Direction deviation proposed burial depth distance d, get the left bank burial depth point P i , P 1 , P 2 , P a .

[0113] B 1 B 2 : Y b =k b x+b b ;

[0114] In the formula, B 1 B 2 is the top axis equation of the panel dam, k b 、b b is the coefficient, and d is the burial depth.

[0115]

[0116] Where N i X,N i Y is N i Coordinate. k b 、b b For B 1 B 2 is the coefficient of the top axis equation of the panel dam, and d is the buried depth.

[0117]

[0118] Where N 1 X,N 1 Y is N 1 Coordinate. k b 、b b For B 1 B 2 is the coefficient of the top axis equation of the panel dam, and d is the buried depth.

[0119]

[0120] Where N 2 X,N 2 Y is N 2 Coordinate. k b 、b b For B 1 B2 is the coefficient of the top axis equation of the panel dam, and d is the buried depth.

[0121]

[0122] Where N a X,N a Y is N a Coordinate. k b 、b b For B 1 B 2 is the coefficient of the top axis equation of the panel dam, and d is the buried depth.

[0123] Similarly, we can get the right bank burial depth W i , W 1 , W 2 , W a .

[0124] Step 5.3: Further, the left bank toe line P i P 1 With the riverbed toeboard line S 1 S 2 Extend the intersection to get the junction point S between the riverbed and the left bank toe line L It can be expressed as:

[0125] Y=k P x+bP;

[0126] P i P 1 :

[0127] Y=k s x+bs;

[0128] S 1 S 2 :

[0129] S L :

[0130]

[0131] where k p 、b p The straight line P i P 1 Coefficient, k s 、b s The straight line S 1 S 2 Coefficient, S L X, S L Y is S L coordinate.

[0132] Similarly, the junction point S between the riverbed and the right bank toe line can be obtainedR .

[0133] Step 6: Based on the toe board line joint points, combined with the left bank burial depth points and the right bank burial depth points, a complete panel dam toe board line is formed, the toe board line terrain data coordinate points are stored, and the toe board line coordinate points are used for 3D modeling.

[0134] Step 6: Further, the left bank burial depth point P in step 5 is a , P 2 , P 1 , P i , the junction point S between the riverbed and the river bank toe line L , S R , right bank burial depth point W i , W 1 , W 2 , W a The toe board line of the panel dam is connected to form a complete one. The toe board line coordinate points can be used to set out according to the slope ratio to model in the 3D model.

[0135] This solution uses a small amount of input data and original terrain data, and obtains a spatial curve based on the terrain data. The panel dam toe line is obtained by fitting the spatial curve data points, and the key coordinate points of the toe line are obtained to complete the fitting of the panel dam toe line. This solution can quickly and automatically complete the drawing of the panel dam toe line. Subsequently, the panel dam toe line based on the fitting can be quickly built into a three-dimensional panel dam model, greatly improving the efficiency of panel dam modeling and design.

[0136] Example 2: Taking a concrete face rockfill dam project as an example, a method for automatically fitting the toe line of a concrete face rockfill dam is combined with a schematic diagram of the drawing process, such as Figure 2 As shown, follow the steps below:

[0137] Step 1 Extract the coordinates (X, Y, Z) of a series of points of the LL space curve. The data point set is a three-dimensional space coordinate point. According to the proposed face plate dam top axis line segment, determine the center point B of the face plate dam top axis. 1 (1033.567, 209.212, 1085), and shift it 10m to the left bank to point B on the left side of the center point of the top axis of the panel dam 2 (1043.856, 219.513, 1085). The data points below the dam crest elevation of 1085 m are filtered out through script language data recognition, and the remaining data points are taken as valid data points {M}.

[0138] Step 2: Based on the valid data points, filter out riverbed data points according to the specified rules from the valid data points, use the remaining valid data points as riverbank data points, and further divide the riverbank data points into left bank data points and right bank data points.

[0139] Step 2 includes the following steps:

[0140] Step 2.1, based on the valid data points {M} in step 1, the script language is used to automatically select the series of coordinate points with the lowest elevation among the valid data points {M} as the riverbed data points {HC}. The valid data points {M}, except the riverbed data points {HC}, are the riverbank data points {HA}.

[0141] {HA}U{HC}={M}; where, valid data point {M}; riverbed data point {HC}; riverbank data point {HA}.

[0142] Step 2.2, the middle river bank data point set {HA} is further split into left bank data points {HL} and right bank data points {HR} by judgment.

[0143] Step 2.3, among the riverbed data points {HC}, identify the lowest point C of the riverbed elevation through the data m (906.047, 296.49)

[0144] Step 2.4, based on the lowest point C of the riverbed elevation in step 2.3 m The center point B of the top axis of the panel dam entered in step 1 1 , connect B 1 C m , calculate the Y value. If Y>0, the river bank data point is divided into the left bank data point; if Y<0, the river bank data point is divided into the right bank data point. The Y value is calculated by the following formula:

[0145] Y=(C m XB 1 X)*(PY-B 1 Y)-(PX-B 1 X)*(C m YB 1 .Y); where C m is the coordinate of the lowest point of the riverbed elevation, and its coordinate is C m (C m X, C m Y); B 1 is the coordinate of the center point of the top axis of the panel dam, and its coordinate is B 1 (B 1 X, B 1 Y); P is the coordinate of any point in the river bank data point.

[0146] Step 3: Find the rightmost point C of the riverbed based on the riverbed data point {M} in step 2 R , the leftmost point of the riverbed C L .

[0147] The steps are as follows:

[0148] Step 3.1, connect any two points of the left bank data points {HL} to form a straight line Q. Draw a perpendicular line from any point of the riverbed data points {M} to the straight line Q, and the perpendicular distance is

[0149]

[0150] Q: Y = kX + b;

[0151] In the formula, k is the slope of the straight line Q, and b is the intercept; (X 0 , Y 0 ) is any point in the riverbed data points {M}. Determine the minimum D Q The corresponding point in {M} is the leftmost point C of the riverbed L Similarly, we can get the rightmost point C of the riverbed R .

[0152] Step 3.2, leftmost point C of the riverbed L Add the left bank data point {HL} and the rightmost point C of the riverbed R Add the right bank data point {HR}.

[0153] Step 3.3, Connect C L C R , since the riverbed assembly points are mostly along C L C R Evenly distributed, so directly connect C L C R Form a riverbed fitting terrain line.

[0154] Specifically, the key data points of the riverbed are shown in Table 1.

[0155] Table 1 Key data points of riverbed

[0156]

[0157] Step 4: Based on the riverbank data point {HA} determined in step 3, use the algorithm to find the feature point in the left bank data point {HL}. The feature point is the point farthest from the distribution in the set. The specific algorithm steps are as follows:

[0158] Step 4.1, the data identifies the left bank data points {HL}, the highest point L on the left bank a The lowest point on the left bank is L i As the starting point of the algorithm. Connect L a L i , L a L i The equation of the straight line is:

[0159] L a L i:Ax+by+C=0;where A, B, C are straight line L a L i The coefficient of .

[0160] Use any point on the left bank data point {HL} (X 1 , Y 1 ) a L i Draw a vertical line with a vertical distance of H until all elements in the left bank data point {HL} are traversed. Where H:

[0161] Where A, B, C are the coefficients of the straight line LaLi; (X 0 , Y 0 ) is a point in the left bank data point {HL}. The point corresponding to the maximum vertical distance is determined as the first feature point L on the left bank. 1 , the corresponding vertical distance is H 1 .

[0162] Step 4.2, based on the feature point L found in step 4.1 1 L a Li is divided into two sections. Connect L a L 1 , L 1 L i Two straight lines are formed. Then draw perpendicular lines from all the left bank data points {HL} to this straight line, and take the point with the largest distance L 2 , the corresponding vertical distance is H 2 And so on until the third feature point L is extracted 3 , the corresponding vertical distance is H 3 The fourth feature point L 4 , the corresponding vertical distance is H 4 .

[0163] Step 4.3, find the feature point L 1 , L 2 , L 3 , L 4 The corresponding distance H 1 , H 2 , H 3 , H 4 . Determine the final feature point, H v Average distance.

[0164] Where L 1 , L 2 , L 3 , L 4 The corresponding distance H 1 , H 2 , H3 , H 4 .

[0165] Select H 1 , H 2 , H 3 , H 4 The feature points corresponding to the values ​​greater than the average value are taken as the feature points of the left bank. The left bank data points are segmented by combining the lowest point of the left bank elevation and the highest point of the left bank elevation to obtain the left bank data point segments. Similarly, the feature points R of the right bank can be obtained. i , R 1 , R 2 , R a , combined with the lowest point of the right bank elevation and the lowest point of the right bank elevation, the right bank data points are segmented to obtain the right bank data point segments.

[0166] Specifically, the key data points of the river bank are shown in Table 2.

[0167] Table 2 Key data points of riverbank

[0168]

[0169] Step 4.4 is based on the three segments L divided into the left bank in step 4.3 i L 1 , L 1 L 2 , L 2 L a The terrain data points are fitted. The three straight lines are:

[0170] L i L 1 :Y 1 =k 1 x+b 1

[0171] L 1 L 2 :Y 2 =k 2 x+b 2

[0172] L 2 L a :Y 2 =k 3 x+b 3 ;

[0173] Where k and b are the coefficients in the straight line equation.

[0174] Step 4.5, L 1 Horizontal axis L 1 X into Y 1 =k 1 x+b 1Fitting point N 1 (L 1 X,k 1 L 1 X+b 1 ). Similarly, we can get the fitting point N i、 N 2 (L 2 X,k 2 L 2 X+b 2 ), N a (L a X,k 3 L a X+b 3 ). Left bank terrain fitting line N i N 1 、N 1 N 2 、N 2 Na. Similarly, the right bank fitting point E can be obtained i 、E 1 、E 2 、E a and the right bank topographic fitting line E i E 1 、E 1 E 2 、E 2 E a .

[0175] Specifically, the coefficients of the left bank straight line equation are shown in Table 3.

[0176] Table 3 Coefficients of the left bank straight line equation

[0177] k1 -0.2819016 b1 612.16219 k2 0.5533752 b2 -272.5850904 k3 -0.8352826 b3 1257.874212

[0178] Specifically, the left bank fitting points are shown in Table 4.

[0179] Table 4 Left bank fitting points

[0180]

[0181] Specifically, the right bank fitting points are shown in Table 5.

[0182] Table 5 Right bank fitting points

[0183]

[0184] Step 5: Based on the riverbed topography line and the left bank burial depth points and the right bank burial depth points, combined with the proposed burial depth distance, the toe line is shifted as a whole in the proposed offset direction to obtain the toe line, and then the toe line junction point is obtained.

[0185] In step 5.1, the riverbed topography line C fitted in step 3 is further LC R , along C L C R The vertical line direction is offset by the burial depth d proposed in step 1, and S is obtained. 1 The deepest point on the right side of the riverbed, S 2 The deepest point on the left side of the riverbed.

[0186] C L C R :Y C =k C x+b C ;

[0187]

[0188]

[0189] In the formula, C L X, C L Y is C L Coordinate, C R X, C R Y is C R Coordinate. k C 、b C For straight line C L C R coefficient.

[0190] Step 5.2 For the left bank terrain fitting line N in step 4 i N 1 、N 1 N 2 、N 2 N a Middle fitting point N i 、N 1 、N 2 、N a Along the top axis B of the concrete face dam 1 B 2 Direction deviation proposed burial depth distance d, get the left bank burial depth point P i , P 1 , P 2 , P a .

[0191] B 1 B 2 : Y b =k b x+b b ;

[0192] In the formula, B 1 B 2 is the top axis equation of the panel dam, k b 、b b is the coefficient, and d is the burial depth.

[0193]

[0194] Where N i X,N i Y is N i Coordinate. k b 、b b For B 1 B 2 is the coefficient of the top axis equation of the panel dam, and d is the buried depth.

[0195]

[0196] Where N 1 X,N 1 Y is N 1 Coordinate. k b 、b b For B 1 B 2 is the coefficient of the top axis equation of the panel dam, and d is the buried depth.

[0197]

[0198] Where N 2 X,N 2 Y is N 2 Coordinate. k b 、b b For B 1 B 2 is the coefficient of the top axis equation of the panel dam, and d is the buried depth.

[0199]

[0200] Where N a X,N a Y is N a Coordinate. k b 、b b For B 1 B 2 is the coefficient of the top axis equation of the panel dam, and d is the buried depth.

[0201] Similarly, we can get the right bank burial depth W i , W 1 , W 2 , W a .

[0202] Specifically, the left bank burial depth points are shown in Table 6. The right bank burial depth points are shown in Table 7.

[0203] Table 6 Left bank burial depth points

[0204]

[0205] Table 7 Right bank burial depth points

[0206]

[0207] Step 5.3: Further, the left bank toe line P i P 1 With the riverbed toeboard line S 1 S 2 Extend the intersection to get the junction point S between the riverbed and the left bank toe line L It can be expressed as:

[0208] Y=k P x+b P ;

[0209] P i P 1 :

[0210] Y=k s x+b s ;

[0211] S 1 S 2 :

[0212] S L :

[0213]

[0214] where k p 、b p The straight line P i P 1 Coefficient, k s 、b s The straight line S 1 S 2 Coefficient, S L X, S L Y is S L coordinate.

[0215] Similarly, the junction point S between the riverbed and the right bank toe line can be obtained R .

[0216] Specifically, the junction point between the riverbed and the riverbank toe line is shown in Table 8.

[0217] Table 8 The junction of riverbed and riverbank toe line

[0218]

[0219] Step 6: Further, the left bank burial depth point P in step 5 is a , P 2 , P 1 , Pi , the junction point S between the riverbed and the river bank toe line L , S R , right bank burial depth point W i , W 1 , W 2 , W a The toe board line of the panel dam is connected to form a complete one. The toe board line coordinate points can be used to set out according to the slope ratio to model in the 3D model.

[0220] The above specific embodiments are preferred embodiments of the present invention, and are not intended to limit the specific implementation scope of the present invention. The scope of the present invention includes but is not limited to the present specific embodiments. All equivalent changes made in accordance with the shape, structure, and method of the present invention are within the protection scope of the present invention.

Claims

1. A method for automatically fitting the toe line of a concrete face rockfill dam, characterized in that: The steps include: S1. Obtain the terrain data within the target design range and the spatial curve of the plane where the toe line is located, determine the offset direction according to the proposed panel dam top axis line segment, propose the buried depth distance and extract the three-dimensional spatial coordinate points of the spatial curve, determine the target dam top elevation and filter the data points to obtain valid data points; S2, dividing the valid data points into riverbed data points and riverbank data points, and dividing the riverbank data points into left bank data points and right bank data points according to a judgment rule; S3, determining the leftmost point and the rightmost point of the riverbed based on the riverbed data points, classifying the riverbed according to the leftmost point and the rightmost point of the riverbed, and performing fitting to obtain a riverbed topographic line; S4, determining left bank feature points and right bank feature points through the left bank data points and the right bank data points, and segmenting the left bank data points and the right bank data points according to the left bank feature points and the right bank feature points, thereby obtaining left bank fitting points and right bank fitting points; S5, based on the riverbed topography line and the left bank burial depth points and the right bank burial depth points, combined with the proposed burial depth distance, the toe plate line is shifted in the proposed shift direction as a whole to obtain the toe plate line, and then the toe plate line junction point is obtained; S6. Based on the toe board line joint point, the left bank buried depth point and the right bank buried depth point are combined to form a complete panel dam toe board line.

2. The method for automatically fitting the toe line of a concrete face rockfill dam according to claim 1, characterized in that: The S1 includes: S11, obtaining topographic coordinate point cloud data and topographic line graphics within a certain range around the target dam axis, using the proposed panel dam axis and upstream surface to cut the original topographic line, and obtaining a spatial curve on the plane where the toe board line is located; S12, according to the proposed face plate dam top axis line segment, determine the center point of the face plate dam top axis, and offset it to the left bank by a custom distance, record it as an offset point, connect the center point and the offset point as the offset direction, propose the buried depth distance, and extract the three-dimensional space coordinate point of the space curve; S13, determine the dam crest elevation, and filter out data points that are less than the dam crest elevation through data recognition to obtain valid data points.

3. The automatic fitting method for the toe line of a concrete face rockfill dam according to claim 1 is characterized in that: The S2 includes: S21, based on the valid data points, using a scripting language to select the three-dimensional spatial coordinate point with the lowest elevation among the valid data points as the riverbed data point, and the remaining valid data points as the riverbank data points; S22, further dividing the river bank data points into left bank data points and right bank data points according to a judgment rule; S23, selecting the lowest point of the riverbed elevation among the riverbed data points through data identification.

4. The method for automatically fitting the toe line of a concrete face rockfill dam according to claim 3, characterized in that: In S22, the judgment rules include, based on the lowest point of the riverbed elevation and the center point of the axis of the panel dam top, connecting the two points to obtain a straight line equation, calculating the Y value, if Y>0, dividing the river bank data points into left bank data points; if Y<0, dividing the river bank data points into right bank data points.

5. The method for automatically fitting the toe line of a concrete face rockfill dam according to claim 4, characterized in that: The calculated Y value is calculated by the following formula: Y=(C m X-B1X)*(PY-B1Y)-(PX-B1X)*(C m Y-B1.Y); In the formula, C m is the coordinate of the lowest point of the riverbed elevation, and its coordinate is C m (C m X, C m Y); B1 is the coordinate of the center point of the dam crest axis, and its coordinate is B1(B1X, B1Y); P is the coordinate of any point in the river bank data point.

6. The method for automatically fitting the toe line of a concrete face rockfill dam according to claim 1, characterized in that: The S3 includes: S31, connecting any two points among the left bank data points to form a first straight line, selecting any point among the riverbed data points to draw a perpendicular line to the first straight line, and obtaining the perpendicular line distance; S32, selecting a point among the riverbed data points corresponding to the smallest vertical distance as the leftmost point of the riverbed, and similarly obtaining the rightmost point of the riverbed; S33, classifying the leftmost point of the riverbed into the left bank data point, and classifying the rightmost point of the riverbed into the right bank data point; S34, connecting the leftmost point of the riverbed and the rightmost point of the riverbed to form a riverbed topographic line.

7. The method for automatically fitting the toe line of a concrete face rockfill dam according to claim 1, characterized in that: The S4 includes: S41, obtaining the highest point of the left bank elevation and the lowest point of the left bank elevation among the left bank data points through data recognition, and taking the highest point of the left bank elevation and the lowest point of the left bank elevation as starting points, connecting the highest point of the left bank elevation and the lowest point of the left bank elevation to obtain a second straight line; S42, based on any point in the left bank data points, draw a perpendicular line to the second straight line to obtain a vertical distance, repeat the above process until all points in the left bank data points are traversed, and construct a second straight line vertical distance set; S43, based on the second straight line perpendicular distance set, selecting a point corresponding to a perpendicular line distance with a maximum value as a first feature point, and the corresponding perpendicular line distance is a first perpendicular line distance; S44, based on the first feature point, the second straight line is divided into two segments, and the highest point of the left bank elevation and the first feature point of the left bank are respectively connected to form two straight lines, so as to obtain a third straight line and a fourth straight line; all points of the left bank data points are respectively perpendicular to the third straight line and the fourth straight line to obtain a set of vertical distances of the third straight line; S45, selecting the point with the largest distance value in the third straight line vertical distance set as the second feature point, and the corresponding vertical line distance is the second vertical line distance; and so on until the third feature point and its corresponding third vertical line distance, the fourth feature point and its corresponding fourth vertical line distance are obtained; S46, based on the first vertical distance, the second vertical distance, the third vertical distance and the fourth vertical distance, an average value is calculated, and the first vertical distance, the second vertical distance, the third vertical distance and the fourth vertical distance are respectively compared with the average value, and characteristic points corresponding to the distances greater than the average value are selected as left bank characteristic points, and the left bank data points are segmented by combining the lowest point of the left bank elevation and the highest point of the left bank elevation to obtain left bank data point segments; S47, obtaining the right bank characteristic points according to the derivation method of S41-S46, combining the right bank lowest elevation point and the right bank lowest elevation point, segmenting the right bank data points to obtain right bank data point segments; S48, based on the left bank data point segments, fitting each data point segment to obtain a corresponding straight line; S49, bringing the horizontal coordinates of the highest point of the left bank elevation and the left bank feature point into the corresponding straight line fitting, respectively, to obtain the left bank fitting point, and fitting between two adjacent points to obtain the left bank terrain fitting line; S410, obtaining right bank feature points according to the derivation method of S48-S49, and then obtaining right bank fitting points and right bank terrain fitting lines.

8. The method for automatically fitting the toe line of a concrete face rockfill dam according to claim 1, characterized in that: The S5 includes: S51, based on the riverbed topographic line, offset the proposed burial depth distance along the perpendicular direction of the riverbed topographic line to obtain the rightmost burial depth point of the riverbed and the leftmost burial depth point of the riverbed; S52, based on the left bank fitting point, the proposed burial depth distance is offset along the axial direction of the panel dam crest to obtain left bank burial depth points, namely, a first left bank burial depth point, a second left bank burial depth point, a third left bank burial depth point and a fourth left bank burial depth point; S53, similarly, the right bank burial depth points can be obtained, which are the first right bank burial depth point, the second right bank burial depth point, the third right bank burial depth point and the fourth right bank burial depth point; S54, connect the first left bank buried depth point and the second left bank buried depth point to obtain the left bank toe line, connect the rightmost buried depth point of the riverbed and the leftmost buried depth point of the riverbed to obtain the riverbed toe line, extend the left bank toe line and the riverbed toe line to intersect them to obtain the junction point of the riverbed and the left bank toe line, and similarly obtain the junction point of the riverbed and the right bank toe line, both of which are the toe line junction points.

9. The method for automatically fitting the toe line of a concrete face rockfill dam according to claim 8, characterized in that: In S51, the rightmost buried depth point of the riverbed is specifically: The leftmost buried depth of the riverbed is: In the formula, C L X, C L Y is C L Coordinate, C R X, C R Y is C R Coordinate. k C , b C For straight line C L C R coefficient.

10. The method for automatically fitting the toe line of a concrete face rockfill dam according to claim 1, characterized in that: The S6 includes connecting the first left bank buried depth point, the second left bank buried depth point, the third left bank buried depth point and the fourth left bank buried depth point, the junction point of the riverbed and the left bank toe board line, the junction point of the riverbed and the right bank toe board line, the first right bank buried depth point, the second right bank buried depth point, the third right bank buried depth point and the fourth right bank buried depth point in S5 to form a complete panel dam toe board line.

Images