Automated quantification of three-dimensional subject motion to monitor image quality in high-resolution peripheral quantitative computed tomography

Bone imaging modality precision and agreement between DXA, pQCT, and HR-pQCT

Quantifying precision error for DXA, peripheral QCT (pQCT), and HR-pQCT is crucial for monitoring longitudinal changes in body composition and musculoskeletal outcomes. Agreement and associations between bone variables assessed using pQCT and second-generation HR-pQCT are unclear. This study aimed to determine the precision of, and agreement and associations between, bone variables assessed via DXA, pQCT, and second-generation HR-pQCT. Thirty older adults (mean age: 64.2 ± 8.0 yr; women: 67%) were recruited. DXA scans were performed at the total hip, lumbar spine, and whole body. Distal (4%) and proximal (30%/33%/66%) skeletal sites at the radius and tibia were scanned with pQCT and/or HR-pQCT. Root-mean-squared coefficients of variation (%CVRMS) were calculated to define precision errors, and Bland-Altman plots assessed agreement between densitometric estimates. Pearson correlations and linear regression explored relationships between bone variables at different skeletal sites and proportional bias, respectively. Precision errors ranged between 0.55% and 1.6% for DXA, 0.40% and 4.8% for pQCT, and 0.13% and 30.7% for HR-pQCT. Systematic bias was identified between pQCT- and HR-pQCT-determined radius and tibia volumetric BMD (vBMD) estimates (all p<.001). Proportional bias was not observed between vBMD measures at any skeletal site (all p>.05). pQCT- and HR-pQCT-determined total, trabecular, and cortical vBMD and estimates of bone strength at the radius and tibia were strongly correlated (all p<.05). Precision error was low for most bone variables and within the expected range for all imaging modalities. We observed significant systematic bias, but no proportional bias, between pQCT- and second-generation HR-pQCT-determined vBMD estimates at the radius and tibia. Nevertheless, measures of bone density and strength were strongly correlated at all skeletal sites. These findings suggest that although bone density and strength estimates from both imaging modalities are not interchangeable, they are strongly related and likely have similar fracture prediction capabilities.

Deep Convolutional Neural Networks Provide Motion Grading for High-Resolution Peripheral Quantitative Computed Tomography of the Scaphoid

In vivo high-resolution peripheral quantitative computed tomography (HR-pQCT) studies on bone characteristics are limited, partly due to the lack of standardized and objective techniques to describe motion artifacts responsible for lower-quality images. This study investigates the ability of such deep-learning techniques to assess image quality in HR-pQCT datasets of human scaphoids. In total, 1451 stacks of 482 scaphoid images from 53 patients, each with up to six follow-ups within one year, and each with one non-displaced fractured and one contralateral intact scaphoid, were independently graded by three observers using a visual grading scale for motion artifacts. A 3D-CNN was used to assess image quality. The accuracy of the 3D-CNN to assess the image quality compared to the mean results of three skilled operators was between 92% and 96%. The 3D-CNN classifier reached an ROC-AUC score of 0.94. The average assessment time for one scaphoid was 2.5 s. This study demonstrates that a deep-learning approach for rating radiological image quality provides objective assessments of motion grading for the scaphoid with a high accuracy and a short assessment time. In the future, such a 3D-CNN approach can be used as a resource-saving and cost-effective tool to classify the image quality of HR-pQCT datasets in a reliable, reproducible and objective way.

Influence of demographic factors on the occurrence of motion artefacts in HR-pQCT

The study shows a high incidence of motion artefacts in a central European population and a significant increase of those artefacts with higher age. These findings may impact on the design and conduct of future in vivo HR-pQCT studies or at least help to estimate the potential number of drop outs due to unusable image quality.

PURPOSE

Motion artefacts in high-resolution peripheral quantitative computed tomography (HR-pQCT) are challenging, as they introduce error into the resulting measurement data. The aim of this study was to assess the general occurrence of motion artefacts in healthy distal radius and to evaluate the influence of demographic factors.

METHODS

The retrospective study is based on 525 distal radius second-generation HR-pQCT scans of 95 patients. All stacks were evaluated by two experienced observers and graded according to the visual grading scale recommended by the manufacturer, ranging from grade 1 (no visible motion artefacts) to grade 5 (severe motion artefacts). Correlations between demographic factors and image quality were evaluated using a linear mixed effects model analysis.

RESULTS

