Geologic body three-dimensional model top and bottom plate extraction method and electronic equipment

Abstract

The invention discloses a geologic body three-dimensional model top and bottom plate extraction method and electronic equipment. The method comprises the steps of generating a geologic body three-dimensional model top and bottom plate boundary contour line set based on coordinate data of each contour point; determining intersection point coordinates of the top plate and the bottom plate of the geologic body three-dimensional model based on the space intersection state of any two edges in the first triangular mesh set, the number of intersection points of any two edges projected to the two-dimensional plane and the intersection point coordinates; generating a vertex coordinate set of the top and bottom plates of the geologic body three-dimensional model based on the vertex coordinates of the top and bottom plates; based on each original triangulation network topological structure line, constructing an original triangulation network topological structure line constraint set of the top and bottom plates of the geologic body three-dimensional model; and performing triangulation based on the intersection point coordinate set, the vertex coordinate set, the geologic body three-dimensional model roof and floor boundary contour line set and the original triangulation network topological structure line constraint set to generate the geologic body three-dimensional model roof and floor. Therefore, the top plate and the bottom plate can be finely and efficiently extracted.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of three-dimensional modeling, in particular to a method for extracting a top / bottom plate of a three-dimensional geological model and an electronic device. BACKGROUND

[0002] A three-dimensional geological model is a geological model that visually expresses the form and structure of a geological body in the ground through digitization. In the three-dimensional geological model, the top plate or bottom plate model of the three-dimensional geological model plays a crucial role in precise mining design, quantitative safety evaluation, and efficient production management in mines.

[0003] Currently, the extraction of the top plate or bottom plate model is usually performed in a human-computer interaction manner, that is, in a three-dimensional model display interface, the triangular facets belonging to the top plate or bottom plate on the surface of the geological body are manually selected and merged to construct the corresponding model. Although this method can meet the actual needs to some extent, it has the following disadvantages: (1) the human-computer interaction method is time-consuming and laborious, and the operation efficiency is low, and after the dynamic update of the three-dimensional geological model, the top plate or bottom plate model cannot be updated in time; (2) the human-computer interaction method is difficult to handle the top / bottom plate model at the junction of multiple geological body models, which is prone to logical errors of spatial intersection at the junction of multiple geological body models.

[0004] To overcome the above problems, the prior art such as the method for extracting a top / bottom plate of a three-dimensional geological model based on a meshing principle, device, equipment, medium and product discloses a method for extracting a top / bottom plate of a three-dimensional model based on a meshing principle, which can effectively solve the above problems, but still has certain limitations: the method reconstructs the triangular mesh topology of the top / bottom plate, which leads to a certain deviation between the extracted top / bottom plate and the original three-dimensional model. SUMMARY

[0005] Therefore, the embodiments of the present application provide a method for extracting a top / bottom plate of a three-dimensional geological model and an electronic device, aiming to solve the problems of low efficiency in the extraction process of the top plate or bottom plate model, lagging model update, logical errors at the junction, and inconsistent topological structure in the prior art, so as to realize more refined and consistent top / bottom plate extraction.

[0006] The technical scheme of the embodiments of the present application is as follows:

[0007] In a first aspect, the embodiments of the present application provide a method for extracting a top / bottom plate of a three-dimensional geological model, which comprises:

[0008] obtaining a polygon set of a first triangular mesh set of a plurality of geologic body three-dimensional models after two-dimensional plane projection, a polygon in the polygon set comprising at least one of the following types: an outer contour, an island, and a hole; emitting a ray along a set direction for each polygon vertex of the polygon set, determining a first intersection point set between the polygon vertex and the first triangular mesh set; generating coordinate data of a contour point based on a highest-elevation first intersection point in the first intersection point set and two-dimensional coordinate data of the polygon vertex; and generating a top / bottom boundary contour line set of the geologic body three-dimensional model based on coordinate data of each of the contour points.

[0009] determining intersection point coordinates of a top / bottom plate of the geologic body three-dimensional model based on a spatial intersection state of any two edges in the first triangular mesh set, a number of intersection points of the any two edges projected onto a two-dimensional plane, and intersection point coordinates; and determining an intersection point coordinate set of the top / bottom plate of the geologic body three-dimensional model based on each of the intersection point coordinates.

[0010] for each triangular mesh vertex in the first triangular mesh set, emitting a ray along a set direction, and determining a vertex coordinate of a top / bottom plate of the geologic body three-dimensional model if the triangular mesh vertex does not have a ray intersection point with other triangular meshes in the first triangular mesh set; and generating a vertex coordinate set of the top / bottom plate of the geologic body three-dimensional model based on each of the vertex coordinates of the top / bottom plate.

[0011] traversing each triangular edge in the first triangular mesh set, determining that the triangular edge is a primitive triangular net topology line of the top / bottom plate of the geologic body three-dimensional model if both vertices of the triangular edge are in the vertex coordinate set and no intersection point in the intersection point coordinate set is on the triangular edge; and constructing a primitive triangular net topology line constraint set of the top / bottom plate of the geologic body three-dimensional model based on each of the primitive triangular net topology lines.

[0012] triangulating based on the intersection point coordinate set, the vertex coordinate set, the top / bottom boundary contour line set of the geologic body three-dimensional model, and the primitive triangular net topology line constraint set of the top / bottom plate of the geologic body three-dimensional model, to generate a top / bottom plate of the geologic body three-dimensional model.

[0013] In some embodiments, the method further comprises:

[0014] obtaining a second triangular mesh set of the first triangular mesh set of the plurality of geologic body three-dimensional models after two-dimensional plane projection, the second triangular mesh set comprising two-dimensional coordinate data of a plurality of triangles;

[0015] Performing a Boolean operation on the two-dimensional coordinate data of each second triangular mesh in the second set of triangular meshes to generate the set of polygons, the set of polygons including at least one of the following: coordinate data of an outer contour, coordinate data of an island, and coordinate data of a hole, the outer contour and the island being counterclockwise polygons, and the hole being a clockwise polygon.

[0016] In some embodiments, the intersection point coordinates of the top and bottom plates of the geological body three-dimensional model are determined based on a spatial intersection state of any two edges in the first set of triangular meshes, a number of intersection points and intersection point coordinates of the any two edges projected onto a two-dimensional plane, the spatial intersection state of the any two edges including non-intersection, and the non-intersection including:

[0017] If it is determined that the spatial intersection state of the any two edges is non-intersection and the number of intersection points of the any two edges projected onto a two-dimensional plane is 0, the any two edges are ignored.

[0018] In some embodiments, the intersection point coordinates of the top and bottom plates of the geological body three-dimensional model are determined based on a spatial intersection state of any two edges in the first set of triangular meshes, a number of intersection points and intersection point coordinates of the any two edges projected onto a two-dimensional plane, the spatial intersection state of the any two edges including intersection, and the intersection including:

[0019] If it is determined that the spatial intersection state of the any two edges is intersection and the number of intersection points of the any two edges projected onto a two-dimensional plane is one, the intersection point is set as a second intersection point, a ray is emitted in a set direction based on the second intersection point, and a second set of intersection points between the second intersection point and the first set of triangular meshes is obtained; if an elevation value of a third intersection point with the highest elevation value in the second set of intersection points is equal to an elevation value of the first intersection point, coordinates of the third intersection point are determined as the intersection point coordinates of the top and bottom plates of the geological body three-dimensional model.

