1. A cross section mapping method in hydraulic engineering measurement is characterized in that: the method comprises the following steps:

s1, importing design route data and automatically drawing a design route;

s2, configuring cross section attributes; the cross-sectional attributes include cross-sectional spacing and cross-sectional length;

s3, sequentially calculating the coordinates of the center point of the section line of each cross section and the coordinates of the left and right base points of the section line according to the cross section distance;

s4, automatically generating a cross section base point file of the RTK identification;

s5, importing data of each cross section point measured by field work;

s6, constructing a triangle by each section point and the left and right base points of each section line respectively, making a perpendicular line to the section lines by taking the section point as the base point, and calculating the coordinates of the foot hanging point, the distance D from the foot hanging point to the center point of the section line and the distance L from the foot hanging point to the section point;

s7, configuring a section line threshold value, and screening out measurement data near each section according to the section line threshold value;

s8, sorting the measured data near each cross section and the distance D from the foot point corresponding to the cross section point represented by the measured data obtained by calculation in the step S6 to the central point of the cross section line into a new data set F according to the data organization format according to the screening result in the step S7;

s9, sorting according to the distance D from the foot drop point to the central point of the cross section line from small to large to form a new data set F';

and S10, sequentially drawing the measurement data in the data set F 'to form a cross-sectional diagram, and generating a cross-sectional data table from the data set F'.

2. A cross-sectional mapping method in hydraulic engineering survey according to claim 1, characterized in that: in step S3, the calculation process calculates the coordinates of the center point of the cross section line and the coordinates of the left and right base points of the cross section line by repeatedly using the mathematical formulas of forward coordinate calculation and backward coordinate calculation.

3. A cross-sectional mapping method in hydraulic engineering survey according to claim 1, characterized in that: in S6, the calculating repeatedly calculates the coordinates of the foot point, the distance D from the foot point to the center point of the cross-sectional line, and the distance L from the foot point to the cross-sectional line using the mathematical formulas of coordinate forward calculation and coordinate backward calculation.

4. A cross-sectional mapping method in hydraulic engineering survey according to claim 1, characterized in that: in step S7, the process of screening out the measurement data near each cross section is to determine the relationship between the distance L from the foot drop point to the cross section point in step S6 and the threshold of the cross section line, that is: when L is less than or equal to the threshold of the section line, the section point is the measurement data near the section; when L is greater than the profile line threshold, then the profile point is not measurement data near the profile.

5. A cross-sectional mapping method in hydraulic engineering survey according to claim 1, characterized in that: in S10, the cross-sectional diagram and the cross-sectional data table are txt documents or documents that can be identified by cad drawing software.

Background

With the comprehensive popularization of a CORS (continuous operation satellite positioning service integrated system) network, the traditional method for using a total station instrument to measure the cross section is gradually eliminated. At present, RTK (abbreviation of English Real-time kinematic) is mainly used by field industry as an important means for section measurement; however, much work is required from field measurement data to a real sectional view: firstly, according to a known design line, generating section line data in batch, then exporting the section line data in a specific format, and performing field lofting by field measurement personnel according to the exported section line data. The section data measured by field workers in the field needs to be converted into two forms of section diagrams and reports in batch. The integration degree of section mapping software in the market is not high at present, a designed line is combined with a section chart generated by actually measured data, and software meeting the requirements of hydraulic engineering does not exist at present.

Disclosure of Invention

The invention aims to provide a cross section mapping method in hydraulic engineering measurement.

In order to achieve the purpose, the invention can adopt the following technical scheme:

the invention relates to a cross section mapping method in hydraulic engineering measurement, which comprises the following steps:

s1, importing design route data and automatically drawing a design route;

s2, configuring cross section attributes; the cross-sectional attributes include cross-sectional spacing and cross-sectional length;

s3, sequentially calculating the coordinates of the center point of the section line of each cross section and the coordinates of the left and right base points of the section line according to the cross section distance;

s4, automatically generating a cross section base point file of the RTK identification;

s5, importing data of each cross section point measured by field work;

s6, constructing a triangle by each section point and the left and right base points of each section line respectively, making a perpendicular line to the section lines by taking the section point as the base point, and calculating the coordinates of the foot hanging point, the distance D from the foot hanging point to the center point of the section line and the distance L from the foot hanging point to the section point;

s7, configuring a section line threshold value, and screening out measurement data near each section according to the section line threshold value;

