CN100569195C - Manufacturing method of fine personalized skull model capable of describing teeth occluding relation - Google Patents

Manufacturing method of fine personalized skull model capable of describing teeth occluding relation Download PDF

Info

Publication number
CN100569195C
CN100569195C CNB2006101131320A CN200610113132A CN100569195C CN 100569195 C CN100569195 C CN 100569195C CN B2006101131320 A CNB2006101131320 A CN B2006101131320A CN 200610113132 A CN200610113132 A CN 200610113132A CN 100569195 C CN100569195 C CN 100569195C
Authority
CN
China
Prior art keywords
data
dimensional
tooth
occluding relation
model
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.)
Expired - Fee Related
Application number
CNB2006101131320A
Other languages
Chinese (zh)
Other versions
CN1919157A (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.)
Individual
Original Assignee
Individual
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 Individual filed Critical Individual
Priority to CNB2006101131320A priority Critical patent/CN100569195C/en
Publication of CN1919157A publication Critical patent/CN1919157A/en
Application granted granted Critical
Publication of CN100569195C publication Critical patent/CN100569195C/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Abstract

Manufacturing method of fine personalized skull model capable of describing teeth occluding relation of the present invention, the skull model that solves CT scan in the prior art can't accurately be described the defective of teeth occluding relation.This method comprises that employing CT tomoscan obtains patient's cranio-maxillofacial faultage image data; The cranio-maxillofacial computerized three-dimensional reconstruction; The Gypsum Fibrosum delivery of tooth and the three dimensional data collection of plaster mold; The dental part three-dimensional data of the 3rd step gained is carried out data fusion and separated with the 3 D graphic data of second step to the CT image reconstruction; With mandibular bone data and the head data of having peeled off mandibular bone respectively rapid shaping be processed into real model, and carry out basic steps such as organic connection, the present invention can obtain fine personalized skull model capable of describing teeth occluding relation, and clearly teeth occluding relation helps the doctor to understand patient's oral cavity anatomical detail intuitively, user can find treatment needed various points, plane and angle on head mould, so that better grasp the clinical operation operant skill.

Description

