Automatic multi-parametric quantification of the proximal femur with quantitative computed tomography

SSDL-an automated semi-supervised deep learning approach for patient-specific 3D reconstruction of proximal femur from QCT images

Deep Learning (DL) techniques have recently been used in medical image segmentation and the reconstruction of 3D anatomies of a human body. In this work, we propose a semi-supervised DL (SSDL) approach utilizing a CNN-based 3D U-Net model for femur segmentation from sparsely annotated quantitative computed tomography (QCT) slices. Specifically, QCT slices at the proximal end of the femur forming ball and socket joint with acetabulum were annotated for precise segmentation, where a segmenting binary mask was generated using a 3D U-Net model to segment the femur accurately. A total of 5474 QCT slices were considered for training among which 2316 slices were annotated. 3D femurs were further reconstructed from segmented slices employing polynomial spline interpolation. Both qualitative and quantitative performance of segmentation and 3D reconstruction were satisfactory with more than 90% accuracy achieved for all of the standard performance metrics considered. The spatial overlap index and reproducibility validation metric for segmentation-Dice Similarity Coefficient was 91.8% for unseen patients and 99.2% for validated patients. An average relative error of 12.02% and 10.75% for volume and surface area, respectively, were computed for 3D reconstructed femurs. The proposed approach demonstrates its effectiveness in accurately segmenting and reconstructing 3D femur from QCT slices.

A new hip fracture risk index derived from FEA-computed proximal femur fracture loads and energies-to-failure

Hip fracture risk assessment is an important but challenging task. Quantitative CT-based patient-specific finite element (FE) analysis (FEA) incorporates bone geometry and bone density in the proximal femur. We developed a global FEA-computed fracture risk index to increase the prediction accuracy of hip fracture incidence.

PURPOSE

Quantitative CT-based patient-specific finite element (FE) analysis (FEA) incorporates bone geometry and bone density in the proximal femur to compute the force (fracture load) and energy necessary to break the proximal femur in a particular loading condition. The fracture loads and energies-to-failure are individually associated with incident hip fracture, and provide different structural information about the proximal femur.

METHODS

We used principal component analysis (PCA) to develop a global FEA-computed fracture risk index that incorporates the FEA-computed yield and ultimate failure loads and energies-to-failure in four loading conditions of 110 hip fracture subjects and 235 age- and sex-matched control subjects from the AGES-Reykjavik study. Using a logistic regression model, we compared the prediction performance for hip fracture based on the stratified resampling.

RESULTS

We referred the first principal component (PC1) of the FE parameters as the global FEA-computed fracture risk index, which was the significant predictor of hip fracture (p-value < 0.001). The area under the receiver operating characteristic curve (AUC) using PC1 (0.776) was higher than that using all FE parameters combined (0.737) in the males (p-value < 0.001).

CONCLUSIONS

The global FEA-computed fracture risk index increased hip fracture risk prediction accuracy in males.

People living with HIV have low trabecular bone mineral density, high bone marrow adiposity, and poor trabecular bone microarchitecture at the proximal femur

People living with HIV (PLWH) have increased risk of osteoporosis and fractures. We assessed the proximal femur of PLWH and age-matched seronegative controls using quantitative computed tomography and magnetic resonance imaging. Results suggest that the trabecular compartment is compromised at fracture-prone regions in the proximal femur of PLWH.

INTRODUCTION

People living with HIV (PLWH) have increased risk of osteoporosis and fractures. However, studies assessing the main determinants of bone strength in the proximal femur exclude this vulnerable population. We assessed the proximal femur of 40 PLWH and 26 age-matched seronegative controls using quantitative computed tomography and magnetic resonance imaging.

METHODS

We examined cortical volumetric bone mineral density (Ct.vBMD), trabecular vBMD (Tb.vBMD), cortical thickness (Ct.Th), bone marrow adiposity (BMA), and trabecular number, separation, and bone volume fraction. Parametric comparisons between the two groups were made for the femoral head, femoral neck, trochanter, and total hip using linear regression adjusting for several covariates, including metrics of body composition. In addition, we investigated the associations of BMA with Tb.vBMD and trabecular microarchitecture with Spearman's rank partial correlations.

