A novel methodology for generating 3D finite element models of the hip from 2D radiographs

2D-3D reconstruction of the proximal femur from DXA scans: Evaluation of the 3D-Shaper software

Introduction: Osteoporosis is currently diagnosed based on areal bone mineral density (aBMD) computed from 2D DXA scans. However, aBMD is a limited surrogate for femoral strength since it does not account for 3D bone geometry and density distribution. QCT scans combined with finite element (FE) analysis can deliver improved femoral strength predictions. However, non-negligible radiation dose and high costs prevent a systematic usage of this technique for screening purposes. As an alternative, the 3D-Shaper software (3D-Shaper Medical, Spain) reconstructs the 3D shape and density distribution of the femur from 2D DXA scans. This approach could deliver a more accurate estimation of femoral strength than aBMD by using FE analysis on the reconstructed 3D DXA. Methods: Here we present the first independent evaluation of the software, using a dataset of 77 ex vivo femora. We extend a prior evaluation by including the density distribution differences, the spatial correlation of density values and an FE analysis. Yet, cortical thickness is left out of this evaluation, since the cortex is not resolved in our FE models. Results: We found an average surface distance of 1.16 mm between 3D DXA and QCT images, which shows a good reconstruction of the bone geometry. Although BMD values obtained from 3D DXA and QCT correlated well (r ² = 0.92), the 3D DXA BMD were systematically lower. The average BMD difference amounted to 64 mg/cm³, more than one-third of the 3D DXA BMD. Furthermore, the low correlation (r ² = 0.48) between density values of both images indicates a limited reconstruction of the 3D density distribution. FE results were in good agreement between QCT and 3D DXA images, with a high coefficient of determination (r ² = 0.88). However, this correlation was not statistically different from a direct prediction by aBMD. Moreover, we found differences in the fracture patterns between the two image types. QCT-based FE analysis resulted mostly in femoral neck fractures and 3D DXA-based FE in subcapital or pertrochanteric fractures. Discussion: In conclusion, 3D-Shaper generates an altered BMD distribution compared to QCT but, after careful density calibration, shows an interesting potential for deriving a standardized femoral strength from a DXA scan.

Femoral strength can be predicted from 2D projections using a 3D statistical deformation and texture model with finite element analysis

Ultimate force of the proximal human femur can be predicted using Finite Element Analysis (FEA), but the models rely on 3D computed tomography images. Landmark-based statistical appearance models (SAM) and B-Spline transformation-based statistical deformation models (SDM) have been used to estimate 3D images from 2D projections, which facilitates model generation and reduces the radiation dose. However, there is no literature on the accuracy of SDM-based FEA models of bones with respect to experimental results. In this study, a methodology for an enhanced SDM with textural information is presented. The statistical deformation and texture models (SDTMs) are based on a set of 37 quantitative CT (QCT) images. They were used to estimate 3D images from two or one projections of the set in a leave-one-out setup. These estimations where then used to create FEA models. The ultimate force predicted by FEA models estimated from two or one projection using the SDTMs were compared to the experimental ultimate force from a previous study on the same femora and to the results of standard QCT-based FEA models. High correlations between predictions and experimental measurements were found for FEA models reconstructed from 2D projections with R²=0.835 when based on two projections and R²=0.724 when using one projection. The correlations were comparable to those reached with standard QCT-based FE-models with the experimental results (R²=0.795). This study shows the high potential of SDTM-based 3D image reconstruction and FEA modelling from 2D projections to predict femoral ultimate force.

Validation of 3D finite element models from simulated DXA images for biofidelic simulations of sideways fall impact to the hip