The average visual grading was 2.7 (SD ± 0.7). Age and severity of motion artefacts significantly correlated (p = 0.026). Patients aged 65 years or above had an average image quality between grades 1 and 3 in 72.7% of cases, while patients younger than 65 had an average image quality between grades 1 and 3 in 91.9% of cases. Gender, smoking behaviour, and handedness had no significant influence on motion artefacts.

CONCLUSION

This study showed a high incidence of motion artefacts in a representative central European population, but also a significant increase of motion artefacts with higher age. This could impact further study designs by planning for a sufficiently large and if possible a more selective study population to gain a representative amount of high-quality image data.

Multi-atlas segmentation and quantification of muscle, bone and subcutaneous adipose tissue in the lower leg using peripheral quantitative computed tomography

Accurate and reproducible tissue identification is essential for understanding structural and functional changes that may occur naturally with aging, or because of a chronic disease, or in response to intervention therapies. Peripheral quantitative computed tomography (pQCT) is regularly employed for body composition studies, especially for the structural and material properties of the bone. Furthermore, pQCT acquisition requires low radiation dose and the scanner is compact and portable. However, pQCT scans have limited spatial resolution and moderate SNR. pQCT image quality is frequently degraded by involuntary subject movement during image acquisition. These limitations may often compromise the accuracy of tissue quantification, and emphasize the need for automated and robust quantification methods. We propose a tissue identification and quantification methodology that addresses image quality limitations and artifacts, with increased interest in subject movement. We introduce a multi-atlas image segmentation (MAIS) framework for semantic segmentation of hard and soft tissues in pQCT scans at multiple levels of the lower leg. We describe the stages of statistical atlas generation, deformable registration and multi-tissue classifier fusion. We evaluated the performance of our methodology using multiple deformable registration approaches against reference tissue masks. We also evaluated the performance of conventional model-based segmentation against the same reference data to facilitate comparisons. We studied the effect of subject movement on tissue segmentation quality. We also applied the top performing method to a larger out-of-sample dataset and report the quantification results. The results show that multi-atlas image segmentation with diffeomorphic deformation and probabilistic label fusion produces very good quality over all tissues, even for scans with significant quality degradation. The application of our technique to the larger dataset reveals trends of age-related body composition changes that are consistent with the literature. Because of its robustness to subject motion artifacts, our MAIS methodology enables analysis of larger number of scans than conventional state-of-the-art methods. Automated analysis of both soft and hard tissues in pQCT is another contribution of this work.

Using 3D image registration to maximize the reproducibility of longitudinal bone strength assessment by HR-pQCT and finite element analysis

We developed and validated a finite element (FE) approach for longitudinal high-resolution peripheral quantitative computed tomography (HR-pQCT) studies using 3D image registration to account for misalignment between images. This reduced variability in longitudinal FE estimates and improved our ability to measure in vivo changes in HR-pQCT studies of bone strength.

INTRODUCTION

We developed and validated a finite element (FE) approach for longitudinal high-resolution peripheral quantitative computed tomography (HR-pQCT) studies using 3D rigid-body registration (3DR) to maximize reproducibility by accounting for misalignment between images.

METHODS

In our proposed approach, we used the full common bone volume defined by 3DR to estimate standard FE parameters. Using standard HR-pQCT imaging protocols, we validated the 3DR approach with ex vivo samples of the distal radius (n = 10, four repeat scans) by assessing whether 3DR can reduce measurement variability from repositioning error. We used in vivo data (n = 40, five longitudinal scans) to assess the sensitivity of 3DR to detect changes in bone strength at the distal radius by the standard deviation of the rate of change (σ), where the ideal value of σ is minimized to define true change. FE estimates by 3DR were compared to estimates by no registration (NR) and slice-matching (SM).

RESULTS

Group-wise comparisons of ex vivo variation (CV_RMS, %) found that FE measurement precision was improved by SM (CV_RMS < 0.80%) and 3DR (CV_RMS < 0.62%) compared to NR (CV_RMS~2%), and 3DR was advantageous as repositioning error increased. Longitudinal in vivo reproducibility was minimized by 3DR for failure load estimates (σ = 0.008 kN/month).

CONCLUSION

Although 3D registration cannot negate motion artifacts, it plays an important role in detecting and reducing variability in FE estimates for longitudinal HR-pQCT data and is well suited for estimating effects of interventions in in vivo longitudinal studies of bone strength.

Motion estimation and correction in SPECT, PET and CT

