Computational study of estimating 3D trabecular bone microstructure for the volume of interest from CT scan data

Inspired by the self-optimizing capabilities of bone, a new concept of bone microstructure reconstruction has been recently introduced by using 2D synthetic skeletal images. As a preliminary clinical study, this paper proposes a topology optimization-based method that can estimate 3D trabecular bone microstructure for the volume of interest (VOI) from 3D computed tomography (CT) scan data with enhanced computational efficiency and phenomenological accuracy. For this purpose, a localized finite element (FE) model is constructed by segmenting a target bone from CT scan data and determining the physiological local loads for the VOI. Then, topology optimization is conducted with multiresolution bone mineral density (BMD) deviation constraints to preserve the patient-specific spatial bone distribution obtained from the CT scan data. For the first time, to our knowledge, this study has demonstrated that 60-μm resolution trabecular bone images can be reconstructed from 600-μm resolution CT scan data (a 62-year-old woman with no metabolic bone disorder) for the 4 VOIs in the proximal femur. The reconstructed trabecular bone includes the characteristic trabecular patterns and has morphometric indices that are in good agreement with the anatomical data in the literature. As for computational efficiency, the localization for the VOI reduces the number of FEs by 99%, compared with that of the full FE model. Compared with the previous single-resolution BMD deviation constraint, the proposed multiresolution BMD deviation constraints enable at least 65% and 47% reductions in the number of iterations and computing time, respectively. These results demonstrate the clinical feasibility and potential of the proposed method.

[1]  M. Bendsøe,et al.  Topology Optimization: "Theory, Methods, And Applications" , 2011 .

[2]  Ralph Müller,et al.  Long-term prediction of three-dimensional bone architecture in simulations of pre-, peri- and post-menopausal microstructural bone remodeling , 2005, Osteoporosis International.

[3]  G S Beaupré,et al.  An approach for time‐dependent bone modeling and remodeling—theoretical development , 1990, Journal of orthopaedic research : official publication of the Orthopaedic Research Society.

[4]  R. Huiskes,et al.  Load distribution in the healthy and osteoporotic human proximal femur during a fall to the side. , 2008, Bone.

[5]  Tomasz Lekszycki,et al.  A 2‐D continuum model of a mixture of bone tissue and bio‐resorbable material for simulating mass density redistribution under load slowly variable in time , 2014 .

[6]  M. Holick,et al.  Vitamin D: importance in the prevention of cancers, type 1 diabetes, heart disease, and osteoporosis. , 2004, The American journal of clinical nutrition.

[7]  Thomas M Link,et al.  osteoporosis imaging : State of the Art and Advanced Imaging 1 , 2022 .

[8]  Christopher Boyle,et al.  Three-dimensional micro-level computational study of Wolff's law via trabecular bone remodeling in the human proximal femur using design space topology optimization. , 2011, Journal of biomechanics.

[9]  P. Parizel,et al.  A brief review of parallel magnetic resonance imaging , 2003, European Radiology.

[10]  Jung Jin Kim,et al.  Fully automated segmentation of a hip joint using the patient-specific optimal thresholding and watershed algorithm , 2018, Comput. Methods Programs Biomed..

[11]  Marco Viceconti,et al.  Subject-specific finite element models implementing a maximum principal strain criterion are able to estimate failure risk and fracture location on human femurs tested in vitro. , 2008, Journal of biomechanics.

[12]  C R Appledorn,et al.  A new approach to the interpolation of sampled data , 1996, IEEE Trans. Medical Imaging.

[13]  K. Svanberg.  The method of moving asymptotes—a new method for structural optimization , 1987 .

[14]  Jung Jin Kim,et al.  Image resolution enhancement for healthy weight-bearing bones based on topology optimization. , 2016, Journal of biomechanics.

[15]  W. C. Hayes,et al.  Stress distributions within the proximal femur during gait and falls: Implications for osteoporotic fracture , 2005, Osteoporosis International.

[16]  L. Feldkamp,et al.  Practical cone-beam algorithm , 1984 .

[17]  Truong Q. Nguyen,et al.  Novel Example-Based Method for Super-Resolution and Denoising of Medical Images , 2014, IEEE Transactions on Image Processing.

[18]  Ralph Müller,et al.  The Clinical Biomechanics Award 2012 - presented by the European Society of Biomechanics: large scale simulations of trabecular bone adaptation to loading and treatment. , 2014, Clinical biomechanics.

[19]  J J Hamilton,et al.  Bone remodeling and structural optimization. , 1994, Journal of biomechanics.

[20]  Sharmila Majumdar,et al.  High-resolution imaging techniques for the assessment of osteoporosis. , 2010, Radiologic clinics of North America.

[21]  R Huiskes,et al.  A theoretical framework for strain-related trabecular bone maintenance and adaptation. , 2005, Journal of biomechanics.

[22]  Il Yong Kim,et al.  Analogy of strain energy density based bone-remodeling algorithm and structural topology optimization. , 2009, Journal of biomechanical engineering.

[23]  In Gwun Jang,et al.  Application of design space optimization to bone remodeling simulation of trabecular architecture in human proximal femur for higher computational efficiency , 2010 .

[24]  S A Goldstein,et al.  The relationship between the structural and orthogonal compressive properties of trabecular bone. , 1994, Journal of biomechanics.

[25]  Jung Jin Kim,et al.  Estimation of Local Bone Loads for the Volume of Interest. , 2016, Journal of biomechanical engineering.

