A novel registration-based methodology for prediction of trabecular bone fabric from clinical QCT: A comprehensive analysis

Osteoporosis leads to hip fractures in aging populations and is diagnosed by modern medical imaging techniques such as quantitative computed tomography (QCT). Hip fracture sites involve trabecular bone, whose strength is determined by volume fraction and orientation, known as fabric. However, bone fabric cannot be reliably assessed in clinical QCT images of proximal femur. Accordingly, we propose a novel registration-based estimation of bone fabric designed to preserve tensor properties of bone fabric and to map bone fabric by a global and local decomposition of the gradient of a non-rigid image registration transformation. Furthermore, no comprehensive analysis on the critical components of this methodology has been previously conducted. Hence, the aim of this work was to identify the best registration-based strategy to assign bone fabric to the QCT image of a patient’s proximal femur. The normalized correlation coefficient and curvature-based regularization were used for image-based registration and the Frobenius norm of the stretch tensor of the local gradient was selected to quantify the distance among the proximal femora in the population. Based on this distance, closest, farthest and mean femora with a distinction of sex were chosen as alternative atlases to evaluate their influence on bone fabric prediction. Second, we analyzed different tensor mapping schemes for bone fabric prediction: identity, rotation-only, rotation and stretch tensor. Third, we investigated the use of a population average fabric atlas. A leave one out (LOO) evaluation study was performed with a dual QCT and HR-pQCT database of 36 pairs of human femora. The quality of the fabric prediction was assessed with three metrics, the tensor norm (TN) error, the degree of anisotropy (DA) error and the angular deviation of the principal tensor direction (PTD). The closest femur atlas (CTP) with a full rotation (CR) for fabric mapping delivered the best results with a TN error of 7.3 ± 0.9%, a DA error of 6.6 ± 1.3% and a PTD error of 25 ± 2°. The closest to the population mean femur atlas (MTP) using the same mapping scheme yielded only slightly higher errors than CTP for substantially less computing efforts. The population average fabric atlas yielded substantially higher errors than the MTP with the CR mapping scheme. Accounting for sex did not bring any significant improvements. The identified fabric mapping methodology will be exploited in patient-specific QCT-based finite element analysis of the proximal femur to improve the prediction of hip fracture risk.

[1]  Paul M. Thompson,et al.  Generalized Tensor-Based Morphometry of HIV/AIDS Using Multivariate Statistics on Deformation Tensors , 2008, IEEE Transactions on Medical Imaging.

[2]  Mauricio Reyes,et al.  Image-based vs. mesh-based statistical appearance models of the human femur: implications for finite element simulations. , 2014, Medical engineering & physics.

[3]  Daniel Rueckert,et al.  Nonrigid registration using free-form deformations: application to breast MR images , 1999, IEEE Transactions on Medical Imaging.

[4]  E. Dall’Ara,et al.  Orthotropic HR-pQCT-based FE models improve strength predictions for stance but not for side-way fall loading compared to isotropic QCT-based FE models of human femurs. , 2014, Journal of the mechanical behavior of biomedical materials.

[5]  Robert C. Bolles,et al.  Random sample consensus: a paradigm for model fitting with applications to image analysis and automated cartography , 1981, CACM.

[6]  Xavier Pennec,et al.  Simultaneous Multiscale Polyaffine Registration by Incorporating Deformation Statistics , 2012, MICCAI.

[7]  Guido Gerig,et al.  Unbiased diffeomorphic atlas construction for computational anatomy , 2004, NeuroImage.

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

[9]  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.

[10]  Mauricio Reyes,et al.  Prediction of Trabecular Bone Anisotropy from Quantitative Computed Tomography Using Supervised Learning and a Novel Morphometric Feature Descriptor , 2015, MICCAI.

[11]  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.

[12]  Keita Ito,et al.  A novel approach to estimate trabecular bone anisotropy from stress tensors , 2015, Biomechanics and modeling in mechanobiology.

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

[14]  Hervé Delingette,et al.  Towards a Statistical Atlas of Cardiac Fiber Structure , 2006, MICCAI.

[15]  William H. Press,et al.  Numerical Recipes 3rd Edition: The Art of Scientific Computing , 2007 .

[16]  Mauricio Reyes,et al.  Statistical analysis of the inter-individual variations of the bone shape, volume fraction and fabric and their correlations in the proximal femur. , 2017, Bone.

[17]  Maxime Sermesant,et al.  Spatio-Temporal Tensor Decomposition of a Polyaffine Motion Model for a Better Analysis of Pathological Left Ventricular Dynamics , 2015, IEEE Transactions on Medical Imaging.

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

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

[20]  Alejandro F. Frangi,et al.  Patient-Specific Biomechanical Modeling of Bone Strength Using Statistically-Derived Fabric Tensors , 2015, Annals of Biomedical Engineering.

[21]  Mauricio Reyes,et al.  Biomechanical Role of Bone Anisotropy Estimated on Clinical CT Scans by Image Registration , 2016, Annals of Biomedical Engineering.

[22]  M. Fréchet.  Les éléments aléatoires de nature quelconque dans un espace distancié , 1948 .

[23]  Dinggang Shen,et al.  Multi-contrast diffusion tensor image registration with structural MRI , 2012, 2012 9th IEEE International Symposium on Biomedical Imaging (ISBI).

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

[25]  Hrvoje Kalinić,et al.  Atlas-based image segmentation: A Survey , 2009 .

[26]  L. R. Dice.  Measures of the Amount of Ecologic Association Between Species , 1945 .

[27]  Xavier Pennec,et al.  Geometry-Aware Multiscale Image Registration via OBBTree-Based Polyaffine Log-Demons , 2011, MICCAI.

[28]  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.

[29]  Hervé Delingette,et al.  A Computational Framework for the Statistical Analysis of Cardiac Diffusion Tensors: Application to a Small Database of Canine Hearts , 2007, IEEE Transactions on Medical Imaging.

[30]  Kemal Tuncali,et al.  Multimodality non-rigid image registration for planning, targeting and monitoring during CT-guided percutaneous liver tumor cryoablation. , 2010, Academic radiology.

[31]  V. Bousson,et al.  Hip Fracture Risk and Proximal Femur Geometry from DXA Scans , 2002, Osteoporosis International.

[32]  P. Zysset,et al.  Assessment of volume fraction and fabric in the distal radius using HR-pQCT. , 2009, Bone.