QCT of the femur: Comparison between QCTPro CTXA and MIAF Femur.

QCT is commonly employed in research studies and clinical trials to measure BMD at the proximal femur. In this study we compared two analysis software options, QCTPro CTXA and MIAF-Femur, using CT scans of the semi-anthropometric European Proximal Femur Phantom (EPFP) and in vivo data from 130 Chinese elderly men and women aged 60-80 years. Integral (Int), cortical (Cort) and trabecular (Trab) vBMD, volume, and BMC of the neck (FN), trochanter (TR), inter-trochanter (IT), and total hip (TH) VOIs were compared. Accuracy was determined in the 5 mm wide central portion of the femoral neck of the EPFP. Nominal values were: cross-sectional area (CSA) 4.9 cm2, cortical thickness (C.Th) 2 mm, CortBMD 723 mg/cm3 and TrabBMD 100 mg/cm3. In MIAF the so-called peeled trabecular VOI was analyzed, which excludes subcortical bone to avoid partial volume artefacts at the endocortical border that artificially increase TrabBMD. For CTXA uncorrected, so called raw cortical values were used for the analysis. QCTPro and MIAF phantom results were: CSA 5.9 cm2 versus 5.1 cm2; C.Th 1.68 mm versus 1.92 mm; CortBMD 578 mg/cm3 versus 569 mg/cm3; and TrabBMD 154 mg/cm3 versus 104 mg/cm3. In vivo correlations (R2) of integral and trabecular bone parameters ranged from 0.63 to 0.96. Bland-Altman analysis for TH and FN TrabBMD showed that lower mean values were associated with higher differences, which means that TrabBMD differences between MIAF and CTXA are larger for osteoporotic than for normal patients, which can be largely explained by the inclusion of subcortical BMD in the trabecular VOI analyzed by CTXA in combination with fixed thresholds used to separate cortical from trabecular bone compartments. Correlations between CTXA corrected CortBMD and MIAF were negative, whereas raw data correlated positively with MIAF measurements for all VOIs questioning the validity of the CTXA corrections. The EPFP results demonstrated higher MIAF accuracy of cortical thickness and TrabBMD. Integral and trabecular bone parameters were highly correlated between CTXA and MIAF. Partial volume artefacts at the endocortical border artificially increased trabecular BMD by CTXA, especially for osteoporosis patients. With respect to volumetric cortical measurements with CTXA, the use raw data is recommended, because corrected data cause a negative correlation with MIAF CortBMD.

Three-dimensional Distribution of Muscle and Adipose Tissue of the Thigh at CT: Association with Acute Hip Fracture.

Purpose To evaluate determinants of hip fracture by assessing soft-tissue composition of the upper thigh at CT. Materials and Methods In this retrospective analysis of prospectively collected data, CT studies in 55 female control participants (mean age, 73.1 years ± 9.3 [standard deviation]) were compared with those in 40 female patients (mean age, 80.2 years ± 11.0) with acute hip fractures. Eighty-seven descriptors of the soft-tissue composition were determined. A multivariable best subsets analysis was used to extract parameters best associated with hip fracture. Results were adjusted for age, height, and weight. Results of soft-tissue parameters were compared with bone mineral density (BMD) and cortical bone thickness. Areas under the receiver operating characteristic curve (AUCs) adjusted for multiple comparisons were determined to discriminate fracture. Results The hip fracture group was characterized by lower BMD, lower cortical thickness, lower relative adipose tissue volume of the upper thigh, and higher extramyocellular lipid (EML) surface density. The relative volume of adipose tissue combined with EML surface density (model S1) was associated with hip fracture (AUC, 0.85; 95% confidence interval [CI]: 0.78, 0.93), as well as trochanteric trabecular BMD combined with neck cortical thickness (model B2) (AUC, 0.84; 95% CI: 0.75, 0.92). The model including all four parameters provided significantly better (P < .01) discrimination (AUC, 0.92; 95% CI: 0.86, 0.97) than model S1 or B2. Conclusion In addition to bone mineral density and geometry of the proximal femur, the amount of adipose tissue of the upper thigh and the distribution of the adipocytes in the muscles are significantly associated with acute hip fracture at CT. © RSNA, 2018 Online supplemental material is available for this article.

