A fast calculation method and system for earthwork volume in multi-scale excavation areas with differentiated lithology

Abstract

The invention discloses a method and system for quickly calculating the earthwork volume of a multi-scale excavation area for distinguishing lithology, and relates to the technical field of earthwork volume calculation; the method comprises: obtaining engineering data of an area to be excavated and selecting a contour line of a calculation area; performing spatial interpolation on point cloud data obtained based on the engineering data to obtain elevation grid data; based on the contour line of the calculation area, intercepting target elevation grid data in a target area, and generating a corresponding mask matrix to calculate and obtain the volume excavation volume of each geological layer; obtaining the excavation volume of each geological layer in each layer according to the mask matrix and the step height; obtaining the excavation volume of each geological layer in each strip according to the point cloud data, a preset strip spacing and a strip division reference line; obtaining the excavation volume of each geological layer in each block according to the point cloud data, a preset block spacing and a block division reference line; the invention improves the calculation efficiency while ensuring the accuracy, and quickly obtains the earthwork volume of each geological layer in different scales of volume, layer, strip and block.

Description

Technical Field

[0001] The present invention relates to the technical field of earthwork volume calculation, and in particular to a method and system for quickly calculating earthwork volume in a multi-scale excavation area for distinguishing lithology. Background Art

[0002] The calculation of earthwork volume is of great significance for the construction of large-scale civil engineering projects such as airports and rockfill dams. It is directly related to the optimization of schemes and the budget of engineering costs. Different filling areas have different requirements for the properties of excavated materials, so it is necessary to calculate the excavation volume of each geological layer. The commonly used method is to create terrain surfaces and geological layer surfaces in Civi13D software, use the Extract Entity from Surface command to build a three-dimensional geological model, and query the volume of each geological layer.

[0003] In the process of earthwork allocation, in order to better meet the progress requirements and material matching requirements, it is necessary to conduct a more fine-grained analysis of the excavation volume, that is, to count the excavation volume of each geological layer in each layer, each strip, and each block; however, when using Civi13D software for fine-grained analysis of excavation volume, multiple sectioning operations are required on the three-dimensional geological model, which is inefficient.

[0004] In the prior art, although there are related technical solutions for using Civil3D to calculate the volume of earth and stone by establishing a curved surface, or constructing a geological stratification model to count the volume of each geological layer, none of them has further realized the calculation of the volume of earth and stone for each geological layer in different scales of body-layer-strip-block; for example, inventions with application numbers CN202010059628.4, CN202110009012.0, and CN201811440239.5. Summary of the invention

[0005] In order to solve the problems existing in the prior art, the present invention provides a method for quickly calculating the earthwork volume of multi-scale excavation areas that distinguishes rock types, which can improve the calculation efficiency and quickly obtain the earthwork volume of each geological layer in different scales of body, layer, strip and block to meet the construction period requirements of large-scale projects.

[0006] To achieve the above object, the present invention provides a method for quickly calculating earthwork volume in a multi-scale excavation area for distinguishing lithology, comprising:

[0007] Acquire engineering data of the area to be excavated and calculate the area contour; the engineering data includes: original terrain data, design surface data and exploration hole data;

[0008] Acquire corresponding point cloud data according to the engineering data, and perform spatial interpolation based on the point cloud data to obtain corresponding elevation grid data; the point cloud data includes: original terrain point cloud, design surface point cloud, and thickness point cloud of each geological layer;

[0009] Based on the calculation area contour, determine the target area in the elevation grid data, intercept the elevation grid data in the target area to obtain target elevation grid data, and generate a mask matrix corresponding to the target elevation grid data;

[0010] Based on the grid volume calculation principle, matrix operations are performed on the target elevation grid data to obtain the volume matrix of each geological layer;

[0011] According to the mask matrix and the volume matrix of each geological layer, the volume excavation volume of each geological layer is obtained; the volume excavation volume of each geological layer is used to characterize the volume of each geological layer within the contour line of the calculation area;

[0012] Using the mask matrix and the step height, the elevation grid data is processed in layers to obtain a plurality of layered elevation grid data, and matrix operations are performed on the elevation grid data of each layer to obtain a volume matrix of each geological layer in each layer, and then the excavation volume of each geological layer in each layer is obtained according to the mask matrix and the volume matrix of each geological layer in each layer;

[0013] The layered elevation grid data is processed into strips using point cloud data, preset strip spacing and strip division reference lines to obtain multiple strips of elevation grid data and corresponding multiple strip mask matrices, and matrix operations are performed on each strip of elevation grid data to obtain a volume matrix of each geological layer in each strip, and then the excavation volume of each geological layer in each strip is obtained according to each strip mask matrix and the corresponding volume matrix of each geological layer in each strip;

[0014] The strip elevation grid data is divided into blocks using point cloud data, preset block spacing and block division reference lines to obtain multiple block mask matrices, and the excavation volume of each geological layer in each block is obtained based on the volume matrix of each geological layer in each strip and the corresponding block mask matrix.

[0015] Optionally, the elevation grid data includes: original terrain elevation grid data, design surface elevation grid data, and interface elevation grid data of adjacent geological layers; the spatial interpolation based on the point cloud data to obtain the corresponding elevation grid data specifically includes:

[0016] Performing spatial interpolation on the original terrain point cloud to obtain corresponding original terrain elevation grid data;

