Can Dental Cone Beam Computed Tomography Assess Bone Mineral Density?
Article information
Abstract
Mineral density distribution of bone tissue is altered by active bone modeling and remodeling due to bone complications including bone disease and implantation surgery. Clinical cone beam computed tomography (CBCT) has been examined whether it can assess oral bone mineral density (BMD) in patient. It has been indicated that CBCT has disadvantages of higher noise and lower contrast than conventional medical computed tomography (CT) systems. On the other hand, it has advantages of a relatively lower cost and radiation dose but higher spatial resolution. However, the reliability of CBCT based mineral density measurement has not yet been fully validated. Thus, the objectives of this review are to discuss 1) why assessment of BMD distribution is important and 2) whether the clinical CBCT can be used as a potential tool to measure the BMD. Brief descriptions of image artefacts associated with assessment of gray value, which has been used to account for mineral density, in CBCT images are provided. Techniques to correct local and conversion errors in obtaining the gray values in CBCT images are also introduced. This review can be used as a quick reference for users who may encounter these errors during analysis of CBCT images.
Why assessment of bone mineral density (BMD) distribution is important?
Bone is a connective tissue designed to bear loading generated by muscle movement and body weight during daily activity.[1,2] Bone in a live skeletal system constantly changes. Approximately 7.5% of bone turnover occurs annually.[3] Bone modeling and remodeling are inherently involved in the alteration of bone. In the process of bone remodeling, a bone resorbing cell, osteoclast, is activated to remove the pre-existing bone tissue and subsequently, new bone is deposited by a bone forming cell, osteoblast.[3,4,5,6,7] The bone modeling is an uncoupled process consisting of activation and resorption or activation and formation. When bone resorption and formation are balanced, the net quantity of bone is maintained. It has been observed that bone quantity is associated with its mechanical strength.[8,9] As bone in a live patient cannot be subject to mechanical testing, many non-invasive imaging technologies have been developed to estimate the bone quantity as a surrogate for a patient's bone fragility. Computed tomography (CT) has been widely used non-invasive device to assess bone mass.[10,11] This X-ray based tool mainly detects mineral, which is the densest component in bone tissue materials, and provides Hounsfield units (HU) as the measures. It is well known that bone is composed of approximately 30% water, 40% mineral, and 30% organic protein (mostly collagen) at the tissue level.[1] In the initial process of bone formation, osteoblasts produce immature collagen and mineral composite.[12] The mineral content rapidly increases up to 70% of full mineralization within one week post-initiation of bone formation.[13] Following the first mineralization, more minerals are added to the collagen-mineral composite during the secondary mineralization process that lasts for years.[13,14] Bone modeling and remodeling processes are activated at different time points. In particular, active bone remodeling is stimulated by bone disease and complications including estrogen deficiency in postmenopausal osteoporosis, osteoarthritis, fracture healing, and implantation surgery.[2,6,7,15,16,17,18] As a result, heterogeneity of bone tissue mineral density (TMD) increases (Fig. 1).

(A) Detailed tissue mineral density (TMD) distribution in vertebral trabecular bone. A darker color represents less TMD. (B) A typical TMD histogram of a micro-computed tomography image (voxel size 16×16×16 µm3). The TMD distribution was different between the control sham surgery (black) and ovariectomized (OVX) (gray) groups. [Reprinted from "Increased variability of bone tissue mineral density resulting from estrogen deficiency influences creep behavior in a rat vertebral body", by Kim DG, Navalgund AR, Tee BC, Noble GJ, Hart RT, Lee HR, 2012, Bone, 51(5), pp. 868-75. Copyright 2012 by the Elsevier. Reprinted with permission].
BMD accounts for bone quantity as the mineral contents within an apparent volume of bone, which includes porosity, bone marrow, as well as bone matrix. On the other hand, the TMD represents mineral contents only in the matrix of bone.[18,19] Bone quantity measurements have been widely used to diagnose bone disease. For example, according to the World Health Organization (WHO) criteria, osteoporosis is defined as 2.5 standard deviations or more below the average value for young healthy women.[20,21] However, it has been indicated that bone quantity alone cannot explain bone fragility.[20,22] On the other hand, the bone quality is assessed by its morphology and TMD distribution.[23,24] Mechanical properties of trabecular bone are explained using three-dimensional (3D) morphological parameters of trabecular bone such as trabecular number, thickness, and anisotropy ratio.[25,26,27] Recently, many studies have also indicated that TMD is responsible for controlling elastic, plastic, and viscoelastic mechanical properties at the tissue level of bone.[28,29,30,31] As such, efforts have been made to develop a new imaging technology that can non-invasively assess the bone quality. Dual energy X-ray absorptiometry (DXA) has been recognized as a standard technique to obtain the BMD of a patient's spine or hip clinically. However, the DXA produces a low resolution (approximately 500 µm [32]) 2D image that only provides areal BMD and cannot precisely delineate cancellous bone structure with the trabecular thickness level of 50 to 200 µm.[25,33,34] A high-resolution (up to less than 1 µm) 3D micro-CT image can describe detailed bone morphology and TMD more accurately than other CT images (Fig. 2). However, much higher radiation doses are generated during the longer scanning time of micro-CT for the same size of a specimen, which limits its use to laboratory research purposes only. On the other hand, cone beam CT (CBCT) provides higher image resolution with less radiation dose and exposure time than the conventional clinical multidetector CT (MDCT) (Table 1). Thus, CBCT has been widely used for the clinical diagnosis of dental complications.[35,36] The 3D image produced by CBCT allows for detailed morphological analysis of bone. However, reliability of CBCT based bone density measurement remains to be determined.[37,38,39] Thus, this review will discuss whether the CBCT can be used to assess BMD at next sections.

