CN112365534A - Large coal pile volume measurement method based on monocular camera three-dimensional reconstruction - Google Patents

Large coal pile volume measurement method based on monocular camera three-dimensional reconstruction Download PDF

Info

Publication number
CN112365534A
CN112365534A CN202011263711.XA CN202011263711A CN112365534A CN 112365534 A CN112365534 A CN 112365534A CN 202011263711 A CN202011263711 A CN 202011263711A CN 112365534 A CN112365534 A CN 112365534A
Authority
CN
China
Prior art keywords
volume
point
coal pile
dimensional reconstruction
monocular camera
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.)
Granted
Application number
CN202011263711.XA
Other languages
Chinese (zh)
Other versions
CN112365534B (en
Inventor
张勇
黄诚
胡臻
朱华明
王思彪
全文涛
王健
董延超
沈润杰
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shanghai Waigaoqiao Power Generation Co ltd
Suzhou Tonggang Uav Research Institute Co ltd
Shanghai Minghua Power Technology Co ltd
Original Assignee
Shanghai Waigaoqiao Power Generation Co ltd
Suzhou Tonggang Uav Research Institute Co ltd
Shanghai Minghua Power Technology Co ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Shanghai Waigaoqiao Power Generation Co ltd, Suzhou Tonggang Uav Research Institute Co ltd, Shanghai Minghua Power Technology Co ltd filed Critical Shanghai Waigaoqiao Power Generation Co ltd
Priority to CN202011263711.XA priority Critical patent/CN112365534B/en
Publication of CN112365534A publication Critical patent/CN112365534A/en
Application granted granted Critical
Publication of CN112365534B publication Critical patent/CN112365534B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/60Analysis of geometric attributes
    • G06T7/62Analysis of geometric attributes of area, perimeter, diameter or volume
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T15/003D [Three Dimensional] image rendering
    • G06T15/10Geometric effects
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00Three dimensional [3D] modelling, e.g. data description of 3D objects
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10004Still image; Photographic image
    • G06T2207/10012Stereo images

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Graphics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Software Systems (AREA)
  • Length Measuring Devices By Optical Means (AREA)

Abstract

The invention relates to a monocular camera three-dimensional reconstruction-based large-scale coal pile volume measurement method, which uses an unmanned aerial vehicle to carry a monocular camera, realizes autonomous positioning and navigation of the unmanned aerial vehicle, samples data of a coal pile and reconstructs the coal pile in a three-dimensional manner to calculate the volume, and arranges a two-dimensional code for coal pile scale reduction in the coal pile. Compared with the prior art, the method has the advantages that the automation of the whole operation process can be realized, more manpower and material resources are saved, the real volume of the coal pile can be obtained in a short time, the measurement precision can be ensured to meet the actual production requirement, and the like.

Description

Large coal pile volume measurement method based on monocular camera three-dimensional reconstruction
Technical Field
The invention relates to a large coal pile volume measuring method, in particular to a large coal pile volume measuring method based on monocular camera three-dimensional reconstruction.
Background
In a thermal power plant, the volume measurement of a large coal pile is often required, and the mapping is difficult. The traditional method adopts a manual measurement mode, a large amount of manpower, material resources and time are needed to shape the coal pile, and then the volume is estimated through manual measurement. At present, some superior inventions can more accurately acquire the volume of the coal pile.
The retrieved Chinese patent publication No. CN1022441C discloses an automatic measuring device for the volume of a large coal pile, which carries out point-by-point sampling and distance measurement on the surface height of the large coal pile, and comprises a mechanical scanning system and a photoelectric detection system, the volume of the large coal pile can be rapidly measured, but the precision of the method changes along with the number of sampling points, if a more accurate true value is to be obtained, multiple times of sampling are needed, and the device carrying a laser and a camera is heavy in machinery and difficult to operate and maintain.
Chinese patent publication No. CN 101858730 a discloses a method for automatically measuring the volume of a coal pile and a special device thereof, wherein a cruise trolley is used to carry a laser range finder, and the laser range finder is used to perform range finding sampling at a plurality of positions on the top of the coal pile, so as to establish a three-dimensional map of the coal pile scene for volume measurement.
Chinese patent publication No. CN 106875440 a discloses a method and an apparatus for calculating a volume of a coal pile, in which a large amount of data is sampled by a camera to perform three-dimensional reconstruction using image data satisfying requirements, so as to accurately obtain a three-dimensional model of the coal pile for volume estimation, but an actual coal pile scale is required when estimating the volume, and when using a monocular camera for sampling reconstruction, a true scale of a reconstruction model cannot be obtained, which causes a certain influence on volume estimation of the coal pile.
Disclosure of Invention
The invention aims to overcome the defects of the prior art and provide a large coal pile volume measuring method based on monocular camera three-dimensional reconstruction.
The purpose of the invention can be realized by the following technical scheme:
a monocular camera three-dimensional reconstruction-based large coal pile volume measurement method is characterized in that an unmanned aerial vehicle is used for carrying a monocular camera, autonomous positioning and navigation of the unmanned aerial vehicle are achieved, data sampling and three-dimensional reconstruction are conducted on a coal pile to calculate the volume, and meanwhile two-dimensional codes for coal pile scale reduction are arranged in the coal pile.
Preferably, the method comprises the following steps:
step 1) data collection;
step 2), obtaining an effective 3D point set;
step 3), performing spatial transformation, namely projecting all three-dimensional point clouds onto a two-dimensional plane to obtain a 2D point set;
step 4), obtaining a Delaunay triangulation network, and triangulating the 2D point set to obtain the Delaunay triangulation network;
and 5) converting the solution of the point cloud volume after three-dimensional reconstruction into the solution of a plurality of triangular prism integrals through the steps, and eliminating the triangular prism volume values of which partial deviations exceed a set threshold value through an algorithm.
Preferably, the step 1) data collection specifically includes:
and continuously acquiring a scene by adopting a monocular camera, taking the acquired image as input, and obtaining point cloud of the scene through three-dimensional reconstruction.
Preferably, the step 2) of obtaining the effective 3D point set specifically includes:
in the process of generating the point cloud, the contents of partial non-interesting areas and some discrete error noise points are added into the point cloud result, and the point cloud result is removed by the following method:
201) directly cutting and removing the non-interested region by a three-dimensional space plane equation aiming at the non-interested region;
202) and (4) rejecting random error noise points through an outlier detection algorithm.
Preferably, 201) is specifically:
in the point cloud after three-dimensional reconstruction, only a 3D point set in an interested region needs to be collected, a space tangent plane equation is solved through the joint constraint of the solved horizontal plane equation and the known three-dimensional point coordinates in the scene to cut the scene, and therefore the non-interested region is directly removed.
Preferably, the outlier detection algorithm in 202) is specifically:
for each point, calculate its average distance to all nearby points, define the points whose average distance is outside the standard range as outliers and remove them from the data.
Preferably, the spatial transformation in step 3) is specifically:
301) acquiring normal vectors P (A, B, C) and Q (0,0,1) of an xy plane equation in a point cloud horizontal plane and a reference coordinate system;
302) set all 3D points O (x)0,y0,z0) Projecting the horizontal plane according to the following formula to obtain a new 3D point set P (x)p,yp,zp)
Figure BDA0002775450560000031
Figure BDA0002775450560000032
Figure BDA0002775450560000033
Wherein D is a fourth parameter of the horizontal plane equation;
simultaneously, the distance d between each three-dimensional space point and the equation of the horizontal plane is calculated according to the following formula
Figure BDA0002775450560000034
303) Calculating a rotation angle: the rotation angle θ between the two vectors is calculated according to the following equation
Figure BDA0002775450560000035
Wherein P and Q are two three-dimensional vectors;
304) calculating a rotating shaft: when the rotation angle is calculated, the plane where the angle is located is a plane formed by P and Q, the rotation axis must be perpendicular to the plane, and the vector before rotation is assumed to be a (a)1,a2,a3) After rotation the vector is b (b)1,b2,b3) The rotation axis c (c) is defined by the cross product1,c2,c3) Is composed of
Figure BDA0002775450560000036
305) Converting the rotation angle, the expression of the rotation axis, into a rotation matrix R according to the rodriger rotation formula;
306) set P (x) of 3D pointsp,yp,zp) And right multiplying the rotation matrix R to obtain a series of 3D point sets on the same plane, and extracting x and y values of all the 3D points to form a 2D point set p (x, y).
Preferably, the step 4) of obtaining the Delaunay triangulation network specifically includes:
401) constructing a super triangle, including all scatter points, and putting the super triangle into a triangle linked list;
402) inserting scattered points in the point set in sequence, finding out a triangle (called as an influence triangle of the point) of which the circumscribed circle comprises an insertion point from the triangle linked list, deleting a common edge of the influence triangle, and connecting the insertion point with all vertexes of the influence triangle, thereby completing the insertion of one point in the Delaunay triangle linked list;
403) optimizing the local newly formed triangles according to an optimization criterion, and putting the formed triangles into a Delaunay triangle linked list;
404) loop through 402) until all scatter insertions have been completed.
Preferably, the set of triangles and the distance corresponding to each vertex in the triangles are obtained through the above process, the average value of the distances of the three vertices is taken as the height of the triangular prism, and the area of each triangle is calculated according to the following formula, so that the volume of all the triangular prisms can be obtained:
Figure BDA0002775450560000041
preferably, the triangular prism volume values with partial deviation exceeding the set threshold are rejected by an algorithm by using a box diagram method, specifically:
need to solve for the numberAccording to 5 characteristic data values, including the lower quartile Q1Upper quartile Q3Four-bit distance IQR ═ Q3-Q1And an upper limit Q3+1.5IQR and lower limit Q11.5IQR, eliminating the content outside the upper and lower limits, finally combining the estimated scale and the volume of all effective triangular prisms according to the following formula to estimate the actual volume value of the interested area in the object,
Figure BDA0002775450560000042
wherein K is a scale factor recovered in the point cloud processing process, djIs the length of the side edge of a triangular prism, here
Figure BDA0002775450560000043
The volume is obtained by taking the average length of the side edges of the triangular prism as the height.
Compared with the prior art, the unmanned aerial vehicle carries the monocular camera, autonomous positioning and navigation of the unmanned aerial vehicle are achieved, data sampling and three-dimensional reconstruction are carried out on the coal pile to calculate the volume, a small number of two-dimensional codes need to be arranged in the coal pile, and scale reduction of the coal pile is facilitated. Through actual scene measurement, the invention can realize the automation of the whole operation process, save more manpower and material resources, can obtain the real volume of the coal pile in a short time, and can ensure that the measurement precision meets the actual production requirement.
Drawings
FIG. 1 is a schematic diagram of a three-dimensional reconstruction point cloud result;
FIG. 2 is a schematic diagram of a three-dimensional reconstructed region of interest;
FIG. 3 is a diagram illustrating a point cloud after statistical filtering;
FIG. 4 is a schematic diagram of a Delaunay triangulation network;
FIG. 5 is a box diagram illustration;
FIG. 6 is a flow chart of the present invention.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are some, not all, embodiments of the present invention. All other embodiments, which can be obtained by a person skilled in the art without any inventive step based on the embodiments of the present invention, shall fall within the scope of protection of the present invention.
In order to solve the problem of volume measurement of the coal pile after three-dimensional reconstruction, the invention provides a set of flow combining computer vision and geometric calculation, and the actual volume of the coal pile can be solved on the basis after the horizontal plane and the proportional scale of the point cloud after three-dimensional reconstruction are obtained.
The innovation points of the invention comprise:
1) when the volume of the coal pile is calculated, the calculation of the point cloud volume is converted into the integral of a plurality of triangular prism volumes
2) In the calculation of the coal pile volume, outliers are presented by line graphs.
The technical scheme adopted by the invention for solving the technical problem is as follows:
firstly, projecting all point cloud results on a horizontal plane, and rejecting abnormal points according to the actual space constraint of a scene and a filtering method.
Further, three-dimensional space points are converted into a 2D point set on a horizontal plane through space transformation.
Further, solving a Delaunay triangulation network for all point sets, and calculating the area and the corresponding horizontal height of each triangle to obtain the volume of the triangular prism
And finally, removing abnormal triangular prisms, and summing the integrals of all effective triangular prisms to obtain the total volume.
As shown in fig. 6, the method for measuring the volume of a large coal pile based on monocular camera three-dimensional reconstruction of the present invention specifically comprises the following steps:
1. data collection: the scene can be continuously collected through a monocular camera, the collected images are used as input, point clouds of the scene are obtained through three-dimensional reconstruction, and the point cloud result is shown in fig. 1.
2. Obtaining a valid 3D point set: the contents of partial non-interesting regions and some discrete error noise points are added to the point cloud result in the process of generating the point cloud. Aiming at the two situations, the invention respectively provides the following corresponding solutions:
1) aiming at the non-interested region, the non-interested region is directly cut and removed through a three-dimensional space plane equation;
2) aiming at random error noise points, the random error noise points are removed through an outlier detection algorithm.
For the scheme 1), in the point cloud after the three-dimensional reconstruction, only a 3D point set in the region of interest needs to be collected, as shown in fig. 2, only the region in the red block diagram is the region of interest, and the scene can be cut by solving a space tangent plane equation through the joint constraint of the solved horizontal plane equation and the known three-dimensional point coordinates in the scene.
For scheme 2), the invention uses a statistical filter to perform a statistical analysis on the neighborhood of each point and prunes off some points that do not meet the standard, and the specific method is to calculate the distance distribution from the point to the adjacent point in the input data, and calculate the average distance from each point to all the adjacent points (assuming that the obtained result is a gaussian distribution, and the shape is determined by the mean and standard deviation), then the points whose average distance is outside the standard range can be defined as outliers and removed from the data. After filtering the discrete points in the point cloud by using statistical filtering, as shown in fig. 3.
3. Spatial transformation: through the steps, an effective 3D point set can be obtained, all three-dimensional point clouds need to be projected to a two-dimensional plane in the subsection to obtain a 2D point set, and the simplified steps are as follows:
1) obtaining the normal vectors of the xy plane equation in the point cloud horizontal plane and the reference coordinate system, P (A, B, C) and Q (0,0,1)
2) Set all 3D points O (x)0,y0,z0) Projecting the horizontal plane according to the following formula to obtain a new 3D point set P (x)p,yp,zp)
Figure BDA0002775450560000061
Figure BDA0002775450560000062
Figure BDA0002775450560000063
Meanwhile, the distance between each three-dimensional space point and the equation of the horizontal plane can be calculated according to the following formula
Figure BDA0002775450560000064
3) Calculating a rotation angle: the rotation angle between the two vectors is calculated as follows
Figure BDA0002775450560000065
4) Calculating a rotating shaft: when the rotation angle is calculated, it is known that the plane in which the angle is located is a plane formed by P and Q, and the rotation axis must be perpendicular to the plane. Assume that the pre-rotation vector is a (a)1,a2,a3) After rotation the vector is b (b)1,b2,b3) The rotation axis c (c) is defined by the cross product1,c2,c3) Is composed of
Figure BDA0002775450560000071
5) Converting the rotation angle, the expression of the rotation axis, into a rotation matrix R according to the Rodrigue rotation formula
6) Set P (x) of 3D pointsp,yp,zp) Multiplying the rotation matrix R to the right can obtain a series of 3D point sets on the same plane (all converted point z values are the same), and extracting the x and y values of all the 3D points to form a 2D point set p (x, y).
The 2D point sets converted from the 3D point sets on the same plane and the distance value D corresponding to each 2D point set can be obtained through the steps.
4. Acquiring a Delaunay triangulation network: triangulating the set of 2D points yields a Delaunay triangulation, which is a set of connected but non-overlapping triangles whose circumcircles do not contain any other point of the face region, as shown in fig. 4. The general steps for obtaining the Delaunay triangular mesh are as follows:
1) constructing a super triangle, including all scatter points, and putting the super triangle into a triangle linked list;
2) inserting scattered points in the point set in sequence, finding out a triangle (called as an influence triangle of the point) of which the circumscribed circle comprises an insertion point from the triangle linked list, deleting a common edge of the influence triangle, and connecting the insertion point with all vertexes of the influence triangle, thereby completing the insertion of one point in the Delaunay triangle linked list;
3) optimizing the local newly formed triangles according to an optimization criterion, and putting the formed triangles into a Delaunay triangle linked list;
4) and (5) circularly executing the step 2 until all scatter points are inserted.
Because the Delaunay subdivision has the characteristics of being closest (a triangle is formed by the nearest three points and all line segments do not intersect) and regional (a vertex is newly added, deleted and moved, only the adjacent triangle is influenced), the method is adopted for fitting the plane formed by the point set. The set of triangles and the distance corresponding to each vertex in the triangles can be obtained through the process, the average value of the distances of the three vertices is taken as the height of the triangular prism, and the area of each triangle is calculated according to the following formula, so that the volume of all the triangular prisms can be obtained.
Figure BDA0002775450560000072
5. Integrating triangular prisms: through the above steps, the solution of the three-dimensional reconstructed point cloud volume can be converted into the solution of the plurality of triangular prism integrals, because all of themThe volume values of the triangular prisms are distributed in a certain interval (not completely obeying positive distribution), and the partial triangular prism volume values with large deviation are removed through an algorithm. The invention adopts a box diagram method to deal with the problem, the box diagram is not influenced by abnormal values, the discrete distribution situation of the data can be accurately and stably described, and the cleaning of the data is also facilitated. The boxplot is shown in FIG. 5, and mainly requires solving 5 characteristic data values in the data, including the lower quartile Q1Upper quartile Q3Four-bit distance IQR ═ Q3-Q1And an upper limit Q3+1.5IQR and lower limit Q1And 1.5IQR, eliminating the content outside the upper limit and the lower limit, and finally combining the estimated scale and the volumes of all effective triangular prisms according to the following formula to estimate the actual volume value of the region of interest in the object.
Figure BDA0002775450560000081
While the invention has been described with reference to specific embodiments, the invention is not limited thereto, and various equivalent modifications and substitutions can be easily made by those skilled in the art within the technical scope of the invention. Therefore, the protection scope of the present invention shall be subject to the protection scope of the claims.

Claims (10)

1. A method for measuring the volume of a large coal pile based on monocular camera three-dimensional reconstruction is characterized in that an unmanned aerial vehicle is used for carrying a monocular camera, autonomous positioning and navigation of the unmanned aerial vehicle are achieved, data sampling and three-dimensional reconstruction are conducted on the coal pile to calculate the volume, and meanwhile a two-dimensional code for coal pile scale reduction is arranged in the coal pile.
2. The large coal pile volume measurement method based on monocular camera three-dimensional reconstruction as recited in claim 1, specifically comprising the steps of:
step 1) data collection;
step 2), obtaining an effective 3D point set;
step 3), performing spatial transformation, namely projecting all three-dimensional point clouds onto a two-dimensional plane to obtain a 2D point set;
step 4), obtaining a Delaunay triangulation network, and triangulating the 2D point set to obtain the Delaunay triangulation network;
and 5) converting the solution of the point cloud volume after three-dimensional reconstruction into the solution of a plurality of triangular prism integrals through the steps, and eliminating the triangular prism volume values of which partial deviations exceed a set threshold value through an algorithm.
3. The method for measuring the volume of the large coal pile based on the monocular camera three-dimensional reconstruction as recited in claim 2, wherein the step 1) data collection specifically comprises:
and continuously acquiring a scene by adopting a monocular camera, taking the acquired image as input, and obtaining point cloud of the scene through three-dimensional reconstruction.
4. The method for measuring the volume of the large coal pile based on the monocular camera three-dimensional reconstruction as recited in claim 2, wherein the step 2) of obtaining the effective 3D point set specifically comprises:
in the process of generating the point cloud, the contents of partial non-interesting areas and some discrete error noise points are added into the point cloud result, and the point cloud result is removed by the following method:
201) directly cutting and removing the non-interested region by a three-dimensional space plane equation aiming at the non-interested region;
202) and (4) rejecting random error noise points through an outlier detection algorithm.
5. The method for measuring the volume of the large coal pile based on the monocular camera three-dimensional reconstruction as recited in claim 4, wherein 201) is specifically:
in the point cloud after three-dimensional reconstruction, only a 3D point set in an interested region needs to be collected, a space tangent plane equation is solved through the joint constraint of the solved horizontal plane equation and the known three-dimensional point coordinates in the scene to cut the scene, and therefore the non-interested region is directly removed.
6. The method for measuring the volume of the large coal pile based on the monocular camera three-dimensional reconstruction as recited in claim 4, wherein the outlier detection algorithm in 202) is specifically as follows:
for each point, calculate its average distance to all nearby points, define the points whose average distance is outside the standard range as outliers and remove them from the data.
7. The method for measuring the volume of the large coal pile based on the monocular camera three-dimensional reconstruction as recited in claim 2, wherein the spatial transformation in the step 3) is specifically as follows:
301) acquiring normal vectors P (A, B, C) and Q (0,0,1) of an xy plane equation in a point cloud horizontal plane and a reference coordinate system;
302) set all 3D points O (x)0,y0,z0) Projecting the horizontal plane according to the following formula to obtain a new 3D point set P (x)p,yp,zp)
Figure FDA0002775450550000021
Figure FDA0002775450550000022
Figure FDA0002775450550000023
Wherein D is a fourth parameter of the horizontal plane equation;
simultaneously, the distance d between each three-dimensional space point and the equation of the horizontal plane is calculated according to the following formula
Figure FDA0002775450550000024
303) Calculating a rotation angle: the rotation angle θ between the two vectors is calculated according to the following equation
Figure FDA0002775450550000025
Wherein P and Q are two three-dimensional vectors;
304) calculating a rotating shaft: when the rotation angle is calculated, the plane where the angle is located is a plane formed by P and Q, the rotation axis must be perpendicular to the plane, and the vector before rotation is assumed to be a (a)1,a2,a3) After rotation the vector is b (b)1,b2,b3) The rotation axis c (c) is defined by the cross product1,c2,c3) Is composed of
Figure FDA0002775450550000026
305) Converting the rotation angle, the expression of the rotation axis, into a rotation matrix R according to the rodriger rotation formula;
306) set P (x) of 3D pointsp,yp,zp) And right multiplying the rotation matrix R to obtain a series of 3D point sets on the same plane, and extracting x and y values of all the 3D points to form a 2D point set p (x, y).
8. The method for measuring the volume of the large coal pile based on the monocular camera three-dimensional reconstruction as recited in claim 2, wherein the step 4) of obtaining the Delaunay triangulation network specifically comprises:
401) constructing a super triangle, including all scatter points, and putting the super triangle into a triangle linked list;
402) inserting scattered points in the point set in sequence, finding out a triangle (called as an influence triangle of the point) of which the circumscribed circle comprises an insertion point from the triangle linked list, deleting a common edge of the influence triangle, and connecting the insertion point with all vertexes of the influence triangle, thereby completing the insertion of one point in the Delaunay triangle linked list;
403) optimizing the local newly formed triangles according to an optimization criterion, and putting the formed triangles into a Delaunay triangle linked list;
404) loop through 402) until all scatter insertions have been completed.
9. The method for measuring the volume of the large coal pile based on the monocular camera three-dimensional reconstruction as recited in claim 8, wherein the set of triangles and the distance corresponding to each vertex in the triangle are obtained through the above process, the average value of the distances of the three vertices is taken as the height of the triangular prism, and the area of each triangle is calculated according to the following formula, so that the volume of all the triangular prisms can be obtained:
Figure FDA0002775450550000031
10. the method for measuring the volume of the large coal pile based on the monocular camera three-dimensional reconstruction as recited in claim 9, wherein the triangular prism volume values with partial deviation exceeding the set threshold are removed by an algorithm and a box chart method is adopted, specifically:
the 5 characteristic data values in the data are required to be solved, including the lower quartile Q1Upper quartile Q3Four-bit distance IQR ═ Q3-Q1And an upper limit Q3+1.5IQR and lower limit Q11.5IQR, eliminating the content outside the upper and lower limits, finally combining the estimated scale and the volume of all effective triangular prisms according to the following formula to estimate the actual volume value of the interested area in the object,
Figure FDA0002775450550000032
wherein K is a scale factor recovered in the point cloud processing process, djIs the length of the side edge of a triangular prism, here
Figure FDA0002775450550000033
The volume is obtained by taking the average length of the side edges of the triangular prism as the height.
CN202011263711.XA 2020-11-12 2020-11-12 Large coal pile volume measurement method based on monocular camera three-dimensional reconstruction Active CN112365534B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011263711.XA CN112365534B (en) 2020-11-12 2020-11-12 Large coal pile volume measurement method based on monocular camera three-dimensional reconstruction

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011263711.XA CN112365534B (en) 2020-11-12 2020-11-12 Large coal pile volume measurement method based on monocular camera three-dimensional reconstruction