[0017] Performing spatial interpolation on the design surface point cloud to obtain corresponding design surface elevation grid data;

[0018] Perform spatial interpolation on the thickness point cloud of each geological layer to obtain the corresponding grid data of the thickness of each geological layer;

[0019] The interface elevation grid data of adjacent geological layers are obtained based on the original terrain elevation grid data and the thickness grid data of each geological layer.

[0020] Optionally, determining the target area in the elevation grid data based on the calculation area contour line includes:

[0021] Extracting contour line point data in the contour line of the calculation area;

[0022] Obtaining the maximum and minimum values ​​of the x-coordinate and the maximum and minimum values ​​of the y-coordinate in the contour line point data;

[0023] Determine the minimum circumscribed rectangular frame of the calculation area contour line in the elevation grid data according to the maximum and minimum values ​​of the x-coordinate and the maximum and minimum values ​​of the y-coordinate;

[0024] The four sides of the minimum circumscribed rectangular frame are expanded according to a preset distance to obtain a target area.

[0025] Optionally, the value in the mask matrix is ​​0 or 1; wherein 0 is used to represent point data in the target area that is outside the contour line of the calculation area, and 1 is used to represent point data in the target area that is inside the contour line of the calculation area.

[0026] Optionally, obtaining the volume excavation volume of each geological layer according to the mask matrix and the volume matrix of each geological layer includes:

[0027] Multiplying the volume matrix of each geological layer by the mask matrix respectively to obtain the multiplication result of each geological layer;

[0028] The sum of the positive values ​​of the multiplication results is calculated respectively to obtain the volume excavation of each geological layer.

[0029] Optionally, the target elevation grid data is subjected to hierarchical processing using the mask matrix and the step height to obtain a plurality of hierarchical elevation grid data, including:

[0030] Based on the target elevation grid data and the mask matrix, obtaining the highest point of the terrain and the lowest point of the design surface within the contour line of the calculation area;

[0031] Determine the layering parameters according to the step height, the highest point of the terrain and the lowest point of the design surface; the layering parameters include the number of layers of the target elevation grid data, the top elevation of each layer, and the bottom elevation of each layer;

[0032] According to the layered parameters, a correction operation is performed on the target elevation grid data to obtain a plurality of layered elevation grid data;

[0033] The performing of the correction operation comprises:

[0034] Based on the number of layers, the original terrain elevation grid data and interface elevation grid data in each layer of the target elevation grid data that are higher than the layer top elevation are corrected to the layer top elevation, and the design surface elevation grid data in each layer of the target elevation grid data that are lower than the layer bottom elevation are corrected to the layer bottom elevation.

[0035] Optionally, the layered elevation grid data is processed into strips using the point cloud data, a preset strip spacing and a strip division reference line to obtain a plurality of strip elevation grid data, including:

[0036] Based on the point cloud data, a target point cloud corresponding to a portion of the layered elevation grid data where the original terrain elevation grid data is higher than the design surface elevation grid data is obtained;

[0037] Determining a plurality of stripe interception areas according to the target point cloud, a preset stripe spacing and a stripe division reference line;

[0038] Based on the plurality of intercepted areas, the layered elevation grid data is intercepted to obtain a corresponding plurality of strips of elevation grid data.

[0039] Optionally, determining a plurality of stripe interception areas according to the target point cloud, a preset stripe spacing and a stripe division reference line includes:

[0040] The data point closest to the dividing reference line in the target point cloud is taken as the first target point, and the data point farthest from the dividing reference line is taken as the second target point;

[0041] Determine a strip division start line passing through the first target point, and a strip division end line passing through the second target point; the strip division start line and the strip division end line are both parallel to the strip division reference line;

[0042] Determine a normal line of the partition reference line that is outside the range of the target point cloud;

[0043] The data point closest to the normal line in the target point cloud is used as the third target point, and the data point farthest from the normal line is used as the fourth target point;

[0044] Determine a first boundary line passing through the third target point, and a second boundary line passing through the fourth target point; the first boundary line and the second boundary line are both parallel to the normal line;

[0045] Calculating the interval distance between the strip division starting line and the strip division ending line;

[0046] Dividing the interval distance by the preset strip spacing to obtain the number of strip divisions;

[0047] Dividing a rectangular area enclosed by the strip division start line, the strip division end line, the first boundary line and the second boundary line into a plurality of strip areas according to the strip division number;

[0048] The minimum circumscribed rectangular frame of each of the strip regions is determined, and the obtained multiple minimum circumscribed rectangular frames are used as corresponding multiple strip interception regions.

[0049] According to the specific embodiments provided by the present invention, the present invention discloses the following technical effects:

[0050] The method for quickly calculating the earthwork volume of multi-scale excavation areas that distinguishes rock types provided by the present invention can quickly obtain the earthwork volume of each geological layer (different rock types) in different scales of body, layer, strip and block. Compared with the use of Civi13D software for sectioning analysis in the prior art, it has the advantages of less time and higher efficiency, especially when conducting excavation engineering volume analysis of fine-grained layers, strips and blocks, the time consumed is greatly reduced; at the same time, the calculation results obtained by the calculation method provided by the present invention are slightly different from the calculation results obtained using Civi13D, can meet the engineering requirements for accuracy, and can be directly applied to actual engineering. BRIEF DESCRIPTION OF THE DRAWINGS

