CN112906124B - Method, device, equipment and medium for calculating earthwork - Google Patents
Method, device, equipment and medium for calculating earthwork Download PDFInfo
- Publication number
- CN112906124B CN112906124B CN202110360582.4A CN202110360582A CN112906124B CN 112906124 B CN112906124 B CN 112906124B CN 202110360582 A CN202110360582 A CN 202110360582A CN 112906124 B CN112906124 B CN 112906124B
- Authority
- CN
- China
- Prior art keywords
- vertexes
- dimensional model
- live
- action
- vertex
- 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
Images
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
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- 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/10—Constructive solid geometry [CSG] using solid primitives, e.g. cylinders, cubes
-
- 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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Software Systems (AREA)
- Computer Graphics (AREA)
- General Engineering & Computer Science (AREA)
- Architecture (AREA)
- Civil Engineering (AREA)
- Structural Engineering (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Image Processing (AREA)
- Processing Or Creating Images (AREA)
Abstract
The application discloses a method, a device, equipment and a medium for calculating earthwork volume, which are used for obtaining a real-scene three-dimensional model of a region to be measured, triangulating the real-scene three-dimensional model and obtaining the vertex of the real-scene three-dimensional model; filtering non-ground points in the vertexes of the live-action three-dimensional model to obtain a filtered live-action three-dimensional model; when the filtered live-action three-dimensional model has a cavity area, interpolating the cavity area, and performing surface fitting on the interpolated cavity area by a least square method to obtain a fitting surface; and extracting the filtered live three-dimensional model and the elevation points of the fitting curved surface, and calculating the earthwork of the region to be measured based on the elevation points. The method improves the technical problems that the existing earthwork calculation method needs to measure real plane coordinate data and elevation data of a large number of field measurement points on original topography in real time, the measurement process is complex, and the calculation efficiency is low.
Description
Technical Field
The present application relates to the field of earthwork computation, and in particular, to an earthwork computation method, apparatus, device, and medium.
Background
The calculation of the earthwork is an important step of engineering construction, and the engineering design stage must budget the earthwork, which is directly related to engineering cost and scheme preference. When the earth and stone amount calculation is carried out, the original topography must be mapped, the mapping scale is determined according to the earth amount calculation precision and the topography condition, and the original topography elevation, the original design elevation and the changed design elevation must be marked.
Disclosure of Invention
The application provides an earthwork calculation method, an earthwork calculation device, earthwork calculation equipment and earthwork calculation media, which are used for improving the technical problems that the prior earthwork calculation method needs to measure real plane coordinate data and elevation data of a large number of field measurement points in real time for original topography, and has complex measurement process and low calculation efficiency.
In view of this, a first aspect of the present application provides an earthmoving mass calculating method including:
acquiring a real-scene three-dimensional model of a region to be measured, triangulating the real-scene three-dimensional model, and acquiring vertexes of triangles to obtain vertexes of the real-scene three-dimensional model;
filtering non-ground points in the vertexes of the live-action three-dimensional model to obtain a filtered live-action three-dimensional model;
when the filtered live-action three-dimensional model has a cavity area, interpolating the cavity area, and performing surface fitting on the interpolated cavity area by a least square method to obtain a fitting surface;
and extracting the filtered live three-dimensional model and the elevation points of the fitting curved surface, and calculating the earthwork of the region to be measured based on the elevation points.
Optionally, the filtering the non-ground points in the vertex of the live-action three-dimensional model to obtain the filtered live-action three-dimensional model includes:
extracting scanning lines from vertexes of the live-action three-dimensional model;
acquiring the lowest vertex and the highest vertex on each scanning line, wherein the lowest vertex is a ground point, and the highest vertex is a non-ground point;
coarse classification is carried out on the vertexes on each scanning line based on the lowest vertexes on each scanning line, so that a non-ground point set is obtained;
setting contour lines with preset sizes by taking the highest vertex on each scanning line as a central point, adding the vertices in the contour lines into the non-ground point set, wherein the contour lines are rectangular or irregular polygons;
and filtering the non-ground points in the vertexes of the live-action three-dimensional model based on the non-ground point set to obtain the filtered live-action three-dimensional model.
Optionally, the coarse classification of the vertices on each scan line based on the lowest vertex on each scan line, to obtain a non-ground point set, includes:
taking the lowest vertex on each scanning line as a target vertex;
calculating gradient values of the target vertexes and adjacent vertexes of the target vertexes on the scanning lines, and judging whether the adjacent vertexes of the target vertexes are non-ground points or not based on the gradient values;
and returning the adjacent vertexes of the target vertexes serving as new target vertexes to calculate gradient values of the target vertexes and adjacent vertexes of the target vertexes on the scanning lines, and judging whether the adjacent vertexes of the target vertexes are non-ground points or not based on the gradient values until the vertexes on the scanning lines are traversed, so as to obtain a non-ground point set.
Optionally, the taking the lowest vertex on each scan line as a target vertex further includes:
detecting whether an outlier exists in the vertexes of each scanning line, and if yes, eliminating the outlier.
Optionally, when the filtered live-action three-dimensional model has a cavity area, interpolating the cavity area includes:
when the filtered live-action three-dimensional model has a cavity area, constructing a triangular mesh by boundary points of the cavity area and vertexes in an area with a preset size around the cavity area;
and inserting a point into the cavity area by a triangular net interpolation method, wherein when the boundary characteristic of the cavity area is mainly a straight line, the elevation value of the inserted point is the elevation average value of the boundary point of the cavity area, and when the boundary characteristic of the cavity area is mainly a curve, the elevation value of the inserted point is obtained by the triangular net interpolation nearest to the point.
Optionally, the calculating the earthwork of the area to be measured based on the elevation point includes:
generating a triangular net based on the elevation points to obtain a plurality of triangular cones; or dividing the region to be measured into squares, obtaining the elevation value of each square vertex through interpolation, and generating a plurality of triangular pyramids based on each square vertex and the elevation point;
and calculating the volume of each triangular pyramid, and counting the volumes of all the triangular pyramids to obtain the earthwork of the region to be measured.
A second aspect of the present application provides an earthwork calculation apparatus, comprising:
the acquisition unit is used for acquiring a real-scene three-dimensional model of the area to be measured, triangulating the real-scene three-dimensional model, acquiring the vertexes of the triangle, and obtaining the vertexes of the real-scene three-dimensional model;
the filtering unit is used for filtering non-ground points in the vertexes of the live-action three-dimensional model to obtain the filtered live-action three-dimensional model;
the interpolation and fitting unit is used for interpolating the cavity area when the filtered live-action three-dimensional model has the cavity area, and performing surface fitting on the interpolated cavity area by a least square method to obtain a fitting surface;
and the calculation unit is used for extracting the filtered live three-dimensional model and the elevation point of the fitting curved surface, and calculating the earthwork quantity of the region to be measured based on the elevation point.
Optionally, the filtering unit specifically includes:
an extraction subunit, configured to extract a scan line from a vertex of the live-action three-dimensional model;
the acquisition subunit is used for acquiring the lowest vertex and the highest vertex on each scanning line, wherein the lowest vertex is a ground point, and the highest vertex is a non-ground point;
the classifying subunit is used for roughly classifying the vertexes on the scanning lines based on the lowest vertexes on the scanning lines to obtain a non-ground point set;
a setting subunit, configured to set a contour line with a preset size with the highest vertex on each scan line as a center point, and add a vertex in the contour line to the non-ground point set, where the contour line is rectangular or irregular polygon;
and the filtering subunit is used for filtering the non-ground points in the vertexes of the live-action three-dimensional model based on the non-ground point set to obtain the filtered live-action three-dimensional model.
A third aspect of the present application provides an earthwork calculation apparatus comprising a processor and a memory;
the memory is used for storing program codes and transmitting the program codes to the processor;
the processor is configured to execute the method for calculating an earthvolume according to any one of the first aspects according to instructions in the program code.
A fourth aspect of the present application provides a computer-readable storage medium for storing program code for performing any one of the earthmoving methods of the first aspect.
From the above technical scheme, the application has the following advantages:
the application provides an earthwork volume calculating method, which comprises the following steps: acquiring a real-scene three-dimensional model of a region to be measured, triangulating the real-scene three-dimensional model, and acquiring vertexes of triangles to obtain vertexes of the real-scene three-dimensional model; filtering non-ground points in the vertexes of the live-action three-dimensional model to obtain a filtered live-action three-dimensional model; when the filtered live-action three-dimensional model has a cavity area, interpolating the cavity area, and performing surface fitting on the interpolated cavity area by a least square method to obtain a fitting surface; and extracting the filtered live three-dimensional model and the elevation points of the fitting curved surface, and calculating the earthwork of the region to be measured based on the elevation points.
In the method, the real-scene three-dimensional model of the area to be measured is obtained, a series of processing is carried out on the real-scene three-dimensional model, so that non-ground points such as surface buildings, digital woods and the like in the real-scene three-dimensional model are filtered, then elevation data are obtained to calculate the earthwork, a large amount of complex measurement work is not needed, and the calculation efficiency is high; further, when the filtered live-action three-dimensional model has a cavity area, the cavity area is interpolated, and the interpolated cavity area is subjected to surface fitting by a least square method to form a digital elevation model conforming to the shape of a real earth surface, so that the accuracy of an earthwork calculation result is ensured, and the technical problems that the existing earthwork calculation method needs to measure real plane coordinate data and elevation data of a large number of field measurement points in real time on the original topography, and is complex in measurement process and low in calculation efficiency are solved.
Drawings
In order to more clearly illustrate the embodiments of the present application or the technical solutions in the prior art, the drawings that are required in the embodiments or the description of the prior art will be briefly described below, it being obvious that the drawings in the following description are only some embodiments of the present application, and that other drawings may be obtained according to these drawings without inventive faculty for a person skilled in the art.
Fig. 1 is a schematic flow chart of an earthwork calculation method according to an embodiment of the present application;
fig. 2 is a schematic structural diagram of an earthwork calculating device according to an embodiment of the present application.
Detailed Description
The application provides an earthwork calculation method, an earthwork calculation device, earthwork calculation equipment and earthwork calculation media, which are used for improving the technical problems that the prior earthwork calculation method needs to measure real plane coordinate data and elevation data of a large number of field measurement points in real time for original topography, and has complex measurement process and low calculation efficiency.
In order to make the present application solution better understood by those skilled in the art, the following description will clearly and completely describe the technical solution in the embodiments of the present application with reference to the accompanying drawings in the embodiments of the present application, and it is apparent that the described embodiments are only some embodiments of the present application, not all embodiments. All other embodiments, which can be made by one of ordinary skill in the art without undue burden from the present disclosure, are within the scope of the present disclosure.
For ease of understanding, referring to fig. 1, an embodiment of a method for calculating earthmoving mass provided in the present application includes:
and 101, acquiring a real-scene three-dimensional model of the region to be measured, triangulating the real-scene three-dimensional model, and acquiring vertexes of triangles to obtain vertexes of the real-scene three-dimensional model.
Image data of the region to be measured can be acquired through a camera and the like, and then a live-action three-dimensional model of the region to be measured is constructed based on the image data. And triangulating the live-action three-dimensional model, extracting triangular vertexes, and obtaining vertexes of the live-action three-dimensional model. The triangle formed by all the vertexes forms the minimum texture mapping unit of the live-action three-dimensional model, and smooth interpolation or thinning of the vertexes of the triangle can be carried out according to actual needs so as to balance calculation accuracy and calculation time. Wherein, when the interpolation of the triangular net is carried out, a slope threshold value a relative to the horizontal plane is set, a plurality of grades are divided according to the threshold value a, and the interpolation calculation process is carried outIn which different classes correspond to triangle vertices V 1 、V 2 、V 3 There are different weights to ensure that the inserted points are close to the real model.
According to the embodiment of the application, the elevation data is obtained by acquiring the image data of the area to be measured to construct the live-action three-dimensional model, and complex measurement work is not needed.
And 102, filtering non-ground points in the vertexes of the live-action three-dimensional model to obtain the filtered live-action three-dimensional model.
Non-ground points in the vertices are filtered before the amount of earth is calculated. The specific process is as follows:
s1, extracting scanning lines from vertexes of the live-action three-dimensional model.
Scanlines may be extracted from vertices of the live-action three-dimensional model by a scanline method.
S2, acquiring the lowest vertex and the highest vertex on each scanning line, wherein the lowest vertex is a ground point, and the highest vertex is a non-ground point.
S3, roughly classifying the vertexes on each scanning line based on the lowest vertexes on each scanning line to obtain a non-ground point set.
Coarse classification is carried out on the vertexes of all the scanning lines based on the lowest vertexes of all the scanning lines, and the specific process for obtaining the non-ground point set is as follows:
s30, taking the lowest vertex on each scanning line as a target vertex;
before the lowest vertex is obtained, whether an outlier exists in the vertices on each scanning line or not can be detected according to an outlier detection algorithm, if yes, the outlier is removed, and if not, the lowest vertex on each scanning line is obtained.
S31, calculating gradient values of a target vertex and adjacent vertices of the target vertex on each scanning line, and judging whether the adjacent vertices of the target vertex are non-ground points or not based on the gradient values;
and calculating gradient values of the target vertex and adjacent vertices of the target vertex on each scanning line, and judging whether the adjacent vertices of the target vertex are ground points or non-ground points by comparing the gradient values with gradient threshold values, wherein the calculation of the gradient values belongs to the prior art, and no repeated description is given here.
S32, taking the adjacent vertexes of the target vertexes as new target vertexes, returning to the step S31 until the vertexes on each scanning line are traversed, and obtaining a non-ground point set.
S4, setting contour lines with preset sizes by taking the highest vertex on each scanning line as a central point, adding the vertices in the contour lines into a non-ground point set, wherein the contour lines are rectangular or irregular polygons;
after the above-mentioned coarse classification, fine filtering may be further performed, and the highest vertex on each scan line is used as a marker point, and this point is usually the center point position of a building or a tree crown, etc., a rectangle or an irregular polygon is defined to be expanded outwards, the inter-line threshold L and the in-circle point set density u are defined according to contour expansion, non-ground points are filtered again, the inter-line threshold L and the density u are adjusted during filtering, and points meeting the condition (i.e., vertices in the contour) are added to the non-ground point set by adding the outwardly expanded threshold.
And S5, filtering the non-ground points in the vertexes of the live-action three-dimensional model based on the non-ground point set to obtain a filtered live-action three-dimensional model.
And filtering the vertexes of the non-ground points in the live-action three-dimensional model according to the acquired non-ground point set, so as to remove the non-ground points such as ground surface buildings, digital woods and the like, and reserving the ground points.
And 103, when the filtered live-action three-dimensional model has a cavity area, interpolating the cavity area, and performing surface fitting on the interpolated cavity area by a least square method to obtain a fitting surface.
The real-scene three-dimensional model with the non-ground points removed can have a cavity area, and interpolation needs to be carried out on the cavity area. When the filtered live-action three-dimensional model has a cavity area, constructing a triangular mesh by boundary points of the cavity area and vertexes in an area with a preset size around the cavity area; and inserting a point into the cavity area by a triangular net interpolation method, wherein when the boundary characteristic of the cavity area is mainly a straight line, the elevation value of the inserted point is the elevation average value of the boundary point of the cavity area, and when the boundary characteristic of the cavity area is mainly a curve, the elevation value of the inserted point is obtained by the triangular net interpolation nearest to the point. There are two methods for interpolating the hole area. Selecting the change type of the elevation in the cavity area according to the boundary characteristics of the cavity area, and considering the elevation value in the cavity area as the average value if the boundary characteristics are mainly straight lines; if the boundary features are mainly curves, the elevation value of the point inserted into the cavity area is obtained by interpolation of the nearest triangular net of the point. And finally, performing surface fitting by adopting a least square method to form a surface which is close to the real topography after the non-ground points are removed.
And 104, extracting the filtered live-action three-dimensional model and the elevation points of the fitting curved surface, and calculating the earthwork of the region to be measured based on the elevation points.
And extracting the elevation points according to the designated grid width according to the storage mode of the three-dimensional model when extracting the filtered live-action three-dimensional model and the Gao Chengdian of the fitting curved surface. In order to improve the extraction efficiency, the elevation points can be extracted in a region-by-region mode in a multithreading mode, and finally the elevation points are combined.
There are two methods for calculating the earth volume of the area to be measured based on the elevation points, in one embodiment, a triangle net is generated based on the elevation points, and a plurality of triangular cones are obtained; and calculating the volume of each triangular pyramid, and counting the volumes of all triangular pyramids to obtain the earth volume of the area to be measured. And generating a triangular net according to the obtained elevation points, so that the topography of the whole area to be measured forms a set formed by triangular pyramid combination. And determining a zero plane according to the design elevation given by a user, performing volume calculation by using "+" to represent filling, performing volume calculation by using "-" to represent excavation, and finally counting the excavation filling amounts of all triangular pyramids to obtain the earthwork amount of the region to be measured.
In another embodiment, dividing the area to be measured into squares, obtaining the elevation value of each square vertex through interpolation, and generating a plurality of triangular pyramids based on each square vertex and the elevation point; and calculating the volume of each triangular pyramid, and counting the volumes of all triangular pyramids to obtain the earth volume of the area to be measured. Dividing the area to be measured according to preset sizes such as 20m or 40m to form squares, generating elevation values of the top points of the squares through interpolation, generating triangular pyramids based on the top points of the squares and the elevation points extracted in the previous steps, performing volume calculation by using "+" to represent filling, performing volume calculation by using "-" to represent cutting, and finally counting the cutting and filling amounts of all the triangular pyramids to obtain the earthwork amount of the area to be measured.
In the embodiment of the application, the real-scene three-dimensional model of the area to be measured is obtained, a series of processing is carried out on the real-scene three-dimensional model, so that non-ground points such as surface building, digital wood and the like in the real-scene three-dimensional model are filtered, then elevation data are obtained to calculate the earthwork, a large amount of complex measurement work is not needed, and the calculation efficiency is high; further, when the filtered live-action three-dimensional model has a cavity area, the cavity area is interpolated, and the interpolated cavity area is subjected to surface fitting by a least square method to form a digital elevation model conforming to the shape of a real earth surface, so that the accuracy of an earthwork calculation result is ensured, and the technical problems that the existing earthwork calculation method needs to measure real plane coordinate data and elevation data of a large number of field measurement points in real time on the original topography, and is complex in measurement process and low in calculation efficiency are solved.
The foregoing is one embodiment of an earthwork calculating method provided herein, and the following is one embodiment of an earthwork calculating apparatus provided herein.
Referring to fig. 2, an apparatus for calculating earthwork according to an embodiment of the present application includes:
the acquisition unit is used for acquiring a real-scene three-dimensional model of the area to be measured, triangulating the real-scene three-dimensional model, acquiring the vertexes of the triangle, and acquiring the vertexes of the real-scene three-dimensional model;
the filtering unit is used for filtering the non-ground points in the vertexes of the live-action three-dimensional model to obtain a filtered live-action three-dimensional model;
the interpolation and fitting unit is used for interpolating the hollow area when the filtered live-action three-dimensional model has the hollow area, and performing surface fitting on the interpolated hollow area by a least square method to obtain a fitting surface;
the calculation unit is used for extracting the filtered live three-dimensional model and the elevation points of the fitting curved surface, and calculating the earthwork of the region to be measured based on the elevation points.
As a further improvement, the filter unit comprises in particular:
an extraction subunit, configured to extract a scan line from a vertex of the live-action three-dimensional model;
the acquisition subunit is used for acquiring the lowest vertex and the highest vertex on each scanning line, wherein the lowest vertex is a ground point, and the highest vertex is a non-ground point;
the classifying subunit is used for roughly classifying the vertexes on each scanning line based on the lowest vertexes on each scanning line to obtain a non-ground point set;
the setting subunit is used for setting a contour line with a preset size by taking the highest vertex on each scanning line as a central point, adding the vertex in the contour line into a non-ground point set, and enabling the contour line to be rectangular or irregular polygonal;
and the filtering subunit is used for filtering the non-ground points in the vertexes of the live-action three-dimensional model based on the non-ground point set to obtain a filtered live-action three-dimensional model.
As a further improvement, the classification subunit is specifically configured to:
taking the lowest vertex on each scanning line as a target vertex;
calculating gradient values of a target vertex and adjacent vertices of the target vertex on each scanning line, and judging whether the adjacent vertices of the target vertex are non-ground points or not based on the gradient values;
and returning the adjacent vertexes of the target vertexes to calculate gradient values of the target vertexes on each scanning line and the adjacent vertexes of the target vertexes, and judging whether the adjacent vertexes of the target vertexes are non-ground points or not based on the gradient values until the vertexes on each scanning line are traversed, so as to obtain a non-ground point set.
As a further improvement, the filter unit further comprises:
and the detection subunit is used for detecting whether the outliers exist in the vertexes of the scanning lines, and if yes, the outliers are removed.
In the embodiment of the application, the real-scene three-dimensional model of the area to be measured is obtained, a series of processing is carried out on the real-scene three-dimensional model, so that non-ground points such as surface building, digital wood and the like in the real-scene three-dimensional model are filtered, then elevation data are obtained to calculate the earthwork, a large amount of complex measurement work is not needed, and the calculation efficiency is high; further, when the filtered live-action three-dimensional model has a cavity area, the cavity area is interpolated, and the interpolated cavity area is subjected to surface fitting by a least square method to form a digital elevation model conforming to the shape of a real earth surface, so that the accuracy of an earthwork calculation result is ensured, and the technical problems that the existing earthwork calculation method needs to measure real plane coordinate data and elevation data of a large number of field measurement points in real time on the original topography, and is complex in measurement process and low in calculation efficiency are solved.
The embodiment of the application also provides an earthwork calculation device, which comprises a processor and a memory;
the memory is used for storing the program codes and transmitting the program codes to the processor;
the processor is configured to execute the earthmoving mass calculation method in the foregoing method embodiment according to instructions in the program code.
The present application also provides a computer-readable storage medium storing program code for executing the earthmoving mass calculating method in the foregoing method embodiment.
It will be clear to those skilled in the art that, for convenience and brevity of description, specific working procedures of the apparatus and units described above may refer to corresponding procedures in the foregoing method embodiments, which are not described herein again.
In the several embodiments provided in this application, it should be understood that the disclosed apparatus and method may be implemented in other ways. For example, the apparatus embodiments described above are merely illustrative, e.g., the division of the units is merely a logical function division, and there may be additional divisions when actually implemented, e.g., multiple units or components may be combined or integrated into another system, or some features may be omitted or not performed. Alternatively, the coupling or direct coupling or communication connection shown or discussed with each other may be an indirect coupling or communication connection via some interfaces, devices or units, which may be in electrical, mechanical or other form.
The units described as separate units may or may not be physically separate, and units shown as units may or may not be physical units, may be located in one place, or may be distributed on a plurality of network units. Some or all of the units may be selected according to actual needs to achieve the purpose of the solution of this embodiment.
In addition, each functional unit in each embodiment of the present application may be integrated in one processing unit, or each unit may exist alone physically, or two or more units may be integrated in one unit. The integrated units may be implemented in hardware or in software functional units.
The integrated units, if implemented in the form of software functional units and sold or used as stand-alone products, may be stored in a computer readable storage medium. Based on such understanding, the technical solution of the present application may be embodied in essence or a part contributing to the prior art or all or part of the technical solution in the form of a software product stored in a storage medium, including several instructions to execute all or part of the steps of the methods described in the embodiments of the present application by a computer device (which may be a personal computer, a server, or a network device, etc.). And the aforementioned storage medium includes: u disk, mobile hard disk, read-Only Memory (ROM), random access Memory (RandomAccess Memory, RAM), magnetic disk or optical disk, etc.
The above embodiments are merely for illustrating the technical solution of the present application, and not for limiting the same; although the present application has been described in detail with reference to the foregoing embodiments, it should be understood by those of ordinary skill in the art that: the technical scheme described in the foregoing embodiments can be modified or some technical features thereof can be replaced by equivalents; such modifications and substitutions do not depart from the spirit and scope of the corresponding technical solutions.
Claims (9)
1. An earthwork calculation method, comprising:
acquiring a real-scene three-dimensional model of a region to be measured, triangulating the real-scene three-dimensional model, and acquiring vertexes of triangles to obtain vertexes of the real-scene three-dimensional model;
filtering non-ground points in the vertexes of the live-action three-dimensional model to obtain a filtered live-action three-dimensional model;
filtering non-ground points in the vertexes of the live-action three-dimensional model to obtain the filtered live-action three-dimensional model, wherein the filtering comprises the following steps:
extracting scanning lines from vertexes of the live-action three-dimensional model;
acquiring the lowest vertex and the highest vertex on each scanning line, wherein the lowest vertex is a ground point, and the highest vertex is a non-ground point;
coarse classification is carried out on the vertexes on each scanning line based on the lowest vertexes on each scanning line, so that a non-ground point set is obtained;
setting contour lines with preset sizes by taking the highest vertex on each scanning line as a central point, adding the vertices in the contour lines into the non-ground point set, wherein the contour lines are rectangular or irregular polygons;
filtering non-ground points in the vertexes of the live-action three-dimensional model based on the non-ground point set to obtain a filtered live-action three-dimensional model;
when the filtered live-action three-dimensional model has a cavity area, interpolating the cavity area, and performing surface fitting on the interpolated cavity area by a least square method to obtain a fitting surface;
and extracting the filtered live three-dimensional model and the elevation points of the fitting curved surface, and calculating the earthwork of the region to be measured based on the elevation points.
2. The method of claim 1, wherein the coarsely classifying vertices on each of the scan lines based on the lowest vertex on each of the scan lines to obtain a set of non-ground points, comprises:
taking the lowest vertex on each scanning line as a target vertex;
calculating gradient values of the target vertexes and adjacent vertexes of the target vertexes on the scanning lines, and judging whether the adjacent vertexes of the target vertexes are non-ground points or not based on the gradient values;
and returning the adjacent vertexes of the target vertexes serving as new target vertexes to calculate gradient values of the target vertexes and adjacent vertexes of the target vertexes on the scanning lines, and judging whether the adjacent vertexes of the target vertexes are non-ground points or not based on the gradient values until the vertexes on the scanning lines are traversed, so as to obtain a non-ground point set.
3. The method of calculating an earthvolume according to claim 2, wherein said setting the lowest vertex on each of said scan lines as a target vertex, further comprises:
detecting whether an outlier exists in the vertexes of each scanning line, and if yes, eliminating the outlier.
4. The earth volume calculation method according to claim 1, wherein, when the filtered live-action three-dimensional model has a void area, interpolating the void area, comprising:
when the filtered live-action three-dimensional model has a cavity area, constructing a triangular mesh by boundary points of the cavity area and vertexes in an area with a preset size around the cavity area;
and inserting a point into the cavity area by a triangular net interpolation method, wherein when the boundary characteristic of the cavity area is mainly a straight line, the elevation value of the inserted point is the elevation average value of the boundary point of the cavity area, and when the boundary characteristic of the cavity area is mainly a curve, the elevation value of the inserted point is obtained by the triangular net interpolation nearest to the point.
5. The earth volume calculating method according to claim 1, characterized in that the calculating the earth volume of the area to be measured based on the elevation point includes:
generating a triangular net based on the elevation points to obtain a plurality of triangular cones; or dividing the region to be measured into squares, obtaining the elevation value of each square vertex through interpolation, and generating a plurality of triangular pyramids based on each square vertex and the elevation point;
and calculating the volume of each triangular pyramid, and counting the volumes of all the triangular pyramids to obtain the earthwork of the region to be measured.
6. An earthwork calculation apparatus, comprising:
the acquisition unit is used for acquiring a real-scene three-dimensional model of the area to be measured, triangulating the real-scene three-dimensional model, acquiring the vertexes of the triangle, and obtaining the vertexes of the real-scene three-dimensional model;
the filtering unit is used for filtering non-ground points in the vertexes of the live-action three-dimensional model to obtain the filtered live-action three-dimensional model;
the filter unit is specifically used for:
extracting scanning lines from vertexes of the live-action three-dimensional model;
acquiring the lowest vertex and the highest vertex on each scanning line, wherein the lowest vertex is a ground point, and the highest vertex is a non-ground point;
coarse classification is carried out on the vertexes on each scanning line based on the lowest vertexes on each scanning line, so that a non-ground point set is obtained;
setting contour lines with preset sizes by taking the highest vertex on each scanning line as a central point, adding the vertices in the contour lines into the non-ground point set, wherein the contour lines are rectangular or irregular polygons;
filtering non-ground points in the vertexes of the live-action three-dimensional model based on the non-ground point set to obtain a filtered live-action three-dimensional model;
the interpolation and fitting unit is used for interpolating the cavity area when the filtered live-action three-dimensional model has the cavity area, and performing surface fitting on the interpolated cavity area by a least square method to obtain a fitting surface;
and the calculation unit is used for extracting the filtered live three-dimensional model and the elevation point of the fitting curved surface, and calculating the earthwork quantity of the region to be measured based on the elevation point.
7. The earth volume computing device of claim 6, wherein the filtering unit specifically comprises:
an extraction subunit, configured to extract a scan line from a vertex of the live-action three-dimensional model;
the acquisition subunit is used for acquiring the lowest vertex and the highest vertex on each scanning line, wherein the lowest vertex is a ground point, and the highest vertex is a non-ground point;
the classifying subunit is used for roughly classifying the vertexes on the scanning lines based on the lowest vertexes on the scanning lines to obtain a non-ground point set;
a setting subunit, configured to set a contour line with a preset size with the highest vertex on each scan line as a center point, and add a vertex in the contour line to the non-ground point set, where the contour line is rectangular or irregular polygon;
and the filtering subunit is used for filtering the non-ground points in the vertexes of the live-action three-dimensional model based on the non-ground point set to obtain the filtered live-action three-dimensional model.
8. An earthwork computation apparatus, characterized in that the apparatus comprises a processor and a memory;
the memory is used for storing program codes and transmitting the program codes to the processor;
the processor is configured to perform the earthmoving computation method of any one of claims 1-5 according to instructions in the program code.
9. A computer-readable storage medium storing a program code for performing the earthmoving mass calculating method as recited in any one of claims 1-5.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110360582.4A CN112906124B (en) | 2021-04-02 | 2021-04-02 | Method, device, equipment and medium for calculating earthwork |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110360582.4A CN112906124B (en) | 2021-04-02 | 2021-04-02 | Method, device, equipment and medium for calculating earthwork |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112906124A CN112906124A (en) | 2021-06-04 |
| CN112906124B true CN112906124B (en) | 2023-06-13 |
Family
ID=76109880
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110360582.4A Active CN112906124B (en) | 2021-04-02 | 2021-04-02 | Method, device, equipment and medium for calculating earthwork |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112906124B (en) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114527461A (en) * | 2021-12-22 | 2022-05-24 | 上海振懋机电科技有限公司 | Space measuring device and measuring method thereof |
| CN115482269B (en) * | 2022-09-22 | 2023-05-09 | 佳都科技集团股份有限公司 | Method and device for calculating earthwork, terminal equipment and storage medium |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109191581A (en) * | 2018-08-29 | 2019-01-11 | 李东帅 | A kind of Earth Volume of Road Engineering amount high-resolution method based on BIM technology |
Family Cites Families (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP3407883B2 (en) * | 2000-09-11 | 2003-05-19 | 国際航業株式会社 | Earthwork design support system using digital map data |
| JP2002203234A (en) * | 2000-12-28 | 2002-07-19 | Inst Of Physical & Chemical Res | Method and device for preparing contour line data |
| CN109658477A (en) * | 2017-10-12 | 2019-04-19 | 西南科技大学 | A kind of DEM generating algorithm based on LIDAR data |
| CN110285792B (en) * | 2019-07-02 | 2021-06-01 | 山东省交通规划设计院集团有限公司 | Fine grid earthwork metering method for unmanned aerial vehicle oblique photography |
| CN110657855B (en) * | 2019-09-23 | 2021-06-04 | 广州市中海达测绘仪器有限公司 | Method and device for measuring earth volume, computer equipment and storage medium |
| CN111325837B (en) * | 2020-01-23 | 2022-08-09 | 江西理工大学 | Side slope DEM generation method based on ground three-dimensional laser point cloud |
-
2021
- 2021-04-02 CN CN202110360582.4A patent/CN112906124B/en active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109191581A (en) * | 2018-08-29 | 2019-01-11 | 李东帅 | A kind of Earth Volume of Road Engineering amount high-resolution method based on BIM technology |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112906124A (en) | 2021-06-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN112906124B (en) | Method, device, equipment and medium for calculating earthwork | |
| US7272264B2 (en) | System and method for hole filling in 3D models | |
| Li et al. | An improved building boundary extraction algorithm based on fusion of optical imagery and LIDAR data | |
| Truong-Hong et al. | Octree-based, automatic building facade generation from LiDAR data | |
| CN108225334B (en) | Positioning method and device based on three-dimensional live-action data | |
| JP2011501301A (en) | Geospatial modeling system and related methods using multiple sources of geographic information | |
| Zhu et al. | Seed point selection method for triangle constrained image matching propagation | |
| CN111784840B (en) | LOD (line-of-sight) level three-dimensional data singulation method and system based on vector data automatic segmentation | |
| CN111323788B (en) | Building change monitoring method and device and computer equipment | |
| Zandbergen | Accuracy considerations in the analysis of depressions in medium-resolution Lidar DEMs | |
| Tseng et al. | Extraction of building boundary lines from airborne lidar point clouds | |
| CN108198194A (en) | Oblique photography single object extraction method based on lateral multi-plane cutting | |
| KR101547943B1 (en) | Apparatus and Method for Watershed Extraction using Polygon and Digital Elevation Model | |
| Pentland | Fractal-Based Description. | |
| CN112687007B (en) | Stereoscopic grid chart generation method based on LOD technology | |
| CN111583406A (en) | Pole tower foot base point coordinate calculation method and device and terminal equipment | |
| CN116051777A (en) | Super high-rise building extraction method, apparatus and readable storage medium | |
| CN112700538B (en) | LOD generation method and system | |
| Yakar et al. | The effect of grid resolution in defining terrain surface | |
| CN114019532A (en) | Project progress checking method and device | |
| Dakowicz et al. | Extracting meaningful slopes from terrain contours | |
| Rutzinger et al. | Object-based analysis of airborne laser scanning data for natural hazard purposes using open source components | |
| KR101673442B1 (en) | The method and apparatus for remeshing visual hull approximation by DBSS(displaced butterfly subdivision surface) | |
| KR102531766B1 (en) | Method and Apparatus for Segmenting Density-based Tetrahedron Model for Level-Of-Detail Generation | |
| CN111354081A (en) | Large-scale three-dimensional forest drawing method |
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 |