Accuracy and reliability of mandibular digital model superimposition based on the morphological characteristics of vessels in extraction adult patients

Abstract

Background

This study aimed to validate the availability of superimposing full-color mandibular digital models (DMs) by the morphological characteristics of vessels in extraction adult patients.

Methods

Twenty-eight adult patients were included, and their DMs were superimposed with pre- and posttreatment cone beam computed tomography (CBCT) and the morphological characteristics of lingual vessels. The measurements of each tooth were compared under the same coordinate system.

Results

The ICC results displayed exceptional agreement in intra- and interrater assessments, with scores exceeding 0.891 in the crown for intrarater agreement and scores surpassing 0.888 in the crown for interrater agreement. Furthermore, no statistically significant differences were found in the 2 superimposition methods (P > 0.05).

Conclusion

The morphological characteristics of vessels under the mucogingival junction in the lingual side of mandible of are stable enough for the superimposition of mandibular DMs in the adult patients undergo orthodontic treatment with premolars extraction.

Background

Dento-maxillofacial function, stability, aesthetics and health can be achieved by orthodontic tooth movement. Therefore, orthodontists must assess three-dimensional tooth position accurately and effectively. Several methods have been reported to evaluate tooth movement: plaster model, radiography such as lateral cephalogram and cone beam computed tomography (CBCT). However, plaster model casts have gradually declined because of their breakable characteristic [1, 2]. The lateral cephalograms reveal the 2-D tooth movement and the airway with defects such as magnification and distortion in the image [3, 4]. As for CBCT images, which are now widespread in the clinic for presenting 3-D images of dental-maxillofacial anatomy without image distortion [5,6,7], it has been reported that CBCT is capable of evaluating tooth movement through voxel-based superimposition [8,9,10,11]. However, for orthodontists, it violates the principle of as low as reasonably achievable (ALARA) to take CBCT repeatedly in a short period only to evaluate tooth movement.

3-D intraoral scanning has become widespread for diagnosis and treatment in dental clinics. This technology allows recording teeth and portions of periodontal tissue without exposing the patient to harmful radiation [12, 13]. The intraoral structures are saved as digital data, which are unbreakable and available for 3-D measurement. Some studies have reported that the palatal rugae and palatal vault could remain stable during orthodontic treatment and be the reference region in maxillary DMs superimposition [14,15,16,17,18]. Regarding mandibular DMs, few anatomical structures can be recorded by a scanner. Therefore, some studies have superimposed mandibular DMs indirectly by CBCT registration or occlusion of the maxilla and mandible [19,20,21]. It is difficult to superimpose the mandibular DMs accurately, even though some researchers considered the mandibular torus could be stable during orthodontic treatment [22]. However, the torus is not prevalent in every patient. Therefore, it is urgent to find a method to superimpose mandibular DMs efficiently and accurately.

Some researchers have confirmed that the branches of vessels in the fundus from different periods are available for image superimposition and further instruct the clinical diagnosis [23, 24]. Thus, we suggested the hypothesis that the morphological characteristics of vessels could also be a potential reference for superimposing mandibular DMs. Our study aims to establish the viability of superimposing full-color mandibular digital models (DMs) by evaluating vessel morphological features in adult patients who have completed orthodontic treatment with premolars extraction. The outcome of this evaluation will be compared to that of CBCT voxel-based DM superimposition, regarded as the benchmark standard.

Methods

All patients were collected from the Department of Orthodontics with the following inclusion criteria: (1) adult patient; (2) CBCT images and DMs data were preserved perfectly; (3) patient finished the orthodontic treatment with 4 premolars extraction. Moreover, the exclusion criteria were as follows: (1) patients with missing teeth (excluding the third molars); (2) patients who had undergone orthognathic surgery; (3) patients with severe periodontal diseases. Twenty-eight adult patients (2 males, 26 females, 22.54 ± 2.19 years) were included, and their CBCT images and DM data at 2 different time points, pretreatment (T1) and posttreatment (T2), were obtained as well. The mean duration of treatment was 26 months.

The department of the medical image acquired the CBCT images with Carestream Health CS 9300 (Carestream Health Inc., Rochester, NY), and the parameters were as follows: Field of View (FOV) 170 mm* 135 mm; 0.3 mm³ pixel size; 90 kVp; 4 mA; and 12–28 seconds of exposure. All patients were required to maintain the intercuspal position with the head fixed.

The department of orthodontics collected the DMs data with TRIOS 3D intraoral scanner (v1.3.4.7, 3-Shape Inc., Copenhagen, Denmark). All scanning operations were conducted under the guidance of the manufacturer. The data we collected were further processed through the procedures shown as follows.

