A Robust Algorithm for Thickness Computation at Low Resolution and Its Application to In Vivo Trabecular Bone CT Imaging

Adult bone diseases, especially osteoporosis, lead to increased risk of fracture which in turn is associated with substantial morbidity, mortality, and financial costs. Clinically, osteoporosis is defined by low bone mineral density; however, increasing evidence suggests that the microarchitectural quality of trabecular bone (TB) is an important determinant of bone strength and fracture risk. Accurate measures of TB thickness and marrow spacing is of significant interest for early diagnosis of osteoporosis or treatment effects. Here, we present a new robust algorithm for computing TB thickness and marrow spacing at a low resolution achievable in vivo. The method uses a star-line tracing technique that effectively deals with partial voluming effects of in vivo imaging with voxel size comparable to TB thickness. Also, the method avoids the problem of digitization associated with conventional algorithms based on sampling distance transform along skeletons. Accuracy of the method was examined using computer-generated phantom images, while the robustness of the method was evaluated on human ankle specimens in terms of stability across a wide range of voxel sizes, repeat scan reproducibility under in vivo conditions, and correlation between thickness values computed at ex vivo and in vivo imaging resolutions. Also, the sensitivity of the method was examined by evaluating its ability to predict the bone strength of cadaveric specimens. Finally, the method was evaluated in a human study involving 40 healthy young-adult volunteers (age: 19-21 years; 20 males and 20 females) and ten athletes (age: 19- 21 years; six males and four females). Across a wide range of voxel sizes, the new method is significantly more accurate and robust as compared to conventional methods. Both TB thickness and marrow spacing measures computed using the new method demonstrated strong associations (R2 ∈ [0.83, 0.87]) with bone strength. Also, the TB thickness and marrow spacing measures allowed discrimination between male and female volunteers (p ∈ [0.01, 0.04]) as well as between athletes and nonathletes (p ∈ [0.005, 0.03]).

[1]  W. Edwards,et al.  Cortical and trabecular bone contribute strength to the osteopenic distal radius , 1994, Journal of orthopaedic research : official publication of the Orthopaedic Research Society.

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

[3]  Punam K. Saha,et al.  Tensor scale: A local morphometric parameter with applications to computer vision and image processing , 2005, Comput. Vis. Image Underst..

[4]  P. Rüegsegger,et al.  In vivo high resolution 3D-QCT of the human forearm. , 1998, Technology and health care : official journal of the European Society for Engineering and Medicine.

[5]  S. Majumdar,et al.  In Vivo High Resolution MRI of the Calcaneus: Differences in Trabecular Structure in Osteoporosis Patients , 1998, Journal of bone and mineral research : the official journal of the American Society for Bone and Mineral Research.

[6]  Patrik Rogalla,et al.  Assessment of trabecular bone structure of the calcaneus using multi-detector CT: correlation with microCT and biomechanical testing. , 2009, Bone.

[7]  Jayaram K. Udupa,et al.  Scale-Based Fuzzy Connected Image Segmentation: Theory, Algorithms, and Validation , 2000, Comput. Vis. Image Underst..

[8]  H. Gundersen,et al.  Marrow space star volume in the iliac crest decreases in osteoporotic patients after continuous treatment with fluoride, calcium, and vitamin D2 for five years. , 1991, Bone.

[9]  Sharmila Majumdar,et al.  Magnetic Resonance Imaging of Trabecular Bone Structure , 2002, Topics in magnetic resonance imaging : TMRI.

[10]  Bidyut Baran Chaudhuri,et al.  Detection of 3-D Simple Points for Topology Preserving Transformations with Application to Thinning , 1994, IEEE Trans. Pattern Anal. Mach. Intell..

[11]  Punam K. Saha,et al.  Quantitative Characterization of Trabecular Bone Micro-architecture Using Tensor Scale and Multi-Detector CT Imaging , 2012, MICCAI.

[12]  Hajime Orimo,et al.  Multi‐Detector Row CT Imaging of Vertebral Microstructure for Evaluation of Fracture Risk , 2005, Journal of bone and mineral research : the official journal of the American Society for Bone and Mineral Research.

[13]  Punam K. Saha,et al.  A new algorithm for trabecular bone thickness computation at low resolution achieved under in vivo condition , 2013, 2013 IEEE 10th International Symposium on Biomedical Imaging.

[14]  J. L. Melton,et al.  Epidemiology of Spinal Osteoporosis , 1997, Spine.

[15]  P. Meunier,et al.  Clinical Use of Bone Biopsy , 2001 .

