1
|
Szyszko JA, Aldieri A, La Mattina AA, Viceconti M. Phantomless calibration of CT scans for hip fracture risk prediction in silico: Comparison with phantom-based calibration. PLoS One 2024; 19:e0305474. [PMID: 38875268 PMCID: PMC11178222 DOI: 10.1371/journal.pone.0305474] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/11/2024] [Accepted: 05/30/2024] [Indexed: 06/16/2024] Open
Abstract
Finite element models built from quantitative computed tomography images rely on element-wise mapping of material properties starting from Hounsfield Units (HU), which can be converted into mineral densities upon calibration. While calibration is preferably carried out by scanning a phantom with known-density components, conducting phantom-based calibration may not always be possible. In such cases, a phantomless procedure, where the scanned subject's tissues are used as a phantom, is an interesting alternative. The aim of this study was to compare a phantom-based and a phantomless calibration method on 41 postmenopausal women. The proposed phantomless calibration utilized air, adipose, and muscle tissues, with reference equivalent mineral density values of -797, -95, and 38 mg/cm3, extracted from a previously performed phantom-based calibration. A 9-slice volume of interest (VOI) centred between the femoral head and knee rotation centres was chosen. Reference HU values for air, adipose, and muscle tissues were extracted by identifying HU distribution peaks within the VOI, and patient-specific calibration was performed using linear regression. Comparison of FE models calibrated with the two methods showed average relative differences of 1.99% for Young's modulus1.30% for tensile and 1.34% for compressive principal strains. Excellent correlations (R2 > 0.99) were identified for superficial maximum tensile and minimum compressive strains. Maximum normalised root mean square relative error (RMSRE) values settled at 4.02% for Young's modulus, 2.99% for tensile, and 3.22% for compressive principal strains, respectively. The good agreement found between the two methods supports the adoption of the proposed methodology when phantomless calibration is needed.
Collapse
Affiliation(s)
- Julia A Szyszko
- Medical Technology Lab, IRCCS Istituto Ortopedico Rizzoli, Bologna, Italy
- Department of Industrial Engineering, Alma Mater Studiorum-University of Bologna, Bologna, Italy
| | - Alessandra Aldieri
- Medical Technology Lab, IRCCS Istituto Ortopedico Rizzoli, Bologna, Italy
- PolitoBIOMed Lab, Department of Mechanical and Aerospace Engineering, Politecnico di Torino, Torino, Italy
| | - Antonino A La Mattina
- Medical Technology Lab, IRCCS Istituto Ortopedico Rizzoli, Bologna, Italy
- Department of Industrial Engineering, Alma Mater Studiorum-University of Bologna, Bologna, Italy
| | - Marco Viceconti
- Medical Technology Lab, IRCCS Istituto Ortopedico Rizzoli, Bologna, Italy
- Department of Industrial Engineering, Alma Mater Studiorum-University of Bologna, Bologna, Italy
| |
Collapse
|
2
|
Maquer G, Mueri C, Henderson A, Bischoff J, Favre P. Developing and Validating a Model of Humeral Stem Primary Stability, Intended for In Silico Clinical Trials. Ann Biomed Eng 2024; 52:1280-1296. [PMID: 38361138 DOI: 10.1007/s10439-024-03452-w] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/31/2023] [Accepted: 01/12/2024] [Indexed: 02/17/2024]
Abstract
In silico clinical trials (ISCT) can contribute to demonstrating a device's performance via credible computational models applied on virtual cohorts. Our purpose was to establish the credibility of a model for assessing the risk of humeral stem loosening in total shoulder arthroplasty, based on a twofold validation scheme involving both benchtop and clinical validation activities, for ISCT applications. A finite element model computing bone-implant micromotion (benchtop model) was quantitatively compared to a bone foam micromotion test (benchtop comparator) to ensure that the physics of the system was captured correctly. The model was expanded to a population-based approach (clinical model) and qualitatively evaluated based on its ability to replicate findings from a published clinical study (clinical comparator), namely that grit-blasted stems are at a significantly higher risk of loosening than porous-coated stems, to ensure that clinical performance of the stem can be predicted appropriately. Model form sensitivities pertaining to surgical variation and implant design were evaluated. The model replicated benchtop micromotion measurements (52.1 ± 4.3 µm), without a significant impact of the press-fit ("Press-fit": 54.0 ± 8.5 µm, "No press-fit": 56.0 ± 12.0 µm). Applied to a virtual population, the grit-blasted stems (227 ± 78µm) experienced significantly larger micromotions than porous-coated stems (162 ± 69µm), in accordance with the findings of the clinical comparator. This work provides a concrete example for evaluating the credibility of an ISCT study. By validating the modeling approach against both benchtop and clinical data, model credibility is established for an ISCT application aiming to enrich clinical data in a regulatory submission.
Collapse
Affiliation(s)
- Ghislain Maquer
- Zimmer Biomet, Sulzerallee 8, 8404, Winterthur, Switzerland.
| | | | - Adam Henderson
- Zimmer Biomet, Sulzerallee 8, 8404, Winterthur, Switzerland
| | - Jeff Bischoff
- Zimmer Biomet, 1800 West Center St., Warsaw, IN, 46580, USA
| | | |
Collapse
|
3
|
Bartenschlager S, Cavallaro A, Pogarell T, Chaudry O, Uder M, Khosla S, Schett G, Engelke K. Opportunistic Screening With CT: Comparison of Phantomless BMD Calibration Methods. J Bone Miner Res 2023; 38:1689-1699. [PMID: 37732678 DOI: 10.1002/jbmr.4917] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 04/17/2023] [Revised: 09/08/2023] [Accepted: 09/17/2023] [Indexed: 09/22/2023]
Abstract
Opportunistic screening is a new promising technique to identify individuals at high risk for osteoporotic fracture using computed tomography (CT) scans originally acquired for an clinical purpose unrelated to osteoporosis. In these CT scans, a calibration phantom traditionally required to convert measured CT values to bone mineral density (BMD) is missing. As an alternative, phantomless calibration has been developed. This study aimed to review the principles of four existing phantomless calibration methods and to compare their performance against the gold standard of simultaneous calibration (ΔBMD). All methods were applied to a dataset of 350 females scanned with a highly standardized CT protocol (DS1) and to a second dataset of 114 patients (38 female) from clinical routine covering a large range of CT acquisition and reconstruction parameters (DS2). Three of the phantomless calibration methods must be precalibrated with a reference dataset containing a calibration phantom. Sixty scans from DS1 and 57 from DS2 were randomly selected for this precalibration. For each phantomless calibration method first the best combination of internal reference materials (IMs) was selected. These were either air and blood or subcutaneous adipose tissue, blood, and cortical bone. In addition, for phantomless calibration a fifth method based on average calibration parameters derived from the reference dataset was applied. For DS1, ΔBMD results (mean ± standard deviation) for the phantomless calibration methods requiring a precalibration ranged from 0.1 ± 2.7 mg/cm3 to 2.4 ± 3.5 mg/cm3 with similar means but significantly higher standard deviations for DS2. Performance of the phantomless calibration method, which does not require a precalibration was worse (ΔBMD DS1: 12.6 ± 13.2 mg/cm3 , DS2: 0.5 ± 8.8 mg/cm3 ). In conclusion, phantomless BMD calibration performs well if precalibrated with a reference dataset. © 2023 The Authors. Journal of Bone and Mineral Research published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research (ASBMR).
Collapse
Affiliation(s)
- Stefan Bartenschlager
- Department of Medicine 3, Friedrich-Alexander-University Erlangen-Nürnberg and Universitätsklinikum Erlangen, Erlangen, Germany
- Institute of Medical Physics, Friedrich-Alexander-University Erlangen-Nürnberg, Erlangen, Germany
| | - Alexander Cavallaro
- Institute of Radiology, Friedrich-Alexander-University Erlangen-Nürnberg and Universitätsklinikum Erlangen, Erlangen, Germany
| | - Tobias Pogarell
- Institute of Radiology, Friedrich-Alexander-University Erlangen-Nürnberg and Universitätsklinikum Erlangen, Erlangen, Germany
| | - Oliver Chaudry
- Department of Medicine 3, Friedrich-Alexander-University Erlangen-Nürnberg and Universitätsklinikum Erlangen, Erlangen, Germany
- Institute of Medical Physics, Friedrich-Alexander-University Erlangen-Nürnberg, Erlangen, Germany
| | - Michael Uder
- Institute of Radiology, Friedrich-Alexander-University Erlangen-Nürnberg and Universitätsklinikum Erlangen, Erlangen, Germany
| | - Sundeep Khosla
- Division of Endocrinology and Kogod Center on Aging, Mayo Clinic, Rochester, MN, USA
| | - Georg Schett
- Department of Medicine 3, Friedrich-Alexander-University Erlangen-Nürnberg and Universitätsklinikum Erlangen, Erlangen, Germany
| | - Klaus Engelke
- Department of Medicine 3, Friedrich-Alexander-University Erlangen-Nürnberg and Universitätsklinikum Erlangen, Erlangen, Germany
- Institute of Medical Physics, Friedrich-Alexander-University Erlangen-Nürnberg, Erlangen, Germany
| |
Collapse
|
4
|
Eliyahu L, Yosibash Z, Avivi I, Cohen YC, Ariel G, Sadovnic O, Sternheim A. On the influence of computed tomography's slice thickness on computer tomography based finite element analyses results. Clin Biomech (Bristol, Avon) 2023; 102:105889. [PMID: 36774735 DOI: 10.1016/j.clinbiomech.2023.105889] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 07/03/2022] [Revised: 12/06/2022] [Accepted: 01/09/2023] [Indexed: 01/19/2023]
Abstract
BACKGROUND Patient-specific autonomous finite element analyses of femurs, based on clinical computed tomography scans may be used to monitor the progression of bone-related diseases. Some CT scan protocols provide lower resolution (slice thickness of 3 mm) that affects the accuracy. To investigate the impact of low-resolution scans on the CT-based finite element analyses results, identical CT raw data were reconstructed twice to generate a 1 mm ("gold standard") and a 3 mm slice thickness scans. METHODS CT-based finite element analyses of twenty-four femurs (twelve patients) under stance and sideways fall loads were performed based on 1 and 3 mm slice thickness scans. Bone volume, load direction, and strains were extracted at different locations along the femurs and differences were evaluated. FINDINGS Average differences in bone volume were 1.0 ± 1.5%. The largest average difference in strains in stance position was in the neck region (11.0 ± 13.4%), whereas in other regions these were much smaller. For sidewise fall loading, the average differences were at most 9.2 ± 16.0%. INTERPRETATION Whole-body low dose CT scans (3 mm-slice thickness) are suboptimal for monitoring strain changes in patient's femurs but may allow longitudinal studies if larger than 5% in all areas and larger than 12% in the upper neck. CT-based finite element analyses with slice thickness of 3 mm may be used in clinical practice for patients with smoldering myeloma to associate changes in strains with progression to active myeloma if above ∼10%.
Collapse
Affiliation(s)
- Leetal Eliyahu
- Computational Mechanics and Experimental Biomechanics Lab, School of Mechanical Engineering, Tel-Aviv University, Israel
| | - Zohar Yosibash
- Computational Mechanics and Experimental Biomechanics Lab, School of Mechanical Engineering, Tel-Aviv University, Israel.
| | - Irit Avivi
- Sackler Faculty of Medicine, Tel Aviv University, Tel Aviv, Israel; Hematology Division, Tel Aviv Sourasky Medical Center, Tel Aviv, Israel
| | - Yael C Cohen
- Sackler Faculty of Medicine, Tel Aviv University, Tel Aviv, Israel; Hematology Division, Tel Aviv Sourasky Medical Center, Tel Aviv, Israel
| | - Gal Ariel
- Sackler Faculty of Medicine, Tel Aviv University, Tel Aviv, Israel; National Unit of Orthopaedic Oncology, Tel Aviv Sourasky Medical Center, Tel Aviv, Israel
| | - Ofer Sadovnic
- Sackler Faculty of Medicine, Tel Aviv University, Tel Aviv, Israel; Radiology Division, Tel Aviv Sourasky Medical Center, Tel Aviv, Israel
| | - Amir Sternheim
- Sackler Faculty of Medicine, Tel Aviv University, Tel Aviv, Israel; National Unit of Orthopaedic Oncology, Tel Aviv Sourasky Medical Center, Tel Aviv, Israel
| |
Collapse
|
5
|
Abstract
PURPOSE OF REVIEW Opportunistic screening is a combination of techniques to identify subjects of high risk for osteoporotic fracture using routine clinical CT scans prescribed for diagnoses unrelated to osteoporosis. The two main components are automated detection of vertebral fractures and measurement of bone mineral density (BMD) in CT scans, in which a phantom for calibration of CT to BMD values is not used. This review describes the particular challenges of opportunistic screening and provides an overview and comparison of current techniques used for opportunistic screening. The review further outlines the performance of opportunistic screening. RECENT FINDINGS A wide range of technologies for the automatic detection of vertebral fractures have been developed and successfully validated. Most of them are based on artificial intelligence algorithms. The automated differentiation of osteoporotic from traumatic fractures and vertebral deformities unrelated to osteoporosis, the grading of vertebral fracture severity, and the detection of mild vertebral fractures is still problematic. The accuracy of automated fracture detection compared to classical radiological semi-quantitative Genant scoring is about 80%. Accuracy errors of alternative BMD calibration methods compared to simultaneous phantom-based calibration used in standard quantitative CT (QCT) range from below 5% to about 10%. The impact of contrast agents, frequently administered in clinical CT on the determination of BMD and on fracture risk determination is still controversial. Opportunistic screening, the identification of vertebral fracture and the measurement of BMD using clinical routine CT scans, is feasible but corresponding techniques still need to be integrated into the clinical workflow and further validated with respect to the prediction of fracture risk.
Collapse
Affiliation(s)
- Klaus Engelke
- Department of Medicine III, Friedrich-Alexander University of Erlangen-Nürnberg, University Hospital Erlangen, Ulmenweg 18, 91054, Erlangen, Germany.
- Institute of Medical Physics (IMP), Friedrich-Alexander-Universität Erlangen-Nürnberg (FAU), Henkestr. 91, 91052, Erlangen, Germany.
| | - Oliver Chaudry
- Department of Medicine III, Friedrich-Alexander University of Erlangen-Nürnberg, University Hospital Erlangen, Ulmenweg 18, 91054, Erlangen, Germany
- Institute of Medical Physics (IMP), Friedrich-Alexander-Universität Erlangen-Nürnberg (FAU), Henkestr. 91, 91052, Erlangen, Germany
| | - Stefan Bartenschlager
- Department of Medicine III, Friedrich-Alexander University of Erlangen-Nürnberg, University Hospital Erlangen, Ulmenweg 18, 91054, Erlangen, Germany
- Institute of Medical Physics (IMP), Friedrich-Alexander-Universität Erlangen-Nürnberg (FAU), Henkestr. 91, 91052, Erlangen, Germany
| |
Collapse
|
6
|
Eggermont F, van der Linden Y, Verdonschot N, Dierselhuis E, Ligthert S, Bitter T, Westhoff P, Tanck E. A Patient-Specific Fracture Risk Assessment Tool for Femoral Bone Metastases: Using the Bone Strength (BOS) Score in Clinical Practice. Cancers (Basel) 2022; 14:cancers14235904. [PMID: 36497388 PMCID: PMC9740241 DOI: 10.3390/cancers14235904] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/13/2022] [Revised: 10/18/2022] [Accepted: 11/24/2022] [Indexed: 12/04/2022] Open
Abstract
Patients with femoral metastases are at risk of fracturing bones. It is important to prevent fractures in order to maintain mobility and quality of life. The BOne Strength (BOS) score is based on a computed tomography (CT)-based patient-specific finite element (FE) computer model that objectively calculates bone strength. In this pilot study, the added clinical value of the BOS score towards treatment-related decision making was assessed. In December 2019, the BOS score was implemented in four radiotherapy centers. The BOS scores and fracture risks of individual patients were calculated and returned to the physician to assist in treatment decisions. The physicians filled out a questionnaire, which was qualitatively analyzed. A follow-up to identify fractures and/or death was performed after six months. Until June 2021, 42 BOS scores were delivered (20 high, 9 moderate, and 13 low fracture risk). In 48%, the BOS score led to an adaptation of treatment plans. Physicians indicated that the BOS score provided objective insight into fracture risk, was reassuring for physicians and patients, and improved multidisciplinary discussions and shared decision making. In conclusion, the BOS score is an objective tool to assess fracture risk in femoral bone metastases and aids physicians and patients in making a more informed decision regarding the most appropriate treatment.
Collapse
Affiliation(s)
- Florieke Eggermont
- Orthopaedic Research Laboratory, Radboud Institute for Health Sciences, Radboud University Medical Center, 6525 GA Nijmegen, The Netherlands
- Correspondence:
| | - Yvette van der Linden
- Department of Radiotherapy, Leiden University Medical Center, 2333 ZA Leiden, The Netherlands
- Netherlands Comprehensive Cancer Organisation (IKNL), 3511 DT Utrecht, The Netherlands
| | - Nico Verdonschot
- Orthopaedic Research Laboratory, Radboud Institute for Health Sciences, Radboud University Medical Center, 6525 GA Nijmegen, The Netherlands
- Laboratory of Biomechanical Engineering, University of Twente, 7522 NB Enschede, The Netherlands
| | - Edwin Dierselhuis
- Department of Orthopedics, Radboud Institute for Health Sciences, Radboud University Medical Center, 6525 GA Nijmegen, The Netherlands
| | - Steven Ligthert
- Orthopaedic Research Laboratory, Radboud Institute for Health Sciences, Radboud University Medical Center, 6525 GA Nijmegen, The Netherlands
| | - Thom Bitter
- Orthopaedic Research Laboratory, Radboud Institute for Health Sciences, Radboud University Medical Center, 6525 GA Nijmegen, The Netherlands
| | - Paulien Westhoff
- Department of Radiation Oncology, Radboud Institute for Health Sciences, Radboud University Medical Center, 6525 GA Nijmegen, The Netherlands
| | - Esther Tanck
- Orthopaedic Research Laboratory, Radboud Institute for Health Sciences, Radboud University Medical Center, 6525 GA Nijmegen, The Netherlands
| |
Collapse
|