[0051] The above and other objects, features and advantages of the present invention will become more apparent through a more detailed description of exemplary embodiments of the present invention in conjunction with the accompanying drawings, wherein like reference numerals generally represent like components throughout the exemplary embodiments of the present invention.

[0052] Figure 1 A schematic diagram of a method flow chart of a method for rapidly calculating earthwork volume in a multi-scale excavation area for distinguishing lithology according to an embodiment of the present invention;

[0053] Figure 2 A schematic diagram of three-dimensional and XOY projection of a point cloud before and after Kriging interpolation according to an embodiment of the present invention;

[0054] Figure 3 A schematic diagram of a matrix format for storing elevation grid data according to an embodiment of the present invention;

[0055] Figure 4 A schematic diagram of intercepting elevation grid data of a target area shown in an embodiment of the present invention;

[0056] Figure 5 A schematic diagram of a process for calculating the volume of each geological layer of a prism using a grid method according to an embodiment of the present invention;

[0057] Figure 6 This is a schematic diagram of the calculation process of the volume of "body" excavation according to an embodiment of the present invention;

[0058] Figure 7 This is a schematic diagram of the calculation process of the “layer” excavation volume according to an embodiment of the present invention;

[0059] Figure 8 A schematic diagram of correction of the i-th layer elevation grid data shown in an embodiment of the present invention;

[0060] Fig. 9 This is a schematic diagram of the calculation process of the "strip" excavation volume according to an embodiment of the present invention;

[0061] Fig.10 A schematic diagram of determining the interception range of the j-th elevation grid data of the i-th layer according to an embodiment of the present invention;

[0062] Fig.11 A schematic diagram of the calculation process of the "block" excavation volume according to an embodiment of the present invention;

[0063] Fig.12 This is a schematic diagram of the module structure of a system for rapidly calculating earthwork volume in a multi-scale excavation area for distinguishing lithology according to an embodiment of the present invention. DETAILED DESCRIPTION

[0064] The following will be combined with the drawings in the embodiments of the present invention to clearly and completely describe the technical solutions in the embodiments of the present invention. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0065] The purpose of the present invention is to provide a method for quickly calculating the earthwork volume of a multi-scale excavation area that distinguishes lithology, and can quickly calculate the earthwork volume of each geological layer in different scales such as body, layer, strip and block.

[0066] See also Figure 1 , Figure 1 Schematic diagram of the method flow for rapid calculation of earthwork volume in multi-scale excavation areas for distinguishing lithology.

[0067] The above earthwork volume calculation method includes:

[0068] Step 101: Obtain engineering data of the area to be excavated and calculate the area contour.

[0069] Among them, engineering data includes: original terrain data, design surface data and exploration hole data. Original terrain data can be obtained by field measurement using surveying and mapping tools, such as scanning the area to be excavated by a high-precision laser scanner or a 3D camera carried by a drone, so as to obtain original terrain data; design surface data can be obtained according to actual engineering needs. The design surface in the calculation of earthwork volume refers to the calculation reference surface of earthwork volume determined according to engineering needs during the design stage, which represents the expected ground morphology or structure in the engineering design; exploration hole data of the area to be excavated can be obtained through geological exploration means.

[0070] The original terrain data and the designed surface data have the same storage structure, and each point data consists of x coordinates, y coordinates and elevation value z. The exploration hole data consists of the exploration point number, x coordinates, y coordinates and multiple geological layer thicknesses.

[0071] Step 102: Obtain corresponding point cloud data according to the engineering data, and perform spatial interpolation based on the point cloud data to obtain corresponding elevation grid data.

[0072] The point cloud data includes: original terrain point cloud, design surface point cloud, and thickness point cloud of each geological layer. In the application, the original terrain point cloud, design surface point cloud, and thickness point cloud of each geological layer can be obtained according to the original terrain data, design surface data, and exploration hole data of the area to be excavated.

[0073] After obtaining the point cloud data, you can use methods such as inverse distance weighted interpolation (IDW), Kriging interpolation (Kriging), natural neighbor interpolation (Natural Neighbor) or radial basis function interpolation (RBF) to perform spatial interpolation on the point cloud data; taking Kriging interpolation as an example, see Figure 2 , Figure 2 Schematic diagram of the 3D and XOY projection of the point cloud before and after Kriging interpolation.

[0074] In the application, the point cloud data can be spatially interpolated according to the preset grid spacing. The preset grid spacing can be flexibly set according to the actual situation. Figure 4 , Figure 4 , the elevation grid data obtained after interpolation according to the 25m grid spacing (preset grid spacing) is shown.

[0075] Corresponding to the point cloud data, the elevation grid data includes the original terrain elevation grid data, the design surface elevation grid data, and the interface elevation grid data of adjacent geological layers; specifically, the corresponding elevation grid data is obtained by spatial interpolation based on the point cloud data, including:

[0076] Perform spatial interpolation on the original terrain point cloud to obtain the corresponding original terrain elevation grid data;

[0077] Perform spatial interpolation on the design surface point cloud to obtain the corresponding design surface elevation grid data;

[0078] Perform spatial interpolation on the thickness point cloud of each geological layer to obtain the corresponding grid data of the thickness of each geological layer;

[0079] The interface elevation grid data of adjacent geological layers are obtained based on the original terrain elevation grid data and the thickness grid data of each geological layer.

