Search

The 3D distribution of the cortical and trabecular bone mass in the proximal femur is a critical component in determining fracture resistance that is not taken into account in clinical routine Dual-energy X-ray Absorptiometry (DXA) examination. In this paper, a statistical shape and appearance model together with a 3D-2D registration approach are used to model the femoral shape and bone density distribution in 3D from an anteroposterior DXA projection. A model-based algorithm is subsequently used to segment the cortex and build a 3D map of the cortical thickness and density. Measurements characterising the geometry and density distribution were computed for various regions of interest in both cortical and trabecular compartments. Models and measurements provided by the "3D-DXA" software algorithm were evaluated using a database of 157 study subjects, by comparing 3D-DXA analyses (using DXA scanners from three manufacturers) with measurements performed by Quantitative Computed Tomography (QCT). The mean point-to-surface distance between 3D-DXA and QCT femoral shapes was 0.93 mm. The mean absolute error between cortical thickness and density estimates measured by 3D-DXA and QCT was 0.33 mm and 72mg/cm(3). Correlation coefficients (R) between the 3D-DXA and QCT measurements were 0.86, 0.93, and 0.95 for the volumetric bone mineral density at the trabecular, cortical, and integral compartments respectively, and 0.91 for the mean cortical thickness. 3D-DXA provides a detailed analysis of the proximal femur, including a separate assessment of the cortical layer and trabecular macrostructure, which could potentially improve osteoporosis management while maintaining DXA as the standard routine modality.

Purpose: Cortical thickness and density are critical components in determining the strength of bony structures. Computed tomography(CT) is one possible modality for analyzing the cortex in 3D. In this paper, a model-based approach for measuring the cortical bone thickness and density from clinical CTimages is proposed. Methods: Density variations across the cortex were modeled as a function of the cortical thickness and density, location of the cortex, density of surrounding tissues, and imaging blur. High resolution micro-CT data of cadaver proximal femurs were analyzed to determine a relationship between cortical thickness and density. This thickness-density relationship was used as prior information to be incorporated in the model to obtain accurate measurements of cortical thickness and density from clinical CT volumes. The method was validated using micro-CT scans of 23 cadaver proximal femurs. Simulated clinical CTimages with different voxel sizes were generated from the micro-CT data. Cortical thickness and density were estimated from the simulated images using the proposed method and compared with measurements obtained using the micro-CT images to evaluate the effect of voxel size on the accuracy of the method. Then, 19 of the 23 specimens were imaged using a clinical CT scanner. Cortical thickness and density were estimated from the clinical CTimages using the proposed method and compared with the micro-CT measurements. Finally, a case-control study including 20 patients with osteoporosis and 20 age-matched controls with normal bone density was performed to evaluate the proposed method in a clinical context. Results: Cortical thickness (density) estimation errors were 0.07 ± 0.19 mm (−18 ± 92 mg/cm3) using the simulated clinical CT volumes with the smallest voxel size (0.33 × 0.33 × 0.5 mm3), and 0.10 ± 0.24 mm (−10 ± 115 mg/cm3) using the volumes with the largest voxel size (1.0 × 1.0 × 3.0 mm3). A trend for the cortical thickness and density estimation errors to increase with voxel size was observed and was more pronounced for thin cortices. Using clinical CT data for 19 of the 23 samples, mean errors of 0.18 ± 0.24 mm for the cortical thickness and 15 ± 106 mg/cm3 for the density were found. The case-control study showed that osteoporotic patients had a thinner cortex and a lower cortical density, with average differences of −0.8 mm and −58.6 mg/cm3 at the proximal femur in comparison with age-matched controls (p-value < 0.001). Conclusions: This method might be a promising approach for the quantification of cortical bone thickness and density using clinical routine imaging techniques. Future work will concentrate on investigating how this approach can improve the estimation of mechanical strength of bony structures, the prevention of fracture, and the management of osteoporosis.

Dual Energy X-ray Absorptiometry (DXA) is the standard exam for osteoporosis diagnosis and fracture risk evaluation at the spine. However, numerous patients with bone fragility are not diagnosed as such. In fact, standard analysis of DXA images does not differentiate between trabecular and cortical bone; neither specifically assess of the bone density in the vertebral body, which is where most of the osteoporotic fractures occur. Quantitative computed tomography (QCT) is an alternative technique that overcomes limitations of DXA-based diagnosis. However, due to the high cost and radiation dose, QCT is not used for osteoporosis management. We propose a method that provides a 3-D subject-specific shape and density estimation of the lumbar spine from a single anteroposterior (AP) DXA image. A 3-D statistical shape and density model is built, using a training set of QCT scans, and registered onto the AP DXA image so that its projection matches it. Cortical and trabecular bone compartments are segmented using a model-based algorithm. Clinical measurements are performed at different bone compartments. Accuracy was evaluated by comparing DXA-derived to QCT-derived 3-D measurements for a validation set of 180 subjects. The shape accuracy was 1.51 mm at the total vertebra and 0.66 mm at the vertebral body. Correlation coefficients between DXA and QCT-derived measurements ranged from 0.81 to 0.97. The method proposed offers an insightful 3-D analysis of the lumbar spine, which could potentially improve osteoporosis and fracture risk assessment in patients who had an AP DXA scan of the lumbar spine without any additional examination.

Fragility fractures are an important public health issue. The aim of this study was to analyze the association of the main osteoporotic risk factors related to fragility fracture in a cohort of women with an indication of bone densitometry (BD).A retrospective cohort was followed-up until a fragile fracture occurred, in a population of women aged 40 to 90 years with a first visit for BD between January 1992 and February 2008. We calculated the incidence rate of fracture per 1000 women-years of follow-up, and the hazard ratio (HR) of fragile fracture using a Cox regression model.A total of 49,735 women were studied. The average age of participants was 57.8 years (SD: 8.5). Of these, 3631 women (7.1%) reported a new fragility fracture in post-baseline visits. Risk factors with higher adjusted HR were age ≥ 75 years compared with age55 years (HR: 3.8; 95% CI: 3.3-4.4) and having a BC result evaluated as osteoporosis compared to normal (HR: 2.0; 95% CI: 1.8-2.2). A personal history of humerus, hip or vertebral fractures had an adjusted HR of 1.2 (95% CI: 1.1-1.3).The main risk factors for fragility fracture were advanced age, BD result and a personal history of fracture, although 74% of fractures were detected with a bone mineral density classified as normal or osteopenia. Other relevant factors were rheumatoid arthritis or having received prolonged corticosteroid therapy.

This work presents a statistical model of both the shape and Bone Mineral Density (BMD) distribution of the proximal femur for fracture risk assessment. The shape and density model was built from a dataset of Quantitative Computed Tomography scans of fracture patients and a control group. Principal Component Analysis and Horn's parallel analysis were used to reduce the dimensionality of the shape and density model to the main modes of variation. The input data was then used to analyze the model parameters for the optimal separation between the fracture and control group. Feature selection using the Fisher criterion determined the parameters with the best class separation, which were used in Fisher Linear Discriminant Analysis to find the direction in the parameter space that best separates the fracture and control group. This resulted in a Fisher criterion value of 6.70, while analyzing the Dualenergy X-ray Absorptiometry derived femur neck areal BMD of the same subjects resulted in a Fisher criterion value of 0.98. This indicates that a fracture risk estimation approach based on the presented model might improve upon the current standard clinical practice.

The diagnosis of osteoporosis and the prevention of femur fractures is a major challenge for our society. However, the diagnosis performed in clinical routine from Dual Energy X-ray Absorptiometry (DXA) images is limited. This paper proposes a 3D reconstruction method of both the shape and the Bone Mineral Density (BMD) distribution of the proximal femur from routinely used DXA images. The reconstruction accuracy that can be obtained from single-view and multi-view DXA devices was assessed. This evaluation, from 20 bone specimens and simulated DXA images, highlighted a mean shape accuracy of 1.3mm and a BMD accuracy of 4.4% from a single-view DXA image. A multi-view configuration with 2 views (frontal-sagittal) appeared as a good compromise (mean shape accuracy of 0.9mm and BMD accuracy of 3.2%). We are currently using this method for in vivo clinical studies in order to improve the diagnosis of osteoporosis and the prevention of femur fractures.

Comunicació presentada al congrés: Medical Imaging 2010: Image Processing (SPIE) celebrat del 13 al 18 de febrer de 2010 a San Diego, California. Area Bone Mineral Density (aBMD) measured by Dual-energy X-ray Absorptiometry (DXA) is an established criterion in the evaluation of hip fracture risk. The evaluation from these planar images, however, is limited to 2D while it has been shown that proper 3D assessment of both the shape and the Bone Mineral Density (BMD) distribution improves the fracture risk estimation. In this work we present a method to reconstruct both the 3D bone shape and 3D BMD distribution of the proximal femur from a single DXA image. A statistical model of shape and a separate statistical model of the BMD distribution were automatically constructed from a set of Quantitative Computed Tomography (QCT) scans. The reconstruction method incorporates a fully automatic intensity based 3D-2D registration process, maximizing the similarity between the DXA and a digitally reconstructed radiograph of the combined model. For the construction of the models, an in vitro dataset of QCT scans of 60 anatomical specimens was used. To evaluate the reconstruction accuracy, experiments were performed on simulated DXA images from the QCT scans of 30 anatomical specimens. Comparisons between the reconstructions and the same subject QCT scans showed a mean shape accuracy of 1.2mm, and a mean density error of 81mg/cm3. The results show that this method is capable of accurately reconstructing both the 3D shape and 3D BMD distribution of the proximal femur from DXA images used in clinical routine, potentially improving the diagnosis of osteoporosis and fracture risk assessments at a low radiation dose and low cost. This research has been financially supported by the 3D-FemOs project (3-Dimensional Femur Reconstruction for Osteoporotic Fracture Risk Assessment) with a grant awarded by ACC1Ó Valtec (Tipus 1). In addition this work has been supported by a grant from the Instituto de Salud Carlos III (FI09/00757) and by a grant from the \Deutsche Forschungsgemeinschaft" (LO 730 / 3-1). The authors acknowledge Benedikt Schuler (Institute for Biomedical Image Analysis, UMIT, Hall in Tirol, Austria) and Dr. E. M. Lochmüuller (Universitäts-Frauenklinik der LMU, Mnchen, Germany) for the data collection.