Skip to main content

Comparison of in vivo 3D cone-beam computed tomography tooth volume measurement protocols

Abstract

Background

The objective of this study is to analyze a set of previously developed and proposed image segmentation protocols for precision in both intra- and inter-rater reliability for in vivo tooth volume measurements using cone-beam computed tomography (CBCT) images.

Methods

Six 3D volume segmentation procedures were proposed and tested for intra- and inter-rater reliability to quantify maxillary first molar volumes. Ten randomly selected maxillary first molars were measured in vivo in random order three times with 10 days separation between measurements. Intra- and inter-rater agreement for all segmentation procedures was attained using intra-class correlation coefficient (ICC).

Results

The highest precision was for automated thresholding with manual refinements.

Conclusions

A tooth volume measurement protocol for CBCT images employing automated segmentation with manual human refinement on a 2D slice-by-slice basis in all three planes of space possessed excellent intra- and inter-rater reliability. Three-dimensional volume measurements of the entire tooth structure are more precise than 3D volume measurements of only the dental roots apical to the cemento-enamel junction (CEJ).

Background

Historically, the in vivo detection of changes to dental root morphology such as those associated with external root resorption (ERR) during the course of orthodontic treatment or related to trauma has been mainly through use of two-dimensional (2D) radiographs, most notably periapical radiographs [1]-[3]. Although histological studies have identified a relatively high incidence of apical ERR, 2D radiographic studies have been less definitive [4],[5] and have in general proven inaccurate for the reliable detection of small ERR defects [6]. In fact, 2D periapical radiographs do not reveal external root resorption to an appreciable extent, except for frank apical root resorption, which appears to be in their realm of identification [7]. In addition, there are geometric limitations associated with 2D imaging of a three-dimensional (3D) phenomenon; therefore, the quantitative value of 2D radiographs to measure ERR is questionable [8]-[10]. When considering panoramic films, the distortion in both tooth position and angulation combined with varying magnification, distortion, superimposition, and imaging artifacts in different parts of the image [11],[12] leads to similar limitations in the use of panoramic films to assess ERR [13]-[15]. Therefore, although 2D radiography may be a good screening tool, its use in the quantification of ERR remains controversial [9].

Advancements into 3D imaging techniques have facilitated volumetric imaging capabilities not previously available on an in vivo basis; however, accurate dental volume measurement procedures are required in order to fully utilize the technology. The resulting use of 3D imaging has enabled the quantification and measurement of ERR to be completed with a high level of diagnostic accuracy and repeatability when compared to periapical radiographs [9],[16]-[18]. The strength of cone-beam computed tomography (CBCT) for accurate dental volume measurements in vivo has been shown not to be statistically significantly different as in vitro measurements in one study [19] and even when comparing its accuracy to in vitro micro-CT imaging methods [20]; however, there may exist machine-specific variations. The feasibility of in vivo dental volume measurements using CBCT imaging has similarly been reported by Liu et al.; however, their use of post-processing surface smoothing has been shown to decrease 3D volume measurements [21]. Conversely, increasing the voxel size has been shown in vitro to actually increase volume measurements [22]. It is intuitively apparent that the accuracy of the 3D segmentation procedure is related to the voxel size during acquisition [23] with 0.25 mm voxel size an appropriate compromise between diagnostic accuracy and patient radiation dose [22]. An additional factor is the development of a clearly defined measurement protocol, which appears lacking in the literature as the study currently employing CBCT as a means of determining root volume loss with maxillary expansion lacks a clearly defined measurement protocol involving incorrectly utilized Hounsfield units with the teeth of interest ‘segmented cautiously’ [24]. There exists a potential limiting factor inherent in the use of CBCT scans to measure accurate volumetric information as the time period required to capture the radiograph as patient movement during scans can reduce the accuracy of measurements [15],[25].

The validation of CBCT as a tool for measurement of both root lengths and volumes has been focused on in numerous studies [15]-[17],[19],[21],[26]-[33]. The investigation of in vivo volumetric determination utilizing CBCT images has been shown to yield slight differences from actual physical volumes within −4% to +7% [21]. However, there is a lack of a clearly defined gold-standard 3D segmentation protocol in the literature.

Given the inconsistencies of the techniques reported in the literature, the development of an appropriate CBCT measurement protocol possessing accuracy and precision in both intra- and inter-rater reliability for in vivo dental volume measurements is desired. Due to the relative infancy of the area of 3D dental volume segmentation, with a lack of a gold-standard technique, the need to employ and evaluate segmentation techniques to identify which measurement protocol is most superior is a necessity. Through the establishment of precise and accurate dental volume measurement protocols, clinicians can more confidently employ the available tools to monitor such phenomena as ERR during the orthodontic process and ERR related to dental trauma.

The objective of this study is therefore to analyze a set of developed and proposed image segmentation protocols for precision in both intra- and inter-rater reliability for in vivo tooth volume measurements using CBCT images.

Methods

Cone-beam CT images

The radiographic data set used for the analysis of dental volume was previously acquired as part of a randomized clinical trial at the University of Alberta, Edmonton, Alberta, Canada. Subjects were recruited during an 18-month period. Inclusion criteria for selection included transverse maxillary deficient adolescents with no previous orthodontic treatment. The age range of patients selected for this study ranged from 11 to 17 years old. Subjects were not excluded based on the presence or absence of coronal restorations. Informed consent from the patients' parents and ethical approval from the Ethics Committee at University of Alberta was obtained.