Patient motion impacts single photon emission computed tomography (SPECT), positron emission tomography (PET) and X-ray computed tomography (CT) by giving rise to projection data inconsistencies that can manifest as reconstruction artifacts, thereby degrading image quality and compromising accurate image interpretation and quantification. Methods to estimate and correct for patient motion in SPECT, PET and CT have attracted considerable research effort over several decades. The aims of this effort have been two-fold: to estimate relevant motion fields characterizing the various forms of voluntary and involuntary motion; and to apply these motion fields within a modified reconstruction framework to obtain motion-corrected images. The aims of this review are to outline the motion problem in medical imaging and to critically review published methods for estimating and correcting for the relevant motion fields in clinical and preclinical SPECT, PET and CT. Despite many similarities in how motion is handled between these modalities, utility and applications vary based on differences in temporal and spatial resolution. Technical feasibility has been demonstrated in each modality for both rigid and non-rigid motion, but clinical feasibility remains an important target. There is considerable scope for further developments in motion estimation and correction, and particularly in data-driven methods that will aid clinical utility. State-of-the-art machine learning methods may have a unique role to play in this context.

Longitudinal bone microarchitectural changes are best detected using image registration

Longitudinal studies of bone using high-resolution medical imaging may result in non-physiological measurements of longitudinal changes. In this study, we determined that three-dimensional image processing techniques best capture realistic longitudinal changes in bone density and should therefore be used with high-resolution imaging when studying bone changes over time.

INTRODUCTION

The purpose of this study was to determine which longitudinal analysis technique (no registration (NR), slice-match (SM) registration, or three-dimensional registration (3DR)) produced the most realistic longitudinal changes in a 3-year study of bone density and structure using high-resolution peripheral quantitative computed tomography (HR-pQCT).

METHODS

We assessed HR-pQCT scans of the distal radius and tibia for men and women (N = 40) aged 55-70 years at baseline and 6, 12, 24, and 36 months. To evaluate which longitudinal analysis technique (NR, SM, or 3DR) best captured physiologically reasonable 3-year changes, we calculated the standard deviation of the absolute rate of change in each bone parameter. The data were compared between longitudinal analysis techniques using repeated measures ANOVA and post hoc analysis.

RESULTS

As expected, both SM and 3DR better captured physiological longitudinal changes than NR. At the tibia, there were no differences between SM and 3DR; however, at the radius where precision was lower, 3DR produced better results for total bone mineral density.

CONCLUSIONS

At least SM or 3DR should be implemented in longitudinal studies using HR-pQCT. 3DR is preferable, particularly at the radius, to ensure that physiological changes in bone density are observed.

Deficits in bone strength, density and microarchitecture in women living with HIV: A cross-sectional HR-pQCT study

PURPOSE

With the advent of combined antiretroviral therapy (cART), life expectancy has increased among persons living with HIV, but so too has risk for comorbidities including osteoporosis and fragility fracture. To explore whether HIV status and cART influence three-dimensional measures of BMD, bone microarchitecture and strength we aimed to compare these outcomes between women living with HIV (WLWH; n = 50; 50.4 ± 1.2 years, 44% postmenopausal) and without HIV (controls; n = 50; 51.8 ± 1.2 years, 52% postmenopausal).

METHODS

Outcomes were lumbar spine, total hip and femoral neck areal BMD by DXA; distal radius and tibia trabecular BMD, thickness and number, and cortical BMD and area by HR-pQCT; and finite element analysis-derived bone strength (failure load). Multivariable regression analysis compared bone outcomes between groups adjusting for known osteoporosis risk factors. Within WLWH, we examined associations between bone outcomes and HIV-related factors including disease severity and cART duration.

RESULTS

WLWH were diagnosed 20 ± 4 years ago, were on cART for 123 ± 37 months and 80% had HIV plasma viral load <40 copies/mL. For women ≥50 years (n = 61), total hip aBMD T-Score was lower among WLWH than controls. Adjusted distal radius trabecular BMD and thickness and distal tibia trabecular BMD and failure load were 8-19% lower in WLWH than controls (p < 0.05). Cortical BMD and area did not differ between groups at either site. Lifetime cART duration and current plasma viral load were not associated with bone outcomes in WLWH; however, previous treatment with tenofovir was negatively associated with distal radius trabecular BMD and trabecular number and LS aBMD T-score.

CONCLUSIONS