Manufacturing method of fine personalized skull model capable of describing teeth occluding relation
Technical field:
The present invention relates to a kind of manufacture method of skull model, particularly a kind of manufacturing method of fine personalized skull model capable of describing teeth occluding relation.
Background technology:
Personalized head mould plays important effect in scheme guidance, the skills training in conceptual design, the art before the art of cranio-maxillofacial plastic surgery, the department of stomatology and other related disciplines.But at present personalized head mould adopts CT data rapid shaping to finish usually, because the spatial resolution of CT data is in the millimeter rank, so meticulous anatomy relationship such as teeth occluding relation can't obtain embodying, and this has just greatly limited the application of head mould in tooth section.And the common refined model that obtains tooth by the mode of Gypsum Fibrosum delivery in the tooth section, but this way only can obtain the last lower jaw data and the tooth data of part, and user can not be determined the anatomy of the cranio-maxillofacial data beyond this part data and the neighbour relation that goes up mandibular bone and cranio-maxillofacial.Therefore existing model production method has limited the use clinically of personalized head mould.
Clearly teeth occluding relation helps the doctor to understand patient's oral cavity anatomical detail intuitively.User can find treatment needed various points, plane and angle on head mould, so that better grasp the clinical operation operant skill.Jaw part device true to nature has sagittal condyle portion consistent with actual patient and condyle road, is fit to the needs of multiple operation on oral cavity planning
Summary of the invention:
The spatial resolution that the present invention solves CT cranio-maxillofacial data in the prior art is low, the anatomical detail in oral cavity can't be described and tooth section model only can obtain the last lower jaw data and the tooth data of part according to its cranio-maxillofacial model of making, can't obtain the technical problem of the neighbour relation of mandibular bone and cranio-maxillofacial, a kind of manufacture method is provided, can produce fine personalized skull model capable of describing teeth occluding relation, can help the doctor to understand patient's oral cavity anatomical detail intuitively.
Technology contents of the present invention is:
A kind of manufacturing method of fine personalized skull model capable of describing teeth occluding relation comprises the steps:
The first step, employing CT tomoscan obtain patient's cranio-maxillofacial faultage image data, become optical file to preserve with standard DICOM format record data.
Second step, cranio-maxillofacial computerized three-dimensional reconstruction: CT cranio-maxillofacial faultage image data is carried out three-dimensional reconstruction, promptly adopt medical image control system software to form the 3 D graphic data of reflection cranio-maxillofacial original shape, the output file form is STL.
The 3rd step, the Gypsum Fibrosum delivery of tooth and the three dimensional data collection of plaster mold: at first the clinician adopts the mode of traditional Gypsum Fibrosum delivery to obtain and the identical plaster model of patient's tooth, adopts the 3-dimensional digital camera to obtain the 3 D graphic data of tooth mould then.
The 4th goes on foot, the dental part three-dimensional data of the 3rd step gained is carried out data fusion and separated with the 3 D graphic data of second step to the CT image reconstruction: at first carry out data fusion, the 3 D graphic data of two kinds of data mode is carried out image registration, make the corresponding three-dimensional point set X of two block graphics data A(x 1, y 1, z 1) and X B(x 2, y 2, z 2) reach the locus and dissect structural in full accord, from above-mentioned two block graphics data, select the anatomical structure point of a plurality of correspondences, and then employing obtains the mapping relations of two kinds of mode based on the method for registering of the singular value decomposition-iterative closest point of contour feature: promptly by following formula definite three-dimensional zoom factor c and translation matrix t:
c = 1 σ x 2 tr ( DS ) - - - ( 1 ) t=μ y-cRμ x (2)
σ in last two formulas xRepresentative
Figure C20061011313200051
The mean square deviation of point coordinates position, tr is a matrix trace, and D is a diagonal matrix, and S is an eigenmatrix, μ yRepresent three-dimensional point set In the average coordinates value of all data, R is a spin matrix, μ xRepresent three-dimensional point set
Figure C20061011313200053
In the average coordinates value of all data,
Then to eclipsed partial data in the two block graphics data, keep the high view data of resolution, the data that resolution is low are used image three-dimensional editor's means and are removed;
Carry out data separating then, promptly the overall data of the data after merging as skull model, again with mandibular joint as separation, mandibular bone is separated, so that the model that post-production goes out just has dynamic interlock analog functuion, and mandibular bone data and the head data of having peeled off mandibular bone are saved as the file of STL form respectively.
The 5th step, mandibular bone data and the head data of having peeled off mandibular bone are input to rapid prototyping system respectively, and then be processed into real model with rapid shaping, adopt the spring of meticulous design that two parts model is connected at the lower jaw joint at last, finally obtain fine personalized skull model capable of describing teeth occluding relation.
Technique effect: the present invention passed through to merge and isolation technics in the 4th step, the three-dimensional data that the three-dimensional data of CT scan and tooth is meticulous merges, and the mandibular bone data are separated, the head data of mandibular bone have been obtained peeling off, and carry out rapid shaping respectively and process, connect then, just can obtain fine personalized skull model capable of describing teeth occluding relation, and clearly teeth occluding relation helps the doctor to understand patient's oral cavity anatomical detail intuitively, user can find the needed various points of treatment on head mould, plane and angle are so that better grasp the clinical operation operant skill.Carry out accurately taking getting owing to jaw portion adopts the dental model three-dimensional digital camera, can obtain relative CT scan jaw portion more true to nature, be fit to the needs of multiple operation on oral cavity planning.
Preferably, the three dimensional data collection of the 3rd step plaster mold of the present invention adopts the white light three-dimensional camera to realize, and adopts software system to obtain 3 d image data.The imaging process of described three-dimensional camera and gatherer process are: laser beam illumination micro optical element produces striated structure light, through the throw light on surface of tested tooth or tooth mould of prism vergence, striated structure light is subjected to the deforming stripe of the modulation formation of tooth three-dimensional appearance, again through image rotation prism and telecentric imaging system imaging on picture receiver, the imaging receptor of telecentric imaging system must be an infrared thermoviewer, to realize the detection to infrared light, the deforming stripe figure that picture receiver receives delivers to computer through video output interface; Software system in the computer is finished following process: the three dimensional display of the analysis automatically of the collection of original image, striped, phase unwrapping, the mapping of degree of depth picture, tooth profile, three-dimensional editor, and the conversion of the multiple general format of tooth three-dimensional model and output.
Adopt above-mentioned homemade three-dimensional camera, the digitized tooth model spatial resolution of collection about 70 nanometers, can clear performance the structure and the occluding relation of tooth; The self-editing software system that employing is complementary, can realize the accurate foundation of 3 D graphic data and the conversion and the output of multiple general format, be convenient to fusion with the CT three-dimensional data, also be convenient to the rapid shaping processing of model, make the model that makes can be in conjunction with the advantage of different imaging modes, the skull model that obtains consistent with the true cranio-maxillofacial of patient and can the fine description teeth occluding relation, have sagittal condyle portion consistent and condyle road with actual patient, the needs that are fit to multiple operation on oral cavity planning, can be used for the phantom of the classical operation of cranio-maxillofacial and tooth section and the formulation of quantification operation plan, can carry out effectively evaluating the postoperative curative effect.The design of movable mandibular joint, except that simple dental articulation motion, the jaw part device can be simulated multiple mandibular movement, as protrusive movement, lateral movement etc., offers the doctor with real oral cavity simulated environment.
In the preferably described second step cranio-maxillofacial three-dimensional reconstruction step, the performing step of medical image control system software comprises: cut apart; Extract area-of-interest; Contour line extracts, follows the tracks of; The fifth wheel profile; Generate tri patch; Photechic effect calculates; Generate 3-D view; Output file, file output format are STL.
Above-mentioned self-editing medical image control system software is based on the software of surface rendering, but in order to obtain the higher 3-D view of rotatable solid sense and degree of accuracy, the applicant adopts the object plotting method based on voxel, also designed and be used to diagnose the three-dimensional reconstruction medical software of working out with operation plan, performing step is: A, cut apart; B, extraction area-of-interest; C, interpolation; D, utilize visual theory that voxel is projected to display plane at last to show; The wherein said multiple pattern of cutting apart that adopts of cutting apart, comprise the pattern of cutting apart based on conditions such as density value and form, proximity relations, described interpolation both can be constructed complete volume data to the shape interpolation of sparse synusia The data based on the chamfering distance, can guarantee the verity of data form again.Adopt the three-dimensional reconstruction software of above-mentioned steps, the three-dimensional reconstruction process control of whole skull is in half a minute, and picture quality is clear, true to nature, has high degree of accuracy.On this basis, this software also can carry out rotation, cutting and complicated measurement of 3-D view.Patient's image of sense of reality can not only be provided for the doctor, more help the formulation of state of an illness diagnosis and treatment and medical operating scheme.
Above-mentioned software based on surface rendering, can from volume data, splice to come the match body surface earlier by geometric units, utilize traditional computer graphics technology that the body surface of rebuilding is drawn then, promptly adopting efficient succinct contour line to link method draws the surface of cranio-maxillofacial, make the threedimensional model surface of generation more smooth than other algorithms, the fairing degree is good, the data file of the STL form of software output can directly be used as the input data of industrial design software and rapidform machine, possesses the favorable compatibility energy, whole three-dimensional reconstruction process need not manual intervention, the surface links the smooth no spring layer sense of nature, and can guarantee very high degree of accuracy.The accuracy of configuration of surface and data has guaranteed the accurate of subsequent design.
Preferably, the CT tomoscan of first step has been write DICOM-SLICE data transaction compressed software at the DICOM image data acquiring, and its specific implementation is: to original DICOM view data be shifted in turn, closely mode is compressed; Next sequential scanning, calculated difference data; The efficient replacement of difference data is represented then, obtains the data of high compression ratio at last.
Description of drawings:
Fig. 1 is the general flow chart of the inventive method;
Fig. 2 is the structural representation that the present invention makes three-dimensional camera by oneself;
Fig. 3 is the FB(flow block) that the present invention is based on the three-dimensional reconstruction software of voxel;
Fig. 4 is the FB(flow block) that the present invention is based on the three-dimensional reconstruction software of surface rendering;
Fig. 5 is the FB(flow block) of compression of images software of the present invention;
Fig. 6 is the cranium jaw face graphics that the present invention's three-dimensional reconstruction software shown in Figure 3 is constructed;
Fig. 7 is the software processing system that the present invention is used for the white light three-dimensional camera.
The specific embodiment:
The present invention is described in further detail below in conjunction with accompanying drawing.
See Fig. 1, described the general flow chart of the inventive method, comprise five basic steps:
The first step, employing CT tomoscan obtain patient's cranio-maxillofacial faultage image data; Second step, cranio-maxillofacial computerized three-dimensional reconstruction; The 3rd step, the Gypsum Fibrosum delivery of tooth and the three dimensional data collection of plaster mold; The 4th goes on foot, the dental part three-dimensional data of the 3rd step gained is carried out data fusion and separated with the 3 D graphic data of second step to the CT image reconstruction; The 5th step, with mandibular bone data and the head data of having peeled off mandibular bone respectively rapid shaping be processed into real model, and carry out organic connection.
Wherein, the first step, employing CT tomoscan obtain patient's cranio-maxillofacial faultage image data, and the data of collection become optical file to preserve with standard DICOM format record.After only collecting view data, just can carry out follow-up medical 3 D shows, the graphical design of the cranio-maxillofacial except that last mandibular bone and manufacturing, in hospital, adopt CT equipment that patient's head is carried out fine scanning, obtain the faultage image of cranio-maxillofacial, faultage image is stored on the storage mediums such as CD with the DICOM form of industry global general-use.Preferably, the applicant is according to the needs of real work, formulated the image file format that the meets dicom standard-SLICE data format of oneself, write DICOM-SLICE data transaction compressed software at the DICOM image data acquiring, this software can be discerned the dateout of all main CT producers at present, and be automatically converted to SLICE packed data form, unified pattern of the input as follow-up work, we can not only finish the CT image data acquiring that all meet dicom standard like this, and can finish the compression of data, help the transmission of online data.
Fig. 5 has shown this DICOM-SLICE data transaction compressed software flow process, for: to original DICOM view data be shifted at first in turn, closely mode is compressed; Next sequential scanning, calculated difference data; Do the efficient replacement of difference data then and represent, obtain the data of high compression ratio at last.The principle that realizes compression is: the DICOM data of CT are 12 bit data, but take two bytes store, will produce 4 spare bits like this, by the mode that is shifted in turn, data are stored in tight mode, no longer include spare bits, adopt the further compressing image data of differential technique then, this is a kind of reversible compression algorithm, the described mode of displacement in turn is: for the gray level image among most of DICOM, gray value between two adjacent pixels differs very little, in computer graphics through the statistics can find all to be distributed in basically-8-+8 between, therefore for every scanning line, if the pixel value of first point is C1, the value of second pixel is C2, if X1=C2-C1, then C2 can be expressed as C1+X1, by that analogy, can obtain X2, X3, X4...... as from the foregoing, these values are distributed in mostly-8-+8 between, definition Y1, Y2, Y3......Y16 represents 0000 respectively, 0001,0010....1111 these 16 binary values, the difference data that obtains is used Y1 respectively, Y2 ... Deng replacement, finally obtain compression ratio than higher view data.
Second step, cranio-maxillofacial computerized three-dimensional reconstruction; Promptly adopt medical image control system software 3 d data field to form the 3 D graphic data (also can be described as 3-D graphic) of reflection cranio-maxillofacial original shape, the output file form is STL.Medical image control system software can adopt conventional Mimics software system, but this system degree of accuracy when being used for that the so complicated system of cranium jaw face set up 3 D graphic data can't satisfy the needs of follow-up hot investment casting fully, therefore, preferably, the medical image control system software that adopts the applicant to work out voluntarily carries out three-dimensional reconstruction.
Above-mentioned self-editing medical image control system software preferably includes two parts, and first is a three-dimensional reconstruction software, called after 3DMSee (3dimension medicine see CT).For patient's cranio-maxillofacial image of sense of reality is provided to the doctor, adopt and realize three-dimensional reconstruction based on the object plotting method of voxel, FB(flow block) is seen Fig. 3, performing step is: at 3 d data field, carry out A, cut apart; B, extraction area-of-interest; C, interpolation; D, utilize visual theory that voxel is projected to display plane at last to show.Wherein A step cuts apart the multiple pattern of cutting apart that adopts, comprise based on density value and form, conditions such as proximity relations cut apart pattern, the interpolation in C step is to the shape interpolation of sparse synusia The data based on the chamfering distance, both can construct complete volume data, can guarantee the verity of data form again, the D step at first composes with visual signature data, comprise to the data of each classification and compose color value and give the data tax portion transparency of each classification, carry out the calculating of photechic effect then, the last image that just can carry out is synthetic, and soon voxel is projected in display plane and carries out stereo display.Adopt the three-dimensional reconstruction software of above-mentioned steps, the three-dimensional reconstruction process control that can make whole skull is in half a minute, and picture quality is clear, true to nature, has high degree of accuracy, sees Fig. 6.This software also can carry out rotation, cutting and complicated measurement of 3-D view.Patient's image of sense of reality not only can be provided for the doctor, more help the formulation of state of an illness diagnosis and treatment and medical operating scheme.
Second portion is the software based on surface rendering, and purpose is in order to be connected mutually with the computer-aided design of repair member, to be used for the manufacturing of cranium jaw face.Because the 3-D view of constructing based on the object plotting method of voxel can not be directly used in the cranio-maxillofacial rapid shaping, setting up the 3 d surface model of describing how much topological structure between each data just can be used for making, therefore, from the rapid shaping angle, the software of first is not the basic step of the inventive method.Because the step that obtains before the volume data is identical with front first, so can from the volume data that first obtains, splice to come the match body surface earlier by geometric units, utilize traditional computer graphics technology that the body surface of rebuilding is drawn then, wherein commonly used have fifth wheel profile method, Marching Cube algorithm and a Dividing Cube algorithm etc., we have adopted efficient succinct contour line connection method that the cranio-maxillofacial surface is drawn, software called after 3DMSR (3 dimension medicine surface rendering).This part software is seen Fig. 4 from the FB(flow block) that scans back gained faultage image data and begin, and flow process is: A, at first the data of three-dimensional space data field are carried out the work of the image segmentation of data reduction and multi-method; B also carries out region of interest and extracts, and above-mentioned steps is identical with A, the B step of first, so can directly utilize these two gained data of first; C, carry out profile then and follow the tracks of to determine the contour line of each region of interest; D, use labelling, methods such as difference shadow etc. to determine the corresponding relation between each layer contour line cleverly thereafter; Thereby E is fifth wheel profile structure three-dimensional object surface data accurately, generate tri patch; F, process photechic effect calculate; G, generate 3-D view at last.Efficient succinct contour line links method the surface of cranio-maxillofacial is drawn, make the threedimensional model surface of generation more smooth than other algorithms, the fairing degree is good, the data file of the STL form of software output can directly be used as the input data of industrial design software and rapidform machine, possesses the favorable compatibility energy, whole three-dimensional reconstruction process need not manual intervention, and the surface links the smooth no spring layer sense of nature, and can guarantee very high degree of accuracy.The accuracy of configuration of surface and data has guaranteed the accurate of subsequent design.
The 3rd step, the Gypsum Fibrosum delivery of tooth and the three dimensional data collection of plaster mold: at first the clinician adopts the mode of traditional Gypsum Fibrosum delivery to obtain and the identical plaster model of patient's tooth, adopt the 3-dimensional digital camera to obtain the 3 D graphic data of tooth mould, the high camera of 3-dimensional digital camera preferable space resolution then.
Common CT data because spatial resolution in the millimeter level, so can't reconstruct the fine structure of tooth clearly, the description of the occluding relation of tooth is not accurate enough, so we wish to finish by advanced person's dentistry digital form the digital collection work of this part.At first the clinician adopts the mode of traditional Gypsum Fibrosum delivery to obtain and the identical plaster model of experimenter's tooth, the tooth mould is all very complicated on geometry and topology, for guaranteeing to keep in the plaster model digitized process precision of himself, preferably, adopt the white light three-dimensional camera of independent research to realize the digitized work of tooth mould.Described white light three-dimensional camera adopts high-precision 3 D digital imaging and the Modeling Technology differentiated more, has merged the theoretical and corresponding sensibility variable multiplephase bitmap method for reconstructing of coded illumination of the variable ability of sensitivity.The principle of this technology is that the striped automatic analysis technology is combined with optics More principle, and mutually synthetic projected fringe figure carries out the 3-dimensional digital picture that digital demodulation is obtained the space object surface by para-position.Infrared laser is selected the laser diode of 850nm for use.The white light three-dimensional camera structure is seen Fig. 2, comprise the telecentric imaging system 2 and the micro optical element 3 that are connected with image rotation prism 1 light respectively, micro optical element 3 can be converted into white light with laser, so that the phase function is directly taken a picture to gather three-dimensional data to tooth, telecentric imaging system 2 links to each other with the input of picture receiver 4, and the outfan 5 output video image data of picture receiver 4 are to the software processing system; 3 of micro optical elements are arranged on the laser optical path of infrared diode 6 generations, infrared diode 6 another termination power interfaces 7.
The imaging process of white light three-dimensional camera and image acquisition process are: laser beam illumination micro optical element produces striated structure light through turning to the surface of the tested tooth of prism 1 turning light, striated structure light is subjected to the deforming stripe of the modulation formation of tooth three-dimensional appearance, be imaged on the picture receiver 4 through image rotation prism 1 and telecentric imaging system 2 again, imaging receptor 4 must be an infrared thermoviewer, to realize the detection to infrared light.The deforming stripe figure that receives delivers to computer through video output interface 4, the flow chart that is used for the software processing system of this white light three-dimensional camera is seen Fig. 7, remove the collection of finishing original image, also finish the three dimensional display of striped analysis automatically, phase unwrapping, the mapping of degree of depth picture, tooth profile, three-dimensional editor, and functions such as the conversion of the multiple general format of tooth three-dimensional model and output.The digitized tooth model spatial resolution of gathering by above-mentioned white light three-dimensional camera about 70 nanometers, structure and occluding relation that can the clear performance tooth.
The 4th goes on foot, the dental part three-dimensional data of the 3rd step gained is carried out data fusion and separated with the 3 D graphic data of second step to the CT image reconstruction:
At first carry out data fusion, the 3 D graphic data of two kinds of data mode is carried out image registration, make the corresponding three-dimensional point set X of two block graphics data A(x 1, y 1, z 1) and X B(x 2, y 2, z 2) reach the locus and dissect structural in full accord.The 3-D view of the 3-D view that the CT data reconstruction the goes out Gypsum Fibrosum tooth mould that scanning obtains with the laser three-D digital camera derives from different imaging modes, in scanning process, there is the misalignment phenomenon on direction and the yardstick, therefore needs us that the 3-D view of two kinds of mode is carried out registration.So-called image registration promptly by seeking certain spatial alternation, makes the corresponding point of two groups of images reach the locus and dissect structural in full accord.The result of registration should make all anatomic points on two width of cloth images, all reaches coupling.Two width of cloth image X A(x 1, y 1, z 1) and X B(x 2, y 2, z 2) registration, seek a kind of mapping relations T:X exactly A→ X B, make X AEach the point at X BOn all have unique point corresponding with it, and these 2 corresponding to same anatomical position, mapping relations T shows as one group of successive spatial alternation.Because all there is error in two groups of images that obtain on direction and yardstick, we take the mapping relations of affine change.We are by empirical mode in the time of registration, select a plurality of corresponding anatomical structure points the 3-D view of the Gypsum Fibrosum tooth mould that obtains from the scanning of CT 3-D view and laser three-D digital camera, and then adopt method for registering based on the singular value decomposition-iterative closest point of contour feature.The advantage that this method combines singular value decomposition optimum solution analysis method and iterative search solves the matching problem of image contour point, adopt singular value decomposition method that the space point range is mated to many groups anatomical features point of selecting, by iterating the optimum registration parameter that obtains profile point.
Singular value decomposition based on contour feature can be described as:
To two three-dimensional point sets
Figure C20061011313200081
With Calculate
μ x = 1 n Σ i = 1 n x i , μ y = 1 n Σ i = 1 n y i , μ wherein xRepresent three-dimensional point set
Figure C20061011313200085
In the average coordinates value of all data, μ yRepresent three-dimensional point set
Figure C20061011313200086
In the average coordinates value of all data,
σ x = 1 n Σ i = 1 n | | x i - μ x | | 2 , σ y = 1 n Σ i = 1 n | | y i - μ y | | 2 σ in the formula x, σ yRepresent the mean square deviation of coordinate position, matrix H = Σ xy = 1 n Σ i = 1 n ( y i - μ y ) ( x i - μ x ) T
H is the covariance matrix of data point coordinate in the formula, and T represents matrix transpose
Then H is made singular value decomposition:
H=UDV T, wherein D is a diagonal matrix, U and V are rule one orthogonal matrix.
Spin matrix R is determined by U, V:
R=VSU T
The mapping relations of two kinds of mode: zoom factor c and translation matrix t are determined by following formula:
c = 1 σ x 2 tr ( DS ) - - - ( 1 ) t=μ y-cRμ x (2)
σ in last two formulas xRepresentative
Figure C20061011313200095
The mean square deviation of point coordinates position, tr is a matrix trace, and D is a diagonal matrix, and S is an eigenmatrix.
Then to eclipsed partial data in the two block graphics data, keep the high view data of resolution, the data that resolution is low are used image three-dimensional editor's means and are removed;
Carry out data separating then, promptly the overall data of the data after merging as skull model, again with mandibular joint as separation, mandibular bone is separated, so that the model that post-production goes out just has dynamic interlock analog functuion, and mandibular bone data and the head data of having peeled off mandibular bone are saved as the file of STL form respectively.
The 5th step, mandibular bone data and the head data of having peeled off mandibular bone are input to rapid prototyping system respectively, and then be processed into real model with rapid shaping, adopt the spring of meticulous design that two parts model is connected at the lower jaw joint at last, finally obtain fine personalized skull model capable of describing teeth occluding relation.
The skull model that obtains by above-mentioned steps combines the advantage of different imaging modes, consistent with the true cranio-maxillofacial of patient and can the fine description teeth occluding relation, can be used for the phantom of the classical operation of cranio-maxillofacial and tooth section and the formulation of quantification operation plan, can carry out effectively evaluating the postoperative curative effect.The design of movable mandibular joint, except that simple dental articulation motion, the jaw part device can be simulated multiple mandibular movement, as protrusive movement, lateral movement etc., offers the doctor with real oral cavity simulated environment.