• 1. Marker placement: In 3-Shape Appliance Analyze (v1.9.2.3, 3-Shape Inc., Copenhagen, DK), the branches of vessels in the bilateral lingual side of T1 mandibular DMs could be easily recognized. The markers (sphere in 1 mm diameter) were placed manually along the distribution of vessels (Fig. 1A). The same procedures were repeated in T2 DMs (Fig. 1B). All the markers were merged with DMs (DM-T1 and DM-T2) and exported as standard triangulated language (STL) format (Fig. 1C, D).

• 2. CBCT registration: T1 and T2 CBCT images were imported into Dolphin Imaging (v11.95, Dolphin Imaging and Management Solutions Inc., Chatsworth, CA, USA) and registered together with the mandibular body and part of the ramus as the reference region (Fig. 2), which was considered as voxel-based registration. The registered CBCT-T1 and CBCT-T2 were then exported in DICOM format and converted to STL format (Fig. 2C) by Mimics Research (version 21.0, Materialise N.V., Technologielann, Leuven, BE)

• 3. Settle the position of DM-T1: DM-T1 acquired from STEP 1 was registered to CBCT-T1 acquired from STEP 2 by the morphology of dentition as a reference area in 3-matic Research (v13.0, Materialise N.V., Technologielann, Leuven, BE). Then, the position of DM-T1 was fixed (Fig. 3A).

• 4. CBCT-based DMs superimposition (CBCT group): Because the position of DM-T1 was fixed in CBCT-T1, and the CBCT-T2 had already been registered over CBCT-T1, we only needed to register the DM-T2 over CBCT-T2 as STEP 3 did (Fig. 3B), and then we could obtain the DMs superimposed by CBCT voxel-based registration: DM-CBCT-T2 (Fig. 3C).

• 5. Vessel-based DMs superimposition (Vessel group): DM-T2 was superimposed on DM-T1 acquired from STEP 3 with the vessel markers as the reference area (Fig. 4A). Finally, we acquired the DMs superimposed by the characteristics of vessels: DM-V-T2. Now, we can visually recognize the differences between the 2 superimposition methods (Fig. 4B, C).

• 6. Segmentation: For further study of tooth movement during the treatment, data acquired from the above needed to be further processed. In the 3-Shape Appliance Analyze, pre- and posttreatment mandibular DMs were segmented to the tooth models (Fig. 5A).

• 7. Establish the coordinate system: The tooth models obtained in STEP 6 were aligned with DM-T1, DM-CBCT-T2, and DM-V-T2. These scanning models were then transformed into dentition models (Fig. 5B). Three spheres were positioned at the mesial-buccal cusp of the bilateral mandibular first molars and the midpoint of the bilateral central incisors on DM-T1. These spheres were duplicated on DM-CBCT-T2 and DM-V-T2 to ensure the consistency of the coordinate system (Fig. 6).

• 8. Measurement: The crown points were set at the midpoint of the incisal ridge, cups of canines and the central fossa of molars and premolars (Fig. 7). The 3-D coordinate values of each tooth crown could be determined as (X, Y, Z) in the system (X represented the horizontal tooth movement, Y represented the vertical tooth movement, Z represented the sagittal tooth movement). In addition, 3-D measurements of bilateral lower first molars, canines, incisors movement between T1 and T2 in 2 groups were also calculated. Six landmarks were placed on the middle of incisal edges of the lower middle incisors, canine cusp tips and the central fossa of first molars as black landmarks in DM-T1 models (Fig. 8A, D) and as red landmarks in DM-CBCT-T2 and DM-V-T2 models (Fig. 8B, C, E and F).

Fig. 1

Marker placement. Large boxed areas show high-magnification views. A Pretreatment DMs (DM-T1) with markers. B Posttreatment DMs (DM-T2) with markers. C DM-T1 in standard triangulated language format (Red). D DM-T2 in standard triangulated language format (yellow)

Fig. 2

A Pretreatment CBCT (CBCT-T1). B Posttreatment CBCT (CBCT-T2). C Voxel-based CBCT registration (CBCT-T1: white, CBCT-T2: green, the red box displays the chosen reference area)

Fig. 3

A Settle the position of DM-T1. B Registered the DM-T2 to CBCT-T2. C Acquisition of the CBCT-based DMs superimposition, DM-T1 (red) and DM-CBCT-T2 (blue)

Fig. 4

A Vessel-based DMs superimposition, registration of DM-V-T2 (yellow) to DM-T1 (red). B DM-T1 (red) and DM-CBCT-T2 (blue) and DM-V-T2 (yellow). C Top views of 2 methods of superimposition