Computed tomography (CT)-derived finite element (FE) models have been proposed as a tool to improve the current clinical assessment of osteoporosis and personalized hip fracture risk by providing an accurate estimate of femoral strength. However, this solution has two main drawbacks, namely: (i) 3D CT images are needed, whereas 2D dual-energy x-ray absorptiometry (DXA) images are more generally available, and (ii) quasi-static femoral strength is predicted as a surrogate for fracture risk, instead of predicting whether a fall would result in a fracture or not. The aim of this study was to combine a biofidelic fall simulation technique, based on 3D computed tomography (CT) data with an algorithm that reconstructs 3D femoral shape and BMD distribution from a 2D DXA image. This approach was evaluated on 11 pelvis-femur constructs for which CT scans, ex vivo sideways fall impact experiments and CT-derived biofidelic FE models were available. Simulated DXA images were used to reconstruct the 3D shape and bone mineral density (BMD) distribution of the left femurs by registering a projection of a statistical shape and appearance model with a genetic optimization algorithm. The 2D-to-3D reconstructed femurs were meshed, and the resulting FE models inserted into a biofidelic FE modeling pipeline for simulating a sideways fall. The median 2D-to-3D reconstruction error was 1.02 mm for the shape and 0.06 g/cm³ for BMD for the 11 specimens. FE models derived from simulated DXAs predicted the outcome of the falls in terms of fracture versus non-fracture with the same accuracy as the CT-derived FE models. This study represents a milestone towards improved assessment of hip fracture risk based on widely available clinical DXA images.

3D Analysis of the Proximal Femur Compared to 2D Analysis for Hip Fracture Risk Prediction in a Clinical Population

Due to the adverse impacts of hip fractures on patients' lives, it is crucial to enhance the identification of people at high risk through accessible clinical techniques. Reconstructing the 3D geometry and BMD distribution of the proximal femur could be beneficial in enhancing hip fracture risk predictions; however, it is associated with a high computational burden. It is also not clear whether it provides a better performance than 2D model analysis. Therefore, the purpose of this study was to compare the 2D and 3D model reconstruction's ability to predict hip fracture risk in a clinical population of patients. The DXA scans and CT scans of 16 cadaveric femurs were used to create training sets for the 2D and 3D model reconstruction based on statistical shape and appearance modeling. Subsequently, these methods were used to predict the risk of sustaining a hip fracture in a clinical population of 150 subjects (50 fractured, and 100 non-fractured) that were monitored for five years in the Canadian Multicentre Osteoporosis Study. 3D model reconstruction was able to improve the identification of patients who sustained a hip fracture more accurately than the standard clinical practice (by 40%). Also, the predictions from the 2D statistical model didn't differ significantly from the 3D ones (p > 0.76). These results indicated that to enhance hip fracture risk prediction in clinical practice implementing 2D statistical modeling has comparable performance with lower associated computational load.

Perspectives on the non-invasive evaluation of femoral strength in the assessment of hip fracture risk

We reviewed the experimental and clinical evidence that hip bone strength estimated by BMD and/or finite element analysis (FEA) reflects the actual strength of the proximal femur and is associated with hip fracture risk and its changes upon treatment.

INTRODUCTION

The risk of hip fractures increases exponentially with age due to a progressive loss of bone mass, deterioration of bone structure, and increased incidence of falls. Areal bone mineral density (aBMD), measured by dual-energy X-ray absorptiometry (DXA), is the most used surrogate marker of bone strength. However, age-related declines in bone strength exceed those of aBMD, and the majority of fractures occur in those who are not identified as osteoporotic by BMD testing. With hip fracture incidence increasing worldwide, the development of accurate methods to estimate bone strength in vivo would be very useful to predict the risk of hip fracture and to monitor the effects of osteoporosis therapies.

METHODS

We reviewed experimental and clinical evidence regarding the association between aBMD and/orCT-finite element analysis (FEA) estimated femoral strength and hip fracture risk as well as their changes with treatment.

RESULTS

Femoral aBMD and bone strength estimates by CT-FEA explain a large proportion of femoral strength ex vivo and predict hip fracture risk in vivo. Changes in femoral aBMD are strongly associated with anti-fracture efficacy of osteoporosis treatments, though comparable data for FEA are currently not available.

CONCLUSIONS

Hip aBMD and estimated femoral strength are good predictors of fracture risk and could potentially be used as surrogate endpoints for fracture in clinical trials. Further improvements of FEA may be achieved by incorporating trabecular orientations, enhanced cortical modeling, effects of aging on bone tissue ductility, and multiple sideway fall loading conditions.

Present and future for technologies to develop patient-specific medical devices: a systematic review approach

The main purpose of this investigation was to systematically review the literature regarding case studies on patient-specific implants and devices, with the goal of analyzing the process of developing custom-made medical devices. A content analysis was performed to identify design processes and methodologies implemented to develop devices such as implants adapted to bone geometries. Reverse engineering, computer-aided design, simulation of assets, and rapid prototyping technologies were selected according to their interoperability in a process framework for developing new products. Finally, results from the case studies and process stages identified in the consulted research were analyzed. These results showed a relationship between the scope and complexity of the process and the stage of technology integration of the patient-specific device development. The analyzed case studies were characterized by technical, scientific, and multidisciplinary components to achieve research goals. Likewise, integration of technologies using patient-specific technologies is needed for product development that converges into designing devices, such as implants, biomodels, and cutting drilling guides.

An analytical method to create patient-specific deformed bone models using X-ray images and a healthy bone model

Generation of patient-specific bone models from X-ray images is useful for various medical applications such as total hip replacement, implant manufacturing, knee kinematic studies and deformity correction. These models may provide valuable information required for a more reliable operation. In this work, we propose a new algorithm for generating patient-specific 3D models of femur and tibia with deformity, using only a generic healthy bone model and some simple measurements taken on the X-ray images of the diseased bone. Using the X-ray measurements, an interpolation function (a polynomial or a cubic spline) is fit to the mid-diaphyseal curve of the actual bone and the generic bone model is deformed in the guidance of this function with free form deformation method. The created models are intended to be used mainly for the visualization of fixation procedure in software-supported external fixation systems. An error measure is defined to quantify the error in this matching procedure. The method is found to be capable of producing deformed tibia models that satisfactorily reflect the actual bones, as confirmed by two orthopaedic surgeons who use software-supported external fixation systems regularly.

Interlaboratory comparison of femur surface reconstruction from CT data compared to reference optical 3D scan

BACKGROUND

The present study contrasts the accuracy of different reconstructed models with distinctive segmentation methods performed by various experts. Seven research groups reconstructed nine 3D models of one human femur based on an acquired CT image using their own computational methods. As a reference model for accuracy assessment, a 3D surface scan of the human femur was created using an optical measuring system. Prior to comparison, the femur was divided into four areas; "neck and greater trochanter", "proximal metaphysis", "diaphysis", and "distal metaphysis". The deviation analysis was carried out in GEOMAGIC studio v.2013 software.

RESULTS

The results revealed that the highest deviation errors occurred in "neck and greater trochanter" area and "proximal metaphysis" area with RMSE of 0.84 and 0.83 mm respectively.

CONCLUSION

In conclusion, this study shows that the average deviation of reconstructed models prepared by experts with various methods, skills and software from the surface 3D scan is lower than 0.79 mm, which is not a significant discrepancy.

An Imaged-Based Three-Dimensional Study of First Metatarsal Protrusion Distance in Women with and Without Hallux Valgus

BACKGROUND

First metatarsal protrusion distance (MPD) has been commonly studied as a characteristic of hallux valgus deformity. To date, the majority of investigations have used radiographic methods, with most reporting first metatarsal (ray) protrusion to be associated with deformity. As an alternative, this study used a three-dimensional (3-D) image acquisition and data analysis method to quantify MPD.

METHODS

Magnetic resonance images were acquired in weightbearing on 29 women (19 with hallux valgus; 10 controls). After the 3-D images were reconstructed into virtual bone models, two examiners measured MPD in relation to the navicular. In addition to a reliability analysis, a t test assessed for group differences in demographics, foot posture (hallux valgus, intermetatarsal angles), and MPD.

RESULTS

Group demographics were not different, while measures of hallux valgus and intermetatarsal angles were different ( P < 0.01) between groups. The measurement of MPD was highly reliable (ICC [Formula: see text] 0.99; SEM [Formula: see text] 0.78 mm). Metatarsal protrusion averaged approximately -2.0 mm in both groups. There was no statistical group difference ( P = 0.89) in MPD.

CONCLUSIONS