[0080] After the elevation grid data is obtained, it can be saved in matrix format to facilitate subsequent matrix operations; Figure 3 As shown, Figure 3 Diagram for saving matrix format for elevation grid data.

[0081] Step 103: Based on the calculated area contour, determine the target area in the elevation grid data, intercept the elevation grid data in the target area to obtain target elevation grid data, and generate a mask matrix corresponding to the target elevation grid data.

[0082] The calculation area contour line is to determine and depict the boundary line of the area to be excavated, which is a closed curve data; the calculation area contour line can be flexibly set according to the actual engineering needs. In the application, the x-coordinate and y-coordinate of the calculation area contour point data can be extracted, and the target area can be determined according to the size range of the x-coordinate and y-coordinate.

[0083] Exemplarily, the above-mentioned determination of the target area in the elevation grid data based on the calculated area contour line may include the following steps:

[0084] Extracting contour point data from contour lines of the calculation area;

[0085] Get the maximum and minimum values ​​of the x-coordinate and the maximum and minimum values ​​of the y-coordinate in the contour line point data;

[0086] According to the maximum and minimum values ​​of the x-coordinate and the maximum and minimum values ​​of the y-coordinate, the minimum circumscribed rectangular frame of the contour line of the calculation area in the elevation grid data is determined;

[0087] The four sides of the minimum bounding rectangle are expanded according to the preset distance to obtain the target area.

[0088] Among them, the preset distance can be flexibly adjusted according to actual conditions, for example, it can be expanded by 30m.

[0089] See further Figure 4 , Figure 4 It is a schematic diagram for intercepting the elevation grid data of the target area; Figure 4The interception range in the target area is also the target area. In the application, the target area is determined and the elevation grid data in the target area is intercepted to reduce the calculation range and improve the calculation efficiency.

[0090] After intercepting the target elevation grid data, a corresponding mask matrix can be made according to the contour line of the calculation area; specifically, a mask matrix composed of 0 and 1 can be generated; wherein 0 represents the point data in the target area that is outside the contour line of the calculation area, and 1 represents the point data in the target area that is inside the contour line of the calculation area.

[0091] Step 104: Based on the volume calculation principle of the grid method, matrix operations are performed on the target elevation grid data to obtain the volume matrix of each geological layer.

[0092] In the application, matrix operations can be performed on the target elevation grid data based on the grid volume calculation principle; see Figure 5 , Figure 5 The following is a schematic diagram of the process of calculating the volume of each geological layer of a prism using the grid method. Figure 5 As shown, there is a one-to-one correspondence between the original terrain elevation grid data, the design surface elevation grid data and the interface elevation grid data of adjacent geological layers in the target elevation grid data, that is, the elevation grid data of the original terrain surface, the design surface and the interface of each geological layer are exactly the same in size, and the grids correspond one to one, together constituting a calculation model of the excavation body composed of a number of prisms.

[0093] When calculating the volume of a prism composed of upper and lower grids, the following formula can be used:

[0094] V r =f(upper layer, lower layer)=S×h r ;

[0095] Where V r represents the volume of the rth prism; S represents the area of ​​the grid. The smaller the grid spacing, the smaller the square area. Therefore, by reducing the grid spacing, the accuracy of the square volume calculation can be improved to a certain extent; h r Indicates the height of the rth prism. When calculating, the average height of the four sides of the prism is taken as the calculated height of the prism. Where r is a positive integer.

[0096] Take the number of geological layers as 3 as an example, see Figure 5 , at this time, the interface elevation grid data of adjacent geological layers specifically includes: the first interface elevation grid data and the second interface elevation grid data; in the application, the following formula can be used to calculate the prism volume:

[0097]

[0098] Where V r，1represents the volume of the rth prismatic geological layer 1, V r，2 represents the volume of the rth prismatic geological layer 2, V r，3 It represents the volume of the rth prismatic geological layer 3; the terrain surface, the design surface, the interface 1 and the interface 2 correspond to the original terrain elevation grid data, the design surface elevation grid data, the first interface elevation grid data and the second interface elevation grid data, respectively.

[0099] Based on the above-mentioned grid volume calculation principle, matrix operations are used to improve the calculation efficiency. Based on the terrain surface and the design surface elevation grid data, the prism volume matrix between the terrain surface and the design surface is calculated, recorded as V_matrix1; based on the interface 1 and the design surface elevation grid data, the prism volume matrix between the interface 1 and the design surface is calculated, recorded as V_matrix2; based on the interface 2 and the design surface elevation grid data, the prism volume matrix between the interface 2 and the design surface is calculated, recorded as V_matrix3. Based on these prism volume matrices, the volume matrices of each geological layer are calculated. The volume matrix of geological layer 1 is V_matrix1-V_matrix2, the volume matrix of geological layer 2 is V_matrix2-V_matrix3, and the volume matrix of geological layer 3 is V_matrix3.

[0100] Step 105: Obtain the volume excavation volume of each geological layer according to the mask matrix and the volume matrix of each geological layer.

[0101] Among them, the volume excavation volume of each geological layer is used to characterize the volume of each geological layer within the contour line of the calculation area. In the application, the volume matrix of each geological layer can be multiplied with the mask matrix respectively to obtain the multiplication result of each geological layer; the sum of the positive values ​​of each multiplication result is calculated respectively to obtain the volume excavation volume of each geological layer.