Claims (6)

1, a kind of manufacture method of describing the fine personalized skull model of teeth occluding relation comprises the steps:
The first step, employing CT tomoscan obtain patient's cranio-maxillofacial faultage image data, become optical file to preserve with standard DICOM format record data;
Second step, cranio-maxillofacial computerized three-dimensional reconstruction: CT cranio-maxillofacial faultage image data is carried out three-dimensional reconstruction, promptly adopt medical image control system software to form the 3 D graphic data of reflection cranio-maxillofacial original shape, the output file form is STL;
The 3rd step, the Gypsum Fibrosum delivery of tooth and the three dimensional data collection of plaster mold: at first the clinician adopts the mode of traditional Gypsum Fibrosum delivery to obtain and the identical plaster model of patient's tooth, adopts the 3-dimensional digital camera to obtain the 3 D graphic data of tooth mould then;
The 4th goes on foot, the 3rd dental part 3 D graphic data that goes on foot gained is carried out data fusion and separates with the 3 D graphic data of second step to the CT image reconstruction: at first carry out data fusion, the 3 D graphic data of two kinds of data mode is carried out image registration, make the corresponding three-dimensional point set X of two block graphics data A(x 1, y 1, z 1) and X B(x 2, y 2, z 2) reach the locus and dissect structural in full accord, from above-mentioned two block graphics data, select the anatomical structure point of a plurality of correspondences, and then employing obtains the mapping relations of two kinds of mode based on the method for registering of the singular value decomposition-iterative closest point of contour feature: promptly by following formula definite three-dimensional zoom factor c and translation matrix t:
c = 1 σ x 2 tr ( DS ) - - - ( 1 ) t=μ y-cRμ x (2)
σ in last two formulas xRepresentative The mean square deviation of point coordinates position, tr is a matrix trace, and D is a diagonal matrix, and S is an eigenmatrix, μ yRepresent three-dimensional point set
Figure C2006101131320002C3
In the average coordinates value of all data, R is a spin matrix, μ xRepresent three-dimensional point set
Figure C2006101131320002C4
In the average coordinates value of all data,
Then to eclipsed partial data in the two block graphics data, keep the high view data of resolution, the data that resolution is low are used image three-dimensional editor's means and are removed;
Carry out data separating then, promptly the overall data of the data after merging as skull model, again with mandibular joint as separation, mandibular bone is separated, so that the model that post-production goes out just has dynamic interlock analog functuion, and mandibular bone data and the head data of having peeled off mandibular bone are saved as the file of STL form respectively;
The 5th step, mandibular bone data and the head data of having peeled off mandibular bone are input to rapid prototyping system respectively, and then be processed into real model with rapid shaping, adopt the spring of meticulous design that two parts model is connected at the lower jaw joint at last, finally obtain describing the fine personalized skull model of teeth occluding relation.
2, the manufacture method of the fine personalized skull model of description teeth occluding relation according to claim 1, the three dimensional data collection that it is characterized in that described the 3rd step plaster mold adopts the white light three-dimensional camera to realize, and adopt software system to obtain 3 D graphic data, the imaging process of described three-dimensional camera and image acquisition process are: laser beam illumination micro optical element produces striated structure light, through the throw light on surface of tested tooth or tooth mould of prism vergence, striated structure light is subjected to the deforming stripe of the modulation formation of tooth three-dimensional appearance, again through image rotation prism and telecentric imaging system imaging on picture receiver, the imaging receptor of telecentric imaging system must be an infrared thermoviewer, to realize the detection to infrared light, the deforming stripe figure that picture receiver receives delivers to computer through video output interface; Software system in the computer is finished following process: the three dimensional display of the analysis automatically of the collection of original image, striped, phase unwrapping, the mapping of degree of depth picture, tooth profile, three-dimensional editor, and the conversion of the multiple general format of tooth three-dimensional model and output.
3, the manufacture method of the fine personalized skull model of description teeth occluding relation according to claim 1 and 2 is characterized in that the performing step of medical image control system software comprises: cut apart in the described second step cranio-maxillofacial three-dimensional reconstruction step; Extract area-of-interest; Contour line extracts, follows the tracks of; The fifth wheel profile; Generate tri patch; Photechic effect calculates; Generate 3-D view; Output file, file output format are STL.
4, the manufacture method of the fine personalized skull model of description teeth occluding relation according to claim 3 is characterized in that above-mentioned medical image control system software, also comprises a system based on object plotting method, and implementation step is: A, cut apart; B, extraction area-of-interest; C, interpolation; D, utilize visual theory that voxel is projected to display plane at last to show.
5, the manufacture method of the fine personalized skull model of description teeth occluding relation according to claim 4, it is characterized in that the described multiple pattern of cutting apart that adopts of cutting apart, comprise the pattern of cutting apart based on density value and form, proximity relations condition, described interpolation is to the shape interpolation of sparse synusia The data based on the chamfering distance.
6, the manufacture method of the fine personalized skull model of description teeth occluding relation according to claim 4, it is characterized in that, the CT tomoscan of first step has been write DICOM-SLICE data transaction compressed software at the DICOM image data acquiring, and its specific implementation is: to original DICOM view data be shifted in turn, closely mode is compressed; Next sequential scanning, calculated difference data; The efficient replacement of difference data is represented then, obtains the data of high compression ratio at last.
CNB2006101131320A 2006-09-15 2006-09-15 Manufacturing method of fine personalized skull model capable of describing teeth occluding relation Expired - Fee Related CN100569195C (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CNB2006101131320A CN100569195C (en) 2006-09-15 2006-09-15 Manufacturing method of fine personalized skull model capable of describing teeth occluding relation

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CNB2006101131320A CN100569195C (en) 2006-09-15 2006-09-15 Manufacturing method of fine personalized skull model capable of describing teeth occluding relation

Publications (2)

Publication Number Publication Date
CN1919157A CN1919157A (en) 2007-02-28
CN100569195C true CN100569195C (en) 2009-12-16

Family

ID=37777098

Family Applications (1)

Application Number Title Priority Date Filing Date
CNB2006101131320A Expired - Fee Related CN100569195C (en) 2006-09-15 2006-09-15 Manufacturing method of fine personalized skull model capable of describing teeth occluding relation

Country Status (1)

Country Link
CN (1) CN100569195C (en)

Families Citing this family (25)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101689309A (en) * 2007-06-29 2010-03-31 3M创新有限公司 The synchronized views of video data and three-dimensional modeling data
EP2229914B1 (en) * 2009-03-20 2018-05-30 Nobel Biocare Services AG System and method for aligning virtual models
CN101604354B (en) * 2009-07-17 2011-05-04 南京航空航天大学 Method for accurately recovering single tooth form on three-dimensional tooth-mandible model
JP5383416B2 (en) * 2009-10-06 2014-01-08 キヤノン株式会社 Image processing apparatus, control method therefor, and program
CN101803958B (en) * 2010-03-12 2013-03-20 北京欧亚瑞康新材料科技有限公司 Integrated fiber post for dental restoration and preparation method thereof
CN101799938B (en) * 2010-04-28 2012-04-25 福州大学 Method for producing real object simulated computer image
CN101964155A (en) * 2010-09-16 2011-02-02 南方医科大学 Method for manufacturing anthropotomy cast specimen model
CA2855645C (en) * 2011-11-15 2020-06-23 Trispera Dental Inc. Method and system for acquiring data from an individual for preparing a 3d model
CN102663818A (en) * 2012-04-16 2012-09-12 上海交通大学医学院附属第九人民医院 Method and device for establishing three-dimensional craniomaxillofacial morphology model
DE102012104282A1 (en) 2012-05-16 2013-11-21 Isra Vision Ag Method and device for inspecting surfaces of a tested object
CN102920534A (en) * 2012-11-07 2013-02-13 上海交通大学医学院附属第九人民医院 Production method of personalized obturator
CN103100792B (en) * 2013-03-12 2015-04-22 合肥知常光电科技有限公司 Laser preprocessing and restoring method and device having on-line detection function and used for optical elements
WO2015123759A1 (en) 2014-02-21 2015-08-27 Trispera Dental Inc. Augmented reality dental design method and system
TWI552729B (en) * 2014-06-13 2016-10-11 Cheng-Xin She A system and method for image correction design of jaw jaw surgery
CN105411716B (en) * 2014-09-22 2017-06-16 北京大学口腔医学院 A kind of edentulous jaw alveolar ridge intercuspal position Direct Determination
EP3353752B1 (en) * 2015-09-26 2021-08-04 Boston Scientific Scimed Inc. Adjustable depth anatomical shell editing
CN105342708B (en) * 2015-12-14 2017-10-24 四川大学 Digitlization bite guide and its method for reconstructing based on CT and CBCT fused datas
CN105528807A (en) * 2016-01-29 2016-04-27 北京正齐口腔医疗技术有限公司 Teeth arrangement design method and device
CN105894445B (en) * 2016-03-31 2019-06-11 北京思创贯宇科技开发有限公司 A kind of coronary artery image processing method and device
CN106934236A (en) * 2017-03-09 2017-07-07 重庆医科大学附属口腔医院 Mandibular motion simulation method based on conical beam CT
RU2652014C1 (en) * 2017-09-20 2018-04-24 Общество с ограниченной ответственностью "Авантис3Д" Method of using a dynamic virtual articulator for simulation modeling of occlusion when designing a dental prosthesis for a patient and a carrier of information
CN108113780B (en) * 2018-02-06 2020-11-10 北京大学口腔医学院 Method for measuring motion trail of mandible condyle for manufacturing mandible condyle prosthesis
CN108992211B (en) * 2018-07-24 2019-12-31 北京大学口腔医学院 Method for manufacturing titanium mesh for alveolar bone defect
CN112022384B (en) * 2020-09-04 2022-12-06 上海交通大学医学院附属第九人民医院 CAD/CAM-based digital resin resolution plate, preparation method and application
CN116671956B (en) * 2023-07-17 2023-10-03 广州医思信息科技有限公司 Oral cavity data acquisition method based on comparison model

Also Published As

Publication number Publication date
CN1919157A (en) 2007-02-28

Similar Documents

Publication Publication Date Title
CN100569195C (en) Manufacturing method of fine personalized skull model capable of describing teeth occluding relation
US20060127854A1 (en) Image based dentition record digitization
KR101744079B1 (en) The face model generation method for the Dental procedure simulation
CN108711187A (en) The method that registration fusion CT and MRI signal establish human body lumbar vertebrae three-dimensional simulation model
CN106570859B (en) Tooth image analysis system
KR20160004862A (en) Patient specific virtual set-up and its simulation and the method manufacturing Orthodontics device or Orthodontics guidance
CN101706971A (en) Automatic division method of dental crowns in dental models
CN110033465A (en) A kind of real-time three-dimensional method for reconstructing applied to binocular endoscope medical image
CN102663818A (en) Method and device for establishing three-dimensional craniomaxillofacial morphology model
CN112270993B (en) Ultrasonic robot online decision-making method and system taking diagnosis result as feedback
CN112151169B (en) Autonomous scanning method and system of humanoid-operation ultrasonic robot
CN106327479A (en) Apparatus and method for identifying blood vessels in angiography-assisted congenital heart disease operation
CN105389444A (en) Gingival margin curve design method for personalized implant tooth
Douros et al. An improved algorithm for reconstruction of the surface of the human body from 3D scanner data using local B-spline patches
Buchaillard et al. 3D statistical models for tooth surface reconstruction
Molnar et al. Idealized landmark-based geometric reconstructionsof poorly preserved fossil material: a case study of an early tetrapod vertebra
CN101533518A (en) Method for re-establishing surface of three dimensional target object by unparallel dislocation image sequence
CN101996415B (en) Three-dimensional modeling method for eyeball
Beylot et al. 3D interactive topological modeling using visible human dataset
CN113057677B (en) Heart image modeling method, system and equipment for fusing ultrasonic image and CT image
CN112132805B (en) Ultrasonic robot state normalization method and system based on human body characteristics
Racasan et al. Integrating the concept of reverse engineering in medical applications
Kim et al. Automatic registration of dental CT and 3D scanned model using deep split jaw and surface curvature
Zhang Virtual reality technology
CN114037744A (en) Oral cavity region-of-interest three-dimensional reconstruction method and device based on convolutional neural network

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20091216

Termination date: 20150915

EXPY Termination of patent right or utility model