CN107993277A - Damage location artificial skelecton patch formation model method for reconstructing based on priori - Google Patents
Damage location artificial skelecton patch formation model method for reconstructing based on priori Download PDFInfo
- Publication number
- CN107993277A CN107993277A CN201711213243.3A CN201711213243A CN107993277A CN 107993277 A CN107993277 A CN 107993277A CN 201711213243 A CN201711213243 A CN 201711213243A CN 107993277 A CN107993277 A CN 107993277A
- Authority
- CN
- China
- Prior art keywords
- layer
- point
- repairing
- priori
- artificial skelecton
- 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
- 238000000034 method Methods 0.000 title claims abstract description 43
- 230000006378 damage Effects 0.000 title claims abstract description 33
- 230000015572 biosynthetic process Effects 0.000 title claims abstract description 31
- 210000000988 bone and bone Anatomy 0.000 claims abstract description 63
- 210000003625 skull Anatomy 0.000 claims abstract description 15
- 230000008439 repair process Effects 0.000 claims abstract description 12
- 238000012952 Resampling Methods 0.000 claims abstract description 8
- 230000000295 complement effect Effects 0.000 claims abstract description 7
- 238000003325 tomography Methods 0.000 claims abstract description 7
- 239000010410 layer Substances 0.000 claims description 141
- 230000008569 process Effects 0.000 claims description 9
- 230000009466 transformation Effects 0.000 claims description 8
- 230000005484 gravity Effects 0.000 claims description 6
- 238000005070 sampling Methods 0.000 claims description 5
- 238000004364 calculation method Methods 0.000 claims description 4
- 238000013459 approach Methods 0.000 claims description 3
- 239000011229 interlayer Substances 0.000 claims description 2
- 238000006243 chemical reaction Methods 0.000 abstract description 4
- 210000002454 frontal bone Anatomy 0.000 abstract description 3
- 210000003455 parietal bone Anatomy 0.000 abstract description 3
- 238000001356 surgical procedure Methods 0.000 abstract description 3
- 239000003814 drug Substances 0.000 abstract description 2
- 238000013467 fragmentation Methods 0.000 abstract description 2
- 238000006062 fragmentation reaction Methods 0.000 abstract description 2
- 208000010392 Bone Fractures Diseases 0.000 description 7
- 238000013461 design Methods 0.000 description 7
- 238000010586 diagram Methods 0.000 description 6
- 238000011160 research Methods 0.000 description 4
- 206010017076 Fracture Diseases 0.000 description 3
- 238000000605 extraction Methods 0.000 description 3
- PXFBZOLANLWPMH-UHFFFAOYSA-N 16-Epiaffinine Natural products C1C(C2=CC=CC=C2N2)=C2C(=O)CC2C(=CC)CN(C)C1C2CO PXFBZOLANLWPMH-UHFFFAOYSA-N 0.000 description 2
- 208000024779 Comminuted Fractures Diseases 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- 208000027418 Wounds and injury Diseases 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 210000005013 brain tissue Anatomy 0.000 description 1
- 125000004122 cyclic group Chemical group 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 229940079593 drug Drugs 0.000 description 1
- 238000003708 edge detection Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 208000014674 injury Diseases 0.000 description 1
- 230000003902 lesion Effects 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 210000000056 organ Anatomy 0.000 description 1
- 230000000717 retained effect Effects 0.000 description 1
- 238000013179 statistical model Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000013519 translation Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T15/00—3D [Three Dimensional] image rendering
- G06T15/005—General purpose rendering architectures
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/22—Matching criteria, e.g. proximity measures
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0012—Biomedical image inspection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30008—Bone
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2215/00—Indexing scheme for image rendering
- G06T2215/06—Curved planar reformation of 3D line structures
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Physics & Mathematics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Data Mining & Analysis (AREA)
- General Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Evolutionary Computation (AREA)
- Bioinformatics & Computational Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Artificial Intelligence (AREA)
- Life Sciences & Earth Sciences (AREA)
- Evolutionary Biology (AREA)
- General Health & Medical Sciences (AREA)
- Medical Informatics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Quality & Reliability (AREA)
- Computer Graphics (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Prostheses (AREA)
Abstract
The invention discloses a kind of damage location artificial skelecton patch formation model method for reconstructing based on priori.First, from medicine tomography CT images extract bone outermost contour, is introduced into repairing template and according to tomography thickness from repair template in resampling contour line；Secondly, alignd conversion using searching out template layer the most matched with layer to be repaired in method of characteristic repairing complementary modulus plate and establish both；Then, the point with layer breakpoint isomorphism to be repaired is obtained in template layer by public angle-off set and calculates the fit equation between isomorphism point, reversely substituted into fit equation and obtain layer interpolation point coordinates to be repaired；Then, internal layer profile is improved people's living condition into according to given bone ectonexine；Finally, artificial skelecton patch formation model is obtained using VTK three-dimensional reconstructions.The invention be applied to cause after the position such as human body frontal bone, parietal bone is inflicted heavy losses on the serious fragmentation of skull can not Direct Surgery reparation, it is necessary to which the actual operation for being implanted into personalized artificial skelecton operates in.
Description
Technical field
The invention belongs to computer graphics techniques field, especially suitable for repairing the artificial skelecton repairing of skull damage location
Model three dimensional design.
Background technology
It is suitable for the comminuted well damage of skull using artificial skelecton repairing skull damage location, and can not be directly by
Have in the clinical operation that broken bone splicing is completed.The design of artificial skelecton three-dimensional patch formation model has to and the position of lesion, size
It is consistent, and the situation of comminuted damage is directed to, it is artificial that doctor can not show that general shape is used as by the reparation to broken bone
The reference of bone design.For template frequently as auxiliary tool in three-dimensional reconstruction, template here can be existing healthy bone
Virtual Reconstruction result after reconstructed results or equalization.There are many documents to be directed to repairing using template auxiliary both at home and abroad
Affected bone profile is mended, such as automatic identification intersection point and breakpoint, is then reference using existing template, connects with reference to most short diagonal
Connection is by the fracture contour line reparation extracted into closed contour.Due to how much being had differences between template and current research case,
Therefore in order to ensure also introducing deformation mechanism in some undistorted documents of Extracting contour, such as in the three-dimensional reconstruction for vertebra
In, it is thus proposed that a set of statistical model is used to represent mean parameter as a result, template is passed through with current research case affine
Alignment, then deforms template, obtains the Real Three-dimensional Reconstruction result for meeting current research case.In addition, there are some texts
Offer during skull threedimensional model is obtained from template, constantly it is deformed using iteration, it is straight to cause to approach currently to grind
Study carefully case, it is final to obtain the complete skull threedimensional model with current research case characteristic.
Using existing same type bone collectively as priori, it is necessary to be found according to feature described function with current layer most
For matched same type bone contours line.Profile starting point registration is a kind of highly effective method, can by cyclic shift
Obtain most like in differently contoured starting point matching one.Alignment conversion can realize characteristic point image space coordinate transform,
Establish the correspondence between current layer and same type bone contours line.
The content of the invention
The present invention is to solve drawbacks described above and deficiency existing in the prior art, there is provided a kind of based on priori
The positions such as damage location artificial skelecton patch formation model method for reconstructing, frontal bone, parietal bone for human body skull occur after being inflicted heavy losses on
Comminuted fracture, skull fragmentation seriously can not Direct Surgery reparation situation, using different human body skull bone have the overall situation
Similitude；There is sequentiality, local similarity between Medical Slice Images；The bone contours of extraction have the features such as continuity, will
The bone that has been averaged obtains the bone threedimensional model of fracture location as priori automatically, is repaiied for post personalization artificial skelecton
Mend operation and data basis is provided.Compared with existing broken bone splicing, all it is improved in the accuracy and timeliness of repairing.
In order to solve the above technical problems, the present invention provides a kind of damage location artificial skelecton repairing mould based on priori
Type method for reconstructing, is guided using having repaired template contours as priori, and automatic obtain is used to repair skeleton pars affecta
The artificial skelecton patch formation model of position, detailed process are as follows：
(1) bone contours extracted from Medical Slice Images are pre-processed, and searches for the bone of all not closed faults
Bone contour line, the step include as follows：
Step S01：Medical Slice Images are read in successively, and the bone contours of Medical Slice Images are obtained using Canny methods
Map, some closures of every Medical Slice Images or inc bone contours line are obtained with reference to edge following algorithm；
Step S02：There is layer where not closing bone contours line for the first time in record, and is labeled as A, similarly will last time
Layer where appearance does not close bone contours line is labeled as Z；
(2) repairing template and the correspondence of not closed contour are established, which includes as follows：
Step S03：The coordinate of the recording layer A points of interruption P, Q；
Step S04：Setting repairing template contour interval, resampling obtain one group of closure bone contours；
Step S05：Found in repairing complementary modulus plate and the most matched repairing template layer A ' -1 of closed contour layer A-1；
(3) to repair template as priori, obtain and be used for the artificial skelecton for repairing skeleton damage location at tomography
Contour curve, the step include as follows：
Step S06：The Optimum Matching being calculated according to step S05, acquisition and the most matched Optimum Matching layer A ' of layer A,
Both alignment transformation equations of structure, then using the transformation equation that aligns, make registration process by layer A, A '；
Step S07：Using public angle, the point P ' coordinates with breakpoint P-isomorphism are tried to achieve in Optimum Matching layer A ', are similarly asked
Point Q ' the coordinates of the isomorphism of another breakpoint Q are obtained, sequence of calculation point P ' arrives the fit equation X of Q '；
Step S08：Interpolation point coordinates in the fit equation X computation layers A obtained using step S07；
Step S09：According to given bone thickness d, the internal layer profile coordinate points of calculating artificial skelecton profile；
Step S10：Layer A+1 is denoted as A；
Step S11：Repeat step S03 to step S11, until A and Z is overlapped.
(4) three-dimensional reconstruction artificial skelecton patch formation model, the step include as follows：
Step S12：Artificial skelecton patch formation model is obtained using VTK three-dimensional reconstructions.
Wherein, in step S01, closure or inc bone contours line refer to that the bone contours line on tomography is in circuit
Or the state of open circuit.
Wherein, in step S02, there is layer where not closing bone contours line for the first time in record, and is labeled as A, similarly will
There is layer where not closing bone contours line and is labeled as Z in last time, it is concretely comprised the following steps：
Step i01：Searched for since top, the degree d (v) of every layer of upper each vertex v is calculated in a manner of nonoriented edge, if d
(v)<2 represent that the point is hitch point, which not closing phenomenon occurs, labeled as A；
Step i02：Searched for since the bottom, the degree d (v) of every layer of upper each vertex v is calculated in a manner of nonoriented edge, if d
(v)<2 represent that the point is hitch point, which not closing phenomenon occurs, labeled as Z.
Wherein, in step S03, the coordinate of the recording layer A points of interruption P, Q, it is concretely comprised the following steps：
Step q01：Searched for since top and position layer where to the contour line labeled as A；
Step q02：It is P, Q to mark by vertex of the degree of vertex in this layer less than 2.
Wherein, in step S04, repairing template refer to the average skeleton data identical with Medical faultage image position and
The average skeleton data can represent complete whole characteristic informations when damage location；Contour refers to scan Medical Slice Images
Spacing of the spacing as repairing template resampling；The resampling refers to repairing template according to Medical Slice Images scanning room
Sampled away from, direction, the bone contours line of acquisition is known as repairing template layer.
Wherein, in step S05, found in repairing complementary modulus plate with the most matched repairing template layer A ' of closed contour layer A-1-
1, it is concretely comprised the following steps：
Step p01：Positioning is to above layer A and the closed contour layer A-1 adjacent with layer A；
Step p02：Make similitude judgement in layer A-1 and repairing template interlayer using method of characteristic, so as to find and layer A-1
Most matched repairing template layer, labeled as A ' -1；Otherwise, first sampling point is adjusted, step S04 is performed again and obtains new sampling
Data.
Wherein, in step S06, using public angle, the point P ' coordinates with breakpoint P-isomorphism are tried to achieve in layer A ', its is specific
Step is：
Step l01：With R, R ' expression layers A, A respectively ' barycentric coodinates, P, Q are the A points of interruption；
Step l02：Calculate the angle theta that the line PR from P to center of gravity R is formed with horizontal line；
Step l03：From R ' do one with horizontal line angle be all θ ray, ray and layer A ' intersect at point P ', then
P ' is the point with breakpoint P-isomorphism.
Wherein, in step S07, fit equation refers to utilize least square curve fit principle, by the number on discrete point
According to collection (P ', P+1 ', P+2 ' ..., P+i ' ..., Q-j ' ..., Q-2 ', Q-1 ', Q ', wherein i<J) analytical function is configured to, and is made
The curve infinite approach discrete point of the function.
Wherein, in step S08, using interpolation point coordinates in fit equation X computation layers A, it is concretely comprised the following steps：
Step w01：N equally distributed point x are inserted on the X-axis coordinate of P and Q1,x2,…,xn；
Step w02：Bring the value of x into fit equation X, y is calculated1,y2,…,ynValue.This n coordinate points (xi,yi)
As artificial skelecton profile point.
Wherein, in step S09, according to given bone thickness d, internal layer profile coordinate points are calculated, it is concretely comprised the following steps：
Step h01：It is respectively the coordinate points on outer layer and internal layer profile to make (x, y), (x ', y '), and (x, y) it is known that need
Ask (x ', y ')；
Step h02：According to thickness average value d inside and outside every layer of skull of profile geographic survey；
Step h03：It is respectively x-x ', y-y ' to make dx, dy；
Step h04：Bring d into equationsWherein make dx=dy, you can coordinate (x ', y ') is calculated
Value；
Step h05：Repeat step h01- step h04, calculate and obtain every layer of all internal layer artificial skelecton profile point coordinate value.
The advantageous effects that the present invention is reached：The present invention provides a kind of damage location based on priori is artificial
There is comminuted fracture, head after being inflicted heavy losses in the position such as skeletal repair model reconstruction method, frontal bone, parietal bone for human body skull
Osteomiosis seriously can not Direct Surgery reparation situation, there is global similitude using different human body skull bone；Medicine is broken
There is sequentiality, local similarity between tomographic image；The bone contours of extraction have the features such as continuity, and the bone that will be averaged is made
Obtain the bone threedimensional model of fracture location automatically for priori, data are provided for post personalization artificial skelecton repairing operation
Basis.Compared with existing broken bone splicing, all it is improved in the accuracy and timeliness of repairing.
Brief description of the drawings
Fig. 1 is the damage location artificial skelecton patch formation model method for reconstructing workflow based on priori in the present invention
Figure；
Fig. 2 is that Canny algorithms extraction bone contours process schematic is utilized in the present invention；
Fig. 3 is the schematic diagram of the degree on vertex in the present invention；
Fig. 4 is current layer and repairing template layer position relationship schematic diagram in the present invention；
Fig. 5 is search Optimum Matching algorithm flow chart in the present invention；
Fig. 6 is Optimum Matching and its feature described function schematic diagram in the present invention；
Fig. 7 is to make alignment conversion schematic diagram in the present invention to repairing template on the basis of current layer；
Fig. 8 is to obtain matched curve schematic diagram using isomorphism point in the present invention；
Fig. 9 is generation internal layer outline drawing in the present invention；
Figure 10 is to obtain damage location artificial skelecton repairing reconstruction model to repair template as priori in the present invention and show
It is intended to.
Embodiment
With reference to specific embodiment, the invention will be further described.Following embodiments are only used for clearly illustrating
Technical scheme, and be not intended to limit the protection scope of the present invention and limit the scope of the invention.
Patent of the present invention is further illustrated with reference to the accompanying drawings and examples.
As shown in Figure 1, the damage location artificial skelecton patch formation model method for reconstructing based on priori, it is characterised in that：
Guided using template contours have been repaired for priori, it is automatic to obtain for the position that skeleton degree of injury is larger
Process for the artificial skelecton patch formation model of repairing.
First, the bone contours that all Medical Slice Images extract are done to the bone of all not closed faults of pretreatment search
Bone contour line, the step include as follows：
Step S01：Medical Slice Images are read in successively, and the bone contours of Medical Slice Images are obtained using Canny methods
Map, some closures of every Medical Slice Images or inc bone contours are obtained with reference to edge following algorithm；
Step S02：There is layer where not closing bone contours line for the first time in record, and is labeled as A, similarly will last time
Layer where appearance does not close bone contours line is labeled as Z；
Secondly, repairing template and the correspondence of current not closed contour are established, which includes as follows：
Step S03：Recording layer A interrupts point coordinates P, Q；
Step S04：Setting repairing template contour interval, resampling obtain one group of closure bone contours；
Step S05：The most matched repairing template layers of the closed contour layer A-1 adjacent with layer A are found in repairing complementary modulus plate
A’-1；
Then, to repair template as priori, obtain this tomography be used for repair artificial skelecton contour curve, the step
Including as follows：
Step S06：The Optimum Matching being calculated according to step S05, it is known that with layer A Optimum Matchings layer it is A ', structure two
Person's alignment transformation equation, then using the transformation equation that aligns, makees registration process by layer A, A '；
Step S07：Using public angle, the point P ' coordinates with breakpoint P-isomorphism are tried to achieve in layer A ', similarly try to achieve another
Point Q ' the coordinates of the isomorphism of breakpoint Q, sequence of calculation point P ' arrive the fit equation X of Q '；
Step S08：Interpolation point coordinates in the fit equation X computation layers A obtained using step S07；
Step S09：According to given bone thickness d, internal layer profile coordinate points are calculated；
Step S10：Layer A+1 is denoted as A；
Step S11：Repeat step S03 to step S09, until A and Z is overlapped.
Finally, three-dimensional reconstruction artificial skelecton patch formation model, the step include as follows：
Step S12：Artificial skelecton patch formation model is obtained using VTK three-dimensional reconstructions.
The process of Medical Slice Images outermost contour line is obtained as shown in Figure 2.Fig. 2-1. it is Medical Slice Images；Fig. 2-2.
After setting binary-state threshold to Medical Slice Images, after having filtered out the organ pixels such as skin in artwork, brain tissue, eye ear, protect
The bone pixel stayed；Fig. 2-it is 3. the bone contours map obtained using Canny edge detection methods, including bone internal layer wheel
It is wide；Fig. 2-4. it is to obtain outer layer bone contours using edge following algorithm.
As shown in figure 3, the side number that degree d (v) refers to and the vertex v is associated.As one embodiment of the present of invention, such as scheme
1. 3- show the 25th layer of bone outermost contour line, which is to occur not closed contour place layer, its midpoint v for the first timejFor even
Continuous point, degree d (vj)=2；Point vi、vkRespectively breakpoint, degree d (vi)=d (vk)=2；Such as Fig. 3-2. show outside the 74th layer of bone
Layer contour line, dotted line frame represent partial enlargement, its midpoint vm、vnRespectively crosspoint, degree d (vm)=d (vn)=3.
As shown in figure 4, represent current layer and repairing template layer position relationship.3. wherein Fig. 4-1., 4- represent first respectively
Layer A (25) and its adjacent closed layer A-1 (24) where secondary appearance does not close bone contours line；4. Fig. 4-2., 4- are represented respectively
Repair template A ' (78) layers and A ' -1 (77) layer.
As shown in figure 5, found in repairing complementary modulus plate with most matched repairing template layer A ' -1 of layer A-1, it is concretely comprised the following steps：
Step p01：Positioning is to above layer A and closed contour layer A-1 adjacent thereto；
Step p02：Make similitude judgement between layer A-1 and repairing template using method of characteristic, so as to find with layer A-1 most
Matched repairing template layer, labeled as A ' -1；Otherwise, first sampling point is adjusted, step S04 is performed again and obtains new hits
According to.
As shown in fig. 6, Fig. 6-1. it is layer A-1；Fig. 6-2. represent in repairing template with the most matched layer A '-1 of layer A-1；Figure
6- is 3., Fig. 6-4. is utilizes the feature described function of method of characteristic computation layer A-1, repairing masterplate layer A '-1：fM(n) and f(A-1)
(n), it is concretely comprised the following steps：
Step u01：Mean error minimum method is utilized to polymerize the profile point repaired template and repaired in template non-significant
Characteristic point, identical quantative attribute point is retained using threshold value in two templates；
Step u02：Both, which are obtained, using method of characteristic matches starting point；
Step u03：It can be found if it can obtain Optimum Matching equidirectional with adjacent non-repairing template contours the most
Matched repairing template contours；Otherwise adjustment repairing template contour first sampling point, continues step S04；
Step u04：Repeat step t01 to step t03, find it is all can with it is non-repairing template matches repairing template.
As shown in fig. 7, utilize alignment conversion aligned layer A and repairing template layer A '.Fig. 7-1. it is layer A '-1；Fig. 7-2. it is to incite somebody to action
Layer A ' -1 rotates clockwise result figure after α centered on center of gravity R '；Fig. 7-3. it is that layer A '-1 is translated into d (dx, dy) results afterwards
Figure；Fig. 7-and 4. it is layer A, wherein center of gravity R is that the vertical of-1 barycentric coodinates R of layer A ' maps；Fig. 7-5. it is to be with center of gravity R ' by layer A '
Center, equally rotates clockwise result figure after α；Fig. 7-6. it is that layer A ' is translated into d (dx, dy) result figures afterwards.Wherein rotate, translate
Concretely comprise the following steps：
Step k01：Computation layer A-1, A ' -1 barycentric coodinates：R, R '；
Step k02：The matching starting point that record is obtained according to step u02 is to coordinate：S, S '；
Step k03：Calculate line segment RS and horizontal line angle theta and R ' S ' and horizontal line angle theta '；
Step k04：Displacement difference dx, dy between calculating R, R '；
Step k04：α=θ-θ ' are brought into spin matrix, dx, dy are brought into translation matrix, obtain alignment transformation equation.
As shown in figure 8, represent to obtain matched curve using isomorphism point.
The affine alignment of Fig. 8-1. represent layer A, A ' is as a result, wherein black is not close layer A where bone contours line, grey
To repair a certain layer A ' in template, and layer A ' is matched the most with layer A.P, Q is the A points of interruption, P ', Q ' be in A ' with P, Q isomorphism
Point, its calculating concretely comprise the following steps：
Step l01：With R, R ' expression layers A, A respectively ' barycentric coodinates, P, Q are the A points of interruption；
Step l02：Calculate the line PR from P to center of gravity R and form angle theta with horizontal line；
Step l03：Sent from R ' do one with horizontal line angle be all θ ray, ray and layer A ' intersect at point P ', then
P ' is the point with breakpoint P-isomorphism.
Fig. 8-2. middle vertex P ', Q ' half interval contours are matched curve X, and concrete methods of realizing is：By the data set on discrete point
(P ', P+1 ', P+2 ' ..., P+i ' ..., Q-j ' ..., Q-2 ', Q-1 ', Q ', wherein i<J) it is configured to analytical function so that the letter
Number curve is as close possible to former discrete point.
The middle curve PQ of Fig. 8-3. is the coordinate for going out interpolation point in layer A by the equation calculation of matched curve X, its specific steps
For：
Step w01：N equally distributed point x are inserted on the X-axis coordinate of P and Q1,x2,…,xn；
Step w02：Bring the value of x into fit equation X, y is calculated1,y2,…,ynValue.This n coordinate points (xi,yi)
As repair artificial skelecton profile point.
As shown in figure 9, according to given bone thickness d, internal layer profile coordinate points are calculated.
Fig. 9-1. middle sequence of points P ..., P+m ..., Q-n ..., Q, m<N is according to step w01- step w02, passes through fitting
The interpolation point set that equation obtains；
Fig. 9-and 2. it is according to existing outermost contour coordinate points and given bone thickness, the internal layer profile being calculated is sat
The result diagram of punctuate, its calculating concretely comprise the following steps：
Step h01：It is respectively the coordinate points on outer layer and internal layer profile to make (x, y), (x ', y '), and (x, y) it is known that need
Ask (x ', y ')；
Step h02：According to thickness average value d inside and outside this layer of skull of profile geographic survey；
Step h03：It is respectively x-x ', y-y ' to make dx, dy；
Step h04：Bring d into equationsDx=dy is wherein made, the value of coordinate (x ', y ') is calculated；
Step h05：Repeat step h01- step h04, calculate and obtain this layer of all internal layer repairing profile point coordinate value.
Fig. 9-3. it is the design sketch being embedded into obtained artificial skelecton ectonexine contour line in layer A.
Priori, the process that damage location artificial skelecton patch formation model is rebuild are based on below by the method for the present invention.
Embodiment one：
As shown in Figure 10, Figure 10-1. is damage skull three-dimensional reconstruction design sketch；Figure 10-2. it is to be obtained by the method for the present invention
The artificial skelecton patch formation model three-dimensional reconstruction design sketch arrived；Figure 10-3., Figure 10-is 4. respectively that artificial skelecton patch formation model is embedding
Enter the different visual angles design sketch of damage location.
The present invention is disclosed with preferred embodiment above, so it is not intended to limiting the invention, all to take equivalent substitution
Or the technical solution that the scheme of equivalent transformation is obtained, all fall within protection scope of the present invention.
Claims (10)
1. the damage location artificial skelecton patch formation model method for reconstructing based on priori, it is characterised in that：To have repaired mould
Web wheel is wide to be guided for priori, automatic to obtain the artificial skelecton patch formation model for being used for repairing skeleton damage location, tool
Body process is as follows：
(1) bone contours extracted from Medical Slice Images are pre-processed, and search for all damage locations, that is, the wheel obtained
Exterior feature is fracture shape.It is inc bone contours line to define this contour line for occurring damaging, which includes as follows：
Step S01：Medical Slice Images are read in successively, and the bone contours map of Medical Slice Images is obtained using Canny methods,
Some closures of every Medical Slice Images or inc bone contours line are obtained with reference to edge following algorithm；
Step S02：There is layer where not closing bone contours line for the first time in record, and is labeled as A, similarly last time occurs
Layer where not closing bone contours line is labeled as Z；
(2) repairing template and the correspondence of not closed contour are established, which includes as follows：
Step S03：The coordinate of the recording layer A points of interruption P, Q；
Step S04：Setting repairing template contour interval, resampling obtain one group of closure bone contours；
Step S05：Found in repairing complementary modulus plate and closed contour layer A-1Most matched repairing template layer A '-1；
(3) to repair template as priori, obtain and be used for the artificial skelecton repairing for repairing skeleton damage location at tomography
Model silhouette curve, the step include as follows：
Step S06：The Optimum Matching being calculated according to step S05, obtains and the most matched Optimum Matching layer A ' of layer A, structure
Both alignment transformation equations, then using the transformation equation that aligns, make registration process by layer A, A '；
Step S07：Using public angle, the point P ' coordinates with breakpoint P-isomorphism are tried to achieve in Optimum Matching layer A ', are similarly tried to achieve another
Point Q ' the coordinates of the isomorphism of one breakpoint Q, sequence of calculation point P ' arrive the fit equation X of Q '；
Step S08：Interpolation point coordinates in the fit equation X computation layers A obtained using step S07；
Step S09：According to given bone thickness d, the internal layer profile coordinate points of calculating artificial skelecton patch formation model contour curve；
Step S10：By layer A+1It is denoted as A；
Step S11：Repeat step S03 to step S11, until A and Z is overlapped.
(4) three-dimensional reconstruction artificial skelecton patch formation model, the step include as follows：
Step S12：Artificial skelecton patch formation model is obtained using VTK three-dimensional reconstructions.
2. the damage location artificial skelecton patch formation model method for reconstructing according to claim 1 based on priori, it is special
Sign is：Closure or inc bone contours line refer to that the bone contours line on tomography is in circuit or disconnected in step S01
The state on road.
3. the damage location artificial skelecton patch formation model method for reconstructing according to claim 1 based on priori, it is special
Sign is：In step S02, there is layer where not closing bone contours line for the first time in record, and is labeled as A, similarly by last
Layer where secondary appearance does not close bone contours line is labeled as Z, it is concretely comprised the following steps：
Step i01：Searched for since top, the degree d (v) of every layer of upper each vertex v is calculated in a manner of nonoriented edge, if d (v)<2
It is hitch point to represent the point, which not closing phenomenon occurs, labeled as A；
Step i02：Searched for since the bottom, the degree d (v) of every layer of upper each vertex v is calculated in a manner of nonoriented edge, if d (v)<2
It is hitch point to represent the point, which not closing phenomenon occurs, labeled as Z.
4. the damage location artificial skelecton patch formation model method for reconstructing according to claim 1 based on priori, it is special
Sign is：In step S03, the coordinate of the recording layer A points of interruption P, Q, it is concretely comprised the following steps：
Step q01：Searched for since top and position layer where to the contour line labeled as A；
Step q02：It is P, Q to mark by vertex of the degree of vertex in this layer less than 2.
5. the damage location artificial skelecton patch formation model method for reconstructing according to claim 1 based on priori, it is special
Sign is：In step S04, repairing template refers to the average skeleton data identical with Medical faultage image position and this is average
Skeleton data can represent complete whole characteristic informations when damage location；Contour refers to make Medical Slice Images sweep span
To repair the spacing of template resampling；The resampling refers to repairing template according to Medical Slice Images sweep span, direction
Sampled, the bone contours line of acquisition is known as repairing template layer.
6. the damage location artificial skelecton patch formation model method for reconstructing according to claim 1 based on priori, it is special
Sign is：In step S05, found in repairing complementary modulus plate and most matched repairing template layer A ' -1 of closed contour layer A-1, its tool
Body step is：
Step p01：Positioning is to above layer A and the closed contour layer A-1 adjacent with layer A；
Step p02：Make similitude judgement in layer A-1 and repairing template interlayer using method of characteristic, so as to find with layer A-1 most
The repairing template layer matched somebody with somebody, labeled as A ' -1；Otherwise, first sampling point is adjusted, step S04 is performed again and obtains new hits
According to.
7. the damage location artificial skelecton patch formation model method for reconstructing according to claim 1 based on priori, it is special
Sign is：In step S06, using public angle, the point P ' coordinates with breakpoint P-isomorphism are tried to achieve in layer A ', it is concretely comprised the following steps：
Step l01：With R, R ' expression layers A, A respectively ' barycentric coodinates, P, Q are the A points of interruption；
Step l02：Calculate the angle theta that the line PR from P to center of gravity R is formed with horizontal line；
Step l03：From R ' do one with horizontal line angle be all θ ray, ray and layer A ' intersect at point P ', then P ' is i.e.
For the point with breakpoint P-isomorphism.
8. the damage location artificial skelecton patch formation model method for reconstructing according to claim 1 based on priori, it is special
Sign is：In step S07, fit equation refers to utilize least square curve fit principle, by the data set on discrete point
(P ', P+1 ', P+2 ' ..., P+i ' ..., Q-j ' ..., Q-2 ', Q-1 ', Q ', wherein i<J) analytical function is configured to, and makes the letter
Several curve infinite approach discrete points.
9. the damage location artificial skelecton patch formation model method for reconstructing according to claim 1 based on priori, it is special
Sign is：In step S08, using interpolation point coordinates in fit equation X computation layers A, it is concretely comprised the following steps：
Step w01：N equally distributed point x are inserted on the X-axis coordinate of P and Q1,x2,…,xn；
Step w02：Bring the value of x into fit equation X, y is calculated1,y2,…,ynValue.This n coordinate points (xi,yi) be
Point on artificial skelecton contour curve.
10. the damage location artificial skelecton patch formation model method for reconstructing according to claim 9 based on priori, it is special
Sign is：In step S09, according to given bone thickness d, internal layer profile coordinate points are calculated, it is concretely comprised the following steps：
Step h01：It is respectively the coordinate points on outer layer and internal layer profile to make (x, y), (x ', y '), and (x, y) it is known that demand
(x’,y’)；
Step h02：According to thickness average value d inside and outside every layer of skull of profile geographic survey；
Step h03：It is respectively x-x ', y-y ' to make dx, dy；
Step h04：Bring d into equationsWherein make dx=dy, you can the value of coordinate (x ', y ') is calculated；
Step h05：Repeat step h01- step h04, calculate and obtain every layer of all internal layer artificial skelecton profile point coordinate value.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711213243.3A CN107993277B (en) | 2017-11-28 | 2017-11-28 | priori knowledge-based reconstruction method of artificial bone repair model of damaged part |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711213243.3A CN107993277B (en) | 2017-11-28 | 2017-11-28 | priori knowledge-based reconstruction method of artificial bone repair model of damaged part |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107993277A true CN107993277A (en) | 2018-05-04 |
| CN107993277B CN107993277B (en) | 2019-12-17 |
Family
ID=62032370
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201711213243.3A Expired - Fee Related CN107993277B (en) | 2017-11-28 | 2017-11-28 | priori knowledge-based reconstruction method of artificial bone repair model of damaged part |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107993277B (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109887077A (en) * | 2019-03-07 | 2019-06-14 | 百度在线网络技术（北京）有限公司 | Method and apparatus for generating threedimensional model |
| WO2020014889A1 (en) * | 2018-07-18 | 2020-01-23 | 中山大学附属第一医院 | Method for establishing bionic prosthesis model for limb trunk nerves |
| CN110766666A (en) * | 2019-10-09 | 2020-02-07 | 河海大学常州校区 | Skeleton fault outline large sample library generation method based on region segmentation and GAN model |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6978188B1 (en) * | 2002-09-30 | 2005-12-20 | Medical Modeling, Llc | Method for contouring bone reconstruction plates |
| CN101882326A (en) * | 2010-05-18 | 2010-11-10 | 广州市刑事科学技术研究所 | Three-dimensional craniofacial reconstruction method based on overall facial structure shape data of Chinese people |
| CN103310432A (en) * | 2013-06-25 | 2013-09-18 | 西安电子科技大学 | Computerized Tomography (CT) image uniformization metal artifact correction method based on four-order total-variation shunting |
| CN104809736A (en) * | 2015-05-12 | 2015-07-29 | 河海大学常州校区 | Medical tomographic image closed skeleton outline calculation method based on priori knowledge |
| CN104881875A (en) * | 2015-06-09 | 2015-09-02 | 河海大学常州校区 | Priori-knowledge-based reconstruction method of closed skeleton contour in medical fault image |
| CN105184861A (en) * | 2015-10-20 | 2015-12-23 | 南京大学 | Automatic human three-dimensional (3D) reconstruction system |
| CN105816232A (en) * | 2016-05-17 | 2016-08-03 | 南方医科大学 | Designing and forming method for anatomical bone plate of individuation bone model |
| CN105982722A (en) * | 2015-02-02 | 2016-10-05 | 北京大璞三维科技有限公司 | Digital skeleton operation repair method and system |
-
2017
- 2017-11-28 CN CN201711213243.3A patent/CN107993277B/en not_active Expired - Fee Related
Patent Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6978188B1 (en) * | 2002-09-30 | 2005-12-20 | Medical Modeling, Llc | Method for contouring bone reconstruction plates |
| CN101882326A (en) * | 2010-05-18 | 2010-11-10 | 广州市刑事科学技术研究所 | Three-dimensional craniofacial reconstruction method based on overall facial structure shape data of Chinese people |
| CN103310432A (en) * | 2013-06-25 | 2013-09-18 | 西安电子科技大学 | Computerized Tomography (CT) image uniformization metal artifact correction method based on four-order total-variation shunting |
| CN105982722A (en) * | 2015-02-02 | 2016-10-05 | 北京大璞三维科技有限公司 | Digital skeleton operation repair method and system |
| CN104809736A (en) * | 2015-05-12 | 2015-07-29 | 河海大学常州校区 | Medical tomographic image closed skeleton outline calculation method based on priori knowledge |
| CN104809736B (en) * | 2015-05-12 | 2017-09-12 | 河海大学常州校区 | Medical Slice Images closure bone contours computational methods based on priori |
| CN104881875A (en) * | 2015-06-09 | 2015-09-02 | 河海大学常州校区 | Priori-knowledge-based reconstruction method of closed skeleton contour in medical fault image |
| CN105184861A (en) * | 2015-10-20 | 2015-12-23 | 南京大学 | Automatic human three-dimensional (3D) reconstruction system |
| CN105816232A (en) * | 2016-05-17 | 2016-08-03 | 南方医科大学 | Designing and forming method for anatomical bone plate of individuation bone model |
Non-Patent Citations (10)
| Title |
|---|
| P. JIANG ET AL.: ""Macroporous/mesoporous bioglasses doped with Ag/TiO2 for dual drug action property and bone repair"", 《IET NANOBIOTECHNOLOGY》 * |
| 丁焕文 等: ""计算机辅助股骨粗隆部肿瘤精确切除和保髋重建的临床研究"", 《中华关节外科杂志(电子版)》 * |
| 孙坚 等: ""应用SimPlant软件精确匹配与缺损外形吻合的髂骨肌瓣行下颌骨重建"", 《组织工程与重建外科杂志》 * |
| 张永弟 等: ""提高CT图像手骨模型重建精度的方法"", 《计算机辅助设计与图形学学报》 * |
| 张霞 等: ""基于熔融沉积成型技术的骨骼成型优化方法"", 《计算机应用研究》 * |
| 李昭 等: ""自体骨植骨修复膝关节置换术中胫骨缺损的临床研究"", 《河北医科大学学报》 * |
| 王琳 等: ""基于模板的骨骼参数自动测量方法"", 《计算机科学》 * |
| 王琳 等: ""基于特征点映射关系的接骨板系列化设计方法"", 《计算机辅助设计与图形学学报》 * |
| 陈小中 等: ""参数化股骨曲面特征模型的构建"", 《计算机辅助设计与图形学学报》 * |
| 高静 等: ""基于模板叠加相异特征匹配的CT图像断点标定"", 《科技通报》 * |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2020014889A1 (en) * | 2018-07-18 | 2020-01-23 | 中山大学附属第一医院 | Method for establishing bionic prosthesis model for limb trunk nerves |
| CN109887077A (en) * | 2019-03-07 | 2019-06-14 | 百度在线网络技术（北京）有限公司 | Method and apparatus for generating threedimensional model |
| CN109887077B (en) * | 2019-03-07 | 2022-06-03 | 百度在线网络技术（北京）有限公司 | Method and apparatus for generating three-dimensional model |
| CN110766666A (en) * | 2019-10-09 | 2020-02-07 | 河海大学常州校区 | Skeleton fault outline large sample library generation method based on region segmentation and GAN model |
| CN110766666B (en) * | 2019-10-09 | 2022-11-04 | 河海大学常州校区 | Skeleton fault outline large sample library generation method based on region segmentation and GAN model |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107993277B (en) | 2019-12-17 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110459301A (en) | Brain neuroblastoma surgical navigation method for registering based on thermodynamic chart and facial key point | |
| CN102880866B (en) | Method for extracting face features | |
| CN107993277A (en) | Damage location artificial skelecton patch formation model method for reconstructing based on priori | |
| CN108053433A (en) | A kind of multi-modal arteria carotis MRI method for registering based on physical alignment and outline | |
| Iqtait et al. | Feature extraction for face recognition via active shape model (ASM) and active appearance model (AAM) | |
| CN101339670B (en) | Computer auxiliary three-dimensional craniofacial rejuvenation method | |
| CN106469465A (en) | A kind of three-dimensional facial reconstruction method based on gray scale and depth information | |
| TW201109001A (en) | Method of fabricating artificial implant | |
| CN106887039A (en) | The three-D imaging method and system of a kind of organ and focus based on medical image | |
| CN110197462A (en) | A kind of facial image beautifies in real time and texture synthesis method | |
| CN102930602B (en) | Tomography-image-based facial skin three-dimensional surface model reconstructing method | |
| Akbari et al. | 3D segmentation of prostate ultrasound images using wavelet transform | |
| CN105849777B (en) | A kind of cerebral cortex electrode merges the method and device of display with magnetic resonance image | |
| CN106898044B (en) | Organ splitting and operating method and system based on medical images and by utilizing VR technology | |
| CN104240180A (en) | Method and device for achieving automatic adjusting of images | |
| CN110175529A (en) | A kind of three-dimensional face features' independent positioning method based on noise reduction autoencoder network | |
| CN109753997A (en) | A kind of liver neoplasm automatic and accurate Robust Segmentation method in CT image | |
| CN107240128A (en) | A kind of X-ray film and photochrome method for registering based on contour feature | |
| Arora et al. | 3D fingerprint phantoms | |
| CN108961278A (en) | The method and its system of abdominal wall muscle segmentation based on image data | |
| Tao et al. | Statistical study on cortical sulci of human brains | |
| Guerroudji et al. | Automatic brain tumor segmentation, and 3d reconstruction and visualization using augmented reality | |
| CN106845138A (en) | Method is previewed before a kind of surgery | |
| Heimann et al. | Model-based segmentation | |
| CN106023094A (en) | Image-based bone tissue microstructure restoration system and restoration method thereof |
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 | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20191217 |