Publications (2)

Publication Number Publication Date
CN112365534A true CN112365534A (en) 2021-02-12
CN112365534B CN112365534B (en) 2023-02-03

Family

ID=74514622

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011263711.XA Active CN112365534B (en) 2020-11-12 2020-11-12 Large coal pile volume measurement method based on monocular camera three-dimensional reconstruction

Country Status (1)

Country Link
CN (1) CN112365534B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114111567A (en) * 2021-09-26 2022-03-01 华能国际电力股份有限公司大连电厂 Method and device for measuring coal pile for raw coal bunker

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1653202A1 (en) * 2004-10-11 2006-05-03 Robert Schmidt Method and device for measurement of a volume flow of a belt conveyor
CN104778720A (en) * 2015-05-07 2015-07-15 东南大学 Rapid volume measurement method based on spatial invariant feature
CN106777604A (en) * 2016-12-02 2017-05-31 重庆市勘测院 A kind of pipe threeway construction method
CN109341626A (en) * 2018-10-24 2019-02-15 中国铁道科学研究院集团有限公司标准计量研究所 The difference calculation method of straightness calculation method, section maximum gauge and minimum diameter
CN109389635A (en) * 2018-09-11 2019-02-26 常州大学 A kind of coal yard excavation amount calculation method based on unmanned plane image sequence
CN109870106A (en) * 2019-02-14 2019-06-11 南京邮电大学 A kind of building volume measurement method based on unmanned plane picture
CN110893617A (en) * 2018-09-13 2020-03-20 深圳市优必选科技有限公司 Obstacle detection method and device and storage device
CN111429504A (en) * 2020-03-02 2020-07-17 武汉大学 Automatic material pile extraction and volume measurement method and system based on three-dimensional point cloud

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1653202A1 (en) * 2004-10-11 2006-05-03 Robert Schmidt Method and device for measurement of a volume flow of a belt conveyor
CN104778720A (en) * 2015-05-07 2015-07-15 东南大学 Rapid volume measurement method based on spatial invariant feature
CN106777604A (en) * 2016-12-02 2017-05-31 重庆市勘测院 A kind of pipe threeway construction method
CN109389635A (en) * 2018-09-11 2019-02-26 常州大学 A kind of coal yard excavation amount calculation method based on unmanned plane image sequence
CN110893617A (en) * 2018-09-13 2020-03-20 深圳市优必选科技有限公司 Obstacle detection method and device and storage device
CN109341626A (en) * 2018-10-24 2019-02-15 中国铁道科学研究院集团有限公司标准计量研究所 The difference calculation method of straightness calculation method, section maximum gauge and minimum diameter
CN109870106A (en) * 2019-02-14 2019-06-11 南京邮电大学 A kind of building volume measurement method based on unmanned plane picture
CN111429504A (en) * 2020-03-02 2020-07-17 武汉大学 Automatic material pile extraction and volume measurement method and system based on three-dimensional point cloud

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114111567A (en) * 2021-09-26 2022-03-01 华能国际电力股份有限公司大连电厂 Method and device for measuring coal pile for raw coal bunker