[0020] In some embodiments, the intersection point coordinates of the top and bottom plates of the geological body three-dimensional model are determined based on a spatial intersection state of any two edges in the first set of triangular meshes, a number of intersection points and intersection point coordinates of the any two edges projected onto a two-dimensional plane, the spatial intersection state of the any two edges including non-intersection, and the non-intersection including:

[0021] If it is determined that the spatial intersection state of the two edges is not intersected, and the two edges projected to a two-dimensional plane have two intersection points, the two intersection points are set as a fourth intersection point and a fifth intersection point respectively; rays are emitted to the fourth intersection point and the fifth intersection point respectively, and a third intersection point set between the fourth intersection point, the fifth intersection point and the first triangular mesh set is obtained; if the elevation value of a sixth intersection point with the highest elevation value in the third intersection point set is the same as the elevation value of the fourth intersection point or the fifth intersection point, the coordinates of the sixth intersection point are determined as the intersection point coordinates of the top and bottom plates of the geological body three-dimensional model.

[0022] In some embodiments, the method further comprises:

[0023] If the triangular mesh vertex has a ray intersection point with other triangular meshes in the first triangular mesh set, the triangular mesh vertex is ignored.

[0024] In some embodiments, the method further comprises:

[0025] Each triangular edge in the first triangular mesh set is traversed, and if it is determined that at least one of the two vertices of the triangular edge is not in the vertex coordinate set, the triangular edge is ignored.

[0026] In some embodiments, the method further comprises:

[0027] Each triangular edge in the first triangular mesh set is traversed, and if it is determined that both vertices of the triangular edge are in the vertex coordinate set, and it is determined that there is an intersection point in the intersection point coordinate set on the triangular edge, the triangular edge is ignored.

[0028] In some embodiments, the triangular mesh is generated based on the intersection point coordinate set, the vertex coordinate set, the boundary contour line set of the top and bottom plates of the geological body three-dimensional model, and the original triangular mesh topological structure line constraint set of the top and bottom plates of the geological body three-dimensional model, including:

[0029] The line-constrained Delaunay triangulation is performed based on the intersection point coordinate set, the vertex coordinate set, the boundary contour line set of the top and bottom plates of the geological body three-dimensional model, and the original triangular mesh topological structure line constraint set of the top and bottom plates of the geological body three-dimensional model, to generate a set of triangulation meshes;

[0030] The other triangular meshes in the set of triangulation meshes that are outside the outer contour of the boundary contour line set of the top and bottom plates of the geological body three-dimensional model are deleted, and the top and bottom plates of the geological body three-dimensional model are generated based on the set of triangulation meshes after deletion.

[0031] In a second aspect, an electronic device is provided, which includes a processor and a memory for storing a computer program capable of running on the processor, wherein the processor is configured to run the computer program to perform the steps of the method in the first aspect of the embodiments.

[0032] The technical scheme provided by the embodiments of the present application, the method comprises: obtaining a polygon set after the first triangular mesh set of a plurality of geologic body three-dimensional models is projected on a two-dimensional plane, and the polygons in the polygon set include at least one of the following types: an outer contour, an island, and a hole; a ray is emitted in a set direction for each polygon vertex of the polygon set, and a first intersection point set between the polygon vertex and the first triangular mesh set is determined; based on the two-dimensional coordinate data of the polygon vertex and the first intersection point with the highest elevation value in the first intersection point set, coordinate data of a contour point is generated; based on the coordinate data of each contour point, a set of top / bottom plate boundary contour lines of the geologic body three-dimensional model is generated; based on the spatial intersection state of any two edges in the first triangular mesh set, the number of intersection points and the intersection point coordinates of the any two edges projected on a two-dimensional plane, the intersection point coordinates of the top / bottom plate of the geologic body three-dimensional model are determined; based on each intersection point coordinate, a set of intersection point coordinates of the top / bottom plate of the geologic body three-dimensional model is determined; for each triangular mesh vertex in the first triangular mesh set, a ray is emitted in a set direction, and if the triangular mesh vertex does not have a ray intersection point with other triangular meshes in the first triangular mesh set, the triangular mesh vertex is determined as a vertex coordinate of the top / bottom plate of the geologic body three-dimensional model; based on each vertex coordinate of the top / bottom plate, a set of vertex coordinates of the top / bottom plate of the geologic body three-dimensional model is generated; each triangular edge in the first triangular mesh set is traversed, and if it is determined that both vertices of the triangular edge are in the set of vertex coordinates and that there is no intersection point on the triangular edge in the set of intersection point coordinates, the triangular edge is determined as a primitive triangular net topology line of the top / bottom plate of the geologic body three-dimensional model; based on each primitive triangular net topology line, a primitive triangular net topology line constraint set of the top / bottom plate of the geologic body three-dimensional model is constructed; based on the set of intersection point coordinates, the set of vertex coordinates, the set of top / bottom plate boundary contour lines of the geologic body three-dimensional model, and the primitive triangular net topology line constraint set of the top / bottom plate of the geologic body three-dimensional model, a triangular mesh is generated for the top / bottom plate of the geologic body three-dimensional model.

[0033] Thus, this application can achieve refined, convenient, and efficient extraction of the top and bottom plates of a three-dimensional geological model, including: (1) solving the problem of time-consuming, laborious, and inefficient human-computer interactive extraction of the top and bottom plates, while ensuring that the top or bottom plate model can be updated in a timely manner after the dynamic update of the three-dimensional geological model. Specifically, the embodiment of this application uses triangular mesh data as input and automatically executes processes such as two-dimensional projection, Boolean operation, intersection calculation, and contour extraction to generate the boundary information of the top and bottom plates without manual intervention; at the same time, since the processing flow is entirely based on the basic mesh structure, it supports rapid recalculation and incremental updates after the geological model is modified. (2) It supports the extraction of the top and bottom plate models at the junction of multiple geological models, while avoiding logical errors of spatial intersection of the top and bottom plate models at the junction of multiple geological models. Specifically, the embodiment of this application judges the intersection relationship between the triangular sides in space and two-dimensional projection, and combines the maximum elevation value to accurately identify the effective intersection points in the junction area, avoiding geometric conflicts such as spatial overlap and interpenetration in advance, and ensuring that the generated top and bottom plate models maintain logical and geometric consistency in the junction area. (3) The extracted top and bottom plates can maintain the same triangular mesh topology as the original 3D model to the greatest extent, thus achieving the goal of refined extraction of the top and bottom plate 3D model. Specifically, in the topology restoration stage of this application, a constraint set is constructed based on the original triangle edge relationship, and line constraint Delaunay triangulation is performed in combination with boundary contours, intersections, vertices and other information, so that the generated top and bottom plate structure can inherit the mesh topology features of the original model to the greatest extent. Attached Figure Description

[0034] Figure 1 A flowchart illustrating the method for extracting the top and bottom plates of a three-dimensional geological body model provided in this application embodiment;

[0035] Figure 2 A schematic diagram of a three-dimensional model of a geological body provided in an embodiment of this application;

[0036] Figure 3 A detailed structural schematic diagram of the three-dimensional model of the geological body provided in the embodiments of this application;

[0037] Figure 4 A schematic diagram of the projection result of the three-dimensional model of the geological body on a two-dimensional plane provided in the embodiments of this application;

[0038] Figure 5 This is a schematic diagram of the top boundary outline of the three-dimensional geological model provided in the embodiments of this application;

[0039] Figure 6 A schematic diagram of the top plate of the three-dimensional geological model provided in the embodiments of this application;

[0040] Figure 7A flowchart of a method for extracting a roof of a geological body three-dimensional model is shown as an application example of the present application.

[0041] Figure 8 A structural diagram of an electronic device is shown as an embodiment of the present application. DETAILED DESCRIPTION

[0042] The present application will be further described in detail below with reference to the accompanying drawings and embodiments.

[0043] Unless otherwise defined, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this application belongs. The terminology used in the description herein is for describing particular embodiments only and is not intended to be limiting of the present application.

[0044] An embodiment of the present application provides a method for extracting a roof and floor of a geological body three-dimensional model, as shown in Figure 1 The method comprises the following steps:

[0045] Step 110: Obtain a polygon set after a first triangular mesh set of a plurality of geological body three-dimensional models is projected on a two-dimensional plane, and the polygons in the polygon set include at least one of the following types: an outer contour, an island, and a hole; shoot a ray along a set direction for each polygon vertex of the polygon set, determine a first intersection set between the polygon vertex and the first triangular mesh set; generate coordinate data of a contour point based on the highest elevation value of the first intersection in the first intersection set and the two-dimensional coordinate data of the polygon vertex; and generate a set of geological body three-dimensional model roof and floor boundary contour lines based on the coordinate data of each contour point.

[0046] In the present embodiment, the geological body three-dimensional model is a mathematical model used to express the three-dimensional morphology of the underground geological structure, and is usually constructed based on drilling data, seismic data, and other geological data, and is widely used in the fields of geological exploration, resource assessment, and environmental impact analysis.

[0047] As shown in Figure 2 , Figure 2 A model diagram of a geological body three-dimensional model is shown. By observing the geological body three-dimensional model, it can be seen that the geological body includes different layered structures and irregular surface morphologies, so that the three-dimensional morphology and internal structural features of the geological body can be intuitively understood.

[0048] It can be understood that the triangular mesh is the basic unit that constitutes the geological body three-dimensional model, and therefore, the geological body three-dimensional model can be divided into a series of triangular meshes to approximately represent its shape, and each triangular vertex has a specific three-dimensional coordinate, which defines the position and morphology of the geological body surface. In the present application, the geological body three-dimensional model includes a first triangular mesh set composed of a plurality of first triangular meshes.

[0049] Exemplarily, a plurality of geological body three-dimensional models are composed of n triangular meshes, which are set as a set T = {t1, t2,..., t n}, i.e., a first triangular mesh set. As shown in Figure 3 , Figure 3 a detailed structure of a geological body three-dimensional model is shown. In Figure 3 , it can be seen that the geological body is subdivided into many small triangular meshes, which together constitute the surface and internal structure of the geological body. Each triangular mesh has three vertices, and each vertex has a specific coordinate position in three-dimensional space, which together define the shape and position of the geological body.

[0050] It can be understood that two-dimensional plane projection is a geometric transformation process for projecting points or figures in three-dimensional space onto a two-dimensional plane. In this process, three-dimensional coordinates are converted into two-dimensional coordinates for analysis and processing on the plane.

[0051] In this embodiment, the polygon set obtained after the first triangular mesh set of a plurality of geological body three-dimensional models is projected onto a two-dimensional plane is obtained, and the polygons in the polygon set include at least one of the following types: an outer contour, an island, and a hole.

[0052] Here, the outer contour is the outermost boundary of the geological body on the two-dimensional plane, usually represented as a counterclockwise polygon. It defines the range and shape of the entire geological body on the plane. The hole is a gap area inside the outer contour, represented as a clockwise polygon. The existence of the hole indicates that there is a discontinuous part inside the geological body, which may be caused by geological structure or other factors. The island is an independent area inside the outer contour but not connected to the outer contour, also represented as a counterclockwise polygon. The existence of the island may represent some special structure or feature inside the geological body.

[0053] Exemplarily, the first triangular mesh set T = {t1, t2,..., t n} corresponding to a plurality of geological body three-dimensional models is obtained, and the m polygon sets G' = {g1', g2',..., g m '} corresponding thereto can be obtained, wherein the outer contour or "island" in G' is defined as a counterclockwise polygon, and the "hole" in G' is a clockwise polygon.

[0054] As shown in Figure 4 , Figure 4 a geological body three-dimensional model and its projection result on a two-dimensional plane are shown. By Figure 4In this way, the two-dimensional projection processing manner can simplify the spatial analysis of the geological body structure, while retaining the key geometric boundary information, and facilitate subsequent processing.

[0055] In this embodiment, a ray is emitted in a set direction for each polygon vertex of the polygon set, and a first intersection set between the polygon vertex and the first triangular mesh set is determined; based on the highest elevation value of the first intersection in the first intersection set and the two-dimensional coordinate data of the polygon vertex, coordinate data of an outline point is generated; and based on the coordinate data of each outline point, a set of top and bottom plate boundary contour lines of the geological body three-dimensional model is generated.

[0056] For example, assuming that the polygon set is a set G' = {g1', g2',..., g m m polygons, for each polygon vertex v in the polygon set G', a ray is made upward, the intersection with the first triangular mesh set T is calculated, and based on the highest elevation value of the first intersection in the first intersection set and the two-dimensional coordinate data of the vertex, the coordinate data of an outline point is generated; and based on the coordinate data of each outline point, a set of top and bottom plate boundary contour lines of the geological body three-dimensional model is generated.

[0057] For example, the elevation of the highest point in the intersection is assigned to v, and a set of top and bottom plate boundary contour lines G = {g1, g2,..., g m m} of the geological body three-dimensional model is obtained. As shown in FIG. 6, taking the top plate as an example, the black lines represent the finally generated top plate boundary contour lines, which accurately outline the top plate or bottom plate boundary of the geological body model in the three-dimensional space, and provide a basis for subsequent construction of the model subdivision structure and boundary constraints. Figure 5

[0058] Step 120: determining the intersection point coordinates of the top and bottom plates of the geological body three-dimensional model based on the spatial intersection state of any two edges in the first triangular mesh set, the number of intersection points and the intersection point coordinates of the projection of any two edges to a two-dimensional plane; and determining a set of intersection point coordinates of the top and bottom plates of the geological body three-dimensional model based on the intersection point coordinates.

[0059] In this embodiment, in order to accurately extract the intersection points of the top and bottom plates of the geological body three-dimensional model, the intersection relationship of any two edges in the first triangular mesh set in the three-dimensional space and the geometric relationship under the two-dimensional plane projection need to be comprehensively analyzed.

[0060] ​Specifically, the so-called "spatial intersection state" refers to whether two triangular mesh edges have an actual intersection point in three-dimensional space, i.e., whether they actually geometrically penetrate in space. The two-dimensional projection intersection point refers to whether the two edges have an intersection point on the plane after being projected onto the XY plane (two-dimensional plane), as well as the number and position of the intersection points.

[0061] Illustratively, assuming that the first triangular mesh set is set T, all edges in set T are combined two by two, the projection intersection point of each pair of edges on the XY plane is first calculated, and the number of two-dimensional projection intersection points of each pair of edges is counted. Then, in combination with the relative positions of the two edges in three-dimensional space, it is determined whether there is an actual intersection point. Only in the case of satisfying certain spatial relationship and projection geometric conditions, the intersection point coordinates of the top and bottom plates of the geological body three-dimensional model are further determined; based on the intersection point coordinates, the intersection point coordinate set of the top and bottom plates of the geological body three-dimensional model is determined.

[0062] Illustratively, the intersection point coordinate set of the top and bottom plates of the geological body three-dimensional model can be defined as set P I .

[0063] It can be understood that after the above steps 110 and 120 are completed, the boundary contour line set for defining the boundary of the top and bottom plates and the intersection point coordinate set for describing the key positions inside the boundary have been obtained. The former determines the peripheral closed shape of the top and bottom plates, and the latter supplements the geometric constraint points inside the boundary.

[0064] Step 130: For each triangular mesh vertex in the first triangular mesh set, a ray is made in a set direction, and if the triangular mesh vertex does not have a ray intersection point with other triangular meshes in the first triangular mesh set, the triangular mesh vertex is determined as a vertex coordinate of the top and bottom plates of the geological body three-dimensional model; based on the vertex coordinates of the top and bottom plates, a vertex coordinate set of the top and bottom plates of the geological body three-dimensional model is generated.

[0065] In this embodiment, in order to identify the original structure vertices directly exposed on the surface of the top and bottom plates in the geological body three-dimensional model, it is necessary to further analyze the visibility or boundary of each triangular mesh vertex in three-dimensional space. The core idea is to determine whether each vertex is blocked by other triangular meshes through ray detection, so as to filter out the explicit vertices that are actually located on the surface of the top and bottom plates.

[0066] Specifically, assuming that the first triangular mesh set is set T, which contains all triangular meshes constituting the geological body three-dimensional model. Define the vertex coordinate set P I for storing all coordinates finally determined as the vertex of the top and bottom plates.

[0067] For each vertex of each triangular mesh in the set T, a ray is emitted from the vertex with a defined direction (e.g. "up" or positive Z-axis direction), and the ray is checked for intersection with other triangular meshes in the set T.

[0068] If the ray does not intersect with any other triangular mesh after leaving the emitting vertex, it means that the vertex is a "bare vertex" on the top / bottom plate exposed surface, i.e. the vertex is not covered by the structure above and can be considered as an actual geometric vertex of the top / bottom plate. At this time, the vertex is added to the vertex coordinate set P I .

[0069] Exemplarily, the vertex coordinate set P R of the top / bottom plate of the geological body three-dimensional model is defined. T For all vertices in T, a ray is emitted upward, and p T is taken as an example to calculate whether there is an intersection point with T. If there is no intersection point, p R is added to the set P i .

[0070] Step 140: Each triangular edge in the first triangular mesh set is traversed, and if it is determined that both vertices of the triangular edge are in the vertex coordinate set and that there is no intersection point on the triangular edge in the intersection point coordinate set, it is determined that the triangular edge is an original triangular net topology line of the top / bottom plate of the geological body three-dimensional model; based on the original triangular net topology lines, an original triangular net topology line constraint set of the top / bottom plate of the geological body three-dimensional model is constructed.

[0071] In this embodiment, after the screening of the vertex coordinate set and the intersection point coordinate set in the above steps 120 and 130, in order to further restore the local topological structure relationship of the top / bottom plate, it is necessary to identify the boundary line segments that meet the conditions from the original triangular mesh structure, i.e. the so-called "original triangular net topology line". These structure lines are not only an important part of the surface morphology of the top / bottom plate, but also provide necessary boundary constraint conditions for subsequent Delaunay triangulation.

[0072] It can be understood that the original triangular net topology line constraint set of the top / bottom plate of the geological body three-dimensional model retains the topological connection information directly related to the top / bottom plate in the original model, which is used to constrain the subsequent triangulation process.

[0073] Exemplarily,

[0074] (1) defining the original triangular net topology line constraint set C of the top / bottom plate of the geological body three-dimensional model;

[0075] (2) traversing the edges of the triangular meshes in the first triangular mesh set T in turn, taking edge e i as an example, to determine whether both vertices of e i are in the set P RIn the case that all are in set P R In the case that all are in set P

[0076] (3) Determine whether there is a point on edge e I In the case that all are in set P i In the case that all are in set P i Add edge e

[0077] Finally, all edges e i that meet the above conditions are collected into set C to form a set of original triangular net topology structure lines of the top and bottom plates of the three-dimensional geological model, which is used for constrained triangulation calculation in subsequent steps to ensure that the reconstructed top and bottom plate structures are consistent with the original topology.

[0078] Step 150: Perform triangulation based on the intersection point coordinate set, the vertex coordinate set, the set of top and bottom plate boundary contour lines of the three-dimensional geological model, and the set of original triangular net topology structure lines of the top and bottom plates of the three-dimensional geological model to generate the top and bottom plates of the three-dimensional geological model.

[0079] In the present embodiment, the extraction of the related geometry and topology features of the top and bottom plates of the geological body has been completed in the foregoing steps 110 to 140, including boundary contour line extraction, intersection point calculation, vertex identification, and original topology structure line extraction. Therefore, to construct the final geological body top and bottom plate grid structure that can be used for modeling and display, the present step takes these generated data as input to perform a constrained triangulation operation.

[0080] It can be understood that the intersection point coordinate set (derived from step 120, used to express edge intersection features); the vertex coordinate set (derived from step 130, representing exposed vertices of the model surface); the set of top and bottom plate boundary contour lines of the three-dimensional geological model (derived from step 110, used to define contour constraints); and the set of original triangular net topology structure lines of the top and bottom plates of the three-dimensional geological model (derived from step 140, used to maintain the original topology structure connection).

[0081] Based on the above input sets, line-constrained Delaunay triangulation is performed on the two-dimensional plane. This algorithm generates a triangular grid that meets the quality requirements while preserving the boundary line segments, ensuring the geometric and topological continuity of the top and bottom plate structures.

[0082] Exemplarily, the line-constrained Delaunay triangulation is performed on the union of P I , P R , G, and C to obtain T U , and T Utriangular meshes outside the outer contour of G, to obtain a three-dimensional model of the top and bottom plate of the geological body, as shown in Figure 6 Figure 6 is a schematic diagram of the top plate of the three-dimensional model of the geological body.

[0083] Thus, the present application can realize fine, convenient and efficient extraction of the top and bottom plate of the three-dimensional model of the geological body, including: (1) solving the problem of time-consuming and laborious, low efficiency of human-computer interaction extraction of the top and bottom plate, while meeting the requirement of timely updating of the top plate or bottom plate model after dynamic updating of the three-dimensional model of the geological body. Specifically, the triangular mesh data is taken as input, and the boundary information of the top and bottom plate is generated without manual intervention through automatic execution of processes such as two-dimensional projection, Boolean operation, intersection calculation and contour extraction; at the same time, since the processing flow is completely based on the basic grid structure, the rapid recalculation and incremental updating after modification of the geological body model are supported. (2) Supporting the extraction of the top and bottom plate model at the junction of multiple geological body models, while avoiding logical errors of spatial intersection of the top and bottom plate model at the junction of multiple geological body models. Specifically, the present application embodiment accurately identifies the effective intersection point of the junction area by judging the intersection relationship between the triangular edges in space and two-dimensional projection, and combining the maximum elevation value, to avoid geometric conflicts such as spatial overlap and insertion in advance, and ensure that the generated top and bottom plate model maintains logical and geometric consistency in the junction area. (3) The extracted top and bottom plate can maintain the same triangular net topology structure as the original three-dimensional model to the greatest extent, achieving the purpose of fine extraction of the top and bottom plate three-dimensional model. Specifically, in the topology structure restoration stage, the present application embodiment constructs a constraint set according to the original triangular edge relationship, and executes line-constrained Delaunay triangulation combining the boundary contour, intersection point, vertex and other information, so that the generated top and bottom plate structure can inherit the grid topology features of the original model to the greatest extent.