Quantitative analysis of skeletal muscle by computed tomography imaging-State of the art.

The radiological assessment of muscle properties-size, mass, density (also termed radiodensity), composition, and adipose tissue infiltration-is fundamental in muscle diseases. More recently, it also became obvious that muscle atrophy, also termed muscle wasting, is caused by or associated with many other diseases or conditions, such as inactivity, malnutrition, chronic obstructive pulmonary disorder, cancer-associated cachexia, diabetes, renal and cardiac failure, and sarcopenia and even potentially with osteoporotic hip fracture. Several techniques have been developed to quantify muscle morphology and function. This review is dedicated to quantitative computed tomography (CT) of skeletal muscle and only includes a brief comparison with magnetic resonance imaging. Strengths and limitations of CT techniques are discussed in detail, including CT scanner calibration, acquisition and reconstruction protocols, and the various quantitative parameters that can be measured with CT, starting from simple volume measures to advanced parameters describing the adipose tissue distribution within muscle. Finally, the use of CT in sarcopenia and cachexia and the relevance of muscle parameters for the assessment of osteoporotic fracture illustrate the application of CT in two emerging areas of medical interest.

Methods for segmentation of rheumatoid arthritis bone erosions in high-resolution peripheral quantitative computed tomography (HR-pQCT).

OBJECTIVE: The comparison between different techniques to quantify the 3-dimensional size of inflammatory bone erosions in rheumatoid arthritis(RA) patients. METHODS: Anti-cyclic citrullinated peptide antibody(ACPA) positive RA patients received high-resolution peripheral quantitative computed tomography (HR-pQCT) scans of the metacarpophalangeal joints (MCP). Erosions were measured by three different segmentation techniques: (1) manual method with calculation by half-ellipsoid formula, (2) semi-automated modified Evaluation Script for Erosions (mESE), and (3) semi-automated Medical Image Analysis Framework (MIAF) software. Bland & Altman plots were used to describe agreement between methods. Furthermore, shape of erosions was classified as regular or irregular and then compared to the sphericity obtained by MIAF. RESULTS: A total of 76 erosions from 65 RA patients (46 females/19 males), median age 57 years, median disease duration 6.1 years and median disease activity score 28 of 2.8 units were analyzed. While mESE and MIAF showed good agreement in the measurement of erosion size, the manual method with calculation by half-ellipsoid formula underestimated erosions size, particularly with larger erosions. Accurate segmentation is particularly important in larger erosions, which are irregularly shaped. In all three segmentation techniques irregular erosions were larger in size than regular erosions (MIAF: 19.7 vs. 3.4mm3; mESE: 15.5 vs. 2.3mm3; manual = 7.2 vs. 1.52mm3; all p < 0.001). In accordance, sphericity of erosions measured by MIAF significantly decreased with their size (p < 0.001). CONCLUSION: MIAF and mESE allow segmentation of inflammatory bone erosions in RA patients with excellent inter reader reliability. They allow calculating erosion volume independent of erosion shape and therefore provide an attractive tool to quantify structural damage in individual joints of RA patients.

A new method to determine cortical bone thickness in CT images using a hybrid approach of parametric profile representation and local adaptive thresholds: Accuracy results.

MOTIVATION: Cortical bone is an important contributor to bone strength and is pivotal to understand the etiology of osteoporotic fractures and the specific mechanisms of antiosteoporotic treatment regimen. 3D computed tomography (CT) can be used to measure cortical thickness, density, and mass in the proximal femur, lumbar vertebrae, and distal forearm. However, the spatial resolution of clinical whole body CT scanners is limited by radiation exposure; partial volume artefacts severely impair the accurate assessment of cortical parameters, in particular in locations where the cortex is thin such as in the lumbar vertebral bodies or in the femoral neck. METHOD: Model-based deconvolution approaches recover the cortical thickness by numerically deconvolving the image along 1D profiles using an estimated scanner point spread function (PSF) and a hypothesized uniform cortical bone mineral density (reference density). In this work we provide a new essentially analytical unique solution to the model-based cortex recovery problem using few characteristics of the measured profile and thus eliminate the non-linear optimization step for deconvolution. Also, the proposed approach allows to get rid of the PSF in the model and reduces sensitivity to errors in the reference density. Additionally, run-time and memory effective computation of cortical thickness was achieved with the help of a lookup table. RESULTS: The method accuracy and robustness was validated and compared to that of a deconvolution approach recently proposed for cortical bone and of the 50% relative threshold technique: in a simulated environment with noise and various error levels in the reference density and using CT acquisitions of the European Forearm Phantom (EFP II), a modification of a widely used anthropomorphic standard of cortical and trabecular bone compartments that was scanned with various scan protocols. CONCLUSION: Results of simulations and of phantom data analysis verified the following properties of the new method: 1) Robustness against errors in the reference density. 2) Excellent accuracy on ground truth data with various noise levels. 3) Very fast computation using a lookup table.

