Adaptive Segmentation of Knee Radiographs for Selecting the Optimal ROI in Texture Analysis

The purposes of this study were to investigate: 1) the effect of placement of region-of-interest (ROI) for texture analysis of subchondral bone in knee radiographs, and 2) the ability of several texture descriptors to distinguish between the knees with and without radiographic osteoarthritis (OA). Bilateral posterior-anterior knee radiographs were analyzed from the baseline of OAI and MOST datasets. A fully automatic method to locate the most informative region from subchondral bone using adaptive segmentation was developed. We used an oversegmentation strategy for partitioning knee images into the compact regions that follow natural texture boundaries. LBP, Fractal Dimension (FD), Haralick features, Shannon entropy, and HOG methods were computed within the standard ROI and within the proposed adaptive ROIs. Subsequently, we built logistic regression models to identify and compare the performances of each texture descriptor and each ROI placement method using 5-fold cross validation setting. Importantly, we also investigated the generalizability of our approach by training the models on OAI and testing them on MOST dataset.We used area under the receiver operating characteristic (ROC) curve (AUC) and average precision (AP) obtained from the precision-recall (PR) curve to compare the results. We found that the adaptive ROI improves the classification performance (OA vs. non-OA) over the commonly used standard ROI (up to 9% percent increase in AUC). We also observed that, from all texture parameters, LBP yielded the best performance in all settings with the best AUC of 0.840 [0.825, 0.852] and associated AP of 0.804 [0.786, 0.820]. Compared to the current state-of-the-art approaches, our results suggest that the proposed adaptive ROI approach in texture analysis of subchondral bone can increase the diagnostic performance for detecting the presence of radiographic OA

ToScA North America (6 – 8 June 2017, The University of Texas, Austin, TX) Program

ToScA North America will address key areas of science, including Multi-modal Imaging, Geosciences, Forensics, Increasing Contrast, Educational Outreach, Data, Materials Science and Medical and Biological Science.University of Texas High-Resolution X-ray CT Facility (UTCT); Jackson School of Geosciences, The University of Texas at Austin; Natural History Museum (London); Royal Microscopical Society (Oxford, UK)Geological Science

A new straightforward method for semi-automated segmentation of trabecular bone from cortical bone in diverse and challenging morphologies

Many physiological, biomechanical, evolutionary and clinical studies that explore skeletal structure and function require successful separation of trabecular from cortical compartments of a bone that has been imaged by X-ray micro-computed tomography (micro-CT) prior to analysis. Separation often involves manual subdivision of these two similarly radio-opaque compartments, which can be time-consuming and subjective. We have developed an objective, semi-automated protocol which reduces user bias and enables straightforward, user-friendly segmentation of trabecular from the cortical bone without requiring sophisticated programming expertise. This method can conveniently be used as a ‘recipe’ in commercial programmes (Avizo herein) and applied to a variety of datasets. Here, we characterize and share this recipe, and demonstrate its application to a range of murine and human bone types, including normal and osteoarthritic specimens, and bones with distinct embryonic origins and spanning a range of ages. We validate the method by testing inter-user bias during the scan preparation steps and confirm utility in the architecturally challenging analysis of growing murine epiphyses. We also report details of the recipe, so that other groups can readily re-create a similar method in open access programmes. Our aim is that this method will be adopted widely to create a reproducible and time-efficient method of segmenting trabecular and cortical bone

Automated 3D trabecular bone structure analysis of the proximal femur—prediction of biomechanical strength by CT and DXA