[0084] In some embodiments, the method further comprises:

[0085] obtaining a second triangular mesh set projected on a two-dimensional plane from a first triangular mesh set of a plurality of three-dimensional models of geological bodies, the second triangular mesh set comprising two-dimensional coordinate data of a plurality of triangles;

[0086] performing Boolean operation on the two-dimensional coordinate data of each second triangular mesh in the second triangular mesh set to generate a polygon set, the polygon set comprising at least one of the following: coordinate data of an outer contour, coordinate data of an island and coordinate data of a hole, the outer contour and the island being counterclockwise polygons, and the hole being a clockwise polygon.

[0087] ​In the present embodiment, the construction of the polygon set is the basis for subsequent roof-floor boundary extraction. First, the first triangular mesh set of the plurality of geological body three-dimensional models needs to be projected in two dimensions to obtain a second triangular mesh set. In this set, each element is a triangle on a two-dimensional plane, and the coordinate data is derived from the projection results of the corresponding triangular mesh vertices in the original three-dimensional model.

[0088] Next, Boolean operations are performed on the two-dimensional coordinate data of each triangle in the second triangular mesh set to calculate their set relationship on the plane, thereby generating a set of closed polygons. These polygons are used to describe the overall boundary and internal structural features of the geological body on the two-dimensional plane. Among them, the outer contour represents the outermost boundary of the model on the plane, defined as a counterclockwise polygon; the island is a closed region that exists independently inside the outer contour, also defined as a counterclockwise polygon; the hole is a cavity region inside the outer contour, represented in a clockwise direction to distinguish it from subsequent geometric processing and topological analysis.