A reproducible semi-automatic method to quantify the muscle-lipid distribution in clinical 3D CT images of the thigh.

Many studies use threshold-based techniques to assess in vivo the muscle, bone and adipose tissue distribution of the legs using computed tomography (CT) imaging. More advanced techniques divide the legs into subcutaneous adipose tissue (SAT), anatomical muscle (muscle tissue and adipocytes within the muscle border) and intra- and perimuscular adipose tissue. In addition, a so-called muscle density directly derived from the CT-values is often measured. We introduce a new integrated approach to quantify the muscle-lipid system (MLS) using quantitative CT in patients with sarcopenia or osteoporosis. The analysis targets the thigh as many CT studies of the hip do not include entire legs The framework consists of an anatomic coordinate system, allowing delineation of reproducible volumes of interest, a robust semi-automatic 3D segmentation of the fascia and a comprehensive method to quantify of the muscle and lipid distribution within the fascia. CT density-dependent features are calibrated using subject-specific internal CT values of the SAT and external CT values of an in scan calibration phantom. Robustness of the framework with respect to operator interaction, image noise and calibration was evaluated. Specifically, the impact of inter- and intra-operator reanalysis precision and addition of Gaussian noise to simulate lower radiation exposure on muscle and AT volumes, muscle density and 3D texture features quantifying MLS within the fascia, were analyzed. Existing data of 25 subjects (age: 75.6 ± 8.7) with porous and low-contrast muscle structures were included in the analysis. Intra- and inter-operator reanalysis precision errors were below 1% and mostly comparable to 1% of cohort variation of the corresponding features. Doubling the noise changed most 3D texture features by up to 15% of the cohort variation but did not affect density and volume measurements. The application of the novel technique is easy with acceptable processing time. It can thus be employed for a comprehensive quantification of the muscle-lipid system enabling radiomics approaches to musculoskeletal disorders.

Advanced Knee Structure Analysis (AKSA): a comparison of bone mineral density and trabecular texture measurements using computed tomography and high-resolution peripheral quantitative computed tomography of human knee cadavers.

BACKGROUND: A change of loading conditions in the knee causes changes in the subchondral bone and may be a cause of osteoarthritis (OA). However, quantification of trabecular architecture in vivo is difficult due to the limiting spatial resolution of the imaging equipment; one approach is the use of texture parameters. In previous studies, we have used digital models to simulate changes of subchondral bone architecture under OA progression. One major result was that, using computed tomography (CT) images, subchondral bone mineral density (BMD) in combination with anisotropy and global homogeneity could characterize this progression. The primary goal of this study was a comparison of BMD, entropy, anisotropy, variogram slope, and local and global inhomogeneity measurements between high-resolution peripheral quantitative CT (HR-pQCT) and CT using human cadaveric knees. The secondary goal was the verification of the spatial resolution dependence of texture parameters observed in the earlier simulations, two important prerequisites for the interpretation of in vivo measurements in OA patients. METHOD: The applicability of texture analysis to characterize bone architecture in clinical CT examinations was investigated and compared to results obtained from HR-pQCT. Fifty-seven human knee cadavers (OA status unknown) were examined with both imaging modalities. Three-dimensional (3D) segmentation and registration processes, together with automatic positioning of 3D analysis volumes of interest (VOIs), ensured the measurement of BMD and texture parameters at the same anatomical locations in CT and HR-pQCT datasets. RESULTS: According to the calculation of dice ratios (>0.978), the accuracy of VOI locations between methods was excellent. Entropy, anisotropy, and global inhomogeneity showed significant and high linear correlation between both methods (0.68 < R 2 < 1.00). The resolution dependence of these parameters simulated earlier was confirmed by the in vitro measurements. CONCLUSION: The high correlation of HR-pQCT- and CT-based measurements of entropy, global inhomogeneity, and anisotropy suggests interchangeability between devices regarding the quantification of texture. The agreement of the experimentally determined resolution dependence of global inhomogeneity and anisotropy with earlier simulations is an important milestone towards their use to quantify subchondral bone structure. However, an in vivo study is still required to establish their clinical relevance.

