Diagnosis of osteoporosis by extraction of trabecular features from hip radiographs using support vector machine: An investigation panorama with DXA

BACKGROUND Lifespan and its quality can be improved by early diagnosis of osteoporosis. Analysis of trabecular boundness on digital hip radiographs could be useful for identifying subjects with low bone mineral density (BMD) or osteoporosis. The main aim of our study was to evaluate the ability of a kernel-based support vector machine (SVM) with respect to diagnosis and add to knowledge about the trabecular features of digital hip radiographs for identifying subjects with low BMD. METHOD In this paper we present an SVM kernel classifier-based computer-aided diagnosis (CAD) system for osteoporotic risk detection using digital hip radiographs. Initially, the original radiograph was intensified, then trabecular features such as boundness, orientation, solidity of spur and delta were evaluated and radial bias function (RBF) based discrimination was manifested. The next step was the evaluation of the diagnostic capability of the proposed method in order to spot subjects with low BMD at the femoral neck in 50 (50.7 ± 14.3 years) South Indian women with no previous history of osteoporotic fracture. Out of 50 subjects, 28 were used to train the classifier and the other 22 were used for testing. RESULTS The proposed system has achieved the highest classification accuracy documented so far by means of a fivefold cross-validation analysis with mean accuracy of 90% (95% confidence interval (CI): 82 to 98%); sensitivity and positive predictive value (PPV) were 90% (95% CI: 82 to 98%) and 89% (95% CI: 81 to 97%), respectively. Pearson's correlation was observed at the level of p<0.001, between extracted image trabecular features with age and BMDs measured by dual energy x-ray absorptiometry (DXA). Extracted image features also demonstrated significant differences between high and low BMD groups at the level of p<0.001. CONCLUSION Our findings suggest that the proposed CAD system with SVM would be useful for spotting women vulnerable to osteoporotic risk.

[1]  H. Kopera,et al.  Briefe an die Redaktion , 1989, Klinische Wochenschrift.

[2]  Michel Couprie,et al.  Isthmus-Based 6-Directional Parallel Thinning Algorithms , 2011, DGCI.

[3]  N. Balakrishna,et al.  Bone status of Indian women from a low-income group and its relationship to the nutritional status , 2005, Osteoporosis International.

[4]  Sankar K. Pal,et al.  Automatic Exact Histogram Specification for Contrast Enhancement and Visual System Based Quantitative Evaluation , 2011, IEEE Transactions on Image Processing.

[5]  T. Jämsä,et al.  Experimental hip fracture load can be predicted from plain radiography by combined analysis of trabecular bone structure and bone geometry , 2008, Osteoporosis International.

[6]  R F Kilcoyne,et al.  Ability of four different techniques of measuring bone mass to diagnose vertebral fractures in postmenopausal women , 1987, Journal of bone and mineral research : the official journal of the American Society for Bone and Mineral Research.

[7]  H. Kopera [Prophylaxis and treatment of osteoporosis]. , 1989, Klinische Wochenschrift.

[8]  Vladimir N. Vapnik,et al.  The Nature of Statistical Learning Theory , 2000, Statistics for Engineering and Information Science.

[9]  Bernhard Schölkopf,et al.  Comparing support vector machines with Gaussian kernels to radial basis function classifiers , 1997, IEEE Trans. Signal Process..

[10]  H. Weinans,et al.  Mechanical Consequences of Bone Loss in Cancellous Bone , 2001, Journal of bone and mineral research : the official journal of the American Society for Bone and Mineral Research.

[11]  R. Rizzoli,et al.  Bone strength and its determinants , 2003, Osteoporosis International.

[12]  B. Ripley,et al.  Pattern Recognition , 1968, Nature.

[13]  Isabelle Guyon,et al.  An Introduction to Variable and Feature Selection , 2003, J. Mach. Learn. Res..

[14]  J. Kanis,et al.  Assessment of fracture risk and its application to screening for postmenopausal osteoporosis: Synopsis of a WHO report , 1994, Osteoporosis International.

[15]  Maria-Grazia Ascenzi,et al.  Variation of trabecular architecture in proximal femur of postmenopausal women. , 2011, Journal of biomechanics.

[16]  B. Minasny The Elements of Statistical Learning, Second Edition, Trevor Hastie, Robert Tishirani, Jerome Friedman. (2009), Springer Series in Statistics, ISBN 0172-7397, 745 pp , 2009 .

[17]  S. Cummings,et al.  BMD at Multiple Sites and Risk of Fracture of Multiple Types: Long‐Term Results From the Study of Osteoporotic Fractures , 2003, Journal of bone and mineral research : the official journal of the American Society for Bone and Mineral Research.

[18]  Vladimir Vapnik,et al.  Statistical learning theory , 1998 .

[19]  Osman N. Uçan,et al.  Early diagnosis of osteoporosis using Artificial Neural Networks and Support Vector Machines , 2012, 2012 20th Signal Processing and Communications Applications Conference (SIU).

[20]  A. Hofman,et al.  Fracture incidence and association with bone mineral density in elderly men and women: the Rotterdam Study. , 2004, Bone.

[21]  N. Tandon,et al.  Bone health in healthy Indian population aged 50 years and above , 2011, Osteoporosis International.

[22]  E. Barrett-Connor,et al.  Identification and fracture outcomes of undiagnosed low bone mineral density in postmenopausal women: results from the National Osteoporosis Risk Assessment. , 2001, JAMA.

[23]  P. Delmas,et al.  Bone quality--the material and structural basis of bone strength and fragility. , 2006, The New England journal of medicine.

[24]  Takio Kurita,et al.  Diagnosis of osteoporosis from dental panoramic radiographs using the support vector machine method in a computer-aided system , 2012, BMC Medical Imaging.

[25]  Zicheng Guo,et al.  Parallel thinning with two-subiteration algorithms , 1989, Commun. ACM.

[26]  M. Anburajan,et al.  Analysis of Texture Patterns in Diagnosing Osteoporosis Using Proximal Femur X-Ray Images , 2011 .

[27]  C. Cooper,et al.  Guidelines for diagnosis and management of osteoporosis , 2005, Osteoporosis International.

[28]  C. Netelenbos,et al.  A new method for automatic recognition of the radiographic trabecular pattern , 1990, Journal of bone and mineral research : the official journal of the American Society for Bone and Mineral Research.

[29]  Alexander J. Smola,et al.  Learning with kernels , 1998 .

[30]  Latifa Hamami,et al.  An Interconnectivity Index for Osteoporosis Assessment Using X-Ray Images , 2013 .

[31]  Trevor Hastie,et al.  The Elements of Statistical Learning , 2001 .

[32]  B. Snyder,et al.  Radiographic trabecular 2D and 3D parameters of proximal femoral bone cores correlate with each other and with yield stress , 2009, Osteoporosis International.

[33]  M. Singh,et al.  Changes in trabecular pattern of the upper end of the femur as an index of osteoporosis. , 1970, The Journal of bone and joint surgery. American volume.

[34]  Ashok Samal,et al.  A simple method for fitting of bounding rectangle to closed regions , 2007, Pattern Recognit..

[35]  W. Leow,et al.  DETECTION OF FEMUR AND RADIUS FRACTURES IN X-RAY IMAGES , 2004 .

[36]  C. Cooper,et al.  Hip fractures in the elderly: A world-wide projection , 1992, Osteoporosis International.