Also Published As

Publication number Publication date
CN112365534B (en) 2023-02-03

Similar Documents

Publication Publication Date Title
CN110717983B (en) Building elevation three-dimensional reconstruction method based on knapsack type three-dimensional laser point cloud data
CN112595258B (en) Ground object contour extraction method based on ground laser point cloud
CN111932688B (en) Indoor plane element extraction method, system and equipment based on three-dimensional point cloud
CN110340891B (en) Mechanical arm positioning and grabbing system and method based on point cloud template matching technology
CN109816664B (en) Three-dimensional point cloud segmentation method and device
Grussenmeyer et al. Comparison methods of terrestrial laser scanning, photogrammetry and tacheometry data for recording of cultural heritage buildings
CN109900338B (en) Method and device for measuring volume of pavement pit
CN107481274B (en) Robust reconstruction method of three-dimensional crop point cloud
CN111553292A (en) Rock mass structural plane identification and occurrence classification method based on point cloud data
CN104809759A (en) Large-area unstructured three-dimensional scene modeling method based on small unmanned helicopter
CN113791400B (en) Stair parameter autonomous detection method based on laser radar
CN112233249A (en) B spline surface fitting method and device based on dense point cloud
CN109919955A (en) The tunnel axis of ground formula laser radar point cloud extracts and dividing method
Altantsetseg et al. 3d surface reconstruction of stone tools by using four-directional measurement machine
CN116152306A (en) Method, device, apparatus and medium for determining masonry quality
CN113763529B (en) Substation modeling method based on three-dimensional scanning
CN112687005A (en) Coal pile volume measurement method based on three-dimensional reconstruction
WO2016207669A2 (en) A method of generating a three dimensional representation of an environment or system
CN112365534B (en) Large coal pile volume measurement method based on monocular camera three-dimensional reconstruction
CN116266365A (en) Point cloud data extraction method and device, electronic equipment and storage medium
CN112687010B (en) Digital metering method of end frame drilling jig
CN117036944B (en) Tree carbon sink amount calculating method and system based on point cloud data and image recognition
CN112161568A (en) Large coal pile volume measuring method based on multi-view common camera
CN112505723B (en) Three-dimensional map reconstruction method based on navigation point selection
Wei et al. A 3d vehicle recognition system based on point cloud library

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