The standard diagnostic technique for assessing osteoporosis is dual X-ray absorptiometry (DXA) measuring bone mass parameters. In this study, a combination of DXA and trabecular structure parameters (acquired by computed tomography [CT]) most accurately predicted the biomechanical strength of the proximal femur and allowed for a better prediction than DXA alone. An automated 3D segmentation algorithm was applied to determine specific structure parameters of the trabecular bone in CT images of the proximal femur. This was done to evaluate the ability of these parameters for predicting biomechanical femoral bone strength in comparison with bone mineral content (BMC) and bone mineral density (BMD) acquired by DXA as standard diagnostic technique. One hundred eighty-seven proximal femur specimens were harvested from formalin-fixed human cadavers. BMC and BMD were determined by DXA. Structure parameters of the trabecular bone (i.e., morphometry, fuzzy logic, Minkowski functionals, and the scaling index method [SIM]) were computed from CT images. Absolute femoral bone strength was assessed with a biomechanical side-impact test measuring failure load (FL). Adjusted FL parameters for appraisal of relative bone strength were calculated by dividing FL by influencing variables such as body height, weight, or femoral head diameter. The best single parameter predicting FL and adjusted FL parameters was apparent trabecular separation (morphometry) or DXA-derived BMC or BMD with correlations up to r = 0.802. In combination with DXA, structure parameters (most notably the SIM and morphometry) added in linear regression models significant information in predicting FL and all adjusted FL parameters (up to R adj = 0.872) and allowed for a significant better prediction than DXA alone. A combination of bone mass (DXA) and structure parameters of the trabecular bone (linear and nonlinear, global and local) most accurately predicted absolute and relative femoral bone strength

Quantifying the Tibiofemoral Joint Space Using X-ray Tomosynthesis

Purpose: Digital x-ray tomosynthesis (DTS) has the potential to provide 3D information about the knee joint in a load-bearing posture, which may improve diagnosis and monitoring of knee osteoarthritis compared with projection radiography, the current standard of care. Manually quantifying and visualizing the joint space width (JSW) from 3D tomosynthesis datasets may be challenging. This work developed a semiautomated algorithm for quantifying the 3D tibiofemoral JSW from reconstructed DTS images. The algorithm was validated through anthropomorphic phantom experiments and applied to three clinical datasets. Methods: A user-selected volume of interest within the reconstructed DTS volume was enhanced with 1D multiscale gradient kernels. The edge-enhanced volumes were divided by polarity into tibial and femoral edge maps and combined across kernel scales. A 2D connected components algorithm was performed to determine candidate tibial and femoral edges. A 2D joint space width map (JSW) was constructed to represent the 3D tibiofemoral joint space. To quantify the algorithm accuracy, an adjustable knee phantom was constructed, and eleven posterior–anterior (PA) and lateral DTS scans were acquired with the medial minimum JSW of the phantom set to 0–5 mm in 0.5 mm increments (VolumeRadTM, GE Healthcare, Chalfont St. Giles, United Kingdom). The accuracy of the algorithm was quantified by comparing the minimum JSW in a region of interest in the medial compartment of the JSW map to the measured phantom setting for each trial. In addition, the algorithm was applied to DTS scans of a static knee phantom and the JSW map compared to values estimated from a manually segmented computed tomography (CT) dataset. The algorithm was also applied to three clinical DTS datasets of osteoarthritic patients. Results: The algorithm segmented the JSW and generated a JSW map for all phantom and clinical datasets. For the adjustable phantom, the estimated minimum JSW values were plotted against the measured values for all trials. A linear fit estimated a slope of 0.887 (R2¼0.962) and a mean error across all trials of 0.34 mm for the PA phantom data. The estimated minimum JSW values for the lateral adjustable phantom acquisitions were found to have low correlation to the measured values (R2¼0.377), with a mean error of 2.13 mm. The error in the lateral adjustable-phantom datasets appeared to be caused by artifacts due to unrealistic features in the phantom bones. JSW maps generated by DTS and CT varied by a mean of 0.6 mm and 0.8 mm across the knee joint, for PA and lateral scans. The tibial and femoral edges were successfully segmented and JSW maps determined for PA and lateral clinical DTS datasets. Conclusions: A semiautomated method is presented for quantifying the 3D joint space in a 2D JSW map using tomosynthesis images. The proposed algorithm quantified the JSW across the knee joint to sub-millimeter accuracy for PA tomosynthesis acquisitions. Overall, the results suggest that x-ray tomosynthesis may be beneficial for diagnosing and monitoring disease progression or treatment of osteoarthritis by providing quantitative images of JSW in the load-bearing knee