CN111415415B - Automatic construction method of three-dimensional map cutting geological section - Google Patents
Automatic construction method of three-dimensional map cutting geological section Download PDFInfo
- Publication number
- CN111415415B CN111415415B CN202010242546.3A CN202010242546A CN111415415B CN 111415415 B CN111415415 B CN 111415415B CN 202010242546 A CN202010242546 A CN 202010242546A CN 111415415 B CN111415415 B CN 111415415B
- Authority
- CN
- China
- Prior art keywords
- dimensional
- point
- geological
- line
- stratum
- 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
- 238000010276 construction Methods 0.000 title claims abstract description 14
- 230000015572 biosynthetic process Effects 0.000 claims description 33
- 238000000034 method Methods 0.000 claims description 17
- 238000001341 grazing-angle X-ray diffraction Methods 0.000 claims description 12
- 238000005070 sampling Methods 0.000 claims description 11
- 230000011218 segmentation Effects 0.000 claims description 3
- 230000000694 effects Effects 0.000 abstract description 3
- 230000001351 cycling effect Effects 0.000 abstract 1
- 239000008186 active pharmaceutical agent Substances 0.000 description 8
- 239000003795 chemical substances by application Substances 0.000 description 2
- 230000019771 cognition Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- ATJFFYVFTNAWJD-UHFFFAOYSA-N Tin Chemical compound [Sn] ATJFFYVFTNAWJD-UHFFFAOYSA-N 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000001983 electron spin resonance imaging Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/05—Geographic models
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Geometry (AREA)
- Software Systems (AREA)
- Remote Sensing (AREA)
- Computer Graphics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Processing Or Creating Images (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention discloses an automatic construction method of a three-dimensional map-cut geological section, which comprises the following steps: (1) Reading the geological section line data, the geological map and the DEM of the map cut, and forming a section line set PL, a geological boundary set GeoLine, a stratum surface set GeoPolygon and a grid data set GeoDEM; (2) Any section line is obtained from the set PL, and based on the set GeoLine, the intersection points of the section line and all geological boundary lines are obtained to form a three-dimensional earth surface intersection point set PP; (3) Based on GeoPolygon, a stratum yield set AT and a stratum attribute set MK of the current section line are obtained; (4) Acquiring bottom intersection points of profile bottom side lines corresponding to the current section lines and all stratum lines based on the set GeoDEM and the set PP to form a three-dimensional bottom intersection point set EP; (5) constructing a three-dimensional cut geological section based on the PP and the set EP; (6) And (5) cycling the steps (2) - (5) to obtain all three-dimensional map-cut geological sections. The invention not only improves the three-dimensional expression effect of the cut geological section, but also has higher degree of automation.
Description
Technical Field
The invention relates to the field of three-dimensional modeling and geology, in particular to an automatic construction method of a three-dimensional map-cut geological section.
Background
The map cut geological section is a geological section which is formed by selecting a certain direction on the geological section, and drawing the geological section according to various geological and geographic elements according to a certain scale by a projection method. The two-dimensional graph cut section view is matched with the geological plane view, and the working experience and the space imagination of geological specialists are combined, so that the spreading rule and the cutting relation of the human cognition geological body and the geological structure in the three-dimensional space can be realized to a certain extent.
However, the expression mode of the geological information is too abstract, is unfavorable for information acquisition and rule cognition of geological space, and is difficult to meet the requirements of three-dimensional geological modeling and expression. The three-dimensional map cut geological section can support the three-dimensional geological information expression to a certain extent while the characteristics of the conventional two-dimensional map cut section are reserved, and can provide a data base for the three-dimensional geological modeling application based on the section. Therefore, the research of the automatic construction method of the three-dimensional map-cut geological section is developed, and the method has important practical value and research significance.
Disclosure of Invention
The invention aims to: aiming at the problems existing in the prior art, the invention provides an automatic construction method of a three-dimensional map-cut geological section, which has high automation degree and good expression effect.
The technical scheme is as follows: the automatic construction method of the three-dimensional cut geological section comprises the following steps:
(1) Reading the geological section line data, the geological map and the DEM of the map cut, and forming a section line set PL, a geological boundary set GeoLine, a stratum surface set GeoPolygon and a grid data set GeoDEM;
(2) Any section line is obtained from the section line set PL, and the intersection points of the section line and all geological boundary lines are obtained based on the geological boundary line set GeoLine to form a three-dimensional earth surface intersection point set PP;
(3) Based on a stratum plane set GeoPolygon, a stratum shape set AT and a stratum attribute set MK of the current section line are obtained;
(4) Based on a grid data set GeoDEM and a three-dimensional earth surface intersection point set PP, acquiring the intersection point of the bottom side line of the section corresponding to the current section line and the bottoms of all stratum lines to form a three-dimensional bottom intersection point set EP;
(5) Constructing a three-dimensional map cut geological section based on the three-dimensional earth surface intersection point set PP and the three-dimensional bottom intersection point set EP;
(6) And (5) circulating the steps (2) - (5) until all section lines are processed, and obtaining all three-dimensional graph cut geological sections.
Further, the step (1) specifically includes:
(1-1) reading the cut geological section line data, and storing the cut geological section line data into a section line set PL= { PL i I=1, 2, …, LI } where pl i Indicating the ith section line, LI indicating the number of section lines;
(1-2) extracting geological boundary data and stratum surface data from the geological map, and respectively storing the geological boundary set GeoLine and the stratum surface set GeoPolygon;
(1-3) extracting raster data from the DEM, and storing the raster data into a raster data set GeoDEM.
Further, the step (2) specifically includes:
(2-1) acquiring any one of the section lines pl i Reading all points on the profile generates a set of profile segmentation line endpoints tp= { TP f (x f ,y f ) |f=1, 2, …, PF }, where tp f PF is the f-th point on the cross-sectional lineThe number of table endpoints;
(2-2) for all endpoints in the set TP, the azimuth of the segment line formed by any two adjacent endpoints is calculated according to the following formula and stored in the azimuth set ag= { β s |s＝1,2,…,PF-1}：
in the formula ,(xa ,y a )、(x b ,y b ) Respectively representing coordinates of two adjacent endpoints, beta s An azimuth representing the segment line s;
(2-3) extracting a section line pl i Intersection points with all geological boundaries in the geological boundary set GeoLine, and storing stratum shapes at corresponding positions of each intersection point as attributes of the stratum shapes into an intersection point set InPoint;
(2-4) taking the intersection point set InPoint as a three-dimensional surface intersection point set PP= { PP j I j=1, 2, …, PJ }, where pp j The j-th intersection number is represented, and PJ represents the number of intersections in the three-dimensional surface intersection set PP.
Further, the step (3) specifically includes:
(3-1) obtaining stratum shapes of all points in the three-dimensional surface intersection point set PP to obtain a shape set AT= { alpha n (ρ,θ,δ)|n＝1,2,…,AN}，α n (ρ, θ, δ) represents the nth formation yield, ρ is the formation dip, θ is the formation dip, δ is the formation strike, AN represents the number of yields;
(3-2) obtaining the current section line pl i Forming a line set SeLine;
(3-3) storing the stratum unique number GID of the stratum plane set GeoPolygon as the attribute of the corresponding line in the line set Seline, and establishing the association between GeoPolygon and Seline;
(3-4) obtaining GIDs of all lines from the Seline, and storing the GIDs in the stratum attribute set MK.
Further, the step (4) specifically includes:
(4-1) at the current section line pl i At a preset distance dSampling at intervals to obtain a sampling point set IP= { IP m M=1, 2, …, PM }, where ip m The m-th point is sampled, and PM is the number of sampling points;
(4-2) according to section line pl i Is ordered by the trend of the section line pl i Is integrated with the sampling point set IP to form a ground point set
wherein ,
gp representing the kth Point k PJ is the number of intersection points in the three-dimensional earth surface intersection point set PP;
(4-3) assigning elevation values to all points in the surface point set GP according to the raster data set GeoDEM;
(4-4) calculating the section line pl in the two-dimensional coordinate system i The corresponding profile bottom side line and the ordinate of the intersection point of all stratum lines form a set ZV;
(4-5) based on the three-dimensional earth surface intersection point set PP, the azimuth angle set AG and the set ZV of the section lines, obtaining three-dimensional bottom intersection points of all stratum lines and the profile bottom side lines to form a three-dimensional bottom intersection point set
Representing the coordinates of the j-th point in the EP.
Further, the step (4-4) specifically includes:
(4-4-1) converting the set of surface points GP from three-dimensional coordinates to two-dimensional planar coordinates using the following formula:
in the formula ,
GP representing the kth point in the surface point set GP k Three-dimensional coordinates of>
Represents gp k Is a two-dimensional coordinate of (2);
(4-4-2) moving all points in the ground surface point set GP downward by a preset distance h to form a point set of bottom side lines
(4-4-3) calculating the apparent tilt angle from the set of occurrence AT using
Where ρ is the formation dip, θ is the formation dip,
azimuth for the corresponding segment line;
(4-4-5) two-dimensional plane coordinates and corresponding apparent inclination angles based on the intersection points of the ground surfaces
Calculating to obtain the bottom point coordinates of the formation line to form a bottom point set +.>
in the formula ,
represents the j-th bottom point of the formation line>
Coordinates of->
Represents the jth point pp in the surface intersection set j PJ represents the number of points in the surface intersection set, and le represents the point dp j And point pp j Is a preset distance from the first end of the first frame;
(4-4-6) Point dp based j And point pp j Construction of formation line sl j Form a stratigraphic line set sl= { SL j |j＝1,2,…,PJ}；
(4-4-7) creating a profile bottom edge cl based on the two-dimensional coordinates of the first point of the set of ground points GP, the set of points BP of the bottom edge, and the two-dimensional coordinates of the last point in the set of ground points GP;
(4-4-8) traversing the stratigraphic line set SL, obtaining the ordinate value of the intersection point of each stratigraphic line and the profile bottom side line cl, and storing the ordinate value in the set ZV= { ZV j I j=1, 2, …, PJ }.
Further, the step (4-5) specifically includes:
(4-5-1) calculating the relative elevation of each point in the surface intersection point set to a corresponding point in the three-dimensional bottom intersection point set according to the following formula:
in the formula ,hej Represents the jth point pp in the surface intersection set j Corresponding point ep j Is used for the relative elevation of the base plate,
representation pp j Z axis coordinate of (zv) j Represents the j-th value in the set ZV;
(4-5-2) calculating the three-dimensional coordinates of the three-dimensional bottom intersection point according to the following
Further, the step (5) specifically includes:
(5-1) moving the ground point set GP downward by a preset distance h according to the following formula to obtain a three-dimensional bottom point set
(5-2) constructing a MultiPatch polyhedron element based on the three-dimensional surface intersection point set PP and the three-dimensional bottom intersection point set EP, and generating a three-dimensional cut geological section set SectionMap;
(5-3) creating an ID field in the setonmap attribute table, and storing a formation attribute value corresponding to the three-dimensional formation face element in the ID field;
and (5-4) transmitting other attribute information of the GeoPolygon into an attribute table of the SectionMap according to the peer-to-peer relationship between the ID field in the SectionMap and the GID field in the GeoPolygon.
Further, the step (5-2) specifically includes:
(5-2-1) respectively inserting the first point and the last point in the ground surface point set GP into the first position and the last position in the three-dimensional ground surface intersection point set PP;
(5-2-2) inserting the first point and the last point in the three-dimensional bottom point set WP into the first and last positions in the three-dimensional bottom intersection point set EP, respectively;
(5-2-3) acquiring any two adjacent points PP in the three-dimensional surface intersection point set PP in the surface point set GP j And pp (pp) j+1 All upper points in between constitute an upper point set up= { UP l |l＝1,2,…,PL}；
(5-2-4) acquiring any two adjacent points EP in the three-dimensional bottom intersection point set EP located in the three-dimensional bottom point set WP j And ep j+1 All lower points in between constitute a lower point set xp= { XP l |l＝1,2,…,PU}；
(5-2-5) arranging UP and XP in clockwise order to form a stratum plane vertex set V= { V r |r＝1，2，…，PL+PU}；
(5-2-6) constructing MultiPatch polyhedron elements sequentially in clockwise order based on the stratum surface vertex set V to form a three-dimensional stratum surface fe w ；
(5-2-7) performing (5-2-3) - (5-2-6) circularly to obtain a three-dimensional stratum surface set FE= { FE w |w＝1，2，…，PJ+1}；
(5-2-8) tessellating each three-dimensional bedding plane FE in the three-dimensional bedding plane set FE w And constructing a single three-dimensional map cut geological section to obtain a three-dimensional map cut geological section set SectionMap.
The beneficial effects are that: according to the method, through links such as data loading, intersection point calculation, stratum occurrence and attribute acquisition, three-dimensional bottom intersection point calculation, three-dimensional graph cut section construction and the like, the three-dimensional expression effect of the graph cut geological section is improved, and the method has high automation degree.
Drawings
FIG. 1 is section line and geological map data employed in the present embodiment;
fig. 2 is DEM data employed in the present embodiment;
FIG. 3 is a flow chart of an embodiment of the present invention;
FIG. 4 is a schematic representation of a two-dimensional cross-section provided by the present invention;
FIG. 5 is a collection of three-dimensional cut geological profiles constructed in the present invention;
fig. 6 is a diagram of the relative positions of a set of three-dimensional cut geological profiles and a DEM provided in the present invention.
Detailed Description
In the following, the technical scheme of the present invention is further described in detail, and experimental data of this embodiment adopts geological data of Shang Shan in south-Beijing (fig. 1) and Shang Shan DEM data of Shang Shan in south-Beijing (fig. 2), and a projection coordinate system adopted by the experimental data is wgs_1984_utm_zone_50n. Further description will be provided by describing a specific embodiment with reference to the accompanying drawings.
As shown in fig. 3, the embodiment provides an automatic construction method for a three-dimensional cut geological section, which specifically includes the following steps:
(1) And reading the geological section line data, the geological map (containing stratum boundary lines and strata) and the DEM to form a section line set PL, a geological boundary line set GeoLine, a stratum plane set GeoPolygon and a grid data set GeoDEM.
The method specifically comprises the following steps:
(1-1) reading the cut geological section line data, and storing the cut geological section line data into a section line set PL= { PL i I=1, 2, …, LI } where pl i Indicating the ith section line, LI indicating the number of section lines; in the present embodiment, li=3;
(1-2) extracting geological boundary data and stratum surface data from the geological map, and respectively storing the geological boundary set GeoLine and the stratum surface set GeoPolygon;
(1-3) extracting raster data from the DEM, and storing the raster data into a raster data set GeoDEM.
(2) Any section line is obtained from the section line set PL, and the intersection points of the section line and all geological boundary lines are obtained based on the geological boundary line set GeoLine, so that a three-dimensional surface intersection point set PP is formed.
The method specifically comprises the following steps:
(2-1) acquiring any one of the section lines pl i Reading all points on the profile generates a set of profile segmentation line endpoints tp= { TP f (x f ，y f ) |f=1, 2, …, PF }, where tp f The f-th point on the section line, PF represents the number of endpoints;
(2-2) for all endpoints in the set TP, the azimuth of the segment line formed by any two adjacent endpoints is calculated according to the following formula and stored in the azimuth set ag= { β s |s＝1,2,…,PF-1}：
in the formula ,(xa ,y a )、(x b ,y b ) Respectively representing coordinates of two adjacent endpoints, beta s An azimuth representing the segment line s; in the present embodiment, when i=1, pf=2, x a ＝153224.48,y a ＝150759.99,x b ＝153193.79,,y b ＝150579.27,β＝3.31。
(2-3) extracting the section line pl based on Arcgis Engine API i Intersection points with all geological boundaries in the geological boundary set GeoLine, and storing stratum shapes at corresponding positions of each intersection point as attributes of the stratum shapes into an intersection point set InPoint;
(2-4) taking the intersection point set InPoint as a three-dimensional surface intersection point set PP= { PP j I j=1, 2, …, PJ }, where pp j The j-th intersection number is represented, PJ represents the number of intersections in the three-dimensional surface intersection set PP, and pj=6 when i=1.
(3) And (5) acquiring a stratum yield set AT and a stratum attribute set MK of the current section line based on the stratum surface set GeoPolygon.
The method specifically comprises the following steps:
(3-1) obtaining stratum shapes of all points in the three-dimensional surface intersection point set PP to obtain a shape set AT= { alpha n (ρ,θ，δ)|n＝1，2，…，AN}，α n (ρ, θ, δ) represents the nth formation yield, ρ is the formation dip, θ is the formation dip, δ is the formation strike, AN represents the number of yields; when i=1, an=6;
(3-2) call Arcgis Engine API, obtain the current section line pl i Forming a line set SeLine;
(3-3) storing the stratum unique code GID of the stratum surface set GeoPolygon as the attribute of the corresponding line in the line set Seline, and establishing the association between GeoPolygon and Seline;
(3-4) obtaining GIDs of all lines from the Seline, and storing the GIDs in the stratum attribute set MK.
(4) And acquiring the bottom intersection points of the profile bottom side line corresponding to the current section line and all stratum lines based on the grid data set GeoDEM and the three-dimensional earth surface intersection point set PP to form a three-dimensional bottom intersection point set EP.
The method specifically comprises the following steps:
(4-1) at the current section line pl i Sampling according to a preset distance d at intervals to obtain a sampling point set IP= { IP m M=1, 2, …, PM }, where ip m The m-th point is sampled, and PM is the number of sampling points; in this embodiment, d=4, and when i=1, pm=158;
(4-2) according to section line pl i Is ordered by the trend of the section line pl i Is integrated with the sampling point set IP to form a ground point set
wherein ,
gp representing the kth Point k PJ is the number of intersection points in the three-dimensional earth surface intersection point set PP; when i=1 in the present embodiment, pm+pj+2=165;
(4-3) loading the GeoDEM and calling Arcgis Engine API, and assigning elevation values to all points in the ground surface point set GP according to the grid data set GeoDEM;
(4-4) calculating the section line pl in the two-dimensional coordinate system i The corresponding profile bottom side line and the ordinate of the intersection point of all stratum lines form a set ZV;
(4-5) based on the three-dimensional earth surface intersection point set PP, the azimuth angle set AG and the set ZV of the section lines, obtaining three-dimensional bottom intersection points of all stratum lines and the profile bottom side lines to form a three-dimensional bottom intersection point set
Representing the coordinates of the j-th point in the EP.
Wherein, the step (4-4) specifically comprises:
(4-4-1) converting the set of surface points GP from three-dimensional coordinates to two-dimensional planar coordinates using the following formula:
in the formula ,
GP representing the kth point in the surface point set GP k Three-dimensional coordinates of>
Represents gp k Is a two-dimensional coordinate of (2);
(4-4-2) moving all points in the ground surface point set GP downward by a preset distance h to form a point set of bottom side lines
Wherein h=20;
(4-4-3) calculating the apparent tilt angle from the set of occurrence AT using
Where ρ is the formation dip, θ is the formation dip,
azimuth for the corresponding segment line;
(4-4-5) two-dimensional plane coordinates and corresponding apparent inclination angles based on the intersection points of the ground surfaces
Calculating to obtain the bottom point coordinates of the formation line to form a bottom point set +.>
/>
in the formula ,
represents the j-th bottom point dp of the formation line j Coordinates of->
Represents the jth point pp in the surface intersection set j PJ represents the number of points in the surface intersection set, and le represents the point dp j And point pp j Is 20;
(4-4-6) Point dp based j And point pp j Construction of formation line sl j Form a stratigraphic line set sl= { SL j |j＝1,2,…,PJ}；
(4-4-7) invoking Arcgis Engine API to create a profile bottom edge cl based on the two-dimensional coordinates of the first point of the set of ground points GP, the set of points BP of the bottom edge, and the two-dimensional coordinates of the last point in the set of ground points GP;
(4-4-8) traversing the stratigraphic line set SL to find each stratigraphic layerThe ordinate value of the intersection of the line and the profile bottom edge cl is saved to the set zv= { ZV j I j=1, 2, …, PJ }. In this embodiment, a schematic diagram of the ordinate of the bottom intersection point is shown in fig. 4.
Wherein, the step (4-5) specifically comprises:
(4-5-1) calculating the relative elevation of each point in the surface intersection point set to a corresponding point in the three-dimensional bottom intersection point set according to the following formula:
in the formula ,hej Represents the jth point pp in the surface intersection set j Corresponding point ep j Is used for the relative elevation of the base plate,
representation pp j Z axis coordinate of (zv) j Represents the j-th value in the set ZV;
(4-5-2) calculating the three-dimensional coordinates of the three-dimensional bottom intersection point according to the following
(5) And constructing a three-dimensional map cut geological section based on the three-dimensional surface intersection point set PP and the three-dimensional bottom intersection point set EP.
The method specifically comprises the following steps:
(5-1) moving the ground point set GP downward by a preset distance h according to the following formula to obtain a three-dimensional bottom point set
/>
(5-2) constructing a MultiPatch polyhedron element according to the three-dimensional surface intersection point set PP and the three-dimensional bottom intersection point set EP based on Arcgis Engine API, and generating a three-dimensional cut geological section set SectionMap;
(5-3) creating an ID field in the setonmap attribute table, and storing a formation attribute value corresponding to the three-dimensional formation face element in the ID field;
and (5-4) transmitting other attribute information of the GeoPolygon into an attribute table of the SectionMap according to the peer-to-peer relationship between the ID field in the SectionMap and the GID field in the GeoPolygon.
Wherein, the step (5-2) specifically comprises:
(5-2-1) respectively inserting the first point and the last point in the ground surface point set GP into the first position and the last position in the three-dimensional ground surface intersection point set PP;
(5-2-2) inserting the first point and the last point in the three-dimensional bottom point set WP into the first and last positions in the three-dimensional bottom intersection point set EP, respectively;
(5-2-3) acquiring any two adjacent points PP in the three-dimensional surface intersection point set PP in the surface point set GP j And pp (pp) j+1 All upper points in between constitute an upper point set up= { UP l L=1, 2, …, PL; in the present embodiment, when j=1, pl=21;
(5-2-4) acquiring any two adjacent points EP in the three-dimensional bottom intersection point set EP located in the three-dimensional bottom point set WP j And ep j+1 All lower points in between constitute a lower point set xp= { XP l L=1, 2, …, PU; in the present embodiment, pu=4 when j=1;
(5-2-5) arranging UP and XP in clockwise order to form a stratum plane vertex set V= { V r |r＝1,2，…，PL+PU}；
(5-2-6) constructing MultiPatch polyhedron elements sequentially in clockwise order based on the stratum surface vertex set V to form a three-dimensional stratum surface fe w ；
(5-2-7) performing (5-2-3) - (5-2-6) circularly to obtain a three-dimensional stratum surface set FE= { FE w |w＝1,2,…,PJ+1}；
(5-2-8) mosaic three-dimensional floor setCombining each three-dimensional ground plane FE in FE w And constructing a single three-dimensional map cut geological section to obtain a three-dimensional map cut geological section set SectionMap.
(6) And (5) circulating the steps (2) - (5) until all section lines are processed, and obtaining all three-dimensional graph cut geological sections. In this embodiment, a three-dimensional set of cut geological profiles is constructed as shown in FIG. 5.
In this embodiment, the relative positions of the three-dimensional cut geological profile set and the DEM are as shown in fig. 6. After the construction is completed, the ESRI TIN three-dimensional cut geological section model can be exported to be a three-dimensional model file in obj, FBX and other formats based on Arcgis Engine API.
In the embodiment of the invention, partial GIS operation is provided based on Arcgis Engine API, and related steps can also use APIs of software such as SuperMap, arcgis Object and the like to perform corresponding GIS operation.
Claims (7)
1. An automatic construction method of a three-dimensional cut geological section is characterized by comprising the following steps:
(1) Reading the geological section line data, the geological map and the DEM of the map cut, and forming a section line set PL, a geological boundary set GeoLine, a stratum surface set GeoPolygon and a grid data set GeoDEM;
(2) Any section line is obtained from the section line set PL, and the intersection points of the section line and all geological boundary lines are obtained based on the geological boundary line set GeoLine to form a three-dimensional earth surface intersection point set PP;
(3) Based on a stratum plane set GeoPolygon, a stratum shape set AT and a stratum attribute set MK of the current section line are obtained;
(4) Based on a grid data set GeoDEM and a three-dimensional earth surface intersection point set PP, acquiring the intersection point of the bottom side line of the section corresponding to the current section line and the bottoms of all stratum lines to form a three-dimensional bottom intersection point set EP; the method specifically comprises the following steps:
(4-1) at the current section line pl i Sampling according to a preset distance d at intervals to obtain a sampling point set IP= { IP m M=1, 2, …, PM }, where ip m The m-th point is sampled, and PM is the number of sampling points;
(4-2) according toSection line pl i Is ordered by the trend of the section line pl i Is integrated with the sampling point set IP to form a ground point set
wherein ,
gp representing the kth Point k PJ is the number of intersection points in the three-dimensional earth surface intersection point set PP;
(4-3) assigning elevation values to all points in the surface point set GP according to the raster data set GeoDEM;
(4-4) calculating the section line pl in the two-dimensional coordinate system i The corresponding profile bottom side line and the ordinate of the intersection point of all stratum lines form a set ZV;
(4-5) based on the three-dimensional earth surface intersection point set PP, the azimuth angle set AG and the set ZV of the section lines, obtaining three-dimensional bottom intersection points of all stratum lines and the profile bottom side lines to form a three-dimensional bottom intersection point set
Representing the coordinates of the j-th point in the EP;
(5) Constructing a three-dimensional map cut geological section based on the three-dimensional earth surface intersection point set PP and the three-dimensional bottom intersection point set EP; the method specifically comprises the following steps:
(5-1) moving the ground point set GP downward by a preset distance h according to the following formula to obtain a three-dimensional bottom point set
(5-2) constructing a MultiPatch polyhedron element based on the three-dimensional surface intersection point set PP and the three-dimensional bottom intersection point set EP, and generating a three-dimensional cut geological section set SectionMap;
(5-3) creating an ID field in the setonmap attribute table, and storing a formation attribute value corresponding to the three-dimensional formation face element in the ID field;
(5-4) transmitting other attribute information of GeoPolygon into an attribute table of the SectionMap according to the peer-to-peer relationship between the ID field in the SectionMap and the GID field in the GeoPolygon;
(6) And (5) circulating the steps (2) - (5) until all section lines are processed, and obtaining all three-dimensional graph cut geological sections.
2. The method for automatically constructing a three-dimensional cut geological section according to claim 1, wherein: the step (1) specifically comprises:
(1-1) reading the cut geological section line data, and storing the cut geological section line data into a section line set PL= { PL i I=1, 2, …, LI } where pl i Indicating the ith section line, LI indicating the number of section lines;
(1-2) extracting geological boundary data and stratum surface data from the geological map, and respectively storing the geological boundary set GeoLine and the stratum surface set GeoPolygon;
(1-3) extracting raster data from the DEM, and storing the raster data into a raster data set GeoDEM.
3. The method for automatically constructing a three-dimensional cut geological section according to claim 1, wherein: the step (2) specifically comprises:
(2-1) acquiring any one of the section lines pl i Reading all points on the profile generates a set of profile segmentation line endpoints tp= { TP f (x f ,y f ) |f=1, 2, …, PF }, where tp f The f-th point on the section line, PF represents the number of endpoints;
(2-2) for all endpoints in the set TP, the azimuth of the segment line formed by any two adjacent endpoints is calculated according to the following formula and stored in the azimuth set ag= { β s |s＝1,2,…,PF-1}：
in the formula ,(xa ,y a )、(x b ,y b ) Respectively representing coordinates of two adjacent endpoints, beta s An azimuth representing the segment line s;
(2-3) extracting a section line pl i Intersection points with all geological boundaries in the geological boundary set GeoLine, and storing stratum shapes at corresponding positions of each intersection point as attributes of the stratum shapes into an intersection point set InPoint;
(2-4) taking the intersection point set InPoint as a three-dimensional surface intersection point set PP= { PP j I j=1, 2, …, PJ }, where pp j The j-th intersection number is represented, and PJ represents the number of intersections in the three-dimensional surface intersection set PP.
4. The method for automatically constructing a three-dimensional cut geological section according to claim 1, wherein: the step (3) specifically comprises:
(3-1) obtaining stratum shapes of all points in the three-dimensional surface intersection point set PP to obtain a shape set AT= { alpha n (ρ,θ,δ)|n＝1,2,…,AN}，α n (ρ, θ, δ) represents the nth formation yield, ρ is the formation dip, θ is the formation dip, δ is the formation strike, AN represents the number of yields;
(3-2) obtaining the current section line pl i Forming a line set SeLine;
(3-3) storing the stratum unique number GID of the stratum plane set GeoPolygon as the attribute of the corresponding line in the line set Seline, and establishing the association between GeoPolygon and Seline;
(3-4) obtaining GIDs of all lines from the Seline, and storing the GIDs in the stratum attribute set MK.
5. The method for automatically constructing a three-dimensional cut geological section according to claim 1, wherein: the step (4-4) specifically comprises:
(4-4-1) converting the set of surface points GP from three-dimensional coordinates to two-dimensional planar coordinates using the following formula:
in the formula ,
GP representing the kth point in the surface point set GP k Three-dimensional coordinates of>
Represents gp k Is a two-dimensional coordinate of (2);
(4-4-2) moving all points in the ground surface point set GP downward by a preset distance h to form a point set of bottom side lines
(4-4-3) calculating the apparent tilt angle from the set of occurrence AT using
Wherein ρ is the formation dip, θ is the formation dip, β' s Azimuth for the corresponding segment line;
(4-4-5) two-dimensional plane coordinates and corresponding apparent inclination angles based on the intersection points of the ground surfaces
Calculating to obtain the bottom point coordinates of the formation line to form a bottom point set +.>
in the formula ,
represents the j-th bottom point dp of the formation line j Coordinates of->
Represents the jth point pp in the surface intersection set j PJ represents the number of points in the surface intersection set, and le represents the point dp j And point pp j Is a preset distance from the first end of the first frame;
(4-4-6) Point dp based j And point pp j Construction of formation line sl j Form a stratigraphic line set sl= { SL j |j＝1,2,…,PJ}；
(4-4-7) creating a profile bottom edge cl based on the two-dimensional coordinates of the first point of the set of ground points GP, the set of points BP of the bottom edge, and the two-dimensional coordinates of the last point in the set of ground points GP;
(4-4-8) traversing the stratigraphic line set SL, obtaining the ordinate value of the intersection point of each stratigraphic line and the profile bottom side line cl, and storing the ordinate value in the set ZV= { ZV j I j=1, 2, …, PJ }.
6. The method for automatically constructing a three-dimensional cut geological section according to claim 1, wherein: the step (4-5) specifically comprises:
(4-5-1) calculating the relative elevation of each point in the surface intersection point set to a corresponding point in the three-dimensional bottom intersection point set according to the following formula:
in the formula ,hej Represents the jth point pp in the surface intersection set j Corresponding point ep j Is used for the relative elevation of the base plate,
representation pp j Z axis coordinate of (zv) j Represents the j-th value in the set ZV;
(4-5-2) calculating the three-dimensional coordinates of the three-dimensional bottom intersection point according to the following
7. The method for automatically constructing a three-dimensional cut geological section according to claim 1, wherein: the step (5-2) specifically comprises:
(5-2-1) respectively inserting the first point and the last point in the ground surface point set GP into the first position and the last position in the three-dimensional ground surface intersection point set PP;
(5-2-2) inserting the first point and the last point in the three-dimensional bottom point set WP into the first and last positions in the three-dimensional bottom intersection point set EP, respectively;
(5-2-3) acquiring any two adjacent points PP in the three-dimensional surface intersection point set PP in the surface point set GP j And pp (pp) j+1 All upper points in between constitute an upper point set up= { UP l |l＝1,2,…,PL}；
(5-2-4) acquiring any adjacent ones of the set of three-dimensional bottom points WP located in the set of three-dimensional bottom intersections EPTwo points ep j And ep j+1 All lower points in between constitute a lower point set xp= { XP l |l＝1,2,…,PU}；
(5-2-5) arranging UP and XP in clockwise order to form a stratum plane vertex set V= { V r |r＝1,2,…,PL+PU}；
(5-2-6) constructing MultiPatch polyhedron elements sequentially in clockwise order based on the stratum surface vertex set V to form a three-dimensional stratum surface fe w ；
(5-2-7) performing (5-2-3) - (5-2-6) circularly to obtain a three-dimensional stratum surface set FE= { FE w |w＝1,2,…,PJ+1}；
(5-2-8) tessellating each three-dimensional bedding plane FE in the three-dimensional bedding plane set FE w And constructing a single three-dimensional map cut geological section to obtain a three-dimensional map cut geological section set SectionMap.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010242546.3A CN111415415B (en) | 2020-03-31 | 2020-03-31 | Automatic construction method of three-dimensional map cutting geological section |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010242546.3A CN111415415B (en) | 2020-03-31 | 2020-03-31 | Automatic construction method of three-dimensional map cutting geological section |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111415415A CN111415415A (en) | 2020-07-14 |
| CN111415415B true CN111415415B (en) | 2023-05-09 |
Family
ID=71493482
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010242546.3A Active CN111415415B (en) | 2020-03-31 | 2020-03-31 | Automatic construction method of three-dimensional map cutting geological section |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111415415B (en) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111968231B (en) * | 2020-08-14 | 2023-05-30 | 成都理工大学 | Three-dimensional stratum modeling method based on geological map cutting section |
| CN112233230B (en) * | 2020-09-08 | 2024-02-27 | 南京师范大学 | Three-dimensional model construction method and device for fault structure in cut geological section |
| CN112233238B (en) * | 2020-10-28 | 2024-02-27 | 南京师范大学 | Bedrock geologic body three-dimensional model construction method and device based on graph cut geological parallel section |
| CN112734926B (en) * | 2021-01-29 | 2024-04-19 | 南京师范大学 | Automatic generation method of cut geological section for loose layer coverage area |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105469443B (en) * | 2014-09-30 | 2018-06-26 | 中国地质调查局发展研究中心 | Based on geotraverse（PRB）The method that process Data fuse generates three-dimensional geological figure |
| CN106934858B (en) * | 2017-03-14 | 2018-01-16 | 中国地质科学院矿产资源研究所 | Three-dimensional geological modeling method and system for scale region of mining area |
| CN109753707B (en) * | 2018-12-25 | 2023-10-24 | 核工业北京地质研究院 | Method for extracting stratum boundary line by utilizing exploration line profile to develop three-dimensional modeling |
-
2020
- 2020-03-31 CN CN202010242546.3A patent/CN111415415B/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN111415415A (en) | 2020-07-14 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111415415B (en) | Automatic construction method of three-dimensional map cutting geological section | |
| CN104574511B (en) | A kind of quick progressive three-dimensional geological modeling method | |
| CN103646423B (en) | A kind of method of three-dimensional geological modeling and device | |
| CN111079217B (en) | BIM-based geotechnical engineering comprehensive investigation information interpretation method and system | |
| CN110189409B (en) | PLAXIS-based rapid true three-dimensional geological modeling method and system | |
| CN102609982B (en) | Topology discovery method of space geological data based on unstructured mode | |
| CN109003330A (en) | A kind of three dimensional contour line method based on basement rock boundary constraint | |
| CN102385067B (en) | Drawing method for isoline containing reverse fault | |
| CN109472866B (en) | A kind of modeling method of big geology Building Information Model | |
| CN111553007B (en) | Method for automatically generating geometric working conditions in slope two-dimensional calculation section | |
| CN116152461B (en) | Geological modeling method, device, computer equipment and computer readable storage medium | |
| CN112150582B (en) | Multi-modal data-oriented geological profile approximate expression method | |
| CN110032771A (en) | A kind of accurate method of cutting out of open coal mine DEM for taking local detail feature into account | |
| CN112464479A (en) | Three-dimensional ore formation prediction method and system based on GOCAD | |
| CN109741449A (en) | A kind of geologic section setting method based on three dimensional space coordinate transformation | |
| CN112084548B (en) | 2.5D visual modeling method and system for slope stability analysis | |
| CN106777117A (en) | A kind of automatic identifying method of horizontal stratum tectonic landform | |
| CN109285221A (en) | Deposit three-dimensional Geological Modeling | |
| CN102903149A (en) | Block forming method and device for geologic model | |
| CN110942510A (en) | Three-dimensional rapid combination method for planar geological map | |
| CN106971422B (en) | Method for compiling stability analysis model of arch dam abutment based on CATIA | |
| CN115035258A (en) | Efficient urban three-dimensional geological modeling method based on CAD (computer-aided design) drilling histogram | |
| CN111429573B (en) | Automatic construction method of three-dimensional geological drilling model | |
| Chorowicz et al. | Dip and strike measured systematically on digitized three-dimensional geological maps | |
| CN114463494A (en) | Automatic topographic feature line extracting algorithm |
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 |