RESULTS

PLWH had lower Tb.vBMD and deteriorated trabecular microarchitecture in the femoral neck, trochanter and total hip, and elevated BMA in the femoral head, femoral neck, and total hip. Ct.vBMD and Ct.Th were not significantly different between the two groups. BMA was significantly associated with lower Tb.vBMD and deteriorated trabecular microarchitecture in both groups albeit at different femoral regions.

CONCLUSIONS

Our findings suggest that the trabecular, and not the cortical, compartment is compromised in the proximal femur of PLWH. The observed impairments in fracture-prone regions in PLWH indicate lower femoral strength and suggest higher fracture risk. The inverse associations of BMA with trabecular bone density and microarchitecture quality agree with findings at other anatomic sites and in other populations, suggesting that excess BMA possibly due to a switch from the osteoblast to the adipocyte lineage may be implicated in the pathogenesis of bone fragility at the femur in PLWH.

Physical activity induced adaptation can increase proximal femur strength under loading from a fall onto the greater trochanter

Physical activity enhances proximal femur bone mass, but it remains unclear whether the benefits translate into an enhanced ability to resist fracture related loading. We recently used baseball pitchers as a within-subject controlled model to demonstrate physical activity induced proximal femur adaptation in regions associated with weight bearing and increased strength under single-leg stance loading. However, there was no measurable benefit to resisting common injurious loading (e.g. a fall onto the greater trochanter). A lack of power and a small physical activity effect size may have contributed to the latter null finding. Softball pitchers represent an alternative within-subject controlled model to explore adaptation of the proximal femur to physical activity, exhibiting greater dominant-to-nondominant (D-to-ND) leg differences than baseball pitchers. The current study used quantitative computed tomography, statistical parametric mapping, and subject-specific finite element (FE) modeling to explore adaptation of the proximal femur to physical activity in female softball pitchers (n = 25). Female cross-country runners (n = 15) were included as symmetrically loaded controls, showing very limited D-to-ND leg differences. Softball pitchers had D-to-ND leg differences in proximal femur, femoral neck, and trochanteric volumetric bone mineral density and content, and femoral neck volume. Voxel-based morphometry analyses and cortical bone mapping showed D-to-ND leg differences within a large region connecting the superior femoral head, inferior femoral neck and medial intertrochanteric region, and within the greater trochanter. FE modeling revealed pitchers had 19.4% (95%CI, 15.0 to 23.9%) and 4.9% (95%CI, 1.7 to 8.2%) D-to-ND leg differences in predicted ultimate strength under single-leg stance loading and a fall onto the greater trochanter, respectively. These data affirm the spatial and strength adaptation of the proximal femur to weight bearing directed loading and demonstrate that the changes can also have benefits, albeit smaller, on resisting loads associated with a sideways fall onto the greater trochanter.

Three-Dimensional Quantification of Bone Mineral Density in the Distal Femur and Proximal Tibia Based on Computed Tomography: In Vitro Evaluation of an Extended Standardization Method

While alterations in bone mineral density (BMD) are of interest in a number of musculoskeletal conditions affecting the knee, their analysis is limited by a lack of tools able to take full advantage of modern imaging modalities. This study introduced a new method, combining computed tomography (CT) and computational anatomy algorithms, to produce standardized three-dimensional BMD quantification in the distal femur and proximal tibia. The method was evaluated on ten cadaveric knees CT-scanned twice and processed following three different experimental settings to assess the influence of different scans and operators. The median reliability (intraclass correlation coefficient (ICC)) ranged from 0.96 to 0.99 and the median reproducibility (precision error (RMSSD)) ranged from 3.97 to 10.75 mg/cc for the different experimental settings. In conclusion, this paper presented a method to standardize three-dimensional knee BMD with excellent reliability and adequate reproducibility to be used in research and clinical applications. The perspectives offered by this novel method are further reinforced by the fact it relies on conventional CT scan of the knee. The standardization method introduced in this work is not limited to BMD and could be adapted to quantify other bone parameters in three dimension based on CT images or images acquired using different modalities.

Heterogeneous Spatial and Strength Adaptation of the Proximal Femur to Physical Activity: A Within-Subject Controlled Cross-Sectional Study

Physical activity (PA) enhances proximal femur bone mass, as assessed using projectional imaging techniques. However, these techniques average data over large volumes, obscuring spatially heterogeneous adaptations. The current study used quantitative computed tomography, statistical parameter mapping, and subject-specific finite element (FE) modeling to explore spatial adaptation of the proximal femur to PA. In particular, we were interested in adaptation occurring at the superior femoral neck and improving strength under loading from a fall onto the greater trochanter. High/long jump athletes (n = 16) and baseball pitchers (n = 16) were utilized as within-subject controlled models as they preferentially load their take-off leg and leg contralateral to their throwing arm, respectively. Controls (n = 15) were included but did not show any dominant-to-nondominant (D-to-ND) leg differences. Jumping athletes showed some D-to-ND leg differences but less than pitchers. Pitchers had 5.8% (95% confidence interval [CI] 3.9%-7.6%) D-to-ND leg differences in total hip volumetric bone mineral density (vBMD), with increased vBMD in the cortical compartment of the femoral neck and trochanteric cortical and trabecular compartments. Voxel-based morphometry analyses and cortical bone mapping showed pitchers had D-to-ND leg differences within the regions of the primary compressive trabeculae, inferior femoral neck, and greater trochanter but not the superior femoral neck. FE modeling revealed pitchers had 4.1% (95% CI 1.4%-6.7%) D-to-ND leg differences in ultimate strength under single-leg stance loading but no differences in ultimate strength to a fall onto the greater trochanter. These data indicate the asymmetrical loading associated with baseball pitching induces proximal femur adaptation in regions associated with weight bearing and muscle contractile forces and increases strength under single-leg stance loading. However, there were no benefits evident at the superior femoral neck and no measurable improvement in ultimate strength to common injurious loading during aging (ie, fall onto the greater trochanter), raising questions as to how to better target these variables with PA. © 2019 American Society for Bone and Mineral Research.

A novel method for estimating nail-tract bone density for intertrochanteric fractures

Summary

A novel method based on voxel-based morphometry was proposed to investigate the average volumetric bone mineral density (vBMD) of femoral head nail tract in patients treated with intramedullary nails—proximal femoral nail antirotation (PFNA) and gamma nail (GN). The results showed that there was no significant difference in average vBMD between the two groups.

Background

For unstable intertrochanteric fractures, poor bone quality might be one of the most important causes of cut-out complications in the femoral head during surgical treatment. Bone quality is generally regarded as an equivalent of BMD. Thus, we develop a novel voxel-based morphometry–based method to quantify vBMD of the femoral head nail tract.

Methods

Automatic calculation of average vBMD of nail tracts requires three main steps. First, we built a standard nail tract in a proximal femur template. Then, we mapped the proximal femur structure of each patient to the template by B-spline and Demons registration so that the anatomical positions of the proximal femur of all patients spatially corresponded to the standard template. Finally, we calculated and visualized the average vBMD distribution of the nail tract of all patients. To verify the feasibility of the method, we enrolled 75 patients (52 women and 23 men) with hip fractures to our study to compare measurements. The root mean square of the standard deviation (RMSSD) was calculated, and the coefficient of variation (CV) of the RMSSD (CV-RMSSD) was used to evaluate the reproducibility of intraoperator and interscan measurements. The Mann–Whitney U test was used to compare the average vBMD of nail tracts for the PFNA and GN.

Results

The CV-RMSSD of intraoperator measurements ranged from 1.0% to 2.0%, and the CV-RMSSD of interscan measurements ranged from 3.6% to 4.5%. There was no significant difference in the average vBMD between patients with PFNAs and those with GNs (p > 0.05).

Conclusions

The proposed method is reproducible for determining the average vBMD, which may provide a reference index for selection of appropriate intramedullary nails for individual patients. The current choice of intramedullary nail based on the experience of a surgeon may be biased.

The translational potential of this article

A novel method was proposed to measure the spatial average vBMD of nail tracts, which has good potential to provide a reference index for surgeons to choose appropriate implants.