[26]  Keita Ito,et al.  Inter-individual variability of bone density and morphology distribution in the proximal femur and T12 vertebra. , 2014, Bone.

[27]  H. Rodrigues,et al.  Numerical modeling of bone tissue adaptation--a hierarchical approach for bone apparent density and trabecular structure. , 2009, Journal of Biomechanics.

[28]  T. Lekszycki.  Modelling of Bone Adaptation Based on an Optimal Response Hypothesis* , 2002 .

[29]  M. Viceconti,et al.  Patient-specific finite element estimated femur strength as a predictor of the risk of hip fracture: the effect of methodological determinants , 2016, Osteoporosis International.

[30]  Jing Li,et al.  A viscoelastic model for random ultrasound propagation in cancellous bone , 2013, Comput. Math. Appl..

[31]  P. Rüegsegger,et al.  Direct Three‐Dimensional Morphometric Analysis of Human Cancellous Bone: Microstructural Data from Spine, Femur, Iliac Crest, and Calcaneus , 1999, Journal of bone and mineral research : the official journal of the American Society for Bone and Mineral Research.

[32]  TOR Hildebrand,et al.  Quantification of Bone Microarchitecture with the Structure Model Index. , 1997, Computer methods in biomechanics and biomedical engineering.

[33]  Jie Tang,et al.  Prior image constrained compressed sensing (PICCS): a method to accurately reconstruct dynamic CT images from highly undersampled projection data sets. , 2008, Medical physics.

[34]  S. Hollister,et al.  Optimization of scaffold design for bone tissue engineering: A computational and experimental study. , 2014, Medical engineering & physics.

[35]  R. Müller,et al.  Load-adaptive bone remodeling simulations reveal osteoporotic microstructural and mechanical changes in whole human vertebrae. , 2016, Journal of biomechanics.

[36]  M. Hestenes,et al.  Methods of conjugate gradients for solving linear systems , 1952 .

[37]  Jing Li,et al.  Numerical investigation of ultrasonic attenuation through 2D trabecular bone structures reconstructed from CT scans and random realizations , 2014, Comput. Biol. Medicine.

[38]  Tomasz Lekszycki.  Optimality conditions in modeling of bone adaptation phenomenon , 1999 .

[39]  B. Kwak,et al.  Evolutionary topology optimization using design space adjustment based on fixed grid , 2006 .

[40]  In Gwun Jang,et al.  Homeostasis-based aging model for trabecular changes and its correlation with age-matched bone mineral densities and radiographs. , 2015, European journal of radiology.

[41]  O. Mäkitie,et al.  Heterozygous Mutations in the LDL Receptor‐Related Protein 5 (LRP5) Gene Are Associated With Primary Osteoporosis in Children , 2005, Journal of bone and mineral research : the official journal of the American Society for Bone and Mineral Research.

[42]  Michal Irani,et al.  Improving resolution by image registration , 1991, CVGIP Graph. Model. Image Process..

[43]  Chung-Ming Chen,et al.  Computer-Aided Detection and Diagnosis in Medical Imaging , 2013, Comput. Math. Methods Medicine.

[44]  Steven A. Goldstein,et al.  Measurement and significance of three-dimensional architecture to the mechanical integrity of trabecular bone , 2005, Calcified Tissue International.

[45]  Hsieh Hou,et al.  Cubic splines for image interpolation and digital filtering , 1978 .

[46]  Francesco dell’Isola,et al.  A visco-poroelastic model of functional adaptation in bones reconstructed with bio-resorbable materials , 2016, Biomechanics and Modeling in Mechanobiology.

[47]  B. Kwak,et al.  Design space optimization using design space adjustment and refinement , 2007 .

[48]  In Gwun Jang,et al.  Computational simulation of trabecular adaptation progress in human proximal femur during growth. , 2009, Journal of biomechanics.

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

[50]  William T. Freeman,et al.  Example-Based Super-Resolution , 2002, IEEE Computer Graphics and Applications.

[51]  Takao Yamamuro,et al.  Electron microscopy of calcification during high-density suspension culture of chondrocytes , 1993, Calcified Tissue International.

[52]  Tomonori Yamada,et al.  Computer simulation of trabecular remodeling in human proximal femur using large-scale voxel FE models: Approach to understanding Wolff's law. , 2009, Journal of biomechanics.

[53]  R. Rizzoli,et al.  Bone strength and its determinants , 2003, Osteoporosis International.

[54]  Bert Van Rietbergen,et al.  Finite Element Analysis Based on In Vivo HR‐pQCT Images of the Distal Radius Is Associated With Wrist Fracture in Postmenopausal Women , 2007, Journal of bone and mineral research : the official journal of the American Society for Bone and Mineral Research.

[55]  Francesco dell’Isola,et al.  A mixture model with evolving mass densities for describing synthesis and resorption phenomena in bones reconstructed with bio‐resorbable materials , 2012 .

[56]  In Gwun Jang,et al.  Computational study on the effect of loading alteration caused by disc degeneration on the trabecular architecture in human lumbar spine. , 2010, Journal of biomechanics.

[57]  In Gwun Jang,et al.  Computational study of Wolff's law with trabecular architecture in the human proximal femur using topology optimization. , 2008, Journal of biomechanics.

[58]  Lei Zhu,et al.  Compressed sensing based cone-beam computed tomography reconstruction with a first-order methoda). , 2010, Medical physics.