Search

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.

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 CT images is proposed.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 CT images 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 CT images 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.Cortical thickness (density) estimation errors were 0.07 ± 0.19 mm (-18 ± 92 mg/cm(3)) using the simulated clinical CT volumes with the smallest voxel size (0.33 × 0.33 × 0.5 mm(3)), and 0.10 ± 0.24 mm (-10 ± 115 mg/cm(3)) using the volumes with the largest voxel size (1.0 × 1.0 × 3.0 mm(3)). 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/cm(3) 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/cm(3) at the proximal femur in comparison with age-matched controls (p-value0.001).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.

Purpose: Dual-energy x-ray absorptiometry (DXA) is used in clinical routine to provide a two-dimensional (2D) analysis of the bone mineral density (BMD). 3D reconstruction methods from 2D DXA images could improve the BMD analysis. To find the optimal configuration that should be used in clinical routine, this paper relies on a 3D reconstruction method from DXA images to compare the accuracy that can be obtained from one single-view and from multiview DXA images (two to four projections). Methods: The 3D reconstruction method uses a statistical model and a nonrigid registration technique to recover in 3D the shape and the BMD distribution of the proximal femur. The accuracy was evaluatedin vivo by comparing 3D reconstructions obtained from simulated DXA images of 30 patients (using between one and four DXA views) with quantitative computed tomographyreconstructions. Results: This comparison showed that the use of one single DXA provides accurate 3D reconstructions (mean shape accuracy of 1.0 mm and BMD distribution errors of 7.0%). Among the multiview configurations, the use of two views (0° and 45°) was the best compromise, increasing the accuracy of pose (mean accuracy of 0.7°/1.2°/0.9° against 1.0°/3.5°/3.3° for the single view), reducing slightly the BMD errors (5.7%) while maintaining the same shape accuracy. Conclusions: The use of two views constitutes an interesting configuration when multiview DXA devices are available in clinical routine. However, the use of only one single view remains an accurate solution to recover the shape and the BMD distribution in 3D, with the advantage of a higher potential for clinical translation.

Dual-energy x-ray absorptiometry (DXA) is used in clinical routine to provide a two-dimensional (2D) analysis of the bone mineral density (BMD). 3D reconstruction methods from 2D DXA images could improve the BMD analysis. To find the optimal configuration that should be used in clinical routine, this paper relies on a 3D reconstruction method from DXA images to compare the accuracy that can be obtained from one single-view and from multiview DXA images (two to four projections).The 3D reconstruction method uses a statistical model and a nonrigid registration technique to recover in 3D the shape and the BMD distribution of the proximal femur. The accuracy was evaluated in vivo by comparing 3D reconstructions obtained from simulated DXA images of 30 patients (using between one and four DXA views) with quantitative computed tomography reconstructions.This comparison showed that the use of one single DXA provides accurate 3D reconstructions (mean shape accuracy of 1.0 mm and BMD distribution errors of 7.0%). Among the multiview configurations, the use of two views (0° and 45°) was the best compromise, increasing the accuracy of pose (mean accuracy of 0.7°/1.2°/0.9° against 1.0°/3.5°/3.3° for the single view), reducing slightly the BMD errors (5.7%) while maintaining the same shape accuracy.The use of two views constitutes an interesting configuration when multiview DXA devices are available in clinical routine. However, the use of only one single view remains an accurate solution to recover the shape and the BMD distribution in 3D, with the advantage of a higher potential for clinical translation.