CN110047151A - A kind of fragment plane of disruption joining method and system based on thermonuclear feature - Google Patents
A kind of fragment plane of disruption joining method and system based on thermonuclear feature Download PDFInfo
- Publication number
- CN110047151A CN110047151A CN201910239060.1A CN201910239060A CN110047151A CN 110047151 A CN110047151 A CN 110047151A CN 201910239060 A CN201910239060 A CN 201910239060A CN 110047151 A CN110047151 A CN 110047151A
- Authority
- CN
- China
- Prior art keywords
- disruption
- matching
- fragment
- characteristic point
- point
- 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
Links
- 239000012634 fragment Substances 0.000 title claims abstract description 115
- 238000000034 method Methods 0.000 title claims abstract description 26
- 239000011159 matrix material Substances 0.000 claims description 42
- 238000013519 translation Methods 0.000 claims description 39
- 238000004364 calculation method Methods 0.000 claims description 22
- 230000009466 transformation Effects 0.000 claims description 18
- 238000010276 construction Methods 0.000 claims description 8
- 238000012216 screening Methods 0.000 claims description 7
- 238000001914 filtration Methods 0.000 claims description 3
- 238000012545 processing Methods 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 5
- 239000000284 extract Substances 0.000 description 3
- 229910000906 Bronze Inorganic materials 0.000 description 2
- 239000010974 bronze Substances 0.000 description 2
- KUNSUQLRTQLHQQ-UHFFFAOYSA-N copper tin Chemical compound [Cu].[Sn] KUNSUQLRTQLHQQ-UHFFFAOYSA-N 0.000 description 2
- 230000001154 acute effect Effects 0.000 description 1
- 230000000295 complement effect Effects 0.000 description 1
- 238000005260 corrosion Methods 0.000 description 1
- 230000007797 corrosion Effects 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000017105 transposition Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T3/00—Geometric image transformations in the plane of the image
- G06T3/14—Transformations for image registration, e.g. adjusting or mapping for alignment of images
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/20—Finite element generation, e.g. wire-frame surface description, tesselation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T19/00—Manipulating 3D models or images for computer graphics
- G06T19/20—Editing of 3D images, e.g. changing shapes or colours, aligning objects or positioning parts
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
- Y02P90/30—Computing systems specially adapted for manufacturing
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Graphics (AREA)
- Software Systems (AREA)
- Geometry (AREA)
- Architecture (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
- Investigating Or Analyzing Materials Using Thermal Means (AREA)
Abstract
The invention discloses a kind of fragment plane of disruption joining methods and system based on thermonuclear feature, belong to Computer Image Processing field, include: the local curvature's maximum point for calculating separately two fragment planes of disruption to be matched, these are put and is used as respective set of characteristic points；Calculate the thermonuclear feature of characteristic point in the set of characteristic points；It is matched based on the characteristic point of the thermonuclear feature of the characteristic point to described two fragment planes of disruption, building matching is to set；Based on the matching to set, described two fragment planes of disruption are spliced.The present invention improves the accuracy of fragments mosaicing.
Description
Technical field
The present invention relates to Computer Image Processing field, in particular to a kind of fragment plane of disruption splicing based on thermonuclear feature
Method and system.
Background technique
Effective feature description is the key that realize fragments mosaicing, in existing research, with curvature, profile, the plane of disruption
It is widely used in the splicing of fragment model Deng the feature for representative.Single contour feature more application is in the spelling of thin-walled fragment
It connects, such as Karasik et al. extracts pottery symmetry axis and profile information from a large amount of pottery sherds and analyzes for pottery；Lin Suzhen
Et al. extract the profile of thin-walled bronze ware two dimensional image fragment, and further calculate its curvature, utilize Local Extremum abutment points
Between chord length construction feature descriptor realize bronze ware effective splicing；Fan Shaorong et al. proposes to extract fragment Internal periphery respectively
Curve and outer profile curve are spliced using complementary realize between fragment curve shape；Week, bright congruent people was mentioned based on fragment contour line
A kind of interactive joining method out, effectively realizes historical relic Virtual restora- tion.However, using contour curve as the joining method of main feature
The disadvantage is that there are difficulty and inaccuracy for splicing when the reasons such as missing or corrosion prevent fragment contour curve from completely presenting
Property.
Many fragments often have certain thickness, at bulk, compared with thin-walled fragments mosaicing, the fragments mosaicing with thickness
It is more challenging, it possibly can not be matched by single contour feature, extract and be properly characterized in the accurate plane of disruption
Solve the problems, such as such key.Such as the complete plane of disruption is divided into several small curved surfaces by vertex curvature by Lie group brightness et al.,
Character pair point is obtained between the similarity these curved surfaces, and then fragment is matched；It is constant that MANAY et al. proposes integral
Measure the matching that concept is used for two-dimensional shapes；This concept is expanded to threedimensional model by Pottmann et al., utilizes integral invariant
The concave-convex degree of curved surface is obtained, effectively judgement characteristic point, and the feature extracted has robustness to the noise in data, it can be more
Realize that the plane of disruption matches well.
The curvature or integral invariant on vertex, these information belong to typical local feature either on the plane of disruption,
The splicing of fragment depends on these local features, and they calculate and are easy.However this feature lacks depth, separating capacity
It is weak.Another category feature considers the distribution of shapes of neighborhood of a point, typically Spin image, and this characteristic information is richer, therefore
Separating capacity is stronger, and its advantages are to remain unchanged under rigid body translation, but its calculating is complicated and time-consuming.
Summary of the invention
The purpose of the present invention is to provide a kind of fragment plane of disruption joining methods and system based on thermonuclear feature, to improve
The accuracy of fragments mosaicing.
In order to achieve the above object, on the one hand, the present invention uses a kind of fragment plane of disruption joining method based on thermonuclear feature,
Include the following steps:
These are put and are used as respective spy by the local curvature's maximum point for calculating separately two fragment planes of disruption to be matched
Levy point set；
Calculate the thermonuclear feature of characteristic point in the set of characteristic points；
It is matched, is constructed based on the characteristic point of the thermonuclear feature of the characteristic point to described two fragment planes of disruption
Matching is to set；
Based on the matching to set, described two fragment planes of disruption are spliced.
Preferably, the local curvature's maximum point for calculating separately two fragment planes of disruption to be matched, by these points
Include: as respective set of characteristic points
The local curvature's maximum point for calculating separately two fragment planes of disruption to be matched, by these points as respective first
Beginning set of characteristic points；
The characteristic point is subjected to triangle gridding；
Point in feature vertex neighborhood after the triangle gridding is added in the initial feature point set and is expanded
Exhibition, obtains set of characteristic points.
Preferably, the thermonuclear feature for calculating characteristic point in the set of characteristic points includes:
Calculate the Laplace-Beltrami operator of characteristic point in the set of characteristic points；
The thermonuclear of characteristic point in the set of characteristic points is calculated based on the Laplace-Beltrami operator of the characteristic point
Feature.
Preferably, the Laplace-Beltrami operator of the characteristic point is calculated using grid Laplace operator.
Preferably, described to be clicked through based on the feature of the thermonuclear feature of the characteristic point to described two fragment planes of disruption
Row matching, building matching include: to set
Thermonuclear feature based on the characteristic point calculates each characteristic point and another fragment in a fragment plane of disruption
The distance of all characteristic points in the plane of disruption；
It will be apart from the smallest two characteristic points as a matching pair；
The matching pair of each characteristic point in one fragment plane of disruption is obtained, building matching is to set.
Preferably, described to be matched based on described to set, carrying out splicing to described two fragment planes of disruption includes:
Based on the matching to the matching pair in set, all rigid body translation matrixes are calculated；
Each matching is calculated to the global error under all rigid body translation matrixes；
Using the rigid body translation matrix under the global error minimum as general transformation matrix；
Under the general transformation matrix, described two fragment planes of disruption are spliced.
Preferably, described to be matched based on described to the matching pair in set, calculating all rigid body translation matrixes includes:
Based on the matching to the matching pair in set, three all not conllinear matchings pair are obtained；
Calculate the rigid body translation matrix of three all not conllinear matchings pair.
Second aspect, the present invention use a kind of fragment plane of disruption splicing system based on thermonuclear feature, including characteristic point to obtain
Modulus block, thermonuclear feature calculation module, matching module and splicing module；
Characteristic point obtains local curvature's maximum point that module is used to obtain two fragment planes of disruption to be matched, by these
Point is used as respective set of characteristic points；
Thermonuclear feature calculation module is used to calculate the thermonuclear feature of characteristic point in the set of characteristic points；
Matching module is for the thermonuclear feature based on the characteristic point to the characteristic point of described two fragment planes of disruption
It is matched, building matching is to set；
Splicing module is used to splice described two fragment planes of disruption set based on the matching.
Preferably, the matching module includes metrics calculation unit, judging unit and construction unit；
Metrics calculation unit is used for the thermonuclear feature based on the characteristic point, and it is special to calculate each in a fragment plane of disruption
Sign point is at a distance from all characteristic points in another fragment plane of disruption；
Judging unit will be made for judging whether the distance between two characteristic points is minimum apart from the smallest two characteristic points
For a matching pair；
Construction unit is used to obtain the matching pair of each characteristic point in one fragment plane of disruption, constructs a pairing set
It closes.
Preferably, the splicing module includes converter unit, error calculation unit, screening unit and concatenation unit；
Converter unit is used to calculate all rigid body translation matrixes to the matching pair in set based on the matching；
Error calculation unit is for calculating each matching to the global error under all rigid body translation matrixes；
Screening unit is for filtering out minimum global error, and by the rigid body translation square under the global error minimum
Battle array is used as general transformation matrix；
Concatenation unit is used under the general transformation matrix, is spliced to described two fragment planes of disruption.
Compared with existing algorithm, the invention has the following advantages:
By calculating local curvature's maximum point of the fragment plane of disruption, it regard these points as characteristic point, it is broken avoids calculating
All points in the piece plane of disruption reduce calculation amount to reduce matching search space；
The fragment plane of disruption is matched using the thermonuclear feature of characteristic point, has fully taken into account the characteristic point institute of the plane of disruption
Different neighborhoods, so that matching result is had more robustness；
Since rigid body translation corresponds to rotation and translation, the smaller matching shown a little of error is more accurate, therefore calculates each
The global error under all rigid body translation matrixes is matched, using the rigid body translation matrix under global error minimum as total
Body transformation matrix splices fragment under the constraint of general transformation matrix, and splicing is made to have more accuracy.
Detailed description of the invention
With reference to the accompanying drawing, specific embodiments of the present invention will be described in detail:
Fig. 1 is the flow diagram of fragments mosaicing method；
Fig. 2 is the calculation method schematic diagram of the area of adjacent triangle；
Fig. 3 is the workflow schematic diagram of fragments mosaicing system；
Fig. 4 is the workflow schematic diagram of matching module in fragments mosaicing system；
Fig. 5 is the workflow schematic diagram of splicing module in fragments mosaicing system.
Specific embodiment
In order to further explain feature of the invention, reference should be made to the following detailed description and accompanying drawings of the present invention.Institute
Attached drawing is only for reference and purposes of discussion, is not used to limit protection scope of the present invention.
Embodiment one
As shown in Figure 1, the present invention uses a kind of fragment plane of disruption joining method based on thermonuclear feature, including walk as follows
It is rapid:
S1, the local curvature's maximum point for calculating separately two fragment planes of disruption to be matched, by these points as respective
Set of characteristic points；
S2, the thermonuclear feature for calculating characteristic point in the set of characteristic points；
S3, it is matched based on the characteristic point of the thermonuclear feature of the characteristic point to described two fragment planes of disruption,
Building matching is to set；
S4, it is matched based on described to set, described two fragment planes of disruption is spliced.
Specifically, the thermonuclear feature for calculating all the points on the fragment plane of disruption is avoided in order to reduce matching search space, point
Two fragment plane of disruption F to be matched are not calculatediAnd FjThese are put and are used as respective characteristic point by upper local curvature's maximum point
Set.
Here with one of fragment plane of disruption FiFor, provide the method for calculating local curvature's maximum point:
Calculate point piGaussian curvature:
Wherein, θjIndicate pjpiWith pj+1piAngle, piFor fragment plane of disruption FiIn point, pjFor fragment plane of disruption FjIn
Point, N (i) indicate piAbutment points, A indicates the hybrid projection area of the adjacent triangle in each vertex, as shown in Figure 2:
If adjacent triangle is acute triangle, point p is takeniThe midpoint on connected both sides and triangle excenter, by this three
Point and piComposed area is A1；If adjacent triangle is obtuse angle or right angled triangle, point p is takeniConnected both sides
Midpoint and piThe midpoint of institute's opposite side, by this 3 points and piComposed area is A2。
Calculate point piAverage curvature:
Wherein, αijIndicate pjpj-1With pipj-1Angle, βijIndicate pjpj+1With pipj+1Angle, NiIndicate Δ pipjpj+1
Normal vector.
Calculate point piMaximum principal curvatures K1(pi) and minimum principal curvatures K2(pi):
By K1(pi) and K2(pi) local curvature's maximum can be calculated.
It should be noted that having had the method for calculating maximum based on maximum value and minimum value in the prior art, here
It repeats no more.
After calculating local curvature's maximum, then calculate the thermonuclear feature of these local curvature's maximum.
Based on the thermonuclear feature of local curvature's maximum to fragment plane of disruption FiAnd FjIt is matched, building matching is to set
P.Use fiuIndicate fragment plane of disruption FiIn characteristic point, use fjvIndicate fragment plane of disruption FjIn characteristic point, then match actually
It exactly finds and each characteristic point fiuMatched characteristic point fjv。
Based on matching to set P, to fragment plane of disruption FiAnd FjSpliced.
Preferably, the local curvature's maximum point for calculating separately two fragment planes of disruption to be matched, by these conducts
Respective set of characteristic points include:
The local curvature's maximum point for calculating separately two fragment planes of disruption to be matched, by these points as respective first
Beginning set of characteristic points；
The characteristic point is subjected to triangle gridding；
Point in feature vertex neighborhood after the triangle gridding is added in the initial feature point set and is expanded
Exhibition, obtains set of characteristic points.
Specifically, reducing stitching error in order to improve splicing accuracy, obtaining initial feature point set firstWith, these characteristic points are then subjected to triangle gridding, the point in the feature vertex neighborhood after these triangle griddings is added to
Initial feature point setWithIn be extended, obtain set of characteristic pointsWith。
It should be noted that had the method that triangle gridding is carried out to characteristic point in the prior art, it is no longer superfluous here
It states.
Preferably, the thermonuclear feature for calculating characteristic point in the set of characteristic points includes:
Calculate the Laplace-Beltrami operator of characteristic point in the set of characteristic points；
The thermonuclear of characteristic point in the set of characteristic points is calculated based on the Laplace-Beltrami operator of the characteristic point
Feature.
Preferably, the Laplace-Beltrami operator of characteristic point is calculated using grid Laplace operator.
Specifically, calculating the formula of the thermonuclear feature of characteristic point are as follows:
Wherein, x indicates characteristic point, λsIndicate s-th of characteristic value of Laplace-Beltrami operator, ζsIt indicates
S-th of feature vector of Laplace-Beltrami operator,Indicate point x in the thermonuclear value of moment t.
It is preferably based on characteristic point progress of the thermonuclear feature to described two fragment planes of disruption of the characteristic point
Match, building matching includes: to set
Thermonuclear feature based on the characteristic point calculates each characteristic point and another fragment in a fragment plane of disruption
The distance of all characteristic points in the plane of disruption；
It will be apart from the smallest two characteristic points as a matching pair；
The matching pair of each characteristic point in one fragment plane of disruption is obtained, building matching is to set.
Specifically, to fragment plane of disruption FiAnd FjIn characteristic point fiuAnd fjvMatched is exactly to calculate when different
Under domain, fiuAnd fjvThe distance between:
Wherein, M indicates fragment plane of disruption Fi, the number of K expression time domain, tkIndicate the value of k-th of time domain.
For each characteristic point fiuIf
Then by (fiu,fjv) as matching pair, to obtain fragment plane of disruption FiAnd FjMatching to set P.
The matching is preferably based on to set, carrying out splicing to described two fragment planes of disruption includes:
Based on the matching to the matching pair in set, all rigid body translation matrixes are calculated；
Each matching is calculated to the global error under all rigid body translation matrixes；
Using the rigid body translation matrix under the global error minimum as general transformation matrix；
Under the general transformation matrix, described two fragment planes of disruption are spliced.
Specifically, in order to keep the global error of splicing minimum, it can be disconnected to fragment by the constraint of rigid body translation matrix
Broken face FiAnd FjSpliced, the present embodiment calculates all rigid body translation matrixes using four-tuple method:
Calculate fragment plane of disruption FiAnd FjThe center of gravity of characteristic point:
Wherein,WithIt indicates to match the number to characteristic point in set, and number is identical.
By set of characteristic pointsWithConstruct covariance matrix:
Wherein, fikIndicate fragment plane of disruption FiIn k-th of characteristic point, fjkIndicate fragment plane of disruption FjIn, with fikIt is matched
Characteristic point, T representing matrix transposition.
Symmetrical matrix is constituted by covariance matrix:
Wherein, the mark of tr (R) representing matrix, I3Indicate 3 × 3 unit matrix, Δ=[R '23R′31R′12]T, R'=R-RT, under
Mark the element in 23,31 and 12 representing matrix R'.
The characteristic value and feature vector of Q (R) are calculated, vector corresponding to maximum eigenvalue is best rotating vector:
qR=[q0q1q2q3]T(12),
It should be noted that have the characteristic value and feature vector of calculating matrix in the prior art, and determine most preferably rotate to
The method of amount, which is not described herein again.
According to best rotating vector qR, calculate spin matrix:
Translation vector is calculated again:
Then rigid body translation matrix:
Calculate global error:
If the corresponding global error of some rigid body translation matrix is minimum, using the rigid body translation matrix as global transformation square
Battle array.
Under the constraint of general transformation matrix, to fragment plane of disruption FiAnd FjSpliced.
The matching is preferably based on to the matching pair in set, calculating all rigid body translation matrixes includes:
Based on the matching to the matching pair in set, three all not conllinear matchings pair are obtained；
Calculate the rigid body translation matrix of three all not conllinear matchings pair.
It should be noted that the accuracy in order to guarantee splicing, at least needs according to three not conllinear matchings to calculating
Rigid body translation matrix.
Embodiment two
As shown in figure 3, the present invention uses a kind of fragment plane of disruption splicing system based on thermonuclear feature, including characteristic point to obtain
Modulus block 10, thermonuclear feature calculation module 20, matching module 30 and splicing module 40；
Characteristic point obtains local curvature's maximum point that module 10 is used to obtain two fragment planes of disruption to be matched, by this
A little points are used as respective set of characteristic points；
Thermonuclear feature calculation module 20 is used to calculate the thermonuclear feature of characteristic point in the set of characteristic points；
Matching module 30 is for the thermonuclear feature based on the characteristic point to the feature of described two fragment planes of disruption
Point is matched, and building matching is to set；
Splicing module 40 is used to splice described two fragment planes of disruption set based on the matching.
Specifically, it is corresponding with method described in embodiment one, when system works, first passes through characteristic point and obtain module
10 obtain two fragment plane of disruption F to be matchediAnd FjThese are put and are used as respective characteristic point by upper local curvature's maximum point
Set；The thermonuclear feature of these local curvature's maximum of the calculating of thermonuclear feature calculation module 20；It is special that matching module 30 is based on thermonuclear
Sign is to fragment plane of disruption FiAnd FjIt is matched, building matching is to set P；Splicing module 40 is based on matching to set P, to fragment
Plane of disruption FiAnd FjSpliced.
Preferably, as shown in figure 4, matching module 30 includes metrics calculation unit 31, judging unit 32 and construction unit 33；
Metrics calculation unit 31 is used for the thermonuclear feature based on the characteristic point, calculates each in a fragment plane of disruption
Characteristic point is at a distance from all characteristic points in another fragment plane of disruption；
Judging unit 32 is used for judging whether the distance between two characteristic points is minimum, and will be apart from the smallest two characteristic points
As a matching pair；
Construction unit 33 is used to obtain the matching pair of each characteristic point in one fragment plane of disruption, building matching pair
Set.
Specifically, corresponding with method described in embodiment one, metrics calculation unit 31 is calculated in different time domains
Under, fiuAnd fjvThe distance between d (fiu,fjv)；Judging unit 32 is for judging d (fiu,fjv) whether minimum, it, will if minimum
(fiu,fjv) as a matching pair；Construction unit 33 by all matchings to be built into matching to set P.
Preferably, as shown in figure 5, splicing module 40 includes converter unit 41, error calculation unit 42,43 and of screening unit
Concatenation unit 44；
Converter unit 41 is used to calculate all rigid body translation matrixes to the matching pair in set based on the matching；
Error calculation unit 42 is for calculating each matching to the global error under all rigid body translation matrixes；
Screening unit 43 is for filtering out minimum global error, and by the rigid body translation under the global error minimum
Matrix is as general transformation matrix；
Concatenation unit 44 is used under the general transformation matrix, is spliced to described two fragment planes of disruption.
Specifically, corresponding with method described in embodiment one, converter unit 41 is used for based on the pairing set
Matching pair in conjunction calculates all rigid body translation matrix Φk；Error calculation unit 42 calculates rigid body translation matrix ΦkIt is corresponding
Global errorScreening unit 43 filters out minimum global error, and will be in global error minimum feelings
Rigid body translation matrix Φ under conditionkAs general transformation matrix；Concatenation unit 44 is in general transformation matrix ΦkUnder constraint, to fragment
Plane of disruption FiAnd FjSpliced.
The present invention calculates separately local curvature's maximum point of two fragment planes of disruption to be matched first, these points are made
For respective set of characteristic points, then calculate the thermonuclear feature of characteristic point in set of characteristic points, the thermonuclear feature pair based on characteristic point
The characteristic point of two fragment planes of disruption is matched, and building matching is spliced to set, and to two fragment planes of disruption, thus
Improve the accuracy of fragments mosaicing.
The foregoing is merely presently preferred embodiments of the present invention, is not intended to limit the invention, it is all in spirit of the invention and
Within principle, any modification, equivalent replacement, improvement and so on be should all be included in the protection scope of the present invention.
Claims (10)
1. a kind of fragment plane of disruption joining method based on thermonuclear feature, which comprises the steps of:
These are put and are used as respective characteristic point by the local curvature's maximum point for calculating separately two fragment planes of disruption to be matched
Set；
Calculate the thermonuclear feature of characteristic point in the set of characteristic points；
It is matched based on the characteristic point of the thermonuclear feature of the characteristic point to described two fragment planes of disruption, building matching
To set；
Based on the matching to set, described two fragment planes of disruption are spliced.
2. the method as described in claim 1, which is characterized in that the office for calculating separately two fragment planes of disruption to be matched
Portion's curvature maximum point, these, which are put the respective set of characteristic points of conduct, includes；
These are put and are used as respective initial spy by the local curvature's maximum point for calculating separately two fragment planes of disruption to be matched
Levy point set；
The characteristic point is subjected to triangle gridding；
Point in feature vertex neighborhood after the triangle gridding is added in the initial feature point set and is extended, is obtained
To set of characteristic points.
3. the method as described in claim 1, which is characterized in that the thermonuclear for calculating characteristic point in the set of characteristic points is special
Sign includes:
Calculate the Laplace-Beltrami operator of characteristic point in the set of characteristic points；
Laplace-Beltrami operator based on the characteristic point calculates the thermonuclear feature of characteristic point in the set of characteristic points.
4. method as claimed in claim 3, which is characterized in that calculate the characteristic point using grid Laplace operator
Laplace-Beltrami operator.
5. the method as described in claim 1, which is characterized in that it is described based on the thermonuclear feature of the characteristic point to described two
The characteristic point of the fragment plane of disruption is matched, and building matching includes: to set
Thermonuclear feature based on the characteristic point calculates each characteristic point and another fragment in a fragment plane of disruption and is broken
The distance of all characteristic points in face；
It will be apart from the smallest two characteristic points as a matching pair；
The matching pair of each characteristic point in one fragment plane of disruption is obtained, building matching is to set.
6. the method as described in claim 1, which is characterized in that described to be matched based on described to set, to described two fragments
The plane of disruption carries out splicing
Based on the matching to the matching pair in set, all rigid body translation matrixes are calculated；
Each matching is calculated to the global error under all rigid body translation matrixes；
Using the rigid body translation matrix under the global error minimum as general transformation matrix；
Under the general transformation matrix, described two fragment planes of disruption are spliced.
7. method as claimed in claim 6, which is characterized in that it is described to be matched based on described to the matching pair in set, it calculates
All rigid body translation matrixes include: out
Based on the matching to the matching pair in set, three all not conllinear matchings pair are obtained；
Calculate the rigid body translation matrix of three all not conllinear matchings pair.
8. a kind of fragment plane of disruption splicing system based on thermonuclear feature, which is characterized in that including characteristic point obtain module (10),
Thermonuclear feature calculation module (20), matching module (30) and splicing module (40)；
Characteristic point obtains local curvature's maximum point that module (10) are used to obtain two fragment planes of disruption to be matched, by these
Point is used as respective set of characteristic points；
Thermonuclear feature calculation module (20) is used to calculate the thermonuclear feature of characteristic point in the set of characteristic points；
Matching module (30) is for the thermonuclear feature based on the characteristic point to the characteristic point of described two fragment planes of disruption
It is matched, building matching is to set；
Splicing module (40) is used to splice described two fragment planes of disruption set based on the matching.
9. system as claimed in claim 8, which is characterized in that the matching module (30) include metrics calculation unit (31),
Judging unit (32) and construction unit (33)；
Metrics calculation unit (31) is used for the thermonuclear feature based on the characteristic point, and it is special to calculate each in a fragment plane of disruption
Sign point is at a distance from all characteristic points in another fragment plane of disruption；
Judging unit (32) will be made for judging whether the distance between two characteristic points is minimum apart from the smallest two characteristic points
For a matching pair；
Construction unit (33) is used to obtain the matching pair of each characteristic point in one fragment plane of disruption, constructs a pairing set
It closes.
10. system as claimed in claim 8, which is characterized in that the splicing module (40) includes converter unit (41), error
Computing unit (42), screening unit (43) and concatenation unit (44)；
Converter unit (41) is used to calculate all rigid body translation matrixes to the matching pair in set based on the matching；
Error calculation unit (42) is for calculating each matching to the global error under all rigid body translation matrixes；
Screening unit (43) is for filtering out minimum global error, and by the rigid body translation square under the global error minimum
Battle array is used as general transformation matrix；
Concatenation unit (44) is used under the general transformation matrix, is spliced to described two fragment planes of disruption.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910239060.1A CN110047151B (en) | 2019-03-27 | 2019-03-27 | Fragment fracture surface splicing method and system based on thermonuclear characteristics |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910239060.1A CN110047151B (en) | 2019-03-27 | 2019-03-27 | Fragment fracture surface splicing method and system based on thermonuclear characteristics |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110047151A true CN110047151A (en) | 2019-07-23 |
| CN110047151B CN110047151B (en) | 2023-05-23 |
Family
ID=67275333
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910239060.1A Active CN110047151B (en) | 2019-03-27 | 2019-03-27 | Fragment fracture surface splicing method and system based on thermonuclear characteristics |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110047151B (en) |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102945569A (en) * | 2012-10-23 | 2013-02-27 | 西北工业大学 | Three-dimensional model symmetry analysis method based on heat kernel signal |
| CN105427252A (en) * | 2015-11-06 | 2016-03-23 | 北京航空航天大学 | Mesh model multi-scale geometric detail repair method based on empirical mode decomposition |
| CN107590772A (en) * | 2017-07-18 | 2018-01-16 | 西北大学 | A kind of cultural relic fragments method for automatically split-jointing based on adaptive neighborhood matching |
-
2019
- 2019-03-27 CN CN201910239060.1A patent/CN110047151B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102945569A (en) * | 2012-10-23 | 2013-02-27 | 西北工业大学 | Three-dimensional model symmetry analysis method based on heat kernel signal |
| CN105427252A (en) * | 2015-11-06 | 2016-03-23 | 北京航空航天大学 | Mesh model multi-scale geometric detail repair method based on empirical mode decomposition |
| CN107590772A (en) * | 2017-07-18 | 2018-01-16 | 西北大学 | A kind of cultural relic fragments method for automatically split-jointing based on adaptive neighborhood matching |
Non-Patent Citations (1)
| Title |
|---|
| THOMAS CAISSARD: "heat kernel laplace-beltrami operator on digital surfaces", 《20TH INTERNATIONAL CONFERENCE ON DISCRETE GEOMETRY FOR COMPUTER IMAGERY》 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110047151B (en) | 2023-05-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106548159A (en) | Reticulate pattern facial image recognition method and device based on full convolutional neural networks | |
| CN104463899B (en) | A kind of destination object detection, monitoring method and its device | |
| Mizoguchi et al. | Quantitative scaling evaluation of concrete structures based on terrestrial laser scanning | |
| CN101315661B (en) | Fast three-dimensional face recognition method for reducing expression influence | |
| Xu et al. | Predicting animation skeletons for 3d articulated models via volumetric nets | |
| CN107766776A (en) | A kind of finger vein identification method and device | |
| CN106096503A (en) | A kind of based on key point with the three-dimensional face identification method of local feature | |
| CN110222642A (en) | A kind of planar architectural component point cloud contour extraction method based on global figure cluster | |
| CN110020652A (en) | The dividing method of Tunnel Lining Cracks image | |
| CN107025449A (en) | A kind of inclination image linear feature matching process of unchanged view angle regional area constraint | |
| Munsell et al. | Evaluating shape correspondence for statistical shape analysis: A benchmark study | |
| CN103268500A (en) | Gait identifying method with robustness to walking gait changes | |
| CN112446301B (en) | Method and system for detecting multi-class targets in remote sensing image | |
| CN106204422A (en) | Super large width image Rapid matching joining method based on block subgraph search | |
| CN109272442A (en) | Processing method, device, equipment and the storage medium of panorama spherical surface image | |
| US7181377B1 (en) | Method of modifying a volume mesh using sheet extraction | |
| CN101551865B (en) | Outward appearance and shape is utilized to detect the system and method with matching anatomical structure | |
| CN110047151A (en) | A kind of fragment plane of disruption joining method and system based on thermonuclear feature | |
| Chambolle et al. | Regularity for the optimal compliance problem with length penalization | |
| CN105046698A (en) | Shoulder joint defect parameter statistics method and system based on left and right symmetry information | |
| CN109544578A (en) | Recombination restored method based on the matched three-dimensional fragment of maximum similarity | |
| Wang et al. | Automated labeling of the airway tree in terms of lobes based on deep learning of bifurcation point detection | |
| CN108334873A (en) | A kind of 3D four-dimension hand data discrimination apparatus | |
| Du et al. | An automatic positioning algorithm for archaeological fragments | |
| Sun | Discrete Curvature Theories and Applications |
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 |