Analysis of Bone Quality on Panoramic Radiograph in Osteoporosis Research by Fractal Dimension

DOI: 10.4236/am.2016.74035   PDF   HTML   XML   1,897 Downloads   2,801 Views   Citations


Objective: To Assess the correlation between different quality analysis parameters of trabecular pattern in digital panoramic radiographies and relations with forearm bone mass density (BMD) performed by DXA. Methods: The study was developed using panoramic and peripheral bone densitometry dual energy X-ray absorptiometry (DXA) of 68 patients, 9 males and 59 females (19 - 73 years old). In the panoramic radiographs, evaluation of the trabecular bone morphology through assessment of fractal dimension (FD), connectivity (C) and total number of “bright” pixels (ET) was performed. In DXA, the exam determines the bone mineral density of the forearm to identify who has a high risk of osteoporosis. Statistics analyzed the relationship of these exams and the contribution of dental radiographs in detecting patients at risk for osteoporosis. Results: The average age of subjects was 43.85. In the analysis of trabecular pattern, a significant correlation between the FD, ET and C factors in level of 5% (Pearson correlation test) was found. Correlation tests showed no significant correlation between DF and BMD. Conclusions: The analysis showed correlations with each other, detecting alterations in the trabecular pattern. It cannot be related to BMD with FD but should be taken into account that examining the bone or trabecular alveolar process, when, for example, diagnostic analysis of pre-implant bone quality, is required.

Share and Cite:

Camargo, A. , Côrtes, A. , Aoki, E. , Baladi, M. , Arita, E. and Watanabe, P. (2016) Analysis of Bone Quality on Panoramic Radiograph in Osteoporosis Research by Fractal Dimension. Applied Mathematics, 7, 375-386. doi: 10.4236/am.2016.74035.

Received 27 November 2015; accepted 14 March 2016; published 17 March 2016

1. Introduction

Osteoporosis is a major health problem and affects a significant number of people. This skeletal disorder is characterized by bone fragility due to deterioration of the bone micro-architecture, which confers strength and bone quality; this can increase the risk of risk and is associated with certain serious complications resulting in death [1] [2] . Screening for osteoporosis is currently recommended for all women aged 65 years or older. The complications of this disease may be prevented by early detection [1] . Diagnosis is currently based mainly on measurement of bone mineral density (BMD) using dual-energy X-ray absorptiometry (DXA) considered as “gold standard”, but bone densitometry is expensive and has limited availability for routine use in population screening [3] [4] .

Several studies have demonstrated that individuals with osteoporosis have altered morphology of the mandible and there is evidence that intraoral and panoramic radiographic findings might be indicators of osteoporosis [1] [5] -[10] . The BMD of the mandible is correlated with that in the lumbar spine and femoral neck in osteoporosis [7] .

Panoramic radiography (PR) is highlighted in this context because it is frequently conducted in the dental office; it is rapid and low-cost and uses low-dose X-radiation. The most commonly studied measure of mandibular morphology in relation to osteoporosis is the integrity of the inferior border. Erosions of the inferior border are typically scored using the mandibular cortical index (MCI) [9] . This index (known as the Klemetti index) [11] developed for PR in a study in which the authors concluded that individuals with osteoporosis is also more likely to show erosions. They suggested that this evaluation is useful for identifying postmenopausal women with undetected low skeletal BMD or osteoporosis.

One of the most common bone patterns is trabecular bone. Trabecular bone has a branching pattern, as seen in this vertebral specimen. The idea that the fractal index or fractal dimension (FD) of trabecular bone might be related to bone strength is an appealing one, since the fractal index is simple to calculate from clinical radiographies [12] [13] . Thus, the aim of this study was to assess the correlation between different quality analysis parameters of trabecular pattern in digital panoramic radiographies and relations with forearm bone density performed by pDXA in Brazilian population.

2. Material and Methods

Initially the sample consist in 78 panoramic radiographs (PR), but only 68 of these patients had the BMD tests, by pDXA, for this reason our sample was limited to 68 patients, 9 males and 59 females (19 - 73 years old). The digital panoramic radiographs had been gotten in a Veraviewepocs device of J. Morita MFG. CORP.

2.1. Trabecular Pattern Analysis

The trabecular bone pattern study was done using the morphological analysis methodology proposed by White et al. (2000) [14] . For this analysis we used the NIH ImageJ (National Institutes of Health, Bethesda, MD) software. Thus, was performing measurements of the morphological factors of trabecular architecture in the digital panoramic radiographs.