Fig. 5

A Segmentation. B Transferring the scanning models to the dentition model

Fig. 6

A Coordinate system of dentition models. The X-axis, a line through the midpoint of bilateral incisors and parallel to the line connecting bilateral first molar mesial-buccal cusps, stands for the horizontal direction. The Y-axis was perpendicular to the X-axis and through the midpoint of bilateral incisors, which stands for the vertical direction. The Z-axis was perpendicular to the X-axis and Y-axis through the midpoint of the bilateral incisors, which represents the sagittal direction. B X-Z horizontal plane, constructed on the mandibular oral plane by 3 spheres. C The X-Y frontal plane was perpendicular to the horizontal plane. D The Y-Z sagittal plane was perpendicular to the horizontal plane and frontal plane

Fig. 7

Measurement. The three-dimensional position of each crown was measured. Crown points were set at the middle fossa of molars or premolars, the cusp of canines, and the midpoint of the incisal edge. A Isometric view. B Top view

Fig. 8

3-D measurement. Six landmarks were placed on the middle of incisal edges of the lower middle incisors, canine cusp tips and the central fossa of first molars as black landmarks in DM-T1 models (A, D) and as red landmarks in DM-CBCT-T2 and DM-V-T2 models (B, E). Measurement of the 3-D distance between pretreatment (DM-T1) and posttreatment (DM-CBCT-T2 and DM-V-T2) were both calculated. Landmarks of DM-T1 might be covered by DM-CBCT-T2 or DM-V-T2 (C, F)

The sample size of this study was calculated by PASS (v15.0, NCSS, LLC. Kaysville, Utah, USA). Wilcoxon’s Sign Rank Test was performed to compare the measurement results obtained from two distinct superimposition methods. Statistical significance was established at a P value less than 0.05, provided the data met the standard distribution criteria and homogeneity of variance. The standard deviation of paired differences was set as 0.45 based on the pre-experiment results, the mean of paired differences was set as 0.3 based on the level of clinical acceptability, the testing power was set as 0.9, α = 0.05, and the nominal size was 26.

SPSS Statistics (v25.0, SPSS Inc., Chicago, IL, USA), was used to process all the statistical analyses. Rater 1 determined the intrarater reliability, repeated all the procedures of superimposition and measurement after 1 week and then tested it using the intraclass correlation coefficient (ICC). This study was based on a single rater (k = 1) using absolute agreement and the 2-way mixed-effects model at a 95% confidence interval. As for the interrater reliability, rater 2 repeated the above steps to register the CBCT and DMs; the measurement results were compared by interclass correlation coefficient, based on a mean rating (k = 2), using absolute agreement and the 2-way random-effects model at a 95% confidence interval.

In addition, two independent T test was used to compare the 3-D measurement of 2 groups (significance established at P = 0.05). The values of root mean square (RMS) between DM-CBCT-T2 and DM-V-T2 were also calculated by part comparison analysis in 3-Matic Research.

Results

No significant differences were found between the measurements of the 2 superimposition methods in the horizontal, vertical, and sagittal directions in the crown (P > 0.05) (Table 1). The median, 25th percentile and 75th percentile of tooth movement in each direction are shown in Table 1. The mean differences between the 2 methods were also calculated (DM-CBCT-T2 subtracted from CBCT-V-T2) and shown in Fig. 9. In the crown, the maximum differences between the 2 methods were + 0.100 mm in the horizontal direction, − 0.179 mm in the vertical direction and + 0.100 mm in the sagittal direction.

Table 1 Wilcoxon’s Sign Rank Test for comparing the tooth movement in three-dimensional direction of each crown point in two superimposition methods (P > 0.05)

Fig. 9

The three-dimensional tooth movement of crown in two superimposition methods was demonstrated, as well as the differences between ‘DM-CBCT-T2’ and ‘DM-V-T2’ (calculated by subtracting ‘DM-CBCT-T2’ from ‘CBCT-V-T2’)

The intraclass and interclass correlation coefficient results showed excellent intrarater reliability (e.g., above 0.891 in the crown, Table 2) and interrater agreement (e.g., above 0.888 in the crown, Table 3).

Table 2 Intrarater agreement of rater 1 for the three-dimensional coordinate values of the crown

Table 3 Interrater agreement of rater 1 and rater 2 for the three-dimensional coordinate values of the crown