WLWH have compromised BMD, bone microarchitecture and strength vs. controls of similar age and reproductive status. Treatment with tenofovir may contribute to bone deficits in WLWH.

Please Don't Move-Evaluating Motion Artifact From Peripheral Quantitative Computed Tomography Scans Using Textural Features

Most imaging methods, including peripheral quantitative computed tomography (pQCT), are susceptible to motion artifacts particularly in fidgety pediatric populations. Methods currently used to address motion artifact include manual screening (visual inspection) and objective assessments of the scans. However, previously reported objective methods either cannot be applied on the reconstructed image or have not been tested for distal bone sites. Therefore, the purpose of the present study was to develop and validate motion artifact classifiers to quantify motion artifact in pQCT scans. Whether textural features could provide adequate motion artifact classification performance in 2 adolescent datasets with pQCT scans from tibial and radial diaphyses and epiphyses was tested. The first dataset was split into training (66% of sample) and validation (33% of sample) datasets. Visual classification was used as the ground truth. Moderate to substantial classification performance (J48 classifier, kappa coefficients from 0.57 to 0.80) was observed in the validation dataset with the novel texture-based classifier. In applying the same classifier to the second cross-sectional dataset, a slight-to-fair (κ = 0.01-0.39) classification performance was observed. Overall, this novel textural analysis-based classifier provided a moderate-to-substantial classification of motion artifact when the classifier was specifically trained for the measurement device and population. Classification based on textural features may be used to prescreen obviously acceptable and unacceptable scans, with a subsequent human-operated visual classification of any remaining scans.

Least-detectable and age-related local in vivo bone remodelling assessed by time-lapse HR-pQCT

We previously developed an image analysis approach for the determination of local sites of bone remodelling using time-lapse in vivo HR-pQCT. The involved image filtering for removing noise was chosen rather aggressively, and also removed some effects of the bone remodelling. In this paper, we quantify these filtering settings using ex vivo reproducibility HR-pQCT images, and determine the least-detectable bone remodelling using in vivo reproducibility HR-pQCT images, as well as testing whether the approach is capable of capturing age-related bone remodelling by use of in vivo long-term HR-pQCT images. We found that a threshold value of 225 mg HA/cm³ for the filtering led to acceptable results with falsely determined bone remodelling of less than 0.5%, and that the least-detectable bone formation and bone resorption are 2.0 ± 1.0% and 2.2 ± 0.7% respectively. We also found that age-related local bone remodelling can be captured satisfactorily in postmenopausal women. The latter revealed new insights into the effect of ageing on bone remodelling, and showed that bone remodelling seems to take place through a few small formation packets and many large resorption volumes leading to a net bone loss. We conclude that local in vivo bone remodelling can be successfully assessed with time-lapse in vivo HR-pQCT capable of assessing age-related changes in bone remodelling.

Microarchitecture and Peripheral BMD are Impaired in Postmenopausal White Women With Fracture Independently of Total Hip T-Score: An International Multicenter Study

Because single-center studies have reported conflicting associations between microarchitecture and fracture prevalence, we included high-resolution peripheral quantitative computed tomography (HR-pQCT) data from five centers worldwide into a large multicenter analysis of postmenopausal women with and without fracture. Volumetric BMD (vBMD) and microarchitecture were assessed at the distal radius and tibia in 1379 white postmenopausal women (age 67 ± 8 years); 470 (34%) had at least one fracture including 349 with a major fragility fracture. Age, height, weight, and total hip T-score differed across centers and were employed as covariates in analyses. Women with fracture had higher BMI, were older, and had lower total hip T-score, but lumbar spine T-score was similar between groups. At the radius, total and trabecular vBMD and cortical thickness were significantly lower in fractured women in three out of five centers, and trabecular number in two centers. Similar results were found at the tibia. When data from five centers were combined, however, women with fracture had significantly lower total, trabecular, and cortical vBMD (2% to 7%), lower trabecular number (4% to 5%), and thinner cortices (5% to 6%) than women without fracture after adjustment for covariates. Results were similar at the radius and tibia. Similar results were observed with analysis restricted to major fragility fracture, vertebral and hip fractures, and peripheral fracture (at the radius). When focusing on osteopenic women, each SD decrease of total and trabecular vBMD was associated with a significantly increased risk of major fragility fracture (OR = 1.55 to 1.88, p < 0.01) after adjustment for covariates. Moreover, trabecular architecture modestly improved fracture discrimination beyond peripheral total vBMD. In conclusion, we observed differences by center in the magnitude of fracture/nonfracture differences at both the distal radius and tibia. However, when data were pooled across centers and the sample size increased, we observed significant and consistent deficits in vBMD and microarchitecture independent of total hip T-score in all postmenopausal white women with fracture and in the subgroup of osteopenic women, compared to women who never had a fracture. © 2016 American Society for Bone and Mineral Research.