[0102] Taking the above three geological layers as an example, the mask matrix is ​​multiplied by the volume matrices of geological layer 1, geological layer 2 and geological layer 3 respectively, and then the sum of their positive values ​​is calculated respectively, so that the volume excavation volume of each geological layer within the contour line can be obtained.

[0103] At this point, the calculation of the volume of excavation of the "body" has been completed. Figure 6 , Figure 6 Schematic diagram of the calculation process for the volume of excavation of the "body".

[0104] Step 106: Use the mask matrix and the step height to perform layered processing on the target elevation grid data to obtain multiple layered elevation grid data, and perform matrix operations on each layered elevation grid data to obtain the volume matrix of each geological layer in each layer. Then, based on the mask matrix and the volume matrix of each geological layer in each layer, obtain the excavation volume of each geological layer in each layer.

[0105] See also Figure 7 , Figure 7 Schematic diagram of the calculation process for “layer” excavation volume.

[0106] The target elevation grid data is processed in layers using the mask matrix and the step height to obtain multiple layers of elevation grid data, including:

[0107] Based on the target elevation grid data and mask matrix, the highest point of the terrain and the lowest point of the design surface within the contour line of the calculation area are obtained;

[0108] Determine the layering parameters according to the step height, the highest point of the terrain and the lowest point of the design surface; the layering parameters include the number of layers of the target elevation grid data, the top elevation of each layer, and the bottom elevation of each layer;

[0109] According to the layering parameters, a correction operation is performed on the target elevation grid data to obtain a plurality of layered elevation grid data.

[0110] The correction operation includes:

[0111] Based on the number of layers, the original terrain elevation grid data and interface elevation grid data that are higher than the layer top elevation in each layer of the target elevation grid data are corrected to the layer top elevation, and the design surface elevation grid data that are lower than the layer bottom elevation in each layer of the target elevation grid data are corrected to the layer bottom elevation.

[0112] See also Figure 7 and Figure 8 After obtaining the highest point of the terrain and the lowest point of the design surface within the contour line of the calculation area, the number of layers layer_num is calculated in combination with the step height; then, the elevation grid data of the original terrain, design surface and adjacent geological layer interface of each layer in the target elevation grid data can be corrected layer by layer from top to bottom.

[0113] Further, matrix operations can be performed on the elevation grid data of each layer to obtain the volume matrix of each geological layer in each layer. When performing matrix operations on the elevation grid data of each layer, the volume calculation principle of the grid method can also be used. The calculation process can be referred to the relevant introduction of the above step 104, which will not be repeated here.

[0114] Finally, the volume matrix of each geological layer in each layer can be multiplied with the mask matrix to obtain the multiplication result of each geological layer in each layer, and then the sum of the positive values ​​of each multiplication result is calculated to obtain the excavation volume of each geological layer in each layer.

[0115] At this point, the calculation of the “layer” excavation volume has been completed.

[0116] Step 107: Use point cloud data, preset strip spacing and strip division reference lines to perform strip processing on the layered elevation grid data to obtain multiple strips of elevation grid data and corresponding multiple strip mask matrices, and perform matrix operations on each strip of elevation grid data to obtain the volume matrix of each geological layer in each strip, and then obtain the excavation volume of each geological layer in each strip based on each strip mask matrix and the corresponding volume matrix of each geological layer in each strip.

[0117] Among them, the preset bar spacing and bar division reference lines can be flexibly set according to actual needs.

[0118] See also Fig. 9 , Fig. 9 Schematic diagram of the calculation process for the volume of excavation for a "strip".

