CN115311132A - Underground cavern flattening method, geological logging method, electronic device and storage medium - Google Patents
Underground cavern flattening method, geological logging method, electronic device and storage medium Download PDFInfo
- Publication number
- CN115311132A CN115311132A CN202210949035.4A CN202210949035A CN115311132A CN 115311132 A CN115311132 A CN 115311132A CN 202210949035 A CN202210949035 A CN 202210949035A CN 115311132 A CN115311132 A CN 115311132A
- Authority
- CN
- China
- Prior art keywords
- point
- cavern
- point cloud
- area
- left wall
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 64
- 238000004590 computer program Methods 0.000 claims description 7
- 230000014509 gene expression Effects 0.000 claims description 4
- 230000006870 function Effects 0.000 claims description 3
- 238000002372 labelling Methods 0.000 claims description 2
- 238000010586 diagram Methods 0.000 description 4
- 238000004422 calculation algorithm Methods 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000013499 data model Methods 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000000638 solvent extraction Methods 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T3/00—Geometric image transformations in the plane of the image
- G06T3/06—Topological mapping of higher dimensional structures onto lower dimensional surfaces
- G06T3/067—Reshaping or unfolding 3D tree structures onto 2D planes
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T3/00—Geometric image transformations in the plane of the image
- G06T3/14—Transformations for image registration, e.g. adjusting or mapping for alignment of images
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Processing Or Creating Images (AREA)
Abstract
The invention discloses an underground cavern flattening method, a geological logging method, electronic equipment and a storage medium, wherein the flattening method comprises the steps of obtaining an original cavern point cloud model and coordinates of points of the original cavern point cloud model; selecting a control point, and dividing an original cavern point cloud model into a crown region, a left wall region and a right wall region according to the control point; rotating the original cavern point cloud model according to the rotation angle; the left wall area is kept unchanged, the top arch area and the right wall area are flattened to the surface where the left wall area is located, and the left wall area, the flattened top arch area and the right wall area are integrated to obtain a point cloud data set of the flattened cavern point cloud model; and reversely rotating the flattened cavern point cloud model according to the rotation angle to ensure that the unchanged area of the flattened cavern point cloud model is completely overlapped with the corresponding area of the original cavern point cloud model, so as to obtain a point cloud data set of the reversely rotated cavern point cloud model. The method can solve the problems of low recording efficiency and the like caused by point cloud shielding of the left wall and the right wall in the model.
Description
Technical Field
The invention belongs to the technical field of point cloud data processing, and particularly relates to an underground cavern flattening method based on three-dimensional point cloud data, a geological logging method, electronic equipment and a storage medium.
Background
The three-dimensional laser scanning technology has the advantages of high modeling speed, high spatial position precision and the like, so that the three-dimensional point cloud data model is more and more widely applied to the geological recording process of hydropower engineering. On the basis of a three-dimensional point cloud model of the cavern, geological elements in the cavern can be accurately and quickly drawn, and geometric parameters and occurrence information of the geological elements are automatically calculated, so that the traditional manual geological logging method is greatly simplified, the problems of low logging efficiency and accuracy are solved, and the safety of geological operators is guaranteed.
According to the geological survey technical standard of the Chinese geological survey bureau, the underground cavern recording technology mainly adopts a two-wall one-top expansion method, namely, the underground cavern is divided into a left wall, a top arch and a right wall, and geological elements are drawn on the relative position of coordinate paper according to the proportion along the central axis of the cavern, so that a geological recording map is obtained. The existing geological logging method based on the three-dimensional cavern model has the problems that logging efficiency is low due to mutual shielding of the left wall and the right wall, geological elements drawn on the three-dimensional point cloud model are not uniform with a traditional two-dimensional geological logging map, and the like.
Disclosure of Invention
The invention aims to provide an underground cavern flattening method, a geological logging method, electronic equipment and a storage medium, and aims to solve the problems that the left wall and the right wall are mutually shielded and the logging efficiency is low in the traditional method, and the problem that the logging result on a three-dimensional cavern model is not matched with a traditional two-dimensional geological logging map.
The invention solves the technical problems through the following technical scheme: a method for flattening an underground cavern comprises the following steps:
acquiring an original cavern point cloud model and coordinates of each point in the original cavern point cloud model on a coordinate system;
selecting a control point on the original cavern point cloud model, and dividing the original cavern point cloud model into a crown area, a left wall area and a right wall area according to the control point;
calculating a rotation angle, and rotating the original cavern point cloud model according to the rotation angle to enable an opening extension line to be parallel to an X axis of the coordinate system;
based on the rotated cavern point cloud model, flattening the top arch area and the right wall area to the surface of the left wall area, and integrating the left wall area, the flattened top arch area and the right wall area to obtain a point cloud data set of the flattened cavern point cloud model under the coordinate system; or the right wall area is kept unchanged, the top arch area and the left wall area are flattened to the surface where the right wall area is located, and the right wall area, the flattened top arch area and the left wall area are integrated to obtain a point cloud data set of the flattened cavern point cloud model under the coordinate system;
and reversely rotating the flattened cavern point cloud model on the XY plane of the coordinate system according to the rotation angle to ensure that the unchanged area of the flattened cavern point cloud model is completely overlapped with the corresponding area of the original cavern point cloud model, thereby obtaining a point cloud data set of the cavern point cloud model after reverse rotation in the coordinate system.
Further, the coordinate system is a geodetic coordinate system or a custom three-dimensional coordinate system.
Further, the control points comprise a left wall upper point, a right wall upper point, a top arch middle point and a right wall lower point; the concrete implementation process of dividing the original cavern point cloud model into a crown region, a left wall region and a right wall region according to the control point is as follows:
dividing the original cavern point cloud model into an arch crown area and other areas according to the left wall upper point and the right wall upper point;
the remaining regions are divided into left and right wall regions according to the crown midpoint.
Further, the specific calculation of the rotation angle is:
calculating a rotation angle according to the opening extension line vector and the unit vector on the X axis of the coordinate system, wherein the opening extension line vector points from the upper point of the left wall to the upper point of the right wall;
when P is ry -P ly >0, the rotation angle θ is:
when P is ry -P ly <0, the rotation angle θ is:
when P is present ry -P ly =0, the rotation angle θ is 0;
wherein, P lx 、P ly Respectively represent a point P on the left wall l X, y values, P, on a coordinate system rx 、P ry Respectively represent points P on the right wall r X, y values in a coordinate system, cos -1 Representing an inverse cosine function.
Further, the specific implementation process of the rotation is as follows:
and rotating the original cavern point cloud model on an XOY plane of the coordinate system to enable the cavern extending line to be parallel to an X axis, so as to obtain a point cloud data set of the rotated cavern point cloud model in the coordinate system.
Further, the specific implementation process of flattening the crown region to the left wall region is as follows:
calculating a fitting curve equation of the section of the dome area after rotation in the XOZ plane of the coordinate system, which specifically comprises the following steps:
wherein x is 2 、z 2 Is the x, z value, P 'of a point on a fitted curve on the coordinate system' lx 、P′ lz Respectively represent a rotated left wall upper point P' l The x and z values on the coordinate system, a is the minor axis radius of the fitting curve, and b is the major axis radius of the fitting curve;
substituting the coordinates of the rotated left wall upper point, the rotated right wall upper point and the rotated top arch midpoint into the fitting curve equation to obtain values a and b;
dividing the rotated vault area into different sections along the central axis of the cavern point cloud model, and using the section fitting curve expressions corresponding to the different sections to obtain the fitting curve equation and the plane equation y 2 ＝P c1y Is formed of P c1y Is a point P on the cross section c1 A y value on the coordinate system;
calculating the point P on each section c1 And dot P' c1 Of geometric distance d, its middle point P' c1 Fitting the center point O and point P of the curve to the corresponding cross section c1 The intersection point of the connecting line and the section fitting curve, and the center point O of the section fitting curve is (P' lx +a,P c1y ,P′ lz )；
Calculating the left starting point to the point P 'according to the calculus principle' c1 The length of the arc section between the two points, wherein the left starting point is the starting point on the section fitting curve on the same side of the left wall area;
determining point P 'according to the arc segment length' c1 Coordinates flattened to left wall region, according to Point P' c1 Determining point P by coordinates and distance d flattened to left wall area c1 Flattening to the coordinates of the left wall region, realizing point P c1 And flattening, so that each point on each section is flattened.
Further, the specific implementation process of flattening the right wall region to the left wall region or flattening the left wall region to the right wall region is as follows:
rotating the right wall area by 180 degrees to ensure that the upper point of the right wall area is superposed with the upper point of the right wall of the flattened top arch area to obtain the flattened right wall area;
or rotating the left wall area by 180 degrees to ensure that the upper point of the left wall area is superposed with the upper point of the left wall of the flattened arch area to obtain the flattened left wall area.
Based on the same invention concept, the invention also provides a geological logging method of the underground cavern, which comprises the following steps:
finding out a geological structure surface and exposed trace positions on the point cloud data set of the cavern point cloud model obtained by the underground cavern flattening method after the reverse rotation in the coordinate system, and selecting a point cloud data set V capable of representing the space position of geological elements;
performing three-dimensional plane fitting on all points of the point cloud data set V by using a least square method to obtain a three-dimensional fitting plane;
projecting the point cloud data set V onto a three-dimensional fitting plane, and connecting projection points into a broken line according to the trend of line elements for the line elements; for the surface elements, acquiring a minimum circumscribed convex polygon of the projection points;
and projecting the line elements and the plane elements onto a projection plane by taking the YOZ plane of the coordinate system as the projection plane, labeling the information of each group of line elements and plane elements, and exporting a document in a DXF format to obtain the geological record sketch.
Preferably, the point cloud data set V capable of representing the spatial position of the geological element is selected in a way of selecting a surface in a box or drawing a line.
Based on the same inventive concept, the invention further provides an electronic device, which comprises a memory and a processor, wherein the memory stores a computer program capable of running on the processor, and the processor executes the steps of the underground cavern flattening method or the underground cavern geological recording method when running the computer program.
Based on the same inventive concept, the present invention also provides a computer-readable storage medium, which is a non-volatile storage medium or a non-transitory storage medium, having stored thereon a computer program, which when executed by a processor, performs the steps of the underground cavern flattening method or the underground cavern geological logging method as described above.
Advantageous effects
Compared with the prior art, the invention has the advantages that:
according to the underground cavern flattening method, the geological logging method, the electronic equipment and the storage medium, the top arch area and the right wall area are rotated to the vertical elevation in the left wall elevation direction, or the top arch area and the left wall area are rotated to the vertical elevation in the right wall elevation direction, so that the problems of overlapping and shielding of point clouds of the left wall and the right wall in the logging process of the three-dimensional cavern model are solved, and operators can conveniently and accurately identify and draw geological elements; the cavern point cloud model flattening process is also a geological element flattening process, a two-dimensional geological record draft can be automatically generated after geological elements are drawn on the cavern point cloud model after flattening, and the problem that the cavern point cloud model is not matched with a traditional two-dimensional geological record draft after being recorded on a three-dimensional cavern model is solved.
The invention firstly rotates and then flattens, ensures that the central axis of the cavern is parallel to the y axis of a coordinate system, and reduces the sequential flattening process of different sections from a three-dimensional space to a two-dimensional space, thereby greatly simplifying the algorithm complexity.
Drawings
In order to more clearly illustrate the technical solution of the present invention, the drawings needed to be used in the description of the embodiments are briefly introduced below, and it is obvious that the drawings in the following description are only one embodiment of the present invention, and it is obvious for those skilled in the art to obtain other drawings based on the drawings without creative efforts.
FIG. 1 is a flow chart of a method for flattening a subterranean cavity in accordance with an embodiment of the present invention;
FIG. 2 is a schematic diagram illustrating selection of an original cavern point cloud model and control points according to an embodiment of the invention;
FIG. 3 is a schematic diagram of flattening a section of an arch region of an original cavern point cloud model in an embodiment of the invention;
FIG. 4 is a schematic diagram illustrating a comparison effect between a flattened cavern point cloud model and an original cavern point cloud model according to an embodiment of the invention;
fig. 5 is a two-dimensional geological record map automatically rendered in an embodiment of the present invention.
Detailed Description
The technical solutions in the present invention are 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 only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be obtained by a person skilled in the art without inventive efforts based on the embodiments of the present invention, shall fall within the scope of protection of the present invention.
The technical means of the present application will be described in detail with specific examples. The following several specific embodiments may be combined with each other, and details of the same or similar concepts or processes may not be repeated in some embodiments.
The embodiment of the invention describes a method for flattening the underground cavern based on three-dimensional live-action model data, the whole model comprises 1697879 point clouds, and the model is subjected to coordinate registration by using a total station, so that the real coordinate of each point cloud in the model in a geodetic coordinate system can be obtained. As shown in FIG. 1, the method for flattening the underground cavern comprises the following steps:
1. and acquiring coordinates of each point in the original cavern point cloud model and the original cavern point cloud model on a coordinate system.
In this embodiment, the coordinate system is a geodetic coordinate system or a custom three-dimensional coordinate system. The self-defined three-dimensional coordinate system is a random three-dimensional coordinate system self-defined in the space where the original cavern point cloud model is located, and exemplarily, the three-dimensional coordinate system is constructed by taking the upper point of the left wall of the original cavern point cloud model as the original point, the upper point of the left wall pointing to the upper point of the right wall as the positive direction of the X axis, the central axis of the cavern pointing from the opening to the inside of the cavern as the direction of the Y axis, and the vertical direction of the vertical plane of the left wall pointing to the top of the cavern as the Z axis.
After a coordinate system is determined (the coordinate system is used as a reference in subsequent steps), the coordinates of each point in the original cavern point cloud model in the coordinate system can be determined, and a point cloud data set of the original cavern point cloud model in the coordinate system is obtained and is marked as S = { p = 1 ,p 2 ,p 3 ,…,p i ,…,p n N is the number of points in the point cloud dataset S, p i Is the ith point in the point cloud data set S.
2. Region partitioning
And selecting a control point on the original cavern point cloud model. In the present embodiment, as shown in FIG. 2, the control points include a point P on the left wall l Right wall upper point P r Middle point of top archP c And the right wall lower point P d . Obtaining basic parameters of the three-dimensional cavern model according to the left wall upper point, the right wall upper point, the top arch middle point and the right wall lower point, wherein the basic parameters comprise wall height, top height, tunnel width and ground elevation; the wall height is equal to the vertical height difference between the upper point of the right wall and the lower point of the right wall, the top height is equal to the vertical height difference between the middle point of the top arch and the upper point of the left wall, the hole width is equal to the horizontal straight-line distance between the upper point of the left wall and the upper point of the right wall, and the ground elevation is equal to the vertical height of the lower point of the right wall. Illustratively, through four control points, the wall height of the original cavern point cloud model can be calculated to be 4.06m, the crown arch to be 7.38m, the cavern width to be 7.68m, the cavern length to be 15.59m, and the ground elevation to be 413.71m.
Dividing an original cavern point cloud model into a crown region, a left wall region and a right wall region according to the control points, and specifically realizing the following steps:
step 2.1: dividing the original cavern point cloud model into a crown region and other regions according to the left wall upper point and the right wall upper point;
step 2.2: the remaining areas are divided into left and right wall areas according to the crown midpoint.
3. Rotation of original cavern point cloud model
Because the original cavern point cloud model has three-dimensional space coordinate values, the complexity of the algorithm is increased by directly fitting a section curve to the arch region. In the flattening process of the original cavern point cloud model, the central axis of the cavern is taken as a symmetry axis, the original cavern point cloud model is rotated and then flattened, the central axis of the cavern can be ensured to be parallel to the Y axis of a coordinate system, the flattening process of different sections is reduced from a three-dimensional space to a two-dimensional plane in sequence, and therefore algorithm complexity is greatly simplified.
In this embodiment, the rotation angle is calculated first, and the original cavern point cloud model is rotated according to the rotation angle, and the specific implementation process is as follows:
step 3.1: and calculating to obtain the rotation angle according to the opening extension line vector and the unit vector on the X axis of the coordinate system, wherein the opening extension line vector points from the upper point of the left wall to the upper point of the right wall.
When P is present ry -P ly >0, the rotation angle is:
when P is present ry -P ly <0, the rotation angle is:
when P is present ry -P ly And =0, the rotation angle is θ =0.
Wherein, P lx 、P ly Respectively represent a point P on the left wall l X, y values, P, on a coordinate system rx 、P ry Respectively represent points P on the right wall r X, y values in a coordinate system, cos -1 Representing an inverse cosine function.
Step 3.2: rotating the original cavern point cloud model on an XOY surface of a coordinate system to enable an opening extension line to be parallel to an X axis, enabling a central axis of the cavern point cloud model to be parallel to a Y axis of the coordinate system, and obtaining a point cloud data set of the rotated cavern point cloud model in the coordinate system, wherein the point cloud data set is marked as S '= { p' 1 ,p′ 2 ,p′ 3 ,…,p′ i ,…,p′ n }, the point cloud data set S 'comprises a left wall region S' l And right wall region S' r And dome region S' c 。
Rotation matrix R of original cavern point cloud model θ Comprises the following steps:
4. area flattening
Based on the rotated cavern point cloud model, under the condition that the left wall area is kept unchanged, the top arch area and the right wall area are flattened to the surface where the left wall area is located, or under the condition that the right wall area is kept unchanged, the top arch area and the left wall area are flattened to the surface where the right wall area is located. In this embodiment, as shown in fig. 3, a specific implementation process of flattening is described by taking a surface where the top arch region and the right wall region are flattened to the left wall region as an example:
step 4.1: calculating a fitting curve equation of the section of the dome area after rotation in the XOZ plane of the coordinate system, which specifically comprises the following steps:
wherein x is 2 、z 2 Is the x, z value, P 'of the point on the fitted curve on the coordinate system' lx 、P′ lz Respectively represent a rotated left wall upper point P' l The x and z values on the coordinate system, wherein a is the minor axis radius of the fitting curve, and b is the major axis radius of the fitting curve;
step 4.2: and substituting the coordinates of the rotated left wall upper point, the rotated right wall upper point and the center point of the top arch into a fitting curve equation to obtain the values of a and b.
Step 4.3: dividing the rotated vault area into different sections along the central axis of the cavern point cloud model, and using the section fitting curve expressions corresponding to the different sections to obtain the fitting curve equation and the plane equation y 2 ＝P c1y The composition is as follows:
wherein, P c1y Is a point P on the cross section c1 A y value on the coordinate system;
step 4.4: calculating the point P on each section c1 And dot P' c1 Of geometric distance d, its middle point P' c1 Fitting the center point O and the point P of the curve to the corresponding section c1 The intersection point of the connecting line and the section fitting curve, and the center point O of the section fitting curve is (P' lx +a,P c1y ,P′ lz ) (ii) a If point P c1 And d is negative above the section fitting curve, and is positive if not.
Step 4.5: calculating the left starting point to the point P 'according to the calculus principle' c1 The length of the arc segment between, i, wherein the left starting point refers toThe starting point on the section fitting curve on the same side with the left wall area and the specific calculation formula of the arc section length are as follows:
wherein the content of the first and second substances,
is a left starting point, a center point O and a point P' c1 The angle of the formed sector.
Step 4.6: determining point P 'according to arc segment length l' c1 Flattened to the coordinates of the left wall region, point P' c1r Of point P' c1r Is (P' lx ，P c1y ,P′ lz +l)。
Step 4.7: according to point P' c1 Coordinates flattened to left wall region (i.e., point P' c1r Coordinates of) and distance d determines point P c1 Coordinates flattened to the left wall region (i.e., point P) c1r Coordinates of) point P), point P c1r Is (P' lx +d，P c1y ,P′ lz + l), the realization point P c1 And flattening the surface of the left wall area.
Step 4.8: traversing each point on each section of the crown region, repeating the steps of 4.3-4.6, flattening the crown region to the surface of the left arm region, and obtaining a point cloud data set S' of the flattened crown region c 。
The specific implementation process of flattening the right wall region to the left wall region or flattening the left wall region to the right wall region is as follows:
rotating the right wall area by 180 degrees to ensure that the upper point of the right wall area coincides with the upper point of the right wall of the flattened crown area to obtain the flattened right wall area and obtain a point cloud data set S' of the flattened right wall area r ；
Or rotating the left wall area by 180 degrees to ensure that the upper point of the left wall area coincides with the upper point of the left wall of the flattened arch area to obtain the flattened left wall area and obtain the point cloud data set of the flattened left wall areaS″ l 。
The left wall region S' l Flattened crown region S ″) c And a right wall region S ″ r Integrating to obtain a point cloud data set S' of the flattened cavern point cloud model in a coordinate system; or the right wall region S' r Flattened crown region S ″) c And a left wall region S ″ l Integrating to obtain a point cloud data set S '= { p' of the flattened cavern point cloud model in a coordinate system 1 ,p″ 2 ,p″ 3 ,…,p″ i ,…,p″ n }。
5. Reverse reversal
And reversely rotating the flattened cavern point cloud model in an XY plane of a coordinate system according to the rotation angle to ensure that the unchanged area of the flattened cavern point cloud model is completely overlapped with the corresponding area of the original cavern point cloud model, thereby obtaining a point cloud data set S' of the cavern point cloud model in the coordinate system after reverse rotation.
Illustratively, under the condition that the left wall area is kept unchanged, the crown arch area and the right wall area are flattened to the surface where the left wall area is located, and then during reverse rotation, the flattened cavern point cloud model is reversely rotated on an XO Y plane of a coordinate system according to the rotation angle, so that the left wall area is completely overlapped with the left wall area of the original cavern point cloud model, as shown in FIG. 4.
Based on the same invention concept, the embodiment of the invention also provides a geological logging method of the underground cavern, which comprises the following steps:
1. finding out the position of a geological structure surface and an exposed trace on the point cloud data set S' obtained by the underground cavern flattening method, and selecting the point cloud data set V capable of representing the space position of the geological element in a surface selection or line drawing mode.
2. And (3) performing three-dimensional plane fitting on all points of the point cloud data set V by using a least square method, wherein the general expression of the three-dimensional plane fitting is as follows:
Ax 3 +By 3 +Cz 3 +D＝0
a, B, C is not 0 at the same time, and n = { A, B, C } forms the normal vector of the plane, D is a constant term, x 3 、y 2 、z 3 And (4) fitting x, y and z values of points on the plane on the coordinate system.
3. Projecting the point cloud data set V onto a three-dimensional fitting plane, and connecting projection points into a broken line according to the trend of line elements for the line elements; for surface elements, the minimum circumscribed convex polygon of the projection points is obtained.
4. The YOZ plane of the coordinate system is used as a projection plane, the line elements and the plane elements are projected onto the projection plane, information such as names of all groups of line elements and plane elements is marked, and a DXF format file is exported, so that a geological record sketch (shown in figure 5) can be obtained, and vector diagram support is provided for operators to carry out geological record drawings.
The above disclosure is only for the specific embodiments of the present invention, but the scope of the present invention is not limited thereto, and any person skilled in the art can easily conceive of changes or modifications within the technical scope of the present invention, and shall be covered by the scope of the present invention.
Claims (10)
1. A method for flattening an underground cavern is characterized by comprising the following steps:
acquiring an original cavern point cloud model and coordinates of each point in the original cavern point cloud model on a coordinate system;
selecting a control point on the original cavern point cloud model, and dividing the original cavern point cloud model into a top arch area, a left wall area and a right wall area according to the control point;
calculating a rotation angle, and rotating the original cavern point cloud model according to the rotation angle to enable an opening extension line to be parallel to an X axis of the coordinate system;
based on the rotated cavern point cloud model, flattening the top arch area and the right wall area to the surface of the left wall area, and integrating the left wall area, the flattened top arch area and the right wall area to obtain a point cloud data set of the flattened cavern point cloud model under the coordinate system;
or the right wall area is kept unchanged, the top arch area and the left wall area are flattened to the surface where the right wall area is located, and the right wall area, the flattened top arch area and the left wall area are integrated to obtain a point cloud data set of the flattened cavern point cloud model under the coordinate system;
and reversely rotating the flattened cavern point cloud model on the XY plane of the coordinate system according to the rotation angle to ensure that the unchanged area of the flattened cavern point cloud model is completely overlapped with the corresponding area of the original cavern point cloud model, thereby obtaining a point cloud data set of the cavern point cloud model after reverse rotation in the coordinate system.
2. The underground cavern flattening method of claim 1, wherein the coordinate system is a geodetic coordinate system or a custom three-dimensional coordinate system.
3. The underground cavern flattening method of claim 1, wherein the control points comprise a left wall upper point, a right wall upper point, a roof arch midpoint and a right wall lower point; the concrete implementation process of dividing the original cavern point cloud model into a crown region, a left wall region and a right wall region according to the control point is as follows:
dividing the original cavern point cloud model into an item arch area and other areas according to the left wall upper point and the right wall upper point;
and dividing the rest area into a left wall area and a right wall area according to the arch crown midpoint.
4. The underground cavern flattening method of claim 1, wherein the rotation angle is specifically calculated as follows:
calculating a rotation angle according to the opening extension line vector and a unit vector on the X axis of the coordinate system, wherein the opening extension line vector points from a left wall upper point to a right wall upper point;
when P is present ry -P ly If > 0, the rotation angle theta is as follows:
when P is present ry -P ly If < 0, the rotation angle theta is:
when P is ry -P ly =0, the rotation angle θ is 0;
wherein, P lx 、P ly Respectively represent a point P on the left wall l X, y values, P, on a coordinate system rx 、P ry Respectively represent points P on the right wall r X, y values, cos, on a coordinate system -1 Representing an inverse cosine function.
5. The underground cavern flattening method of claim 1, wherein the rotation is realized by the following specific processes:
and rotating the original cavern point cloud model on an XOY plane of the coordinate system to enable the cavern extending line to be parallel to an X axis, so as to obtain a point cloud data set of the rotated cavern point cloud model in the coordinate system.
6. A method for flattening a subterranean cavern as claimed in any one of claims 1 to 5, wherein the specific implementation process of flattening the crown region to the left wall region is as follows:
calculating a fitting curve equation of the section of the dome area after rotation in the XOZ plane of the coordinate system, which specifically comprises the following steps:
wherein x is 2 、z 2 Is the x, z value, P 'of a point on a fitted curve on the coordinate system' lx 、P′ lz Respectively represent a rotated left wall upper point P' l The x and z values on the coordinate system, a is the minor axis radius of the fitting curve, and b is the major axis radius of the fitting curve;
substituting the coordinates of the rotated left wall upper point, the rotated right wall upper point and the rotated top arch midpoint into the fitting curve equation to obtain values a and b;
dividing the rotated vault area into different sections along the central axis of the cavern point cloud model, and using the section fitting curve expressions corresponding to the different sections as the fitting curve equation and the plane equation y 2 ＝P c1y Is formed of P c1y Is a point P on the cross section c1 A y value on the coordinate system;
calculating the point P on each section c1 And dot P' c1 A geometric distance d therebetween, wherein point P' c1 Fitting the center point O and point P of the curve to the corresponding cross section c1 The intersection point of the connecting line and the section fitting curve, and the center point O of the section fitting curve is (P' lx +a，P c1y ，P′ lz )；
Calculating the left starting point to the point P 'according to the calculus principle' c1 The length of the arc section between the two points, wherein the left starting point is the starting point on the section fitting curve on the same side of the left wall area;
determining point P 'according to the arc segment length' c1 Coordinates flattened to left wall region, according to Point P' c1 Determining point P by flattening coordinates and distance d to left wall region c1 Flattening to the coordinates of the left wall region, realizing point P c1 And flattening, so that each point on each section is flattened.
7. A method for flattening a subterranean cavern as recited in any one of claims 1 to 5, wherein the right wall region is flattened to the left wall region or the left wall region is flattened to the right wall region by:
rotating the right wall area by 180 degrees to ensure that the upper point of the right wall area is superposed with the upper point of the right wall of the flattened top arch area to obtain the flattened right wall area;
or rotating the left wall area by 180 degrees to ensure that the upper point of the left wall area is superposed with the upper point of the left wall of the flattened arch area to obtain the flattened left wall area.
8. A geological logging method for underground caverns is characterized by comprising the following steps:
finding out a geological structure surface and an exposed trace position on a point cloud data set under a coordinate system by using the reversely rotated cavern point cloud model obtained by the underground cavern flattening method according to any one of claims 1 to 7, and selecting a point cloud data set V capable of representing the space position of geological elements;
performing three-dimensional plane fitting on all points of the point cloud data set V by using a least square method to obtain a three-dimensional fitting plane;
projecting the point cloud data set V onto a three-dimensional fitting plane, and connecting projection points into a broken line according to the trend of line elements for the line elements; for the surface elements, acquiring a minimum circumscribed convex polygon of a projection point;
projecting the line elements and the plane elements onto a projection plane by taking a YOZ plane of a coordinate system as the projection plane, labeling information of each group of line elements and plane elements, and exporting a document in a DXF format to obtain a geological record sketch;
preferably, the point cloud data set V capable of representing the spatial position of the geological element is selected in a way of selecting a surface in a box or drawing a line.
9. An electronic device comprising a memory and a processor, the memory having stored thereon a computer program operable on the processor, wherein: the processor, when executing the computer program, performs the steps of the method for flattening a subterranean cavern as defined in any one of claims 1 to 7 or the method for geological logging of a subterranean cavern as defined in claim 8.
10. A computer-readable storage medium, the computer-readable storage medium being a non-volatile storage medium or a non-transitory storage medium having a computer program stored thereon, characterized in that: the computer program when executed by a processor performs the steps of the subterranean cavern flattening method as defined in any one of claims 1 to 7 or the subterranean cavern geological logging method as defined in claim 8.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210949035.4A CN115311132A (en) | 2022-08-09 | 2022-08-09 | Underground cavern flattening method, geological logging method, electronic device and storage medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210949035.4A CN115311132A (en) | 2022-08-09 | 2022-08-09 | Underground cavern flattening method, geological logging method, electronic device and storage medium |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN115311132A true CN115311132A (en) | 2022-11-08 |
Family
ID=83860133
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210949035.4A Pending CN115311132A (en) | 2022-08-09 | 2022-08-09 | Underground cavern flattening method, geological logging method, electronic device and storage medium |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115311132A (en) |
-
2022
- 2022-08-09 CN CN202210949035.4A patent/CN115311132A/en active Pending
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103761397B (en) | Three-dimensional model slice for surface exposure additive forming and projection plane generating method | |
| KR100916474B1 (en) | Methods of partitioning a region represented by contours into smaller polygonal zones and calculating data for digital elevation model and data for constructing geographical surfaces of the region, and a recording medium in which programs for implementing the methods are recorded | |
| CN110489778B (en) | Graph segmentation method and laser etching control system for laser etching processing | |
| US8994726B1 (en) | Systems and methods for preparing a model of an environment for display | |
| JPH05342310A (en) | Method and device for three-dimensional conversion of linear element data | |
| JPH0550030B2 (en) | ||
| CN108665530B (en) | Three-dimensional modeling implementation method based on single picture | |
| CN104361632A (en) | Triangular mesh hole-filling method based on Hermite radial basis function | |
| CN107146287B (en) | Two-dimensional projection image to threedimensional model mapping method | |
| CN110032771A (en) | A kind of accurate method of cutting out of open coal mine DEM for taking local detail feature into account | |
| CN107729707A (en) | Engineering construction setting out method based on mobile augmented reality technology and BIM | |
| CN113591300A (en) | 3D printing file generation method and device, computer equipment and storage medium | |
| CN113421338A (en) | Method for three-dimensional rendering of two-dimensional vector | |
| CN110132168A (en) | A kind of three-dimensional laser point cloud data processing system | |
| CN116822160A (en) | Cartesian grid generation method, device, equipment and medium | |
| CN115311132A (en) | Underground cavern flattening method, geological logging method, electronic device and storage medium | |
| CN103617650A (en) | Displaying method for complex three-dimensional terrain | |
| CN115035233B (en) | Global multi-resolution space grid spherical ray projection visualization method and device | |
| CN108986212B (en) | Three-dimensional virtual terrain LOD model generation method based on crack elimination | |
| CN113420460A (en) | Urban building height limit rapid analysis method and system based on OSG data skyline | |
| CN118097339B (en) | Deep learning sample enhancement method and device based on low-altitude photogrammetry | |
| CN110335353B (en) | Method, device, equipment and storage medium for automatically subdividing blocks | |
| JPH07113928B2 (en) | Figure layout display method | |
| CN114329575B (en) | Contour line encryption method, system, medium, and device based on interpolated ridge line | |
| CN115311426A (en) | Underground cavern automatic flattening method, recording method, electronic equipment and storage medium |
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 |