Prediction of Trabecular Bone Anisotropy from Quantitative Computed Tomography Using Supervised Learning and a Novel Morphometric Feature Descriptor

Patient-specific biomechanical models including local bone mineral density and anisotropy have gained importance for assessing musculoskeletal disorders. However the trabecular bone anisotropy captured by high-resolution imaging is only available at the peripheral skeleton in clinical practice. In this work, we propose a supervised learning approach to predict trabecular bone anisotropy that builds on a novel set of pose invariant feature descriptors. The statistical relationship between trabecular bone anisotropy and feature descriptors were learned from a database of pairs of high resolution QCT and clinical QCT reconstructions. On a set of leave-one-out experiments, we compared the accuracy of the proposed approach to previous ones, and report a mean prediction error of 6% for the tensor norm, 6% for the degree of anisotropy and 19i¾? for the principal tensor direction. These findings show the potential of the proposed approach to predict trabecular bone anisotropy from clinically available QCT images.

[1]  Willi A Kalender,et al.  An anatomic coordinate system of the femoral neck for highly reproducible BMD measurements using 3D QCT. , 2005, Computerized Medical Imaging and Graphics.

[2]  F. Taddei,et al.  A novel approach to estimate trabecular bone anisotropy using a database approach. , 2013, Journal of biomechanics.

[3]  Steven K Boyd,et al.  Mapping anisotropy of the proximal femur for enhanced image based finite element analysis. , 2014, Journal of biomechanics.

[4]  Raphaël Marée,et al.  Fast Multi-class Image Annotation with Random Subwindows and Multiple Output Randomized Trees , 2009, VISAPP.

[5]  Alejandro F. Frangi,et al.  Statistical estimation of femur micro-architecture using optimal shape and density predictors. , 2015, Journal of biomechanics.

[6]  Dieter H. Pahr,et al.  From high-resolution CT data to finite element models: development of an integrated modular framework , 2009 .

[7]  P. Zysset,et al.  Bone Volume Fraction and Fabric Anisotropy Are Better Determinants of Trabecular Bone Stiffness Than Other Morphological Variables , 2015, Journal of bone and mineral research : the official journal of the American Society for Bone and Mineral Research.

[8]  Raphaël Marée,et al.  Extremely Randomized Trees and Random Subwindows for Image Classification, Annotation, and Retrieval , 2013 .

[9]  Max A. Viergever,et al.  elastix: A Toolbox for Intensity-Based Medical Image Registration , 2010, IEEE Transactions on Medical Imaging.

[10]  Xavier Pennec,et al.  A Riemannian Framework for Tensor Computing , 2005, International Journal of Computer Vision.

[11]  Robert M. Haralick,et al.  Textural Features for Image Classification , 1973, IEEE Trans. Syst. Man Cybern..

[12]  F. Kainberger,et al.  A nonlinear QCT-based finite element model validation study for the human femur tested in two configurations in vitro. , 2013, Bone.

[13]  R. Mann,et al.  Characterization of microstructural anisotropy in orthotropic materials using a second rank tensor , 1984 .

[14]  Yan Kang,et al.  A new accurate and precise 3-D segmentation method for skeletal structures in volumetric CT data , 2003, IEEE Transactions on Medical Imaging.

[15]  Andrés Ortiz,et al.  Segmentation of Brain MRI Using SOM-FCM-Based Method and 3D Statistical Descriptors , 2013, Comput. Math. Methods Medicine.

[16]  Mauricio Reyes,et al.  Including the trabecular anisotropy from registered microCT data in homogenized FE model improves the bone's mechanical predictions , 2014 .

[17]  V. Gudnason,et al.  Assessment of incident spine and hip fractures in women and men using finite element analysis of CT scans , 2014, Journal of bone and mineral research : the official journal of the American Society for Bone and Mineral Research.