|Home | About | Journals | Submit | Contact Us | Français|
The study aim is to investigate the influence of scan field, mouth opening, voxel size, and segmentation threshold selections on the quality of the three-dimensional (3D) surface models of the dental arches from cone beam computed tomography (CBCT). 3D models of 25 patients scanned with one image intensifier CBCT system (NewTom 3G, QR SLR, Verona, Italy) using three field sizes in open- and closed-mouth positions were created at different voxel size resolutions. Two observers assessed the quality of the models independently on a five-point scale using specified criteria. The results indicate that large-field selection reduced the visibility of the teeth and the interproximal space. Also, large voxel size reduced the visibility of the occlusal surfaces and bone in the anterior region in both maxilla and mandible. Segmentation threshold was more variable in the maxilla than in the mandible. Closed-mouth scan complicated separating the jaws and reduced teeth surfaces visibility. The preliminary results from this image-intensifier system indicate that the use of medium or small scan fields in an open-mouth position with a small voxel is recommended to optimize quality of the 3D surface model reconstructions of the dental arches from CBCT. More research is needed to validate the results with other flat-panel detector-based CBCT systems.
Computerized three-dimensional (3D) models of the craniofacial region are a recent trend in dentistry. 3D models of the hard-tissue dental arches including the teeth and the jawbones have been used for various clinical applications including pre-operative treatment planning of dental implants and craniofacial surgical procedures, fabrication of dental and craniofacial prosthesis, analysis of arch dimensions, virtual treatment simulation in orthodontics, and postoperative treatment outcome assessment [1–6]. 3D models of the dental arches are typically obtained by digitizing the dental cast with high-resolution surface laser scanning [7–10]. Those models can be used to evaluate the dentition and inter-occlusal space; however, they do not provide information regarding the amount of alveolar bone available and the relation of the tooth root to its socket. In addition, the accuracy of the digitized model is limited by the accuracy of the dental impression and cast, which could be variable and inconsistent overtime depending on several factors [11–13].
3D models of the dental arches can also be obtained from computed tomography (CT) scans. Conventional CT technology provides accurate reconstructions of the alveolar bone but not the teeth due to limited spatial resolution and strong streak artifacts caused by metallic dental restorations and orthodontic brackets when present. Therefore, “composite models” were developed where the high-quality 3D reconstructions of the alveolar bone and the roots of the teeth obtained from CT were combined with the high-quality 3D reconstructions of the teeth crowns obtained from laser surface scanning of the dental cast. And, while this technique produced satisfactory results, it was deemed unsuitable for the clinical practice due to time and cost constraints associated with complicated setup and the extensive user experience required for correct implementation [2, 14–20].
Cone beam CT (CBCT) specifically developed for the maxillofacial region provides comparable images to conventional CT at reduced radiation dose and cost [21–26]. CBCT reconstructions at present have smaller voxels in comparison with conventional CT , which could be advantageous in obtaining more accurate 3D surface models of the teeth crowns. However, CBCT subjective image quality is still inferior in comparison with conventional CT. Several artifacts including beam hardening and inhomogeneity and truncation influence image contrast and bone border definition in CBCT [28–31]. The Feldkamp filtered back-projection reconstruction algorithm employed by most CBCT manufacturers is exactly identical to the Radon inverse-transform algorithm used in conventional CT in the mid-plane . As such, it does not inherently reduce the effect of streak artifacts caused by metallic dental filling and orthodontic brackets. Moreover, since CBCT is inferior to CT in terms of contrast resolution, the effect of metal streak artifacts may be even more pronounced . Several reconstruction algorithms and techniques were proposed to reduce those artifacts in CBCT images, but they remain to date computationally expensive and not yet widely adopted by CBCT manufacturers .
There are currently many CBCT systems which are commercially available. Those systems are categorized according to detector design technology into: (1) image intensifier tube/charged coupled devices (IIT/CCD) combination or (2) flat panel detector (FPD) . It has been reported that IIT/CCD suffers from more artifacts and increased noise levels compared to FPD systems . Also, in practice several patient scanning and data reconstruction parameters have influence on CBCT subjective image quality [35, 36]. All those factors combined could thus directly influence the quality of the 3D surface models reconstructions of the dental arches from CBCT. It is necessary to assess the quality of those models and the influence those parameters might have in order to optimize patient scanning, data reconstruction, and 3D surface model creation protocols. This is important in order to assess whether it is possible to obtain accurate 3D surface models of the dental arches from CBCT or not. Few studies based on dry skull samples assessed the accuracy of 3D models reconstructions from CBCT [37, 38]. However, 3D surface models created from phantoms, dry skulls, or even formalin-fixed cadavers are not realistic and do not represent actual patient. The quality of the segmented model could deviate largely from what is observed clinically. The objective of this study is to assess the quality of 3D surface models of the dental arches with respect to the influence of different scanning and reconstruction parameters in one IIT-based CBCT system (NewTom 3G) in an in vivo sample of patients.
Twenty-five datasets scanned with the IIT/CCD system Newtom 3G CBCT (QR SLR, Verona, Italy) were selected from a larger database of CBCT patients. The datasets were divided into three groups: Group (A) consisted of ten patients scanned with the 12-in. (large) detector scan field of view (FoV). Group (B) consisted of ten patients scanned with the 9-in. (medium) scan FoV. Group (C) consisted of five patients scanned with the 6-in. (small) scan FoV. Groups (A) and (B) were scanned in a closed-mouth position with the teeth in maximum intercuspid relation, while group (C) was scanned in an open-mouth position. Additional selection criteria were that (1) the upper and the lower jaws are both visible in the scan, (2) no more than four teeth are missing for both jaws excluding the third molars, and that (3) there are no orthodontic brackets or large metal restorations. Informed consent was obtained from the patients to use their data for research purposes. Datasets were exported according to the manufacturer's default settings for each scan field in DICOM 3 file format at the isotropic voxel size of (0.3, 0.25, and 0.2 mm3) for the 12-, 9-, and 6-in. scan FoVs, respectively. The datasets were 12 bits in depth, and the gray values (212=4,096) range was (−1,000 to 3,095). The datasets were imported into 3D analysis software (Amira v4.2, Visage Imaging, Carlsbad, CA, USA) for making the 3D models.
The histogram of a CBCT dataset is composed of a wide range of gray-scale values that represent the X-ray attenuation profiles of the different soft and hard tissues. It is generally more difficult to specify the correct threshold value to separate the bone from soft tissue and background in CBCT than in conventional CT due to inherent inconsistencies in the histogram. A single threshold value was specified to segment the bone and teeth from the background and soft tissue for each 3D model. The surface models were created using the marching cube algorithm . To determine an optimal threshold value, the histogram of each model was approximated as a mixture of Gaussians by using the stochastic expectation maximization (SEM) algorithm. The threshold value was then specified in the region of the intersection of the two Gaussians representing soft tissue and bone. To eliminate the stochastic part of the SEM, the average of three threshold values was taken [29, 40]. Threshold values were determined using this method separately for the maxilla and for the mandible when the jaws were separated (Fig. 1).
For groups (A), (B), and (C), a smaller region of interest (ROI) limited only to the dental arches was selected, and the rest was digitally removed (Fig. 2). In each group, the upper arch was separated from the lower arch by using a cubic ROI selection. For this experiment, the isotropic voxel size of 0.3 mm3 was chosen for each model. Fifty 3D surface models were created in total (25 upper and 25 lower arches) for the three groups (A, B, C).
Group (C) was selected to assess the influence of the voxel size on 3D model quality. This was done in order to assess the influence of voxel size selection on the definition of the occlusal surfaces since groups (A) and (B) were scanned in closed-mouth position that the occlusal surfaces were difficult to visualize. To assess the influence of each level of the “voxel size” factor independently, voxels were categorized into two components: (a) pixel resolution in the (x,y) scan plane (PRxy) and (b) axial slice thickness (AST) in the (z) plane. First, PRxy was fixed at 0.3 mm, and AST was increased from 0.3 to 1.2 mm using 0.3 mm steps. Then, AST was fixed at 0.3 mm and PRxy was increased from 0.3 to 1.2 mm using 0.3 mm steps. A 3D bi-cubic resampling filter was used to manipulate the voxel size to change the pixel resolution and the axial slice thickness. The bi-cubic resampling filter uses information from 16 adjacent points to each voxel for resampling the data. This has the added advantage of smoother resampling and less interpolation artifact compared to bilinear or the nearest neighbor filtering. In total, eight voxel combinations were obtained per model (AST 0.3, 0.6, 0.9, 1.2 and PRxy 0.3, 0.6, 0.9, 1.2). A total of 80 models were created (40 upper and 40 lower) dental arches.
A single investigator who did not participate in the observations created the models. All 3D models were coded, and two observers (one maxillofacial radiologist and one maxillofacial surgeon) assessed the quality of all the models independently. The observers were blind to the scan field and voxel size selections. The models were presented to the observers in a random order. All models were viewed on a 19-in. flat panel screen (1,280×1,024, Philips Brilliance 200WP, Brussels, Belgium). The observers were allowed to rotate, scale (zoom), and translate (move) the models to improve visibility of certain structures, but adjusting the threshold value was not permitted. The observers were asked to assess the subjective visibility of the following structures for each model in the maxilla and the mandible on a five-point scale (1 = very poor, 2 = poor, 3 = satisfactory, 4 = good, 5 = very good):
The observation data were entered and analyzed using SPSS software (SPSS v.15, SPSS, Chicago, IL, USA). The observers' ratings were entered into a multivariate analysis of variance (MANOVA) with the Wilk's lambda test to assess the influence of scan field, voxel size, and segmentation threshold on visibility of the different structures of the 3D models. The three scan field groups and the different voxel combinations were compared to each other using Helmert planned contrasts. Inter-observer agreement was determined using Cohen kappa, and observers' interaction with scan FoV selection was calculated in the MANOVA analysis. Alpha level was set to 0.05
Mean threshold values were (660±44.72, 850±51.63, and 784±60.82) for groups (A), (B), and (C), respectively. There was a statistically significant difference in the threshold value among the three groups (p=0.001). Mean threshold and standard deviation value were (847±33.21) and (662±121.2) for the mandible and maxilla, respectively.
There was an overall statistically significant difference among the three groups (A, B, and C) in 3D surface model quality (p=0.0001). Inter-observer agreement was moderate (kappa=0.53). There was no statistically significant interaction between the observers and the scan field groups (p=0.12). Specifically, the visibility of the external surfaces of the teeth in the maxilla and the interproximal space between the teeth in the anterior region in the maxilla and the mandible was better in group (C) than in group (A) (p=0.0001) with no significant differences from group (B) (p=0.48; Fig. 3). Image noise was significantly less in group (A) than in groups (B) and (C) (p=0.0001, 0.001), respectively.
There was an overall statistically significant difference in 3D model quality in group (C) due to manipulating pixel resolution in the scan plane (PRxy; p=0.0001) but not due to axial slice thickness manipulation (p=0.87). There was no significant interaction between the observers and PRxy (p=0.08). Specifically, the visibility of the external and occlusal surfaces of the teeth, anterior alveolar bone, and interproximal space anteriorly in the maxilla and the mandible was significantly improved by the small pixel resolution PRxy=0.3 mm in comparison with other selections (p=0.001; Fig. 4). There were no significant differences in subjective image quality between (PRxy=0.6 mm) and (PRxy=0.9 mm) selections (p=0.17) but a significant difference between (PRxy=0.6 mm) and (PRxy=1.2 mm) selections (p=0.03).
The results indicate that scan field selection has significant influence on the quality of the 3D models. The small scan field selection provided best visibility for the different structures of the 3D models. Interestingly, there was no significant difference between the small and the medium scan fields in the quality of the 3D models. This means that the medium field, which covers the temporomandibular joints bilaterally and several other important anatomical regions including the maxillary sinus, can be used without loss in quality. The large field reduced the visibility of the teeth surfaces and the interproximal space between the teeth. This is significant because the standard scanning protocol of orthodontic patients with CBCT is with the large scan field. The results here show, however, that there is a significant loss in quality of the 3D models of the dental arches when the large field is used. Conversely, there were more image artifacts associated with the small scan field than with the large scan field. This could be due to cone truncation for small field selections.
An important factor when comparing different scan FoVs is the local tomography problem. When a small FoV is selected, the anatomical structures outside the region of interest are imaged as well since the line integrals sampled by the detector pass through them. Due to the fact that those tissues are only sampled for a small angular range, the resulting image reconstruction is inconsistent [41–43]. Smaller CBCT scan FoVs suffer from greater variability in the density gray values as well (i.e., more inconsistencies) compared to larger scan FoVs [35, 44].
The different scan fields have different voxel sizes, and various CBCT scanners also differ from each other in voxel size selections for each scan field. The importance of voxel size stems from a practical observation that very small voxels (e.g., 0.2 or 0.3 mm3 isotropic) result in an extremely large surface mesh model, which is difficult to process to create an accurate 3D surface model for pre-operative treatment planning and simulation. The results here show that large pixel resolution of 0.6, 0.9, or 1.2 mm in the x,y plane significantly reduce the visibility of the occlusal surfaces of the teeth, interproximal space between the teeth, and alveolar bone. Interestingly, increasing the axial slice thickness in the z plane up to 0.9 mm did not significantly reduce the quality of the 3D models. That means that anisotropic voxels with small pixel area in the x,y plane and larger slice thickness can be used instead of the “standard” small isotropic voxels. This has the added advantage of significantly reducing the model size to facilitate processing while maintaining image quality and also to reduce image noise. The choice of larger voxels reduces image noise, which is caused by photon count statistics, by averaging the gray-level values across slices, which in turn, reduces the overall noise level in the image .
Segmentation threshold value was automatically determined and observer independent . There was more variation in the threshold value due to scan field selection than due to the different samples within the same field. Variation in the threshold value was less in the mandible than in the maxilla. This can be explained that cortical bone in the mandible is thick enough to keep the attenuation profile uniform across the entire bone surface, while in the maxilla, the varied thin cortical bone especially in the palate and tuberosity regions creates significant “bone dehiscence and fenestration” artifacts effect in the 3D model. Due to limited contrast, the roots could not be separated from their sockets using binary thresholding alone without creating artifacts. No attempt was made to systematically evaluate the visibility of the roots since they were very difficult to visualize.
A closed-mouth position made it extremely difficult to separate the jaws using a cubic ROI selection to visualize the occlusal surfaces of the teeth and to assess the inter-occlusal relationship. A closed-mouth position necessitates manual segmentation to separate the teeth from each other, which was a tedious, time-consuming, and user-dependent procedure that resulted in unsatisfactory results. An open-mouth position facilitates separating the teeth to visualize the interproximal space. The mandible can be then virtually rotated into maximum occlusion using software tools.
This study was conducted to evaluate the influence of several scanning and reconstruction parameters on the quality of 3D model reconstructions of the dental arches from one IIT/CCD CBCT scanner (NewTom 3G). The required surface description includes clear display of the teeth surfaces and the occlusal contacts, exact reconstruction of the alveolar bone, plus separation of the teeth from each other and each tooth root from its socket. The study results show that it is still difficult to meet all those demands since separating the root from its socket is still difficult to achieve. Also, to produce a “usable” 3D surface model with sufficient quality, several scanning and reconstruction parameters need to be optimized first, and an optimized segmentation approach must be adhered to.
This study was limited in that no interaction between the scan field selection and the voxel size was assumed. For practical time and resources constraints, it was not possible to assess the influence of voxel size selection on all fields. Instead, the small field was selected a priori based on the hypothesis that it provides better subjective image quality. That hypothesis was supported by the study results. Also, the visibility of the occlusal surfaces with the large and medium fields were not assessed since it was not possible to obtain scan material with the 9- and the 12-in. scan fields with open-mouth position. The study is limited only to one IIT/CCD CBCT system (NewTom 3G). Other FPD based systems are different with respect to image quality and scan and reconstruction parameters. Therefore, caution should be exercised when extrapolating the results to other CBCT systems. Lastly, the influence of streak artifacts was not investigated to limit the number of assessed factors on subjective image quality. However, it is expected that streak artifacts will have a significant influence on the quality of the 3D models.
In conclusion, within the limitations of the current study, several scanning and reconstruction parameters need to be optimized first before good-quality models can be created. The use of the small or medium field in an open-mouth scan position coupled with an anisotropic voxel of 0.3–0.4 mm pixel resolution in the scan plane and an axial slice thickness of 0.6–0.7 mm is recommended when creating 3D models of the dental arches from NewTom 3G CBCT (Fig. 5).
Conflict of interest The authors declare that they have no conflict of interest.
Open Access This article is distributed under the terms of the Creative Commons Attribution Noncommercial License which permits any noncommercial use, distribution, and reproduction in any medium, provided the original author(s) and source are credited.