Earthwork calculation method and system

Abstract

The invention relates to the technical field of surveying and mapping, and discloses an earthwork calculation method and system, and the method comprises the following steps: opening a graphic file, selecting a field area and a ridge line, generating a plurality of fields through a left-turn algorithm, and marking the number, area and elevation point of each field; generating a natural terrain surface according to all elevation points on the drawing surface; enabling the user to input the total net volume of the field block area, and calculating the net volume of each field block in combination with the input total net volume of the field block area and the proportion of each field block in the total area; judging the position relationship between the natural terrain surface and the field, and generating a closed polyline conforming to the terrain; and on the basis of the closed polylines, calculating by adopting an earthwork calculation method to obtain the design elevation of each field parcel, and completing earth-rock calculation. The method solves the problems of large calculation amount and high cost of the existing earth-rock calculation technology, and has the characteristic of being capable of accurately reflecting local topographic features.

Description

Technical Field

[0001] The present invention relates to the technical field of mapping and drawing, and more specifically, to an earthwork calculation method and system. Background Art

[0002] Earthwork calculation is an important task in land consolidation, which directly affects the calculation of project quantities and costs. This article mainly discusses the earthwork calculation problem of terraced fields in mountainous and hilly areas. In mountain terraced fields, each field block needs to be calculated separately and finally summarized and statistically analyzed.

[0003] The prior art is an earthwork measurement and construction system, including a true three-dimensional engineering model establishment module for generating a true three-dimensional engineering model; an excavation plan module for formulating or adjusting an excavation plan for earthwork; an earthwork volume measurement module for measuring personnel to measure the earthwork volume excavated according to a set time; a three-dimensional engineering model display module for combining the excavation plan module and the earthwork volume measurement module in the true three-dimensional engineering model establishment module; a construction plan module for providing construction plan information for construction personnel according to the information in the three-dimensional engineering model display module; and a construction animation module for real-time simulation display of the driving path of transport vehicles and the excavation volume.

[0004] However, the existing earthwork calculation technology has problems of large calculation amount and high cost. Therefore, how to invent an earthwork calculation method that can reduce the calculation amount is a technical problem that urgently needs to be solved in this technical field. Summary of the Invention

[0005] In order to solve the problems of large calculation amount and high cost of the existing earthwork calculation technology, the present invention provides an earthwork calculation method and system, which has the characteristic of being able to accurately reflect local terrain features.

[0006] To achieve the above object of the present invention, the following technical solutions are adopted: An earthwork calculation method includes the following steps: Open a graphic file, select a field block area and a ridge line, generate multiple field blocks through a left-turn algorithm, and mark the number, area, and elevation points of each field block; Generate a natural terrain surface according to all elevation points on the drawing; Let the user input the total net volume of the field block area, and calculate the net volume of each field block by combining the input total net volume of the field block area and the proportion of each field block in the total area; Judge the positional relationship between the natural terrain surface and the field block, and generate a closed polyline that conforms to the terrain; Based on the closed polyline, calculate the design elevation of each field block by using an earthwork calculation method, and the earthwork calculation is completed.

[0007] Preferably, when generating the natural terrain surface, the Delaunay triangulation construction algorithm is specifically used to divide the entire calculation area into multiple triangular pyramids, obtain the point set of each triangulation network, and generate the natural terrain surface M based on the point set and draw it on the drawing surface.

[0008] Furthermore, calculate the net volume of each field plot, specifically: t1 = T * (S1 / S) Where, t1 is the net volume of each field plot, T is the total net volume, S1 is the area of each field plot, and S is the total area.

[0009] Even further, the specific steps for determining the positional relationship between the natural terrain surface and the field plot and generating a closed polyline that conforms to the terrain are as follows: According to the existing natural terrain surface M, from M1 to M n , determine the positional relationship between M n and the field plot P1: Record the natural terrain surface M that intersects and is contained within the field plot P1 s , use M s and P1 to generate a closed polyline M that conforms to the terrain using the left-turn algorithm sp , determine whether the center point of the polyline M sp is within the field plot P1. If it is, draw it on the drawing surface according to the point set. If not, it means the closed line is outside the field plot P1 and does not need to be drawn; Repeat the drawing until the positional relationship between all field plots P n and the natural terrain surface M n is determined.