The reconstructed image datasets captured the 3-D spatial relationship of the anatomy. Measurements of MPD were reliable. The first ray measured 2 mm shorter than the second ray in both the hallux valgus and control groups. Though unexpected, this result may prompt future study of the pathokinematics associated with hallux valgus that include the quantification of metatarsal protrusion with 3-D methods, instead of relying solely on single-plane radiograph reports.

Hip joint geometry effects on cartilage contact stresses during a gait cycle

The cartilage surface geometry of natural human hip joint is commonly regarded as sphere. It has been widely applied in computational simulation and hip joint prosthesis design. Some new geometry models have been developed and the sphere assumption has been questioned recently. The objective of this study was to analyze joint geometry effects on cartilage contact stress distribution and investigate contact patterns during a whole gait cycle. Hip surface was reconstructed from CT data of a healthy volunteer. Three finite element (FE) models of hip joint were developed from different cartilage geometries: natural geometry, sphere and rotational ellipsoid. Loads at ten instants of gait cycle were applied to these models based on published in-vivo data. FE predictions of peak contact pressure during gait of natural hip were compared with sphere and rotational ellipsoid replaced hip joint. Contact occurs mainly in upper anterior region of both acetabulum and femur distributing along sagittal plane of human body. It moves towards inferolateral aspect as the resultant joint reaction force changes during walking for natural hip. Peak pressures at the instant with maximum contact force were 7.48 MPa, 14.97 MPa and 13.12 MPa for models with natural hip surface, sphere replaced and rotational ellipsoid replaced surface respectively. During the whole gait cycle, contact pressure of natural hip ranked lowest in most of the instants, followed by rotational ellipsoid replaced and sphere replaced hip. The results indicate that rotational ellipsoid is more consistent with natural hip cartilage geometry than sphere during normal walking. This means rotational ellipsoid prosthesis could give a better description of physiological structure compared with standard sphere prosthesis. Therefore, rotational ellipsoid would be a better choice for prosthesis design.

Prediction of femoral strength using 3D finite element models reconstructed from DXA images: validation against experiments

Computed tomography (CT)-based finite element (FE) models may improve the current osteoporosis diagnostics and prediction of fracture risk by providing an estimate for femoral strength. However, the need for a CT scan, as opposed to the conventional use of dual-energy X-ray absorptiometry (DXA) for osteoporosis diagnostics, is considered a major obstacle. The 3D shape and bone mineral density (BMD) distribution of a femur can be reconstructed using a statistical shape and appearance model (SSAM) and the DXA image of the femur. Then, the reconstructed shape and BMD could be used to build FE models to predict bone strength. Since high accuracy is needed in all steps of the analysis, this study aimed at evaluating the ability of a 3D FE model built from one 2D DXA image to predict the strains and fracture load of human femora. Three cadaver femora were retrieved, for which experimental measurements from ex vivo mechanical tests were available. FE models were built using the SSAM-based reconstructions: using only the SSAM-reconstructed shape, only the SSAM-reconstructed BMD distribution, and the full SSAM-based reconstruction (including both shape and BMD distribution). When compared with experimental data, the SSAM-based models predicted accurately principal strains (coefficient of determination >0.83, normalized root-mean-square error <16%) and femoral strength (standard error of the estimate 1215 N). These results were only slightly inferior to those obtained with CT-based FE models, but with the considerable advantage of the models being built from DXA images. In summary, the results support the feasibility of SSAM-based models as a practical tool to introduce FE-based bone strength estimation in the current fracture risk diagnostics.

Finite Element-Based Mechanical Assessment of Bone Quality on the Basis of In Vivo Images

Beyond bone mineral density (BMD), bone quality designates the mechanical integrity of bone tissue. In vivo images based on X-ray attenuation, such as CT reconstructions, provide size, shape, and local BMD distribution and may be exploited as input for finite element analysis (FEA) to assess bone fragility. Further key input parameters of FEA are the material properties of bone tissue. This review discusses the main determinants of bone mechanical properties and emphasizes the added value, as well as the important assumptions underlying finite element analysis. Bone tissue is a sophisticated, multiscale composite material that undergoes remodeling but exhibits a rather narrow band of tissue mineralization. Mechanically, bone tissue behaves elastically under physiologic loads and yields by cracking beyond critical strain levels. Through adequate cell-orchestrated modeling, trabecular bone tunes its mechanical properties by volume fraction and fabric. With proper calibration, these mechanical properties may be incorporated in quantitative CT-based finite element analysis that has been validated extensively with ex vivo experiments and has been applied increasingly in clinical trials to assess treatment efficacy against osteoporosis.

In-vivo assessment of femoral bone strength using Finite Element Analysis (FEA) based on routine MDCT imaging: a preliminary study on patients with vertebral fractures

Purpose

To experimentally validate a non-linear finite element analysis (FEA) modeling approach assessing in-vitro fracture risk at the proximal femur and to transfer the method to standard in-vivo multi-detector computed tomography (MDCT) data of the hip aiming to predict additional hip fracture risk in subjects with and without osteoporosis associated vertebral fractures using bone mineral density (BMD) measurements as gold standard.

Methods

One fresh-frozen human femur specimen was mechanically tested and fractured simulating stance and clinically relevant fall loading configurations to the hip. After experimental in-vitro validation, the FEA simulation protocol was transferred to standard contrast-enhanced in-vivo MDCT images to calculate individual hip fracture risk each for 4 subjects with and without a history of osteoporotic vertebral fractures matched by age and gender. In addition, FEA based risk factor calculations were compared to manual femoral BMD measurements of all subjects.

Results

In-vitro simulations showed good correlation with the experimentally measured strains both in stance (R² = 0.963) and fall configuration (R² = 0.976). The simulated maximum stress overestimated the experimental failure load (4743 N) by 14.7% (5440 N) while the simulated maximum strain overestimated by 4.7% (4968 N). The simulated failed elements coincided precisely with the experimentally determined fracture locations. BMD measurements in subjects with a history of osteoporotic vertebral fractures did not differ significantly from subjects without fragility fractures (femoral head: p = 0.989; femoral neck: p = 0.366), but showed higher FEA based risk factors for additional incident hip fractures (p = 0.028).

Conclusion

FEA simulations were successfully validated by elastic and destructive in-vitro experiments. In the subsequent in-vivo analyses, MDCT based FEA based risk factor differences for additional hip fractures were not mirrored by according BMD measurements. Our data suggests, that MDCT derived FEA models may assess bone strength more accurately than BMD measurements alone, providing a valuable in-vivo fracture risk assessment tool.

High-spatial-resolution bone densitometry with dual-energy X-ray absorptiometric region-free analysis

PURPOSE

To outline the conceptual development of dual-energy absorptiometric (DXA) region-free analysis, quantify its precision, and evaluate its application to quantify the change in longitudinal femoral periprosthetic bone mineral density (BMD) in patients during the 12 months after total hip arthroplasty.

MATERIALS AND METHODS

All subjects had undergone total hip arthroplasty for idiopathic arthritis, and the scans were collected as part of previous ethically approved studies (1998-2005) for which informed consent was provided. Contemporary image processing approaches were used to develop a region of interest-free DXA analysis method with increased spatial resolution for assessment of proximal femoral BMD. The method was calibrated, and its accuracy relative to a proprietary algorithm was assessed by using a hip phantom. The precision of the method was examined by using repeat DXA acquisitions in 29 patients, and its ability to allow spatial resolution of localized periprosthetic BMD change at the hip was assessed in an independent group of 19 patients who were measured throughout a 12-month period. Differences were evaluated with t tests (P < .05).

RESULTS

The method allowed spatial resolution of more than 10 000 individual BMD data points on a typical archived prosthetic hip scan. The median data point-level error of the method after calibration was -1.9% (interquartile range, -7.2% to 3.5%) relative to a proprietary algorithm. The median data point-level precision, expressed as a coefficient of variation, was 1.4% (interquartile range, 1.2%-1.6%). Evaluation of BMD change in a model of periprosthetic bone loss demonstrated large but highly focal changes in BMD that would not be resolved by using traditional region of interest-based analysis approaches.

CONCLUSION

The proposed approach provides a quantitative, precise method for extracting high-spatial-resolution BMD data from existing DXA datasets without the limitations imposed by region of interest-based analysis.