CN109584357B - Three-dimensional modeling method, device and system based on multiple contour lines and storage medium - Google Patents
Three-dimensional modeling method, device and system based on multiple contour lines and storage medium Download PDFInfo
- Publication number
- CN109584357B CN109584357B CN201811421209.XA CN201811421209A CN109584357B CN 109584357 B CN109584357 B CN 109584357B CN 201811421209 A CN201811421209 A CN 201811421209A CN 109584357 B CN109584357 B CN 109584357B
- Authority
- CN
- China
- Prior art keywords
- constraint
- distance
- interpolation
- point
- sampling
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
Abstract
The invention discloses a three-dimensional modeling method, a device, a system and a storage medium based on multiple contour lines, wherein the method comprises the following steps: discretizing and sampling the plurality of contour lines to obtain a first constraint condition and a second constraint condition; constructing a signed distance field for characterizing the initialized three-dimensional model according to the first constraint and the second constraint; modifying the signed distance field to obtain a distance constraint value for a constraint point for spatial interpolation; carrying out interpolation calculation on the plurality of contour lines by utilizing an interpolation function according to the distance constraint value to obtain an implicit function distance field for representing the idealized three-dimensional model; and sampling and calculating the implicit function distance field based on a reference isosurface extraction algorithm to obtain a three-dimensional model. The invention effectively improves the intellectualization and self-adaptive level of three-dimensional modeling.
Description
Technical Field
The invention relates to the field of three-dimensional modeling, in particular to a three-dimensional modeling method, a three-dimensional modeling device, a three-dimensional modeling system and a storage medium based on multiple contour lines.
Background
Three-dimensional models have been widely used in the fields of geological exploration, medical inspection, movie and television production, and in order to meet the requirements for system analysis of various inspection data and image data, three-dimensional models are often generated by adopting a three-dimensional modeling technology according to the contour lines of structures. For example, in the field of three-dimensional geological modeling, a geological engineer may need to outline an ore body for a region of interest on an exploratory line profile based on borehole sampling information and reconstruct a three-dimensional ore body surface model from a series of intersecting ore body contours.
In the prior art, contour line splicing is mainly adopted for contour line-based three-dimensional reconstruction, and various contour line splicing algorithms such as a maximum volume method, a minimum surface area method, a side length minimum method, a shortest diagonal method, a synchronous advancement method, an incision stitching method and the like are generated in order to adapt to various complex situations.
In the process of implementing the technical scheme of the embodiment of the application, the inventor finds that the algorithms have different adaptability, but does not have any theoretical guidance on the adaptability of the algorithms, and usually needs to perform manual interactive modeling in a heuristic mode in actual operation, so that the existing three-dimensional reconstruction method based on the contour line is not intelligent enough and cannot meet the adaptability of various complex scenes.
Disclosure of Invention
In view of this, embodiments of the present invention provide a method, an apparatus, a system and a storage medium for three-dimensional modeling based on multiple contour lines, and aim to improve the intelligence and adaptive level of three-dimensional modeling.
The technical scheme of the embodiment of the invention is realized as follows:
according to a first aspect of the embodiments of the present invention, there is provided a multi-contour-line-based three-dimensional modeling method, including:
discretizing and sampling a plurality of contour lines to obtain a first constraint condition and a second constraint condition, wherein the first constraint condition comprises sampling points for representing the boundary of each contour line real model and is used for limiting the constraint of the contour lines on the three-dimensional model; the second constraint condition comprises internal sampling points and/or external sampling points which are positioned inside the closed contour line formed by the contour lines and used for controlling the internal and external properties and the extrapolation tendency of the three-dimensional model;
constructing a signed distance field for characterizing the initialized three-dimensional model according to the first constraint and the second constraint;
modifying the signed distance field to obtain a distance constraint value for a constraint point for spatial interpolation;
carrying out interpolation calculation on the plurality of contour lines by utilizing an interpolation function according to the distance constraint value to obtain an implicit function distance field for representing the idealized three-dimensional model;
and sampling and calculating the implicit function distance field based on a reference isosurface extraction algorithm to obtain a three-dimensional model.
According to a second aspect of embodiments of the present invention, there is provided a multi-contour-line-based three-dimensional modeling apparatus, the apparatus including:
the sampling module is used for carrying out discretization sampling on a plurality of contour lines to obtain a first constraint condition and a second constraint condition, wherein the first constraint condition comprises sampling points for representing the boundary of each contour line real model and is used for limiting the constraint of the contour lines on the three-dimensional model; the second constraint condition comprises internal sampling points and/or external sampling points which are positioned inside the closed contour line formed by the contour lines and used for controlling the internal and external properties and the extrapolation tendency of the three-dimensional model;
a first distance field construction module that constructs a signed distance field characterizing the initialized three-dimensional model according to the first constraint and the second constraint;
a correction module for correcting the signed distance field to obtain a distance constraint value of a constraint point for spatial interpolation;
a second distance field construction module, configured to perform interpolation calculation on the plurality of contour lines by using an interpolation function according to the distance constraint value to obtain an implicit function distance field for representing the idealized three-dimensional model;
and the model construction module is used for carrying out sampling calculation on the implicit function distance field based on a reference isosurface extraction algorithm to obtain a three-dimensional model.
According to a third aspect of the embodiments of the present invention, there is provided a multi-contour-line-based three-dimensional modeling system, including a memory, a processor, and an executable program stored on the memory and capable of being executed by the processor, where the processor executes the executable program to perform the multi-contour-line-based three-dimensional modeling method according to any one of the above embodiments.
According to a fourth aspect of the embodiments of the present invention, a storage medium having an executable program stored thereon is characterized in that the executable program, when executed by a processor, implements the multi-contour-based three-dimensional modeling method according to any one of the foregoing embodiments.
According to the three-dimensional modeling method, the three-dimensional modeling device, the three-dimensional modeling system and the storage medium based on the multi-contour line, the distance constraint value of the constraint point for spatial interpolation is obtained by correcting the signed distance field, the implicit function distance field for representing the idealized three-dimensional model is obtained by performing interpolation calculation on the plurality of contour lines according to the distance constraint value through the interpolation function, the three-dimensional model is obtained by performing sampling calculation on the implicit function distance field based on the reference isosurface extraction algorithm, three-dimensional modeling can be performed according to the implicit function distance field obtained after the interpolation calculation, a large number of manual interaction operation processes in the modeling process are avoided, and the intelligence and self-adaption level of the three-dimensional modeling are improved.
Drawings
FIG. 1 is a schematic flow chart of a multi-contour-based three-dimensional modeling method according to an embodiment of the present invention;
FIG. 2 is a schematic structural diagram of adaptive sampling according to an embodiment of the present invention;
FIG. 3 is a schematic diagram of distance field correction according to an embodiment of the invention;
FIG. 4 is a second schematic diagram of distance field correction according to an embodiment of the invention;
FIG. 5 is a schematic flow chart of an iterative closest point correction algorithm according to an embodiment of the present invention;
FIG. 6 is a diagram illustrating one of the effects of parallel contour modeling according to an embodiment of the present invention;
FIG. 7 is a second diagram illustrating the effect of parallel contour modeling according to an embodiment of the present invention;
FIG. 8 is one of the effects of cross-contour modeling using an embodiment of the present invention;
FIG. 9 is a second schematic diagram illustrating the effect of cross-contour modeling according to an embodiment of the present invention;
FIG. 10 is a schematic structural diagram of a multi-contour-based three-dimensional modeling apparatus according to an embodiment of the present invention;
FIG. 11 is a schematic structural diagram of a multi-contour-line-based three-dimensional modeling system according to an embodiment of the present invention.
Detailed Description
The technical scheme of the invention is further elaborated by combining the drawings and the specific embodiments in the specification. It should be understood that the examples provided herein are merely illustrative of the present invention and are not intended to limit the present invention. In addition, the following embodiments are provided as partial embodiments for implementing the present invention, not all embodiments for implementing the present invention, and the technical solutions described in the embodiments of the present invention may be implemented in any combination without conflict.
Unless defined otherwise, 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 invention belongs. The terminology used herein in the description of the invention is for the purpose of describing particular embodiments only and is not intended to be limiting of the application.
It should be noted that, in the embodiments of the present invention, the terms "comprises", "comprising" or any other variation thereof are intended to cover a non-exclusive inclusion, so that a method or apparatus including a series of elements includes not only the explicitly recited elements but also other elements not explicitly listed or inherent to the method or apparatus. Without further limitation, an element defined by the phrase "comprising an … …" does not exclude the presence of other related elements in a method or apparatus including the element (e.g., steps in a method or elements in an apparatus, such as a unit may be part of a circuit, part of a processor, part of a program or software, etc.).
Fig. 1 is a schematic flow chart of a multi-contour-based three-dimensional modeling method according to an embodiment of the present invention. Referring to fig. 1, an embodiment of the present invention provides a three-dimensional modeling method based on multiple contour lines, including:
Here, contour lines refer to a plurality of segments defining the boundaries of a model, for example, in the field of mineral exploration, and a plurality of contour lines refer to the contours of a mineral body drawn by a geological engineer from borehole sampling information over the section of the mineral body under exploration. Discretizing and sampling each contour line according to set sampling precision to form a first constraint condition comprising a plurality of boundary sampling points, wherein the first constraint condition is a hard constraint for limiting the boundary of the contour line of the real three-dimensional model. In this embodiment, the hard constraint is a deterministic point constraint that must be interpolated accurately, such as a point constraint sampled on a contour with a determined distance value.
In this embodiment, in order to perform interpolation constraint on the contours, a soft constraint is introduced, which refers to an incrementally sampled point constraint with approximate distance values, and this constraint is used for distance correction in constructing an implicit function distance field of an idealized three-dimensional model. In an alternative embodiment, the soft constraints include constraints on the incremental sample points inside and outside the closed contour (such as the second constraint below). In another alternative embodiment, the soft constraint may further add an additional boundary constraint, which is manually specified, to the second constraint for controlling the extrapolation tendency of the three-dimensional model.
In this embodiment, discretizing the plurality of contour lines to obtain a second constraint condition includes:
carrying out coordinate transformation on each contour line by taking the corresponding tangent plane as a projection plane to obtain two-dimensional data corresponding to the closed contour line;
carrying out self-adaptive sampling on the closed contour line according to the set sampling precision to obtain sampling data;
and constructing an internal sampling point and/or an external sampling point according to the sampling data.
In an optional implementation manner, a tangent plane where each contour line is located is used as a projection plane to perform coordinate transformation, the projection plane is converted into two-dimensional data corresponding to a closed contour line, and then a quadtree structure is used to perform two-dimensional adaptive sampling on the closed contour line until the closed contour line is subdivided to a pre-specified subdivision level or subdivision size. The pre-specified number of subdivision levels or subdivision sizes may be set or adjusted according to the sampling accuracy requirements. And then, constructing an internal sampling point inside the closed contour line and/or an external sampling point outside the closed contour line by taking the self-adaptive sampled quad-tree unit as a center. In the actual modeling process, a plurality of contour lines can be defined according to a plurality of parallel flat sections, and at the moment, the plurality of contour lines are subjected to coordinate transformation through a projection plane to form a closed contour line. In this case, the plurality of contour lines may be formed into a plurality of closed contour lines by coordinate conversion using a plurality of projection planes. Constructing internal sampling points inside the closed contour line by taking the self-adaptive sampled quad-tree unit as a center; or constructing an external sampling point outside the closed contour line by taking the self-adaptive sampled quad-tree unit as the center; and constructing an internal sampling point inside the closed contour line and an external sampling point outside the closed contour line by taking the self-adaptive sampled quad-tree unit as a center. Fig. 2 is a schematic structural diagram of adaptive sampling in the embodiment of the present invention, which shows a schematic diagram of constructing internal sampling points inside a closed contour line.
In this embodiment, a Signed Distance Function (SDF) field may be constructed by calculating the closest Distance from an internal sample point to a closed contour; the signed distance field may also be constructed by calculating the closest distance from the external sample point to the closed contour line, or by calculating the closest distance from the internal sample point and the external sample point to the closed contour line. The signed distance field of this embodiment forms a spatial constraint using the boundary sample points, the internal sample points, and/or the external sample points, and the relationship between the function values f (x) of the signed distance field and the model boundary can be expressed as:
wherein X is any point in three-dimensional space, R3And representing a three-dimensional space, wherein the function value of the internal sampling point is a negative value, the function value of the external sampling point is a positive value, and a zero level set extracted according to the condition that the function value is zero is the boundary of the model.
In order to avoid repeated, abnormal and ambiguous distance function values, the embodiment verifies the distance function values of the discrete points, and through a distance field correction process, the closest distance from a point to a corresponding tangent plane contour line can be improved to be used for estimating the closest distance from a sampling point to a curved surface.
In this embodiment, the modifying the signed distance field to obtain a distance constraint value for a constraint point used for spatial interpolation includes: and correcting the constraint points inside or outside the closed contour line based on the definition of Euclidean space distance and the principle of minimum distance between the constraint points and the closed contour line to obtain a distance constraint value of the constraint points for spatial interpolation.
In distance fields, the signed implicit function s (x) should satisfy the principle of minimum boundary distance and the definition of distance in Euclidean space. The triangle inequality condition defined by distance d (x, z) + d (z, y) ≧ d (x, y) is combined with the definition of the signed distance field
A and B are arbitrary interpolation central points in the spatial domain, and the positive and negative signs corresponding to the function values are respectively SignAAnd SignBA 'and B' are A, B, respectively, the closest point to the curved surface, d (A, A ') is the Euclidean distance between A and A', and d (B, B ') is the Euclidean distance between B and B'. According to s (x)A) And s (x)B) The difference in symbols can be classified into the following two cases:
referring to FIG. 3, when s (x)A) And s (x)B) The symbols are the same, have
Referring to FIG. 4, when s (x)A) And s (x)B) When the symbols are opposite, there are
Therefore, the above formula should be satisfied for each out-of-plane constraint, otherwise, a correction should be made. Taking the case of the same sign as an example, if the soft constraint occurs, | s (x)A)|＞|s(xB) If | d (A, B), the distance value of A can be | s (x)A)|＝|s(xB) L + d (A, B). Soft constraints and hard constraints are considered simultaneously when the distance field is corrected, and the addition of the hard constraints is beneficial to better evaluating the distribution trend of the distance field. In this embodiment, the distance field modification can iteratively modify the distance values for soft constraints, which are used for distance modification calculations (and do not participate in distance modification).
And 104, carrying out interpolation calculation on the plurality of contour lines by utilizing an interpolation function according to the distance constraint value to obtain an implicit function distance field for representing the idealized three-dimensional model.
In an optional embodiment, interpolating the plurality of contour lines by an interpolation function according to the distance constraint value to obtain an implicit function distance field for characterizing the idealized three-dimensional model includes:
determining the interpolation function according to the constraint points which accord with the distance constraint value;
carrying out interpolation calculation on the plurality of contour lines according to the interpolation function to obtain interpolation constraint points;
and updating the signed distance field according to the interpolation constraint point to obtain the implicit function distance field.
Here, determining the interpolation function from the constraint points that conform to the distance constraint value includes: solving a large linear equation set, obtaining an unknown coefficient of an interpolation function, and establishing an interpolation function formula as follows:
wherein x ═ (x, y, z) is any interpolation constraint point to be solved; f (x) is the radial basis function (abbreviated RBF); p (x) is a low order polynomial, typically a low order polynomial; lambda [ alpha ]iCoefficients that are radial basis functions; phi is a strictly positive or conditionally positive kernel function; c. CiIs the interpolated center point of the radial basis function.
Optionally, the plurality of contour lines are interpolated according to the interpolation function to obtain interpolation constraint points, and the interpolation data may be divided into near-field data and far-field data according to a fast multi-polarization method according to a far-near hierarchy from an interpolation center x equal to (x, y, z). The near field data is solved by a direct method, and the far field data is solved by an approximate method within a certain error range, so that the speed of interpolation calculation can be greatly increased, high-efficiency interpolation and curved surface reconstruction can be performed on large-scale constraint data, and the problem of efficiency of curved surface reconstruction after incremental sampling of the contour line is solved.
In another optional embodiment, the interpolating the plurality of contour lines according to the distance constraint value by using an interpolation function to obtain an implicit function distance field for characterizing the idealized three-dimensional model includes:
determining the interpolation function according to the constraint points which accord with the distance constraint value;
carrying out interpolation calculation on the plurality of contour lines according to the interpolation function to obtain interpolation constraint points;
constructing additional boundary constraints, wherein the additional boundary constraints are used for controlling the extrapolation trend of the three-dimensional model;
and updating the signed distance field according to the interpolation constraint point and the additional boundary constraint to obtain the implicit function distance field.
For sparse contour data, local adjustments to the automated modeling process may be made, such as creating local ore model constraint lines according to particular geological rules, forming additional boundary constraints, and updating the signed distance field according to the interpolation constraint points and the additional boundary constraints to obtain the implicit function distance field.
In order to construct the additional boundary constraint, a discretization method can be adopted to convert the interpolation constraint line into an interpolation constraint point. And constructing normal information by limiting the tangent plane direction of the constraint line by adopting the constraint line of the additional normal constraint section. The tangent plane of the constraint line defines the plane where the normal direction of the constraint line is located, and the normal direction information of the constraint line is determined in sections according to the section directions of the constraint line; and determining the normal direction (pointing to the outside of the model) by combining the internal and external relations of the model cut by the tangent plane. Discretizing and sampling the interpolation constraint lines according to a certain interval, and converting the additional constraint lines into interpolation constraint points so as to update and generate an implicit function distance field representing the three-dimensional model.
And 105, performing sampling calculation on the implicit function distance field based on a reference isosurface extraction algorithm to obtain a three-dimensional model.
In this embodiment, the performing sampling computation on the implicit function distance field based on a reference isosurface extraction algorithm to obtain a three-dimensional model includes:
extracting a minimum outer package for defining a model boundary according to the interpolation constraint point;
determining the resolution of three-dimensional reconstruction of the model according to the minimum outsourcing;
carrying out space segmentation on the minimum outer package according to the resolution ratio to construct a discretized space rule data field;
and carrying out sampling calculation on the implicit function distance field according to the discretized space rule data field, and reconstructing the three-dimensional model.
Here, the minimum outsourcing refers to a solid frame surrounding the three-dimensional model, and the initial resolution of the three-dimensional reconstruction of the ore body model is determined according to 0.1 times of the minimum side length of the minimum outsourcing in length, width and height. Optionally, the resolution may also be adjusted according to actual needs. And carrying out space segmentation on the minimum outsourcing according to the selected size of the cubic unit to construct a discretized space rule data field. And (3) sampling and calculating the implicit function distance field according to the mobile cube by adopting a voxel growing method, and reconstructing a three-dimensional model.
In this embodiment, by constructing the initial voxel seed point and tracking the curved surface reconstruction process by using a certain voxel growth rule, the evaluation of the function value by the radial basis function on all voxel grid points can be avoided, and the model reconstruction process is accelerated.
In the three-dimensional modeling method based on multiple contour lines, the distance constraint value of the constraint point for spatial interpolation is obtained by correcting the signed distance field, the implicit function distance field for representing the idealized three-dimensional model is obtained by performing interpolation calculation on the plurality of contour lines according to the distance constraint value by using the interpolation function, the three-dimensional model is obtained by performing sampling calculation on the implicit function distance field based on the reference isosurface extraction algorithm, and a large number of manual interaction operation processes in the modeling process are avoided, so that the intellectualization and the adaptability level of the three-dimensional modeling are improved.
In an optional implementation manner, referring to fig. 5, the constraint points inside or outside the closed contour line are modified based on the definition of the euclidean spatial distance and the principle that the distance from the constraint points to the closed contour line is the minimum, so as to obtain the distance constraint value of the constraint point used for spatial interpolation, and an iterative closest point modification algorithm is adopted. The iterative closest point correction algorithm comprises the following steps:
and step 501, sorting constraint points inside or outside the closed contour line from small to large according to the distance value between the constraint points and the closed contour line to obtain a data set.
Here, since the iterative process is a process of gradually reducing the distance value, in order to avoid a large number of repeated comparison processes, each soft constraint point may be sorted from small to large according to the distance value, and the data set is S.
Constructing a current minimum distance constraint value dminAs a starting value of the distance constraint value, an initial value thereof may be set to zero.
Here, d is never smaller thanminSelects a constraint point from the subset of (a) as a first constraint point xAFinding the first constraint point x by adopting a kd-tree spatial index algorithmAThe closest constraint point is taken as a second constraint point x for distance comparisonB。
Determining the first constraint point x according to a distance determination formulaAChecking to determine | s (x)A)|＞|s(xB) If | + d (a, B), it is determined that the correction condition is satisfied, and if the correction condition is satisfied, step 505 is executed, and if the correction condition is not satisfied, step 507 is executed.
For the first constraint point xAIs corrected so that | s (x)A)|＝|s(xB) L + d (a, B), but without changing the sign of its distance value. A, B is an arbitrary interpolation center point, | s (x) in the spatial domainA) L is xAIs given by the distance value, | s (x)B) L is xBD (a, B) is the euclidean distance between a and B.
Reordering the data set S from small to large according to the distance value, and according to the updated first constraint point xAUpdate the distance constraint value dminAnd returns to step 503.
And step 507, ending the iteration.
And when the first constraint point is judged not to meet the correction condition, judging that the iteration process is converged, and ending the iteration closest point correction algorithm.
The three-dimensional modeling method based on the multiple contour lines can be applied to the field of ore body exploration, can automatically model according to the parallel contour lines of the ore body flat section, and can also automatically model contour lines of a mixed flat section comprising the ore body flat section and the ore body longitudinal section. The embodiment avoids a large number of manual interactive operation processes in modeling. By using an iterative closest point correction algorithm, iterative correction of the signed distance field can be performed based on constraints and model internal and external relationships. In addition, the efficiency of the algorithm on the interpolation and reconstruction process can be improved by adopting a rapid multi-polarization method and a curved surface tracking method, the high-efficiency interpolation and curved surface reconstruction can be carried out on large-scale constraint data, and the problem of the efficiency of curved surface reconstruction after incremental sampling of the contour line is solved.
It should be noted that the three-dimensional modeling method based on multiple contour lines in this embodiment may be applied to other similar contour line modeling fields, such as computer graphics, medical CT slice modeling, and the like, in addition to the ore body contour line modeling field of geological modeling.
In order to test the universality of the multi-contour-line-based three-dimensional modeling method of the present embodiment, fig. 6 to 9 show. The results show that for dense parallel sections (as shown in fig. 6 and 7), the results of sampling the contours reconstruct the original model well without adding additional interpolation constraints. For sparse cross-sections (as shown in fig. 8 and 9), additional boundary constraints may be interactively added to limit model boundaries in order to control the interpolated extrapolated boundaries. For non-parallel cross sections, on one hand, the contour lines have better constraint action on the internal and external relations of the model, and can better restore the original model; on the other hand, ambiguity constraint is easy to occur in sampling between intersecting contour lines, and the reconstructed model generates abnormity such as distortion and recess. In the embodiment, the ambiguity problem of the soft constraint distance value between the cross sections can be effectively solved by adding the additional interpolation constraint condition, and the defects are avoided, so that the modeling quality is improved, and the subsequent engineering application and analysis are facilitated.
Referring to fig. 10, an embodiment of the present invention further provides a three-dimensional modeling apparatus based on multiple contour lines, including:
the sampling module 110 is configured to perform discretization sampling on a plurality of contour lines to obtain a first constraint condition and a second constraint condition, where the first constraint condition includes a sampling point for representing a boundary of each contour line real model, and is used to limit constraint of the contour lines on the three-dimensional model; the second constraint condition comprises internal sampling points and/or external sampling points which are positioned inside the closed contour line formed by the contour lines and used for controlling the internal and external properties and the extrapolation tendency of the three-dimensional model;
a first distance field construction module 120 that constructs a signed distance field characterizing the initialized three-dimensional model according to the first constraint and the second constraint;
a correction module 130 for correcting the signed distance field to obtain a distance constraint value for a constraint point for spatial interpolation;
a second distance field construction module 140, configured to perform interpolation calculation on the plurality of contour lines by using an interpolation function according to the distance constraint value to obtain an implicit function distance field for representing the idealized three-dimensional model;
and the model construction module 150 is used for performing sampling calculation on the implicit function distance field based on a reference isosurface extraction algorithm to obtain a three-dimensional model.
In this embodiment, the sampling module 110 is specifically configured to perform discretization sampling on each contour line according to a set sampling precision, so as to form a first constraint condition including a plurality of boundary sampling points.
The sampling module 110 performs discretization sampling on the plurality of contour lines to obtain a second constraint condition, which specifically includes:
carrying out coordinate transformation on each contour line by taking the corresponding tangent plane as a projection plane to obtain two-dimensional data corresponding to the closed contour line;
carrying out self-adaptive sampling on the closed contour line according to the set sampling precision to obtain sampling data;
and constructing an internal sampling point and/or an external sampling point according to the sampling data.
In this embodiment, the correcting module 130 is configured to correct the constraint point inside or outside the closed contour based on the definition of the euclidean spatial distance and the rule that the distance from the constraint point to the closed contour is the minimum, so as to obtain a distance constraint value of the constraint point for spatial interpolation.
In one embodiment, the second distance field construction module 140 is specifically configured to:
determining the interpolation function according to the constraint points which accord with the distance constraint value;
carrying out interpolation calculation on the plurality of contour lines according to the interpolation function to obtain interpolation constraint points;
and updating the signed distance field according to the interpolation constraint point to obtain the implicit function distance field.
In another embodiment, the second distance field construction module 140 is specifically configured to:
determining the interpolation function according to the constraint points which accord with the distance constraint value;
carrying out interpolation calculation on the plurality of contour lines according to the interpolation function to obtain interpolation constraint points;
constructing additional boundary constraints, wherein the additional boundary constraints are used for controlling the extrapolation trend of the three-dimensional model;
and updating the signed distance field according to the interpolation constraint point and the additional boundary constraint to obtain the implicit function distance field.
In one embodiment, the model building module 150 is specifically configured to:
extracting a minimum outer package for defining a model boundary according to the interpolation constraint point;
determining the resolution of three-dimensional reconstruction of the model according to the minimum outsourcing;
carrying out space segmentation on the minimum outer package according to the resolution ratio to construct a discretized space rule data field;
and carrying out sampling calculation on the implicit function distance field according to the discretized space rule data field, and reconstructing the three-dimensional model.
It should be noted that: in the three-dimensional modeling apparatus based on multiple contour lines provided in the above embodiments, only the division of the program modules is illustrated, and in practical applications, the processing may be distributed to different program modules according to needs, that is, the internal structure of the three-dimensional modeling apparatus may be divided into different program modules to complete all or part of the processing described above. In addition, the multi-contour-line-based three-dimensional modeling apparatus provided in the above embodiments and the multi-contour-line-based three-dimensional modeling method embodiments belong to the same concept, and specific implementation processes thereof are described in the method embodiments and are not described herein again.
In practical applications, each of the program modules may be implemented by a Central Processing Unit (CPU) on the server, a microprocessor Unit (MPU), a Digital Signal Processor (DSP), a Field Programmable Gate Array (FPGA), or the like.
In order to realize the three-dimensional modeling method based on the multiple contour lines, the embodiment of the invention also provides a hardware structure of the three-dimensional modeling system based on the multiple contour lines. A multi-profile based three-dimensional modeling system that may be implemented in the form of various types of processors, such as clients, servers, etc., that implement the present invention will now be described with reference to the accompanying drawings. In the following, the hardware structure of the multi-contour based three-dimensional modeling system according to the embodiment of the present invention is further described, it is to be understood that fig. 11 only shows an exemplary structure of the multi-contour based three-dimensional modeling system, and not a whole structure, and a part of the structure or the whole structure shown in fig. 11 may be implemented as needed.
Referring to fig. 11, fig. 11 is a schematic diagram of a hardware structure of a multi-contour-based three-dimensional modeling system according to an embodiment of the present invention, which may be applied to the aforementioned server running an application program in practical applications, where the multi-contour-based three-dimensional modeling system 1100 shown in fig. 11 includes: at least one processor 1101, memory 1102, a user interface 1103, and at least one network interface 1104. The various components in the three-dimensional modeling system 1100 are coupled together by a bus system 1105. It will be appreciated that the bus system 1105 is used to enable communications among the components. The bus system 1105 may include a power bus, a control bus, and a status signal bus in addition to a data bus. For clarity of illustration, however, the various buses are labeled in fig. 11 as the bus system 1105.
The user interface 1103 may include, among other things, a display, a keyboard, a mouse, a trackball, a click wheel, keys, buttons, a touch pad, or a touch screen.
It will be appreciated that the memory 1102 can be either volatile memory or nonvolatile memory, or can include both volatile and nonvolatile memory.
The memory 1102 in embodiments of the present invention is used to store various types of data to support the operation of the multi-profile based three-dimensional modeling system 1100. Examples of such data include: any computer program for operating on a multi-profile based three-dimensional modeling system, such as executable program 11021 and operating system 11022, a program implementing the multi-profile based three-dimensional modeling method of an embodiment of the present invention may be included in executable program 11021.
The multi-contour-line-based three-dimensional modeling method disclosed by the embodiment of the invention can be applied to the processor 1101 or realized by the processor 1101. The processor 1101 may be an integrated circuit chip having signal processing capabilities. In implementation, the steps of the above-described multi-contour-line-based three-dimensional modeling method may be implemented by integrated logic circuits of hardware or instructions in the form of software in the processor 1101. The processor 1101 described above may be a general purpose processor, a DSP, or other programmable logic device, discrete gate or transistor logic device, discrete hardware components, or the like. The processor 1101 may implement or perform the multi-contour line based three-dimensional modeling methods, steps and logic blocks provided in the embodiments of the present invention. A general purpose processor may be a microprocessor or any conventional processor or the like. The steps of the multi-contour-line-based three-dimensional modeling method provided by the embodiment of the invention can be directly implemented by a hardware decoding processor or implemented by combining hardware and software modules in the decoding processor. The software modules may be located in a storage medium located in the memory 1102, and the processor 1101 reads the information in the memory 1102, and completes the steps of the multi-contour line based three-dimensional modeling method provided by the embodiment of the present invention in combination with the hardware thereof.
In this embodiment, the multi-contour-line-based three-dimensional modeling system 1100 includes a memory 1102, a processor 1101, and an executable program 11021 stored in the memory 1102 and capable of being executed by the processor 1101, and the processor 1101 executes the executable program 11021 to implement the multi-contour-line-based three-dimensional modeling method according to the foregoing embodiment.
In an exemplary embodiment, an embodiment of the present invention further provides a storage medium, which may be a storage medium such as a removable storage device, a Read Only Memory (ROM), an optical disc, a flash Memory, or a magnetic disc, and may be selected as a non-transitory storage medium. The storage medium stores an executable program 11021, and when the executable program 11021 is executed by a processor, the multi-contour-line-based three-dimensional modeling method of the foregoing embodiment is realized.
As will be appreciated by one skilled in the art, embodiments of the present invention may be provided as a method, system, or executable program product. Accordingly, the present invention may take the form of a hardware embodiment, a software embodiment, or an embodiment combining software and hardware aspects. Furthermore, the present invention may take the form of an executable program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, optical storage, and the like) having computer-usable program code embodied therein.
The present invention has been described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems), and executable program products according to embodiments of the invention. It will be understood that each flow and/or block of the flow diagrams and/or block diagrams, and combinations of flows and/or blocks in the flow diagrams and/or block diagrams, can be implemented by executable program instructions. These executable program instructions may be provided to a general purpose computer, special purpose computer, embedded processor, or processor with reference to a programmable data processing apparatus to produce a machine, such that the instructions, which execute via the computer or processor with reference to the programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These executable program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These executable program instructions may also be loaded onto a computer or reference programmable data processing apparatus to cause a series of operational steps to be performed on the computer or reference programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or reference programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
The above description is only for the specific embodiments of the present invention, but the scope of the present invention is not limited thereto, and any person skilled in the art can easily conceive of the changes or substitutions within the technical scope of the present invention, and all the changes or substitutions should be covered within the scope of the present invention. Therefore, the protection scope of the present invention shall be subject to the protection scope of the appended claims.
Claims (8)
1. A three-dimensional modeling method based on multiple contour lines is characterized by comprising the following steps:
discretizing and sampling a plurality of contour lines to obtain a first constraint condition and a second constraint condition, wherein the first constraint condition comprises sampling points for representing the boundary of each contour line real model and is used for limiting the constraint of the contour lines on the three-dimensional model; the second constraint condition comprises internal sampling points and/or external sampling points which are positioned inside the closed contour line formed by the contour lines and used for controlling the internal and external properties and the extrapolation tendency of the three-dimensional model;
constructing a signed distance field for characterizing the initialized three-dimensional model according to the first constraint and the second constraint;
modifying the signed distance field to obtain a distance constraint value for a constraint point for spatial interpolation;
carrying out interpolation calculation on the plurality of contour lines by utilizing an interpolation function according to the distance constraint value to obtain an implicit function distance field for representing the idealized three-dimensional model;
sampling and calculating the implicit function distance field based on a reference isosurface extraction algorithm to obtain a three-dimensional model;
the modifying the signed distance field to obtain a distance constraint value for a constraint point for spatial interpolation, comprising:
correcting the constraint points inside or outside the closed contour line based on the definition of Euclidean space distance and the principle of minimum distance between the constraint points and the closed contour line to obtain a distance constraint value of the constraint points for spatial interpolation;
the step of correcting the constraint points inside or outside the closed contour line based on the definition of the Euclidean spatial distance and the principle that the distance between the constraint points and the closed contour line is the minimum to obtain the distance constraint value of the constraint points for spatial interpolation comprises the following steps:
sorting constraint points inside or outside the closed contour line from small to large according to the distance value between the constraint points and the closed contour line to obtain a data set;
initializing the distance constraint value;
traversing the data set, and selecting a constraint point from the subsets not less than the distance constraint value as a first constraint point
Searching for the first constraint point by using a spatial index algorithm
The closest constraint point is taken as the second constraint point of the distance comparison
；
Determining the first constraint point according to a distance determination formula
Checking to determine
For the first constraint point
Is corrected so that
Wherein A, B is an arbitrary interpolation center point in the spatial domain,
is composed of
The value of (a) is determined,
is composed of
The value of (a) is determined,
is the Euclidean distance between A and B;
reordering the data set according to the updated first constraint point
Updating the distance constraint value;
returning the traversal of the data set until the iterative process converges.
2. The multi-contour-line-based three-dimensional modeling method according to claim 1, wherein discretizing the plurality of contour lines to obtain the second constraint comprises:
carrying out coordinate transformation on each contour line by taking the corresponding tangent plane as a projection plane to obtain two-dimensional data corresponding to the closed contour line;
carrying out self-adaptive sampling on the closed contour line according to the set sampling precision to obtain sampling data;
and constructing an internal sampling point and/or an external sampling point according to the sampling data.
3. The multi-contour-line-based three-dimensional modeling method of claim 1, wherein interpolating the plurality of contour lines using an interpolation function based on the distance constraint values to obtain an implicit function distance field for characterizing an idealized three-dimensional model comprises:
determining the interpolation function according to the constraint points which accord with the distance constraint value;
carrying out interpolation calculation on the plurality of contour lines according to the interpolation function to obtain interpolation constraint points;
and updating the signed distance field according to the interpolation constraint point to obtain the implicit function distance field.
4. The multi-contour-line-based three-dimensional modeling method of claim 1, wherein interpolating the plurality of contour lines using an interpolation function based on the distance constraint values to obtain an implicit function distance field for characterizing an idealized three-dimensional model comprises:
determining the interpolation function according to the constraint points which accord with the distance constraint value;
carrying out interpolation calculation on the plurality of contour lines according to the interpolation function to obtain interpolation constraint points;
constructing additional boundary constraints, wherein the additional boundary constraints are used for controlling the extrapolation trend of the three-dimensional model;
and updating the signed distance field according to the interpolation constraint point and the additional boundary constraint to obtain the implicit function distance field.
5. The multi-contour-based three-dimensional modeling method of claim 4, wherein said sampling computation of said implicit function distance field based on a reference iso-surface extraction algorithm to obtain a three-dimensional model comprises:
extracting a minimum outer package for defining a model boundary according to the interpolation constraint point;
determining the resolution of three-dimensional reconstruction of the model according to the minimum outsourcing;
carrying out space segmentation on the minimum outer package according to the resolution ratio to construct a discretized space rule data field;
and carrying out sampling calculation on the implicit function distance field according to the discretized space rule data field, and reconstructing the three-dimensional model.
6. A three-dimensional modeling apparatus based on multiple contour lines, comprising:
the sampling module is used for carrying out discretization sampling on a plurality of contour lines to obtain a first constraint condition and a second constraint condition, wherein the first constraint condition comprises sampling points for representing the boundary of each contour line real model and is used for limiting the constraint of the contour lines on the three-dimensional model; the second constraint condition comprises internal sampling points and/or external sampling points which are positioned inside the closed contour line formed by the contour lines and used for controlling the internal and external properties and the extrapolation tendency of the three-dimensional model;
a first distance field construction module that constructs a signed distance field characterizing the initialized three-dimensional model according to the first constraint and the second constraint;
a correction module for correcting the signed distance field to obtain a distance constraint value for a constraint point for spatial interpolation, comprising: correcting the constraint points inside or outside the closed contour line based on the definition of Euclidean space distance and the principle of minimum distance between the constraint points and the closed contour line to obtain a distance constraint value of the constraint points for spatial interpolation;
the step of correcting the constraint point inside or outside the closed contour line based on the definition of the Euclidean spatial distance and the principle that the distance between the constraint point and the closed contour line is the minimum to obtain the distance constraint value of the constraint point for spatial interpolation comprises the following steps:
sorting constraint points inside or outside the closed contour line from small to large according to the distance value between the constraint points and the closed contour line to obtain a data set;
initializing the distance constraint value;
go throughSelecting a constraint point from the subset not less than the distance constraint value as a first constraint point for the data set
Searching for the first constraint point by using a spatial index algorithm
The closest constraint point is taken as the second constraint point of the distance comparison
；
Determining the first constraint point according to a distance determination formula
Checking to determine
For the first constraint point
Is corrected so that
Wherein A, B is an arbitrary interpolation center point in the spatial domain,
is composed of
The value of (a) is determined,
is composed of
The value of (a) is determined,
is the Euclidean distance between A and B;
reordering the data set according to the updated first constraint point
Updating the distance constraint value;
returning the traversal of the data set until the iterative process converges;
a second distance field construction module, configured to perform interpolation calculation on the plurality of contour lines by using an interpolation function according to the distance constraint value to obtain an implicit function distance field for representing the idealized three-dimensional model;
and the model construction module is used for carrying out sampling calculation on the implicit function distance field based on a reference isosurface extraction algorithm to obtain a three-dimensional model.
7. A multi-profile based three-dimensional modeling system comprising a memory, a processor and an executable program stored on the memory and executable by the processor, wherein the processor executes the executable program to perform the multi-profile based three-dimensional modeling method according to any one of claims 1 to 5.
8. A storage medium having stored thereon an executable program, wherein the executable program, when executed by a processor, implements the multi-profile based three-dimensional modeling method according to any one of claims 1 to 5.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811421209.XA CN109584357B (en) | 2018-11-27 | 2018-11-27 | Three-dimensional modeling method, device and system based on multiple contour lines and storage medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811421209.XA CN109584357B (en) | 2018-11-27 | 2018-11-27 | Three-dimensional modeling method, device and system based on multiple contour lines and storage medium |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109584357A CN109584357A (en) | 2019-04-05 |
| CN109584357B true CN109584357B (en) | 2022-05-17 |
Family
ID=65924658
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811421209.XA Active CN109584357B (en) | 2018-11-27 | 2018-11-27 | Three-dimensional modeling method, device and system based on multiple contour lines and storage medium |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109584357B (en) |
Families Citing this family (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110033519B (en) * | 2019-04-23 | 2022-06-28 | 中南大学 | Three-dimensional modeling method, device and system based on implicit function and storage medium |
| CN111968231B (en) * | 2020-08-14 | 2023-05-30 | 成都理工大学 | Three-dimensional stratum modeling method based on geological map cutting section |
| CN112036030B (en) * | 2020-08-31 | 2021-06-01 | 中南大学 | Ore body combination constraint modeling method, device, equipment and storage medium |
| CN112935575A (en) * | 2021-01-12 | 2021-06-11 | 大族激光科技产业集团股份有限公司 | Cutting path optimization method and device and computer readable storage medium |
| CN112927350B (en) * | 2021-02-23 | 2021-11-02 | 中南大学 | Multi-domain geologic body model construction method, device, equipment and storage medium |
| CN113011040B (en) * | 2021-04-01 | 2022-07-08 | 中南大学 | Pulse geologic body modeling method, device, equipment and storage medium |
| CN113376953B (en) * | 2021-05-20 | 2022-09-27 | 达闼机器人股份有限公司 | Object projection reconstruction system |
| CN113592711B (en) * | 2021-07-22 | 2024-01-12 | 中山大学 | Three-dimensional reconstruction method, system, equipment and storage medium for non-uniform point cloud data |
| CN114357683B (en) * | 2022-01-11 | 2022-08-02 | 中南大学 | Ore body modeling method and device based on cross contour line normal dynamic estimation |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104715508A (en) * | 2015-02-13 | 2015-06-17 | 中南大学 | Grid modeling method and device |
| CN104966317A (en) * | 2015-06-04 | 2015-10-07 | 中南大学 | Automatic three-dimensional modeling method based on contour line of ore body |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6867784B2 (en) * | 2003-04-30 | 2005-03-15 | Landmark Graphics Corporation | Methods, systems, and computer readable media for using closed contours of a modeled surface |
| US9613424B2 (en) * | 2013-09-23 | 2017-04-04 | Beihang University | Method of constructing 3D clothing model based on a single image |
| US10678959B2 (en) * | 2017-03-23 | 2020-06-09 | Autodesk, Inc. | Creating gradients of different materials for three-dimensional models in computer aided design applications |
-
2018
- 2018-11-27 CN CN201811421209.XA patent/CN109584357B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104715508A (en) * | 2015-02-13 | 2015-06-17 | 中南大学 | Grid modeling method and device |
| CN104966317A (en) * | 2015-06-04 | 2015-10-07 | 中南大学 | Automatic three-dimensional modeling method based on contour line of ore body |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109584357A (en) | 2019-04-05 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109584357B (en) | Three-dimensional modeling method, device and system based on multiple contour lines and storage medium | |
| US7443393B2 (en) | Method, system, and program product for re-meshing of a three-dimensional input model using progressive implicit approximating levels | |
| CN104268934B (en) | Method for reconstructing three-dimensional curve face through point cloud | |
| Song et al. | A progressive point cloud simplification algorithm with preserved sharp edge data | |
| Di Angelo et al. | A new mesh-growing algorithm for fast surface reconstruction | |
| CN110033519B (en) | Three-dimensional modeling method, device and system based on implicit function and storage medium | |
| US11704802B2 (en) | Multi-dimensional model merge for style transfer | |
| JP2020115337A (en) | Set of neural networks | |
| Chen et al. | Point cloud resampling using centroidal Voronoi tessellation methods | |
| WO2017203528A1 (en) | Systems and methods for generating volumetric models | |
| WO2013123636A1 (en) | Method and apparatus for mesh simplification | |
| JP2020109661A (en) | Learning of neural network for inference of solid cad features | |
| US20150127301A1 (en) | Updating A CAD Model To Reflect Global Or Local Shape Changes | |
| Liu et al. | Quality improvement of surface triangular mesh using a modified Laplacian smoothing approach avoiding intersection | |
| Daniels et al. | Semi‐regular quadrilateral‐only remeshing from simplified base domains | |
| Tur et al. | Stabilized method of imposing Dirichlet boundary conditions using a recovered stress field | |
| Yu et al. | Geodesics on point clouds | |
| JP2021082297A (en) | Processing 3d signal of shape attribute over real object | |
| Pan et al. | Q-MAT+: An error-controllable and feature-sensitive simplification algorithm for medial axis transform | |
| Foucault et al. | Generalizing the advancing front method to composite surfaces in the context of meshing constraints topology | |
| Li et al. | On surface reconstruction: A priority driven approach | |
| US10943037B2 (en) | Generating a CAD model from a finite element mesh | |
| WO2020023811A1 (en) | 3d object design synthesis and optimization using existing designs | |
| Vančo et al. | Surface reconstruction from unorganized point data with quadrics | |
| Quadros et al. | 3 D discrete skeleton generation by wave propagation on PR-octree for finite element mesh sizing |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |