Biomechanical Role of Bone Anisotropy Estimated on Clinical CT Scans by Image Registration

Image-based modeling is a popular approach to perform patient-specific biomechanical simulations. Accurate modeling is critical for orthopedic application to evaluate implant design and surgical planning. It has been shown that bone strength can be estimated from the bone mineral density (BMD) and trabecular bone architecture. However, these findings cannot be directly and fully transferred to patient-specific modeling since only BMD can be derived from clinical CT. Therefore, the objective of this study was to propose a method to predict the trabecular bone structure using a µCT atlas and an image registration technique. The approach has been evaluated on femurs and patellae under physiological loading. The displacement and ultimate force for femurs loaded in stance position were predicted with an error of 2.5% and 3.7%, respectively, while predictions obtained with an isotropic material resulted in errors of 7.3% and 6.9%. Similar results were obtained for the patella, where the strain predicted using the registration approach resulted in an improved mean squared error compared to the isotropic model. We conclude that the registration of anisotropic information from of a single template bone enables more accurate patient-specific simulations from clinical image datasets than isotropic model.

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

[2]  E. Itoi,et al.  Patellar morphology and femoral component geometry influence patellofemoral contact stress in total knee arthroplasty without patellar resurfacing , 2012, Knee Surgery, Sports Traumatology, Arthroscopy.

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

[4]  U. Whitcher Boundary Conditions , 2011 .

[5]  P. Zysset,et al.  An anisotropic elastic-viscoplastic damage model for bone tissue , 2012, Biomechanics and Modeling in Mechanobiology.

[6]  P Rüegsegger,et al.  Non-invasive bone biopsy: a new method to analyse and display the three-dimensional structure of trabecular bone. , 1994, Physics in medicine and biology.

[7]  M Doblaré,et al.  Application of an anisotropic bone-remodelling model based on a damage-repair theory to the analysis of the proximal femur before and after total hip replacement. , 2001, Journal of biomechanics.

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

[9]  Yoshikazu Nakajima,et al.  A Parallel Implementation of 2-D/3-D Image Registration for Computer-Assisted Surgery , 2005, 11th International Conference on Parallel and Distributed Systems (ICPADS'05).

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

[11]  David Larsson,et al.  Measurement of structural anisotropy in femoral trabecular bone using clinical-resolution CT images. , 2013, Journal of biomechanics.

[12]  Z. Yosibash,et al.  Patient-specific Finite-element Analyses of the Proximal Femur with Orthotropic Material Properties Validated , 2022 .

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

[14]  Christian Hellmich,et al.  Can the diverse elastic properties of trabecular and cortical bone be attributed to only a few tissue-independent phase properties and their interactions? , 2004, Biomechanics and modeling in mechanobiology.

[15]  Yoon-Hyuk Kim,et al.  A new patient-specific planning method based on joint contact force balance with soft tissue release in total knee arthroplasty , 2013 .

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

[17]  C Kober,et al.  Consideration of anisotropic elasticity minimizes volumetric rather than shear deformation in human mandible , 2006, Computer methods in biomechanics and biomedical engineering.

[18]  Yifei Dai,et al.  Robust QCT/FEA Models of Proximal Femur Stiffness and Fracture Load During a Sideways Fall on the Hip , 2011, Annals of Biomedical Engineering.

[19]  Xavier Pennec,et al.  Capturing the multiscale anatomical shape variability with polyaffine transformation trees , 2012, Medical Image Anal..

[20]  M. Bouxsein,et al.  In vivo assessment of trabecular bone microarchitecture by high-resolution peripheral quantitative computed tomography. , 2005, The Journal of clinical endocrinology and metabolism.

[21]  Berthold K. P. Horn,et al.  Closed-form solution of absolute orientation using unit quaternions , 1987 .

[22]  Philippe K. Zysset,et al.  An alternative model for anisotropic elasticity based on fabric tensors , 1995 .

[23]  Z. Tabor,et al.  Quantifying anisotropy of trabecular bone from gray-level images. , 2007, Bone.

[24]  Yoshikazu Nakajima,et al.  A parallel implementation of 2-D/3-D image registration for computer-assisted surgery , 2006 .

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

[26]  Uwe Wolfram,et al.  Vertebral trabecular main direction can be determined from clinical CT datasets using the gradient structure tensor and not the inertia tensor--a case study. , 2009, Journal of biomechanics.

[27]  Sébastien Ourselin,et al.  Block Matching: A General Framework to Improve Robustness of Rigid Registration of Medical Images , 2000, MICCAI.

[28]  David Larsson,et al.  Assessment of Transverse Isotropy in Clinical-Level CT Images of Trabecular Bone Using the Gradient Structure Tensor , 2014, Annals of Biomedical Engineering.

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

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

[31]  C. Hellmich,et al.  Micromechanics-Based Conversion of CT Data into Anisotropic Elasticity Tensors, Applied to FE Simulations of a Mandible , 2008, Annals of Biomedical Engineering.

[32]  T San Antonio,et al.  Orientation of orthotropic material properties in a femur FE model: a method based on the principal stresses directions. , 2012, Medical engineering & physics.

[33]  Thomas M. Link,et al.  High-resolution Computed Tomography for Clinical Imaging of Bone Microarchitecture , 2011, Clinical orthopaedics and related research.

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