First of all was selected the region of interest (ROI). It was selected the rectangular selection tool, a fixed manner, in the dimension of 230 × 130 pixels located in the region of the mandibular body below the lower canine tooth apex, always on the right, avoiding overlapping anatomical structures, as the mental foramen, Genis apophyses and regions of masticatory stress (premolars and molars) [15] . Thus, it was cut the same size rectangular pictures in all PR [16] (Figure 1).

In Image J software (National Institutes of Health, Bethesda, MD) to analyze the images and possibilities for changes in trabecular pattern of the ROI. The following steps were followed [14] : 1) ROI image duplication; 2) Image blurring (From Gaussian) with 33 radius (pixels)―this step remove all the fine and medium scale structures and retains only large variations in density (low-pass filter); 3) Image subtraction―The blurred image was subtracted from the original of the same patient; 4) Adding―This operation adding a constant to each pixel of the result subtracted image. This generates an image with an average value of about 128, regardless of the original image intensity, discarding gross variations in the intensity RP; 5) Binary transformation―The previous image was transformed into binary image that is black and white; 6) Erode―The previous result was “eroded”. Each pixel is placed with the minimum value of your neighborhood; 7) Dilate―Previous result was “dilated”.

Figure 1. The region of interest (ROI) selected the images of digital panoramic radiographs, in the region of the mandibular body, before the mental foramen and after the mental area, usually in the apical region of the tooth 33/43.

Each pixel is loaded with the maximum value in the vicinity 3 × 3; 8) Inversion―The image of the previous result was reversed. Similar to a photographic negative; 9) Skeletonized―It was removed edges of the pixels of the binary-object-dilated-eroded-inverted images until they were ducted to a single wide skeleton of a single pixel. The object (trabecular) is black and background is white (Table 1).

The skeletonized images from each patient were then saved in TIFF format. Using the same software Image J performed the analysis of the images as follows: Analyze particles-Analysis of the trabecular area checking the average size of trabeculae (particles), TMP, and the number of trabeculae different sizes (particles), DTP. These values represent the connectivity of the trabecular meshwork. The higher the value of TMP and the smaller the value of the DTP, greater trabecular bone connectivity of ROI. To facilitate the analysis of those parameters did a ratio using the following equation:

The second analysis in the same image was by histogram. Thus, evaluating the total number of “bright” pixels (black) in the binary image of the ROI in terms of percentage, called ET.

The third analysis was Fractal Box Count. This analysis provides us with the value of “D” called FD or fractal dimension. Boxes were used in the following sizes 2, 3, 4, 6, 8, 12, 16, 32, 64.

2.2. Bone Mineral Density Measure (BMD)

The peripheral bone densitometry measurements of absorption by dual energy X-ray absorptiometry (DXA) were carried out in the apparatus of Norland, pDXA determining the bone mineral density of the forearm to identify who has a high risk of osteoporosis. The region of choice for scanning and harvesting of data for risk assessment osteoporosis was the distal forearm [17] . The scanning resolution was 1.00 × 1.00 mm , prior to scanning. The radiation dose was less than 0.03 mSv for each examination.

2.3. Statistical Analysis

The statistical analysis was performed by the GMC program, version 2002, Biological Research. In addition to verifying the normality tests were performed Pearson.

3. Results

The sample was 78 panoramic radiographs, but only 68 of these patients had the BMD tests, by pDXA, for this reason our sample was limited to 68 patients, 9 males and 59 females. It should be noted that WHO has no reference data for densitometry measurements for man even so they were considered.

BMD measurements of the total sample based on T-score value for the three sites and the average of these

Table 1. Methods of digital image processing, for skeletal samples from each patient.

values for each patient. It may be noted that only 26 patients have correlation between the three values: Proximal radio = osteopenia, proximal radius plus ulna = osteopenia and distal of radius plus ulna = osteopenia, or 38.2% of the sample (Table 2). HORNER et al. 2002 [10] considered that any measure with standard deviation below or equal to −1 classify the patient in the group of patients with low bone mass, or possessed of some risk for osteoporosis, and consequently the fracture risk.

The correlation was made between the tested parameters, performing the Pearson correlation tests to verify the correlation coefficient between the variables associated and the statistical significance of the possible relationship (Table 3). For the FD values, higher values, the greater the complexity of the structure analyzed in this work, the mandibular ROI [18] .

