1
|
Dudle A, Gugler Y, Pretterklieber M, Ferrari S, Lippuner K, Zysset P. 2D-3D reconstruction of the proximal femur from DXA scans: Evaluation of the 3D-Shaper software. Front Bioeng Biotechnol 2023; 11:1111020. [PMID: 36937766 PMCID: PMC10014626 DOI: 10.3389/fbioe.2023.1111020] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/29/2022] [Accepted: 02/15/2023] [Indexed: 03/05/2023] Open
Abstract
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 2 = 0.92), the 3D DXA BMD were systematically lower. The average BMD difference amounted to 64 mg/cm3, more than one-third of the 3D DXA BMD. Furthermore, the low correlation (r 2 = 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 2 = 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.
Collapse
Affiliation(s)
- Alice Dudle
- ARTORG Center for Biomedical Engineering Research, University of Bern, Bern, Switzerland
- *Correspondence: Alice Dudle, ; Yvan Gugler,
| | - Yvan Gugler
- ARTORG Center for Biomedical Engineering Research, University of Bern, Bern, Switzerland
- *Correspondence: Alice Dudle, ; Yvan Gugler,
| | - Michael Pretterklieber
- Division of Anatomy, Gottfried Schatz Research Center, Medical University of Graz, Graz, Austria
- Division of Anatomy, Center for Anatomy and Cell Biology, Medical University of Vienna, Vienna, Austria
| | - Serge Ferrari
- Division of Bone Diseases, Geneva University Hospitals (HUG), Geneva, Switzerland
| | - Kurt Lippuner
- Department of Osteoporosis, Inselspital, Bern University Hospital, University of Bern, Bern, Switzerland
| | - Philippe Zysset
- ARTORG Center for Biomedical Engineering Research, University of Bern, Bern, Switzerland
| |
Collapse
|
2
|
Femoral strength can be predicted from 2D projections using a 3D statistical deformation and texture model with finite element analysis. Med Eng Phys 2021; 93:72-82. [PMID: 34154777 DOI: 10.1016/j.medengphy.2021.05.012] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/18/2020] [Revised: 04/11/2021] [Accepted: 05/18/2021] [Indexed: 11/22/2022]
Abstract
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 R2=0.835 when based on two projections and R2=0.724 when using one projection. The correlations were comparable to those reached with standard QCT-based FE-models with the experimental results (R2=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.
Collapse
|
3
|
Grassi L, Fleps I, Sahlstedt H, Väänänen SP, Ferguson SJ, Isaksson H, Helgason B. Validation of 3D finite element models from simulated DXA images for biofidelic simulations of sideways fall impact to the hip. Bone 2021; 142:115678. [PMID: 33022451 DOI: 10.1016/j.bone.2020.115678] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 07/03/2020] [Revised: 09/11/2020] [Accepted: 09/30/2020] [Indexed: 10/23/2022]
Abstract
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/cm3 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.
Collapse
Affiliation(s)
- Lorenzo Grassi
- Department of Biomedical Engineering, Lund University, Lund, Sweden.
| | - Ingmar Fleps
- Institute for Biomechanics, ETH Zürich, Zürich, Switzerland
| | | | - Sami P Väänänen
- Department of Applied Physics, University of Eastern Finland, Kuopio, Finland; Diagnostic Imaging Center, Kuopio University Hospital, Kuopio, Finland
| | | | - Hanna Isaksson
- Department of Biomedical Engineering, Lund University, Lund, Sweden
| | | |
Collapse
|
4
|
Jazinizadeh F, Quenneville CE. 3D Analysis of the Proximal Femur Compared to 2D Analysis for Hip Fracture Risk Prediction in a Clinical Population. Ann Biomed Eng 2020; 49:1222-1232. [PMID: 33123827 DOI: 10.1007/s10439-020-02670-2] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/20/2020] [Accepted: 10/20/2020] [Indexed: 01/22/2023]
Abstract
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.
Collapse
Affiliation(s)
- Fatemeh Jazinizadeh
- Department of Mechanical Engineering, McMaster University, ABB-C308, 1280 Main St. West, Hamilton, ON, L8S 4L8, Canada
| | - Cheryl E Quenneville
- Department of Mechanical Engineering, McMaster University, ABB-C308, 1280 Main St. West, Hamilton, ON, L8S 4L8, Canada.
- School of Biomedical Engineering, McMaster University, Hamilton, ON, Canada.
| |
Collapse
|
5
|
Bouxsein ML, Zysset P, Glüer CC, McClung M, Biver E, Pierroz DD, Ferrari SL. Perspectives on the non-invasive evaluation of femoral strength in the assessment of hip fracture risk. Osteoporos Int 2020; 31:393-408. [PMID: 31900541 DOI: 10.1007/s00198-019-05195-0] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 08/23/2019] [Accepted: 10/04/2019] [Indexed: 10/25/2022]
Abstract
UNLABELLED 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.
Collapse
Affiliation(s)
- M L Bouxsein
- Center for Advanced Orthopedic Studies, Beth Israel Deaconess Medical Center, and Department of Orthopedic Surgery, Harvard Medical School, Boston, MA, USA
| | - P Zysset
- ARTORG Center for Biomedical Engineering Research, University of Bern, Bern, Switzerland
| | - C C Glüer
- Section of Biomedical Imaging, Department of Radiology and Neuroradiology, University Medical Center of Schleswig-Holstein, Campus Kiel, Kiel, Germany
| | - M McClung
- Oregon Osteoporosis Center, Portland, OR, USA
- Mary MacKillop Institute for Health Research, Australian Catholic University, Melbourne, VIC, Australia
| | - E Biver
- Division of Bone Disease, Department of Internal Medicine Specialties, Faculty of Medicine, Geneva University Hospital, Geneva, Switzerland
| | - D D Pierroz
- International Osteoporosis Foundation (IOF), Nyon, Switzerland
| | - S L Ferrari
- Division of Bone Disease, Department of Internal Medicine Specialties, Faculty of Medicine, Geneva University Hospital, Geneva, Switzerland.
| | | |
Collapse
|
6
|
López Gualdrón CI, Bravo Ibarra ER, Murillo Bohórquez AP, Garnica Bohórquez I. Present and future for technologies to develop patient-specific medical devices: a systematic review approach. MEDICAL DEVICES-EVIDENCE AND RESEARCH 2019; 12:253-273. [PMID: 31496840 PMCID: PMC6689557 DOI: 10.2147/mder.s215947] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/16/2019] [Accepted: 07/08/2019] [Indexed: 01/25/2023] Open
Abstract
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.
Collapse
Affiliation(s)
| | - Edna-Rocío Bravo Ibarra
- Industrial and Business Studies School, Universidad Industrial de Santander, Bucaramanga, Colombia
| | | | - Israel Garnica Bohórquez
- Industrial and Business Studies School, Universidad Industrial de Santander, Bucaramanga, Colombia
| |
Collapse
|
7
|
Ün MK, Avşar E, Akçalı İD. An analytical method to create patient-specific deformed bone models using X-ray images and a healthy bone model. Comput Biol Med 2018; 104:43-51. [PMID: 30423529 DOI: 10.1016/j.compbiomed.2018.11.003] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/11/2018] [Revised: 10/20/2018] [Accepted: 11/03/2018] [Indexed: 10/27/2022]
Abstract
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.
Collapse
Affiliation(s)
- M Kerem Ün
- Department of Electrical and Electronics Engineering, Çukurova University, 01330, Sarıçam, Adana, Turkey; Department of Biomedical Engineering, Çukurova University, 01330, Sarıçam, Adana, Turkey.
| | - Ercan Avşar
- Department of Electrical and Electronics Engineering, Çukurova University, 01330, Sarıçam, Adana, Turkey.
| | - İbrahim D Akçalı
- MACTIMARUM Research and Application Center, Çukurova University, 01330, Sarıçam, Adana, Turkey.
| |
Collapse
|
8
|
Soodmand E, Kluess D, Varady PA, Cichon R, Schwarze M, Gehweiler D, Niemeyer F, Pahr D, Woiczinski M. Interlaboratory comparison of femur surface reconstruction from CT data compared to reference optical 3D scan. Biomed Eng Online 2018; 17:29. [PMID: 29495963 PMCID: PMC5833145 DOI: 10.1186/s12938-018-0461-0] [Citation(s) in RCA: 23] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/17/2017] [Accepted: 02/23/2018] [Indexed: 11/10/2022] Open
Abstract
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.
Collapse
Affiliation(s)
- Ehsan Soodmand
- Biomechanics and Implant Technology Research Laboratory, Department of Orthopaedics, University Medicine Rostock, Doberaner Strasse 142, 18057, Rostock, Germany.
| | - Daniel Kluess
- Biomechanics and Implant Technology Research Laboratory, Department of Orthopaedics, University Medicine Rostock, Doberaner Strasse 142, 18057, Rostock, Germany
| | - Patrick A Varady
- Trauma Center Murnau Institute of Biomechanics, Professor-Küntscher-Str. 882418, Murnau am Staffelsee, Germany
| | - Robert Cichon
- Chair of Mechanics and Robotics, University Duisburg-Essen, Lotharstrasse 1, 47057, Duisburg, Germany
| | - Michael Schwarze
- Laboratory for Biomechanics and Biomaterials of the Orthopaedic Clinic, Hannover Medical School, Anna-von-Borries-Strasse 1-7, 30625, Hannover, Germany
| | - Dominic Gehweiler
- Department of Trauma, Hand and Reconstructive Surgery, University Hospital Münster, Albert-Schweitzer-Campus 1, 48149, Münster, Germany
| | - Frank Niemeyer
- Fraunhofer Research Institution for Large Structures in Production Engineering (IGP), Albert-Einstein-Str. 30, 18059, Rostock, Germany
| | - Dieter Pahr
- Institute of Lightweight Design and Structural Biomechanics, TU Vienna, Getreidemarkt 9, 1060, Vienna, Austria
| | - Matthias Woiczinski
- Department of Orthopedic Surgery, Physical Medicine and Rehabilitation, University Hospital of Munich (LMU), Marchioninistr. 15, 81377, Munich, Germany
| |
Collapse
|
9
|
Rodriguez MW, Pioske KE, Glasoe WM. An Imaged-Based Three-Dimensional Study of First Metatarsal Protrusion Distance in Women with and Without Hallux Valgus. J Am Podiatr Med Assoc 2017; 107:531-537. [PMID: 28768118 DOI: 10.7547/16-104] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Indexed: 02/03/2023]
Abstract
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.
Collapse
Affiliation(s)
| | - Kara E. Pioske
- Division in Physical Therapy, University of Minnesota, Minneapolis, Minnesota
| | - Ward M. Glasoe
- Division in Physical Therapy, University of Minnesota, Minneapolis, Minnesota
| |
Collapse
|
10
|
Taylor M, Perilli E, Martelli S. Development of a surrogate model based on patient weight, bone mass and geometry to predict femoral neck strains and fracture loads. J Biomech 2017; 55:121-127. [DOI: 10.1016/j.jbiomech.2017.02.022] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/09/2016] [Revised: 11/18/2016] [Accepted: 02/16/2017] [Indexed: 10/20/2022]
|
11
|
Hip joint geometry effects on cartilage contact stresses during a gait cycle. ANNUAL INTERNATIONAL CONFERENCE OF THE IEEE ENGINEERING IN MEDICINE AND BIOLOGY SOCIETY. IEEE ENGINEERING IN MEDICINE AND BIOLOGY SOCIETY. ANNUAL INTERNATIONAL CONFERENCE 2017; 2016:6038-6041. [PMID: 28269629 DOI: 10.1109/embc.2016.7592105] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
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.
Collapse
|
12
|
Grassi L, Väänänen SP, Ristinmaa M, Jurvelin JS, Isaksson H. Prediction of femoral strength using 3D finite element models reconstructed from DXA images: validation against experiments. Biomech Model Mechanobiol 2016; 16:989-1000. [PMID: 28004226 PMCID: PMC5422489 DOI: 10.1007/s10237-016-0866-2] [Citation(s) in RCA: 25] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/30/2016] [Accepted: 12/10/2016] [Indexed: 12/01/2022]
Abstract
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.
Collapse
Affiliation(s)
- Lorenzo Grassi
- Department of Biomedical Engineering, Lund University, BMC D13, 221 84, Lund, Sweden.
| | - Sami P Väänänen
- Department of Applied Physics, University of Eastern Finland, Kuopio, Finland
- Department of Orthopaedics, Traumatology and Hand Surgery, Kuopio University Hospital, Kuopio, Finland
| | | | - Jukka S Jurvelin
- Department of Applied Physics, University of Eastern Finland, Kuopio, Finland
- Diagnostic Imaging Center, Kuopio University Hospital, Kuopio, Finland
| | - Hanna Isaksson
- Department of Biomedical Engineering, Lund University, BMC D13, 221 84, Lund, Sweden
| |
Collapse
|
13
|
Abstract
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.
Collapse
Affiliation(s)
- Dieter H Pahr
- Institute of Lightweight Design and Structural Biomechanics, Vienna University of Technology, Vienna, Austria
| | - Philippe K Zysset
- Institute for Surgical Technology and Biomechanics, University of Bern, Bern, Switzerland.
| |
Collapse
|
14
|
Generation of 3D shape, density, cortical thickness and finite element mesh of proximal femur from a DXA image. Med Image Anal 2015; 24:125-134. [DOI: 10.1016/j.media.2015.06.001] [Citation(s) in RCA: 59] [Impact Index Per Article: 6.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/03/2014] [Revised: 06/03/2015] [Accepted: 06/11/2015] [Indexed: 11/19/2022]
|
15
|
Liebl H, Garcia EG, Holzner F, Noel PB, Burgkart R, Rummeny EJ, Baum T, Bauer JS. 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. PLoS One 2015; 10:e0116907. [PMID: 25723187 PMCID: PMC4344329 DOI: 10.1371/journal.pone.0116907] [Citation(s) in RCA: 26] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/21/2014] [Accepted: 12/16/2014] [Indexed: 01/26/2023] Open
Abstract
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 (R2 = 0.963) and fall configuration (R2 = 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.
Collapse
Affiliation(s)
- Hans Liebl
- Department of Radiology, Klinikum rechts der Isar, Technische Universitaet Muenchen, Muenchen, Germany
| | - Eduardo Grande Garcia
- Department of Radiology, Klinikum rechts der Isar, Technische Universitaet Muenchen, Muenchen, Germany; Department of Orthopaedic Surgery, Klinikum rechts der Isar, Technische Universitaet Muenchen, Muenchen, Germany
| | - Fabian Holzner
- Department of Orthopaedic Surgery, Klinikum rechts der Isar, Technische Universitaet Muenchen, Muenchen, Germany
| | - Peter B Noel
- Department of Radiology, Klinikum rechts der Isar, Technische Universitaet Muenchen, Muenchen, Germany
| | - Rainer Burgkart
- Department of Orthopaedic Surgery, Klinikum rechts der Isar, Technische Universitaet Muenchen, Muenchen, Germany
| | - Ernst J Rummeny
- Department of Radiology, Klinikum rechts der Isar, Technische Universitaet Muenchen, Muenchen, Germany
| | - Thomas Baum
- Department of Radiology, Klinikum rechts der Isar, Technische Universitaet Muenchen, Muenchen, Germany
| | - Jan S Bauer
- Section of Neuroradiology, Klinikum rechts der Isar, Technische Universitaet Muenchen, Muenchen, Germany
| |
Collapse
|
16
|
Morris RM, Yang L, Martín-Fernández MA, Pozo JM, Frangi AF, Wilkinson JM. High-spatial-resolution bone densitometry with dual-energy X-ray absorptiometric region-free analysis. Radiology 2014; 274:532-9. [PMID: 25222069 DOI: 10.1148/radiol.14140636] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Abstract
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.
Collapse
Affiliation(s)
- Richard M Morris
- From the Academic Unit of Bone Metabolism, University of Sheffield, Northern General Hospital, Beech Hill Road, Sheffield S10 2RX, England (R.M.M., L.Y., J.M.W.); and Centre for Computational Imaging & Simulation Technologies in Biomedicine (CISTIB), Department of Mechanical Engineering, University of Sheffield, Sheffield, England (M.A.M.F., J.M.P., A.F.F.)
| | | | | | | | | | | |
Collapse
|