All CBCT images were taken with the NewTom 3G (QR SRL, Verona, Italy) device at 110 kV, 6.19 mAs, and 8-mm aluminum filtration with the patient in maximum intercuspation following common CBCT imaging protocols. Images were converted to DICOM format by using the NewTom software to a voxel size of 0.25 mm. The DICOM-formatted images were volume rendered with Avizo 3D analysis software (Visualization Sciences Group, Berlin, Germany) [34]. Patient images were acquired at two time points during the trial: T1 (before treatment) and T2 (after treatment, approximately 12 months).

Tooth volume measurement protocols

Three protocols for dental tooth volume determination were investigated using Avizo 3D analysis software:

  1. 1.

    Manual human segmentation on a repeated 2D basis;

  2. 2.

    Automated segmentation without human refinement; and

  3. 3.

    Automated segmentation with manual human refinement on a repeated 2D basis.

In addition, two methods for tooth volume selection were simultaneously investigated. These involved the entire tooth structure including the crown and only the dental root structure apical to the cemento-enamel junction (CEJ). All three protocols and two methods were applied to determine the technique producing greatest intra- and inter-rater reliability. The dental pulp chamber and canals were included in the volume measurements. The investigator was blinded to whether they were T1 (before treatment) or T2 (after treatment) radiographs. In all, a total of six different approaches (combination of three protocols and two tooth volumes) to tooth volume segmentation were investigated. Ten randomly selected maxillary first molars (selected from both T1 and T2 patient images) were measured in vivo in random order three times with 10 days separation between measurements.

The threshold value for image segmentation was set for each tooth. This same threshold value was used in all protocols to assess the particular tooth of interest to limit variability between methods. The first protocol did not require a threshold value to be explicitly set as the protocol was strictly manual human tracing of the image on a 2D slice-by-slice basis.

No image orientation adjustments were completed prior to testing of the protocols. The sagittal plane was utilized for initial evaluations for each technique, as it appeared most useful in the visualization and evaluation of the tooth structure of the crown and root simultaneously.

1. Manual human segmentation on a repeated 2D basis

The first protocol involved manual image segmentation procedures on a 2D slice-by-slice basis through the use of Avizo's ‘lasso’ tool, which allows one to define an area freehand by generating a closed contour curve in 2D. The delineation of tooth structure from the surrounding alveolar and cortical bone was first determined on a slice-by-slice basis in the YZ (sagittal) plane (Figure 1) based upon visual inspection only. Refinements in the XY (axial) plane (Figure 2) were then manually completed for the observation of the tooth anatomy from a different perspective. An axial view enabled root structure and interproximal contact point refinements. Finally, additional refinements in the XZ (coronal) plane (Figure 3) were again manually completed. A coronal view enabled refinements to root structure that was in close proximity to the buccal and palatal cortical plates. The 3D resultant tooth was evaluated for approximately normal maxillary first molar dental anatomy to limit gross misidentification of dental structures (Figure 4). Once segmentation was completed, the software automatically computed the tooth's radiographic volume. No smoothing functions were applied to the 3D tooth structure to prevent smoothing of minor root defects/imperfections or possible resorption lacunae. Both the complete tooth volume (Figure 5A) and the dental root volume, defined as the anatomical root apical to the CEJ, (Figure 5B) were measured.

2. Automated segmentation without human refinement

The second protocol involved the use of the ‘magic wand’ tool in Avizo 3D imaging software as a ‘region-growing’ tool. The ‘magic wand’ tool allows one to perform the so-called ‘region-growing’ in either 2D or 3D. Selecting an individual ‘seed voxel’ of a tooth root or crown selects the largest connected area (either 2D or 3D) that contains the voxel itself and all voxels with gray values contained within a user-specified range. The range can be chosen to represent absolute gray values or gray values relative to that of the seed voxel. For the purposes of our investigation, absolute gray values were chosen to limit variability in selection of the seed voxel gray value. Segmentation was performed using strictly an automated approach after minor operator input to the selection of the seed voxel in the enamel of the tooth without focused manual refinements in an attempt to test an efficient measurement procedure. The user input to select the seed voxel proved to be a necessity given the software. The rest of the procedure required no operator input for the actual segmentation procedure. A visually defined optimal threshold value was set for each tooth in the YZ (sagittal) plane (Figure 1). The threshold level was set to most clearly show the tooth anatomy with minimal interference from the surrounding bone and adjacent structures. The 3D resultant tooth was evaluated for approximately normal maxillary first molar dental anatomy to limit gross misidentification of dental structures. Once segmentation was completed, the software automatically computed the tooth's radiographic volume. As in protocol 1, no smoothing functions were applied and both the complete tooth volume and dental root volume were measured.

3. Automated segmentation with manual human refinement on a repeated 2D basis

The third protocol also involved the use of the ‘magic wand’ tool in Avizo 3D imaging software as a ‘region-growing’ tool, similar to that utilized in the second protocol; however, in this case, segmentation was performed using a mixture of an automated approach with manual localized visual refinements to the tooth structure. For the purposes of our investigation, absolute gray values were chosen to limit variability in selection of the seed voxel gray value. The same absolute gray value range was selected as in the second protocol to limit variability between methods for each tooth. Segmentation was performed using a mixture of an automated approach with manual localized visual refinements to the tooth structure. A visually defined optimal threshold value was set for each tooth in the YZ (sagittal) plane (Figure 1). The threshold level was set to most clearly show the tooth anatomy with minimal interference from the surrounding bone and adjacent structures. Manual refinements were processed on a slice-by-slice basis to enhance accuracy by correcting for over- and under-contoured voxels in the tooth volume. Initial refinements occurred in the YZ (sagittal) plane. Secondary refinements were performed in the XY (axial) plane (Figure 2) to refine root structure and interproximal dental contact points. Tertiary refinements were performed in the XZ (coronal) plane (Figure 3) to verify tooth anatomy and focus on the delineation of dental root structure from the buccal and palatal cortical plates. The 3D resultant tooth was evaluated for approximately normal maxillary first molar dental anatomy to limit gross misidentification of dental structures. Once segmentation was completed, the software automatically computed the tooth's radiographic volume. As in protocols 1 and 2, no smoothing functions were applied and both the complete tooth volume and dental root volume were measured.