Hip Fracture Discrimination Based on Statistical Multi-parametric Modeling (SMPM)

Studies using quantitative computed tomography (QCT) and data-driven image analysis techniques have shown that trabecular and cortical volumetric bone mineral density (vBMD) can improve the hip fracture prediction of dual-energy X-ray absorptiometry areal BMD (aBMD). Here, we hypothesize that (1) QCT imaging features of shape, density and structure derived from data-driven image analysis techniques can improve the hip fracture discrimination of classification models based on mean femoral neck aBMD (Neck.aBMD), and (2) that data-driven cortical bone thickness (Ct.Th) features can improve the hip fracture discrimination of vBMD models. We tested our hypotheses using statistical multi-parametric modeling (SMPM) in a QCT study of acute hip fracture of 50 controls and 93 fragility fracture cases. SMPM was used to extract features of shape, vBMD, Ct.Th, cortical vBMD, and vBMD in a layer adjacent to the endosteal surface to develop hip fracture classification models with machine learning logistic LASSO. The performance of these classification models was evaluated in two aspects: (1) their hip fracture classification capability without Neck.aBMD, and (2) their capability to improve the hip fracture classification of the Neck.aBMD model. Assessments were done with 10-fold cross-validation, areas under the receiver operating characteristic curve (AUCs), differences of AUCs, and the integrated discrimination improvement (IDI) index. All LASSO models including SMPM-vBMD features, and the majority of models including SMPM-Ct.Th features performed significantly better than the Neck.aBMD model; and all SMPM features significantly improved the hip fracture discrimination of the Neck.aBMD model (Hypothesis 1). An interesting finding was that SMPM-features of vBMD also captured Ct.Th patterns, potentially explaining the superior classification performance of models based on SMPM-vBMD features (Hypothesis 2). Age, height and weight had a small impact on model performances, and the model of shape, vBMD and Ct.Th consistently yielded better performances than the Neck.aBMD models. Results of this study clearly support the relevance of bone density and quality on the assessment of hip fracture, and demonstrate their potential on patient and healthcare cost benefits.

Adaptation of the proximal humerus to physical activity: A within-subject controlled study in baseball players

The proximal humerus is a common, yet understudied site for osteoporotic fracture. The current study explored the impact of prolonged physical activity on proximal humerus bone health by comparing bone properties between the throwing and nonthrowing arms within professional baseball players. The proximal humerus in throwing arms had 28.1% (95% CI, 17.8 to 38.3%) greater bone mass compared to nonthrowing arms, as assessed using dual-energy x-ray absorptiometry. At the level of the surgical neck, computed tomography revealed 12.0% (95% CI, 8.2 to 15.8%) greater total cross-sectional area and 31.0% (95% CI, 17.8 to 44.2%) greater cortical thickness within throwing arms, which contributed to 56.8% (95% CI, 44.9 to 68.8%) greater polar moment of inertia (i.e., estimated ability to resist torsional forces) compared to nonthrowing arms. Within the humeral head and greater tubercle regions, throwing arms had 3.1% (95% CI, 1.1 to 5.1%) more trabecular bone, as assessed using high-resolution magnetic resonance imaging. Three-dimensional mapping of voxel- and vertex-wise differences between arms using statistical parametric mapping techniques revealed throwing arms had adaptation within much of the proximal diaphysis, especially the posterolateral cortex. The pattern of proximal diaphysis adaptation approximated the pattern of strain energy distribution within the proximal humerus during a fastball pitch derived from a musculoskeletal and finite element model in a representative player. These data demonstrate the adaptive ability of the proximal humerus to physical activity-related mechanical loads. It remains to be established how they translate to exercise prescription to improve bone health within the proximal humerus; however, they provide unique insight into the relationship between prolonged loading and skeletal adaptation at a clinically relevant osteoporotic site.

MRI in idiopathic, stable, slipped capital femoral epiphysis: evaluation of contralateral pre-slip

PURPOSE

Early diagnosis and treatment of slipped capital femoral epiphysis (SCFE) is important to prevent slip progression and avoid complications. We sought to determine if MRI findings in patients with unilateral SCFE could indicate 'pre-slip' or predict future SCFE in the contralateral hip.