Prediction of Hip Failure Load: In Vitro Study of 80 Femurs Using Three Imaging Methods and Finite Element Models-The European Fracture Study (EFFECT).

Purpose To evaluate the performance of three imaging methods (radiography, dual-energy x-ray absorptiometry [DXA], and quantitative computed tomography [CT]) and that of a numerical analysis with finite element modeling (FEM) in the prediction of failure load of the proximal femur and to identify the best densitometric or geometric predictors of hip failure load. Materials and Methods Institutional review board approval was obtained. A total of 40 pairs of excised cadaver femurs (mean patient age at time of death, 82 years ± 12 [standard deviation]) were examined with (a) radiography to measure geometric parameters (lengths, angles, and cortical thicknesses), (b) DXA (reference standard) to determine areal bone mineral densities (BMDs), and (c) quantitative CT with dedicated three-dimensional analysis software to determine volumetric BMDs and geometric parameters (neck axis length, cortical thicknesses, volumes, and moments of inertia), and (d) quantitative CT-based FEM to calculate a numerical value of failure load. The 80 femurs were fractured via mechanical testing, with random assignment of one femur from each pair to the single-limb stance configuration (hereafter, stance configuration) and assignment of the paired femur to the sideways fall configuration (hereafter, side configuration). Descriptive statistics, univariate correlations, and stepwise regression models were obtained for each imaging method and for FEM to enable us to predict failure load in both configurations. Results Statistics reported are for stance and side configurations, respectively. For radiography, the strongest correlation with mechanical failure load was obtained by using a geometric parameter combined with a cortical thickness (r(2) = 0.66, P < .001; r(2) = 0.65, P < .001). For DXA, the strongest correlation with mechanical failure load was obtained by using total BMD (r(2) = 0.73, P < .001) and trochanteric BMD (r(2) = 0.80, P < .001). For quantitative CT, in both configurations, the best model combined volumetric BMD and a moment of inertia (r(2) = 0.78, P < .001; r(2) = 0.85, P < .001). FEM explained 87% (P < .001) and 83% (P < .001) of bone strength, respectively. By combining (a) radiography and DXA and (b) quantitative CT and DXA, correlations with mechanical failure load increased to 0.82 (P < .001) and 0.84 (P < .001), respectively, for radiography and DXA and to 0.80 (P < .001) and 0.86 (P < .001) , respectively, for quantitative CT and DXA. Conclusion Quantitative CT-based FEM was the best method with which to predict the experimental failure load; however, combining quantitative CT and DXA yielded a performance as good as that attained with FEM. The quantitative CT DXA combination may be easier to use in fracture prediction, provided standardized software is developed. These findings also highlight the major influence on femoral failure load, particularly in the trochanteric region, of a densitometric parameter combined with a geometric parameter. (©) RSNA, 2016 Online supplemental material is available for this article.

Comparison of proximal femur and vertebral body strength improvements in the FREEDOM trial using an alternative finite element methodology.

