CN119004603B - Rapid calculation method and system for soil and stone quantity of multi-scale excavation area for distinguishing lithology - Google Patents
Rapid calculation method and system for soil and stone quantity of multi-scale excavation area for distinguishing lithology Download PDFInfo
- Publication number
- CN119004603B CN119004603B CN202411075554.8A CN202411075554A CN119004603B CN 119004603 B CN119004603 B CN 119004603B CN 202411075554 A CN202411075554 A CN 202411075554A CN 119004603 B CN119004603 B CN 119004603B
- Authority
- CN
- China
- Prior art keywords
- data
- volume
- strip
- grid data
- layer
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/13—Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/05—Geographic models
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/20—Finite element generation, e.g. wire-frame surface description, tesselation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- Software Systems (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Physics (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Computer Graphics (AREA)
- Computer Hardware Design (AREA)
- Data Mining & Analysis (AREA)
- General Engineering & Computer Science (AREA)
- Databases & Information Systems (AREA)
- Computing Systems (AREA)
- Remote Sensing (AREA)
- Architecture (AREA)
- Civil Engineering (AREA)
- Algebra (AREA)
- Structural Engineering (AREA)
- Evolutionary Computation (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention discloses a rapid calculation method and a rapid calculation system for soil and stone quantities of a multi-scale excavation area for distinguishing lithology, which relate to the technical field of soil and stone quantity calculation and comprise the steps of obtaining engineering data of an area to be excavated, selecting a contour line of the calculation area, carrying out spatial interpolation on point cloud data obtained based on the engineering data to obtain elevation grid data, intercepting target height Cheng Wangge data in a target area based on the contour line of the calculation area, generating a corresponding mask matrix, calculating and obtaining volume excavation quantities of each geological layer, obtaining the excavation quantities of each geological layer in each layering according to the mask matrix and the height of a ladder section, obtaining the excavation quantities of each geological layer in each layering according to the point cloud data, preset strip spacing and strip division reference lines, obtaining the excavation quantities of each geological layer in each strip according to the point cloud data, the preset block spacing and the strip division reference lines, and obtaining the excavation quantities of each geological layer in each block.
Description
Technical Field
The invention relates to the technical field of earth and stone volume calculation, in particular to a method and a system for rapidly calculating the earth and stone volume of a multiscale excavation area for distinguishing lithology.
Background
The calculation of the earth and rock mass is of great significance to the construction of large civil engineering such as airports, rock-fill dams and the like, the calculation is directly related to scheme optimization and engineering cost budget, and different filling areas have different requirements on the material property of the excavated materials, so that the excavation engineering mass of each geological layer is necessary to be calculated. The current common method is to create a terrain surface and each geologic layer surface in Civi D software, construct a three-dimensional geologic model using the extract entity commands from the surfaces, and query the volume of each geologic layer.
In the process of mixing the earthwork, in order to better meet the requirements of progress and material matching, the excavation engineering quantity needs to be analyzed with finer granularity, namely, the excavation engineering quantity of each geological layer in each layer, each strip and each block is counted, however, when the fine granularity analysis of the excavation engineering quantity is carried out by using Civi D software at present, a plurality of sub-sectioning operations are needed to be carried out on the three-dimensional geological model, and the calculation efficiency of the mode is low.
In the prior art, although a related technical scheme of calculating the earth and stone volume by establishing a curved surface calculation volume by utilizing Civi D or calculating the volume of each geological layer by establishing a geological layering model exists, the earth and stone volume calculation of each geological layer in different scales of body-layer-strip-block is not further realized, for example, the inventions with application numbers of CN202010059628.4, CN202110009012.0, CN201811440239.5 and the like are disclosed.
Disclosure of Invention
In order to solve the problems in the prior art, the invention provides a rapid calculation method for the soil and stone quantity of a multi-scale excavation area for distinguishing lithology, which can improve the calculation efficiency, rapidly obtain the soil and stone quantity of each geological layer in different scales of body, layer, strip and block, and meet the construction period requirement of large-scale engineering.
In order to achieve the above purpose, the invention provides a rapid calculation method for the earth and stone quantity of a multiscale excavation region for distinguishing lithology, which comprises the following steps:
Acquiring engineering data of an area to be excavated and calculating an area contour line, wherein the engineering data comprises original terrain data, design surface data and exploration hole data;
acquiring corresponding point cloud data according to the engineering data, and performing spatial interpolation based on the point cloud data to acquire corresponding elevation grid data, wherein the point cloud data comprises original terrain point cloud, design surface point cloud and thickness point cloud of each geological layer;
determining a target area in the elevation grid data based on the calculated area contour line, intercepting the elevation grid data in the target area to obtain target Cheng Wangge data, and generating a mask matrix corresponding to the target elevation grid data;
based on a square grid method volume calculation principle, performing matrix operation on the target elevation grid data to obtain a volume matrix of each geological layer;
obtaining the volume of each geological layer according to the mask matrix and the volume matrix of each geological layer, wherein the volume of each geological layer is used for representing the volume of each geological layer in the outline of the calculation region;
Carrying out layering treatment on the elevation grid data by utilizing the mask matrix and the bench height to obtain a plurality of layering height Cheng Wangge data, carrying out matrix operation on each layering height Cheng Wangge data to obtain a volume matrix of each geological layer in each layering, and then obtaining the excavation square quantity of each geological layer in each layering according to the mask matrix and the volume matrix of each geological layer in each layering;
Performing striping processing on the layered high Cheng Wangge data by utilizing point cloud data, preset stripe intervals and stripe division reference lines to obtain a plurality of stripe high Cheng Wangge data and a plurality of corresponding stripe mask matrixes, performing matrix operation on each stripe high Cheng Wangge data to obtain a volume matrix of each geological layer in each stripe, and then obtaining the excavation amount of each geological layer in each stripe according to each stripe mask matrix and the corresponding volume matrix of each geological layer in each stripe;
And carrying out block division processing on the strip elevation grid data by utilizing the point cloud data, the preset block spacing and the block division reference line to obtain a plurality of block mask matrixes, and obtaining the excavation amount of each geological layer in each block according to the volume matrixes of each geological layer in each strip and the corresponding block mask matrixes.
Optionally, the elevation grid data comprises original terrain elevation grid data, design surface height Cheng Wangge data and interface elevation grid data of adjacent geological layers, and the method comprises the steps of performing spatial interpolation based on point cloud data to obtain corresponding elevation grid data, and specifically comprises the following steps:
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 height Cheng Wangge data;
performing spatial interpolation on thickness point clouds of each geological layer to obtain corresponding geological layer thickness grid data;
And acquiring interface elevation grid data of adjacent geological layers according to the original terrain elevation grid data and the geological layer thickness grid data.
Optionally, the determining, based on the calculated region contour line, a target region in the elevation grid data includes:
extracting contour line point data in the contour line of the calculation region;
Obtaining the maximum value and the minimum value of the x coordinate and the maximum value and the minimum value of the y coordinate in the profile line point data;
determining a minimum circumscribed rectangular frame of the calculated area contour line in the elevation grid data according to the maximum value and the minimum value of the x coordinate and the maximum value and the minimum value of the y coordinate;
and expanding four sides of the minimum circumscribed rectangular frame according to a preset distance to obtain a target area.
Optionally, the value in the mask matrix is 0 or 1, wherein 0 is used for representing point data in the target area, which is located outside the outline of the calculation area, and 1 is used for representing point data in the target area, which is located inside the outline of the calculation area.
Optionally, obtaining the volume excavation square quantity of each geological layer according to the mask matrix and the volume matrix of each geological layer includes:
multiplying the volume matrix of each geological layer with the mask matrix to obtain multiplication results of each geological layer;
And respectively obtaining the sum of positive values of the multiplication results to obtain the volume excavation amount of each geological layer.
Optionally, using the mask matrix and the bench height, performing hierarchical processing on the target elevation grid data to obtain a plurality of hierarchical height Cheng Wangge data, including:
Acquiring a terrain highest point and a design surface lowest point in the contour line of the calculation region based on the target elevation grid data and the mask matrix;
Determining layering parameters according to the bench height, the highest point of the terrain and the lowest point of the design surface, wherein the layering parameters comprise layering quantity of the target elevation grid data, layering top elevation of each layering and layering bottom elevation of each layering;
performing correction operation on the target elevation grid data according to the layering parameters to obtain a plurality of layering height Cheng Wangge data;
the performing a corrective action includes:
And correcting original terrain elevation grid data and interface elevation grid data which are higher than the layer top elevation in each layer of the target elevation grid data into layer top elevation based on the number of layers, and correcting design surface elevation Cheng Wangge data which are lower than the layer bottom elevation in each layer of the target elevation grid data into layer bottom elevation.
Optionally, performing striping processing on the layered height Cheng Wangge data by using the point cloud data, a preset stripe pitch and a stripe division reference line to obtain a plurality of stripe height Cheng Wangge data, including:
acquiring target point clouds corresponding to parts of the original terrain elevation grid data higher than the design surface height Cheng Wangge data in the layered elevation grid data based on the point cloud data;
Determining a plurality of strip interception areas according to the target point cloud, the preset strip intervals and the strip division reference lines;
And intercepting the hierarchical elevation grid data based on the plurality of strip intercepting areas to obtain a plurality of corresponding strip height Cheng Wangge data.
Optionally, determining a plurality of stripe intercepting areas according to the target point cloud, a preset stripe interval and a stripe dividing reference line includes:
Taking a data point closest to the dividing reference line in the target point cloud as a first target point, and a data point farthest from the dividing reference line as a second target point;
Determining a strip division starting line passing through the first target point and a strip division ending line passing through the second target point, wherein the strip division starting line and the strip division ending line are parallel to the strip division reference line;
determining a normal line of the dividing reference line which is positioned outside the target point cloud range;
Taking a data point closest to the normal line in the target point cloud as a third target point, and a data point farthest from the normal line as a fourth target point;
determining a first boundary line passing through the third target point and a second boundary line passing through the fourth target point, wherein the first boundary line and the second boundary line are parallel to the normal;
calculating the interval distance between the strip division starting line and the strip division ending line;
dividing the interval distance by the preset strip interval to obtain the strip dividing number;
Dividing a rectangular area surrounded by the strip division starting line, the strip division ending line, the first boundary line and the second boundary line into a plurality of strip areas according to the strip division quantity;
and determining the minimum circumscribed positive rectangular frame of each strip area, and taking the obtained multiple minimum circumscribed positive rectangular frames as a corresponding multiple strip interception areas.
According to the specific embodiment provided by the invention, the invention discloses the following technical effects:
The rapid calculation method for the soil and stone amount of the multi-scale excavation area for distinguishing lithology provided by the invention can rapidly obtain the soil and stone amount of each geological layer (different lithology) in different scales of a body, a layer, a strip and a block, has the advantages of less time and high efficiency compared with the prior art that the process of carrying out the segmentation analysis by using Civi D software, particularly greatly reduces the time spent in carrying out the excavation engineering amount analysis of the fine granularity of the layer, the strip and the block, and meanwhile, the calculation result obtained by the calculation method provided by the invention has smaller phase difference compared with the calculation result obtained by Civi D, can meet the requirement of engineering on precision, and can be directly applied to practical engineering.
Drawings
The foregoing and other objects, features and advantages of the invention will be apparent from the following more particular descriptions of exemplary embodiments of the invention as illustrated in the accompanying drawings wherein like reference numbers generally represent like parts throughout the exemplary embodiments of the invention.
Fig. 1 is a flow chart of a method for rapidly calculating the earth and rock mass of a multi-scale excavated area for distinguishing lithology according to an embodiment of the present invention;
FIG. 2 is a schematic view of three-dimensional and XOY projections of point clouds before and after Kriging interpolation according to an embodiment of the present invention;
FIG. 3 is a schematic diagram of an elevation grid data preservation matrix format according to an embodiment of the present invention;
FIG. 4 is a schematic diagram of capturing elevation grid data of a target area according to an embodiment of the present invention;
FIG. 5 is a schematic diagram illustrating a process for calculating the volume of each geological formation of a prism using a square method according to an embodiment of the present invention;
FIG. 6 is a schematic diagram of a calculation flow of the "body" excavation amount shown in the embodiment of the present invention;
FIG. 7 is a schematic diagram of a calculation flow of "layer" excavation amount according to an embodiment of the present invention;
FIG. 8 is a schematic diagram of an i-th layer elevation grid data correction according to an embodiment of the present invention;
FIG. 9 is a schematic diagram of a calculation flow of the "strip" excavation amount shown in the embodiment of the present invention;
FIG. 10 is a schematic diagram illustrating determination of a data interception range of an ith layer and a jth elevation grid according to an embodiment of the present invention;
FIG. 11 is a schematic diagram of a calculation flow of "block" excavation amount shown in an embodiment of the present invention;
fig. 12 is a schematic block diagram of a rapid calculation system for the earth and rock mass in a multi-scale excavated area for distinguishing lithology according to an embodiment of the present invention.
Detailed Description
The following description of the embodiments of the present invention will be made clearly and completely with reference to the accompanying drawings, in which it is apparent that the embodiments described are only some embodiments of the present invention, but not all embodiments. All other embodiments, which can be made by those skilled in the art based on the embodiments of the invention without making any inventive effort, are intended to be within the scope of the invention.
The invention aims to provide a rapid calculation method for the earthwork quantity of a multi-scale excavated area for distinguishing lithology, which can rapidly calculate the earthwork quantity of each geological layer in different scales of a body, a layer, a strip and a block.
Referring to fig. 1, fig. 1 is a flow chart of a method for rapidly calculating the earth and rock mass of a multi-scale excavated area for distinguishing lithology.
The method for calculating the earthwork quantity comprises the following steps:
and 101, acquiring engineering data of an area to be excavated and calculating an area contour line.
The engineering data comprises original terrain data, design surface data and exploration hole data. The original terrain data can be obtained by measuring in the field by using a mapping tool, for example, the original terrain data can be obtained by scanning an area to be excavated through a high-precision laser scanner or a three-dimensional camera carried by an unmanned aerial vehicle, and the like, the design surface data can be obtained according to actual engineering requirements, the design surface in the calculation of the earth and stone volume refers to a calculation reference surface of the earth and stone volume determined according to engineering requirements in the design stage, and represents the expected ground shape or structure in engineering design, and the exploration hole data of the area to be excavated can be obtained through geological exploration means.
The original topography data and the design surface data have the same storage structure, and each point data is composed of an x-coordinate, a y-coordinate, and an elevation value z. The survey hole data consists of a survey point number, an x coordinate, a y coordinate, and a plurality of geological layer thicknesses.
Step 102, obtaining corresponding point cloud data according to engineering data, and performing spatial interpolation based on the point cloud data to obtain corresponding elevation grid data.
The point cloud data comprise original terrain point cloud, design surface point cloud and thickness point cloud of each geological layer. In the application, corresponding 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.
After the point cloud data is obtained, spatial interpolation can be performed on the point cloud data by adopting methods such as inverse distance weight Interpolation (IDW), kriging, natural neighborhood interpolation (Natural Neighbor), radial basis function interpolation (RBF) and the like, wherein the Kriging interpolation is taken as an example, and fig. 2 is a schematic diagram of three-dimensional and XOY projection of the point cloud before and after the Kriging interpolation.
In application, the point cloud data can be spatially interpolated according to a preset grid interval, wherein the preset grid interval can be flexibly set according to practical situations, and fig. 4 shows elevation grid data obtained by interpolation according to a 25m grid interval (preset grid interval) in fig. 4.
The elevation grid data comprises original terrain elevation grid data, design surface height Cheng Wangge data and interface elevation grid data of adjacent geological layers, wherein the elevation grid data corresponds to the point cloud data, and concretely, the elevation grid data is obtained by carrying out spatial interpolation based on the point cloud data, and the method specifically comprises the following steps:
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 height Cheng Wangge data;
performing spatial interpolation on thickness point clouds of each geological layer to obtain corresponding geological layer thickness grid data;
And acquiring interface elevation grid data of adjacent geological layers according to the original terrain elevation grid data and the geological layer thickness grid data.
The elevation grid data can be stored as a matrix format after being obtained so as to facilitate subsequent matrix operation, as shown in fig. 3, fig. 3 is a schematic diagram of a high Cheng Wangge data storage matrix format.
Step 103, determining a target area in the elevation grid data based on the calculated area contour line, intercepting the elevation grid data in the target area to obtain target Cheng Wangge data, and generating a mask matrix corresponding to the target Cheng Wangge data.
The calculation region contour line is a boundary line for determining and describing the region to be excavated and is closed curve data, and can be flexibly set according to actual engineering requirements. In the application, the x coordinate and the y coordinate of the contour line point data of the calculation region can be extracted, and the target region is determined according to the size range of the x coordinate and the y coordinate.
Illustratively, determining the target region in the elevation grid data based on the calculated region contour line may include the steps of:
Extracting contour line point data in a contour line of a calculation region;
obtaining the maximum value and the minimum value of an x coordinate and the maximum value and the minimum value of a y coordinate in the profile line point data;
determining a minimum circumscribed rectangular frame of the contour line of the calculated area in the elevation grid data according to the maximum value and the minimum value of the x coordinate and the maximum value and the minimum value of the y coordinate;
and expanding four sides of the minimum circumscribed rectangular frame according to the preset distance to obtain a target area.
The preset distance can be flexibly adjusted according to practical conditions, for example, the distance can be expanded by 30m.
Further referring to fig. 4, fig. 4 is a schematic diagram of capturing elevation grid data of a target area, and in the application of the capturing range in fig. 4, namely, the target area, the calculation range is reduced and the calculation efficiency is improved by determining the target area and capturing the elevation grid data in the target area.
After intercepting and obtaining the target Cheng Wangge data, a corresponding mask matrix can be manufactured according to the calculated area contour line, and specifically, a mask matrix composed of 0 and 1 can be generated, wherein 0 represents point data positioned outside the calculated area contour line in the target area, and 1 represents point data positioned inside the calculated area contour line in the target area.
And 104, performing matrix operation on the target elevation grid data based on the square grid volumetric calculation principle to obtain a volumetric matrix of each geological layer.
In application, the target elevation grid data can be subjected to matrix operation based on the principle of square grid method volume calculation, and referring to fig. 5, fig. 5 is a schematic diagram of a process for calculating the volumes of geological layers of prisms by using a square grid method. As shown in fig. 5, the original terrain elevation grid data, the design surface height Cheng Wangge data and the elevation grid data of the interfaces of adjacent geological layers in the target elevation grid data are in one-to-one correspondence, that is, the original terrain surface, the design surface and the elevation grid data of the interfaces of each geological layer are completely the same in size, and the grids are in one-to-one correspondence, so that an excavation body calculation model consisting of a plurality of prisms is formed together.
In calculating the prism volume composed of the upper and lower grids, the following formula can be used:
V r =f (upper layer, lower layer) =s×h r;
Wherein V r represents the volume of the r prism, S represents the area of the square lattice, the smaller the grid spacing is, the smaller the square lattice area is, therefore, the accuracy of square calculation can be improved to a certain extent by reducing the grid spacing, h r represents the height of the r prism, and the average height of four sides of the prism is taken as the calculated height of the prism during calculation. Wherein r is a positive integer.
Taking the geological layer number of 3 as an example, referring to fig. 5, the interface elevation grid data of the adjacent geological layers specifically includes first interface height Cheng Wangge data and second interface elevation grid data, and in application, the following formula can be adopted to calculate the prism volume:
Wherein V r,1 represents the volume of the r prism geological layer 1, V r,2 represents the volume of the r prism geological layer 2, V r,3 represents the volume of the r prism geological layer 3, and the terrain surface, the design surface, the interface 1 and the interface 2 correspond to the original terrain elevation grid data, the design surface height Cheng Wangge data, the first interface height Cheng Wangge data and the second sub-interface elevation grid data respectively.
Based on the square grid method volume calculation principle, matrix operation is adopted to improve calculation efficiency, a prism volume matrix between a terrain surface and a design surface is calculated based on data of the terrain surface and the design surface, and is marked as V_matrix1, a prism volume matrix between the interface 1 and the design surface is calculated based on data of the interface 1 and the design surface, and is marked as V_matrix2, and a prism volume matrix between the interface 2 and the design surface is calculated based on data of the interface 2 and the design surface, and is marked as V_matrix3. Based on the prism volume matrices, the volume matrix of each geological layer is 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.
And 105, obtaining the volume excavation square quantity of each geological layer according to the mask matrix and the volume matrix of each geological layer.
The volume excavation amount of each geological layer is used for representing the volume of each geological layer in the contour line of the calculated area. In the application, the volume matrix of each geological layer is multiplied by the mask matrix to obtain the multiplication result of each geological layer, and the sum of the positive values of the multiplication results is calculated to obtain the volume excavation square quantity of each geological layer.
Taking the above 3 geologic layers as an example, the mask matrix is multiplied by the volume matrices of geologic layer 1, geologic layer 2 and geologic layer 3 respectively, and the positive values are summed up respectively, so that the volume excavation amount of each geologic layer in the contour line can be obtained.
Thus, the calculation of the "body" excavation amount is completed, and fig. 6 is a schematic diagram of the calculation flow of the "body" excavation amount.
And 106, layering the target height Cheng Wangge data by using the mask matrix and the bench height to obtain a plurality of layered height Cheng Wangge data, performing matrix operation on each layered height Cheng Wangge data to obtain a volume matrix of each geological layer in each layered, and then obtaining the excavation amount of each geological layer in each layered according to the mask matrix and the volume matrix of each geological layer in each layered.
Referring to fig. 7, fig. 7 is a schematic diagram of a calculation flow of the "layer" excavation amount.
The layering processing is performed on the target high Cheng Wangge data by using the mask matrix and the bench height to obtain a plurality of layered high Cheng Wangge data, including:
Acquiring a highest point of the terrain and a lowest point of the design surface in the contour line of the calculated area based on the target Cheng Wangge data and the mask matrix;
Determining layering parameters according to the height of the bench, the highest point of the terrain and the lowest point of the design surface, wherein the layering parameters comprise layering quantity of target Cheng Wangge data, layering top elevation of each layering and layering bottom elevation of each layering;
And performing correction operation on the target Cheng Wangge data according to the layering parameters to obtain a plurality of layering high Cheng Wangge data.
Wherein performing the corrective action includes:
And correcting the original terrain elevation grid data and the interface elevation grid data which are higher than the layer top elevation in each layer of the target elevation grid data into the layer top elevation based on the number of layers, and correcting the design surface elevation Cheng Wangge data which are lower than the layer bottom elevation in each layer of the target elevation grid data into the layer bottom elevation.
Referring to fig. 7 and 8, after the highest point of the topography and the lowest point of the design surface in the contour line of the calculation region are obtained, the layering number layer_num is calculated by combining the bench height, and then the original topography of each layering in the target elevation grid data, the elevation grid data of the design surface and the adjacent geological layer interface can be corrected layer by layer from top to bottom.
Further, matrix operations may be performed on the data of each of the strata height Cheng Wangge to obtain a volumetric matrix for each of the geological layers in each of the strata. When matrix operation is performed on each hierarchical data Cheng Wangge, the calculation process can be referred to the related description of step 104 based on the principle of square grid method volume calculation, which is not described herein.
Finally, the volume matrix of each geological layer in each layering can be multiplied by the mask matrix to obtain the multiplication results of each geological layer in each layering, and then the sum of positive values of each multiplication result is calculated to obtain the excavation amount of each geological layer in each layering.
Thus, the calculation of the "layer" excavation amount is completed.
Step 107, performing striping processing on the layered elevation grid data by utilizing the point cloud data, the preset stripe intervals and the stripe division reference lines to obtain a plurality of stripe height Cheng Wangge data and a plurality of corresponding stripe mask matrixes, performing matrix operation on each stripe height Cheng Wangge data to obtain a volume matrix of each geological layer in each stripe, and then obtaining the excavation amount of each geological layer in each stripe according to each stripe mask matrix and the corresponding volume matrix of each geological layer in each stripe.
The preset strip spacing and the strip dividing reference line can be flexibly set according to actual requirements.
Referring to fig. 9, fig. 9 is a schematic diagram of a calculation flow of the "strip" excavation amount.
Taking the layered elevation grid data of the ith layer as an example, as shown in fig. 10, screening out the target point cloud of which the original topography of the ith layer is higher than the design surface from the point cloud data, calculating the distance from the point cloud of the ith layer to the strip division reference line according to the x and y ranges of the point cloud of the ith layer to obtain the nearest and farthest tangent point coordinates, respectively obtaining strip division starting lines and termination lines, calculating the distance between the two lines and dividing the distance by a preset strip distance to obtain strip division number strip_num, and similarly obtaining two intersecting lines of the normal line of the strip division reference line and the x and y ranges of the point cloud, determining the contour line of the jth strip of the ith layer by the intersecting point of the normal line of the strip division reference line and the point cloud x and y ranges, determining the strip interception area (the interception range shown in fig. 10) based on the contour line so as to intercept the layered elevation Cheng Wangge data of the ith layer, and further reducing the calculation range to obtain the strip elevation grid data of the jth strip. Then, based on the stripe elevation grid data of the jth stripe, acquiring a volume matrix of each geological layer in the jth stripe of the ith layer, and manufacturing a corresponding stripe mask matrix of the jth stripe of the ith layer; and finally, obtaining the excavation amount of each geological layer in the j th layer based on the volume matrix and the bar mask matrix of each geological layer in the j th layer.
It should be noted that, the method for obtaining the volume matrix of each geological layer in each stripe is similar to the method for obtaining the volume matrix of each geological layer, and can be specifically described with reference to step 104, and the method for producing each stripe mask matrix is similar to the mask matrix, and the stripe mask matrix consisting of 0 and 1 is produced according to the contour line of each stripe and the corresponding stripe interception area.
The method can obtain the excavation amount of each geological layer in each layered strip, and thus the calculation of the excavation amount of the strip is completed.
Referring to fig. 10, the obtaining a plurality of stripe height Cheng Wangge data in the layered height Cheng Wangge data according to the point cloud data, the preset stripe pitch and the stripe division reference line specifically includes:
Acquiring target point clouds corresponding to parts of original terrain elevation grid data higher than design surface height Cheng Wangge data in layered height Cheng Wangge data based on the point cloud data;
Determining a plurality of strip interception areas according to the target point cloud, the preset strip intervals and the strip division reference lines;
and based on the plurality of strip interception areas, intercepting the hierarchical height Cheng Wangge data to obtain a corresponding plurality of strip height Cheng Wangge data.
Further, determining a plurality of stripe intercepting areas according to the target point cloud, the preset stripe interval and the stripe dividing reference line comprises:
taking the data point closest to the dividing reference line in the target point cloud as a first target point, and the data point farthest from the dividing reference line as a second target point;
Determining a strip division starting line passing through a first target point and a strip division ending line passing through a second target point, wherein the strip division starting line and the strip division ending line are parallel to a strip division reference line;
determining a normal line of a dividing reference line which is positioned outside the cloud range of the target point;
Taking the data point closest to the normal line in the target point cloud as a third target point and the data point farthest from the normal line as a fourth target point;
Determining a first boundary line passing through the third target point and a second boundary line passing through the fourth target point, wherein the first boundary line and the second boundary line are parallel to the normal;
calculating the interval distance between the strip division starting line and the strip division ending line;
Dividing the interval distance by a preset strip interval to obtain the strip dividing number;
Dividing a rectangular area surrounded by a strip division starting line, a strip division ending line, a first boundary line and a second boundary line into a plurality of strip areas according to the strip division quantity;
And determining the minimum circumscribed positive rectangular frame of each strip area, and taking the obtained multiple minimum circumscribed positive rectangular frames as a corresponding multiple strip interception areas.
And 108, performing block division processing on the stripe height Cheng Wangge data by utilizing the point cloud data, the preset block spacing and the block division reference line to obtain a plurality of block mask matrixes, and obtaining the excavation amount of each geological layer in each block according to the volume matrixes of each geological layer in each stripe and the corresponding block mask matrixes.
The preset block distance and the block division reference line can be flexibly set according to actual requirements.
Referring to fig. 11, fig. 11 is a schematic diagram of a calculation flow of the "block" excavation amount.
Firstly, the volume matrix of each geological layer in the ith layer can be obtained through the steps, secondly, the principle of calculating the contour line of the jth layer in the ith layer in the step 108 is referred to, the contour line of the jth block in the ith layer is calculated based on the preset block spacing and block dividing reference line, the block mask matrix of the jth block in the ith layer is manufactured, finally, the excavation amount of each geological layer in the jth block in the ith layer is obtained based on the volume matrix of each geological layer in the jth layer and the block mask matrix of the jth block in the ith layer, and the excavation amount of each geological layer in each block in each layer can be obtained by utilizing the method, so that the calculation of the 'block' excavation amount is completed.
Corresponding to the embodiment of the application function implementation method, the invention also provides a rapid computing system for the soil and stone quantity of the multiscale excavation area for distinguishing lithology and a corresponding embodiment.
Referring to fig. 12, fig. 12 is a schematic block diagram of the rapid computing system for the earthwork volume of a multiscale excavated area for distinguishing lithology, where the system includes:
The data acquisition module 131 is used for acquiring engineering data of an area to be excavated and calculating an area contour line, wherein the engineering data comprises original topography data, design surface data and exploration hole data;
A data processing module 132 for:
acquiring corresponding point cloud data according to the engineering data, and performing spatial interpolation based on the point cloud data to acquire corresponding elevation grid data, wherein the point cloud data comprises original terrain point cloud, design surface point cloud and thickness point cloud of each geological layer;
determining a target area in the elevation grid data based on the calculated area contour line, intercepting the elevation grid data in the target area to obtain target Cheng Wangge data, and generating a mask matrix corresponding to the target elevation grid data;
the body excavation square calculation module 133 is configured to:
based on a square grid method volume calculation principle, performing matrix operation on the target elevation grid data to obtain a volume matrix of each geological layer;
obtaining the volume of each geological layer according to the mask matrix and the volume matrix of each geological layer, wherein the volume of each geological layer is used for representing the volume of each geological layer in the outline of the calculation region;
Layer excavation amount calculation module 134 for:
Carrying out layering treatment on the elevation grid data by utilizing the mask matrix and the bench height to obtain a plurality of layering height Cheng Wangge data, carrying out matrix operation on each layering height Cheng Wangge data to obtain a volume matrix of each geological layer in each layering, and then obtaining the excavation square quantity of each geological layer in each layering according to the mask matrix and the volume matrix of each geological layer in each layering;
a strip excavation amount calculation module 135 for:
Performing striping processing on the layered high Cheng Wangge data by utilizing point cloud data, preset stripe intervals and stripe division reference lines to obtain a plurality of stripe high Cheng Wangge data and a plurality of corresponding stripe mask matrixes, performing matrix operation on each stripe high Cheng Wangge data to obtain a volume matrix of each geological layer in each stripe, and then obtaining the excavation amount of each geological layer in each stripe according to each stripe mask matrix and the corresponding volume matrix of each geological layer in each stripe;
a block excavation amount calculation module 136 for:
And carrying out block division processing on the strip elevation grid data by utilizing the point cloud data, the preset block spacing and the block division reference line to obtain a plurality of block mask matrixes, and obtaining the excavation amount of each geological layer in each block according to the volume matrixes of each geological layer in each strip and the corresponding block mask matrixes.
The specific manner in which the various modules perform the operations in relation to the systems of the above embodiments have been described in detail in relation to the embodiments of the method and will not be described in detail herein.
The foregoing description of embodiments of the invention has been presented for purposes of illustration and description, and is not intended to be exhaustive or limited to the embodiments disclosed. Many modifications and variations will be apparent to those of ordinary skill in the art without departing from the scope and spirit of the various embodiments described. The terminology used herein was chosen in order to best explain the principles of the embodiments, the practical application, or the improvement of technology in the marketplace, or to enable others of ordinary skill in the art to understand the embodiments disclosed herein.
Claims (8)
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202411075554.8A CN119004603B (en) | 2024-08-06 | 2024-08-06 | Rapid calculation method and system for soil and stone quantity of multi-scale excavation area for distinguishing lithology |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202411075554.8A CN119004603B (en) | 2024-08-06 | 2024-08-06 | Rapid calculation method and system for soil and stone quantity of multi-scale excavation area for distinguishing lithology |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN119004603A CN119004603A (en) | 2024-11-22 |
| CN119004603B true CN119004603B (en) | 2025-03-14 |
Family
ID=93472047
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202411075554.8A Active CN119004603B (en) | 2024-08-06 | 2024-08-06 | Rapid calculation method and system for soil and stone quantity of multi-scale excavation area for distinguishing lithology |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN119004603B (en) |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111191307A (en) * | 2019-12-16 | 2020-05-22 | 福建建工集团有限责任公司 | Earthwork virtual construction method based on BIM + GIS technology |
| CN111898869A (en) * | 2020-07-07 | 2020-11-06 | 舜元建设(集团)有限公司 | Earth and stone measurement and calculation method combining oblique photography, RTK (real time kinematic) and BIM (building information modeling) technologies |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112734929B (en) * | 2020-12-28 | 2022-06-03 | 三峡大学 | Method for calculating excavation volume of complex earth and rockfill dam earth stock ground based on grid subdivision algorithm |
-
2024
- 2024-08-06 CN CN202411075554.8A patent/CN119004603B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111191307A (en) * | 2019-12-16 | 2020-05-22 | 福建建工集团有限责任公司 | Earthwork virtual construction method based on BIM + GIS technology |
| CN111898869A (en) * | 2020-07-07 | 2020-11-06 | 舜元建设(集团)有限公司 | Earth and stone measurement and calculation method combining oblique photography, RTK (real time kinematic) and BIM (building information modeling) technologies |
Also Published As
| Publication number | Publication date |
|---|---|
| CN119004603A (en) | 2024-11-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US10795053B2 (en) | Systems and methods of multi-scale meshing for geologic time modeling | |
| US8050892B2 (en) | Method, device, computer program and data carrier for modeling a multidimensional heterogeneous structure, using a digital processing unit | |
| CN109508508B (en) | Surface mine governance investigation design method | |
| CN104952107A (en) | Three-dimensional bridge reconstruction method based on vehicle-mounted LiDAR point cloud data | |
| CN111595403B (en) | Engineering earthwork measuring method based on point cloud measuring technology | |
| CN116152461A (en) | Geological modeling method, device, computer equipment and computer readable storage medium | |
| CN111797454A (en) | Foundation pit earth volume calculation method based on digital informatization technology | |
| CN112734929B (en) | Method for calculating excavation volume of complex earth and rockfill dam earth stock ground based on grid subdivision algorithm | |
| CN117351165B (en) | Construction and calculation method of three-dimensional stratum model based on ArcGIS | |
| US20080270031A1 (en) | Geospatial modeling system providing data thinning of geospatial data points and related methods | |
| CN107122522B (en) | Digital elevation model precision evaluation method based on fill-cut analysis | |
| CN114332413A (en) | A geological body modeling method and device based on sliding kriging interpolation | |
| CN115797288B (en) | Method for calculating filling and excavating amount based on ground point cloud data | |
| CN107886573B (en) | A three-dimensional finite element mesh generation method for slopes under complex geological conditions | |
| US7590512B2 (en) | System and method for deformation analysis using inverse pre-deformation of finite element mesh | |
| CN114491721B (en) | Photovoltaic module arrangement method and device | |
| CN119004603B (en) | Rapid calculation method and system for soil and stone quantity of multi-scale excavation area for distinguishing lithology | |
| CN107507179A (en) | Rock And Soil quantitative analysis method based on GOCAD | |
| CN117392335B (en) | Conversion method from three-dimensional surface element geologic model to three-dimensional surface element geologic model | |
| CN116976112B (en) | River channel large section data extraction method, system, electronic equipment and storage medium | |
| CN101281662A (en) | A three-dimensional stratigraphic subdivision method of triangular prism and excavation earthwork estimation method | |
| CN117760514A (en) | Three-dimensional model earthwork calculation method, device, medium and equipment based on Boolean operation | |
| CN115630478A (en) | Discrete fracture network generation method based on iterative function system | |
| Tonini et al. | Three-dimensional reconstruction of the Carrara Syncline (Apuane Alps, Italy): An approach to reconstruct and control a geological model using only field survey data | |
| CN106875480A (en) | A kind of method of city three-dimensional data tissue |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |