it is important to estimate the changes in tooth position before and after orthodontic treatment to determine whether the teeth have been moved according to the anchorage requirements or to the planned tooth movement of the treatment objectives .
traditionally , the evaluation of tooth movement has been possible by superimposition of lateral cephalograms .
the method can be technique sensitive and time - consuming while variable head positions in the cephalostat between serial radiographs may lead to errors in the comparison of treatment results .
recent advances in computer technologies such as 3-dimensional ( 3d ) computed tomography ( ct ) , 3d surface scanning technology , computer - aided design , computer - aided manufacturing , and associated software have contributed not only to the precise diagnosis , treatment planning , and simulations , but also to the 3d evaluation of treatment results . although these methods involve additional radiation exposure and time - consuming procedure , more accurate 3d data can be acquired . on the maxillary dental arch , although little is known about the stability of identifiable landmarks on dental casts , palatal rugae has been suggested as relatively stable structures for registration of serial maxillary models.1 the shape of the palatal vault and the medial portions of the palatal rugae are fairly stable throughout the development of the dentition.2 the palatal rugae retain their shape and pattern throughout a person 's lifetime.3 thus , they have been used for identification purposes in forensics.4 almeida et al.5 and bailey et al.6 have studied the potential use of the palatal rugae for the superimposition of serial models and both studies concluded that specific parts of the palatal rugae ( e.g. , medial ) may be sufficiently stable to serve as anatomical references for superimposing serial maxillary models with headgear or premolar extraction treatment .
hoggan and sadowsky7 suggested that certain landmarks on the palatal rugae can be used as reliably as cephalometric superimpositions to assess anteroposterior molar movement .
technical developments technology now permit the 3d superimposition of the maxillary dental arch using 3d model scanning technology .
ashmore et al.8 reported that the method developed for the superimposition of digital configurations of serial dental models allowed accurate measurements of translational movement of the maxillary first - molar in 3 dimensions for both headgear and untreated groups .
miller et al.9 reported that the digital superimposition of the maxillary arch was reproducible and that the error associated with using palatal rugae as reference landmarks may be similar to or less than that of current 2d cephalometric analyses.10 cha
et al.11 reported that use of a 3d digital orthodontic model superimposition technique on the maxillary dental arch was clinically as reliable as cephalometric superimposition for assessing orthodontic tooth movement .
however , none of these studies were performed on the mandibular dental arch since no well - known stable reference points or areas for superimposition have been established on the mandibular dental arch . no accurate and reproducible method has been introduced to date to evaluate the degree of 3d orthodontic tooth movement on the mandibular dental arch despite the available technology . the purpose of this study was to introduce a method for mandibular dental arch superimposition in a non - gorwing patient by the combination of 3d cone beam ct data and 3d digital model data .
the method was designed to use the basal bone area as the reference for superimposition while soft tissue or alveolar bone areas were avoided since they can be altered throughout orthodontic treatment .
two different registration methods for superimposing the mandible were designed to evaluate their reproducibility as well as practical applicability .
one adult patient with no growth potency who had visited the department of orthodontics , kyung hee dental hospital and undergone orthodontic treatment , was examined .
cone beam computed tomography ( cbct ) scan and alginate impression were taken on the same day prior to treatment , and in accordance with the patient 's request , 1 day after debonding .
the superimposion system consisted of 6 main procedures : ( 1 ) acquisition of cbct data , ( 2 ) acquisition of 3d digital model data , ( 3 ) integration of cbct data and digital model data , ( 4 ) superimposition of pre - treatment and post - treatment integrated cbct image , ( 5 ) extraction of the mandibular dentition portion , and ( 6 ) evaluation of tooth movement on the mandibular arch .
craniofacial skeleton 3d cbct image data was acquired with an i - cat imaging device ( imaging sciences international , hatfield , pa , usa ) .
the patient underwent cbct scanning in facial mode in 14-bit gray scale with 1.0 mm slice thickness , 1.0 mm voxel size , and image acquisition by single 360 degree rotation with a 40 s total scan time .
the reconstructed digital data was directly transferred from the cbct scanner to a personal computer and stored as digital imaging and communications in medicine ( dicom ) files .
both pre - treatment and post - treatment dicom data were transferred to rapidform 2006 ( inus technology inc . ,
the mandible was segmented by rapidform 2006 from the craniofacial skeleton image ( figure 2 ) .
orapix scanner ; a 3d surface scanning system with a slit laser beam and non - contact measurement method ( kod-300 , orapix co. ltd . ,
seoul , korea ) was used to obtain 3d data from the dental cast model .
both pre - treatment and post - treatment orthodontic 3d digital models are shown in figure 3 . to create a preliminary fusion model ,
the dentition reconstructed by cbct data was superimposed on the 3d digital model by a regional registration method .
this function of rapidform 2006 , designated as 3d surface - to - surface matching ( best - fit - method ) , employs a least - mean - squared algorithm.12,13 according to the coordination between the tooth image from cbct and the digitally scanned dental cast , the position , size , and posture of each tooth on the 3d digital model was automatically fitted to one of the subject 's images of dentition in the cbct data ( figure 4 ) . finally , both pre- and post - treatment digital models were integrated into each mandible model , and the bases of the digital model superimposed on each image were removed .
thus , we produced individual integrated 3d cbct images ( integrated 3d cbct image [ ici]t1 , icit2 ) of the mandible , before and after the treatment , respectively , with an anatomical structure that included a digital model of higher resolution than the 3d cbct image ( figure 5 ) .
two different registration methods of superimposing icit1 and icit2 were used to evaluate their practical applicability and to estimate their reproducibility ( figure 6 ) . the first method designated as ' surface superimposition ' was based on the surface of the basal bone structure of the mandible by surface - to - surface matching ( best - fit method ) .
the registration area was comprised of the basal bone structures of the mandible , which had not been altered by orthodontic treatment .
the fully automated registration was computed by rapidform 2006 , with the limited manual process of selection and de - selection of the registration area .
the inferior portion of mandibular body and the posterior portion of ramus were selected as the registration area . the procedure for surface superimposition
the centers of both sides of mental foramen and the center of mandibular lingual foramen ( a consistent arterial foramen in the middle of the mandible)14,15 were selected as the landmarks .
the center of each foramen was selected with a grid on the monitor by picking 4 points of the foramen - the most superior , inferior and both lateral points - with the mandible rotated to give a perpendicular view to the foramen ( figures 8 and 9 ) .
the procedure for plane superimposition was shown in figure 10 . before the 2 superimposition procedures
the reference points were the center of right and left mental foramen , the center of lingual foramen , the facial axis points16 of right and left lower central incisors , canines , and 2nd premolars ( figure 11 ) . specially in this patient the mesiobuccal tip of the lower right 1st molar was also marked since it was found to be an ankylosed tooth .
rapidform 2006 has a coordinate system , of which its origin is defined with the 0-point value by the software itself .
the landmarks were established and selected on the 3d digital model image by 2 investigators using the rapidform program .
the points assignment was repeated 5 times each by each investigator , and the mean coordinate values of each landmark were used to define the final coordinate values .
the final coordinate values ( x , y , and z ) of each landmark were transferred to microsoft office excel 2007 ( microsoft , redmond , wa , usa ) .
the dentition was extracted from the superimposed image ( figure 12 ) , and the base of the model previously removed was added again on the dentition .
the software easily performed the re - addition procedure at the same position since the bases were previously constructed on the digital model .
finally , the superimposed 3d mandibular dental arch was ready to be used for evaluation of the orthodontic tooth movement that occurred throughout treatment .
the initially established landmarks consisted of 10 points on both pre- and post - treatment images . after applying the 2 methods of superimposition
, the coordinate values ( x , y , z ) of the same landmarks were transferred to microsoft office excel 2007 .
each method of superimposition was performed 30 times to estimate reproducibility ( steps 4 to 6 ) and all the coordinate values of post - treatment landmarks after each superimposition were transferred to microsoft office excel 2007 .
first , the mean values and variances of each x- , y- , and z - coordinate values of the post - treatment landmarks were calculated with both the surface superimposition and plane superimposition methods .
the equality of variance test was carried out for each x- , y- , and z - coordinate values as a scalar with matlab 7.5 ( minitab inc . , state college , pa , usa ) . the null hypothesis and alternative hypothesis
were expressed as h0 : p < s and ha : p > s ( h0 : null hypothesis , ha : alternative hypothesis , p : variance of plane superimposition , s : variance of surface superimposition )
. the difference in variance between the 2 superimposition methods was 0 at a significance level of 0.05 .
then the mean distance ( mean value of the distance from the mean coordinate value of the 30 trials to the coordinate values of each trial ) was calculated for both methods with matlab .
student 's t - test was performed to determine whether a significant difference existed between the mean distances of the 2 superimposition methods for each landmark using sas 9.1 ( sas institute inc . ,
the null hypothesis was rejected at a significant level of 0.05 . to compare the stability of the reproduced results from both surface and plane superimposition ,
the following 2 methods were used , since the data used for the analysis involved 3d coordinates , and no appropriate method for the precise comparison of 2 different superimposition methods has been established . at first , 3d coordinates were divided into x , y , z axes and the variance for results from each superimposition method was calculated in each axis .
subsequently , f - test for equal variances was used to compare variance values , and the difference in distance between values from each trial and the mean values of the 2 superimposition methods were compared .
craniofacial skeleton 3d cbct image data was acquired with an i - cat imaging device ( imaging sciences international , hatfield , pa , usa ) .
the patient underwent cbct scanning in facial mode in 14-bit gray scale with 1.0 mm slice thickness , 1.0 mm voxel size , and image acquisition by single 360 degree rotation with a 40 s total scan time .
the reconstructed digital data was directly transferred from the cbct scanner to a personal computer and stored as digital imaging and communications in medicine ( dicom ) files . both pre - treatment and post - treatment dicom data
the mandible was segmented by rapidform 2006 from the craniofacial skeleton image ( figure 2 ) .
orapix scanner ; a 3d surface scanning system with a slit laser beam and non - contact measurement method ( kod-300 , orapix co. ltd . ,
seoul , korea ) was used to obtain 3d data from the dental cast model .
both pre - treatment and post - treatment orthodontic 3d digital models are shown in figure 3 . to create a preliminary fusion model ,
the dentition reconstructed by cbct data was superimposed on the 3d digital model by a regional registration method .
this function of rapidform 2006 , designated as 3d surface - to - surface matching ( best - fit - method ) , employs a least - mean - squared algorithm.12,13 according to the coordination between the tooth image from cbct and the digitally scanned dental cast , the position , size , and posture of each tooth on the 3d digital model was automatically fitted to one of the subject 's images of dentition in the cbct data ( figure 4 ) .
finally , both pre- and post - treatment digital models were integrated into each mandible model , and the bases of the digital model superimposed on each image were removed .
thus , we produced individual integrated 3d cbct images ( integrated 3d cbct image [ ici]t1 , icit2 ) of the mandible , before and after the treatment , respectively , with an anatomical structure that included a digital model of higher resolution than the 3d cbct image ( figure 5 ) .
two different registration methods of superimposing icit1 and icit2 were used to evaluate their practical applicability and to estimate their reproducibility ( figure 6 ) . the first method designated as ' surface superimposition '
was based on the surface of the basal bone structure of the mandible by surface - to - surface matching ( best - fit method ) .
the registration area was comprised of the basal bone structures of the mandible , which had not been altered by orthodontic treatment .
the fully automated registration was computed by rapidform 2006 , with the limited manual process of selection and de - selection of the registration area .
the inferior portion of mandibular body and the posterior portion of ramus were selected as the registration area . the procedure for surface superimposition
the centers of both sides of mental foramen and the center of mandibular lingual foramen ( a consistent arterial foramen in the middle of the mandible)14,15 were selected as the landmarks .
the center of each foramen was selected with a grid on the monitor by picking 4 points of the foramen - the most superior , inferior and both lateral points - with the mandible rotated to give a perpendicular view to the foramen ( figures 8 and 9 ) .
the procedure for plane superimposition was shown in figure 10 . before the 2 superimposition procedures ,
the reference points were the center of right and left mental foramen , the center of lingual foramen , the facial axis points16 of right and left lower central incisors , canines , and 2nd premolars ( figure 11 ) . specially in this patient the mesiobuccal tip of the lower right 1st molar was also marked since it was found to be an ankylosed tooth .
rapidform 2006 has a coordinate system , of which its origin is defined with the 0-point value by the software itself .
the landmarks were established and selected on the 3d digital model image by 2 investigators using the rapidform program .
the points assignment was repeated 5 times each by each investigator , and the mean coordinate values of each landmark were used to define the final coordinate values .
the final coordinate values ( x , y , and z ) of each landmark were transferred to microsoft office excel 2007 ( microsoft , redmond , wa , usa ) .
the dentition was extracted from the superimposed image ( figure 12 ) , and the base of the model previously removed was added again on the dentition .
the software easily performed the re - addition procedure at the same position since the bases were previously constructed on the digital model .
finally , the superimposed 3d mandibular dental arch was ready to be used for evaluation of the orthodontic tooth movement that occurred throughout treatment .
craniofacial skeleton 3d cbct image data was acquired with an i - cat imaging device ( imaging sciences international , hatfield , pa , usa ) .
the patient underwent cbct scanning in facial mode in 14-bit gray scale with 1.0 mm slice thickness , 1.0 mm voxel size , and image acquisition by single 360 degree rotation with a 40 s total scan time .
the reconstructed digital data was directly transferred from the cbct scanner to a personal computer and stored as digital imaging and communications in medicine ( dicom ) files . both pre - treatment and post - treatment dicom data
the mandible was segmented by rapidform 2006 from the craniofacial skeleton image ( figure 2 ) .
orapix scanner ; a 3d surface scanning system with a slit laser beam and non - contact measurement method ( kod-300 , orapix co. ltd . ,
seoul , korea ) was used to obtain 3d data from the dental cast model .
both pre - treatment and post - treatment orthodontic 3d digital models are shown in figure 3 .
to create a preliminary fusion model , the dentition reconstructed by cbct data was superimposed on the 3d digital model by a regional registration method .
this function of rapidform 2006 , designated as 3d surface - to - surface matching ( best - fit - method ) , employs a least - mean - squared algorithm.12,13 according to the coordination between the tooth image from cbct and the digitally scanned dental cast , the position , size , and posture of each tooth on the 3d digital model was automatically fitted to one of the subject 's images of dentition in the cbct data ( figure 4 ) .
finally , both pre- and post - treatment digital models were integrated into each mandible model , and the bases of the digital model superimposed on each image were removed .
thus , we produced individual integrated 3d cbct images ( integrated 3d cbct image [ ici]t1 , icit2 ) of the mandible , before and after the treatment , respectively , with an anatomical structure that included a digital model of higher resolution than the 3d cbct image ( figure 5 ) .
two different registration methods of superimposing icit1 and icit2 were used to evaluate their practical applicability and to estimate their reproducibility ( figure 6 ) . the first method designated as ' surface superimposition '
was based on the surface of the basal bone structure of the mandible by surface - to - surface matching ( best - fit method ) .
the registration area was comprised of the basal bone structures of the mandible , which had not been altered by orthodontic treatment .
the fully automated registration was computed by rapidform 2006 , with the limited manual process of selection and de - selection of the registration area .
the inferior portion of mandibular body and the posterior portion of ramus were selected as the registration area . the procedure for surface superimposition
the centers of both sides of mental foramen and the center of mandibular lingual foramen ( a consistent arterial foramen in the middle of the mandible)14,15 were selected as the landmarks .
the center of each foramen was selected with a grid on the monitor by picking 4 points of the foramen - the most superior , inferior and both lateral points - with the mandible rotated to give a perpendicular view to the foramen ( figures 8 and 9 ) .
the procedure for plane superimposition was shown in figure 10 . before the 2 superimposition procedures ,
the reference points were the center of right and left mental foramen , the center of lingual foramen , the facial axis points16 of right and left lower central incisors , canines , and 2nd premolars ( figure 11 ) .
specially in this patient the mesiobuccal tip of the lower right 1st molar was also marked since it was found to be an ankylosed tooth .
rapidform 2006 has a coordinate system , of which its origin is defined with the 0-point value by the software itself .
the landmarks were established and selected on the 3d digital model image by 2 investigators using the rapidform program .
the points assignment was repeated 5 times each by each investigator , and the mean coordinate values of each landmark were used to define the final coordinate values .
the final coordinate values ( x , y , and z ) of each landmark were transferred to microsoft office excel 2007 ( microsoft , redmond , wa , usa ) . the dentition was extracted from the superimposed image ( figure 12 ) , and the base of the model previously removed was added again on the dentition .
the software easily performed the re - addition procedure at the same position since the bases were previously constructed on the digital model .
finally , the superimposed 3d mandibular dental arch was ready to be used for evaluation of the orthodontic tooth movement that occurred throughout treatment .
the first method designated as ' surface superimposition ' was based on the surface of the basal bone structure of the mandible by surface - to - surface matching ( best - fit method ) .
the registration area was comprised of the basal bone structures of the mandible , which had not been altered by orthodontic treatment .
the fully automated registration was computed by rapidform 2006 , with the limited manual process of selection and de - selection of the registration area .
the inferior portion of mandibular body and the posterior portion of ramus were selected as the registration area . the procedure for surface superimposition
the second method , designated as plane superimposition was based on anatomical structure . the centers of both sides of mental foramen and the center of mandibular lingual foramen ( a consistent arterial foramen in the middle of the mandible)14,15 were selected as the landmarks .
the center of each foramen was selected with a grid on the monitor by picking 4 points of the foramen - the most superior , inferior and both lateral points - with the mandible rotated to give a perpendicular view to the foramen ( figures 8 and 9 ) .
the reference points were the center of right and left mental foramen , the center of lingual foramen , the facial axis points16 of right and left lower central incisors , canines , and 2nd premolars ( figure 11 ) .
specially in this patient the mesiobuccal tip of the lower right 1st molar was also marked since it was found to be an ankylosed tooth .
rapidform 2006 has a coordinate system , of which its origin is defined with the 0-point value by the software itself .
the landmarks were established and selected on the 3d digital model image by 2 investigators using the rapidform program .
the points assignment was repeated 5 times each by each investigator , and the mean coordinate values of each landmark were used to define the final coordinate values .
the final coordinate values ( x , y , and z ) of each landmark were transferred to microsoft office excel 2007 ( microsoft , redmond , wa , usa ) .
the dentition was extracted from the superimposed image ( figure 12 ) , and the base of the model previously removed was added again on the dentition .
the software easily performed the re - addition procedure at the same position since the bases were previously constructed on the digital model .
finally , the superimposed 3d mandibular dental arch was ready to be used for evaluation of the orthodontic tooth movement that occurred throughout treatment .
the initially established landmarks consisted of 10 points on both pre- and post - treatment images . after applying the 2 methods of superimposition
, the coordinate values ( x , y , z ) of the same landmarks were transferred to microsoft office excel 2007 .
each method of superimposition was performed 30 times to estimate reproducibility ( steps 4 to 6 ) and all the coordinate values of post - treatment landmarks after each superimposition were transferred to microsoft office excel 2007 .
first , the mean values and variances of each x- , y- , and z - coordinate values of the post - treatment landmarks were calculated with both the surface superimposition and plane superimposition methods .
the equality of variance test was carried out for each x- , y- , and z - coordinate values as a scalar with matlab 7.5 ( minitab inc . , state college , pa , usa ) .
the null hypothesis and alternative hypothesis were expressed as h0 : p < s and ha : p > s ( h0 : null hypothesis , ha : alternative hypothesis , p : variance of plane superimposition , s : variance of surface superimposition ) .
the difference in variance between the 2 superimposition methods was 0 at a significance level of 0.05 .
then the mean distance ( mean value of the distance from the mean coordinate value of the 30 trials to the coordinate values of each trial ) was calculated for both methods with matlab .
student 's t - test was performed to determine whether a significant difference existed between the mean distances of the 2 superimposition methods for each landmark using sas 9.1 ( sas institute inc
the null hypothesis was rejected at a significant level of 0.05 . to compare the stability of the reproduced results from both surface and plane superimposition ,
the following 2 methods were used , since the data used for the analysis involved 3d coordinates , and no appropriate method for the precise comparison of 2 different superimposition methods has been established . at first , 3d coordinates were divided into x , y , z axes and the variance for results from each superimposition method was calculated in each axis .
subsequently , f - test for equal variances was used to compare variance values , and the difference in distance between values from each trial and the mean values of the 2 superimposition methods were compared .
final superimposed 3d digital models were produced ( figure 13 ) to permit the evaluation of tooth movement visually as well as numerically .
twenty - eight coordinate values demonstrated statistically significant differences ( p < 0.05 ) between the variance of each method , indicating that the plane superimposition method generated more variability in coordinate values after each trial .
figure 14 shows a scatter plot of coordinate values on the x - y , y - z , and z - x axis of the lower left incisor , which demonstrated a large variation in the coordinate values ( x : f = 59.3 , p < 0.001 , y : f = 11.61 , p < 0.001 , z : f = 491.16 , p < 0.001 ) after each trial between the 2 methods .
student 's t - test revealed a statistically significant difference ( p < 0.001 ) in the mean distance from the mean coordinate value of each landmark to the corresponding coordinate value of each trial , indicating that the surface superimposition method had relatively more consistent coordinate values closer to the mean on all landmarks ( table 2 ) .
many 3d instruments enable orthodontists to acquire 3d images of the facial structures , including both hard and soft tissues .
not only the 3d visualization of the structures but also the measurements and comparison of different evaluation times have been made possible . with advances in technology , 3d superimposition with 3d digital models on the maxillary arch has been developed with reliable reference areas for superimposition . in recent years
ct has been adopted for 3d virtual surgical planning and simulation of post - operative outcomes in orthognathic surgery with various registration methods for superimposition.17 - 20 this study aimed to introduce a reproducible superimposition technique for the mandibular dental arch by the generating integrated 3d cbct images based on cbct scans and 3d surface scanning data with improved imaging quality of the dentition .
this study was limited to a non - growing patient who did not undergo an orthognathic surgical procedure to exclude the effects of mandibular growth or surgical change from the evaluation and comparison of the reproducibility of the 2 designed methods .
cbct was selected as the device for ct scanning for several reasons ; it provides orthodontists with a new modality for research and permits vertical scanning of patients in a normal seated position .
ferrario et al.21 reported that 3d measurement values were not sensitive to head posture . in this study ,
the patient underwent cbct scanning with a slice thickness of 1.0 mm , which improved the visual quality of 3d reconstruction .
further improvement is possible with a smaller slice thickness , but this change would also result in increased image size , requiring greater computational power and longer user interaction time .
recently , a computerized composite skull model was introduced to orthognathic surgery,22 which eliminated the need for plaster casts to create digital models since a laser surface scanner subsequently scans the impressions of the dentition .
thus , one potential source of errors that reduces the accuracy of the complete system appears to have been removed . however , specialized radiolucent dental impressions trays with fiducial markers are needed to allow for the registration of the dental arches assessed by ct and laser scans .
this technique introduces new sources of error , because several additional steps accompanying the data acquisition and registration procedures are required .
therefore , in this study we chose a simpler approach using a plaster cast scanned with an orapix scanner for the registration of the cbct scans and the 3d dental surface images . extremely accurate modeling of the facial structures
however , it is not possible to ascertain whether dental morphology can be acquired as accurately with this technique due to artifacts from metallic restorations .
moreover , intercuspation during occlusion leads to overlapping images , resulting in inaccurate morphology of the individual tooth .
the resultant effect on the ct images appears as pronounced dark and bright streaks , non - linear edge gradients , and sampling errors arising from the surface restoration.23 swennen et al.24 discussed the drawbacks of ct imaging ; it does not provide detailed surface dental morphology or precise intercuspation data due to its limited resolution , and the interocclusal relationship is often obscured by radio - opaque dental restorations or orthodontic brackets . combining 3d cbct images and the 3d virtual model data can compensate for these issues .
gateno et al.22 introduced a technique to integrate digital models into a virtual 3d ct bone model of the skull .
measurements acquired from using 3 mm titanium spheres as markers demonstrated high accuracy ( 0.1 - 0.5 mm ) in the study of a single dry skull .
nkenke et al.25 suggested a considerable progress in the fusion of images from different imaging modalities using software - based approaches in which 2 modalities are combined and mounted in a single coordinate system .
combining ct with optical 3d imaging seems to be a reasonable approach to correct metallic artifacts , allowing precise assessment of the tooth surface and accurate simulation of the post - operative dental occlusion . in a single patient ,
accuracy levels of 0.66 and 0.56 mm for the mandible and maxilla , respectively , with a modified double ct scan procedure resulted in accurate registration , with registration error of 0.1355  0.0323 mm , and analysis of variance of 0.0564 mm.24 two known methods for the registration of ct data and scanned optical dental images have been setablished : point - based registration of external fiducial markers22,26 and surface - based registration of anatomical structures.25 we used surface - based registration referred to as regional registration with surface - to - surface matching ( best - fit method ) , which included a least square technique .
this approach is advantageous in that it can be applied with no additional fiducial markers , which in turn require more additional time - consuming procedures that may also introduce new sources of error .
although their study did not assess a combination technique for ct and dental images , cha et al.11 evaluated the accuracy of measurements by calculating the superimposition discrepancies between pre- and post - treatment 3d models in the maxillary arch .
surface - to - surface matching ( best - fit method ) resulted in a mean error of 0.0399 mm , standard error of 0.0289 mm , and standard deviation of 0.1583 mm , which the authors concluded as a satisfactory accuracy level .
subsol et al.27 proposed methods that have guided the progress of our study , in which semi - landmarks on the surface were used to incorporate information about vectors in the vicinity of the landmark . the landmarks were established and selected on the 3d digital model images 5 times each by 2 investigators .
the mean coordinate values of each landmark were used to define the final coordinate values .
the determined landmarks were then fixed to the image so as not to be changed during the superimposition procedure , and the coordinate values were then determined after every superimposition procedure for a total of 30 repetitions with each method . in this study
, we introduced 2 methods for superimposing ici , designated as surface superimposition and plane superimposition .
plane superimposition was designed to include more information concerning anatomical structures in the superimposition procedure .
the landmarks used - mental foramen and lingual foramen - are known as relatively stable structures on the mandible . in a study of 314 dried mandibles , mcdonnell et al.28 reported that the lingual foramen present in 311 specimens ( 99.04% ) , and that the wall of the canal , not the genial tubercles , produced the radio - opacity peripheral to the foramen seen on radiographs . although the authors remarked that the foramen was not seen on many radiographs of the lower incisor region , a change in orientation of the x - ray beam could account for this discrepancy . in this study
, the registration plane was oriented to the 3 points of these structures that were considered stable .
a lower variance of samples from one method compared to the other , indicates that the method is relatively more efficient in terms of reproducibility . however , in the present study each sample was a series of 3d - vectors , resulting in a covariance matrix for each method instead of a scaler variance .
therefore , we decomposed these series of vectors into 3 series of scalars using x , y , and z coordinates and then used the f - test to determine whether 2 samples from these 2 methods had the same variance .
the surface method demonstrated a lower sample variance than the plane method within a 5% significance level at 28 coordinates .
we attempted to express the magnitude of the 3d vector while excluding the directional component .
the mean distance from the mean coordinate value of each landmark was measured for both methods .
student 's t - test revealed a statistically significant difference in the mean distance between the 2 methods and also indicated that the surface superimposition method had a relatively more consistent distribution closer to the mean coordinate values over 30 trials .
therefore , the surface method was considered to provide more reproducible results after repeated trials .
this reproducibility of the surface superimposition method was helped by our negligible observer variability , which allowed image analysis procedures largely independent of observer errors . by contrast
, the relative inconsistency of plane superimposition is thought to result from the procedure , which requires manual selection of the anatomical structures that could generate inter- or intraobserver variability .
cevidanes et al.29 introduced a fully automated superimposition method using voxel - wise rigid registration of the cranial base , as the cranial base structures are not altered by surgery .
this technique represents an advance from the processe described by kawamata et al.,30 who used an observer - dependent method to superimpose and rotate the post - surgery ct until anatomical landmarks overlapped these same structures in the pre - surgery semi - transparent model .
our results also show improved reproducibility with the automated superimposition technique ( surface superimposition ) , which minimizing the possibility of observer errors .
because it is almost fully computerized except for the selection process of the registration area , it avoids human error , thereby simplifying the process and producing integrated images with a high degree of reproducibility in positioning .
therefore , surface superimposition method may be more precise and convenient method for mandibular surface model superimposition .
the data used for the analysis involved 3d coordinates , and no appropriate method for the precise comparison has been established . as alternatives , the following 2 methods were used .
first , 3d coordinates were divided into x , y , and z axes and the variances for each result from both superimposition methods was calculated on each axis . the f - test for equal variances
was then used to compare the values of variance , and a conclusion was drawn by comparing this calculated f value to the critical value of f distribution ( alpha = 0.05 ) since the f value is known to demonstrate f distribution under conditions of the null hypothesis .
second , the difference in distance between the results and mean points that were generated by the 2 superimposition methods were compared . for this aim ,
mean x , y , and z values were deduced from 3d coordinates , and the distances of each result from the corresponding mean values were calculated
. then the mean distance distribution and mean values that were generated from the 2 superimposition methods were analyzed .
the null hypothesis was that the mean distance calculated by plane superimposition was less than or equal to that calculated by surface superimposition .
the alternative hypothesis was that the mean distance calculated by plane superimposition was greater than that calculated by surface superimposition .
standard deviation of the test statistic t is known to demonstrate t - distribution under the conditions of the null hypothesis .
therefore , a conclusion was drawn by comparing the t value calculated as described above and the ctitical value of t - distribution ( alpha = 0.05 ) .
future studies should include more patient samples for the superimposition procedures for reproducibility and also a comparison with the measurement of tooth movement detected on the conventional cephalographs .
this study was limited to a non - growing patient , and more studies may be needed to reflect the growth changes of the mandible during superimposition .
the true methodological error of the superimposition system was not estimated in this study , and such an estimation would be challenging because the mandible is an independent structure with a joint on the craniofacial skeleton .
however , the error of registration and the reproducibility of the 2 methods were determined . because significant changes in mandibular position can occur during orthodontic treatment , such as autorotation of the mandible , condylar displacement , and centric occluson / centric relation discrepancy , the registration reference in this study was designed within the mandible itself , not on craniofacial structures .
additionally , the registration reference was based on the basal bone of the mandible , which , in contrast to soft tissues or alveolar bone , had not been altered by orthodontic tooth movement . therefore , this method is assumed to meet the submillimetric accuracy requirements of orthodontic tooth movement in non - growing patients .
it is expected to gain acceptance as the new method helps patients to better understand the 3d changes in the mandibular dental arch throughout the treatment process .
the visualization of 3d superimposition and measurements can help orthodontists to better evaluate treatment outcomes . here
, we introduced 2 methods for 3d mandibular dental arch superimposition ; surface superimposition and plane superimposition .
the surface superimposition method demonstrated more reliable and reproducible results based on its minimization of the possibility of observer errors .
these findings suggest that 3d mandibular dental arch superimposition can be achieved using the surface superimposition method .
 objectivethe purpose of this study was to develop superimposition method on the lower arch using 3-dimensional ( 3d ) cone beam computed tomography ( cbct ) images and orthodontic 3d digital modeling.methodsintegrated 3d cbct images were acquired by substituting the dental portion of 3d cbct images with precise dental images of an orthodontic 3d digital model . 
 images were acquired before and after treatment . 
 for the superimposition , 2 superimposition methods were designed . 
 surface superimposition was based on the basal bone structure of the mandible by surface - to - surface matching ( best - fit method ) . 
 plane superimposition was based on anatomical structures ( mental and lingual foramen ) . for the evaluation , 
 10 landmarks including teeth and anatomic structures were assigned , and 30 times of superimpositions and measurements were performed to determine the more reproducible and reliable method.resultsall landmarks demonstrated that the surface superimposition method produced relatively more consistent coordinate values . 
 the mean distances of measured landmarks values from the means were statistically significantly lower with the surface superimpositions method.conclusionsbetween the 2 superimposition methods designed for the evaluation of 3d changes in the lower arch , surface superimposition was the simpler , more reproducible , reliable method . 
INTRODUCTION
MATERIALS AND METHODS
Superimposition procedure
Step 1. Acquisition of 3D CBCT skeletal data and reconstruction of 3D mandibular bone images
Step 2. Acquisition of 3D digital model data by surface scanning
Step 3. Initial integration of the mandible and mandibular dentition model
Step 4. Superimposition of pre- and post-treatment ICI
Surface superimposition
Plane superimposition
Establishment of landmarks on dentition
Step 5 and 6. Extraction of mandibular dentition
Statistical analysis
RESULTS
DISCUSSION
CONCLUSION