[0010] Even further, calculate the design elevation based on the closed polyline, and the specific steps are as follows: According to the generated M sp , obtain the two-dimensional point set pt2dArr of each polyline through the secondary development interface; Extract the Z values of all points in pt2dArr and store them as pt3dArr. At the same time, record the maximum elevation value dMaxZ and the minimum elevation value dMinZ of each natural terrain surface M; Based on pt3dArr, split M sp into at least one triangle M spd ; Based on M spd , dMaxZ, and dMinZ, use the earthwork calculation method Fun() to calculate the design elevation dBg.

[0011] Further, to extract the Z values of all points in pt2dArr, the specific steps are as follows: Take any point coordinate pt1 in pt2dArr and determine the relationship between the point and the closed line: Find out which natural terrain surface M the point pt1 is inside. If it is inside, according to the set of point coordinates (XYZ) of the natural terrain surface M, use the three-point interpolation method to calculate the Z value of this point pt1; Repeat this step to calculate the Z values of all points in pt2dArr and store them as pt3dArr.

[0012] Further, based on pt3dArr, split M sp into at least one triangle M spd , and the specific steps are as follows: Determine the number of points in pt3dArr. If the number of points is greater than 3, split the concave and convex polygon M sp into multiple triangles M spd , and if it is equal to 3, no processing is required.

[0013] Further, based on M spd , dMaxZ, and dMinZ, use the earthwork calculation method Fun() to calculate the designed elevation dBg. The specific steps are as follows: S1: In each M spd , set the initial designed elevation as dBg = (dMaxZ + dMinz) / 2, and calculate the current earthwork volume t p ; If t p is equal to t1, then dBg is the required designed elevation, and the dichotomy ends; S2: If t p is less than t1, then dMinZ = dBg, and continue to execute S1; S3: If t p is less than t1, then dMaxZ = dBg, and continue to execute S1.

[0014] Further, the earthwork calculation method Fun() is specifically as follows: Use '+' to represent filling and '-' to represent excavation, and mark the triangles; Calculate the volume of each irregular triangle; Respectively count the total volume of the shapes with volume '+' and the total volume of the shapes with volume '-', and obtain the total excavation volume and the total filling volume.

[0015] An earthwork calculation system includes the following modules: A graphics processing module, which is used to select the field area and the ridge line, generate multiple fields through the left-turn algorithm, and mark the number, area, and elevation points of each field; A terrain modeling module, which generates a natural terrain surface according to all the elevation points on the drawing; The net volume distribution module is used to enable the user to input the total net volume of the field area, and calculate the net volume of each field by combining the input total net volume of the field area and the proportion of each field in the total area. The designed elevation calculation module is used to determine the positional relationship between the natural terrain surface and the fields, and generate a closed polyline that conforms to the terrain; based on the closed polyline, the designed elevation of each field is calculated using the earthwork calculation method, and the earthwork calculation is completed.

[0016] The beneficial effects of the present invention are as follows: The earthwork calculation method of the present invention realizes field division, natural terrain surface generation and earthwork calculation through an automated tool, greatly improving work efficiency. Based on the triangulation method and the earthwork calculation method, the accuracy of the calculation results is ensured, and it is especially suitable for mountainous and hilly areas with discontinuous terrain changes, and can meet the planning requirements of areas with discontinuous changes such as terraced fields. Description of the Drawings

[0017] Figure 1 It is a schematic flow chart of an earthwork calculation method of the present invention.

[0018] Figure 2 It is the graphic file imported into the SouthMap software in Example 2.

[0019] Figure 3 It is a schematic diagram of generating multiple fields in Example 2.

[0020] Figure 4 It is a schematic diagram of generating a natural terrain surface according to elevation points in Example 2.