In Table 4 displays have skeletonized ROI samples for each patient and the respective FD, ET and TMP/DTP values. Each patient has a representative color classification by DXA, according to HORNER et al. (2002) [10] , who classified the patients in healthy (green) and with osteoporosis possibility (yellow and red).

Table 2. Classification of patients according to the measurements of bone mineral density (BMD) by pDXA.

Color significance according to pDXA graphics exam: green = healthy; yellow = medium risk to osteoporosis or osteopenia and red = high risk to osteoporosis.

Table 3. Summary results of correlation tests.

FD: fractal dimension; ET: percentage of black pixels (trabeculae); TMP/DTP: average size ratio and different particle sizes (connectivity); BMD: bone mineral density average.

It was also performed the Cohen’s kappa test to evaluate the correlation between the sample data (trabecular morphology, FD-ET-TMP/DTP and BMD. This analysis was based on the visual observation of the ROI and the mean values of the parameters analyzed. We consider the following parameters: FD lower than 1.4500―unsatis- factory; ET less than 8.00%―unsatisfactory; TMP/DTP lower than 150―not satisfactory. It should be emphasized that using the mean values of the ratio of TMP/DTP, does not reach 500, considered approximately 1/3 of this value, so there is no risk of overestimation, since it represents the connectivity samples. Thus, this parameter had significant weight on visual observation of the ROI.

4. Discussion

Bone quality may be underestimated when based on trabecular bone density due to increase of radiographic density and which has no correlation with the state of bone mineral skeleton. The contributions of soft tissue as much as kVp, exposure time, mA and a different image receptors affect the radiographic outcome. The method used to evaluate the cortex morphology are less susceptible to these types of errors [19] , as recommended morpho-digital studies of trabecular bone pattern [9] [14] [20] -[28] . Therefore, it may be more reliable, perform the evaluation considering the appearance of the cortex in the panoramic radiography, examining the trabecular bone or alveolar process, when, for example, diagnostic analysis of pre-implant bone quality, is required.

White & Rudolf (1999) found sufficient evidence that the trabecular pattern is affected in individuals with osteoporosis, either medium or high risk [29] . We found a significant correlation between the FD factors, ET and TMP/DTP (or connectivity) level of 5% in the Pearson correlation test. Despite the highly significant correlation of FD with ET, however correlation tests showed no significant correlation of BMD, with the analysis parameters of the trabecular morphology.

WHO defines osteoporosis as “a disease characterized by low bone mass and deterioration of bone microarchitecture, leading to increased bone fragility and a consequent increase in fracture risk”. Thus, evaluating the connectivity of trabecular bone is of paramount importance to the dentist, the analysis of the mandibular bone quality. Observation of ROI numbers: 16, 17, 35, 41, 42, 43, 45, 52, 56, 60, 63, 65, 69, 70, 71, 72 e 73, clearly demonstrates a lack of connectivity in the region, therefore, values TMP/DTP very low (less than 100), 72.4% correspond to patients with low bone mass. The no statistical correlation between BMD and TMP/DTP may be due to the samples being not normal or need to improve the size of sample.

SLEMENDA et al. (1996) also had recognized to have modest correspondence enter the values of BMD, in the different skeletal small sites. When consider in set the factor age, the results had been still less satisfactory due to great overlapping between sick people and not sick [30] .

The morphological filter used in this study was similar to that used by Kumasaka (1997) for extracting skeletal pattern of trabecular bone in panoramic digital images, and transformed into binary images with enhanced geometric components [9] [14] [20] . The authors found that the morphology of trabecular pattern in anterior mandible based on radiographic studies, was indistinct, and the digital extraction site by morphological filter has a very clear observation and revealed details before scarce. For this reason the election site to the digital

Table 4. Skeletonized areas of trabecular structures (patients: 1-78).

FD: fractal dimension; % E.T.: percentage of black pixelss; TMP/DTP: proportion between average size of trabeculae (particles) and number of trabeculae different sizes (particles). Green = healthy, yellow = ostepenia, red = osteoporosis, by pDEXA, according to HORNER et al. (2002).

processing of panoramic radiographs images, and used in this study. The ROI chosen because it is free of image artifacts showed in processed images much more details. Thus, the results suggest that images “improved” radiographs with “poor” quality, can show details of trabecular skeletal pattern. It is necessary to define the best binary processing threshold for adequate results. In this study was followed the technique which can be further enhanced. This technique was possible to see the direction, continuation or termination of trabeculae and are extremely useful for the analysis of bone structures, computer-assisted, to determine, for example, the prognosis following surgery, or to assess certain medical treatment [20] [31] .

5. Conclusion

In conclusion, due to easy access to the panoramic radiograph by a dentist, low cost exam, and frequency, this exam has been recommended as routine in dental care for initial evaluation; the professional should conduct visual examinations of the cortex and the standard morpho-digital of their patients, relating them to factors such as body bone mass index, age and sex, for in case of finding predictors of osteoporosis, and discuss with other medical professionals, the patient’s indication for examination of bone densitometry by DXA. The no statistical correlation between BMD and fractal dimension and TMP/DTP may be due to the samples being not normal or need to improve the size of sample.


*Corresponding author.

Conflicts of Interest

The authors declare no conflicts of interest.


[1] Erdogan, O., Incki, K.K., Benlidayi, M.E., Seydaoglu, G. and Kelekci, S. (2009) Dental and Radiographic Findings as Predictors of Osteoporosis in Postmenopausal Women. Geriatrics & Gerontology International, 9, 155-164.
[2] Hua, Y., Nackaerts, O., Duyck, J., Maes, F. and Jacobs, R. (2009) Bone Quality Assessment Based on Cone Beam Computed Tomography Imaging. Clinical Oral Implants Research, 20, 767-771.
[3] Reginster, J.Y. and Burlet, N. (2006) Osteoporosis: A Still Increasing Prevalence. Bone, 38, S4-S9.
[4] Nakamoto, T., Taguchi, A., Ohtsuka, M., Suei, Y., Fujita, M., Tanimoto, K., Tsuda, M., Sanada, M., Ohama, K., Takahashi, J. and Rohlin, M. (2003) Dental Panoramic Radiograph as a Tool to Detect Postmenopausal Women with Low Bone Mineral Density: Untrained General Dental Practitioners’ Diagnostic Performance. Osteoporosis International, 14, 659-664.
[5] Alonso, M.B.C.C., Cortes, A.R.G., Camargo, A.J., Arita, E.S., Haiter-Neto, F. and Watanabe, P.C.A. (2011) Assessment of Panoramic Radiomorphometric Indices of the Mandible in a Brazilian Population. ISRN Rheumatology, 2011, Article ID: 854287.
[6] Camargo, A.J., Arita, E.S., Fenandez, M.C.C. and Watanabe, P.C.A. (2015) Comparison of Two Radiological Methods for Evaluation of Bone Density in Postmenopausal Women. International Journal of Morphology, 33, 732-736.
[7] Koh, K.J. and Kim, K.A. (2011) Utility of the Computed Tomography Indices on Cone Beam Computed Tomography Images in the Diagnosis of Osteoporosis in Women. Imaging Science in Dentistry, 41, 101-106.
[8] Dagistan, S. and Bilge, O.M. (2010) Comparison of Antegonial Index, Mental Index, Panoramic Mandibular Index and Mandibular Cortical Index Values in the Panoramic Radiographs of Normal Males and Male Patients with Osteoporosis. Dentomaxillofacial Radiology, 39, 290-294.
[9] White, S.C. (2002) Oral Radiographic Predictors of Osteoporosis. Dentomaxillofacial Radiology, 31, 84-92.
[10] Horner, K., Devlin, H. and Harvey, L. (2002) Detecting Patients with Low Skeletal Bone Mass. Journal of Dentistry, 30, 171-175.
[11] Klemetti, E., Kolmakov, S. and Kroger, H. (1994) Pantomography in Assessment of the Osteoporosis Risk Group. Scandinavian Journal of Dental Research, 102, 68-72.
[12] Sánchez, I. and Uzcátegui, G. (2011) Fractals in Dentistry. Journal of Dentistry, 394, 273-292.
[13] Oliveira, M.L., Pedrosa, E.F.N.C., Cruz, A.D., Haiter-Neto, F., Paula, F.J.A. and Watanabe, P.C.A. (2013) Relationship between Bone Mineral Density and Trabecular Bone Pattern in Postmenopausal Osteoporotic Brazilian Women. Clinical Oral Investigations, 17, 1847-1853.
[14] White, S.C., Cohen, J.M. and Mourshed, F.A. (2000) Digital Analysis of Trabecular Pattern in Jaws of Patients with Sickle Cell Anemia. Dentomaxillofacial Radiology, 29, 119-124.
[15] Klemetti, E. (1996) A Review of Residual Ridge Resorption and Bone Density. Journal of Prosthetic Dentistry, 75, 512-514.
[16] Shrout, M.K., Hildebolt, C.F. and Potter, B.J. (1997) The Effect of Varying the Region of Interest on Calculations of Fractal Index. Dentomaxillofacial Radiology, 26, 295-298.
[17] Notelovitz, M. (2008) Osteoporosis: Prevention, Diagnosis, and Management. Professional Communications.
[18] Nascimento, M.P.B. (2001) Fractal Dimension. Dissertation. Maths and Statistics Institute, Rio de Janeiro.
[19] Klemetti, E., Kolmakov, S., Heiskanen, P., Vainio, P. and Lassila, V. (1993) Panoramic Mandibular Index and Bone Mineral Densities in Postmenopausal Women. Oral Surgery, Oral Medicine, Oral Pathology, Oral Radiology, and Endodontology, 75, 774-779.
[20] Kumasaka, S. and Kashima, I. (1997) Initial Investigation of Mathematical Morphology for the Digital Extraction of the Skeletal Characteristics of Trabecular Bone. Dentomaxillofacial Radiology, 26, 161-168.
[21] Link, T.M., Majumdar, S., Grampp, S., Guglielmi, G., Van Kuijk, C., Imhof, H. and Adams, J.E. (1999) Imaging of Trabecular Bone Structure in Osteoporosis. European Radiology, 9, 1781-1788.
[22] Geraets, W.G. and Van Der Stelt, P.F. (2000) Fractal Properties of Bone. Dentomaxillofacial Radiology, 29, 144-153.
[23] Bollen, A.M., Taguchi, A., Hujoel, P.P. and Hollender, L.G. (2001) Fractal Dimension on Dental Radiographs. Dentomaxillofacial Radiology, 30, 270-275.
[24] Feltrin, G.P., Macchi, V., Saccavini, C., Tosi, E., Dus, C., Fassina, A., Parenti, A. and De Caro, R. (2001) Fractal Analysis of Lumbar Vertebral Cancellous Bone Architecture. Clinical Anatomy, 14, 414-417.
[25] Gurdal, P., Hildebolt, C.F. and Akdeniz, B.G. (2001) The Effects of Different Image File Formats and Image-Analysis Software Programs on Dental Radiometric Digital Evaluations. Dentomaxillofacial Radiology, 30, 50-55.
[26] Jonasson, G., Bankvall, G. and Kiliaridis, S. (2001) Estimation of Skeletal Bone Mineral Density by Means of the Trabecular Pattern of the Alveolar Bone, Its Interdental Thickness, and the Bone Mass of the Mandible. Oral Surgery, Oral Medicine, Oral Pathology, Oral Radiology, and Endodontology, 92, 346-352.
[27] Wojtowicz, A., Chaberek, S., Kryst, L., Urbanowska, E., Ciechowicz, K. and Ostrowski, K. (2003) Fourier and Fractal Analysis of Maxillary Alveolar Ridge Repair Using Platelet Rich Plasma (PRP) and Inorganic Bovine Bone. International Journal of Oral & Maxillofacial Surgery, 321, 84-86.
[28] Leichter, I.S., Simkin, A., Neeman, V., Jabschinsky, C., Schoenfeld, D., Foldes, A.J. and Liebergall, M. (2001) Optical Processing of Radiographic Trabecular Pattern versus Bone Mineral Density of Proximal Femur as Measures of Bone Strength. Journal of Clinical Densitometry, 4, 121-129.
[29] White, S.C. and Rudolph, D.J. (1999) Alterations of the Trabecular Pattern of the Jaws in Patients with Osteoporosis. Oral Surgery, Oral Medicine, Oral Pathology, Oral Radiology, and Endodontology, 88, 628-635.
[30] Slemenda, C., Longcope, C., Peacock, M., Hui, S. and Johnston, C.C. (1996) Sex Steroids, Bone Mass, and Bone Loss. A Prospective Study of Pre-, Peri-, and Postmenopausal Women. Journal of Clinical Investigation, 97, 14.
[31] Oka, K., Kumasaka, S. and Kashima, I. (2002) Assessment of Bone Feature Parameters from Lumbar Trabecular Skeletal Using Mathematical Morphology Image Processing. Journal of Bone and Mineral Metabolism, 20, 201-208.

comments powered by Disqus

Copyright © 2020 by authors and Scientific Research Publishing Inc.

Creative Commons License

This work and the related PDF file are licensed under a Creative Commons Attribution 4.0 International License.