s8, sorting the measured data near each cross section and the distance D from the foot point corresponding to the cross section point represented by the measured data obtained by calculation in the step S6 to the central point of the cross section line into a new data set F according to the data organization format according to the screening result in the step S7;

s9, sorting according to the distance D from the foot drop point to the central point of the cross section line from small to large to form a new data set F';

and S10, sequentially drawing the measurement data in the data set F 'to form a cross-sectional diagram, and generating a cross-sectional data table from the data set F'.

Preferably, in step S3, the calculating process calculates coordinates of the center point of the section line and coordinates of the left and right base points of the section line for each cross section by repeatedly using a mathematical formula of coordinate forward calculation and coordinate backward calculation.

Preferably, in S6, the calculating repeatedly calculates the coordinates of the drop foot point, the distance D from the drop foot point to the center point of the cross-sectional line, and the distance L from the drop foot point to the cross-sectional point by using the mathematical formulas of coordinate forward calculation and coordinate backward calculation.

Preferably, in step S7, the process of selecting the measurement data near each cross section is to determine the relationship between the distance L from the foot point to the cross section point in step S6 and the threshold of the cross section line, that is: when L is less than or equal to the threshold of the section line, the section point is the measurement data near the section; when L is greater than the profile line threshold, then the profile point is not measurement data near the profile.

Preferably, in S10, the cross-sectional diagram or the cross-sectional data table is a txt document or a document that can be identified by cad drawing software.

The invention has the advantages that cross section lines are automatically generated according to design line data, points meeting the point position precision requirement are intelligently selected as cross section points or longitudinal section points according to field measurement data, and longitudinal and cross section maps meeting the water conservancy requirement are automatically generated according to the cross section points, so that the automation of the cross section mapping is realized, the time cost is saved, and the working efficiency is improved.

Drawings

FIG. 1 is a flow chart of the method of the present invention.

FIG. 2 is a cross-sectional view produced by the method of the present invention.

FIG. 3 is a table of results generated by the method of the present invention.

FIG. 4 is a schematic representation of the process of the present invention.

Detailed Description

The technical solutions in the embodiments of the present invention will be described clearly and completely below, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.

As shown in fig. 1 and 4, the cross-sectional mapping method in hydraulic engineering measurement according to the present invention includes the following steps:

s1, importing design route data and automatically drawing a design route 1; the design route data is provided by a designer and is in the format of: serial number, X, Y, turning radius; wherein X and Y represent coordinate values;

for example, certain design line data is as follows:

1,3697394.274,473468.902,0.000

2,3708706.744,484701.259,1500.000

3,3706191.516,495454.010,800.000

4,3712680.396,506384.580,2000.000

5,3708591.133,518374.143,0.000

………………

wherein the first row represents the meaning: the point with the sequence number of 1 on the designed route has the X coordinate of 3697394.274, the Y coordinate of 473468.902 and the turning radius of 0;

s2, configuring cross section attributes; the cross-sectional attributes include cross-sectional spacing and cross-sectional length;

s3, calculating the coordinates of the cross section line center point 2 of the cross section and the coordinates of the left base point 3 and the right base point 4 of the cross section line according to the cross section distance; the whole design route 1 is provided with a plurality of cross sections according to the cross section intervals configured in the step S2, and each cross section is provided with a section line 5; the section line 5 is perpendicular to the design route 1 and is positioned in the same horizontal plane with the design route 1; the length of the section line 5 is equal to the length of the cross section configured in the step S2, the intersection point of the section line and the design route 1 is a section line central point 2, and the section line central point 2 equally divides the section line into two sections; wherein, according to the route advancing direction (the direction indicated by the dotted arrow in fig. 4 is the route advancing direction), the section line end point on the left side of the section line central point 2 is the left base point 3 of the section line; the end point of the section line on the right side of the central point 2 of the section line is a right base point 4 of the section line; the route advancing direction is the route direction indicated by the sequence number in the designed route data from small to large in the step S1;

according to the known point data, the cross section distance and the cross section length on the design route 1 in the step S1, calculating the coordinates of the left base point 3 and the right base point 4 of each cross section and the coordinates of the cross section line central point 2 of each cross section by repeatedly using a coordinate forward calculation and a coordinate backward calculation mathematical formula;

s4, automatically generating a cross section base point file for the southern RTK instrument to identify; the format of the cross section base point file is as follows: ZHZ, N, E, H, ZHY, N, E, H, starting point mileage, section number;