[0021] Figure 5 It is a schematic diagram of generating a closed polyline that conforms to the terrain in Example 2.

[0022] Figure 6 It is a schematic diagram of the designed elevation calculation result in Example 2.

[0023] Figure 7 It is a schematic diagram of the earthwork calculation details table output in Example 2. Detailed Embodiments

[0024] The present invention will be described in detail below with reference to the drawings and specific embodiments.

[0025] Example 1 As Figure 1 shown, an earthwork calculation method includes the following steps: Open the graphic file, select the field area and the field ridge line, generate multiple fields through the left-turn algorithm, and mark the number, area, and elevation points of each field; Generate a natural terrain surface based on all elevation points on the drawing; Prompt the user to input the total net volume of the field area, and calculate the net volume of each field by combining the input total net volume of the field area and the proportion of each field in the total area; Determine the positional relationship between the natural terrain surface and the fields, and generate a closed polyline that conforms to the terrain; Based on the closed polyline, calculate the design elevation of each field using the earthwork calculation method, and the earthwork calculation is completed.

[0026] Embodiment 2 In this embodiment, the earthwork calculation method is implemented in the SouthMap software.

[0027] In this embodiment, the opened graphic file is as Figure 2 shown.

[0028] In this embodiment, first select the field area Y, then select the ridge line L, and use the left-turn algorithm to perform the minimum surface construction on the area Y and the ridge line L together to generate multiple fields P, and label the number, area, and elevation points of each field, as Figure 3 shown.

[0029] In this embodiment, as Figure 4 shown, when generating the natural terrain surface, specifically use the Delaunay triangulation algorithm to divide the entire calculation area into multiple triangular pyramids, obtain the point set of each triangular mesh, and generate the natural terrain surface M for the corresponding through the point set and draw it on the drawing.

[0030] In a specific embodiment, calculating the net volume of each field is specifically: t1 = T * (S1 / S) where t1 is the net volume of each field, T is the total net volume, S1 is the area of each field, and S is the total area.

[0031] In a specific embodiment, the specific steps for determining the positional relationship between the natural terrain surface and the fields and generating a closed polyline that conforms to the terrain are: According to the existing natural terrain surface M, from M1 to M n , determine the positional relationship between M n and the field P1: Record the natural terrain surface M that intersects and is contained inside the field P1 s , and use M s and P1 to construct a surface using the left-turn algorithm to generate a closed polyline M sp that conforms to the terrain, and determine whether the center point of the polyline M sp is inside the field P1. If it is, draw it on the drawing according to the point set. If not, it means the closed line is outside the field P1 and does not need to be drawn; Redraw repeatedly until all the fields P n and the natural terrain surface M n The positional relationship has been determined.

[0032] In this embodiment, the generated closed polyline M sp such as Figure 5 shown by the white polyline in

[0033] In a specific embodiment, calculating the design elevation based on the closed polyline, the specific steps are as follows: According to the generated M sp , obtain the two-dimensional point set pt2dArr of each polyline through the secondary development interface; Extract the Z values of all the points in pt2dArr, store and record them as pt3dArr. At the same time, record the maximum elevation value dMaxZ and the minimum elevation value dMinZ of each natural terrain surface M; Based on pt3dArr, split M sp into at least one triangle M spd ; Based on M spd , dMaxZ, dMinZ, use the earthwork calculation method Fun() to obtain the design elevation dBg.

[0034] In a specific embodiment, the steps of extracting the Z values of all the points in pt2dArr are as follows: Take any point coordinate pt1 in pt2dArr, judge the relationship between the point and the closed line: find out which natural terrain surface M the point pt1 is inside, if so, according to the set of point coordinates (XYZ) of the natural terrain surface M, use the three-point interpolation method to find the Z value of this point pt1; repeat this step to find the Z values of all the points in pt2dArr, store and record them as pt3dArr.

[0035] In a specific embodiment, based on pt3dArr, split M sp into at least one triangle M spd , the specific steps are as follows: Judge the number of points in pt3dArr, if the number of points is greater than 3, then split the concave and convex polygon M sp into multiple triangles M spd , if it is equal to 3, no processing is required.