Figure 1
figure 1

YZ (sagittal) plane.

Figure 2
figure 2

XY (axial) plane.

Figure 3
figure 3

XZ (coronal) plane.

Figure 4
figure 4

3D sample view of maxillary first molar volume without smoothing.

Figure 5
figure 5

Three-dimensional tooth volume. (A) Three-dimensional complete tooth volume. (B) Three-dimensional tooth volume apical to the cemento-enamel junction.

Statistical analysis

The volume data was manually entered into Microsoft Excel 2011 for MAC (Microsoft, Redmond, WA, USA). SPSS for MAC (version 21, IBM, Armonk, New York, USA) was used to run all statistical tests. For all tests, statistical significance was set at an α value of 0.05.

Intra-class correlation coefficient (ICC) was used to measure the agreement between the measurements for the continuous dependent variable (dental tooth volume) taken on the three separate days. A single measure with consistency under two-way mixed model was chosen, thus removing the rater's variation, and the subjects/teeth were chosen randomly with the rater fixed. The technique that produced the highest ICC value and lower bound of the 95% confidence interval was chosen as the preferred measurement protocol.

To assess inter-rater reliability for the two approaches, the second rater (S.N.) was trained directly by the initial rater step-by-step in the use of the software and chosen measurement technique as determined from the intra-examiner reliability assessment. The general use of the software, visualization of the tooth of interest in all three planes of space, automated segmentation procedures, manual refinements, and 3D visualization of the resultant volume were reviewed in training. The second rater (S.N.), who possessed a dental background and knowledge of normal dental anatomy, measured the same ten randomly selected maxillary first molars as measured by the principal investigator (both the whole tooth method and the dental root apical to the CEJ method). ICC was used to measure the agreement between the principal investigator's second measurement, as determined randomly, and the additional investigator's single measurement. A single measure with absolute agreement under two-way mixed model was chosen to account for rater variation, and the subjects/teeth were chosen randomly with the raters fixed.

Results

Intra-rater reliability

Protocol 1: manual human segmentation on a repeated 2D basis

The ICC demonstrated agreement, ICC (single measures) = 0.885, 95% CI (0.707, 0.967), within rater for the whole tooth measurement. The ICC demonstrated agreement, ICC (single measures) = 0.904, 95% CI (0.749, 0.973), within rater for the root measurement apical to the CEJ.

Protocol 2: automated segmentation without human refinement

The ICC demonstrated agreement, ICC (single measures) = 0.826, 95% CI (0.697, 0.952), within rater for the whole tooth measurement. The ICC demonstrated agreement, ICC (single measures) = 0.899, 95% CI (0.742, 0.953), within rater for the root measurement apical to the CEJ.

Protocol 3: automated segmentation with manual human refinement on a repeated 2D basis

The ICC demonstrated excellent agreement, ICC (single measures) = 0.996, 95% CI (0.989, 0.999), within rater for the whole tooth measurement. The ICC demonstrated agreement, ICC (single measures) = 0.904, 95% CI (0.751, 0.973), within rater for the root measurement apical to the CEJ.

Therefore, the whole tooth measurement utilizing protocol 3 was selected as the measurement method as it possessed the highest ICC value and lower bound of the confidence interval (ICC (single measures) = 0.996, 95% CI (0.989, 0.999)) when compared to all other measurement protocols investigated.

The summary of intra-rater reliability via the ICC is presented in Table 1. In addition, a summary of the largest volume differences for intra-rater repeated measures is presented in Table 2.

Table 1 Intra-rater intra-class correlation coefficient values
Table 2 Largest volume differences for intra-rater repeated measures

Inter-rater reliability

Looking strictly at the variability on an intra-rater basis was the focus of our determination of appropriate methods to be evaluated on an inter-rater basis. Therefore, the method with the highest intra-rater reliability was chosen to further address inter-rater reliability. Given that protocol 3, automated segmentation with manual human refinement on a repeated 2D basis, yielded the highest intra-rater reliability statistics, the inter-rater reliability was computed utilizing measurement protocol 3. The ICC demonstrates excellent agreement, ICC (single measures) = 0.990, 95% CI (0.961, 0.998), between raters for the whole tooth measurement. However, the ICC demonstrates less powerful agreement, ICC (single measures) = 0.728, 95% CI (0.198, 0.926), between raters for the root measurement apical to the CEJ. The inter-rater analysis results are in agreement with the chosen measurement protocol as determined via intra-rater ICC values.

The reliability readings for protocol 3 are included in Table 3. It serves to display the differences in absolute volume measurements for the repeated measures and inter-rater values.

Table 3 Reliability readings for protocol 3 (all units in mm 3 )

One subject had coronal restorative material present in the evaluated tooth. The presence of this restorative material did not have significant effects on the segmentation results, as it was not an outlier in the data set.

Discussion