[0089] For example, projecting each element in the set T to a two-dimensional plane obtains a second triangular mesh set T' = {t1', t2',..., t n m polygons G' = {g1', g2',..., g m m polygons G' = {g1', g2',..., g Figure 4 As shown in FIG. 6, the distribution patterns of the outer contour, island, and hole can clearly reflect the segmentation results of the geological body on the two-dimensional plane, providing accurate geometric input for subsequent outline point extraction, intersection calculation, and triangular subdivision.

[0090] In some embodiments, based on the spatial intersection state of any two edges in the first triangular mesh set, the number of intersection points of any two edges projected onto a two-dimensional plane, and the intersection point coordinates, the intersection point coordinates of the roof and floor of the geological body three-dimensional model are determined. The spatial intersection state of any two edges includes no intersection, including:

[0091] If it is determined that the spatial intersection state of any two edges is no intersection, and the number of intersection points of any two edges projected onto a two-dimensional plane is 0, then the any two edges are ignored.

[0092] In the present embodiment, when the intersection state of any two edges in three-dimensional space is no intersection, and the number of intersection points of any two edges on a two-dimensional projection plane is 0, it means that the two edges have no intersection in space and projection, and cannot constitute a geometric intersection point of the roof and floor, so this pair of edges can be directly ignored to avoid invalid calculation.

[0093] For example, define a set PI For all elements in set T, calculate the intersection point of all edges. In space, calculate the intersection point of two edges, that is, calculate the corresponding point of the intersection point of the two edges projected onto the two-dimensional plane. Taking edges e1 and e2 as an example, if edges e1 and e2 do not intersect in three-dimensional space and do not have an intersection point on the two-dimensional projection plane, they are directly ignored and not included in the calculation range of the subsequent intersection point set.

[0094] In some embodiments, the intersection coordinates of the top and bottom plates of the geological body's three-dimensional model are determined based on the spatial intersection state of any two edges in the first triangular mesh set, the number of intersection points projected onto the two-dimensional plane by any two edges, and the coordinates of the intersection points. The spatial intersection state of any two edges includes intersection, including:

[0095] If the spatial intersection state of any two edges is determined to be intersection, and the number of intersection points when any two edges are projected onto the two-dimensional plane is one, then the intersection point is set as the second intersection point. Based on the second intersection point, a ray is emitted along a set direction, and the second intersection point set between the second intersection point and the first triangular mesh set is obtained. If the elevation value of the third intersection point with the highest elevation value in the second intersection point set is equal to the elevation value of the first intersection point, then the coordinates of the third intersection point are determined as the intersection coordinates of the top and bottom plates of the three-dimensional geological body model.

[0096] In this embodiment, when any two edges intersect in space and have only one intersection point on the two-dimensional projection plane, it indicates that the two edges have a real geometric intersection in three-dimensional space. In this case, this unique projection intersection point is defined as the second intersection point. And at that point Starting from a given position, emit a ray along a predetermined direction (usually perpendicular), and calculate the set of intersection points between this ray and the first triangular mesh set. The intersection point with the highest elevation value in this set is considered the third intersection point. If the elevation value is equal to the pre-calculated elevation value of the first intersection point, it indicates that the intersection point is actually located on the surface of the top and bottom plates. This intersection point can be confirmed as valid intersection point coordinates and added to the intersection point set P of the top and bottom plates of the geological body's 3D model. I Conversely, ignore it.

[0097] For example, define a set P I For all elements in set T, calculate the intersection point of all edges. In space, calculating the intersection point of two edges means calculating the corresponding points on the two edges that represent the intersection point projected onto the two-dimensional plane. Taking edges e1 and e2 as an example, edges e1 and e2 have one intersection point, meaning they intersect in space. Let the intersection point be... right Draw a ray upwards and calculate its intersection with point T. Let the point with the highest elevation among the intersection points be... like Elevation equal to Then the Add the set P I , otherwise ignore.

[0098] In some embodiments, based on the spatial intersection state of any two edges in the first triangular mesh set, the number of intersection points of any two edges projected onto a two-dimensional plane, and the intersection point coordinates, the intersection point coordinates of the top and bottom plates of the geological body three-dimensional model are determined, and the spatial intersection state of any two edges includes no intersection, including:

[0099] If it is determined that the spatial intersection state of any two edges is no intersection, and any two edges projected onto a two-dimensional plane have two intersection points, then the two intersection points are respectively set as a fourth intersection point and a fifth intersection point; respectively emitting a ray from the fourth intersection point and the fifth intersection point, obtaining a third intersection point set between the fourth intersection point and the fifth intersection point and the first triangular mesh set; if the elevation value of a sixth intersection point with the highest elevation value in the third intersection point set is the same as the elevation value of the fourth intersection point or the elevation value of the fifth intersection point, then the coordinates of the sixth intersection point are determined as the intersection point coordinates of the top and bottom plates of the geological body three-dimensional model.

[0100] In this embodiment, when the spatial intersection state of any two edges is no intersection, but the number of intersection points of any two edges on the two-dimensional projection plane is two, it indicates that there is an intersection phenomenon in the plane projection, but they are actually separated in space. At this time, the two intersection points are respectively set as a fourth intersection point and a fifth intersection point and respectively emit a ray from the two points in the set direction, and calculate the intersection point set with the first triangular mesh set. If the elevation value of the intersection point with the highest elevation value in the intersection point set (i.e. the sixth intersection point) is the same as the elevation value of the fourth intersection point or the fifth intersection point , then the sixth intersection point is confirmed as the valid top and bottom plate intersection point coordinates, and is added to the intersection point set P I ; otherwise, it is ignored.

[0101] Exemplarily, define a set P I , calculate the intersection points of all edges of all elements in the set T, calculate the intersection points of two edges in space, i.e. calculate the corresponding points of the intersection points of the two edges projected onto a two-dimensional plane on the two edges. Taking edges e1 and e2 as an example, edges e1 and e2 have two intersection points, i.e. edges e1 and e2 do not intersect in space, but intersect when projected onto a two-dimensional plane. Let the intersection points be and , calculate the intersection points with T by making a ray upward, and let the point with the highest elevation in the intersection points be If and have a point with the same elevation as , then ​Add the set P I Otherwise, ignore.

[0102] In some embodiments, the method further comprises:

[0103] If the triangular mesh vertex has a ray intersection with other triangular meshes in the first set of triangular meshes, ignore the triangular mesh vertex.

[0104] In this embodiment, if the ray intersects with other triangular meshes, it means that the vertex is occluded by other structures and does not participate in the construction of the top and bottom plate boundary, and can be ignored. Specifically, when a triangular mesh vertex in the first set of triangular meshes emits a ray in a set direction (for example, a vertical upward direction), if the ray intersects with other triangular meshes in the set except the triangular mesh where the vertex is located, it is determined that the vertex is occluded by other triangular faces and does not belong to the exposed vertex of the top and bottom plate. Such a vertex is located inside or below in the geometric sense and does not directly participate in the construction of the top and bottom plate boundary or surface, so it can be removed from the top and bottom plate vertex candidate set to ensure the accuracy and simplicity of subsequent model construction.

[0105] Exemplarily, the vertex coordinate set P R of the top and bottom plate of the three-dimensional model of the geological body is defined T For example, it is calculated whether there is an intersection with T, and if there is an intersection, the triangular mesh vertex p T is ignored.

[0106] In some embodiments, the method further comprises:

[0107] In some embodiments, the method further comprises:

[0108] In this embodiment, each triangular edge in the first set of triangular meshes is traversed, and it is determined whether both vertices of the triangular edge exist in the determined vertex coordinate set. If at least one vertex is found not to exist in the set, it means that the triangular edge is not composed of effective vertices of the top and bottom plate, and therefore can be directly ignored and not included in the candidate set of the top and bottom plate original triangular mesh topology line.

[0109] Exemplarily, the edges of the triangular meshes in the set T are traversed in turn, and the edge e i is taken as an example. i It is determined whether both vertices of e R are in the set P R , and if not, it is ignored.

[0110] Therefore, by the vertex validity filtering, irrelevant triangle edges can be excluded in advance, the calculation amount of subsequent boundary judgment and topology construction can be reduced, and false edges can be prevented from participating in the generation of the top and bottom plate model, so that the generated topology structure line set is consistent with the actual top and bottom plate surface.

[0111] In some embodiments, the method further comprises:

[0112] Each triangle edge in the first triangle mesh set is traversed, and if it is determined that both vertices of the triangle edge are in the vertex coordinate set, and it is determined that there is an intersection point in the intersection point coordinate set on the triangle edge, the triangle edge is ignored.

[0113] For the triangle edges that remain after the first screening, it is further determined whether they have a coincidence relationship with the intersection point coordinate set. Specifically, if both vertices of an edge are in the vertex coordinate set, but there is an intersection point on the edge in the intersection point coordinate set, it indicates that the edge is geometrically cut by the intersection point, and should not be used as the original continuous topology structure line, so the edge needs to be ignored.

[0114] In this way, the edge cut by the intersection point can be prevented from being directly used as a topology constraint line, to prevent the occurrence of false connection relationships or geometric conflicts in the subsequent subdivision process, and to maintain the correctness and consistency of the subdivision mesh.

[0115] For example, the edges of the triangle meshes in the set T are traversed in turn, and the edge e i is taken as an example to determine whether both vertices of e i are in the set P R . If both vertices are in the set P R , the next step is performed, that is, it is determined whether there is a point in P I on the edge e i in the two-dimensional projection plane, and if there is, the edge e i is ignored.

[0116] In some embodiments, the triangle subdivision is performed based on the intersection point coordinate set, the vertex coordinate set, the boundary contour line set of the top and bottom plate of the geological body three-dimensional model, and the original triangle mesh topology structure line constraint set of the top and bottom plate of the geological body three-dimensional model to generate the top and bottom plate of the geological body three-dimensional model, including:

[0117] The line constraint Delaunay triangle subdivision is performed based on the intersection point coordinate set, the vertex coordinate set, the boundary contour line set of the top and bottom plate of the geological body three-dimensional model, and the original triangle mesh topology structure line constraint set to generate a set of subdivided triangle meshes;

[0118] triangular meshes belonging to the outer contour of the boundary contour line set of the top / bottom plate of the 3D geological model are deleted; and the top / bottom plate of the 3D geological model is generated based on the set of the deleted dissected triangular meshes.

[0119] In the set of the generated dissected triangular meshes, there may be redundant meshes located outside the top / bottom plate. Therefore, in this embodiment, the triangular meshes located outside the boundary contour line of the top / bottom plate in the dissected result are deleted by the outer contour judgment, and only the effective mesh area inside the outer contour is reserved.

[0120] After the above dissecting and screening operations, the final top / bottom plate mesh structure of the 3D geological model can be obtained based on the set of the deleted dissected triangular meshes, which not only retains the geometric and topological characteristics of the original model, but also has complete dissected continuity, thereby providing a reliable data basis for subsequent visualization, analysis and calculation.

[0121] Next, the present application will be described in detail in combination with an application example.

[0122] Based on the above problems, therefore, it is urgent to introduce a method for extracting the top / bottom plate of a 3D geological model to achieve fine, convenient and efficient extraction of the top / bottom plate of the 3D geological model, and the technical effects include:

[0123] 1. The problem of time-consuming and laborious and low efficiency in extracting the top / bottom plate by human-computer interaction is solved, and the top plate or bottom plate model can be updated in time after the 3D geological model is dynamically updated;

[0124] 2. The top / bottom plate model extraction at the junction of multiple geological models can be supported, and the logical error of spatial intersection of the top / bottom plate model at the junction of multiple geological models can be avoided.

[0125] 3. The extracted top / bottom plate can maintain the triangular mesh topological structure consistent with the original 3D model to the greatest extent, achieving the purpose of fine extraction of the top / bottom plate 3D model.

[0126] This application example takes the extraction of the top plate as an example, Figure 7 the flowchart for extracting the top plate of the 3D geological model is shown in FIG. 7, and the technical solution of the application example will be explained in combination with Figure 7 the flowchart for extracting the top plate of the 3D geological model is shown in FIG. 7, and the technical solution of the application example will be explained in combination with

[0127] Step 701: Extracting the boundary contour line of the top plate of the 3D geological model.

[0128] In actual application, the input data for extracting the top / bottom plate of the 3D geological model includes multiple 3D geological models, as shown in FIG. 6. It is assumed that the multiple 3D geological models are composed of n triangular meshes, as shown in FIG. 7. Figure 2 Figure 3 ​As shown, set T = {t1, t2,..., t n} (i.e. the first triangular mesh set mentioned above).

[0129] Obtain a first triangular mesh set T = {t1, t2,..., t n} corresponding to the plurality of geological body three-dimensional models; the first triangular mesh set includes three-dimensional coordinate data of a plurality of triangular meshes;

[0130] Project each element in set T to a two-dimensional plane to obtain set T' = {t1', t2',..., t n '}; perform a Boolean operation on all elements in set T' on the plane to obtain a union of m polygon sets G' = {g1', g2',..., g m '}; define the outer contour or "island" in G' as a counterclockwise polygon, and the "hole" in G' as a clockwise polygon, as shown in Figure 4 .

[0131] For each vertex v of each polygon in G', make a ray upward, calculate the intersection point with T, and assign the elevation of the highest point in the intersection point to v to obtain a roof boundary contour line G = {g1, g2,..., g m} of the geological body three-dimensional model (i.e. the set of roof and floor boundary contour lines of the geological body three-dimensional model), as shown in Figure 5 .

[0132] Step 702: geological body three-dimensional model intersection point processing.

[0133] In actual application, define set P I (i.e. the intersection point coordinate set of the roof and floor of the geological body three-dimensional model mentioned above); calculate the intersection point of all edges of all elements in set T; in space, calculate the intersection point of two edges, i.e. calculate the corresponding point of the intersection point of the two edges projected onto the two-dimensional plane on the two edges. Take edges e1 and e2 as an example, there are three cases for the intersection point of edges e1 and e2:

[0134] (1) Edges e1 and e2 have no intersection point, which is ignored;

[0135] (2) Edges e1 and e2 have one intersection point, i.e. edges e1 and e2 intersect in space, and the intersection point is set as For , make a ray upward, calculate the intersection point with T, and set the point with the highest elevation in the intersection point as If the elevation of is equal to , add to set P I , otherwise ignore it;

[0136] (3) Edges e1 and e2 have two intersection points, i.e. edges e1 and e2 are not intersected in space, but intersected in the two-dimensional plane, and the intersection points are and For , a ray is made upward to calculate the intersection point with T, and the point with the highest elevation in the intersection point is If and have a point with the same elevation as , then is added to the set P I , otherwise it is ignored.

[0137] Step 703: Processing of the vertices of the three-dimensional geological body model.

[0138] In practical applications, the set P R (i.e. the set of vertex coordinates of the top and bottom plates of the three-dimensional geological body model) is defined, and a ray is made upward for all vertices in T, taking p T as an example, whether there is an intersection point with T is calculated, if there is no intersection point, then p T is added to the set P R , otherwise it is ignored.

[0139] Step 704: Construction of the original triangular mesh topological structure line constraint of the top plate of the three-dimensional geological body model.

[0140] In practical applications, the specific construction steps are as follows:

[0141] (1) Define the original triangular mesh topological structure line constraint set C of the top plate of the three-dimensional geological body model (i.e. the original triangular mesh topological structure line constraint set of the top and bottom plates of the three-dimensional geological body model);

[0142] (2) Traverse the edges of the triangular mesh in set T in turn, taking edge e i as an example, whether both vertices of e i are in set P R is determined, if both are in set P R , the next step is performed, otherwise it is ignored;

[0143] (3) Determine whether there is a point in P I on edge e i in the two-dimensional projection plane, if there is not, then edge e i is added to set C, otherwise it is ignored.

[0144] Step 705: Construction of the top plate of the three-dimensional geological body model.

[0145] In practical applications, the union of P I , P R , G and C is subjected to line-constrained Delaunay triangulation to obtain TU Delete T U The triangular mesh outside the outer contour of G is used to obtain the top plate of the 3D model of the geological body, such as... Figure 6 As shown.

[0146] Based on the hardware implementation of the above program modules, and in order to implement the method of the embodiments of this application, the embodiments of this application also provide an electronic device. Figure 8 The diagram shows only an exemplary structure of the electronic device, not the entire structure; implementation is possible as needed. Figure 8 The diagram shows part or all of the structure. For example... Figure 8 As shown, the electronic device 800 provided in this application embodiment includes: at least one processor 801, a memory 802, a user interface 803, and at least one network interface 804. The various components in the electronic device 800 are coupled together via a bus system 805. It can be understood that the bus system 805 is used to implement communication between these components. In addition to a data bus, the bus system 805 also includes a power bus, a control bus, and a status signal bus. However, for clarity, in... Figure 8 The various buses are all labeled as bus system 805. Among them, the user interface 803 may include a display, keyboard, mouse, trackball, click wheel, buttons, touchpad, or touch screen, etc.

[0147] The memory 802 in this embodiment is used to store various types of data to support the operation of the electronic device. Examples of such data include any computer program used to operate on the electronic device.

[0148] The method for extracting the top and bottom plates of the three-dimensional model of the geological body of the electronic device disclosed in the embodiments of the present application can be applied in the processor 801 or implemented by the processor 801. The processor 801 can be an integrated circuit chip with a signal processing capability. In the implementation process, each step of the method for optimizing the arrangement of the mine early warning device of the electronic device can be completed by the integrated logic circuit of the hardware in the processor 801 or the instructions in the form of software. The processor 801 mentioned above can be a general processor, a digital signal processor (DSP), or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, etc. The processor 801 can implement or execute the methods, steps, and logic block diagrams disclosed in the embodiments of the present application. The general processor can be a microprocessor or any conventional processor, etc. In combination with the steps of the method disclosed in the embodiments of the present application, the hardware decoding processor can be directly embodied to execute the steps, or the hardware and software modules in the decoding processor can be combined to execute the steps. The software module can be located in the storage medium, and the storage medium is located in the memory 802. The processor 801 reads the information in the memory 802, and in combination with the hardware, the steps of the method for extracting the top and bottom plates of the three-dimensional model of the geological body of the electronic device provided in the embodiments of the present application are completed.

[0149] In the example embodiments, the electronic device can be implemented by one or more application specific integrated circuits (ASICs), DSPs, programmable logic devices (PLDs), complex programmable logic devices (CPLDs), field programmable gate arrays (FPGAs), general-purpose processors, controllers, microcontrollers (MCUs), microprocessors (Microprocessors), or other electronic elements, for executing the foregoing method.

[0150] It is to be understood that the memory 802 can be volatile or nonvolatile memory, or both. In one embodiment, the nonvolatile memory can be read only memory (ROM), programmable read only memory (PROM), erasable programmable read only memory (EPROM), electrically erasable programmable read only memory (EEPROM), a magnetic random access memory (FRAM), flash memory, a magnetic surface memory, an optical disk, or a compact disk read only memory (CD-ROM). The magnetic surface memory can be a magnetic disk memory or the volatile memory can be a random access memory (RAM) used as an external cache. By way of example, and not limitation, many forms of RAM can be used, such as static random access memory (SRAM), synchronous static random access memory (SSRAM), dynamic random access memory (DRAM), synchronous dynamic random access memory (SDRAM), double data rate synchronous dynamic random access memory (DDR SDRAM), enhanced synchronous dynamic random access memory (ESDRAM), synchronous link dynamic random access memory (SLDRAM), and direct rambus random access memory (DRRAM). The memory described herein is intended to include, without being limited to, these and any other suitable types of memory.

[0151] In an example embodiment, the embodiments of the present application further provide a computer storage medium, specifically a computer readable storage medium, which stores a computer program executable by a processor to complete the steps of the method of the embodiments of the present application. The computer readable storage medium can be a ROM, a PROM, an EPROM, an EEPROM, a flash memory, a magnetic surface memory, an optical disc, or a CD-ROM memory, etc.

[0152] In an example embodiment, the embodiments of the present application further provide a computer program product comprising a computer program executable by the processor 801 of an electronic device to complete the steps described in the method of the embodiments of the present application.

[0153] It should be noted that "first", "second", etc. are used to distinguish similar objects, and do not necessarily mean a specific order or sequence.

[0154] In addition, the technical solutions described in the embodiments of the present application can be combined arbitrarily without conflict. The above is merely a specific implementation of the present application, but the protection scope of the present application is not limited thereto, and any person skilled in the art can easily think of changes or replacements within the technical range disclosed in the present application, which should be covered within the protection scope of the present application.

Claims

1. A method for extracting the top and bottom plates of a three-dimensional geological body model, characterized in that, The method includes: A polygon set is obtained by projecting a first triangular mesh set of multiple geological body 3D models onto a 2D plane. The polygons in the polygon set include at least one of the following types: outer contour, island, and hole. Rays are emitted along a set direction for each polygon vertex in the polygon set to determine a first set of intersection points between the polygon vertex and the first triangular mesh set. Based on the first intersection point with the highest elevation value in the first set of intersection points and the 2D coordinate data of the polygon vertex, coordinate data of contour points are generated. Based on the coordinate data of each contour point, a set of boundary contour lines of the top and bottom plates of the geological body 3D model is generated. Based on the spatial intersection state of any two sides in the first triangular mesh set, the number of intersection points and the coordinates of the intersection points when the two sides are projected onto the two-dimensional plane, the coordinates of the intersection points of the top and bottom plates of the geological body's three-dimensional model are determined; based on the coordinates of each intersection point, the set of intersection point coordinates of the top and bottom plates of the geological body's three-dimensional model is determined. For each vertex of the first triangular mesh set, a ray is drawn along a set direction. If the vertex of the triangular mesh does not have a ray intersection with any other triangular mesh in the first triangular mesh set, then the vertex of the triangular mesh is determined as the vertex coordinates of the top and bottom plates of the geological body 3D model. Based on the vertex coordinates of each top and bottom plate, a vertex coordinate set of the top and bottom plates of the geological body 3D model is generated. Traverse each triangle edge in the first triangular mesh set. If it is determined that both vertices of the triangle edge are in the vertex coordinate set and that there is no intersection point on the triangle edge in the intersection coordinate set, then the triangle edge is determined to be the original triangular mesh topology line of the top and bottom plates of the geological body three-dimensional model. Based on each of the original triangular mesh topology lines, construct the original triangular mesh topology line constraint set of the top and bottom plates of the geological body three-dimensional model. Triangulation is performed based on the set of intersection coordinates, the set of vertex coordinates, the set of boundary contour lines of the top and bottom plates of the three-dimensional geological body model, and the set of original triangular network topology constraints of the top and bottom plates of the three-dimensional geological body model to generate the top and bottom plates of the three-dimensional geological body model.

2. The method according to claim 1, characterized in that, The method further includes: Obtain a second triangular mesh set after projecting the first triangular mesh set of multiple geological body 3D models onto a 2D plane. The second triangular mesh set includes the 2D coordinate data of multiple triangles. Boolean operations are performed on the two-dimensional coordinate data of each second triangular mesh in the second triangular mesh set to generate the polygon set. The polygon set includes at least one of the following: coordinate data of the outer contour, coordinate data of the island, and coordinate data of the hole. The outer contour and the island are counterclockwise polygons, and the hole is a clockwise polygon.

3. The method according to claim 1, characterized in that, The spatial intersection state of any two edges in the first triangular mesh set, the number of intersection points projected onto the two-dimensional plane by the two edges, and the coordinates of the intersection points are used to determine the intersection coordinates of the top and bottom plates of the three-dimensional model of the geological body. The spatial intersection state of any two edges includes non-intersection, including: If it is determined that the spatial intersection state of any two edges is non-intersection, and the number of intersection points of any two edges projected onto the two-dimensional plane is 0, then the two edges are ignored.

4. The method according to claim 1, characterized in that, The spatial intersection state of any two edges in the first triangular mesh set, the number of intersection points projected onto the two-dimensional plane by the arbitrary two edges, and the coordinates of the intersection points are used to determine the intersection coordinates of the top and bottom plates of the three-dimensional model of the geological body. The spatial intersection state of any two edges includes intersection, including: If it is determined that the spatial intersection state of any two sides is intersecting, and the number of intersection points of any two sides projected onto the two-dimensional plane is one, then the intersection point is set as the second intersection point. Based on the second intersection point, a ray is emitted along a set direction, and the second intersection point set between the second intersection point and the first triangular mesh set is obtained. If the elevation value of the third intersection point with the highest elevation value in the second intersection point set is equal to the elevation value of the first intersection point, then the coordinates of the third intersection point are determined as the intersection coordinates of the top and bottom plates of the three-dimensional geological body model.

5. The method according to claim 1, characterized in that, The spatial intersection state of any two edges in the first triangular mesh set, the number of intersection points projected onto the two-dimensional plane by the two edges, and the coordinates of the intersection points are used to determine the intersection coordinates of the top and bottom plates of the three-dimensional model of the geological body. The spatial intersection state of any two edges includes non-intersection, including: If it is determined that the spatial intersection state of any two edges is non-intersecting, and the projection of any two edges onto the two-dimensional plane has two intersection points, then the two intersection points are set as the fourth intersection point and the fifth intersection point, respectively; rays are emitted from the fourth intersection point and the fifth intersection point to obtain the third intersection point set between the fourth intersection point and the fifth intersection point and the first triangular mesh set; if the elevation value of the sixth intersection point with the highest elevation value in the third intersection point set is the same as the elevation value of the fourth intersection point or the elevation value of the fifth intersection point, then the coordinates of the sixth intersection point are determined as the intersection coordinates of the top and bottom plates of the three-dimensional geological body model.

6. The method according to claim 1, characterized in that, The method further includes: If a vertex of the triangular mesh intersects with another triangular mesh in the first triangular mesh set via a ray, then the vertex of the triangular mesh is ignored.

7. The method according to claim 1, characterized in that, The method further includes: Traverse each triangle edge in the first triangular mesh set. If it is determined that at least one of the two vertices of the triangle edge is not in the vertex coordinate set, then the triangle edge is ignored.

8. The method according to claim 1, characterized in that, The method further includes: Traverse each triangle edge in the first triangular mesh set. If it is determined that both vertices of the triangle edge are in the vertex coordinate set and that there is an intersection point on the triangle edge in the intersection coordinate set, then ignore the triangle edge.

9. The method according to claim 1, characterized in that, The process of triangulation based on the set of intersection coordinates, the set of vertex coordinates, the set of boundary contour lines of the top and bottom plates of the 3D geological body model, and the set of line constraints of the original triangular network topology of the top and bottom plates of the 3D geological body model to generate the top and bottom plates of the 3D geological body model includes: Based on the set of intersection coordinates, the set of vertex coordinates, the set of boundary contours of the top and bottom plates of the three-dimensional geological body model, and the set of line constraints of the original triangular mesh topology of the top and bottom plates of the three-dimensional geological body model, line constraint Delaunay triangulation is performed to generate a set of triangular meshes. Delete other triangular meshes in the triangular mesh set that are outside the outer contour of the boundary contour line set of the top and bottom plates of the geological body 3D model; and generate the top and bottom plates of the geological body 3D model based on the deleted triangular mesh set.

10. An electronic device, characterized in that, include: A processor and memory for storing computer programs that can run on the processor, wherein, The processor, when running a computer program, performs the steps of the method according to any one of claims 1 to 9.

Images