[0036] In a specific embodiment, based on M spd , dMaxZ, dMinZ, use the earthwork calculation method Fun() to obtain the design elevation dBg, the specific steps are as follows: S1: In each M spd , set the initial design elevation as dBg = (dMaxZ + dMinz) / 2, and calculate the current earthwork volume t p; If t p is equal to t1, then dBg is the required design elevation, and the dichotomy method ends; S2: If t p is less than t1, then dMinZ = dBg, and continue to execute S1; S3: If t p is less than t1, then dMaxZ = dBg, and continue to execute S1.

[0037] In a specific embodiment, the earthwork calculation method Fun() is specifically as follows: Use \"+\" to represent filling and \"-\" to represent excavation, and mark triangles; Calculate the volume of each irregular triangle; Respectively count the total volume of the forms with volume \"+\" and volume \"-\", and obtain the total excavation volume and total filling volume.

[0038] In this embodiment, the calculated design elevation dBg is also marked on the drawing in sequence, such as Figure 6 shown in the red characters. Each field is marked with its field block number, area, and elevation. Then, click a point on the drawing according to the command line prompt information; combine the field block number, area, elevation point, and net volume, and use the AcDbTable entity class to mark the earthwork calculation details table on the drawing according to the self-set template. Store and output the earthwork calculation results of each field block obtained during the calculation process according to the classification.

[0039] In this embodiment, the field block method of the present invention is applicable to terrain conditions where the terrain changes discontinuously, such as earthwork calculations for regular planning of farmland terraces, etc. The existing elevation point data, plot boundary line and other information should be input first. First, divide the plot into corresponding field blocks according to the actual plot boundary line and ridge line. Based on the input net volume and measured elevation points, calculate the design elevation and earthwork volume of each field block in reverse, and can quickly and accurately obtain the regular information of each field block, and organize it in tabular form for convenient statistics and management. Embodiment 3 An earthwork calculation system includes the following modules: A graphic processing module, which is used to select the field block area and the ridge line, generate multiple field blocks through the left-turn algorithm, and mark the number, area, and elevation point of each field block; A terrain modeling module, which generates a natural terrain surface according to all the elevation points on the drawing; A net volume distribution module, which is used to enable the user to input the total net volume of the field block area, and calculate the net volume of each field block in combination with the input total net volume of the field block area and the proportion of each field block in the total area; The design elevation calculation module is used to judge the positional relationship between the natural terrain surface and the field block, and generate a closed polyline that conforms to the terrain; based on the closed polyline, the design elevation of each field block is calculated by using the earthwork calculation method, and the earthwork calculation is completed; The result output module is used to output the earthwork calculation result.

[0040] Obviously, the above-mentioned embodiments of the present invention are merely examples for clearly illustrating the present invention, rather than limitations on the implementation manners of the present invention. Any modifications, equivalent replacements, improvements, etc. made within the spirit and principle of the present invention shall be included in the protection scope of the claims of the present invention.

Claims

1. A method for calculating earthwork and stonework, characterized in that, It includes the following steps: Open the graphic file, select the field area and the ridge line, generate multiple fields through the left-turn algorithm, and mark the number, area, and elevation points of each field; Generate the natural terrain surface based on all the elevation points on the drawing; Let the user input the total net volume of the field area, and combine the input total net volume of the field area and the proportion of each field in the total area to calculate the net volume of each field; Judge the positional relationship between the natural terrain surface and the field, and generate a closed polyline that conforms to the terrain; Based on the closed polyline, use the earthwork calculation method to calculate the design elevation of each field, and the earthwork calculation is completed.

2. The earthwork calculation method according to claim 1, characterized in that, When generating the natural terrain surface, specifically use the Delaunay triangulation construction algorithm to divide the entire calculation area into multiple triangular pyramids, obtain the point set of each triangulation network, and generate the natural terrain surface M for each point set and draw it on the drawing.

