Assessment of trabecular bone strength at in vivo CT imaging with space-variant hysteresis and finite element modelling

Osteoporosis is a common bone disease associated with reduced bone strength and increased fracture risk. Finite element modelling (FEM) is a powerful tool to assess bone strength. Fast acquisition and highly reduced radiation in latest multi-row detector CT (MDCT) put it as a frontline imaging modality to assess in vivo trabecular bone (TB) microarchitecture and strength. In current CT image resolution, conventional segmentation methods fail to maintain TB network connectivity limiting reliable assessment of TB strength using FEM. In this paper, we present a new space-variant hysteresis approach to maintain TB connectivity while preserving marrow pores and a high quality mesh generator for FEM. We examine the effectiveness in estimating TB strength. The reproducibility and the ability to predict actual bone strength were examined on MDCT images of cadaveric ankle specimens under in vivo conditions. An intra-class correlation coefficient of 0.97 was observed in computed Young's modulus from repeat scans, and a high linear correlation (R2=0.92) were found between computed and experimental Young's modulus.

[1]  J. Torner,et al.  Automated cortical bone segmentation for multirow-detector CT imaging with validation and application to human studies. , 2015, Medical physics.

[2]  Jayaram K. Udupa,et al.  Optimum Image Thresholding via Class Uncertainty and Region Homogeneity , 2001, IEEE Trans. Pattern Anal. Mach. Intell..

[3]  Patrik Rogalla,et al.  Trabecular Bone Structure Analysis in the Osteoporotic Spine Using a Clinical In Vivo Setup for 64‐Slice MDCT Imaging: Comparison to μCT Imaging and μFE Modeling , 2009, Journal of bone and mineral research : the official journal of the American Society for Bone and Mineral Research.

[4]  Felix W Wehrli,et al.  In Vivo μMRI‐Based Finite Element and Morphological Analyses of Tibial Trabecular Bone in Eugonadal and Hypogonadal Men Before and After Testosterone Treatment , 2008, Journal of bone and mineral research : the official journal of the American Society for Bone and Mineral Research.

[5]  F. Eckstein,et al.  Prediction of bone strength by μCT and MDCT-based finite-element-models: how much spatial resolution is needed? , 2014, European journal of radiology.

[6]  Punam K. Saha,et al.  Volumetric Topological Analysis: A Novel Approach for Trabecular Bone Classification on the Continuum Between Plates and Rods , 2010, IEEE Transactions on Medical Imaging.

[7]  Sundeep Khosla,et al.  Population‐Based Study of Age and Sex Differences in Bone Volumetric Density, Size, Geometry, and Structure at Different Skeletal Sites , 2004, Journal of bone and mineral research : the official journal of the American Society for Bone and Mineral Research.

[8]  J. Shewchuk,et al.  Isosurface stuffing: fast tetrahedral meshes with good dihedral angles , 2007, SIGGRAPH 2007.

[9]  Steve Oudot,et al.  Provably good sampling and meshing of surfaces , 2005, Graph. Model..

[10]  John F. Canny,et al.  A Computational Approach to Edge Detection , 1986, IEEE Transactions on Pattern Analysis and Machine Intelligence.

[11]  M. Kleerekoper,et al.  The role of three-dimensional trabecular microstructure in the pathogenesis of vertebral compression fractures , 1985, Calcified Tissue International.

[12]  J. Shewchuk What Is a Good Linear Finite Element? Interpolation, Conditioning, Anisotropy, and Quality Measures , 2002 .