(A) Micro-computed tomography (CT) image (27×27×27 µm3 voxel size) and (B) cone beam CT image (200×200×200 µm3 voxel size) of the same human condyle.
1. How a CBCT image is created?
The CBCT is relatively recent technology, with the first commercial scanner being introduced in 1998 for dental imaging.[40] This is a type of CT that uses a cone X-ray beam instead of the conventional fan beam (Fig. 3).[35] For standard full field of view for CBCT, polychromatic X-rays are generated from a source tube excited by 100-120 kV tube potential with 1.5-10 mA current.[41] The emitted X-ray photon energy is absorbed by materials located between the X-ray source and detector. The remaining energy after absorption is recorded on a detector as expressed by intensity (I) that represents energy per area and time. For example, the flat-panel detector consists of a pixel array of hydrogenated amorphous silicon thin-film transistors.[35] A scintillating material (e.g., Gd2O2S) detects the X-rays and converts them into visible light that is charged in the pixel transistors. A CT attenuation coefficient value of the material can be obtained by the Lambert-Beer law (Eq. 1).[39]

X-ray beam projection scheme comparing acquisition geometry of conventional or "fan" beam (right) and "cone" beam (left) imaging geometry and resultant image production. The amount of scatter generated (sinusoidal lines) and recorded by cone-beam image acquisition is substantially higher, reducing image contrast and increasing image noise. [Reprinted from "What is cone-beam CT and how does it work?", by Scarfe WC, Farman AG, 2008, Dent Clin North Am, 52(4), pp.707-30. Copyright 2008 by the Elsevier. Reprinted with permission].
where I0 is the intensity emitted by the source, µ represents the attenuation coefficient, and l is a distance over which the attenuation is integrated. A denser material absorbs more energy resulting in a higher attenuation value. The attenuation values are converted to gray values in a digital (8-bit to 16-bit) image while the scanned image slices are reconstructed. As such, the gray value in the CBCT image is theoretically equivalent to the density of the material. However, in order to obtain a reliable density value using a CBCT image, multiple factors must be considered.
2. CBCT artefacts from machine factors
Imaging artefacts associated with general CT systems also occur in CBCT imaging.[39,42,43] Referring to the artefacts listed in previous literatures, this review summarizes the artefacts that can influence the CBCT based bone density measurements and related issues.
1) Noise
This artefact is appeared as inconsistent gray values with large standard deviations. This result arises from low signal-to-noise ratio of intensity, which needs to maintain the low radiation dose. The noise level can be reduced by increasing the excitation potential and current.
2) Scatter artefacts
Scattered X-ray photons from the original path can be added to the primary intensity giving rise to an underestimation of attenuation value in Eq. 1. Larger detectors have a greater chance of encountering the scattered X-ray photons leading to streak artefacts during reconstruction process for the CBCT image.
3) Beam hardening
As the energy levels of polychromatic X-ray beams used in the CBCT are not identical, the lower energy photons can be easily absorbed at the edges of the scanned subject resulting in hardening of the X-ray beam, which produces lower gray values toward the center of the subject (cupping artifact) even if the density of the subject is homogenous.
4) Ring artefacts
Defects or uncalibrated components in the detector may cause ring artefacts with concentric rings in the CBCT image. The ring voxels have inconsistent gray values that can increase overall errors in the assessment of bone density.
5) Partial volume effects
The cubic or rectangular voxels cannot completely delineate the irregular shapes of scanned subjects. Thus, the gray value of voxels at the border between different density materials contains averaged attenuations. If the voxel size of CBCT image increases, it has more incorrect partial volume gray values.
3. CBCT artefacts from software factors
1) Reconstruction algorithm
A sinogram is constructed using digital signals from the pixel transistors. This composite image combines each row of each projection in the CBCT image.[35] Then, a mathematical algorithm converts the sinogram to 2D image slices and reconstructs them to a 3D CBCT image. The most predominantly used algorithm is a convolution-backprojection formula introduced by Feldkamp et al. [35,44] It was indicated that this backprojection algorithm causes unavoidable distortion in the transverse direction and resolution degradation in the longitudinal direction.[35] The reconstruction process also converts the raw intensity values in the sinogram to the gray values in the scale range of image data. Thus, absolute gray values may be different depending on the data size of the CBCT image utilized by the operation software of a CBCT machine company.
2) Field of view (FOV)
CBCT can improve the image resolution by reducing the FOV while decreasing overall radiation dose. However, variability of gray values is increased by using a smaller FOV (5 cm) compared to larger FOVs (10 to 20 cm).[38] Many other studies have indicated the FOV associated effects on gray value variability.[45,46]
3) HU
The attenuation coefficients of the same material scanned by different CT systems can vary if the scanning conditions are not identical. As the HU can be computed relative to the attenuation coefficient of water (Eq. 2), it has been widely used to compare material density between different CT systems.
It was suggested that the HU should be calibrated in order to obtain a consistent density value when the same material is scanned using different CT systems.[47,48] However, accuracy of the CBCT based HU measurement still remains controversial.[37,38,39]
4. CBCT artefacts from patient factors
1) Streak artefact
Scanning dense metallic materials can cause severe streaking artifacts when their gray values exceed the maximum level of operation that the software can handle.[43] This artefact limits analysis of local gray values surrounding dental restorations and metal implants in patients. Other sources to cause streak artefacts include beam hardening, noise, and photon starvation that can appear when insufficient photons reach the detector.
Can the clinical CBCT be used as a potential tool to measure the BMD?
Noise, ring and streak artefacts can be corrected by improving reconstruction algorithms.[39,42,49] Recently, a scatter-correction algorithm was developed to significantly improve the consistency of CBCT gray values including compressed sensing reconstruction algorithms to create a scatter-free CBCT image.[50,51] As the beam hardening effect is the most frequently observed artefact in CT scanning, many efficient correction algorithms have been developed and included in commercial CT software.[52,53] The FOV associated variations may arise because the backprojection algorithm includes gray values of tissues outside the small FOV during reconstruction.[39] Taken together, the recently improved reconstruction process likely corrects the local errors of gray value in the CBCT image. However, as these advanced algorithms need to be verified and calibrated to individual CT systems prior to being utilized in the clinical settings, the artefacts may be still observed in most of the currently used CT machines.
While the local errors can be visually detected, there are some systematic complications to be considered for CBCT based mineral density measurement. The most debated aspect is that the HU values of subjects are not consistent between different CT systems and between different times scanned even using the same CT system. These discrepancies can arise from the non-uniform process of scaling the HU values during reconstruction. It was indicated that a manufacturer's software imports a gray value range of 1 to 4,096 and rescales from -1,024 to 3,072.[47] This image would have 12-bit data size (212=4,096). The gray value of an image with 8-bit data size is scaled from 0 to 255 (28=256). The HU of water and air should be close to 0 and -1,000, respectively. However, the different scaling processes may provide completely different absolute values of HU for the same material. As such, the value of each voxel in a CBCT image may better represent the gray value rather than the HU.
Many studies found that strong positive correlations of the gray values of CBCT image with known density of reference materials and gray values obtained from the conventional clinical CT.[47,48,54,55,56,57,58,59] For example, phantoms of bone materials (hydroxyapatite) with 3 different densities (1,000, 1,250, and 1,750 mg/cm3) were scanned using 3 different resolutions (200, 300, and 400 µm) of CBCT (Fig. 4). As a result, the CBCT based gray values of phantoms had strong positive correlations with the density values. The gray value of scanned specimens can be directly converted to the corresponding mineral density value using this calibration curve when the same CBCT scanner is used. Alternatively, a previous study found that the CBCT gray value strongly correlated with the CT attenuation coefficient (µ) of standard materials, which include aluminium (µ=0.678), outer bone equivalent (µ=0.523), inner bone equivalent (µ=0.255), polymethylmethacrylate (PMMA) (µ=0.220), muscle equivalent (µ=0.205), water equivalent (µ=0.202), adipose equivalent (µ=0.176), and air (µ=0.000).[47] Then, the corrected HU value of scanned materials was derived by estimating the attenuation values in the HU equation (Eq. 2) using the correlation equation with the CBCT gray values. The corrected HU values were comparable between different CBCT scanners. Recently, clinical applicability of this HU correction was validated using intraoral phantoms for patients.[48]