The method involving automated segmentation with manual human refinement on a repeated 2D basis for the whole tooth proved to be the most reliable measurement protocol both within and between observers. Essentially, the intra- and inter-rater analysis results are in agreement with measurement protocol 3 as determined via ICC values. For excellent agreement, the ICC 95% confidence interval should be above 0.750 [35],[36], which is the case for the results obtained for protocol 3 using the entire tooth volume. It is of value to note that the protocol developed possesses similarities to studies investigating the accuracy of dental volume measurement in vivo using CBCT [21] and condylar head volume measurement [37], and hence lends to our segmentation technique's credibility.

The greatest difference across intra-rater repeated measures for the whole tooth approach utilizing protocol 3 was 17.76 mm3 (approximately 1.50% of the average whole tooth volume measured), whereas for the roots, only approach utilizing protocol 3 was 64.79 mm3 (approximately 11.45% of the average root volume apical to the CEJ measured). The intra-rater variability was thus approximately 3.6 times greater in absolute volume (64.79 versus 17.76 mm3) and 7.6 times greater in proportion of structure measured (11.45% versus 1.50%) when measuring roots only as compared to measuring the entire tooth volume.

The visualization of the respective maximum volume difference as displayed in Table 3 applied to a single tooth in various scenarios is displayed in Figure 6. The maximum volumetric discrepancy between repeated intra-observer measurements for the whole tooth was 17.76 mm3. The maximum volumetric discrepancy for inter-observer measurements for the whole tooth was 30.39 mm3. The effects of these measurement variations can be shown visually in a number of ways as displayed in Figure 6. Removal of the maximum inter-observer volume difference (30.39 mm3) strictly from the most apical portion of the palatal root (Figure 6B), from the apical portions of all three roots (Figure 6C), and from the buccal surfaces of the mesiobuccal and distobuccal roots are displayed (Figure 6D). Nearly imperceptible changes when differences are distributed across all roots and on the buccal surfaces of the mesiobuccal and distobuccal roots visually display the inter-observer errors with which tooth volumes may be determined. Visualizing the maximum volumetric discrepancy between repeated intra-observer measurements for the whole tooth of 17.76 mm3 is displayed visually in Figure 7. The volume displayed is the mesiobuccal cusp tip of the maxillary right first molar. The occlusal-apical dimension of the cusp tip volume measures only 1.2 mm, thus providing an approximate clinical interpretation and visualization of the volume differences. In addition, when considering the measurements of ERR with tooth-anchored maxillary expander (TAME) for the maxillary first molars completed by one rater, average maxillary first molar ERR volume changes of 42.67 mm3 have been previously reported in the literature [24], which is approximately 2.4 times greater than the intra-rater reliability protocol established for this technique.

Figure 6
figure 6

Three-dimensional visualization of inter-observer volume differences for the whole tooth measurement. (A) Entire tooth volume. (B) Entire tooth volume with maximum inter-observer variability removed from palatal root. (C) Entire tooth volume with maximum inter-observer variability removed from apical portion of all three roots. (D) Entire tooth volume with maximum inter-observer variability removed from buccal surfaces of the mesiobuccal and distobuccal roots.

Figure 7
figure 7

Three-dimensional visualization of the intra-observer volume difference for the maxillary right first molar mesio-buccal cusp tip.

The method resulting in the worst reliability was automated segmentation without human refinement. There are numerous reasons why this protocol was flawed. The determination of the boundaries between the tooth roots and the buccal and palatal cortical plates is sometimes indistinct given the very close proximity of the roots. The furcation area of the tooth possesses a large surface area of lamina dura, the dense surrounding cortical bone, which can lead to unclear tooth furcation anatomy. The proximity of the erupting second molar in some patients, as well as the interproximal contacts with adjacent teeth, often led to over-contouring of the volume of the crown of the maxillary first molar. In addition, the presence of dense bone islands of increased radiopacity can also result in misidentification of the proper root morphology. Given the limitations associated with a strictly automated method, it is still not possible, at least at a 0.25 mm voxel size, to automate the segmentation procedure. To be precise, the process must still involve manual human refinements with proper knowledge and interpretation of the 3D anatomy. As such, the process is extremely labor intensive given the slice-by-slice refinements that are required in all three planes of space.

In the approaches investigated, the pulpal tissue was included in the volumes as additional errors in delineating dentin from pulpal tissue would be an added source of variation. The additional dentin/soft tissue border, which is likely more challenging anatomy than the tooth to surrounding bony support to identify, due to intricate pulpal canal architecture of small dimension, would require identification. Therefore, since our area of interest is only ERR, internal pulpal changes are irrelevant. Consequently, both the hard and soft tissues within the cementum of the tooth were calculated as a part of the total tooth volume.

Due to the retrospective nature of the study, certain limitations were imposed on our ability to verify the accuracy of the volume measurements. To address this concern, focus was turned to the precision of volume measurements from both an intra- and inter-examiner perspective. The limitations from the retrospective nature of the study are twofold. Firstly, the CBCT machine used to capture the initial images was no longer functional or available for additional measurements such as in vitro dental volume comparisons. Secondly, due to the non-extraction orthodontic treatment of these patients, and that investigation of maxillary first molar volume was desired, the true value of the molar volumes is unknown and is unlikely to be known in future studies as maxillary first molars are not commonly extracted for orthodontic purposes except in rare circumstances. However, there appears no obvious reason for not being able to extrapolate the measurement protocol identified to other teeth. Therefore, with the inability to focus on the validity of the data, the approach was chosen to identify a measurement protocol to give highly precise results, both in intra- and inter-rater conditions.

