Premium
Structural analysis of cortical porosity applied to HR‐pQCT data
Author(s) -
Tjong Willy,
Nirody Jasmine,
Burghardt Andrew J.,
CarballidoGamio Julio,
Kazakia Galateia J.
Publication year - 2014
Publication title -
medical physics
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.473
H-Index - 180
eISSN - 2473-4209
pISSN - 0094-2405
DOI - 10.1118/1.4851575
Subject(s) - porosity , medical imaging , materials science , biomedical engineering , nuclear medicine , medicine , radiology , composite material
Purpose: The investigation of cortical porosity is an important aspect of understanding biological, pathoetiological, and biomechanical processes occurring within the skeleton. With the emergence of HR‐pQCT as a noninvasive tool suitable for clinical use, cortical porosity at appendicular sites can be directly visualized in vivo . The aim of this study was to introduce a novel topological analysis of the cortical pore network for HR‐pQCT data and determine the influence of resolution on measures of cortical pore network microstructure and topology.Methods: Cadaveric radii were scanned using HR‐pQCT at two different voxel sizes (41 and 82 μ m) and also using μ CT at a voxel size of 18 μ m. HR‐pQCT and μ CT image sets were spatially coregistered. Segmentation and quantification of cortical porosity (Ct.Po) and mean pore diameter (Ct.Po.Dm) were achieved using an established extended cortical analysis technique. Topological classification of individual pores was performed using topology‐preserving skeletonization and multicolor dilation algorithms. Based on the pore skeleton topological classification, the following parameters were quantified: total number of planar surface‐skeleton canals (N.Slabs), tubular curve‐skeleton canals (N.Tubes), and junction elements (N.Junctions), mean slab volume (Slab.Vol), mean tube volume (Tube.Vol), mean slab orientation (Slab.θ), mean tube orientation (Tube.θ), N.Slabs/N.Tubes, and integral (total) slab volume/integral tube volume (iSlab.Vol/iTube.Vol). An in vivo reproducibility study was also conducted to assess short‐term precision of the topology parameters. Precision error was characterized using root mean square coefficient of variation (RMSCV%).Results: Correlations to μ CT values for Ct.Po were significant for both the 41 and 82 μ m HR‐pQCT data (41: r 2 = 0.82, p < 0.001, 82: r 2 = 0.75, p < 0.001). For Ct.Po.Dm, only the 41 μ m data were significantly predictive of μ CT values (r 2 = 0.72, p < 0.01) Data at both HR‐pQCT voxel sizes were strongly predictive of the μ CT values for N.Slabs (41: r 2 = 0.93, p < 0.001; 82: r 2 = 0.84, p < 0.001), N.Tubes (41: r 2 = 0.94, p < 0.001; 82: r 2 = 0.84, p < 0.001), and N.Junctions (41: r 2 = 0.93, p < 0.001; 82: r 2 = 0.78, p < 0.001), though proportional bias was evident in these correlations. Weak correlations were seen for iSlab.Vol/iTube.Vol at both voxel sizes (41: r 2 = 0.52, p < 0.01; 82: r 2 = 0.39, p < 0.05). Slab.Vol was significantly correlated to μ CT data at 41 μ m (r 2 = 0.60, p < 0.01) but not at 82 μ m, while Tube.Vol was significantly correlated at both voxel sizes (41: r 2 = 0.79, p < 0.001; 82: r 2 = 0.68, p < 0.01). In vivo precision error for these parameters ranged from 2.31 to 9.68 RMSCV%.Conclusions: Strong correlations between μ CT‐ and HR‐pQCT‐derived measurements were found, particularly in HR‐pQCT images obtained at 41 μ m. These data are in agreement with our previous study investigating the effect of voxel size on standard HR‐pQCT metrics of trabecular and cortical microstructure, and extend our previous findings to include topological descriptors of the cortical pore network.