Challenges in longitudinal measurements with HR-pQCT: evaluation of a 3D registration method to improve bone microarchitecture and strength measurement reproducibility

Definition of identical regions between repeated computed tomography (CT) scans is a key factor to monitor changes in bone microarchitecture. In longitudinal studies, accurate determination of the volume of interest (VOI), using three dimensional (3D) registration may improve precision. Therefore, the aim of our study was to investigate the short-term reproducibility of bone geometry, density, microstructure and biomechanical parameters assessed by HR-pQCT and micro-finite element (μFE) derived analyses, using the cross-sectional area (CSA) registration method in comparison with the use of 3D registration, to find overlapping regions between scans. Fifteen healthy individuals (aged 21-47 years) underwent 3 separate scans at the distal radius and tibia, within a one-month interval. Reproducibility was assessed after double contouring the cortical compartment and after applying three different methods to determine the common region between repeated scans: (i) the VOI was determined with no registration, i.e., on 110 slices, (ii) the VOI was determined after CSA-based registration, and (iii) the VOI was determined after 3D registration. Both pre- and post-registration short-term reproducibility for each subject was determined. With no registration, CVrms of geometry parameters ranged from 0.5 to 3.7%, showing a slight variation in the CSA between scans. When the CSA registration method was employed, the variability of geometry (CVrms<1.8%) and density parameters (CVrms<1.8%), was better than that obtained without registration. By removing the effect of repositioning, the 3D registration further improved the reproducibility of cortical bone measurements compared to other methods. Indeed, significant improvements were found for cortical geometry and microstructure measurements (CVrms ranged from 0.4% to 10.7% at both sites; p<0.05), whereas the impact on trabecular bone measurements was restricted to its geometry parameter. The repositioning error was significantly reduced, most markedly at the radius compared to the tibia. For μFE measures, the impact of 3D registration on whole bone stiffness was negligible, indicating adequate assessment of longitudinal changes in estimated biomechanical properties, even without registration. In conclusion, we have shown that the 3D registration improved the identification of the common region retained for longitudinal analysis, contributing to improve the reproducibility of cortical bone parameter measurements. We also quantified the minimally detectable bone changes to help designing future studies with HR-pQCT.

Quantitative computer tomography in children and adolescents: the 2013 ISCD Pediatric Official Positions

In 2007, International Society of Clinical Densitometry Pediatric Positions Task Forces reviewed the evidence for the clinical application of peripheral quantitative computed tomography (pQCT) in children and adolescents. At that time, numerous limitations regarding the clinical application of pQCT were identified, although its use as a research modality for investigation of bone strength was highlighted. The present report provides an updated review of evidence for the clinical application of pQCT, as well as additional reviews of whole body QCT scans of the central and peripheral skeletons, and high-resolution pQCT in children. Although these techniques remain in the domain of research, this report summarizes the recent literature and evidence of the clinical applicability and offers general recommendations regarding the use of these modalities in pediatric bone health assessment.

Standardizing evaluation of pQCT image quality in the presence of subject movement: qualitative versus quantitative assessment

Peripheral quantitative computed tomography (pQCT) is an essential tool for assessing bone parameters of the limbs, but subject movement and its impact on image quality remains a challenge to manage. The current approach to determine image viability is by visual inspection, but pQCT lacks a quantitative evaluation. Therefore, the aims of this study were to (1) examine the reliability of a qualitative visual inspection scale and (2) establish a quantitative motion assessment methodology. Scans were performed on 506 healthy girls (9-13 years) at diaphyseal regions of the femur and tibia. Scans were rated for movement independently by three technicians using a linear, nominal scale. Quantitatively, a ratio of movement to limb size (%Move) provided a measure of movement artifact. A repeat-scan subsample (n = 46) was examined to determine %Move's impact on bone parameters. Agreement between measurers was strong (intraclass correlation coefficient = 0.732 for tibia, 0.812 for femur), but greater variability was observed in scans rated 3 or 4, the delineation between repeat and no repeat. The quantitative approach found ≥95% of subjects had %Move <25 %. Comparison of initial and repeat scans by groups above and below 25% initial movement showed significant differences in the >25 % grouping. A pQCT visual inspection scale can be a reliable metric of image quality, but technicians may periodically mischaracterize subject motion. The presented quantitative methodology yields more consistent movement assessment and could unify procedure across laboratories. Data suggest a delineation of 25% movement for determining whether a diaphyseal scan is viable or requires repeat.