METHODS

A prospective study evaluated patients with unilateral SCFE over a two-year period. MRI of the asymptomatic hip was performed within the perioperative period. Patients were followed with radiographs until a contralateral slip occurred or until physeal closure. Demographics, clinical stability, severity, posterior slope angle (PSA), modified Oxford Bone Score (mOBS) and patency of the triradiate cartilage were recorded and statistical analysis performed.

RESULTS

In all, 33 of 54 patients with unilateral SCFE were enrolled into the study. In all, 29 (87.8%) had complete follow-up. Five of the enrolled patients (15.2%) developed a sequential slip requiring in situ pinning. Six of 33 (18.2%) patients had positive MRI findings: four of which proceeded to sequential SCFE and two which did not. One sequential slip had a negative MRI. PSA predicted 1/11 sequential slips (sensitivity 9.09%, specificity 81.4%, positive predictive value (PPV) 11.1%, negative predictive value (NPV) 77.8%) and mOBS predicted 5/11 sequential slips (sensitivity 45.5%, specificity 93%, PPV 62.5%, NPV 87%). An open triradiate cartilage was present in 8/11 patients with sequential slips (sensitivity 72.7%, specificity 81.4%, PPV 50%, NPV 92.1%).

CONCLUSION

MRI findings consistent with 'pre-slip' were present in 66.7% of patients who developed a sequential SCFE. Further study on the utility/sensitivity of MRI in predicting sequential SCFE is warranted.

LEVEL OF EVIDENCE

II, diagnostic.

Cortical Bone Mapping: Measurement and Statistical Analysis of Localised Skeletal Changes

PURPOSE OF REVIEW

Cortical bone mapping (CBM) is a technique for measuring localised skeletal changes from computed tomography (CT) images. It can provide measurements with accuracy surpassing the underlying imaging resolution. CBM can detect changes in several properties of the cortex, with no prior assumptions about the likely location of said changes. This paper summarises the theory behind CBM, discusses its strengths and limitations, and reviews some studies in which it has been applied.

RECENT FINDINGS

CBM has revealed associations between fracture risk and cortical properties in specific regions of the proximal femur which present feasible therapeutic targets. Analyses of several pharmaceutical and exercise interventions quantify effects that are distinct both in location and in the nature of the micro-architectural changes. CBM has illuminated age-related changes in the proximal femur and has recently been applied to other bones, as well as to the assessment of cartilage. The CBM processing pipeline is designed primarily for large cohort studies. Its main impact thus far has not been in the realm of clinical practice, but rather to improve our fundamental understanding of localised bone structure and changes.

Sex differences in the spatial distribution of bone in relation to incident hip fracture: Findings from the AGES-Reykjavik study

In this case-cohort study, we used data-driven computational anatomy approaches to assess within and between sex spatial differences in proximal femoral bone characteristics in relation to incident hip fracture. One hundred male and 234 female incident hip fracture cases, and 1047 randomly selected noncase subcohort participants (562 female) were chosen from the population-based AGES-Reykjavik study (mean age of 77 years). The baseline -i.e. before hip fracture- hip quantitative computed tomography scans of these subjects were analyzed using voxel-based morphometry, tensor-based morphometry, and surface-based statistical parametric mapping to assess the spatial distribution of volumetric bone mineral density (vBMD), internal structure, and cortical bone properties (thickness, vBMD and trabecular vBMD adjacent to the endosteal surface) of the proximal femur, respectively, in relation to incident hip fracture. Results showed that in both men and women: 1) the superior aspect of the femoral neck and the trochanteric region (except for cortical bone thickness) were consistently identified as being associated with incident hip fracture, and 2) differences in bone properties between noncases and incident hip fracture cases followed similar trends, were located at compatible regions, and manifested heterogeneity in the spatial distribution of their magnitude with focal regions showing larger differences. With respect to sex differences, most of the regions with a significant interaction between fracture group and sex showed: 1) differences of greater magnitude in men between noncases and incident hip fracture cases with different spatial distributions for all bone properties with the exception of cortical bone thickness, and 2) that while most of these regions showed better bone quality in male cases than in female cases, female cases showed higher vBMD in the principal compressive group and higher endotrabecular vBMD at several regions including the anterior, posterior, and lateral aspects of the proximal femur. These findings indicate the value of these image analysis techniques by providing unique information about the specific patterns of bone deterioration associated with incident hip fracture and their sex differences, highlighting the importance of looking to men and women separately in the assessment of hip fracture risk.