Denosumab reduced the incidence of new fractures in postmenopausal women with osteoporosis by 68% at the spine and 40% at the hip over 36 months compared with placebo in the FREEDOM study. This efficacy was supported by improvements from baseline in vertebral (18.2%) strength in axial compression and femoral (8.6%) strength in sideways fall configuration at 36 months, estimated in Newtons by an established voxel-based finite element (FE) methodology. Since FE analyses rely on the choice of meshes, material properties, and boundary conditions, the aim of this study was to independently confirm and compare the effects of denosumab on vertebral and femoral strength during the FREEDOM trial using an alternative smooth FE methodology. Unlike the previous FE study, effects on femoral strength in physiological stance configuration were also examined. QCT data for the proximal femur and two lumbar vertebrae were analyzed by smooth FE methodology at baseline, 12, 24, and 36 months for 51 treated (denosumab) and 47 control (placebo) subjects. QCT images were segmented and converted into smooth FE models to compute bone strength. L1 and L2 vertebral bodies were virtually loaded in axial compression and the proximal femora in both fall and stance configurations. Denosumab increased vertebral body strength by 10.8%, 14.0%, and 17.4% from baseline at 12, 24, and 36 months, respectively (p<0.0001). Denosumab also increased femoral strength in the fall configuration by 4.3%, 5.1%, and 7.2% from baseline at 12, 24, and 36 months, respectively (p<0.0001). Similar improvements were observed in the stance configuration with increases of 4.2%, 5.2%, and 5.2% from baseline (p≤0.0007). Differences between the increasing strengths with denosumab and the decreasing strengths with placebo were significant starting at 12 months (vertebral and femoral fall) or 24 months (femoral stance). Using an alternative smooth FE methodology, we confirmed the significant improvements in vertebral body and proximal femur strength previously observed with denosumab. Estimated increases in strength with denosumab and decreases with placebo were highly consistent between both FE techniques.

Automated three-dimensional registration of high-resolution peripheral quantitative computed tomography data to quantify size and shape changes of arthritic bone erosions.

OBJECTIVE: To monitor size and shape changes of bone erosions and changes in BMD in the vicinity of the erosion and in the periarticular trabecular compartment of patients with RA using high-resolution peripheral quantitative CT (HR-pQCT) imaging and to compare an automated three-dimensional (3D) image processing technique with manual measurements of erosion width and depth. METHODS: The shape of 40 bone erosions and composition of bone around the erosions were analysed in the MCP joints of 22 RA patients both manually and by semi-automated 3D image processing at two different time points. Periosteal segmentation was performed using volume growing and morphological operations. Image registration was applied for transfer of baseline segmentations to follow-up datasets. RESULTS: Eight erosions decreased in size, 6 increased and 28 remained stable. Increasing erosions were more spherical and smaller at baseline compared with decreasing or stable erosions. BMD in the vicinity of shrinking erosions increased, while it decreased next to expanding erosions. There was moderate agreement in the determination of erosion volume between semi-automated and manual measurements, but agreement was poor when assessing changes in volume over time. CONCLUSION: Longitudinal changes in erosion size and shape and of BMD in the vicinity of an erosion can be measured. BMD changes are associated with progression and regression of erosions. However, the semi-automated and manual approaches did not classify longitudinal changes of erosion volume in the same way. Further research is necessary to define the nature of these differences.

Registration of 2D histological sections with 3D micro-CT datasets from small animal vertebrae and tibiae.

The aim of this study was the registration of digitized thin 2D sections of mouse vertebrae and tibiae used for histomorphometry of trabecular bone structure into 3D micro computed tomography (µCT) datasets of the samples from which the sections were prepared. Intensity-based and segmentation-based registrations (SegRegs) of 2D sections and 3D µCT datasets were applied. As the 2D sections were deformed during their preparation, affine registration for the vertebrae was used instead of rigid registration. Tibiae sections were additionally cut on the distal end, which subsequently undergone more deformation so that elastic registration was necessary. The Jaccard distance was used as registration quality measure. The quality of intensity-based registrations and SegRegs was practically equal, although precision errors of the elastic registration of segmentation masks in tibiae were lower, while those in vertebrae were lower for the intensity-based registration. Results of SegReg significantly depended on the segmentation of the µCT datasets. Accuracy errors were reduced from approximately 64% to 42% when applying affine instead of rigid transformations for the vertebrae and from about 43% to 24% when using B-spline instead of rigid transformations for the tibiae. Accuracy errors can also be caused by the difference in spatial resolution between the thin sections (pixel size: 7.25 µm) and the µCT data (voxel size: 15 µm). In the vertebrae, average deformations amounted to a 6.7% shortening along the direction of sectioning and a 4% extension along the perpendicular direction corresponding to 0.13-0.17 mm. Maximum offsets in the mouse tibiae were 0.16 mm on average.