(A) Strong positive correlations in the calibration curves of gray values for (B) phantoms of bone materials (hydroxyapatite) with 3 different densities (1,000, 1,250, and 1,750 mg/cm3) scanned using 3 different resolutions (200, 300, and 400 µm) of cone beam computed tomography. HU, Hounsfield units.
The aforementioned methods need external standard phantoms to be scanned to obtain the calibration curve. Another methodology was introduced using an internal reference to compare the mineral density between patients.[60,61] In the previous study using CBCT images of human cadaveric mandibles,[61] teeth were digitally separated from the mandible segment using imaging software (Livewire®,[62] Institute of Computing, State University of Campinas, Brazil[61] and ImageJ, National Institutes of Health, Bethesda, Maryland, USA) (Fig. 5). The alveolar bone (AB) region was digitally isolated from within 1 mm of the root surface and the basal control bone (CB) was determined at 0.6 mm from periosteal and endosteal surfaces of the basal bone. The gray value histogram was used to obtain mean, standard deviation (SD), coefficient of variation (COV=SD/Mean), low and high gray values (Low5 and High5) at the 5th and the 95th percentiles of voxel counts in the histogram. The gray value parameters of the CB region were used as the internal references to compare those of AB region between subjects. It was assumed that a core region of basal cortical bone has normal bone remodeling while the AB and marginal basal bone regions have a rapid bone turnover due to active bone formation.[17,63] Percentage difference (%) of these gray level parameters between AB and basal cortical bone ([AB-CB]/[(AB+CB)/2]×100) was then computed for each patient CBCT image. The same image process was performed using the micro-CT images of the same human mandibles. As a result, the relative values of percentage differences were comparable between CBCT and micro-CT images based analyses. Those procedures to assess the gray value parameters were also applied to examine CBCT images collected from routine patients in clinic. An advantage of this method is that the relative value of the percentage difference is computed using the gray values of different regions in the same CBCT image avoiding the effects of different scanning conditions. The gray value histogram analyses were successfully used for clinical investigation of the success and failure of dental grafting and implant systems.[64,65]