Trabecular bone microstructure is impaired in the proximal femur of human immunodeficiency virus-infected men with normal bone mineral density

Background

There is evidence that human immunodeficiency virus (HIV) infection and antiretroviral therapy (ART) are independent risk factors for osteoporosis and fracture which is not solely explained by changes in bone mineral density. Thus, we hypothesized that the assessment of trabecular microstructure might play an important role for bone quality in this population and might explain the increased fracture risk. In this study, we have assessed bone microstructure in the proximal femur using high-resolution magnetic resonance imaging (MRI) as well as in the extremities using high resolution peripheral quantitative computed tomography (HR-pQCT) in HIV-infected men and healthy controls and compared these findings to those based on areal bone mineral density (aBMD) derived from dual X-ray absorptiometry (DXA) which is the standard clinical parameter for the diagnosis of osteoporosis.

Methods

Eight HIV-infected men and 11 healthy age-matched controls were recruited and informed consent was obtained before each scan. High-resolution MRI of the proximal femur was performed using fully balanced steady state free precession (bSSFP) on a 3T system. Three volumes of interest at corresponding anatomic locations across all subjects were defined based on registrations of a common template. Four MR-based trabecular microstructural parameters were analyzed at each region: fuzzy bone volume fraction (f-BVF), trabecular number (Tb.N), thickness (Tb.Th), and spacing (Tb.Sp). In addition, the distal radius and distal tibia were imaged with HR-pQCT. Four HR-pQCT-based microstructural parameters were analyzed: trabecular bone volume fraction (BV/TV), Tb.N, Tb.Th, and Tb.Sp. Total hip and spine aBMD were determined from DXA.

Results

Microstructural bone parameters derived from MRI at the proximal femur and from HR-pQCT at the distal tibia showed significantly lower bone quality in HIV-infected patients compared to healthy controls. In contrast, DXA aBMD data showed no significant differences between HIV-infected patients and healthy controls.

Conclusions

Our results suggest that high-resolution imaging is a powerful tool to assess trabecular bone microstructure and can be used to assess bone health in HIV-infected men who show no differences to healthy males by DXA aBMD. Advances in MRI technology have made microstructural imaging at the proximal femur possible. Further studies in larger patient cohorts are clearly warranted.

Spatial Differences in the Distribution of Bone Between Femoral Neck and Trochanteric Fractures

There is little knowledge about the spatial distribution differences in volumetric bone mineral density and cortical bone structure at the proximal femur between femoral neck fractures and trochanteric fractures. In this case-control study, a total of 93 women with fragility hip fractures, 72 with femoral neck fractures (mean ± SD age: 70.6 ± 12.7 years) and 21 with trochanteric fractures (75.6 ± 9.3 years), and 50 control subjects (63.7 ± 7.0 years) were included for the comparisons. Differences in the spatial distributions of volumetric bone mineral density, cortical bone thickness, cortical volumetric bone mineral density, and volumetric bone mineral density in a layer adjacent to the endosteal surface were investigated using voxel-based morphometry (VBM) and surface-based statistical parametric mapping (SPM). We compared these spatial distributions between controls and both types of fracture, and between the two types of fracture. Using VBM, we found spatially heterogeneous volumetric bone mineral density differences between control subjects and subjects with hip fracture that varied by fracture type. Interestingly, femoral neck fracture subjects, but not subjects with trochanteric fracture, showed significantly lower volumetric bone mineral density in the superior aspect of the femoral neck compared with controls. Using surface-based SPM, we found that compared with controls, both fracture types showed thinner cortices in regions in agreement with the type of fracture. Most outcomes of cortical and endocortical volumetric bone mineral density comparisons were consistent with VBM results. Our results suggest: 1) that the spatial distribution of trabecular volumetric bone mineral density might play a significant role in hip fracture; 2) that focal cortical bone thinning might be more relevant in femoral neck fractures; and 3) that areas of reduced cortical and endocortical volumetric bone mineral density might be more relevant for trochanteric fractures in Chinese women. © 2017 American Society for Bone and Mineral Research.