Advanced CT based in vivo methods for the assessment of bone density, structure, and strength

Based on spiral 3D tomography a large variety of applications have been developed during the last decade to asses bone mineral density, bone macro and micro structure, and bone strength. Quantitative computed tomography (QCT) using clinical whole body scanners provides separate assessment of trabecular, cortical, and subcortical bone mineral density (BMD) and content (BMC) principally in the spine and hip, although the distal forearm can also be assessed. Further bone macrostructure, for example bone geometry or cortical thickness can be quantified. Special high resolution peripheral CT (hr-pQCT) devices have been introduced to measure bone microstructure for example the trabecular architecture or cortical porosity at the distal forearm or tibia. 3D CT is also the basis for finite element analysis (FEA) to determine bone strength. QCT, hr-pQCT, and FEM are increasingly used in research as well as in clinical trials to complement areal BMD measurements obtained by the standard densitometric technique of dual x-ray absorptiometry (DXA). This review explains technical developments and demonstrates how QCT based techniques advanced our understanding of bone biology.

High-resolution peripheral quantitative computed tomography for the assessment of bone strength and structure: a review by the Canadian Bone Strength Working Group

Bone structure is an integral determinant of bone strength. The availability of high resolution peripheral quantitative computed tomography (HR-pQCT) has made it possible to measure three-dimensional bone microarchitecture and volumetric bone mineral density in vivo, with accuracy previously unachievable and with relatively low-dose radiation. Recent studies using this novel imaging tool have increased our understanding of age-related changes and sex differences in bone microarchitecture, as well as the effect of different pharmacological therapies. One advantage of this novel tool is the use of finite element analysis modelling to non-invasively estimate bone strength and predict fractures using reconstructed three-dimensional images. In this paper, we describe the strengths and limitations of HR-pQCT and review the clinical studies using this tool.

Multicenter precision of cortical and trabecular bone quality measures assessed by high-resolution peripheral quantitative computed tomography

High-resolution peripheral quantitative computed tomography (HR-pQCT) has recently been introduced as a clinical research tool for in vivo assessment of bone quality. The utility of this technology to address important skeletal health questions requires translation to standardized multicenter data pools. Our goal was to evaluate the feasibility of pooling data in multicenter HR-pQCT imaging trials. Reproducibility imaging experiments were performed using structure and composition-realistic phantoms constructed from cadaveric radii. Single-center precision was determined by repeat scanning over short-term (<72 hours), intermediate-term (3-5 months), and long-term intervals (28 months). Multicenter precision was determined by imaging the phantoms at nine different HR-pQCT centers. Least significant change (LSC) and root mean squared coefficient of variation (RMSCV) for each interval and across centers was calculated for bone density, geometry, microstructure, and biomechanical parameters. Single-center short-term RMSCVs were <1% for all parameters except cortical thickness (Ct.Th) (1.1%), spatial variability in cortical thickness (Ct.Th.SD) (2.6%), standard deviation of trabecular separation (Tb.Sp.SD) (1.8%), and porosity measures (6% to 8%). Intermediate-term RMSCVs were generally not statistically different from short-term values. Long-term variability was significantly greater for all density measures (0.7% to 2.0%; p < 0.05 versus short-term) and several structure measures: cortical thickness (Ct.Th) (3.4%; p < 0.01 versus short-term), cortical porosity (Ct.Po) (15.4%; p < 0.01 versus short-term), and trabecular thickness (Tb.Th) (2.2%; p < 0.01 versus short-term). Multicenter RMSCVs were also significantly higher than short-term values: 2% to 4% for density and micro-finite element analysis (µFE) measures (p < 0.0001), 2.6% to 5.3% for morphometric measures (p < 0.001), whereas Ct.Po was 16.2% (p < 0.001). In the absence of subject motion, multicenter precision errors for HR-pQCT parameters were generally less than 5%. Phantom-based multicenter precision was comparable to previously reported in in vivo single-center precision errors, although this was approximately two to five times worse than ex vivo short-term precision. The data generated from this study will contribute to the future design and validation of standardized procedures that are broadly translatable to multicenter study designs.

