CN105931297A - Data processing method applied to three-dimensional geological surface model - Google Patents

Abstract

The invention discloses a data processing method applied to a three-dimensional geological surface model. The data processing method includes the following steps that: a sampling point set is obtained; interpolation is performed on the sampling point set through using an interpolation algorithm, and a gain type sampling point set can be generated, so that a geological three-dimensional point cloud structure can be recovered; the gain type sampling point set is processed by using Delaunay triangulation algorithm, so that a triangle set can be generated; malformed triangles are searched in the triangle set; the malformed triangles are processed; the triangle set of which the malformed triangles are processed is optimized through using the windowed sinc interpolation kernel function; the optimized triangle set is processed according to a surface extending algorithm and a surface intersection algorithm, so that a surface boundary line can be generated; and a spatial surface set at the surface boundary line is generated according to a boundary filtering algorithm so as to be used for characterizing a spatial section. The spatial surface set is corresponding to one independent spatial section in space, so that formation shape can be recovered, and the quantitative study of a geological structure can be benefitted.

Description

Data processing method in three-dimensional geological surface model

Technical field

The application relates to GIS-Geographic Information System field, particularly relates to the data in a kind of three-dimensional geological surface model Processing method.

Background technology

When GIS-Geographic Information System carries out in existing three-dimensional geological surface model data process, stressing of research Point is to carry out geologic structure analogue simulation, i.e. lays particular emphasis on and renders geologic structure and visual place Reason.

During realizing prior art, inventor finds that in prior art, at least there are the following problems:

In the three-dimensional geological surface model that GIS-Geographic Information System builds, the standardization of dough sheet unit is poor, precision is low, Geologic structure cannot be carried out quantitative research.Such as, split in geologic structure according to three-dimensional geological surface model Seam, tomography etc. are studied.

Summary of the invention

The embodiment of the present application provides the edge filtering data processing method in a kind of three-dimensional geological surface model, should Data processing method dough sheet unit standardization is strong, precision is high, may be used for the quantitative research of geologic structure.Tool Body, the data processing method in a kind of three-dimensional geological surface model, comprise the following steps:

Obtain sampled point set；

Use interpolation algorithm, sampled point set is carried out interpolation and generates gain-type sampled point set, to recover ground The three-dimensional point cloud structure of matter；

Use delaunay subdivision algorithm, processing gain type sampled point set, generate triangle sets；

Search the lopsided triangle in triangle sets；

Process deformity triangle；

Using windowed sinc interpolation kernel function, optimization processes the triangle sets after deformity triangle；

According to surface extending algorithm and surface intersection algorithm, process the triangle sets after optimizing, generate curved surface Boundary line；

According to edge filtering algorithm, generate the space curved surface set at surface boundary line, break in order to characterize space Face.

The method that the data that the embodiment of the present application provides process, at least has the advantages that

Space curved surface is integrated into a spatially corresponding independent space section, such that it is able to renwing stratum shape, Be conducive to the quantitative research to geologic structure.

Accompanying drawing explanation

Accompanying drawing described herein is used for providing further understanding of the present application, constitutes the part of the application, The schematic description and description of the application is used for explaining the application, is not intended that the improper limit to the application Fixed.In the accompanying drawings:

Interpolated data processing method stream in the three-dimensional geological surface model that Fig. 1 provides for the embodiment of the present application Cheng Tu.

Gridded data processing method in the three-dimensional geological surface model that Fig. 2 provides for the embodiment of the present application Flow chart.

Lopsided triangle data in the three-dimensional geological surface model that Fig. 3 provides for the embodiment of the present application processes Data processing method flow chart.

Surface extending data process side in the three-dimensional geological surface model that Fig. 4 provides for the embodiment of the present application Method data processing data process flow figure.

Surface intersection data process side in the three-dimensional geological surface model that Fig. 5 provides for the embodiment of the present application Method data processing data process flow figure.

Edge filtering data process side in the three-dimensional geological surface model that Fig. 6 provides for the embodiment of the present application Method flow chart.