No statistical significance was found in the 3-D measurement results between the 2 groups (Table 4, P > 0.05). The RMS values between DM-CBCT-T2 and DM-V-T2 are shown in Table 5.All the RMS values were lower than 0.5 mm. The results of part comparison analysis showed in Fig. 10, the color of DM-V-T2 surface, which ranged from green to red, presented for the differences of two parts.

Table 4 Two independent samples T test showed that there was no statistical significance of 3-D measurements between CBCT group and Vessel group (P > 0.05)

Table 5 The mean, standard deviation and root mean square values between DM-CBCT-T2 and DM-V-T2

Fig. 10

Part comparison analysis. A Top view of DM-CBCT-T2 (blue) and DM-V-T2 (yellow). B, C Top view and frontal view of part comparison analysis of two parts, the color bar (range from green to red) shows the differences

A flow chart depicting the steps for the method of mandibular DM superimposition is provided in Fig. 11.

Fig. 11

Flowchart of the steps for the mandibular digital model (DM) registration method

Discussion

Three-dimensional digital models are widely used for clinical practice because it can record intraoral tissue without radiation efficiently [25]. Meanwhile, DMs are potential tools for sequential orthodontic treatment assessment. However, researchers have not reached a definite conclusion about the evident and stable anatomical structure in mandibular DMs thus far. Thus, it is urgent to find a method to superimpose mandibular DMs efficiently and accurately, especially in extraction cases.

The mucogingival line was considered as a stable region, and its morphological features were used as a reference for mandibular DMs in non-extraction cases [26,27,28,29]. Some researchers have reported that the vessels of the fundus in the different periods were available for superimposing the images [23, 24]. Our study also uncovered an extensive network of vessels on both bilateral lingual sides of the mandibular DMs. Therefore, we proposed the hypothesis of superimposing mandibular DMs by the characteristics of vessels. This study aimed to validate the availability and reliability of superimposing full-color mandibular digital models through the morphological characteristics of vessels.

Presently, few studies have reported methods of superimposing mandibular DMs. A previous study proved that the mandibular torus was a potential reference area of superimposition in the extraction case, and the mandibular alveolar surfaces seem unstable because of active remodeling [22]. Dai attempted to preliminarily superimpose the mandibular DMs through the palatal rugae of the maxillary and the occlusion relationship, then adjust the mandibular rotation and translation by cephalometric film created from CBCT images, and compared the results with the CBCT surface-based mandibular DMs superimposition [19]. No significant differences were found except for the horizontal movement of incisors and canines. However, these two methods could not adjust the rotation and translation in the horizontal direction in the light of cephalometric films. Besides cephalometric radiographs have several disadvantages, such as overlapping structures, magnification and distortion of the image. In our study, the outcome of tooth 3-D movement was compared to that of CBCT voxel-based DM superimposition, which was regarded as the benchmark standard.

We proposed that the method can superimpose mandibular models directly with the lingual vessels as anatomical reference areas, which is simpler and more effective than the other methods. This method not only reduces radioactive damage to patients but also accurately measures the three-dimensional movement of each tooth in different time periods. The past researches have confirmed that the mucogingival junction could stay relatively stable during orthodontic treatment [26, 30, 31]. In some cases, the high frenal attachment of buccal frenum might tract the periodontal tissue and lower the stability of this area [32, 33], which were not suitable for digital model superimposition. Thus, we chose the vessels in the lingual sides of mandible, near the level of mucogingival junction as the potential reference region for superimposition. To further illustrate the accuracy of the method, we compared each pair of DMs from the CBCT group and Vessel group by the function of part comparison analysis in 3-Matic Research, which visually displayed the differences (Fig. 10). The RMS values were all lower than 0.5 mm (Table 5). The minor differences between the Vessel group and the CBCT group were clinically acceptable and showed no statistical significance. However, the active reconstruction of alveolar bone and significant movement of teeth in extraction cases posed challenges for superimposition. To reduce the errors, all reference points avoided settling at the anterior or posterior active remodeling areas. All markers were positioned in the lingual sides, between the premolars and first molar, near the level of mucogingival junction, and dispersed to the greatest extent possible.

There were some inadequacies for improvement of this method. First, scanning the lingual tissue in the mandible may prove difficult due to the unique shape of the tongue and the depth of the oral cavity in individual patients. Second, the vessel markers were placed manually and limited in the region of the premolars to the first molar, so the markers should be dispersed as widely as possible.

Conclusions

The courses of vessels in the bilateral lingual sides of full-color mandibular DMs can remain stable in adult cases during orthodontic extraction treatment and be a potential reference region for mandibular digital model superimposition. This method provides an efficient way to evaluate tooth movement and decreases the risk of radiation exposure to patients.