[16]  P. Croucher,et al.  Assessment of cancellous bone structure: Comparison of strut analysis, trabecular bone pattern factor, and marrow space star volume , 1996, Journal of bone and mineral research : the official journal of the American Society for Bone and Mineral Research.

[17]  Punam K. Saha,et al.  Fuzzy Distance Transform: Theory, Algorithms, and Applications , 2002, Comput. Vis. Image Underst..

[18]  Bidyut Baran Chaudhuri,et al.  A new shape preserving parallel thinning algorithm for 3D digital images , 1997, Pattern Recognit..

[19]  Gunilla Borgefors,et al.  On Digital Distance Transforms in Three Dimensions , 1996, Comput. Vis. Image Underst..

[20]  Bhabatosh Chanda,et al.  Topology preservation in 3D digital space , 1994, Pattern Recognit..

[21]  P. Rüegsegger,et al.  A new method for the model‐independent assessment of thickness in three‐dimensional images , 1997 .

[22]  R. Recker,et al.  Architecture and vertebral fracture , 1993, Calcified Tissue International.

[23]  Steven Boonen,et al.  Osteoporosis management: impact of fracture type on cost and quality of life in patients at risk for fracture I , 2008, Current medical research and opinion.

[24]  G. Beaupré,et al.  The influence of bone volume fraction and ash fraction on bone strength and modulus. , 2001, Bone.

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

[26]  R. Recker,et al.  Bone Microstructure in Osteoporosis: Transilial Biopsy and Histomorphometry , 2002, Topics in magnetic resonance imaging : TMRI.

[27]  Punam K. Saha,et al.  Measurement of trabecular bone thickness in the limited resolution regime of in vivo MRI by fuzzy distance transform , 2004, IEEE Transactions on Medical Imaging.

[28]  Yan Yang,et al.  Quantification of coronary arterial stenoses in CTA using fuzzy distance transform , 2012, Comput. Medical Imaging Graph..

[29]  Gabriella Sanniti di Baja,et al.  Distance-Driven Skeletonization in Voxel Images , 2011, IEEE Transactions on Pattern Analysis and Machine Intelligence.

[30]  Örjan Smedby,et al.  Estimation of trabecular thickness in gray-scale images through granulometric analysis , 2012, Medical Imaging.

[31]  P. Rüegsegger,et al.  A microtomographic system for the nondestructive evaluation of bone architecture , 2006, Calcified Tissue International.

[32]  H J Gundersen,et al.  Star volume of marrow space and trabeculae of the first lumbar vertebra: sampling efficiency and biological variation. , 1989, Bone.

[33]  I H Parkinson,et al.  Trabecular spacing in post-menopausal Australian women with and without vertebral fractures. , 1992, Australian and New Zealand journal of medicine.

[34]  M. Kleerekoper,et al.  Relationships between surface, volume, and thickness of iliac trabecular bone in aging and in osteoporosis. Implications for the microanatomic and cellular mechanisms of bone loss. , 1983, The Journal of clinical investigation.

[35]  A. Parfitt Implications of architecture for the pathogenesis and prevention of vertebral fracture. , 1992, Bone.

[36]  Steven K Boyd,et al.  Improved reproducibility of high-resolution peripheral quantitative computed tomography for measurement of bone quality. , 2008, Medical engineering & physics.

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

[38]  David G. Kirkpatrick,et al.  Linear Time Euclidean Distance Algorithms , 1995, IEEE Trans. Pattern Anal. Mach. Intell..

[39]  A. Wright,et al.  Deterioration of trabecular architecture in hypogonadal men. , 2003, The Journal of clinical endocrinology and metabolism.

[40]  K. J. Macleod,et al.  High resolution computed tomographic assessment of airway wall thickness in chronic asthma: reproducibility and relationship with lung function and severity , 2002, Thorax.

[41]  F Peyrin,et al.  A method for the automatic characterization of bone architecture in 3D mice microtomographic images. , 2003, Computerized Medical Imaging and Graphics.

[42]  Felix W. Wehrli,et al.  Imaging of Trabecular Bone , 2007 .

[43]  L. Mosekilde Age-related changes in vertebral trabecular bone architecture--assessed by a new method. , 1988, Bone.

[44]  Steven D. Kugelmass,et al.  Relationship between NMR transverse relaxation, trabecular bone architecture, and strength. , 1993, Proceedings of the National Academy of Sciences of the United States of America.

[45]  Milan Sonka,et al.  Topomorphologic Separation of Fused Isointensity Objects via Multiscale Opening: Separating Arteries and Veins in 3-D Pulmonary CT , 2010, IEEE Transactions on Medical Imaging.