Characterization of knee osteoarthritis-related changes in trabecular bone using texture parameters at various levels of spatial resolution-a simulation study.

Articular cartilage and subchondral bone are the key tissues in osteoarthritis (OA). The role of the cancellous bone increasingly attracts attention in OA research. Because of its fast adaptation to changes in the loading distribution across joints, its quantification is expected to improve the diagnosis and monitoring of OA. In this study, we simulated OA progression-related changes of trabecular structure in a series of digital bone models and then characterized the potential of texture parameters and bone mineral density (BMD) as surrogate measures to quantify trabecular bone structure. Five texture parameters were studied: entropy, global and local inhomogeneity, anisotropy and variogram slope. Their dependence on OA relevant structural changes was investigated for three spatial resolutions typically used in micro computed tomography (CT; 10 µm), high-resolution peripheral quantitative CT (HR-pQCT) (90 µm) and clinical whole-body CT equipment (250 µm). At all resolutions, OA-related changes in trabecular bone architecture can be quantified using a specific (resolution dependent) combination of three texture parameters. BMD alone is inadequate for this purpose but if available reduces the required texture parameter combination to anisotropy and global inhomogeneity. The results are summarized in a comprehensive analysis guide for the detection of structural changes in OA knees. In conclusion, texture parameters can be used to characterize trabecular bone architecture even at spatial resolutions below the dimensions of a single trabecula and are essential for a detailed classification of relevant OA changes that cannot be achieved with a measurement of BMD alone.

Characterization and quantification of angiogenesis in rheumatoid arthritis in a mouse model using µCT.

BACKGROUND: Angiogenesis is an important pathophysiological process of chronic inflammation, especially in inflammatory arthritis. Quantitative measurement of changes in vascularization may improve the diagnosis and monitoring of arthritis. The aim of this work is the development of a 3D imaging and analysis framework for quantification of vascularization in experimental arthritis. METHODS: High-resolution micro-computed tomography (µCT) was used to scan knee joints of arthritic human tumor necrosis factor transgenic (hTNFtg) mice and non-arthritic wild-type controls previously perfused with lead-containing contrast agent Microfil MV-122. Vessel segmentation was performed by combination of intensity-based (local adaptive thresholding) and form-based (multi-scale method) segmentation techniques. Four anatomically defined concentric spherical shells centered in the knee joint were used as analysis volumes of interest. Vessel density, density distribution as well as vessel thickness, surface, spacing and number were measured. Simulated digital vessel tree models were used for validation of the algorithms. RESULTS: High-resolution µCT allows the quantitative assessment of the vascular tree in the knee joint during arthritis. Segmentation and analysis were highly automated but occasionally required manual corrections of the vessel segmentation close to the bone surfaces. Vascularization was significantly increased in arthritic hTNFtg mice compared to wild type controls. Precision errors for the morphologic parameters were smaller than 3% and 6% for intra- and interoperator analysis, respectively. Accuracy errors for vessel thickness were around 20% for vessels larger than twice the resolution of the scanner. CONCLUSIONS: Arthritis-induced changes of the vascular tree, including detailed and quantitative description of the number of vessel branches, length of vessel segments and the bifurcation angle, can be detected by contrast-enhanced high-resolution µCT.

Finite element analyses of human vertebral bodies embedded in polymethylmethalcrylate or loaded via the hyperelastic intervertebral disc models provide equivalent predictions of experimental strength.

Quantitative computer tomography (QCT)-based finite element (FE) models of vertebral body provide better prediction of vertebral strength than dual energy X-ray absorptiometry. However, most models were validated against compression of vertebral bodies with endplates embedded in polymethylmethalcrylate (PMMA). Yet, loading being as important as bone density, the absence of intervertebral disc (IVD) affects the strength. Accordingly, the aim was to assess the strength predictions of the classic FE models (vertebral body embedded) against the in vitro and in silico strengths of vertebral bodies loaded via IVDs. High resolution peripheral QCT (HR-pQCT) were performed on 13 segments (T11/T12/L1). T11 and L1 were augmented with PMMA and the samples were tested under a 4° wedge compression until failure of T12. Specimen-specific model was generated for each T12 from the HR-pQCT data. Two FE sets were created: FE-PMMA refers to the classical vertebral body embedded model under axial compression; FE-IVD to their loading via hyperelastic IVD model under the wedge compression as conducted experimentally. Results showed that FE-PMMA models overestimated the experimental strength and their strength prediction was satisfactory considering the different experimental set-up. On the other hand, the FE-IVD models did not prove significantly better (Exp/FE-PMMA: R²=0.68; Exp/FE-IVD: R²=0.71, p=0.84). In conclusion, FE-PMMA correlates well with in vitro strength of human vertebral bodies loaded via real IVDs and FE-IVD with hyperelastic IVDs do not significantly improve this correlation. Therefore, it seems not worth adding the IVDs to vertebral body models until fully validated patient-specific IVD models become available.