[0119] Take the hierarchical elevation grid data of the i-th layer as an example, Fig.10 As shown in FIG. 1 , the target point cloud whose original terrain of the i-th layer is higher than the design surface is selected from the point cloud data. According to the x and y range of the i-th layer point cloud, the distance from the i-th layer point cloud to the strip division reference line is calculated to obtain the coordinates of the nearest and farthest tangent points, and the strip division start line and end line are obtained respectively. The spacing between the two lines is calculated and divided by the preset strip spacing to obtain the strip division number strip_num. Similarly, the two intersection lines of the normal of the strip division reference line and the x and y range of the point cloud can be obtained. The intersection point with the two edge lines of the j-th strip determines the contour line of the j-th strip of the i-th layer. Based on this contour line, the strip interception area (such as Fig.10 The interception range shown in is used to intercept the layered elevation grid data of the i-th layer, further reduce the calculation range, and obtain the strip elevation grid data of the j-th strip. After that, based on the strip elevation grid data of the j-th strip, the volume matrix of each geological layer in the i-th layer and the j-th strip is obtained, and the corresponding strip mask matrix of the i-th layer and the j-th strip is made; finally, based on the volume matrix and the strip mask matrix of each geological layer in the i-th layer and the j-th strip, the excavation volume of each geological layer in the i-th layer and the j-th strip is obtained.

[0120] It should be noted that the method for obtaining the volume matrix of each geological layer in each strip is similar to the method for obtaining the volume matrix of each geological layer. For details, please refer to the relevant introduction in step 104; the method for making each strip mask matrix is ​​similar to the mask matrix. According to the contour lines of each strip and the corresponding strip interception area, a strip mask matrix composed of 0 and 1 is made.

[0121] The above method can be used to obtain the excavation volume of each geological layer in each strip of each layer; thus, the calculation of the excavation volume of the "strip" is completed.

[0122] See also Fig.10 The above-mentioned method of obtaining multiple strips of elevation grid data in the layered elevation grid data according to the point cloud data, the preset strip spacing and the strip division reference line specifically includes:

[0123] Based on the point cloud data, a target point cloud corresponding to a portion of the original terrain elevation grid data that is higher than the design surface elevation grid data is obtained;

[0124] Determine multiple stripe interception areas according to the target point cloud, the preset stripe spacing and the stripe division reference line;

[0125] Based on the multiple interception areas, the layered elevation grid data is intercepted to obtain the corresponding multiple strips of elevation grid data.

[0126] Furthermore, a plurality of stripe interception areas are determined according to the target point cloud, the preset stripe spacing and the stripe division reference line, including:

[0127] The data point closest to the dividing reference line in the target point cloud is taken as the first target point, and the data point farthest from the dividing reference line is taken as the second target point;

[0128] Determine a strip division start line passing through the first target point and a strip division end line passing through the second target point; the strip division start line and the strip division end line are both parallel to the strip division reference line;

[0129] Determine a normal line of a partition reference line that is outside the range of the target point cloud;

[0130] The data point closest to the normal line in the target point cloud is taken as the third target point, and the data point farthest from the normal line is taken as the fourth target point;

[0131] Determine a first boundary line passing through the third target point and a second boundary line passing through the fourth target point; the first boundary line and the second boundary line are both parallel to the normal line;

[0132] Calculate the interval distance between the strip division start line and the strip division end line;

[0133] Divide the interval distance by the preset bar spacing to obtain the number of bar divisions;

[0134] Dividing a rectangular area surrounded by a strip division start line, a strip division end line, a first boundary line and a second boundary line into a plurality of strip areas according to the number of strip divisions;

[0135] The minimum circumscribed rectangular frame of each area is determined, and the obtained multiple minimum circumscribed rectangular frames are used as the corresponding multiple intercepted areas.

[0136] Step 108: Use point cloud data, preset block spacing and block division reference lines to block the strip elevation grid data to obtain multiple block mask matrices, and obtain the excavation volume of each geological layer in each block based on the volume matrix of each geological layer in each strip and the corresponding block mask matrix.

[0137] Among them, the preset block spacing and block division reference line can be flexibly set according to actual needs.

[0138] See also Fig.11 , Fig.11 Schematic diagram of the calculation process for the excavation volume of a "block".

[0139] First, the volume matrix of each geological layer in the i-th layer and the j-th line can be obtained through the above steps; secondly, referring to the principle of calculating the contour line of the i-th layer and the j-th line in step 108, based on the preset block spacing and block division reference line, the contour line of the i-th layer, the j-th line and the k-th block is calculated, and the block mask matrix of the i-th layer, the j-th line and the k-th block is made; finally, based on the volume matrix of each geological layer in the i-th layer and the block mask matrix of the i-th layer, the j-th line and the k-th block, the excavation volume of each geological layer in the i-th layer, the j-th line and the k-th block is obtained; in this way, the excavation volume of each geological layer in each layer, each line and each block can be obtained, and the calculation of the "block" excavation volume is completed.

[0140] Corresponding to the aforementioned application function realization method embodiment, the present invention also provides a multi-scale excavation area earthwork volume rapid calculation system for distinguishing rock types and corresponding embodiments.

[0141] See also Fig.12 , Fig.12 This is a schematic diagram of the module structure of the multi-scale excavation area earthwork volume rapid calculation system for distinguishing lithology; the system includes:

[0142] The data acquisition module 131 is used to acquire the engineering data of the area to be excavated and calculate the area contour line; the engineering data includes: original terrain data, design surface data and exploration hole data;

[0143] The data processing module 132 is used to:

[0144] Acquire corresponding point cloud data according to the engineering data, and perform spatial interpolation based on the point cloud data to obtain corresponding elevation grid data; the point cloud data includes: original terrain point cloud, design surface point cloud, and thickness point cloud of each geological layer;

[0145] Based on the calculation area contour, determine the target area in the elevation grid data, intercept the elevation grid data in the target area to obtain target elevation grid data, and generate a mask matrix corresponding to the target elevation grid data;

[0146] The volume excavation calculation module 133 is used to:

[0147] Based on the grid volume calculation principle, matrix operations are performed on the target elevation grid data to obtain the volume matrix of each geological layer;

[0148] According to the mask matrix and the volume matrix of each geological layer, the volume excavation volume of each geological layer is obtained; the volume excavation volume of each geological layer is used to characterize the volume of each geological layer within the contour line of the calculation area;

[0149] The layer excavation volume calculation module 134 is used to:

[0150] Using the mask matrix and the step height, the elevation grid data is processed in layers to obtain a plurality of layered elevation grid data, and matrix operations are performed on the elevation grid data of each layer to obtain a volume matrix of each geological layer in each layer, and then the excavation volume of each geological layer in each layer is obtained according to the mask matrix and the volume matrix of each geological layer in each layer;

[0151] The excavation volume calculation module 135 is used to:

[0152] The layered elevation grid data is processed into strips using point cloud data, preset strip spacing and strip division reference lines to obtain multiple strips of elevation grid data and corresponding multiple strip mask matrices, and matrix operations are performed on each strip of elevation grid data to obtain a volume matrix of each geological layer in each strip, and then the excavation volume of each geological layer in each strip is obtained according to each strip mask matrix and the corresponding volume matrix of each geological layer in each strip;

[0153] The block excavation volume calculation module 136 is used to:

[0154] The strip elevation grid data is divided into blocks using point cloud data, preset block spacing and block division reference lines to obtain multiple block mask matrices, and the excavation volume of each geological layer in each block is obtained based on the volume matrix of each geological layer in each strip and the corresponding block mask matrix.

[0155] Regarding the system in the above embodiment, the specific manner in which each module performs operations has been described in detail in the embodiment of the method, and will not be elaborated again here.

[0156] The embodiments of the present invention have been described above, and the above description is exemplary, not exhaustive, and is not limited to the disclosed embodiments. Many modifications and variations will be apparent to those of ordinary skill in the art without departing from the scope and spirit of the described embodiments. The selection of terms used herein is intended to best explain the principles of the embodiments, practical applications, or improvements to the technology in the market, or to enable other persons of ordinary skill in the art to understand the embodiments disclosed herein.

Claims

1. A method for quickly calculating earthwork volume in a multi-scale excavation area for distinguishing lithology, characterized in that: include: Obtain engineering data of the area to be excavated and calculate the area contour; The engineering data includes: original terrain data, design surface data and exploration hole data; Acquire corresponding point cloud data according to the engineering data, and perform spatial interpolation based on the point cloud data to obtain corresponding elevation grid data; the point cloud data includes: original terrain point cloud, design surface point cloud, and thickness point cloud of each geological layer; Based on the calculation area contour line, determine the target area in the elevation grid data, intercept the elevation grid data in the target area to obtain the target elevation grid data, and generate a mask matrix corresponding to the target elevation grid data; the value in the mask matrix is ​​0 or 1; wherein 0 is used to represent the point data in the target area that is outside the calculation area contour line, and 1 is used to represent the point data in the target area that is inside the calculation area contour line; Based on the grid volume calculation principle, matrix operations are performed on the target elevation grid data to obtain the volume matrix of each geological layer; According to the mask matrix and the volume matrix of each geological layer, the volume excavation volume of each geological layer is obtained; the volume excavation volume of each geological layer is used to characterize the volume of each geological layer within the contour line of the calculation area; Using the mask matrix and the step height, the elevation grid data is processed in layers to obtain a plurality of layered elevation grid data, and matrix operations are performed on the elevation grid data of each layer to obtain a volume matrix of each geological layer in each layer, and then the excavation volume of each geological layer in each layer is obtained according to the mask matrix and the volume matrix of each geological layer in each layer; The layered elevation grid data is processed into strips using point cloud data, preset strip spacing and strip division reference lines to obtain multiple strips of elevation grid data and corresponding multiple strip mask matrices, and matrix operations are performed on each strip of elevation grid data to obtain a volume matrix of each geological layer in each strip, and then the excavation volume of each geological layer in each strip is obtained according to each strip mask matrix and the corresponding volume matrix of each geological layer in each strip; The strip elevation grid data is divided into blocks using point cloud data, preset block spacing and block division reference lines to obtain multiple block mask matrices, and the excavation volume of each geological layer in each block is obtained based on the volume matrix of each geological layer in each strip and the corresponding block mask matrix.

2. The multi-scale excavation area earthwork volume rapid calculation method for distinguishing lithology according to claim 1 is characterized in that: The elevation grid data includes: original terrain elevation grid data, design surface elevation grid data, and interface elevation grid data of adjacent geological layers; the spatial interpolation based on point cloud data to obtain corresponding elevation grid data specifically includes: Performing spatial interpolation on the original terrain point cloud to obtain corresponding original terrain elevation grid data; Performing spatial interpolation on the design surface point cloud to obtain corresponding design surface elevation grid data; Perform spatial interpolation on the thickness point cloud of each geological layer to obtain the corresponding grid data of the thickness of each geological layer; The interface elevation grid data of adjacent geological layers are obtained based on the original terrain elevation grid data and the thickness grid data of each geological layer.

3. The method for rapidly calculating earthwork volume in a multi-scale excavation area for distinguishing lithology according to claim 1 is characterized in that: The step of determining the target area in the elevation grid data based on the calculation area contour line includes: Extracting contour line point data in the contour line of the calculation area; Obtaining the maximum and minimum values ​​of the x-coordinate and the maximum and minimum values ​​of the y-coordinate in the contour line point data; Determine the minimum circumscribed rectangular frame of the calculation area contour line in the elevation grid data according to the maximum and minimum values ​​of the x-coordinate and the maximum and minimum values ​​of the y-coordinate; The four sides of the minimum circumscribed rectangular frame are expanded according to a preset distance to obtain a target area.

4. The method for rapidly calculating earthwork volume in a multi-scale excavation area for distinguishing lithology according to claim 1 is characterized in that: According to the mask matrix and the volume matrix of each geological layer, the volume excavation volume of each geological layer is obtained, including: Multiplying the volume matrix of each geological layer by the mask matrix respectively to obtain the multiplication result of each geological layer; The sum of the positive values ​​of the multiplication results is calculated respectively to obtain the volume excavation of each geological layer.

5. The method for rapidly calculating earthwork volume in a multi-scale excavation area for distinguishing lithology according to claim 1 is characterized in that: The target elevation grid data is processed in layers using the mask matrix and the step height to obtain a plurality of layered elevation grid data, including: Based on the target elevation grid data and the mask matrix, obtaining the highest point of the terrain and the lowest point of the design surface within the contour line of the calculation area; Determine the layering parameters according to the step height, the highest point of the terrain and the lowest point of the design surface; the layering parameters include the number of layers of the target elevation grid data, the top elevation of each layer, and the bottom elevation of each layer; According to the layered parameters, a correction operation is performed on the target elevation grid data to obtain a plurality of layered elevation grid data; The performing of the correction operation comprises: Based on the number of layers, the original terrain elevation grid data and interface elevation grid data in each layer of the target elevation grid data that are higher than the layer top elevation are corrected to the layer top elevation, and the design surface elevation grid data in each layer of the target elevation grid data that are lower than the layer bottom elevation are corrected to the layer bottom elevation.

6. The method for rapidly calculating earthwork volume in a multi-scale excavation area for distinguishing lithology according to claim 5 is characterized in that: The layered elevation grid data is processed into strips using the point cloud data, the preset strip spacing and the strip division reference line to obtain a plurality of strip elevation grid data, including: Based on the point cloud data, a target point cloud corresponding to a portion of the layered elevation grid data where the original terrain elevation grid data is higher than the design surface elevation grid data is obtained; Determining a plurality of stripe interception areas according to the target point cloud, a preset stripe spacing and a stripe division reference line; Based on the plurality of intercepted areas, the layered elevation grid data is intercepted to obtain a corresponding plurality of strips of elevation grid data.

7. The method for rapidly calculating earthwork volume in a multi-scale excavation area for distinguishing lithology according to claim 6 is characterized in that: Determining a plurality of stripe interception areas according to the target point cloud, the preset stripe spacing and the stripe division reference line includes: The data point closest to the dividing reference line in the target point cloud is taken as the first target point, and the data point farthest from the dividing reference line is taken as the second target point; Determine a strip division start line passing through the first target point, and a strip division end line passing through the second target point; the strip division start line and the strip division end line are both parallel to the strip division reference line; Determine a normal line of the partition reference line that is outside the range of the target point cloud; The data point closest to the normal line in the target point cloud is used as the third target point, and the data point farthest from the normal line is used as the fourth target point; Determine a first boundary line passing through the third target point, and a second boundary line passing through the fourth target point; the first boundary line and the second boundary line are both parallel to the normal line; Calculating the interval distance between the strip division starting line and the strip division ending line; Dividing the interval distance by the preset strip spacing to obtain the number of strip divisions; Dividing a rectangular area enclosed by the strip division start line, the strip division end line, the first boundary line and the second boundary line into a plurality of strip areas according to the strip division number; The minimum circumscribed rectangular frame of each of the strip regions is determined, and the obtained multiple minimum circumscribed rectangular frames are used as corresponding multiple strip interception regions.

8. A multi-scale excavation area earthwork volume rapid calculation system for distinguishing lithology, characterized in that: include: A data acquisition module is used to obtain engineering data of the area to be excavated and calculate the area contour line; The engineering data includes: original terrain data, design surface data and exploration hole data; Data processing module for: Acquire corresponding point cloud data according to the engineering data, and perform spatial interpolation based on the point cloud data to obtain corresponding elevation grid data; the point cloud data includes: original terrain point cloud, design surface point cloud, and thickness point cloud of each geological layer; Based on the calculation area contour line, determine the target area in the elevation grid data, intercept the elevation grid data in the target area to obtain the target elevation grid data, and generate a mask matrix corresponding to the target elevation grid data; the value in the mask matrix is ​​0 or 1; wherein 0 is used to represent the point data in the target area that is outside the calculation area contour line, and 1 is used to represent the point data in the target area that is inside the calculation area contour line. Volume excavation calculation module, used for: Based on the grid volume calculation principle, matrix operations are performed on the target elevation grid data to obtain the volume matrix of each geological layer; According to the mask matrix and the volume matrix of each geological layer, the volume excavation volume of each geological layer is obtained; the volume excavation volume of each geological layer is used to characterize the volume of each geological layer within the contour line of the calculation area; Layer excavation volume calculation module is used for: Using the mask matrix and the step height, the elevation grid data is processed in layers to obtain a plurality of layered elevation grid data, and matrix operations are performed on the elevation grid data of each layer to obtain a volume matrix of each geological layer in each layer, and then the excavation volume of each geological layer in each layer is obtained according to the mask matrix and the volume matrix of each geological layer in each layer; The excavation volume calculation module is used to: The layered elevation grid data is processed into strips using point cloud data, preset strip spacing and strip division reference lines to obtain multiple strips of elevation grid data and corresponding multiple strip mask matrices, and matrix operations are performed on each strip of elevation grid data to obtain a volume matrix of each geological layer in each strip, and then the excavation volume of each geological layer in each strip is obtained according to each strip mask matrix and the corresponding volume matrix of each geological layer in each strip; Block excavation volume calculation module, used for: The strip elevation grid data is divided into blocks using point cloud data, preset block spacing and block division reference lines to obtain multiple block mask matrices, and the excavation volume of each geological layer in each block is obtained based on the volume matrix of each geological layer in each strip and the corresponding block mask matrix.

Images