The desire and ability to verify the true volumes or accuracy for this particular CBCT machine brings into question the capability to replicate an in vivo scenario. Numerous factors could not be addressed in a post hoc replicated model including the lack of a periodontal ligament, cortical bone, soft tissues, and patient movement to name a few. In addition, the imaging of a model as opposed to an in vivo dental volume followed by extraction and in vitro dental volume measurement would lead to the introduction of several errors and inaccuracies.

The validation of CBCT as a tool for measurement of both root lengths and volumes has been addressed in a number of studies with a multitude of image segmentation techniques. [16],[17],[19]-[21],[27]-[33] The weaknesses of the studies include the lack of investigation into more than one image segmentation protocol to provide the most precise experimental data. As an example, the study assessing ERR with maxillary expansion using 3D CBCT images [24] utilizes a segmentation procedure employing a root-only approach (apical to the CEJ for maxillary premolars and apical to the deepest point of the furcation for maxillary molars). The results of our study yielded the greatest intra-rater variability when using a similar method. Therefore, assessment of the root volume only appears to be wrought with errors in identification of the desired volume. Although an identical CBCT machine was not utilized in our own study, the voxel size and imaging parameters were similar to that of another [21]. In general, a change of software or CBCT machine appears to not be significantly clinically important given voxel sizes are held constant. There exist other 3D software programs for analyzing CBCT data with similar functions as the software is being utilized only as a tool to compute a volume. What does appear important however is voxel size and segmentation protocol, not the particular software used, as long as there is segmentation functionality.

Due to the uniqueness of our data set and limited access to the original CBCT machine because of the retrospective nature of the study, validation of our method, was sought in the literature. After the independent development, reliability testing, and subsequent comparison with existing published literature employing image segmentation protocols for dental volumes, some conclusions were reached. With numerous segmentation protocols in the literature, the volume measurement techniques in one CBCT volume validation study [21] were identical (in so far as can be determined from their reported methods) to our protocol 3 (whole tooth), which possessed the most precise volume segmentation results. Given the similarity of our measurement protocols, we feel confident in the validity of our results obtained to the study that verified the accuracy and validity of the dental tooth volume to within −4% to 7% [21].

Traditionally, as reported in the literature, bicuspids were routinely measured in vivo and subsequently extracted for physical volumetric determination [21]. There is an inherent tooth type limitation likely to be present in all studies due to the rarity of maxillary first molar extractions associated with orthodontic treatment. In our study, investigation of the maxillary first molar was chosen for a number of reasons including its complex root anatomy, early eruption and completion of root development in orthodontic-aged adolescent patients, and its use as an anchorage unit for initial phase orthodontic care. The potential for incomplete root development would be a limitation in evaluating any permanent tooth in adolescents, but given the comparatively early eruption of the permanent first molars, this limitation is mitigated as much as possible. An additional reason for this decision is because utilizing CBCT to assess ERR associated with maxillary expansion appliances is the ultimate goal. With maxillary first molars being the most commonly banded teeth that are attached to a maxillary expansion appliance, it makes inherent sense to assess the resorptive changes occurring within the anchor teeth themselves.

A strength of our reliability investigation lies not only in the numerous image measurement protocols investigated but also in the investigation of the entire maxillary first molar tooth volume versus the volume of the roots apical to the CEJ for each protocol as different segmentation cutoff points have been reported in the literature with no justification [24] or mention of technique accuracy or reproducibility. The reason two volumes were investigated is numerous. If the entire tooth volume was used, this adds the potential for patient coronal tooth volume variability between time points due to possible attrition, decay, loss of coronal tooth structure, and the placement or adjustment of new fillings or occlusion, in addition to other unidentified sources. This coronal tooth structure variability was not directly investigated in this particular paper, as the volume measurements for reliability were repeated measures on the same teeth at the same point in time. However, the identification and attention to the possible sources of variability aids in deciding whether the variability present in the root volume only measurements apical to the CEJ approach is more favorable due to the absence of the coronal variability. To address this from a visual perspective, Figure 8 displays the coronal changes that would have to be present on the cusp tips of the maxillary first molars to represent the maximum difference in additional variability of the root versus whole tooth (64.79 mm3 − 30.39 mm3 = 34.40 mm3) in the repeated measures. With visualization of the hypothetical coronal changes, it appears clear that measuring the whole tooth appears superior when compared to the roots only given the additional variability associated with measuring the roots only is more than can be expected from coronal changes over the short term (1 to 2 years). Support for the use of the whole tooth measurement protocol comes from a study exploring crown and root length of teeth using CBCT images. The study found a wider range of limits, and hence more variability, in measuring root lengths as opposed to crown lengths [26]. This increased error can perhaps be extended to root volumes due to difficulties in determining the CEJ location as the enamel is at its thinnest in this area. The CEJ is anatomically not a straight line; however, in many images, the lack of definition of the apical extent of the enamel resulted in a nearly straight-line resultant 3D segmentation. The best suggested method considers all the tooth structure not just the root in order to eliminate CEJ identification. For instance, an error in identification of one axial 2D slice was found to introduce root volume changes between 40 to 65 mm3 depending on tooth size. The increased variability appears to occur due to the fact that the CEJ represents quite a large 2D axial volume. In contrast, attrition at the molar cusp tips leads to almost imperceptible changes in tooth volume. The differences between a cusp tip axial slice area versus a CEJ slice area are displayed in Figure 9. Therefore, due to the increased difficulties in CEJ identification, the role of external coronal volume changes from T1 to T2 were judged to be minimal compared to the effects of an error in CEJ identification.