Quantitative Computed Tomography (QCT) derived Bone Mineral Density (BMD) in finite element studies: a review of the literature

BACKGROUND

Finite element modeling of human bone provides a powerful tool to evaluate a wide variety of outcomes in a highly repeatable and parametric manner. These models are most often derived from computed tomography data, with mechanical properties related to bone mineral density (BMD) from the x-ray energy attenuation provided from this data. To increase accuracy, many researchers report the use of quantitative computed tomography (QCT), in which a calibration phantom is used during image acquisition to improve the estimation of BMD. Since model accuracy is dependent on the methods used in the calculation of BMD and density-mechanical property relationships, it is important to use relationships developed for the same anatomical location and using the same scanner settings, as these may impact model accuracy. The purpose of this literature review is to report the relationships used in the conversion of QCT equivalent density measures to ash, apparent, and/or tissue densities in recent finite element (FE) studies used in common density-modulus relationships. For studies reporting experimental validation, the validation metrics and results are presented.

RESULTS

Of the studies reviewed, 29% reported the use of a dipotassium phosphate (K2HPO4) phantom, 47% a hydroxyapatite (HA) phantom, 13% did not report phantom type, 7% reported use of both K2HPO4 and HA phantoms, and 4% alternate phantom types. Scanner type and/or settings were omitted or partially reported in 31% of studies. The majority of studies used densitometric and/or density-modulus relationships derived from different anatomical locations scanned in different scanners with different scanner settings. The methods used to derive various densitometric relationships are reported and recommendations are provided toward the standardization of reporting metrics.

CONCLUSIONS

This review assessed the current state of QCT-based FE modeling with use of clinical scanners. It was found that previously developed densitometric relationships vary by anatomical location, scanner type and settings. Reporting of all parameters used when referring to previously developed relationships, or in the development of new relationships, may increase the accuracy and repeatability of future FE models.

Statistical Parametric Mapping of HR-pQCT Images: A Tool for Population-Based Local Comparisons of Micro-Scale Bone Features

HR-pQCT enables in vivo multi-parametric assessments of bone microstructure in the distal radius and distal tibia. Conventional HR-pQCT image analysis approaches summarize bone parameters into global scalars, discarding relevant spatial information. In this work, we demonstrate the feasibility and reliability of statistical parametric mapping (SPM) techniques for HR-pQCT studies, which enable population-based local comparisons of bone properties. We present voxel-based morphometry (VBM) to assess trabecular and cortical bone voxel-based features, and a surface-based framework to assess cortical bone features both in cross-sectional and longitudinal studies. In addition, we present tensor-based morphometry (TBM) to assess trabecular and cortical bone structural changes. The SPM techniques were evaluated based on scan-rescan HR-pQCT acquisitions with repositioning of the distal radius and distal tibia of 30 subjects. For VBM and surface-based SPM purposes, all scans were spatially normalized to common radial and tibial templates, while for TBM purposes, rescans (follow-up) were spatially normalized to their corresponding scans (baseline). VBM was evaluated based on maps of local bone volume fraction (BV/TV), homogenized volumetric bone mineral density (vBMD), and homogenized strain energy density (SED) derived from micro-finite element analysis; while the cortical bone framework was evaluated based on surface maps of cortical bone thickness, vBMD, and SED. Voxel-wise and vertex-wise comparisons of bone features were done between the groups of baseline and follow-up scans. TBM was evaluated based on mean square errors of determinants of Jacobians at baseline bone voxels. In both anatomical sites, voxel- and vertex-wise uni- and multi-parametric comparisons yielded non-significant differences, and TBM showed no artefactual bone loss or apposition. The presented SPM techniques demonstrated robust specificity thus warranting their application in future clinical HR-pQCT studies.