Trabecular bone structure analysis in the limited spatial resolution regime of in vivo MRI.

RATIONALE AND OBJECTIVES To develop a method for processing and visualization of trabecular bone networks on the basis of magnetic resonance (MR) images acquired in the limited spatial resolution regime of in vivo imaging at which trabecular thickness is comparable to voxel size. MATERIALS AND METHODS A sequence of processing steps for analyzing the topologic structure of trabecular bone networks is presented and evaluated using three types of datasets: images of synthetic structures with various levels of superimposed Gaussian noise, micro-computed tomographic images of human trabecular bone downsampled to in vivo resolution, and in vivo micro-MR images from a prior longitudinal study investigating the structural implications of testosterone treatment of hypogonadal men. The simulated images were analyzed at a voxel size of 150 microm(3), the clinical MR image data had been acquired with 137 x 137 x 410 microm(3) voxel size. The technique is a modification to the virtual bone biopsy processing chain that involves a sinc convolution step immediately preceding binarization, and employs the Manzanera-Bernard thinning algorithm for obtaining the three-dimensional skeleton before topologic classification. The detectability of plate and rod bone elements was also analyzed theoretically. RESULTS As compared with previously published techniques, the approach produced a more accurate bone skeleton in the micro-computed tomographic and simulation experiments, with clear improvement in preservation of rod and plate elements. Simulations suggest that rods are detectable down to a diameter of approximately 50% of the MR image voxel length, whereas plates can be detected at thicknesses of 20% or more of voxel length. For in vivo studies, it was shown that the method could recover the treatment response in terms of the ensuing topologic changes in patients undergoing antiresorptive treatment. CONCLUSIONS The algorithm for processing of in vivo micro-MR images of trabecular bone is superior to prior approaches in preserving the topology of the network in the presence of noise.

[1]  Felix W Wehrli,et al.  Subvoxel processing: A method for reducing partial volume blurring with application to in vivo MR images of trabecular bone , 2002, Magnetic resonance in medicine.

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

[3]  Antoine Manzanera,et al.  N-dimensional Skeletonization: a Unified Mathematical Framework , 2002, J. Electronic Imaging.

[4]  Steven D. Kugelmass,et al.  Potential role of nuclear magnetic resonance for the evaluation of trabecular bone quality , 2005, Calcified Tissue International.

[5]  C. Nahmias,et al.  In vivo assessment of trabecular bone structure at the distal radius from high-resolution magnetic resonance images. , 1997, Medical physics.

[6]  S. Majumdar,et al.  Correlation of Trabecular Bone Structure with Age, Bone Mineral Density, and Osteoporotic Status: In Vivo Studies in the Distal Radius Using High Resolution Magnetic Resonance Imaging , 1997, Journal of bone and mineral research : the official journal of the American Society for Bone and Mineral Research.

[7]  H. Song,et al.  In vivo micro‐imaging using alternating navigator echoes with applications to cancellous bone structural analysis , 1999, Magnetic resonance in medicine.

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

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

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

[11]  J. Kanis,et al.  Secondary osteoporosis and the microanatomy of trabecular bone , 1989, Clinical Rheumatology.

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

[13]  F. Wehrli Structural and functional assessment of trabecular and cortical bone by micro magnetic resonance imaging , 2007, Journal of magnetic resonance imaging : JMRI.

[14]  Punam K. Saha,et al.  Three‐dimensional digital topological characterization of cancellous bone architecture , 2000 .

[15]  P. Levitz,et al.  A new method for three-dimensional skeleton graph analysis of porous media: application to trabecular bone microarchitecture. , 2000, Journal of microscopy.

[16]  H. Genant,et al.  Quantitative Computed Tomography of Vertebral Spongiosa: A Sensitive Method for Detecting Early Bone Loss After Oophorectomy , 1982, Annals of internal medicine.

[17]  Branimir Vasilic,et al.  Effect of Testosterone Replacement on Trabecular Architecture in Hypogonadal Men , 2005, Journal of bone and mineral research : the official journal of the American Society for Bone and Mineral Research.

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

[19]  J. Adachi,et al.  In vivo assessment of trabecular bone structure at the distal radius from high-resolution computed tomography images. , 1996, Physics in medicine and biology.

[20]  Punam K. Saha,et al.  Noninvasive assessment of bone architecture by magnetic resonance micro-imaging-based virtual bone biopsy , 2003, Proc. IEEE.

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

[22]  H. Song,et al.  Fast 3D large‐angle spin‐echo imaging (3D FLASE) , 1996, Magnetic resonance in medicine.

[23]  Jeremy Magland,et al.  In Vivo Magnetic Resonance Detects Rapid Remodeling Changes in the Topology of the Trabecular Bone Network After Menopause and the Protective Effect of Estradiol , 2008, Journal of bone and mineral research : the official journal of the American Society for Bone and Mineral Research.

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

[25]  S. Cowin,et al.  Wolff's law of trabecular architecture at remodeling equilibrium. , 1986, Journal of biomechanical engineering.

[26]  É. Legrand,et al.  Trabecular Bone Microarchitecture, Bone Mineral Density, and Vertebral Fractures in Male Osteoporosis , 2000, Journal of bone and mineral research : the official journal of the American Society for Bone and Mineral Research.

[27]  C. Rubin,et al.  Quantity and Quality of Trabecular Bone in the Femur Are Enhanced by a Strongly Anabolic, Noninvasive Mechanical Intervention , 2002, Journal of bone and mineral research : the official journal of the American Society for Bone and Mineral Research.

[28]  G. A. Ladinsky,et al.  Image metric‐based correction (autofocusing) of motion artifacts in high‐resolution trabecular bone imaging , 2007, Journal of magnetic resonance imaging : JMRI.

[29]  A. Chines,et al.  Risedronate preserves bone architecture in postmenopausal women with osteoporosis as measured by three-dimensional microcomputed tomography. , 2004, Bone.

[30]  N. Hylton,et al.  Impact of section doubling on MR angiography. , 1992, Radiology.

[31]  J Hale,et al.  Applications of voxel shifting in magnetic resonance imaging. , 1990, Investigative radiology.

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

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

[34]  Felix W. Wehrli,et al.  A novel local thresholding algorithm for trabecular bone volume fraction mapping in the limited spatial resolution regime of in vivo MRI , 2005, IEEE Transactions on Medical Imaging.

[35]  G. A. Ladinsky,et al.  Trabecular Structure Quantified With the MRI‐Based Virtual Bone Biopsy in Postmenopausal Women Contributes to Vertebral Deformity Burden Independent of Areal Vertebral BMD , 2007, Journal of bone and mineral research : the official journal of the American Society for Bone and Mineral Research.

[36]  Yiping P. Du,et al.  Reduction of partial‐volume artifacts with zero‐filled interpolation in three‐dimensional MR angiography , 1994, Journal of magnetic resonance imaging : JMRI.

[37]  M. Bouxsein,et al.  The effects of parathyroid hormone and alendronate alone or in combination in postmenopausal osteoporosis , 2004 .

[38]  Scott N. Hwang,et al.  Digital Topological Analysis of In Vivo Magnetic Resonance Microimages of Trabecular Bone Reveals Structural Implications of Osteoporosis , 2001, Journal of bone and mineral research : the official journal of the American Society for Bone and Mineral Research.

[39]  Scott N. Hwang,et al.  Estimating voxel volume fractions of trabecular bone on the basis of magnetic resonance images acquired in vivo , 1999 .