Figure 8
figure 8

Cusp tip attrition. This would have to be present to represent the maximum difference in additional variability between roots only and whole tooth in the repeated measures. (A) Normal 3D volume. (B) 3D volume representing maximum difference in additional variability of the roots versus the whole tooth (64.79 mm3 − 30.39 mm3 = 34.40 mm3) in the repeated measures removed from cusp tips.

Figure 9
figure 9

Axial slice area at the cement-enamel junction (CEJ) versus cusp tip. (A) Axial slice area at the molar CEJ. (B) Axial slice area at the molar cusp tips.

In general, the presence of radiodense restorative materials, such as amalgam and some highly filled composites, which greatly inhibit the passage of electromagnetic radiation, has the potential to introduce further variability into dental volume determination and result in imaging artifacts. These artifacts can result in the inability to predictably identify the true extent of the radiodense material or adjacent structures and thus affect volume segmentation.

There exists the obvious issue of resolution of a CBCT image in determining the volume of a tooth. Given the voxel size of 0.25 mm in each dimension, there are concerns regarding the potential that the border of a tooth versus bone could be contained within a voxel. A limitation in computed tomography imaging is the so-called ‘partial volume effect.’ In essence, this phenomenon can present issues in the differentiation of different tissue types. For example, a large amount of periodontal ligament (PDL) space and a thin layer of compact bone such as the lamina dura can cause the same attenuation in a voxel as the dentin of a tooth alone. The issue of resolution is complex as improved resolution can be acquired, but at the expense of increased patient radiation dose [30]. Although improved volumetric determination and ERR detection can be obtained with smaller voxel sizes and increased scan times [28], there exists a limit, which has to be established between patient radiation dose and resolution required for appropriate diagnostics. However, early detection of root changes may modify treatment mechanics and thus limit the progression of ERR and the long-term impact on the affected teeth. Using a voxel size of 0.125 mm has the potential to yield in vivo volume measurement of teeth comparable to micro-CT in vitro analysis [20], but understandably has the disadvantage of increased patient radiation dose. A study, which investigated the influence of voxel size on the diagnostic ability of CBCT to evaluate ERR, concluded CBCT to be a reliable method of ERR detection with a voxel size of 0.3 mm as the ‘best protocol’ when balancing patient dose and diagnostic performance [38]. Given the voxel size of 0.25 mm used in this study, there appears to be more than adequate resolution required to measure dental tooth volumes while balancing patient radiation dose.

A limitation of the study involves the use of only one CBCT model to capture the patient images. There exists the potential issue of variation between different CBCT models that may possibly possess varying image quality and gray value distributions for the aforementioned segmentation methods.

Through the establishment of precise and accurate dental volume measurement protocols, clinicians can more confidently employ the available tools to monitor such phenomena as ERR at various stages throughout the orthodontic process and ERR related to dental trauma. However, patient radiation exposure and diagnostic imaging needs require a careful balance to be established. The uses of CBCT imaging are to maximize the diagnostic information available to the clinician while limiting patient radiation exposure to make individualized treatment decisions while considering as many patient specific factors as possible.

Conclusions

The proposed maxillary first molar dental volume measurement protocol for CBCT images employing automated segmentation with manual human refinement on a repeated 2D slice-by-slice basis in all three planes of space possesses excellent intra- and inter-rater reliability and precision.

Maxillary first molar 3D volume measurements of the entire tooth structure are more precise than 3D volume measurements of only the dental roots apical to the CEJ.