3. The earthwork calculation method according to claim 2, wherein, Calculate the net volume of each field, specifically: t1 = T * (S1 / S) Where t1 is the net volume of each field, T is the total net volume, S1 is the area of each field, and S is the total area.

4. The earthwork calculation method according to claim 3, characterized in that, The specific steps for judging the positional relationship between the natural terrain surface and the field and generating a closed polyline that conforms to the terrain are: According to the existing natural terrain surface M, from M1 to M n , judge the positional relationship between M n and the field block P1: Record the natural terrain surfaces M that intersect and are contained within the field block P1 s , and use M s and P1 to generate a terrain-conforming closed polyline M using the left-turn algorithm sp , and judge whether the center point of the polyline M sp is inside the field block P1. If it is, draw it on the drawing according to the point set. If not, it means the closed line is outside the field block P1 and does not need to be drawn; Redraw repeatedly until all the fields P n and the natural terrain surface M n The positional relationship has been determined.

5. The earthwork calculation method according to claim 4, wherein Calculate the design elevation based on the closed polyline, and its specific steps are: According to the generated M sp , obtain the 2D point set pt2dArr of each polyline through the secondary development interface; Extract the Z values of all points in pt2dArr, store them and record them as pt3dArr. At the same time, record the maximum elevation value dMaxZ and the minimum elevation value dMinZ of each natural terrain surface M; Split M based on pt3dArr sp into at least one triangle M spd ; Based on M spd , dMaxZ, dMinZ, the earthwork calculation method Fun() is used to obtain the designed elevation dBg.

6. The earthwork calculation method according to claim 5, characterized in that Extract the Z values of all points in pt2dArr. The specific steps are: take any point coordinate pt1 in pt2dArr, judge the relationship between the point and the closed line: find out which natural terrain surface M the pt1 is inside, if it is, then according to the point coordinate (XYZ) set of the natural terrain surface M, use the three-point interpolation method to find the Z value of this point pt1; repeat this step to find the Z values of all points in pt2dArr, and store them and record them as pt3dArr.

7. The earthwork calculation method according to claim 5, wherein Split M based on pt3dArr sp into at least one triangle M spd , and the specific steps are as follows: judge the number of points in pt3dArr. If the number of points is greater than 3, split the concave-convex polygon M sp into multiple triangles M spd . If it is equal to 3, no processing is required.

8. The earthwork calculation method according to claim 7, wherein Based on M spd , dMaxZ, dMinZ, the design elevation dBg is obtained by using the earthwork calculation method Fun(). The specific steps are as follows: S1: At each M spd , set the initial design elevation as dBg = (dMaxZ + dMinz) / 2, and use the earthwork calculation method Fun() to calculate the current earthwork volume t p ; If t p is equal to t1, then dBg is the required design elevation, and the bisection method ends; S2: If t p is less than t1, then dMinZ = dBg, and continue to execute S1; S3: If t p is less than t1, then dMaxZ = dBg, and continue to execute S1.

9. The earthwork calculation method according to claim 8, characterized in that, The specific earthwork calculation method Fun() is: Use "+" to represent filling and "-" to represent excavation, and mark the triangles; Calculate the volume of each irregular triangle; Respectively count the total volume of the shapes with a volume of "+" and a volume of "-" to obtain the total excavation volume and the total filling volume.

10. An earthwork calculation system, characterized in that, It includes the following modules: Graphic processing module, which is used to select the field area and the ridge line, generate multiple fields through the left-turn algorithm, and mark the number, area, and elevation points of each field; Terrain modeling module, which generates the natural terrain surface based on all the elevation points on the drawing; Net volume distribution module, which is used to let the user input the total net volume of the field area, and combine the input total net volume of the field area and the proportion of each field in the total area to calculate the net volume of each field; Design elevation calculation module, which is used to judge the positional relationship between the natural terrain surface and the field and generate a closed polyline that conforms to the terrain; Based on the closed polyline, use the earthwork calculation method to calculate the design elevation of each field, and the earthwork calculation is completed.

Images