(A) Degree of bone mineralization parameters determined using a grey level histogram, (B) comparison of grey level histograms between alveolar bone (AB, grey line) and basal cortical bone (CB, black line) regions using a three-dimensional (3D) cone beam computed tomography (CT) image (200×200×200 µm3 voxel size), and (C) using 3D micro-CT image (20×20×20 µm3 voxel size) for the same specimen. COV, coefficient of variation; Highs, grey level at the 95th percentile; LOWs, grey level at the 5th percentile; SD, standard deviation; AB, alveolar bone; CB, control bone. [Modified from "Comparison of micro-CT and cone beam CT-based assessments for relative difference of grey level distribution in a human mandible" by Taylor TT, Gans SI, Jones EM, Firestone AR, Johnston WM, Kim DG, 2013, Dentomaxillofac Radiol, 42(3), pp. 25117764. Copyright 2013 by British Institute of Radiology. Reprinted with permission].
To date, CBCT based BMD measurement has been used mainly to estimate bone properties for dental implantation.[66,67,68,69,70] The BMD measured from CBCT images showed a better predicting power to estimate cortical bone fracture than assessed by the conventional DXA images.[71] More applications are performed for the diagnosis of oral bone disease, improvement of orthodontic treatment, and maxillofacial surgery.[36,72,73,74] It is clear that higher resolution of 3D images can provide more useful information. For example, the partial volume effect can be reduced using a smaller voxel size in the image. However, if longer scanning time is required to obtain higher image resolution, the risk of high radiation exposure likely increases against the principle of "as low as reasonably achievable" (ALARA) radiation dose. As complete segmentation of bone voxels from non-bone voxels in the CT image with the current resolution can tremendously improve the capability of CT based analyses for both morphology and density (Fig. 2), numerous studies have been attempted to develop innovative algorithms.[75,76,77,78] Most commercialized CBCT software uses a global thresholding gray value to segment voxels. However, this simple method does not provide sufficient quality of images.[79] Further studies are still required to develop an optimal segmentation technology specific to individual CT images.
CONCLUSION
Mineral density distribution of bone tissue reflects the result of biological activity, which is altered due to bone complications. Clinical CT is an indirect method to measure mineral density distribution based on the X-ray attenuation coefficient of the materials, mainly the mineral, in bone tissue. Hence, CT based density measurement is not as accurate as direct measurement using biopsy. However, it is a very powerful non-destructive tool that allows for longitudinal diagnoses of patients' bone disease. Most of the local artefacts can be corrected by improving the CBCT image reconstruction algorithm and alternative corrective methods can be applied as visualized in the image. However, the image processing errors during converting the attenuation coefficient values to the HUs and gray values should be accurately addressed in order to compare the mineral density in CBCT images scanned under different conditions. This review briefly introduces the types of artefacts that occur and solutions to correct them as a quick reference for users who may find these errors in CBCT images.
Notes
No potential conflict of interest relevant to this article was reported.
The author thanks his student, Trenton B. Johnson, for helpful discussions about the subjects described in this review.