Automated threshold-independent cortex segmentation by 3D-texture analysis of HR-pQCT scans

The quantitative assessment of metabolic bone diseases relies on tissue properties such as bone mineral density (BMD) and bone microarchitecture. In spite of an increasing number of publications using high-resolution peripheral quantitative computed-tomography (HR-pQCT), the accurate and reproducible separation of cortical and trabecular bone remains challenging. In this paper, we present a novel, fully automated, threshold-independent technique for the segmentation of cortical and trabecular bone in HR-pQCT scans. This novel post-processing method is based on modeling appearance characteristics from manually annotated cases. In our experiments the algorithm automatically selected texture features with high differentiating power and trained a classifier to separate cortical and trabecular bone. From this mask, cortical thickness and tissue volume could be calculated with high accuracy. The overlap between the proposed threshold-independent segmentation tool (TIST) and manual contouring was 0.904±0.045 (Dice coefficient). In our experiments, TIST obtained higher overall accuracy in our measurements than other techniques.

Automated motion correction for in vivo optical projection tomography

In in vivo optical projection tomography (OPT), object motion will significantly reduce the quality and resolution of the reconstructed image. Based on the well-known Helgason-Ludwig consistency condition (HLCC), we propose a novel method for motion correction in OPT under parallel beam illumination. The method estimates object motion from projection data directly and does not require any other additional information, which results in a straightforward implementation. We decompose object movement into translation and rotation, and discuss how to correct for both translation and general motion simultaneously. Since finding the center of rotation accurately is critical in OPT, we also point out that the system's geometrical offset can be considered as object translation and therefore also calibrated through the translation estimation method. In order to verify the algorithm effectiveness, both simulated and in vivo OPT experiments are performed. Our results demonstrate that the proposed approach is capable of decreasing movement artifacts significantly thus providing high quality reconstructed images in the presence of object motion.

Quality control for bone quality parameters affected by subject motion in high-resolution peripheral quantitative computed tomography

Subject motion during high-resolution peripheral quantitative computed tomography (HR-pQCT) causes image artifacts that affect morphological analysis of bone quality. The aim of our study was to determine effectiveness of techniques for quality control in the presence of motion in vivo including automated and manual approaches. First, repeatability of manual grading was determined within and between laboratories. Given proper training using a standardized scale and training images (provided by the manufacturer), we found that manual grading is suitable for repeatable image quality grading within and across sites (ICC>0.7). Both a new automated technique providing motion measures based on projection moments, and traditional manual grading (1=best, 5=worst) were subsequently used to assess subject data for motion in N=137 image pairs (scan/re-scan) from the Canadian Multicentre Osteoporosis Study (CaMos) Calgary cohort. High quality image pairs were selected and measurement precision was estimated by calculating the coefficient of variation (CV). Consistent with previous data, density parameters (e.g. total bone mineral density) are more robust than structural (e.g. trabecular number) or finite element parameters (e.g. failure load). To obtain acceptable measurement precision, images should not exceed a manual grading of 3 (on a scale from 1 to 5) or an automatic (ε(T)) grading of 1.2. Automatic and manual grading provide comparable quality control, but the advantage of the automated technique is its ability to provide a motion value at scan time (providing a basis for real time decision regarding re-scan requirements), and the assessment is objective. Notably, automatic motion measurement can be performed retrospectively based on original scan data, and is therefore well suited for multi-center studies as well as any research where objective quality control is paramount.

Noninvasive imaging of bone microarchitecture

The noninvasive quantification of peripheral compartment-specific bone microarchitecture is feasible with high-resolution peripheral quantitative computed tomography (HR-pQCT) and high-resolution magnetic resonance imaging (HR-MRI). In addition to classic morphometric indices, both techniques provide a suitable basis for virtual biomechanical testing using finite element (FE) analyses. Methodical limitations, morphometric parameter definition, and motion artifacts have to be considered to achieve optimal data interpretation from imaging studies. With increasing availability of in vivo high-resolution bone imaging techniques, special emphasis should be put on quality control including multicenter, cross-site validations. Importantly, conclusions from interventional studies investigating the effects of antiosteoporotic drugs on bone microarchitecture should be drawn with care, ideally involving imaging scientists, translational researchers, and clinicians.