CN108230452A - A kind of model filling-up hole method based on textures synthesis - Google Patents

A kind of model filling-up hole method based on textures synthesis Download PDF

Info

Publication number
CN108230452A
CN108230452A CN201711364637.9A CN201711364637A CN108230452A CN 108230452 A CN108230452 A CN 108230452A CN 201711364637 A CN201711364637 A CN 201711364637A CN 108230452 A CN108230452 A CN 108230452A
Authority
CN
China
Prior art keywords
model
sample block
vertex
information
grid
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.)
Granted
Application number
CN201711364637.9A
Other languages
Chinese (zh)
Other versions
CN108230452B (en
Inventor
刘斌
徐晓影
许建文
黄常标
江开勇
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Huaqiao University
Original Assignee
Huaqiao University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Huaqiao University filed Critical Huaqiao University
Priority to CN201711364637.9A priority Critical patent/CN108230452B/en
Publication of CN108230452A publication Critical patent/CN108230452A/en
Application granted granted Critical
Publication of CN108230452B publication Critical patent/CN108230452B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T19/00Manipulating 3D models or images for computer graphics
    • G06T19/20Editing of 3D images, e.g. changing shapes or colours, aligning objects or positioning parts
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00Three dimensional [3D] modelling, e.g. data description of 3D objects
    • G06T17/20Finite element generation, e.g. wire-frame surface description, tesselation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T3/00Geometric image transformations in the plane of the image
    • G06T3/06Topological mapping of higher dimensional structures onto lower dimensional surfaces
    • G06T3/067Reshaping or unfolding 3D tree structures onto 2D planes

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Graphics (AREA)
  • Software Systems (AREA)
  • Architecture (AREA)
  • Computer Hardware Design (AREA)
  • General Engineering & Computer Science (AREA)
  • Geometry (AREA)
  • Image Generation (AREA)
  • Processing Or Creating Images (AREA)

Abstract

The invention discloses the model method for repairing and mending based on textures synthesis, and the textured triangle grid model containing hole is selected first, not considered input model the conventional method filling-up hole of texture information, and carry out fairing processing as input;Secondly, characteristic displacement vector is isolated as matched indicia information by the use of multi-level technique;Then parameterized treatment is carried out to the base surface mesh after fairing and obtains two dimensional surface, completed the sampling on each vertex in parametrization plane, obtain the sample set of blocks being made of each vertex sampling area;Then the best match sample block of sample block to be repaired is searched for, and the new texture information in vertex is added on base surface mesh, spirally updates inner vertex information in region to be repaired successively to be similar to;Finally according to updated grid vertex information, addition constraints solves system of linear equations, and the new coordinate in reverse vertex completes target gridding reconstruct.So as to solve the problems, such as that textural characteristics lack after the textured triangle grid model repairing containing hole.

Description

A kind of model filling-up hole method based on textures synthesis
Technical field
The invention belongs to the methods that 3-dimensional digital geometric manipulations technical field more particularly to grid model are repaired.
Background technology
The hole repair of grid model belongs to one kind in 3-dimensional digital geometric manipulations technology, three-dimensional scanning measurement equipment, The factors such as measuring method and human error cause model data to lack, and after point cloud pretreatment or trigonometric ratio, are formed not Complete surface defect.The formation of hole influences threedimensional model visual effect and later stage geometric manipulations, such as surface modeling and editor Deng, and the development of model repairing technique also provides new means and tool for cultural relic conservation and research.The repairing of general hole Research main method has radial basis functions, diffusion method, B-spline curves interpolation method, template deformation method etc., and hole repair is in itself It is an ill-conditioning problem.Various symmetrical structures can be carried out by extracting symmetry axis, plane of symmetry etc. in the form of symmetry transformation It repairs, and for universal model, how research tendency mainly efficiently uses its peripheral information and realizes that Porous Characteristic details is extensive It is multiple.
Geometric error modeling is commonly used to represent the fine structure on rough geometric model surface, such as fur, fabric surface, synthesizes skill Art also becomes the important means of the enhancing sense of reality;As image texture synthetic technology reaches its maturity, it is extended to three-dimensional geometry The technology of textures synthesis also plays a very important role in fields such as Design and machining of products, analogue simulation, biological informations.Texture Synthetic technology is broadly divided into:It is generated according to the multiple parameters of geometric error modeling parameter model;Printenv synthesizes a pixel every time Point;Based on the multiple pixels of sample single sintering.Because most practical textures have coarse morphology, image line Reason synthetic technology cannot simply meet the synthesis of three-D grain;Most of research all concentrates on image texture and geometric error modeling closes Into theory and method improve on, such as select different geometric error modeling synthetic methods or sample matches criterion could be adjusted to Innovatory algorithm, but geometric error modeling details synthetic effect is not satisfactory.
In conclusion due to existing grid model hole repair method, there are still precision is high, model original geometry texture Minutia cannot obtain the problems such as complete correction, therefore a kind of new it is necessary to propose for the hole of textured grid model Scheme.
Invention content
The purpose of the present invention is to overcome the deficiency in the prior art, proposes a kind of model method for repairing and mending based on textures synthesis. Of less demanding for the little texture model of common Curvature varying itself, geometric error modeling synthetic technology had both been widened this method and had been applicable in The popularity of model, the adaptively sampled method independent of grid vertex neighborhood point also enhance the operability of model repairing And accuracy;At the same time, the texture characteristics after repairing and archetype realization can be connected and transition naturally, compare There is more preferably minutia retentivity in existing model method for repairing and mending.
The technical solution adopted by the present invention to solve the technical problems is:
A kind of model method for repairing and mending based on textures synthesis, flow include following steps:
Step 1) pre-processes;
Textured model M containing hole is inputted, model is carried out at filling-up hole using the conventional method for not considering texture information Reason obtains model Mo, then to MoFairing is carried out to operate to obtain basal plane grid model MS, refer again to multiresolution thought and isolate table Levy the motion vector of minutia information;
Step 2) basal plane Grid Sampling;
To basal plane grid model MSParameterized treatment is carried out, two dimensional surface P is obtained, each vertex sample block is carried out in P Sampling, specifies the range size of single sampling block, and calculates the characteristic information that single sampled point contains;
Step 3) Exemplar Matching;
For the sample block B of sample range covering area to be repaired, its part effective information B is utilized0As between sample block Measuring similarity standard;To in the similar target search region of missing texture, being carried out to each vertex sample block one by one Match, obtain best match sample block A;
Step 4) feature migrates;
It is arrived because the direction of normal on model different location feature vertex is different, therefore by best match sample block A information transfers Unknown message B in sample block B1In the process, it needs to be based on grid model MsThe local frame of foundation is migrated for medium, is protected Demonstrate,prove original attribute of feature to be repaired;
Step 5) spirally updates Porous Characteristic;
After sample block B completes matching and feature migration, basal plane grid model MsForamen primum hole bounds are gradually reduced, and are formed New sample block B ' to be repaired;Then repeat step 3 and step 4 finds fresh target sample block A ' and completes feature migration; Continue to repeat step 3 and step 4 so that area to be repaired information also be similar to spirally to hole border inner more Newly, until its information, all update completion iteration obtains Mn
Step 6) reconstructs target gridding.
Because of model MoWith MsWith identical topological structure, therefore utilize model MoThe Xin Lapu on former perforated Nei Ge vertex Lars coordinate by solving the new cartesian coordinate in system of linear equations reverse target gridding vertex, is finally completed target gridding reconstruct Mt
In a preferred embodiment:The step 1 specifically includes:
A textured triangle grid model containing hole is inputted, model surface hole is closed for simply connected, hole portion Dividing has geometric error modeling missing;Conventional method filling-up hole of texture information is not considered it, and carries out fairing processing, utilizes more points Resolution technology isolates characteristic displacement vector.
In a preferred embodiment:The conventional method filling-up hole method is:
The textured triangle grid model M of the cap holes containing simply connected is inputted, by identifying hole boundary feature, finds institute The region to be repaired;Delaunay Triangulation is carried out to repairing area boundary curve dimensionality reduction, and to the closed area, with entire The net lattice control length of side is that reference judges whether to need the newly-increased point of insertion, and thin according to the round characteristic of sky and side displacement principle progress grid Change, meet the density that newly-increased mesh-density is not less than periphery grid, and avoid generating long-narrow triangular mesh as far as possible;Solve linear side After journey group carries out Mesh Fusion, the model M of missing textural characteristics is obtainedo
In a preferred embodiment:The method of the fairing processing is:
To isolate the motion vector of characterization minutia information, need, by attribute feature is accumulate in grid vertex, to adopt The basal plane grid model M without any texture feature information is obtained with Taubin fairingS, there is no any textural characteristics grids On, it is approximately zero that each apex Laplce, which sits target value,.
In a preferred embodiment:It is described isolate characteristic displacement vector using multi-level technique method be:
High frequency detail feature from the grid model with textural characteristics is separated, will be contained using multi-level technique Laplce's coordinate of textured feature grid does difference operation with Laplce's coordinate without textural characteristics grid;Because of model Mo And MSI-th of vertex viThe topological structure at place is identical, that is, the displacement of characterization grain details geometric properties information separated to Amount.
In a preferred embodiment:The step 2 includes:
Parameterized treatment is carried out to basal plane grid model, each vertex sample block sampling is carried out in parametrization plane, is adjusted Whole parameter determines the range size of single sampling block, and calculates the characteristic information that single sampled point contains.
In a preferred embodiment:It is described to basal plane grid model carry out parameterized treatment method be:
It is after each vertex parameter the result is that according to shortest path in dijkstra's algorithm in mesh parameterization DEM algorithms Previous neighborhood point carries out what diffusive transport obtained on diameter;Added by the relevant information that parametrization vertex has been completed around opposite vertexes Weight average obtains the geodetic vector of each apex and coordinate substrate, then realizes improved dispersion index mapping parameters successively Change method obtains parametrization two dimensional surface P.
In a preferred embodiment:The sample block method of sampling includes:
1. the average side length L of initial voids model M is obtained respectively0, the UV of each apex is sat on grid in parameter field P Mark;
2. single sample block, using vertex as the center position of sampling area, sampled point is whole with spiral array aligned transfer The square distribution of body, and sampling point distributions not necessarily overlap on grid with grid vertex;
3. set SnFor the dimension of sample block square area sampled point array, SlBetween sampled point either vertically or horizontally The unit distance in direction, it is necessary to meet SnFor odd number, SlNo more than L0
4. single sample block square area culminating point considers not as sampled point, the sample centered on current vertex The textural characteristics that block must cover signal period feature or area to be repaired is lacked, sampled point according to spiral array shape Formula is regular to be stored;
5. by Laplce's coordinate information on gore on piece vertex where sampled point, interpolation calculation goes out the sampled point Laplce's coordinate, curvature and method arrow etc. textural characteristics relevant information;
6. calculating all sampling point informations of single sample block, after repeating 5., deposited as a separate unit Storage;
5. and 6. 7. traversing all vertex on grid, repeat, to match for subsequent characteristics and important evidence is provided.
In a preferred embodiment:The step 3 includes:
The sample block B of texture information to be repaired is covered for sample range, its part effective information B can be utilized0As sample Similitude judgment criteria between this block;During searching for the missing most like target sample block of texture sample block, traversal is each A sample block is matched one by one, the minimum as best match sample block A of deviation.
In a preferred embodiment:The similitude judgment criteria is:
Effective information B0That is the sampling point information of known textural characteristics in sample block B, target sample block i.e. its all adopt The characteristic information of sampling point is known;Similitude judgment criteria, that is, sample block B and effective information sampled point in target sample block It is sequentially matched one by one, one group of data of determining method arrow and curvature information deviation minimum.
In a preferred embodiment:The step 4 includes:
It is arrived because the direction of normal on model different location feature vertex is different, therefore by best match sample block A information transfers Unknown message B in sample block B1In the process, it needs to be based on grid model MsThe local frame of foundation carries out feature for medium and moves It moves, it is ensured that original attribute of feature to be repaired.
In a preferred embodiment:The local frame of the foundation is:
1. according to that need to establish local frame grid vertex and parameterize the relationship between source point method arrow, spin matrix is calculated M;
2. a base vector n as local frame is sweared using the method that need to establish local frame grid vertex;
3. upper one of Selecting All Parameters source point tangent plane cuts arrow as benchmark, rotation transformation is carried out using M, obtains office Another base vector e of portion's frame1
4. the base vector n that need to establish local frame grid vertex and another base vector e1 multiplication crosses are obtained secondary method swears e2
5. base vector n, e1, the e2 on current grid vertex have respectively constituted its local frame;
6. repeat above operation the local frame that can obtain other vertex.
In a preferred embodiment:The feature migration includes:
1. the feature of best match sample block is moved to sample block to be repaired, operation object is divided into the central network of sample block Lattice vertex and sampled point;
2. the feature migration between grid vertex, can be by establishing respective local frame, then component carries out corresponding pass It is assignment operation;
3. the feature migration in sample block between sampled point, can be as needed updated by the corresponding sequence of sampling, Then interpolation obtains the texture feature information on patch Intranet lattice vertex to be repaired;
It 4., need to be into row information again for the sample block with the sample block there are overlapping region after current patch update to be repaired Update migrates convenient for subsequent match and feature;
5. after the completion of characteristic matching, other patch texture information updates to be repaired can be completed by continuing to repeat above operation.
In a preferred embodiment:The step 5 is specially:The initial block of samples on coverage hole boundary completes matching and spy It relocates residents from locations to be used for construction of new buildings or factories after shifting, closing of the frontier curve ranges in foramen primum hole are gradually reduced, and form new sample block to be repaired;Then it is repeated in carrying out Exemplar Matching and feature migration operation so that area to be repaired information also be similar to spirally to hole border inner more Newly, until its information, all update completion iteration obtains grid model Mn
In a preferred embodiment:The step 6 is specially:Because of model MoWith MsTopological structure is identical, therefore utilizes MoForamen primum Laplce's coordinate new after each apex feature update, new by solving vertex on system of linear equations reverse target gridding in the region of hole Cartesian coordinate is finally completed target gridding reconstruct Mt
It is of the invention compared with existing the relevant technologies, have the advantages that:
(1) for of less demanding, the applications expanding of geometric error modeling synthetic technology of the little texture model of Curvature varying itself The popularity of this method applicable models, for the model of different texture feature, no matter whether textural characteristics have periodically, no It is whether uniform by its mesh-density, it can realize considerable effect;
(2) the adaptively sampled method independent of grid vertex neighborhood point, enhance model repairing operability and Accuracy reduces the influence error that grid vertex variation causes sample information, spatial sampling problem is converted to two-dimensional geometry Problem, computational efficiency are high, intuitive and easy to understand;
(3) effect model finally realized not only retains and efficiently utilizes the textural characteristics of initial model, but also after repairing The realization of texture characteristics and archetype be connected naturally and transition, compared to existing model method for repairing and mending with more excellent Minutia retentivity.
The present invention is described in further detail with reference to the accompanying drawings and embodiments, but one kind of the present invention is closed based on texture Into model filling-up hole method be not limited to embodiment.
Description of the drawings
Fig. 1 is the flow diagram of the present invention;
Fig. 2 is textured grid model of the embodiment of the present invention containing hole;
Fig. 3 is that the embodiment of filling-up hole effect is carried out to model using the conventional method for not considering texture information;
Fig. 4 is embodiment basal plane grid model of the present invention;
Fig. 5 is the sample block schematic diagram for parameterizing grid vertex in plane;
Fig. 6 is feature principle of displacement schematic diagram between the block to be repaired with target sample on basal plane grid model;
Fig. 7 is that basal plane grid model figure 5 above and Fig. 6 constantly repeat newer path profile, and curve represents spiral rail Mark, dotted line representative sample block more new direction;
Fig. 8 is that embodiment Fig. 7 of the present invention completes the updated basal plane grid model of feature;
Fig. 9 is to have updated characteristic shown in Fig. 8 to carry out the model after Mesh Fusion with initial texture model.
Specific embodiment
As shown in Figure 1, a kind of model filling-up hole method based on textures synthesis, specifically includes following steps:
A kind of model method for repairing and mending based on textures synthesis, specifically includes following steps:
Step 1) pre-processes
A textured triangle grid model containing hole is inputted, model surface hole is closed for simply connected, hole portion Dividing has geometric error modeling missing;Conventional method filling-up hole of texture information is not considered it, and carries out fairing processing, utilizes more points Resolution technology isolates characteristic displacement vector.
It is in the specific practice of step 1:Textured model M containing hole is inputted, using the routine for not considering texture information Method carries out filling-up hole to model and handles to obtain model Mo, then to MoFairing is carried out to operate to obtain basal plane grid model MS, refer again to Multiresolution thought isolates the motion vector of characterization minutia information.
Wherein, the conventional method filling-up hole method for not considering texture information is:
The textured triangle grid model M of the cap holes containing simply connected is inputted, by identifying hole boundary feature, finds institute The region to be repaired;Delaunay Triangulation is carried out to repairing area boundary curve dimensionality reduction, and to the closed area, with entire The net lattice control length of side is that reference judges whether to need the newly-increased point of insertion, and thin according to the round characteristic of sky and side displacement principle progress grid The minimum angle for changing the tri patch formed is maximum;Meanwhile need to meet the density that newly-increased mesh-density is not less than periphery grid, and It avoids generating long-narrow triangular mesh as far as possible;After solving system of linear equations progress Mesh Fusion, the model of missing textural characteristics is obtained Mo
The method for fairing is:
To isolate the motion vector of characterization minutia information, need, by attribute feature is accumulate in grid vertex, to adopt The basal plane grid model M without any texture feature information is obtained with Taubin fairingS, there is no any textural characteristics grids On, it is approximately zero that each apex Laplce, which sits target value,.
The method of the separation characteristic motion vector is:
High frequency detail feature from the grid model with textural characteristics is separated, will be contained using multi-level technique Laplce's coordinate of textured feature grid does difference operation with Laplce's coordinate without textural characteristics grid;Because of model Mo And MSI-th of vertex viThe topological structure at place is identical, that is, the displacement of characterization grain details geometric properties information separated to Amount.
Step 2) basal plane Grid Sampling
Parameterized treatment is carried out to basal plane grid model, each vertex sample block sampling is carried out in parametrization plane, is adjusted Whole parameter determines the range size of single sampling block, and calculates the characteristic information that single sampled point contains.
The specific practice of the step 2 is:To basal plane grid model MSParameterized treatment is carried out, obtains two dimensional surface P, Each vertex sample block sampling is carried out in P, specifies the range size of single sampling block, and calculate the feature that single sampled point contains Information.
Wherein, it is to the method for basal plane grid model progress parameterized treatment:
It is after each vertex parameter the result is that according to shortest path in dijkstra's algorithm in mesh parameterization DEM algorithms Previous neighborhood point carries out what diffusive transport obtained on diameter;Added by the relevant information that parametrization vertex has been completed around opposite vertexes Weight average obtains the geodetic vector of each apex and coordinate substrate, then realizes improved dispersion index mapping parameters successively Change method obtains parametrization two dimensional surface P.
The specific method of sample block sampling is:
1. the average side length L of initial voids model M is obtained respectively0, the UV of each apex is sat on grid in parameter field P Mark;
2. single sample block, using vertex as the center position of sampling area, sampled point is whole with spiral array aligned transfer The square distribution of body, and sampling point distributions not necessarily overlap on grid with grid vertex;
3. set SnFor the dimension of sample block square area sampled point array, SlBetween sampled point either vertically or horizontally The unit distance in direction, it is necessary to meet SnFor odd number, SlNo more than L0
4. single sample block square area culminating point considers not as sampled point, the sample centered on current vertex The textural characteristics that block must cover signal period feature or area to be repaired is lacked, sampled point according to spiral array shape Formula is regular to be stored;
5. by Laplce's coordinate information on gore on piece vertex where sampled point, interpolation calculation goes out the sampled point Laplce's coordinate, curvature and method arrow etc. textural characteristics relevant information;
6. calculating all sampling point informations of single sample block, after repeating 5., deposited as a separate unit Storage;
5. and 6. 7. traversing all vertex on grid, repeat, to match for subsequent characteristics and important evidence is provided.
Step 3) Exemplar Matching
For the sample block of sample range covering area to be repaired, by the use of its part effective information as the phase between sample block Like degree measurement standard;In the target search region similar to missing texture, each vertex sample block is matched one by one, it can Obtain best match sample block.
The specific practice of the step 3 is:
The sample block B of texture information to be repaired is covered for sample range, its part effective information B can be utilized0As sample Similitude judgment criteria between this block;During searching for the missing most like target sample block of texture sample block, traversal is each The characteristic information of a sample block is matched one by one, the minimum as best match sample block A of deviation.
The similitude judgment criteria is:
Effective information B0That is the sampling point information of known textural characteristics in sample block B, target sample block A i.e. its all adopt The characteristic information of sampling point is known;Similitude judgment criteria, that is, sample block B and effective information sampled point in target sample block It is sequentially matched one by one, determining method arrow and curvature information deviation MatchcnOne group of minimum data.
Matchcn=t*Matchc+(1-t)*Matchn (3)
Wherein,WithThe song of s-th of sampled point Laplce's coordinate in difference representative sample block B and target sample block A Rate component information;WithS-th of sampling in representative sample block B and target sample block A respectively The component of curvature information of point Laplce's coordinate;Parameter t represents user-defined weight factor, unified value 0.5.
Step 4) feature migrates
Believe because the direction of normal on model different location feature vertex is different, therefore by best match sample block known features Breath is moved in sample block during unknown message, is needed using the local frame established based on basal plane grid model vertex as medium It is migrated, it is ensured that original attribute of feature to be repaired.
The specific practice of the step 4 is:
It is arrived because the direction of normal on model different location feature vertex is different, therefore by best match sample block A information transfers Unknown message B in sample block B1In the process, it needs to be based on grid model MsThe local frame of foundation carries out feature for medium and moves It moves, it is ensured that original attribute of feature to be repaired.
The local frame of the foundation is:
1. according to that need to establish local frame grid vertex and parameterize the relationship between source point method arrow, spin matrix is calculated M;
2. a base vector n as local frame is sweared using the method that need to establish local frame grid vertex;
3. upper one of Selecting All Parameters source point tangent plane cuts arrow as benchmark, rotation transformation is carried out using M, obtains office Another base vector e of portion's frame1
4. by the base vector n and another base vector e of local frame grid vertex need to be established1Multiplication cross obtains secondary method arrow e2
5. base vector n, e on current grid vertex1、e2Its local frame is respectively constituted;
6. repeat above operation the local frame that can obtain other vertex.
Carry out feature migration specific method be:
1. the feature of best match sample block is moved to sample block to be repaired, operation object is divided into the central network of sample block Lattice vertex and sampled point;
2. the feature migration between grid vertex, can be by establishing respective local frame, then component carries out corresponding pass It is assignment operation;
3. the feature migration in sample block between sampled point, can be as needed updated by the corresponding sequence of sampling, Then interpolation obtains the texture feature information on patch Intranet lattice vertex to be repaired;
It 4., need to be into row information again for the sample block with the sample block there are overlapping region after current patch update to be repaired Update migrates convenient for subsequent match and feature;
5. after the completion of characteristic matching, other patch texture information updates to be repaired can be completed by continuing to repeat above operation.
Step 5) spirally updates Porous Characteristic
After initial block of samples completes matching and feature migration, basal plane grid model foramen primum hole bounds are gradually reduced, shape The sample block to be repaired of Cheng Xin, then carry out step 3 and 4 find fresh target sample block complete feature migration, and repeat step 3 and 4 so that area to be repaired information is also spirally updated with being similar to hole border inner, and until its information, all update is completed Iteration obtains perforated feature and completes the basal plane grid model repaired.
The specific practice of the step 5 is:
After the initial block of samples on coverage hole boundary completes matching and feature migration, foramen primum hole closing of the frontier curve ranges by It is tapered small, form new sample block to be repaired;Then it is repeated in carrying out Exemplar Matching and feature migration operation so as to be repaired Multiple area information is also spirally updated with being similar to hole border inner, and until its information, all iteration is completed in update, is obtained Perforated feature completes the basal plane grid model M repairedn
Step 6) reconstruct target gridding
Because there is identical topological structure with basal plane grid model using model after the conventional filling-up hole for not considering texture information, Therefore new Laplce's coordinate on each vertex of model after conventional filling-up hole is utilized, go out target gridding vertex by solving system of linear equations reverse New cartesian coordinate is finally completed target gridding reconstruct.
The specific practice of the step 6 is:
Because of model M after conventional filling-up holeoWith basal plane grid model MsWith identical topological structure, therefore after utilizing conventional filling-up hole Model MoLaplce's coordinate new after each apex feature update, goes out target by solving system of linear equations reverse in former perforated The new cartesian coordinate in vertex on grid is finally completed target gridding reconstruct Mt
Embodiment 1
Demonstration citing is hereafter carried out to above-mentioned step by specific example:
Step 1, preprocessing process.
Specifically, the model M containing hole (as shown in Figure 2) that input is textured, using the routine side for not considering texture information Method carries out filling-up hole to model and handles to obtain model Mo(as shown in Figure 3), then to MoFairing is carried out to operate to obtain basal plane grid model MS(as shown in Figure 4) refers again to the motion vector that multiresolution thought isolates characterization minutia information.
Step 2, basal plane Grid Sampling.
Specifically, to basal plane grid model MSParameterized treatment is carried out, by having completed grid DEM algorithms around opposite vertexes Vertex relevant information weighted average after parametrization obtains the geodetic vector of each apex and coordinate substrate successively, then real Existing improved dispersion index mapping parameters method obtains parametrization two dimensional surface P, each vertex sample block is carried out in P and is adopted Sample, specifies the range size of single sampling block, and calculates the characteristic information that single sampled point contains.
When sample block samples, the average side length L of initial voids model M is obtained respectively first0, it is each on grid in parameter field P The UV coordinates of a apex, secondly using vertex as the center position in sample block region, sampled point with spiral array rule sequentially Arrangement, whole square distribution, finally, using Laplce's coordinate information on gore on piece vertex where sampled point, Interpolation calculation goes out the textural characteristics relevant informations such as Laplce's coordinate, curvature and the method arrow of each sampled point, repeats above-mentioned behaviour Sampling process can be completed.
Step 3, Exemplar Matching.
Specifically, covering the sample block B (as shown in Figure 5) of texture information to be repaired for sample range, its portion can be utilized Divide effective information B0As similitude judgment criteria between sample block;In search and the missing most like target sample of texture sample block During block, the characteristic information for traversing each sample block is matched one by one, the minimum as best match sample block of deviation A。
Step 4, feature migrates.
It is arrived because the direction of normal on model different location feature vertex is different, therefore by best match sample block A information transfers Unknown message B in sample block B1In the process, it needs to be based on grid model MsThe local frame of foundation carries out feature migration for medium (as shown in Figure 6), it is ensured that original attribute of feature to be repaired.
Step 5, Porous Characteristic is spirally updated.
After the initial block of samples on coverage hole boundary completes matching and feature migration, foramen primum hole closing of the frontier curve ranges by It is tapered small, form new sample block to be repaired;Then it is repeated in carrying out Exemplar Matching and feature migration operation so as to be repaired Multiple area information also spirally updates (as shown in Figure 7) to be similar to hole border inner, until its information has all updated Into iteration, obtain perforated feature and complete the basal plane grid model M repairedn(as shown in Figure 8).
Step 6, target gridding is reconstructed.
Because of model M after conventional filling-up holeoWith basal plane grid model MsWith identical topological structure, therefore after utilizing conventional filling-up hole Model MoLaplce's coordinate new after each apex feature update, goes out target by solving system of linear equations reverse in former perforated The new cartesian coordinate in vertex on grid is finally completed target gridding reconstruct Mt(as shown in Figure 9).
Fig. 2, Fig. 3, Fig. 4, Fig. 8, Fig. 9 are embodiment of the present invention:Fig. 2 is textured grid model of the input containing hole, The STL models include 27321 vertex and 53615 tri patch altogether;Fig. 3 is using the conventional method for not considering texture information Model after filling-up hole, Fig. 4 are embodiment basal plane grid models of the present invention, and Fig. 8 is that embodiment Fig. 7 of the present invention completes feature update Basal plane grid model afterwards, the STL models include 33109 vertex and 65457 tri patch;Implement the model of the present invention After filling-up hole method, obtained texture model is as shown in figure 9, comprising 33109 vertex and 65457 tri patch, whole process It is automatically performed.
The foregoing is merely present pre-ferred embodiments, therefore the technical scope of the present invention cannot be limited according to this, therefore Fan Yiben The equivalent changes and modifications that the technical spirit and description of invention are made, in the range of should all belonging to technical solution of the present invention.

Claims (15)

1. a kind of model method for repairing and mending based on textures synthesis, it is characterised in that:Its flow includes following steps:
Step 1) pre-processes;
Textured model M containing hole is inputted, carrying out filling-up hole to model using the conventional method for not considering texture information is handled To model Mo, then to MoFairing is carried out to operate to obtain basal plane grid model MS, refer again to multiresolution thought and isolate characterization carefully Save the motion vector of characteristic information;
Step 2) basal plane Grid Sampling;
To basal plane grid model MSParameterized treatment is carried out, obtains two dimensional surface P, each vertex sample block sampling is carried out in P, The range size of specified single sample block, and calculate the characteristic information that single sampled point contains;
Step 3) Exemplar Matching;
For the sample block B of sample range covering area to be repaired, its part effective information B is utilized0As similar between sample block Spend measurement standard;In the target search region similar to missing texture, each vertex subsamples block is matched one by one, Obtain best match sample block A;
Step 4) feature migrates;
Because the direction of normal on different location feature vertex on model is different, therefore by best match sample block A information transfers to sample Unknown message B in this block B1In the process, it needs to be based on grid model MsThe local frame of foundation is migrated for medium, is ensured Original attribute of feature to be repaired;
Step 5) spirally updates Porous Characteristic;
After sample block B completes matching and feature migration, basal plane grid model MsForamen primum hole bounds are gradually reduced, and are formed new Sample block B ' to be repaired;Then repeat step 3 and step 4 finds fresh target sample block A ' and completes feature migration;Continue Repeat step 3 and step 4 so that area to be repaired information is also spirally updated to be similar to hole border inner, directly To its information, all update completion iteration obtains Mn
Step 6) reconstructs target gridding.
Because of model MoWith MsWith identical topological structure, therefore utilize model MoThe new Laplce on former perforated Nei Ge vertex Coordinate by solving the new cartesian coordinate in system of linear equations reverse target gridding vertex, is finally completed target gridding reconstruct Mt
2. the model method for repairing and mending according to claim 1 based on textures synthesis, which is characterized in that the step 1 is specific Including:
A textured triangle grid model containing hole is inputted, model surface hole is simply connected closing, and hole part has Geometric error modeling lacks;Conventional method filling-up hole of texture information is not considered it, and carries out fairing processing, utilizes multiresolution Technology isolates characteristic displacement vector.
3. the model method for repairing and mending according to claim 2 based on textures synthesis, which is characterized in that the conventional method is mended Hole method is:
The textured triangle grid model M of the cap holes containing simply connected is inputted, by identifying hole boundary feature, searching to be repaiied The region of benefit;Delaunay Triangulation is carried out to repairing area boundary curve dimensionality reduction, and to the closed area, with entire grid Average side length is that reference judges whether to need to be inserted into newly-increased point, and carry out mesh refinement according to the round characteristic of sky and side displacement principle, Meet condition of the newly-increased mesh-density not less than the density of periphery grid, and avoid generating long-narrow triangular mesh as far as possible;Solve line Property equation group carry out Mesh Fusion after, obtain missing textural characteristics model Mo
4. the model method for repairing and mending according to claim 2 based on textures synthesis, which is characterized in that the fairing processing Method is:
To isolate the motion vector of characterization minutia information, need, by attribute feature is accumulate in grid vertex, to use Taubin fairing obtains the basal plane grid model M without any texture feature informationS, and there is no on any textural characteristics grid It is approximately zero that each apex Laplce, which sits target value,.
5. the model method for repairing and mending according to claim 2 based on textures synthesis, which is characterized in that described to be differentiated using more The method that rate technology isolates characteristic displacement vector is:
High frequency detail feature from the grid model with textural characteristics is separated, will contain line using multi-level technique The Laplce's coordinate and Laplce's coordinate without textural characteristics grid for managing feature grid do difference operation;Because of model MoAnd MS I-th of vertex viThe topological structure at place is identical, that is, the motion vector of characterization grain details geometric properties information separated.
6. the model method for repairing and mending according to claim 1 based on textures synthesis, which is characterized in that the step 2 includes:
Parameterized treatment is carried out to basal plane grid model, each vertex sample block sampling, adjustment ginseng are carried out in parametrization plane The range size of the determining single sample block of number, and calculate the characteristic information that single sampled point contains.
7. the model method for repairing and mending according to claim 6 based on textures synthesis, which is characterized in that described to basal plane grid Model carry out parameterized treatment method be:
It is after each vertex parameter the result is that according on shortest path in dijkstra's algorithm in mesh parameterization DEM algorithms Previous neighborhood point carries out what diffusive transport obtained;It is flat by the relevant information weighting that parametrization vertex has been completed around opposite vertexes , the geodetic vector of each apex and coordinate substrate are obtained successively, then realize improved dispersion index mapping parameters side Method obtains parametrization two dimensional surface P.
8. the model method for repairing and mending according to claim 6 based on textures synthesis, which is characterized in that the sample block sampling Method includes:
1. the average side length L of initial voids model M is obtained respectively0, in parameter field P on grid each apex UV coordinates;
2. single sample block, using vertex as the center position of sampling area, sampled point is in integrally with spiral array aligned transfer Square profile, and sampling point distributions not necessarily overlap on grid with grid vertex;
3. set SnFor the dimension of sample block square area sampled point array, SlThe either vertically or horizontally direction between sampled point Unit distance, it is necessary to meet SnFor odd number, SlNo more than L0
4. single sample block square area culminating point considers that the sample block centered on current vertex must not as sampled point The textural characteristics that palpus covering signal period feature or area to be repaired are lacked, sampled point have according to the form of spiral array Rule is stored;
5. by Laplce's coordinate information on gore on piece vertex where sampled point, interpolation calculation goes out the drawing of the sampled point The textural characteristics relevant informations such as this coordinate of pula, curvature and method arrow;
6. calculating all sampling point informations of single sample block, after repeating 5., stored as a separate unit,
5. and 6. 7. traversing all vertex on grid, repeat, to match for subsequent characteristics and important evidence is provided.
9. the model method for repairing and mending according to claim 1 based on textures synthesis, which is characterized in that the step 3 includes:
The sample block B of texture information to be repaired is covered for sample range, its part effective information B can be utilized0As sample block it Between similitude judgment criteria;During searching for the missing most like target sample block of texture sample block, each sample is traversed Block is matched one by one, the minimum as best match sample block A of deviation.
10. the model method for repairing and mending according to claim 9 based on textures synthesis, which is characterized in that the similitude is commented Sentencing standard is:
Effective information B0That is the sampling point information of known textural characteristics in sample block B, its i.e. all sampled point of target sample block Characteristic information is known;In similitude judgment criteria, that is, sample block B and target sample block effective information sampled point sequentially into Row matches one by one, one group of data of determining method arrow and curvature information deviation minimum.
11. the model method for repairing and mending according to claim 1 based on textures synthesis, which is characterized in that the step 4 is wrapped It includes:
Because the direction of normal on different location feature vertex on model is different, therefore by best match sample block A information transfers to sample Unknown message B in this block B1In the process, it needs to be based on grid model MsThe local frame of foundation carries out feature migration for medium, It can ensure original attribute of feature to be repaired.
12. the model method for repairing and mending according to claim 11 based on textures synthesis, which is characterized in that the office of the foundation Portion's frame is:
1. according to that need to establish local frame grid vertex and parameterize the relationship between source point method arrow, spin matrix M is calculated;
2. a base vector n as local frame is sweared using the method that need to establish local frame grid vertex;
3. upper one of Selecting All Parameters source point tangent plane cuts arrow as benchmark, rotation transformation is carried out using M, obtains local mark Another base vector e of frame1
4. by the base vector n and another base vector e of local frame grid vertex need to be established1Multiplication cross obtains secondary method arrow e2
5. base vector n, e on current grid vertex1、e2Its local frame is respectively constituted;
6. repeat above operation the local frame that can obtain other vertex.
13. the model method for repairing and mending according to claim 11 based on textures synthesis, which is characterized in that the feature migration Including:
1. the feature of best match sample block is moved to sample block to be repaired, operation object is divided into the central gridding top of sample block Point and sampled point;
2. the feature migration between grid vertex, can be by establishing respective local frame, then component carries out correspondence tax Value Operations;
3. the feature migration in sample block between sampled point, can be as needed updated, then by the corresponding sequence of sampling Interpolation obtains the texture feature information on patch Intranet lattice vertex to be repaired;
4. after current patch update to be repaired, for the sample block with the sample block there are overlapping region, need to again be updated into row information, It is migrated convenient for subsequent match and feature;
5. after the completion of characteristic matching, other patch texture information updates to be repaired can be completed by continuing to repeat above operation.
14. the model method for repairing and mending according to claim 1 based on textures synthesis, which is characterized in that the step 5 is specific For:After the initial block of samples on coverage hole boundary completes matching and feature migration, closing of the frontier curve ranges in foramen primum hole gradually contract It is small, form new sample block to be repaired;Then it is repeated in carrying out Exemplar Matching and feature migration operation so that area to be repaired Domain information is also spirally updated with being similar to hole border inner, and until its information, all update completion iteration obtains grid mould Type Mn
15. the model method for repairing and mending according to claim 1 based on textures synthesis, which is characterized in that the step 6 is specific For:Because of model MoWith MsTopological structure is identical, therefore utilizes MoLaplce new after each apex feature update sits in former perforated Mark by solving the new cartesian coordinate in vertex on system of linear equations reverse target gridding, is finally completed target gridding reconstruct Mt
CN201711364637.9A 2017-12-18 2017-12-18 Model hole filling method based on texture synthesis Active CN108230452B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711364637.9A CN108230452B (en) 2017-12-18 2017-12-18 Model hole filling method based on texture synthesis

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711364637.9A CN108230452B (en) 2017-12-18 2017-12-18 Model hole filling method based on texture synthesis

Publications (2)

Publication Number Publication Date
CN108230452A true CN108230452A (en) 2018-06-29
CN108230452B CN108230452B (en) 2021-06-29

Family

ID=62649732

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711364637.9A Active CN108230452B (en) 2017-12-18 2017-12-18 Model hole filling method based on texture synthesis

Country Status (1)

Country Link
CN (1) CN108230452B (en)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109086492A (en) * 2018-07-11 2018-12-25 大连理工大学 A kind of wire frame representation of body structure threedimensional model and deformation method and system
CN109118589A (en) * 2018-09-07 2019-01-01 北京航空航天大学青岛研究院 Geometrical model full range details restorative procedure
CN109615702A (en) * 2018-10-23 2019-04-12 深圳市华讯方舟太赫兹科技有限公司 Grid filling-up hole method, image processing equipment and the device with store function
CN109872394A (en) * 2019-01-10 2019-06-11 重庆大学 Long-narrow triangular mesh grid optimization method based on least square method supporting vector machine
CN112381950A (en) * 2020-11-05 2021-02-19 新拓三维技术(深圳)有限公司 Grid hole repairing method, electronic equipment and computer readable storage medium

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006348574A (en) * 2005-06-15 2006-12-28 Matsushita Electric Works Ltd Moisture absorbing-releasing material with decoration
CN101853525A (en) * 2010-05-19 2010-10-06 北京航空航天大学 Mesh segmentation based simplification method for preserving details of textured model
CN104361632A (en) * 2014-11-03 2015-02-18 北京航空航天大学 Triangular mesh hole-filling method based on Hermite radial basis function
CN105069833A (en) * 2015-07-29 2015-11-18 苏州华漫信息服务有限公司 Three-dimensional model repair method capable of retaining texture information data
CN105427252A (en) * 2015-11-06 2016-03-23 北京航空航天大学 Mesh model multi-scale geometric detail repair method based on empirical mode decomposition

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006348574A (en) * 2005-06-15 2006-12-28 Matsushita Electric Works Ltd Moisture absorbing-releasing material with decoration
CN101853525A (en) * 2010-05-19 2010-10-06 北京航空航天大学 Mesh segmentation based simplification method for preserving details of textured model
CN104361632A (en) * 2014-11-03 2015-02-18 北京航空航天大学 Triangular mesh hole-filling method based on Hermite radial basis function
CN105069833A (en) * 2015-07-29 2015-11-18 苏州华漫信息服务有限公司 Three-dimensional model repair method capable of retaining texture information data
CN105427252A (en) * 2015-11-06 2016-03-23 北京航空航天大学 Mesh model multi-scale geometric detail repair method based on empirical mode decomposition

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
LIONEL C.GONTARD 等: "Photogrammetry of the three-dimensional shape and texture of a nanoscale particle using scanning electron microscopy and free software", 《ULTRAMICROSCOPY》 *
马爽: "纹理图像智能修补关键技术研究", 《万方数据库》 *

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109086492A (en) * 2018-07-11 2018-12-25 大连理工大学 A kind of wire frame representation of body structure threedimensional model and deformation method and system
CN109086492B (en) * 2018-07-11 2022-12-13 大连理工大学 Wire frame representation and deformation method and system for three-dimensional model of vehicle body structure
CN109118589A (en) * 2018-09-07 2019-01-01 北京航空航天大学青岛研究院 Geometrical model full range details restorative procedure
CN109615702A (en) * 2018-10-23 2019-04-12 深圳市华讯方舟太赫兹科技有限公司 Grid filling-up hole method, image processing equipment and the device with store function
CN109615702B (en) * 2018-10-23 2023-09-22 深圳市重投华讯太赫兹科技有限公司 Grid hole filling method, image processing equipment and device with storage function
CN109872394A (en) * 2019-01-10 2019-06-11 重庆大学 Long-narrow triangular mesh grid optimization method based on least square method supporting vector machine
CN112381950A (en) * 2020-11-05 2021-02-19 新拓三维技术(深圳)有限公司 Grid hole repairing method, electronic equipment and computer readable storage medium

Also Published As

Publication number Publication date
CN108230452B (en) 2021-06-29

Similar Documents

Publication Publication Date Title
CN108230452A (en) A kind of model filling-up hole method based on textures synthesis
CN110516388B (en) Harmonic mapping-based curved surface discrete point cloud model circular cutter path generation method
CN109377561A (en) A kind of digital-to-analogue surface mesh generation method based on conformal geometry
CN108171780A (en) A kind of method that indoor true three-dimension map is built based on laser radar
CN109360260A (en) A kind of cut-away restructing algorithm of triangle gridding threedimensional model
CN108595858A (en) Lightweight advanced treating working method based on BIM
CN103903061B (en) Information comprehensive processing device and method in three-dimensional mineral resource prediction evaluation
CN109872394A (en) Long-narrow triangular mesh grid optimization method based on least square method supporting vector machine
CN109472870B (en) Model matching method based on grid reconstruction and multi-influence-domain correction
CN104835202A (en) Quick three-dimensional virtual scene constructing method
CN104572924B (en) Multi-scale expression information generating method for GIS vector building polygons
CN109584357A (en) Three-dimensional modeling method, device, system and storage medium based on more contour lines
CN107610228A (en) Curved surface increment topology rebuilding method based on massive point cloud
Zheng et al. A morphologically preserved multi-resolution TIN surface modeling and visualization method for virtual globes
Pfeifer A subdivision algorithm for smooth 3D terrain models
CN106875484A (en) A kind of geology accumulation body Fast Fitting modeling method based on dimensional topography
CN104537254B (en) A kind of drafting method that becomes more meticulous based on social statistics data
CN105931297A (en) Data processing method applied to three-dimensional geological surface model
CN105869210A (en) Interpolation data processing method in three-dimensional geological surface model
CN110738726B (en) Robot vision-guided three-dimensional object reconstruction method based on octree
CN112380306A (en) Adaptive correction method for Kergin spatial interpolation considering distribution balance
CN104463924B (en) Digital elevation terrain model generation method based on scattered point elevation sample data
Auerbach et al. Computer-aided geometric design of geologic surfaces and bodies
CN110163973A (en) A kind of zone boundary ring method for fairing
CN114817854B (en) Rapid multi-point simulation method oriented to continuous value variable and based on linear regression

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