Segmentation and quantification of bone erosions in high-resolution peripheral quantitative computed tomography datasets of the metacarpophalangeal joints of patients with rheumatoid arthritis.

OBJECTIVE: To develop a precise three-dimensional (3D) segmentation technique for bone erosions in high-resolution peripheral quantitative CT (HR-pQCT) datasets to measure their volume, surface area and shape parameters. Assessment of bone erosions in patients with RA is important for diagnosis and evaluation of treatment efficacy. HR-pQCT allows quantifying periarticular bone loss in arthritis. METHODS: HR-pQCT scans with a spatial resolution of about 120 µm of the second to fourth metacarpophalangeal joints were acquired in patients with RA. Erosions were identified by placing a seed point in each of them. After applying 3D segmentation, the volume, surface area and sphericity of erosions were calculated. Results were compared with an approximation method using manual measurements. Intra- and interoperator precision analysis was performed for both the 3D segmentation and the manual technique. RESULTS: Forty-three erosions were assessed in 18 datasets. Intra- and interoperator precisions (RMSCV/RMSSD) for erosion volume were 5.66%/0.49 mm(3) and 7.76%/0.76 mm(3), respectively. The correlation between manual measurements and their simulation using segmentation volumes was r = 0.87. Precision errors for the manual method were 15.39% and 0.36 mm(3), respectively. CONCLUSION: We developed a new precise 3D segmentation technique for quantification of bone erosions in HR-pQCT datasets that correlates to the volume, shape and surface area of the erosion. The technique allows fast and effective measurement of the erosion size and could therefore be helpful for rapid and quantitative assessment of erosion size.

Binary segmentation masks can improve intrasubject registration accuracy of bone structures in CT images.

Registration of bone structures is a common problem in medical research as well as in clinical applications. Intrasubject rigid 3D monomodality registration of segmented bone structures of CT images and multimodality registration of muMR and segmented muCT bone images were performed with the multiresolution intensity-based technique implemented in ITK. The registration results for binary volumes of interest (VOI) masks and for segmented gray value VOIs were compared. To determine the registration quality, in the monomodality case the sum of squared difference, the sum of absolute differences, and the normalized symmetric difference of binary masks and in the multimodality case Mattes mutual information were applied. The use of binary VOI masks was significantly superior to the use of gray value VOIs.

Comparison of anatomic coordinate systems with rigid multi-resolution 3D registration for the reproducible positioning of analysis volumes of interest in QCT.

In this study we compared two approaches that have recently been used to minimize precision errors in 3D quantitative computed tomography (QCT) images of the hip and the spine in order to optimize the detection of longitudinal changes in bone mineral density (BMD). In 30 subjects we obtained baseline and 1 year follow-up 3D CT scans of the proximal femur and the spine. QCT analysis was applied to a variety of volumes of interest (VOIs) automatically positioned relative to anatomic coordinate systems (ACS). In the first approach (A1) baseline and follow-up scans were analyzed independently. In the second approach (A2) a 3D versor-based rigid intensity registration method was applied to match baseline and follow-up images, and the baseline ACS was mapped on the follow-up image using the registration transformation. Afterwards, the analysis VOIs were again independently calculated for baseline and follow-up images. There were no significant differences of percent BMD changes between baseline and follow-up images between A1 and A2 for any of the VOIs investigated. With advanced image processing methods a time-consuming 3D registration between baseline and follow-up images before the analysis does not improve analysis precision compared to the use of anatomical coordinate systems.