Combined radiogrammetry and texture analysis for early diagnosis of osteoporosis using Indian and Swiss data

Osteoporosis is a bone disorder characterized by bone loss and decreased bone strength. The most widely used technique for detection of osteoporosis is the measurement of bone mineral density (BMD) using dual energy X-ray absorptiometry (DXA). But DXA scans are expensive and not widely available in low-income economies. In this paper, we propose a low cost pre-screening tool for the detection of low bone mass, using cortical radiogrammetry of third metacarpal bone and trabecular texture analysis of distal radius from hand and wrist radiographs. An automatic segmentation algorithm to automatically locate and segment the third metacarpal bone and distal radius region of interest (ROI) is proposed. Cortical measurements such as combined cortical thickness (CCT), cortical area (CA), percent cortical area (PCA) and Barnett Nordin index (BNI) were taken from the shaft of third metacarpal bone. Texture analysis of trabecular network at the distal radius was performed using features obtained from histogram, gray level Co-occurrence matrix (GLCM) and morphological gradient method (MGM). The significant cortical and texture features were selected using independent sample t-test and used to train classifiers to classify healthy subjects and people with low bone mass. The proposed pre-screening tool was validated on two ethnic groups, Indian sample population and Swiss sample population. Data of 134 subjects from Indian sample population and 65 subjects from Swiss sample population were analysed. The proposed automatic segmentation approach shows a detection accuracy of 86% in detecting the third metacarpal bone shaft and 90% in accurately locating the distal radius ROI. Comparison of the automatic radiogrammetry to the ground truth provided by experts show a mean absolute error of 0.04 mm for cortical width of healthy group, 0.12 mm for cortical width of low bone mass group, 0.22 mm for medullary width of healthy group, and 0.26 mm for medullary width of low bone mass group. Independent sample t-test was used to select the most discriminant features, to be used as input for training the classifiers. Pearson correlation analysis of the extracted features with DXA-BMD of lumbar spine (DXA-LS) shows significantly high correlation values. Classifiers were trained with the most significant features in the Indian and Swiss sample data. Weighted KNN classifier shows the best test accuracy of 78% for Indian sample data and 100% for Swiss sample data. Hence, combined automatic radiogrammetry and texture analysis is shown to be an effective low cost pre-screening tool for early diagnosis of osteoporosis.