the serial number of the section is a serial number of the cross section and is an identification code of the cross section; the starting point mileage represents a starting value of a left base point 3 of a southern RTK instrument measuring section line 5, that is, when the value is 0, the measurement data represents the distance from the measuring instrument to the left base point 3; when the value is 10, the measured data represents the distance from the measuring instrument to the left base point 3 plus the calculated value 10; ZHZ represents the pile number 3 of the left base point; ZHY represents the number of 4 piles on the right base point; n and E represent the coordinate values of X and Y of the base point, respectively; h represents the elevation value of the base point;

for example, the specific data of a base point file is as follows:

0+000L,3844273.715,401663.8468,0,,0+000R,3844136.866,401809.6969,0,,0, 1

0+100L,3844311.368,401701.3737,0,, 0+100R,3844169.977,401842.8258,0,,0, 2

0+200L,3844346.731,401736.7213,0,, 0+200R,3844205.34,401878.1735,0,,0, 3

………………

wherein the first row represents the meaning: the starting mileage of a left base point 3 of the first cross section is 0, wherein the pile number of the left base point 3 is 0+000L, the X coordinate is 3844273.715, the Y coordinate is 401663.8468, the elevation H is 0, the pile number of the right base point 4 is 0+000R, the X coordinate is 3844136.866, the Y coordinate is 401809.6969, and the elevation H is 0;

wherein the second row represents the meaning: in the second cross section, the starting point mileage of the left base point 3 is 0, wherein the pile number of the left base point 3 is 0+100L, the X coordinate is 3844311.368, the Y coordinate is 401701.3737, the elevation H is 0, the pile number of the right base point 4 is 0+100R, the X coordinate is 3844169.977, the Y coordinate is 401842.8258, and the elevation H is 0;

s5, importing field measurement data, wherein the data import format is as follows: serial number, attribute of section point 6, Y, X and H;

for example, the imported field measurement data is specifically as follows:

1, waterside, 491919.651,3706209.009,30.050

2, waterside, 491927.131,3706234.953,31.050

3, right bank top, 491956.240,3706335.923,30.300

4, left bank angle, 491925.111,3706271.803,27.247

5 center point, 491924.088,3706269.796,27.331

……………………………………

Wherein the first row represents the meaning: the serial number of the section point 6 is 1, the attribute of the section point 6 is water edge, the X coordinate of the section point 6 is 491919.651, the Y coordinate is 3706209.009, and the elevation H is 30.050;

s6, constructing a triangle (as shown by the two-dot chain line in fig. 4) between each section point 6 and the corresponding left base point 3 and right base point 4 of the section line 5, making a perpendicular line to the section line 5 with the section point 6 as the base point, and repeatedly calculating the coordinates of the foot point 7, the distance D from the foot point 7 to the center point 2 of the section line, and the distance L from the foot point 7 to the section point 6 by using the mathematical formulas of coordinate forward calculation and coordinate backward calculation;

s7, configuring a section line threshold value; screening out measurement data near each section according to the section line threshold;

for example, when the section line threshold is configured to be 0.5 meter, and the measurement data near the first cross section is screened, firstly, each section point 6 and the left base point 3 and the right base point 4 of the first section line 5 form a triangle, and the coordinates of each section point 6 and the foot hanging point 7 of the section line 5, the distance D from the foot hanging point 7 to the section line central point 2, and the distance L from the foot hanging point 7 to the section point 6 are obtained; then judging the relation between L and the threshold value of the section line, namely when L is less than or equal to 0.5 m, the section point 6 is the measurement data near the first cross section; when L is more than 0.5 m, the section point 6 is not the measurement data near the first cross section;

repeating the process of screening the measurement data near the first cross section when screening the measurement data near the second and third cross sections;

s8, sorting the measured data near each cross section and the distance D from the foot point 7 corresponding to the cross section point 6 represented by the measured data obtained by calculation in the step S1.6 to the cross section line central point 2 into a new data set F according to the data organization format according to the screening result in the step S1.7; the data organization format is as follows: the distance D from the foot drop point 7 to the central point 2 of the section line and the attributes X, Y and H of the section point 6 in the measured data;

s9, sorting according to the distance D from the foot drop point 7 to the central point 2 of the section line from small to large to form a new data set F';

s10, sequentially drawing the measurement data in the data set F 'to form a cross-sectional diagram, as shown in FIG. 2, and generating a cross-sectional data table from the data set F', as shown in FIG. 3;

s11, sharing the result file to related workers, realizing result sharing and improving working efficiency; the result file is a txt document or a document which can be identified by cad drawing software.