References

  1. de Freitas MR, Beltrao RT, Janson G, Henriques JF, Chiqueto K: Evaluation of root resorption after open bite treatment with and without extractions. Am J Orthod Dentofacial Orthop. 2007, 132(2):143.e15–143.e22.

    Article  Google Scholar 

  2. Levander E, Malmgren O, Eliasson S: Evaluation of root resorption in relation to two orthodontic treatment regimes. A clinical experimental study. Eur J Orthod. 1994, 16(3):223–8. 10.1093/ejo/16.3.223

    Article  PubMed  Google Scholar 

  3. Levander E, Bajka R, Malmgren O: Early radiographic diagnosis of apical root resorption during orthodontic treatment: a study of maxillary incisors. Eur J Orthod. 1998, 20(1):57–63. 10.1093/ejo/20.1.57

    Article  PubMed  Google Scholar 

  4. Brezniak N, Wasserstein A: Root resorption after orthodontic treatment: part 1. Literature review. Am J Orthod Dentofacial Orthop. 1993, 103(1):62–6. 10.1016/0889-5406(93)70106-X

    Article  PubMed  Google Scholar 

  5. Mirabella AD, Artun J: Prevalence and severity of apical root resorption of maxillary anterior teeth in adult orthodontic patients. Eur J Orthod. 1995, 17(2):93–9. 10.1093/ejo/17.2.93

    Article  PubMed  Google Scholar 

  6. Chapnick L: External root resorption: an experimental radiographic evaluation. Oral Surg Oral Med Oral Pathol. 1989, 67(5):578–82. 10.1016/0030-4220(89)90276-4

    Article  PubMed  Google Scholar 

  7. Forst D, Nijjar S, Khaled Y, Manuel L, Flores-Mir C: Radiographic assessment of external root resorption associated with jackscrew-based maxillary expansion therapies: a systematic review. Eur. J. Orthod. 2013, 36: 576–85. 10.1093/ejo/cjt090

    Article  PubMed  Google Scholar 

  8. Sameshima GT, Sinclair PM: Predicting and preventing root resorption: part I. Diagnostic factors. Am J Orthod Dentofacial Orthop. 2001, 119(5):505–10. 10.1067/mod.2001.113409

    Article  PubMed  Google Scholar 

  9. Chan EK, Darendeliler MA: Exploring the third dimension in root resorption. Orthod Craniofac Res. 2004, 7(2):64–70. 10.1111/j.1601-6343.2004.00280.x

    Article  PubMed  Google Scholar 

  10. Katona TR: Flaws in root resorption assessment algorithms: role of tooth shape. Am J Orthod Dentofacial Orthop. 2006, 130(6):698.e19–698.e27.

    Article  Google Scholar 

  11. Mckee IW, Williamson PC, Lam EW, Heo G, Glover KE, Major PW: The accuracy of 4 panoramic units in the projection of mesiodistal tooth angulations. Am J Orthod Dentofacial Orthop. 2002, 121(2):166–75. quiz 192 quiz 192 10.1067/mod.2002.119435

    Article  PubMed  Google Scholar 

  12. Van Elslande DC, Russett SJ, Major PW, Flores-Mir C: Mandibular asymmetry diagnosis with panoramic imaging. Am J Orthod Dentofacial Orthop. 2008, 134(2):183–92. 10.1016/j.ajodo.2007.07.021

    Article  PubMed  Google Scholar 

  13. Darendeliler MA, Kharbanda OP, Chan EK, Srivicharnkul P, Rex T, Swain MV, Jones AS, Petocz P: Root resorption and its association with alterations in physical properties, mineral contents and resorption craters in human premolars following application of light and heavy controlled orthodontic forces. Orthod Craniofac Res. 2004, 7(2):79–97. 10.1111/j.1601-6343.2004.00281.x

    Article  PubMed  Google Scholar 

  14. Garcia-Figueroa MA, Raboud DW, Lam EW, Heo G, Major PW: Effect of buccolingual root angulation on the mesiodistal angulation shown on panoramic radiographs. Am J Orthod Dentofacial Orthop. 2008, 134(1):93–9. 10.1016/j.ajodo.2006.07.034

    Article  PubMed  Google Scholar 

  15. Ponder SN, Benavides E, Kapila S, Hatch NE: Quantification of external root resorption by low- vs high-resolution cone-beam computed tomography and periapical radiography: a volumetric and linear analysis. Am J Orthod Dentofacial Orthop. 2013, 143(1):77–91. 10.1016/j.ajodo.2012.08.023

    Article  PubMed  Google Scholar 

  16. Ericson S, Kurol J: Incisor root resorptions due to ectopic maxillary canines imaged by computerized tomography: a comparative study in extracted teeth. Angle Orthod. 2000, 70(4):276–83.

    PubMed  Google Scholar 

  17. Lund H, Grondahl K, Grondahl HG: Cone beam computed tomography for assessment of root length and marginal bone level during orthodontic treatment. Angle Orthod. 2010, 80(3):466–73. 10.2319/072909-427.1

    Article  PubMed  Google Scholar 

  18. Patel S, Dawood A, Wilson R, Horner K, Mannocci F: The detection and management of root resorption lesions using intraoral radiography and cone beam computed tomography - an in vivo investigation. Int Endod J. 2009, 42(9):831–8. 10.1111/j.1365-2591.2009.01592.x

    Article  PubMed  Google Scholar 

  19. Li W, Chen F, Zhang F, Ding W, Ye Q, Shi J, Fu B: Volumetric measurement of root resorption following molar mini-screw implant intrusion using cone beam computed tomography. PLoS One 2013, 8(4):e60962. 10.1371/journal.pone.0060962

    Article  PubMed Central  PubMed  Google Scholar 

  20. Wang Y, He S, Yu L, Li J, Chen S: Accuracy of volumetric measurement of teeth in vivo based on cone beam computer tomography. Orthod Craniofac Res 2011, 14(4):206–12. 10.1111/j.1601-6343.2011.01525.x

    Article  PubMed  Google Scholar 

  21. Liu Y, Olszewski R, Alexandroni ES, Enciso R, Xu T, Mah JK: The validity of in vivo tooth volume determinations from cone-beam computed tomography. Angle Orthod. 2010, 80(1):160–6. 10.2319/121608-639.1

    Article  PubMed  Google Scholar 

  22. Ye N, Jian F, Xue J, Wang S, Liao L, Huang W, Yang X, Zhou Y, Lai W, Li J, Wang J: Accuracy of in-vitro tooth volumetric measurements from cone-beam computed tomography. Am J Orthod Dentofacial Orthop. 2012, 142(6):879–87. 10.1016/j.ajodo.2012.05.020

    Article  PubMed  Google Scholar 

  23. Maret D, Telmon N, Peters OA, Lepage B, Treil J, Inglese JM, Peyre A, Kahn JL, Sixou M: Effect of voxel size on the accuracy of 3D reconstructions with cone beam CT. Dentomaxillofac Radiol. 2012, 41(8):649–55. 10.1259/dmfr/81804525

    Article  PubMed Central  PubMed  Google Scholar 

  24. Baysal A, Karadede I, Hekimoglu S, Ucar F, Ozer T, Veli I, Uysal T: Evaluation of root resorption following rapid maxillary expansion using cone-beam computed tomography. Angle Orthod. 2012, 82(3):488–94. 10.2319/060411-367.1

    Article  PubMed  Google Scholar 

  25. Ramires T, Maia RA, Barone JR: Nasal cavity changes and the respiratory standard after maxillary expansion. Rev Bras Otorrinolaringol (Engl Ed). 2008, 74(5):763–9. 10.1590/S0034-72992008000500019

    Article  Google Scholar 

  26. Kim SY, Lim SH, Gang SN, Kim HJ: Crown and root lengths of incisors, canines, and premolars measured by cone-beam computed tomography in patients with malocclusions. Korean J Orthod. 2013, 43(6):271–8. 10.4041/kjod.2013.43.6.271

    Article  PubMed Central  PubMed  Google Scholar 

  27. Wang Y, He S, Guo Y, Wang S, Chen S: Accuracy of volumetric measurement of simulated root resorption lacunas based on cone beam computed tomography. Orthod Craniofac Res. 2013, 16(3):169–76. 10.1111/ocr.12016

    Article  PubMed  Google Scholar 

  28. Dalili Z, Taramsari M, Mousavi Mehr SZ, Salamat F: Diagnostic value of two modes of cone-beam computed tomography in evaluation of simulated external root resorption: an in vitro study. Imaging Sci Dent. 2012, 42(1):19–24. 10.5624/isd.2012.42.1.19

    Article  PubMed Central  PubMed  Google Scholar 

  29. Shokri A, Mortazavi H, Salemi F, Javadian A, Bakhtiari H, Matlabi H: Diagnosis of simulated external root resorption using conventional intraoral film radiography, CCD, PSP, and CBCT: a comparison study. Biomed J. 2013, 36(1):18–22. 10.4103/2319-4170.107156

    Article  PubMed  Google Scholar 

  30. Ren H, Chen J, Deng F, Zheng L, Liu X, Dong Y: Comparison of cone-beam computed tomography and periapical radiography for detecting simulated apical root resorption. Angle Orthod. 2013, 83(2):189–95. 10.2319/050512-372.1

    Article  PubMed  Google Scholar 

  31. Castro IO, Alencar AH, Valladares-Neto J, Estrela C: Apical root resorption due to orthodontic treatment detected by cone beam computed tomography. Angle Orthod. 2013, 83(2):196–203. 10.2319/032112-240.1

    Article  PubMed  Google Scholar 

  32. Xie XY, Zhang ZY: [Diagnostic accuracy of cone beam computed tomography and eight-slice computed tomography for evaluation of external root reabsorption]. Beijing Da Xue Xue Bao. 2012, 44(4):628–32.

    PubMed  Google Scholar 

  33. Lund H, Grondahl K, Hansen K, Grondahl HG: Apical root resorption during orthodontic treatment. A prospective study using cone beam CT. Angle Orthod. 2012, 82(3):480–7. 10.2319/061311-390.1

    Article  PubMed  Google Scholar 

  34. Lagravere MO, Carey J, Heo G, Toogood RW, Major PW: Transverse, vertical, and anteroposterior changes from bone-anchored maxillary expansion vs traditional rapid maxillary expansion: a randomized clinical trial. Am J Orthod Dentofacial Orthop. 2010, 137(3):304.e1–304.12. discussion 304–5 discussion 304–5

    Google Scholar 

  35. Bland JM, Altman DG: Statistical methods for assessing agreement between two methods of clinical measurement. Lancet 1986, 1(8476):307–10. 10.1016/S0140-6736(86)90837-8

    Article  PubMed  Google Scholar 

  36. Fleiss JL: Reliability of measurement. In The Design and Analysis of Clinical Experiments. John Wiley & Sons, Inc., New York; 1999:1–32. 10.1002/9781118032923

    Chapter  Google Scholar 

  37. Xi T, van Loon B, Fudalej P, Berge S, Swennen G, Maal T: Validation of a novel semi-automated method for three-dimensional surface rendering of condyles using cone beam computed tomography data. Int. J. Oral Maxillofac. Surg. 2013, 42(8):1023–9. 10.1016/j.ijom.2013.01.016

    Article  PubMed  Google Scholar 

  38. Liedke GS, da Silveira HE, da Silveira HL, Dutra V, de Figueiredo JA: Influence of voxel size in the diagnostic ability of cone beam tomography to evaluate simulated external root resorption. J. Endod. 2009, 35(2):233–5. 10.1016/j.joen.2008.11.005

    Article  PubMed  Google Scholar 

Download references

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Manuel Lagravere.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

DF: Conception, design, data acquisition, analysis, and interpretation of data. Developed data acquisition methods and completed intra-examiner reliability readings. Drafted the written work. Final approval for publishing. SN: Design, data acquisition, analysis, and interpretation of data. Completed inter-examiner reliability readings. Assisted in drafting the written work. Critical revision of the written work. Final approval for publishing. CFM, JPC and MS: Conception, design and interpretation of data. Critical revision of the written work. Final approval for publishing. MOL: Conception, design, assistance with data acquisition and interpretation of data. Critical revision of the written work. Final approval for publishing. All authors read and approved the final manuscript.

Authors’ original submitted files for images

Rights and permissions

Open Access  This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made.

The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.

To view a copy of this licence, visit https://creativecommons.org/licenses/by/4.0/.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Forst, D., Nijjar, S., Flores-Mir, C. et al. Comparison of in vivo 3D cone-beam computed tomography tooth volume measurement protocols. Prog Orthod. 15, 69 (2014). https://doi.org/10.1186/s40510-014-0069-2

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s40510-014-0069-2

Keywords