Data processing method flow chart in the three-dimensional geological surface model that Fig. 7 provides for the embodiment of the present application.

Detailed description of the invention

For making the purpose of the application, technical scheme and advantage clearer, specifically real below in conjunction with the application Execute example and technical scheme is clearly and completely described by corresponding accompanying drawing.Obviously, described Embodiment is only some embodiments of the present application rather than whole embodiments.Based on the enforcement in the application Example, the every other enforcement that those of ordinary skill in the art are obtained under not making creative work premise Example, broadly falls into the scope of the application protection.

In geological research, three-dimensional geological model can be set up, in order to visualize geologic structure, such that it is able to Visualization represents geologic structure.

The geographical location information that geological prospecting is obtained, and insert in geographical location information according to certain rule Enter some geographical location information with reference to character, the Visualization Model being made up of some dot matrix can be set up.This Plant model and can substantially sketch the contours of the general profile of geologic structure, it is impossible to carry out further quantitative research.

Further, can be by geographical location information, i.e. several spatial point couple together formation dough sheet, also It is exactly the planar slice formed by some spatial point, it is thus possible to the expression geologic structure that becomes more meticulous further.

The application provide embodiment in, geological prospecting can be obtained geographical location information, i.e. several Spatial point, divides and forms some sealene triangle faces unit, or be called sealene triangle dough sheet.So After, can will divide the sealene triangle dough sheet formed, join together to get up the surface of matching geologic structure.

But, owing to geological prospecting obtains the randomness of geographical location information, the sealene triangle face of formation Differing greatly between sheet, be in particular in, the triangle length of side, triangle area fluctuating margin are relatively big, unfavorable Process in follow-up.

In the embodiment of the present application, it is provided that a kind of geology inserting some reference character in geographical location information The method of positional information, so that sealene triangle dough sheet meets certain standardization, consequently facilitating follow-up Process.

Refer to Fig. 1, the embodiment of the present application provides the data processing method in a kind of three-dimensional geological surface model, Comprise the following steps:

S100: obtain sampled point set.

Sampled point set, specifically can show as the geographical location information obtained by geological prospecting.Geology is surveyed Visit the geographical location information obtained, concrete, can be borehole data, geological data, slice of data etc.. Sampled point set PS={p can be remembered₁, p₂…p_n, p_i=(x_i, y_i, z_i), i ∈ [1, n], n represents sampled point Number.

S200: calculate the curvature of each sampled point in sampled point set.

Sampled point set PS is discrete spatial point data acquisition system, does not comprise the topology between any sampled point Relation information.Can be by sampled point set being carried out two dimension Delaunay subdivision, it is thus achieved that between sampled point Open up benefit relation.

In a kind of embodiment that the application provides, calculate the curvature of each sampled point in sampled point set, tool Body includes:

Calculate the Gaussian curvature of each sampled point in sampled point set.

Geographical location information according to sampled point set and its topological structure, the curvature of each sampled point can be led to Cross following equation to obtain:

α(p_i)=K (p_i)₁×K(p_i)₂。

Wherein, α (p_i) represent p_iThe Gaussian curvature of point, K (p_i)₁With K (p_i)₂Represent p respectively_iTwo masters of point Curvature.Wherein p_iThe principal curvatures of point refers to, at p_iThe normal plane of the different directions at Dian and surface intersection curve Curvature in maximum and minima.

If intersection curve is given with parametric form:

S (t)=(x (t), y (t), z (t)), then K (p_i)₁With K (p_i)₂It is exactly that this curve is at a p_iThe curvature at place is Big value and minima, be denoted as F (x, y, z).

$F(x,y,z)=\frac{\sqrt{{(\frac{{\∂}^{2}z\∂y}{{\∂}^{2}t\∂t}-\frac{{\∂}^{2}y\∂z}{{\∂}^{2}t\∂t})}^2+{(\frac{{\∂}^{2}x}{{\∂}^{2}t}\frac{\∂z}{\∂t}-\frac{{\∂}^{2}z}{{\∂}^{2}t}\frac{\∂x}{\∂t})}^2+{(\frac{{\∂}^{2}y\∂x}{{\∂}^{2}t\∂t}-\frac{{\∂}^{2}x\∂y}{{\∂}^{2}t\∂t})}^2}}{{(\frac{\∂x^2}{\∂t}+\frac{\∂y^2}{\∂t}+\frac{\∂z^2}{\∂t})}^{3/2}}.$

Thus can obtain the sampled point set cuvPS={p of subsidiary Gaussian curvature₁, p₂…p_n, wherein

p_i=(x_i, y_i, z_i, α_i), i ∈ [1, n], n represent the number of sampled point, α_iRepresent some p_iPlace Gauss is bent Rate.

S300: sampled point set is projected to horizontal plane.

S400: the control point generating the horizontal plane projection for covering sampled point set at horizontal plane is gathered.

Control point CP_i, i=1,2 ... N.N represents the number at control point.Control point is at prime stratum The point of horizontal plane, these control point will control the interpolation of surface of stratum.

The coordinate at control point is calculated respectively according to below equation:

${\mathrm{XCP}}_i=x_{min}+\frac{x_{max}-x_{min}}{x_{revol}+1}\×i;$

${\mathrm{YCP}}_j={\mathrm{\γ}}_{min}+\frac{y_{max}-y_{min}}{y_{revol}+1}\×j;$

Wherein, i=1,2 ... M；J=1,2 ... N.

XCP_iRepresent the X-coordinate of i-th point in the X direction；

x_minRepresent the minima of X-coordinate in control point；

x_maxRepresent the maximum of X-coordinate in control point；

x_revolRepresent the resolution of X-direction；

YCP_iRepresent the Y coordinate of jth point in the Y direction；

y_minRepresent the minima of Y coordinate in control point；

y_maxRepresent the maximum of Y coordinate in control point；

y_revolRepresent the resolution of Y-direction.

Control point CP_iFor controlling at rectangular area rect_iInterior interpolation point.Rectangular area rect_iWith control point CP_iCentered by point, respectively with one step x_interAnd y_interFor the length of side.

In order to the interpolation obtained, the one step of X-direction and the one step of Y-direction should be equal, I.e. x_inter=y_inter, say, that described parameter meets following equation:

$y_{revol}=\frac{(y_{max}-y_{min})\×(x_{revol}+1)}{x_{max}-x_{min}}-1.$

S500: according to the curvature at the control point that the curvature estimation control point of sampled point is gathered.

According to interpolation method, it is possible to obtain control point set CPCS={p₁, p₂... p_nGaussian curvature, Wherein, p_i=(x_i, y_i, α_i), i ∈ [1, n], n represent the number at control point, α_iRepresent some p_iThe Gauss at place Curvature.

Concrete interpolation method can use distance inverse ratio interpolation algorithm, Kriging regression algorithm etc..

S600: calculate the interpolation density at control point according to curvature density transformation model.

Interpolation density according to below equation calculating control point:

Wherein, δ_iRepresent control point CP_iThe interpolation density at place；

A is conversion factor, in order to control density and the dispersion degree of entirety；

α_iRepresent CP_iThe curvature at place；

α_minAnd α_maxRepresent maximum and the minima of the curvature at all control point.

Such that it is able to obtain having the control point set of two dimensional surface coordinate and interpolation density CPDS={p₁, p₂... p_n, wherein, p_i=(x_i, y_i, δ_i), i ∈ [1, n], n represent the individual of control point Number, δ_iRepresent some p_iThe interpolation density at place.

S700: determine interpolation point set according to the interpolation density at control point.

Further, in the another embodiment that the application provides, determine slotting according to the interpolation density at control point Value point set, specifically includes:

Determine the interpolation point in perspective plane；

Interpolation point in perspective plane is mapped the interpolation point generating surface model.

Further, in the another embodiment that the application provides, the interpolation point in perspective plane is determined, specifically Including:

Determine the area of perspective plane dough sheet unit；

Area according to described dough sheet unit and interpolation density, calculate the interpolation point being positioned at perspective plane dough sheet unit；

Interpolation point is uniformly distributed in dough sheet unit, it is thus achieved that interpolation point set.

Further, in the another embodiment that the application provides, the interpolation point in perspective plane is mapped and generates The interpolation point of surface model, specifically includes:

Determine the third dimension coordinate of the interpolation point of surface model according to below equation, thus obtain surface model Interpolation point set:

$z\left(p\right(x,y\left)\right)=\frac{{\mathrm{\Σ}}_{i=1}^n\frac{z_i}{{\[d_i(x,y)\]}^{\mathrm{\μ}}}}{{\mathrm{\Σ}}_{i=1}^n\frac{1}{{\[d_i\left(xy\right)\]}^{\mathrm{\μ}}}}.$

Wherein,(x, y) to sampled point set to represent interpolation point p Middle i-th point p_iDistance；

The power exponent μ value of weight coefficient is 2, represents that interpolation point and sampled point are several in the Europe of horizontal plane In obtain distance.

Concrete, for each control point in control point set CPDS, its have two-dimensional coordinate information and Interpolation density information, can be centered by control point, with x_interAnd y_interA little square is constructed for the length of side Shape region.At each little rectangle inside, according to the control point density in four end points of rectangle and rectangular area The method utilizing linear interpolation generates the stratum projection plane two dimension point set of rectangle inside RP_i={ p₁, p₂... p_ni, wherein, RP_iRepresent the two-dimentional point set in the little rectangle of i-th. p_j=(x_j, y_j) represent the coordinate of this little rectangle inside jth point.Ni=δ_i ²In representing the little rectangle of i-th The number of the point in portion.It can thus be appreciated that the calculating of the two-dimentional point set RPS generated in the projection plane of all of stratum Method is shown in formula

Then, use distance inverse ratio interpolation and Kriging regression algorithm, be combined into sample with sampling point set, use RPS point set carries out interpolation, recovers three-dimensional point cloud structure.

Determine the third dimension coordinate of the interpolation point of surface model according to below equation, thus obtain surface model Interpolation point set:

$z\left(p\right(x,y\left)\right)=\frac{{\mathrm{\Σ}}_{i=1}^n\frac{z_i}{{\[d_i(x,y)\]}^{\mathrm{\μ}}}}{{\mathrm{\Σ}}_{i=1}^n\frac{1}{{\[d_i(x,y)\]}^{\mathrm{\μ}}}}.$

Wherein,(x, y) to sampled point set to represent interpolation point p Middle i-th point p_iDistance；

The power exponent μ value of weight coefficient is 2, represents that interpolation point and sampled point are several in the Europe of horizontal plane In obtain distance.

In the embodiment of the present application, interpolation point generates according to the interpolation density at control point, and for sampling Point supplements, such that it is able to improve the quality of the dough sheet unit generated, and then improves the specification of dough sheet unit Property and precision, in order to for the quantitative research of geologic structure.

In the embodiment of the present application, further, also provide for one and geological prospecting obtained geographical location information, Several spatial point i.e., divide and form some sealene triangle faces unit, or be called sealene triangle Dough sheet.It is then possible to the sealene triangle dough sheet formed will be divided, join together to get up matching geologic structure Surface.

Refer to Fig. 2, the embodiment of the present application provides at the gridded data in a kind of three-dimensional geological surface model Reason method, comprises the following steps:

S101: generate one and comprise the initial boundary triangle of all sampled points in gain-type sampled point set.

Gain-type sampled point set comprises original sample point and interpolation point.

S201: select each sampled point in gain-type sampled point set one by one, searches circumscribed circle and comprises described The triangle sets of sampled point.

S301: delete triangle sets, forms cavity in initial boundary triangle.

S401: connect described sampled point and each the limit in cavity, form gridding triangle sets.

S501: have the gridding three of same vertices in Delete Mesh triangle sets with initial boundary triangle Dihedral.

In the embodiment of the present application, gain-type sampled point set comprises original sample point and interpolation point, thus, Standardization and the precision of the blocking of geological surface are significantly improved, and the quantization that may be used for geologic structure is ground Study carefully.

Further, in the another embodiment that the application provides,

Select each sampled point in gain-type sampled point set one by one, search circumscribed circle and comprise described sampled point Triangle sets, specifically include:

The guide for use skill of handling needles searches the triangle comprising described sampled point；

Search the adjacent triangle of described triangle；

When all of its neighbor triangle does not comprises described sampled point, form circumscribed circle and comprise described sampled point Triangle sets.

Use the method can significantly improve the efficiency that data process.

In the embodiment of the present application, although the sealene triangle dough sheet total quality formed above is higher, but If the interpolation density difference between the rectangular area that two control point are formed is excessive, it will cause these two Mesh quality between rectangle can not meet real needs.

Refer to Fig. 3, the embodiment of the present application provides the lopsided triangle number in a kind of three-dimensional geological surface model According to processing method, comprise the following steps:

S102: search deformity triangle.

S202: the center of circle of deformity triangle circumscribed circle is added in sampled point set.

S302: search the adjacent sampled point of each sampled point in sampled point set.

S402: use adjacent sample point coordinate to recalculate the coordinate of described sampled point according to interpolation algorithm.

Further, in the another embodiment that the application provides, search deformity triangle, specifically include:

According to formula

Search deformity triangle；

Wherein, δ is threshold parameter.

Further, in the another embodiment that the application provides, use adjacent sampled point according to interpolation algorithm Coordinate recalculates the coordinate of described sampled point, specifically includes:

Windowed sine interpolation kernel function is used to recalculate the coordinate of sampled point.

Use interpolation algorithm that deformity triangle is modified, therefore, the standardization of the blocking of geological surface It is significantly improved with precision, may be used for the specification of the blocking of the quantitative research geological surface of geologic structure Property and precision are significantly improved, and may be used for the quantitative research of geologic structure.

Refer to Fig. 4, the embodiment of the present application provides the surface extending data in a kind of three-dimensional geological surface model Processing method, comprises the following steps:

S103: find out all boundary edge set and edge triangles set in surface of stratum triangle gridding.

Find out all boundary edge set BES and edge triangles BTS set in surface of stratum triangle gridding.

S203: select the boundary edge in boundary edge set one by one, this boundary edge belongs to an edge triangles collection An edge triangles in conjunction.

S303: determine the starting point of extension according to two end points of described boundary edge.

Further, in the another embodiment that the application provides, two end points according to described boundary edge are true The fixed starting point extended, specifically includes:

When two end points of boundary edge areWithDetermine rising of extension Point is

S403: determine the step-length of extension according to described edge triangles.

Further, in the another embodiment that the application provides, extension is determined according to described edge triangles Step-length, specifically include:

Work as e₁、e₂、e₃When being respectively the length of side of described edge triangles, determine that step-length is

S503: determine song according to the direction of described boundary edge, described step-length, the normal direction of described edge triangles The direction that face extends.

Further, in the another embodiment that the application provides, according to the direction of described boundary edge, described Step-length, the normal direction of described edge triangles determine the direction of surface extending, specifically include:

Direction according to below equation calculating surface extending:

Wherein, definition triangle T₀Adjacent triangle be and T₀The triangle on certain limit altogether, defines T₀ K to adjoin triangle be T₀Through less than or equal to k the accessible triangle sets of syntopy, it is designated as φ_k, the n of set_k；

α is the affecting parameters of current triangle；

For when the normal vector of front piece；

K for current border face adjoins the normal vector of the i-th triangle in triangle.

S603: calculate terminal according to described starting point, described step-length, described surface extending telegoniometer.

Further, in the another embodiment that the application provides, according to described starting point, described step-length, institute State surface extending telegoniometer and calculate terminal, specifically include:

Terminal is calculated according to below equation:

S703: described terminal is added in the some set of three-dimensional geological surface model.

S803: generate the triangle gridding of topology to carrying out two dimension Delaunay Triangulation in described set.

Further, in the another embodiment that the application provides, carry out two dimension in described set Delaunay Triangulation generates the triangle gridding of topology, specifically includes:

The initial boundary triangle of all sampled points in the some set of one three-dimensional geological surface model of generation；

Select each sampled point one by one, search the triangle sets that circumscribed circle comprises described sampled point；

Delete triangle sets, initial boundary triangle is formed cavity；

Connect described sampled point and each the limit in cavity, form gridding triangle sets；

Delete Mesh triangle sets has the gridding triangle of same vertices with initial boundary triangle.

In the embodiment of the present application, the method for surface extending, for processing due to the restriction of various objective condition And make the problem that sampling point set can not meet surface intersection condition.Concrete, in actual landform, bent Crossing between face, but, due to gather spatial point the dough sheet making to be made up of spatial point between not There is overlapping relation, then use this surface extending method so that surface intersection, thus, it is impossible to directly Carry out curved surface computing with renwing stratum and the shape of tomography.

Refer to Fig. 5, the embodiment of the present application provides the surface intersection data in a kind of three-dimensional geological surface model Processing method, comprises the following steps:

S104: each first tri patch in selected first surface.

S204: each second tri patch in selected second curved surface.

S304: judge whether the first tri patch and the second tri patch intersect.

Further, in the another embodiment that the application provides, it is judged that the first tri patch and the second triangle Whether dough sheet intersects, and specifically includes:

Judge whether the bounding box of the first tri patch and the bounding box of the second tri patch intersect.

Further, in the another embodiment that the application provides, it is judged that the first tri patch and the second triangle Whether dough sheet intersects, and specifically includes:

When the first tri patch isSecond tri patch isTime, According to formulaWhen t ∈ [0,1] then, first Tri patch and the second tri patch intersect.

S404: when the first tri patch and the second tri patch intersect, determine the first tri patch and second The intersection point of tri patch and intersection section.

S504: connect intersection section and generate intersection.

In the embodiment of the present application, surface intersection data processing method, with the intersection in the hope of two curved surfaces, Thus carry out curved surface and split to recover geologic structure pattern.

Refer to Fig. 6, the embodiment of the present application provides the edge filtering data in a kind of three-dimensional geological surface model Processing method, comprises the following steps:

S105: input all triangle sets and all boundary sections set.

The all boundary points of S205: labelling.

S305: the arbitrary triangle in selected described triangle sets.

S405: judge whether described vertex of a triangle is labeled.

S505: when described vertex of a triangle is not labeled, described triangle is included into first surface set, Described triangle is not processed.

S605: the adjacent triangle of traversal first surface set, when described adjacent vertex of a triangle is not marked Note, and when the adjacent side of described adjacent triangle and described triangle is not belonging to boundary sections set, by institute State adjacent triangle and be included into first surface set.

S705: repeat previous step, until obtaining stable first surface set.

Further, in the another embodiment that the application provides, described method also includes:

Travel through all trianglees in described triangle sets；

Obtain the curved surface set for each triangle.

In the embodiment of the present application, first surface is integrated into a spatially corresponding independent space section, from And the quantitative research to geologic structure can be conducive to renwing stratum shape.

Meanwhile, all triangle sets mixed in together are split as some independent space section so that Stratum shape is more accurate.

To sum up, in the embodiment of the present application, refer to Fig. 7, it is provided that in a kind of three-dimensional geological surface model Data processing method, comprise the following steps:

S110: obtain sampled point set；

S210: use interpolation algorithm, carries out interpolation and generates gain-type sampled point set sampled point set, with Recover the three-dimensional point cloud structure of geology；

S310: use delaunay subdivision algorithm, processing gain type sampled point set, generate triangle sets；

S410: search the lopsided triangle in triangle sets；

S510: process deformity triangle；

S610: use windowed sinc interpolation kernel function, optimization processes the triangle after deformity triangle Set；

S710: according to surface extending algorithm and surface intersection algorithm, processes the triangle sets after optimizing, raw Become surface boundary line；

S810: according to edge filtering algorithm, generate the space curved surface set at surface boundary line, in order to characterize Space section.

Further, in the another embodiment that the application provides, use interpolation algorithm, to sampled point set Carry out interpolation and generate gain-type sampled point set, specifically include:

Calculate the curvature of each sampled point in sampled point set；

Sampled point set is projected to horizontal plane；

The control point generating the horizontal plane projection for covering sampled point set at horizontal plane is gathered；

The curvature at the control point of the curvature estimation control point set according to sampled point；

The interpolation density at control point is calculated according to curvature density transformation model；

Interpolation density according to control point determines interpolation point set；

Collect sampled point set and interpolation point set, form described gain-type sampled point set.

Further, in the another embodiment that the application provides, processing gain type sampled point set, generate Triangle sets, specifically includes:

Generate one and comprise the initial boundary triangle of all sampled points in gain-type sampled point set；

Select each sampled point in gain-type sampled point set one by one, search circumscribed circle and comprise described sampled point Triangle sets；

Delete triangle sets, initial boundary triangle is formed cavity；

Connect described sampled point and each the limit in cavity, form gridding triangle sets；

Delete Mesh triangle sets has the gridding triangle of same vertices with initial boundary triangle.

Further, in the another embodiment that the application provides, search the lopsided triangle in triangle sets Shape, specifically includes:

According to formula

Search deformity triangle；

Wherein, a, b, c are three sides of a triangle, and δ is threshold parameter.

Further, in the another embodiment that the application provides, according to surface extending algorithm and surface intersection Algorithm, processes the triangle sets after optimizing, and generates surface boundary line, specifically includes:

Find out all boundary edge set and edge triangles set in surface of stratum triangle gridding；

Selecting the boundary edge in boundary edge set one by one, this boundary edge belongs in an edge triangles set One edge triangles；

Selecting the boundary edge in boundary edge set one by one, this boundary edge belongs in an edge triangles set One edge triangles；

The step-length of extension is determined according to described edge triangles；

Direction according to described boundary edge, described step-length, the normal direction of described edge triangles determine surface extending Direction；

Terminal is calculated according to described starting point, described step-length, described surface extending telegoniometer；

Described terminal is added in the some set of three-dimensional geological surface model；

The triangle gridding of topology is generated to carrying out two dimension Delaunay Triangulation in described set.

Further, in the another embodiment that the application provides, described method also includes:

Each first tri patch in selected first surface；

Each second tri patch in selected second curved surface；

Judge whether the first tri patch and the second tri patch intersect；

When the first tri patch and the second tri patch intersect, determine the first tri patch and the second triangular facet The intersection point of sheet and intersection section；

Connect intersection section and generate intersection.

Further, in the another embodiment that the application provides, according to edge filtering algorithm, generate curved surface Space curved surface set at boundary line, specifically includes:

Input all triangle sets and all boundary sections set；

The all boundary points of labelling；

Arbitrary triangle in selected described triangle sets；

Judge whether described vertex of a triangle is labeled；

When described vertex of a triangle is not labeled, described triangle is included into first surface set, otherwise Described triangle is not processed；

The adjacent triangle of traversal first surface set, when described adjacent vertex of a triangle is not labeled, and And the adjacent side of described adjacent triangle and described triangle is when being not belonging to boundary sections set, will described adjoin Triangle is included into first surface set；

Repeat previous step, until obtaining stable first surface set.

The foregoing is only embodiments herein, be not limited to the application.For this area skill For art personnel, the application can have various modifications and variations.All institutes within spirit herein and principle Any modification, equivalent substitution and improvement etc. made, within the scope of should be included in claims hereof.

Claims (7)

1. the data processing method in a three-dimensional geological surface model, it is characterised in that include following step Rapid:

Obtain sampled point set；

Use interpolation algorithm, sampled point set is carried out interpolation and generates gain-type sampled point set, to recover ground The three-dimensional point cloud structure of matter；

Use delaunay subdivision algorithm, processing gain type sampled point set, generate triangle sets；

Search the lopsided triangle in triangle sets；

Process deformity triangle；

Using windowed sinc interpolation kernel function, optimization processes the triangle sets after deformity triangle；

According to surface extending algorithm and surface intersection algorithm, process the triangle sets after optimizing, generate curved surface Boundary line；

According to edge filtering algorithm, generate the space curved surface set at surface boundary line, break in order to characterize space Face.

2. the method for claim 1, it is characterised in that use interpolation algorithm, to sampling point set Conjunction carries out interpolation and generates gain-type sampled point set, specifically includes:

Calculate the curvature of each sampled point in sampled point set；

Sampled point set is projected to horizontal plane；

The control point generating the horizontal plane projection for covering sampled point set at horizontal plane is gathered；

The curvature at the control point of the curvature estimation control point set according to sampled point；

The interpolation density at control point is calculated according to curvature density transformation model；

Interpolation density according to control point determines interpolation point set；

Collect sampled point set and interpolation point set, form described gain-type sampled point set.

3. the method for claim 1, it is characterised in that processing gain type sampled point set, raw Triangularity set, specifically includes:

Generate one and comprise the initial boundary triangle of all sampled points in gain-type sampled point set；

Select each sampled point in gain-type sampled point set one by one, search circumscribed circle and comprise described sampled point Triangle sets；

Delete triangle sets, initial boundary triangle is formed cavity；

Connect described sampled point and each the limit in cavity, form gridding triangle sets；

Delete Mesh triangle sets has the gridding triangle of same vertices with initial boundary triangle.

4. the method for claim 1, it is characterised in that search the deformity three in triangle sets Dihedral, specifically includes:

According to formula

Search deformity triangle；

Wherein, a, b, c are three sides of a triangle, and δ is threshold parameter.

5. the method for claim 1, it is characterised in that ask according to surface extending algorithm and curved surface Hand over algorithm, process the triangle sets after optimizing, generate surface boundary line, specifically include:

Find out all boundary edge set and edge triangles set in surface of stratum triangle gridding；

Selecting the boundary edge in boundary edge set one by one, this boundary edge belongs in an edge triangles set One edge triangles；

Selecting the boundary edge in boundary edge set one by one, this boundary edge belongs in an edge triangles set One edge triangles；

The step-length of extension is determined according to described edge triangles；

Direction according to described boundary edge, described step-length, the normal direction of described edge triangles determine surface extending Direction；

Terminal is calculated according to described starting point, described step-length, described surface extending telegoniometer；

Described terminal is added in the some set of three-dimensional geological surface model；

The triangle gridding of topology is generated to carrying out two dimension Delaunay Triangulation in described set.

6. method as claimed in claim 5, it is characterised in that described method also includes:

Each first tri patch in selected first surface；

Each second tri patch in selected second curved surface；

Judge whether the first tri patch and the second tri patch intersect；

When the first tri patch and the second tri patch intersect, determine the first tri patch and the second triangular facet The intersection point of sheet and intersection section；

Connect intersection section and generate intersection.

7. the method for claim 1, it is characterised in that according to edge filtering algorithm, generates song Space curved surface set at boundary line, face, specifically includes:

Input all triangle sets and all boundary sections set；

The all boundary points of labelling；

Arbitrary triangle in selected described triangle sets；

Judge whether described vertex of a triangle is labeled；

When described vertex of a triangle is not labeled, described triangle is included into first surface set, otherwise Described triangle is not processed；

The adjacent triangle of traversal first surface set, when described adjacent vertex of a triangle is not labeled, and And the adjacent side of described adjacent triangle and described triangle is when being not belonging to boundary sections set, will described adjoin Triangle is included into first surface set；

Repeat previous step, until obtaining stable first surface set.