1
|
Koh CWY, Lew KS, Wibawa A, Master Z, Yeap PL, Chua CGA, Lee JCL, Tan HQ, Park SY. First clinical experience following the consensus guide for calibrating a proton stopping power ratio curve in a new proton centre. Phys Med 2024; 120:103341. [PMID: 38554639 DOI: 10.1016/j.ejmp.2024.103341] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 08/22/2023] [Revised: 03/07/2024] [Accepted: 03/25/2024] [Indexed: 04/02/2024] Open
Abstract
BACKGROUND AND PURPOSE This work introduces the first assessment of CT calibration following the ESTRO's consensus guidelines and validating the HLUT through the irradiation of biological material. METHODS Two electron density phantoms were scanned with two CT scanners using two CT scan energies. The stopping power ratio (SPR) and mass density (MD) HLUTs for different CT scan energies were derived using Schneider's and ESTRO's methods. The comparison metric in this work is based on the Water-Equivalent Thickness (WET) difference between the treatment planning system and biological irradiation measurement. The SPR HLUTs were compared between the two calibration methods. To assess the accuracy of using MD HLUT for dose calculation in the treatment planning system, MD vs SPR HLUT was compared. Lastly, the feasibility of using a single SPR HLUT to replace two different energy CT scans was explored. RESULTS The results show a WET difference of less than 3.5% except for the result in the Bone region between Schneider's and ESTRO's methods. Comparing MD and SPR HLUT, the results from MD HLUT show less than a 3.5% difference except for the Bone region. However, the SPR HLUT shows a lower mean absolute percentage difference as compared to MD HLUT between the measured and calculated WET difference. Lastly, it is possible to use a single SPR HLUT for two different CT scan energies since both WET differences are within 3.5%. CONCLUSION This is the first report on calibrating an HLUT following the ESTRO's guidelines. While our result shows incremental improvement in range uncertainty using the ESTRO's guideline, the prescriptional approach of the guideline does promote harmonization of CT calibration protocols between different centres.
Collapse
Affiliation(s)
| | - Kah Seng Lew
- Division of Radiation Oncology, National Cancer Centre Singapore, Singapore; Nanyang Technological University Singapore, Singapore
| | - Andrew Wibawa
- Division of Radiation Oncology, National Cancer Centre Singapore, Singapore
| | - Zubin Master
- Division of Radiation Oncology, National Cancer Centre Singapore, Singapore
| | - Ping Lin Yeap
- Division of Radiation Oncology, National Cancer Centre Singapore, Singapore; Department of Oncology, University of Cambridge, United Kingdom
| | | | - James Cheow Lei Lee
- Division of Radiation Oncology, National Cancer Centre Singapore, Singapore; Nanyang Technological University Singapore, Singapore
| | - Hong Qi Tan
- Division of Radiation Oncology, National Cancer Centre Singapore, Singapore; Oncology Academic Clinical Programme, Duke-NUS Medical School, Singapore.
| | - Sung Yong Park
- Division of Radiation Oncology, National Cancer Centre Singapore, Singapore; Oncology Academic Clinical Programme, Duke-NUS Medical School, Singapore
| |
Collapse
|
2
|
Pettersson E, Thilander-Klang A, Bäck A. Prediction of proton stopping power ratios using dual-energy CT basis material decomposition. Med Phys 2024; 51:881-897. [PMID: 38194501 DOI: 10.1002/mp.16929] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/21/2023] [Revised: 12/04/2023] [Accepted: 12/15/2023] [Indexed: 01/11/2024] Open
Abstract
BACKGROUND Proton radiotherapy treatment plans are currently restricted by the range uncertainties originating from the stopping power ratio (SPR) prediction based on single-energy computed tomography (SECT). Various studies have shown that multi-energy CT (MECT) can reduce the range uncertainties due to medical implant materials and age-related variations in tissue composition. None of these has directly applied the basis material density (MD) images produced by projection-based MECT systems for SPR prediction. PURPOSE To present and evaluate a novel proton SPR prediction method based on MD images from dual-energy CT (DECT), which could reduce the range uncertainties currently associated with proton radiotherapy. METHODS A theoretical basis material decomposition into water and iodine material densities was performed for various pediatric and adult human reference tissues, as well as other non-tissue materials, by minimizing the root-mean-square relative attenuation error in the energy interval from 40 to 140 keV. A model (here called MD-SPR) mapping predicted MDs to theoretically calculated reference SPRs was created with locally weighted scatterplot smoothing (LOWESS) data-fitting. The goodness of fit of the MD-SPR model was evaluated for the included reference tissues. MD images of two electron density phantoms, combined to form a head- and an abdomen-sized phantom setup, were acquired with a clinical projection-based fast-kV switching DECT scanner. The MD images were compared to the theoretically predicted MDs of the tissue surrogates and other non-tissue materials in the phantoms, as well as used for input to the MD-SPR model for generation of SPR images. The SPR images were subsequently compared to theoretical reference SPRs of the materials in the phantoms, as well as to SPR images from a commercial algorithm (DirectSPR, Siemens Healthineers, Forchheim, Germany) using image-based consecutive scan DECT for the same phantom setups. RESULTS The predicted SPRs of the tissue surrogates were similar for MD-SPR and DirectSPR, where the adipose and bone tissue surrogates were within 1% difference to the reference SPRs, while other non-adipose soft tissue surrogates (breast, brain, liver, muscle) were all underestimated by between -0.7% and -1.8%. The SPRs of the non-tissue materials (polymethyl methacrylate (PMMA), polyether ether ketone (PEEK), graphite and Teflon) were within 2.8% for MD-SPR images, compared to 6.8% for DirectSPR. CONCLUSIONS The MD-SPR model performed similar compared to other published methods for the human reference tissues. The SPR prediction for tissue surrogates was similar to DirectSPR and showed potential to improve SPR prediction for non-tissue materials.
Collapse
Affiliation(s)
- Erik Pettersson
- Department of Medical Radiation Sciences, Institute of Clinical Sciences, Sahlgrenska Academy, University of Gothenburg, Gothenburg, Sweden
- Therapeutic Radiation Physics, Department of Medical Physics and Biomedical Engineering, Sahlgrenska University Hospital, Gothenburg, Sweden
| | - Anne Thilander-Klang
- Department of Medical Radiation Sciences, Institute of Clinical Sciences, Sahlgrenska Academy, University of Gothenburg, Gothenburg, Sweden
- Diagnostic Radiation Physics, Department of Medical Physics and Biomedical Engineering, Sahlgrenska University Hospital, Gothenburg, Sweden
| | - Anna Bäck
- Department of Medical Radiation Sciences, Institute of Clinical Sciences, Sahlgrenska Academy, University of Gothenburg, Gothenburg, Sweden
- Therapeutic Radiation Physics, Department of Medical Physics and Biomedical Engineering, Sahlgrenska University Hospital, Gothenburg, Sweden
| |
Collapse
|
3
|
Taasti VT, Wohlfahrt P. From computed tomography innovation to routine clinical application in radiation oncology - A joint initiative of close collaboration. Phys Imaging Radiat Oncol 2024; 29:100550. [PMID: 38390587 PMCID: PMC10881422 DOI: 10.1016/j.phro.2024.100550] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/24/2024] Open
Affiliation(s)
- Vicki Trier Taasti
- Danish Centre for Particle Therapy, Aarhus University Hospital, Aarhus, Denmark
| | - Patrick Wohlfahrt
- Siemens Healthineers, Varian, Cancer Therapy Imaging, Forchheim, Germany
| |
Collapse
|
4
|
Wong YM, Koh CWY, Lew KS, Chua CGA, Yeap PL, Andrew W, Zubin M, Poh SS, Lew WS, Lee JCL, Park SY, Tan HQ. Effects of modern aesthetic dental fillings on proton therapy. Phys Imaging Radiat Oncol 2024; 29:100552. [PMID: 38405428 PMCID: PMC10891317 DOI: 10.1016/j.phro.2024.100552] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/11/2023] [Revised: 01/30/2024] [Accepted: 02/08/2024] [Indexed: 02/27/2024] Open
Abstract
Background and purpose High-density dental fillings pose a non-negligible impact on head and neck cancer treatment. For proton therapy, stopping power ratio (SPR) prediction will be significantly impaired by the associated image artifacts. Dose perturbation is also inevitable, compromising the treatment plan quality. While plenty of work has been done on metal or amalgam fillings, none has touched on composite resin (CR) and glass ionomer cement (GIC) which have seen an increasing usage. Hence, this work aims to provide a detailed characterisation of SPR and dose perturbation in proton therapy caused by CR and GIC. Materials and methods Four types of fillings were used: CR, Fuji Bulk (FB), Fuji II (FII) and Fuji IX (FIX). The latter three belong to GIC category. Measured SPR were compared with SPR predicted using single-energy computed tomography (SECT) and dual-energy computed tomography (DECT). Dose perturbation of proton beams with lower- and higher-energy levels was also quantified using Gafchromic films. Results The measured SPR for CR, FB, FII and FIX were 1.68, 1.77, 1.77 and 1.76, respectively. Overall, DECT could predict SPR better than SECT. The lowest percentage error achieved by DECT was 19.7 %, demonstrating the challenge in estimating SPR, even for fillings with relatively lower densities. For both proton beam energies and all four fillings of about 4.5 mm thickness, the maximum dose perturbation was 3 %. Conclusion This study showed that dose perturbation by CR and GIC was comparatively small. We have measured and recommended the SPR values for overriding the fillings in TPS.
Collapse
Affiliation(s)
- Yun Ming Wong
- Division of Physics and Applied Physics, Nanyang Technological University, Singapore
| | | | - Kah Seng Lew
- Division of Physics and Applied Physics, Nanyang Technological University, Singapore
- Division of Radiation Oncology, National Cancer Centre Singapore, Singapore
| | | | - Ping Lin Yeap
- Division of Radiation Oncology, National Cancer Centre Singapore, Singapore
- Department of Oncology, University of Cambridge, United Kingdom
| | - Wibawa Andrew
- Division of Radiation Oncology, National Cancer Centre Singapore, Singapore
| | - Master Zubin
- Division of Radiation Oncology, National Cancer Centre Singapore, Singapore
| | - Sharon Shuxian Poh
- Division of Radiation Oncology, National Cancer Centre Singapore, Singapore
| | - Wen Siang Lew
- Division of Physics and Applied Physics, Nanyang Technological University, Singapore
| | - James Cheow Lei Lee
- Division of Physics and Applied Physics, Nanyang Technological University, Singapore
- Division of Radiation Oncology, National Cancer Centre Singapore, Singapore
| | - Sung Yong Park
- Division of Radiation Oncology, National Cancer Centre Singapore, Singapore
- Oncology Academic Clinical Programme, Duke-NUS Medical School, Singapore
| | - Hong Qi Tan
- Division of Radiation Oncology, National Cancer Centre Singapore, Singapore
- Oncology Academic Clinical Programme, Duke-NUS Medical School, Singapore
| |
Collapse
|
5
|
Huang C, Xu Z, Zhao Z, Yin Y, Hu Z, She Q, Mao R, Wei K, Yang H, Tang K, Lu Z. Carbon ion radiography with a composite ionization chamber detector. Appl Radiat Isot 2024; 203:111072. [PMID: 37897938 DOI: 10.1016/j.apradiso.2023.111072] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/24/2023] [Revised: 10/05/2023] [Accepted: 10/16/2023] [Indexed: 10/30/2023]
Abstract
Range uncertainty in carbon ion therapy can diminish treatment efficacy because it may cause deviation from the planned dose distribution. The precise and accurate determination of relative stopping power (RSP) maps of carbon ions in the patient is a direct solution to this problem. To obtain RSP maps in patients undergoing carbon ion radiography, our team developed a preliminary prototype of a composite ionization chamber detection (CICD) system. The CICD prototype employs synchronously gated integral electronics with the ability to measure the depth-to-dose curve and the beam profile simultaneously. Carbon ion radiography experiments were performed on hemispherical, sloped, and stepped phantoms using the Heavy Ion Medical Machine (HIMM) beam. The beam energy was 190.19 MeV/μ and the beam spot full width at half maximum (FWHM) was 7.42 mm. The radiographic image of the sloped phantom, the thickness prediction accuracy of each pixel (2 mm) is 88.25%, its absolute mean error (AME) is 1.07 mm, and the maximum absolute deviation (MAD) is 2.64 mm. The prediction accuracy of the CICD prototype is mainly affected by electronic noise, with a noise-to-signal ratio (NSR) of about 14.36 dB. Carbon ion radiography simulations were performed in this study using Geant4 software to eliminate the effect of the electronic noise. The thickness prediction accuracy is 98.54%, 98.62%, and 99.07% per pixel for hemispherical, sloped and stepped phantoms, respectively, with AME of 0.09 mm, 0.27 mm, and 0.48 mm. Carbon ion radiography utilizing the CICD prototype scheme has the ability to refine the accuracy and resolution of radiographic images, consequently establishing a scientific foundation for diminishing the effects of range uncertainty and fully exploiting the advantages of precision particle therapy.
Collapse
Affiliation(s)
- Chuan Huang
- Institute of Modern Physics, Chinese Academy of Sciences, Lanzhou, 730000, China; School of Nuclear Science and Technology, Lanzhou University, Lanzhou, Gansu, 730000, China
| | - Zhiguo Xu
- Institute of Modern Physics, Chinese Academy of Sciences, Lanzhou, 730000, China.
| | - Zulong Zhao
- Institute of Modern Physics, Chinese Academy of Sciences, Lanzhou, 730000, China
| | - Yongzhi Yin
- School of Nuclear Science and Technology, Lanzhou University, Lanzhou, Gansu, 730000, China
| | - Zhengguo Hu
- Institute of Modern Physics, Chinese Academy of Sciences, Lanzhou, 730000, China
| | - Qianshun She
- Institute of Modern Physics, Chinese Academy of Sciences, Lanzhou, 730000, China
| | - Ruishi Mao
- Institute of Modern Physics, Chinese Academy of Sciences, Lanzhou, 730000, China
| | - Kun Wei
- Wuwei Occupational College, Wuwei, Gansu, 730000, China
| | - Herun Yang
- Institute of Modern Physics, Chinese Academy of Sciences, Lanzhou, 730000, China
| | - Kai Tang
- Institute of Modern Physics, Chinese Academy of Sciences, Lanzhou, 730000, China
| | - Ziwei Lu
- Institute of Modern Physics, Chinese Academy of Sciences, Lanzhou, 730000, China
| |
Collapse
|
6
|
Edvardsson A, Gorgisyan J, Andersson KM, Vallhagen Dahlgren C, Dasu A, Gram D, Björk-Eriksson T, Munck af Rosenschöld P. Robustness and dosimetric verification of hippocampal-sparing craniospinal pencil beam scanning proton plans for pediatric medulloblastoma. Phys Imaging Radiat Oncol 2024; 29:100555. [PMID: 38405431 PMCID: PMC10891325 DOI: 10.1016/j.phro.2024.100555] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/16/2023] [Revised: 02/06/2024] [Accepted: 02/08/2024] [Indexed: 02/27/2024] Open
Abstract
Background and Purpose Hippocampal-sparing (HS) is a method that can potentially reduce late cognitive complications for pediatric medulloblastoma (MB) patients treated with craniospinal proton therapy (PT). The aim of this study was to investigate robustness and dosimetric plan verification of pencil beam scanning HS PT. Materials and Methods HS and non-HS PT plans for the whole brain part of craniospinal treatment were created for 15 pediatric MB patients. A robust evaluation of the plans was performed. Plans were recalculated in a water phantom and measured field-by-field using an ion chamber detector at depths corresponding to the central part of hippocampi. All HS and non-HS fields were measured with the standard resolution of the detector and in addition 16 HS fields were measured with high resolution. Measured and planned dose distributions were compared using gamma evaluation. Results The median mean hippocampus dose was reduced from 22.9 Gy (RBE) to 8.9 Gy (RBE), while keeping CTV V95% above 95 % for all nominal HS plans. HS plans were relatively robust regarding hippocampus mean dose, however, less robust regarding target coverage and maximum dose compared to non-HS plans. For standard resolution measurements, median pass rates were 99.7 % for HS and 99.5 % for non-HS plans (p < 0.001). For high-resolution measurements, median pass rates were 100 % in the hippocampus region and 98.2 % in the surrounding region. Conclusions A substantial reduction of dose in the hippocampus region appeared feasible. Dosimetric accuracy of HS plans was comparable to non-HS plans and agreed well with planned dose distribution in the hippocampus region.
Collapse
Affiliation(s)
- Anneli Edvardsson
- Radiation Physics, Department of Hematology, Oncology and Radiation Physics, Skåne University Hospital, Sweden
- Medical Radiation Physics, Department of Clinical Sciences Lund, Lund University, Lund, Sweden
| | - Jenny Gorgisyan
- Radiation Physics, Department of Hematology, Oncology and Radiation Physics, Skåne University Hospital, Sweden
- Medical Radiation Physics, Department of Clinical Sciences Lund, Lund University, Lund, Sweden
| | | | | | - Alexandru Dasu
- The Skandion Clinic, Uppsala, Sweden
- Medical Radiation Sciences, Department of Immunology, Genetics and Pathology, Uppsala University, Uppsala, Sweden
| | - Daniel Gram
- Department of Clinical Oncology and Palliative Care, Radiotherapy, Zealand University Hospital, Næstved, Denmark
- Niels Bohr Institute, University of Copenhagen, Copenhagen, Denmark
- Department of Oncology – Section of Radiotherapy, Rigshospitalet, Copenhagen, Denmark
| | - Thomas Björk-Eriksson
- Department of Oncology, Institute of Clinical Sciences, Sahlgrenska Academy at the University of Gothenburg, Gothenburg, Sweden
- Regional Cancer Centre West, Western Sweden Healthcare Region, Gothenburg, Sweden
| | - Per Munck af Rosenschöld
- Radiation Physics, Department of Hematology, Oncology and Radiation Physics, Skåne University Hospital, Sweden
- Medical Radiation Physics, Department of Clinical Sciences Lund, Lund University, Lund, Sweden
| |
Collapse
|
7
|
Decabooter E, Hilgers GC, De Rouck J, Salvo K, Van Wingerden J, Bosmans H, van der Heyden B, Qamhiyeh S, Papalazarou C, Kaatee R, Pittomvils G, Bogaert E. Survey on fan-beam computed tomography for radiotherapy: Imaging for dose calculation and delineation. Phys Imaging Radiat Oncol 2024; 29:100522. [PMID: 38152701 PMCID: PMC10750173 DOI: 10.1016/j.phro.2023.100522] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/30/2023] [Revised: 11/26/2023] [Accepted: 11/28/2023] [Indexed: 12/29/2023] Open
Abstract
Background and purpose To obtain an understanding of current practice, professional needs and future directions in the field of fan-beam CT in RT, a survey was conducted. This work presents the collected information regarding the use of CT imaging for dose calculation and structure delineation. Materials and methods An online institutional survey was distributed to medical physics experts employed at Belgian and Dutch radiotherapy institutions to assess the status, challenges, and future directions of QA practices for fan-beam CT. A maximum of 143 questions covered topics such as CT scanner availability, CT scanner specifications, QA protocols, treatment simulation workflow, and radiotherapy dose calculation. Answer forms were collected between 1-Sep-2022 and 10-Oct-2022. Results A 66 % response rate was achieved, yielding data on a total of 58 CT scanners. For MV photon therapy, all single-energy CT scans are reconstructed in Hounsfield Units for delineation or dose calculation, and a direct- or stoichiometric method was used to convert CT numbers for dose calculation. Limited use of dual-energy CT is reported for photon (N = 3) and proton dose calculations (N = 1). For brachytherapy, most institutions adopt water-based dose calculation, while approximately 26 % of the institutions take tissue heterogeneity into account. Commissioning and regular QA include eleven tasks, which are performed by two or more professions (29/31) with varying frequencies. Conclusions Dual usage of a planning CT limits protocol optimization for both tissue characterization and delineation. DECT has been implemented only gradually. A variation of QA testing frequencies and tests are reported.
Collapse
Affiliation(s)
- Esther Decabooter
- Department of Radiation Oncology (Maastro Clinic), GROW School for Oncology, Maastricht University Medical Centre+, Maastricht, The Netherlands
| | | | - Joke De Rouck
- Department of Radiotherapy, AZ Sint Lucas, Ghent, Belgium
| | - Koen Salvo
- Department of Radiotherapy, AZ Sint-Maarten, Mechelen, Belgium
| | - Jacobus Van Wingerden
- Department of Medical Physics, Haaglanden Medical Centre, Leidschendam, The Netherlands
| | - Hilde Bosmans
- Department of Medical Radiation Physics, University Hospital Leuven, Belgium
| | - Brent van der Heyden
- IBiTech-MEDISIP, Department of Electronics and Information Systems, Ghent University, Ghent, Belgium
- Department of Oncology, Laboratory of Experimental Radiotherapy, KU Leuven, Leuven, Belgium
| | - Sima Qamhiyeh
- University Hospitals Leuven, Department of Radiation Oncology, Leuven, Belgium
| | - Chrysi Papalazarou
- Department of Radiotherapy, Leiden University Medical Center, Leiden, The Netherlands
| | - Robert Kaatee
- Radiotherapy Institute Friesland, Leeuwarden, The Netherlands
| | - Geert Pittomvils
- Department of Radiation-Oncology, Ghent University Hospital, Ghent, Belgium
| | - Evelien Bogaert
- Department of Radiation-Oncology, Ghent University Hospital, Ghent, Belgium
| |
Collapse
|
8
|
Aehle M, Alme J, Gábor Barnaföldi G, Blühdorn J, Bodova T, Borshchov V, van den Brink A, Eikeland V, Feofilov G, Garth C, Gauger NR, Grøttvik O, Helstrup H, Igolkin S, Keidel R, Kobdaj C, Kortus T, Kusch L, Leonhardt V, Mehendale S, Ningappa Mulawade R, Harald Odland O, O'Neill G, Papp G, Peitzmann T, Pettersen HES, Piersimoni P, Pochampalli R, Protsenko M, Rauch M, Ur Rehman A, Richter M, Röhrich D, Sagebaum M, Santana J, Schilling A, Seco J, Songmoolnak A, Sudár Á, Tambave G, Tymchuk I, Ullaland K, Varga-Kofarago M, Volz L, Wagner B, Wendzel S, Wiebel A, Xiao R, Yang S, Zillien S. Exploration of differentiability in a proton computed tomography simulation framework. Phys Med Biol 2023; 68:244002. [PMID: 37949060 DOI: 10.1088/1361-6560/ad0bdd] [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: 05/16/2023] [Accepted: 11/10/2023] [Indexed: 11/12/2023]
Abstract
Objective.Gradient-based optimization using algorithmic derivatives can be a useful technique to improve engineering designs with respect to a computer-implemented objective function. Likewise, uncertainty quantification through computer simulations can be carried out by means of derivatives of the computer simulation. However, the effectiveness of these techniques depends on how 'well-linearizable' the software is. In this study, we assess how promising derivative information of a typical proton computed tomography (pCT) scan computer simulation is for the aforementioned applications.Approach.This study is mainly based on numerical experiments, in which we repeatedly evaluate three representative computational steps with perturbed input values. We support our observations with a review of the algorithmic steps and arithmetic operations performed by the software, using debugging techniques.Main results.The model-based iterative reconstruction (MBIR) subprocedure (at the end of the software pipeline) and the Monte Carlo (MC) simulation (at the beginning) were piecewise differentiable. However, the observed high density and magnitude of jumps was likely to preclude most meaningful uses of the derivatives. Jumps in the MBIR function arose from the discrete computation of the set of voxels intersected by a proton path, and could be reduced in magnitude by a 'fuzzy voxels' approach. The investigated jumps in the MC function arose from local changes in the control flow that affected the amount of consumed random numbers. The tracking algorithm solves an inherently non-differentiable problem.Significance.Besides the technical challenges of merely applying AD to existing software projects, the MC and MBIR codes must be adapted to compute smoother functions. For the MBIR code, we presented one possible approach for this while for the MC code, this will be subject to further research. For the tracking subprocedure, further research on surrogate models is necessary.
Collapse
Affiliation(s)
- Max Aehle
- Chair for Scientific Computing, University of Kaiserslautern-Landau, D-67663 Kaiserslautern, Germany
| | - Johan Alme
- Department of Physics and Technology, University of Bergen, NO-5007 Bergen, Norway
| | | | - Johannes Blühdorn
- Chair for Scientific Computing, University of Kaiserslautern-Landau, D-67663 Kaiserslautern, Germany
| | - Tea Bodova
- Department of Physics and Technology, University of Bergen, NO-5007 Bergen, Norway
| | | | | | - Viljar Eikeland
- Department of Physics and Technology, University of Bergen, NO-5007 Bergen, Norway
| | | | - Christoph Garth
- Scientific Visualization Lab, University of Kaiserslautern-Landau, D-67663 Kaiserslautern, Germany
| | - Nicolas R Gauger
- Chair for Scientific Computing, University of Kaiserslautern-Landau, D-67663 Kaiserslautern, Germany
| | - Ola Grøttvik
- Department of Physics and Technology, University of Bergen, NO-5007 Bergen, Norway
| | - Håvard Helstrup
- Department of Computer Science, Electrical Engineering and Mathematical Sciences, Western Norway University of Applied Sciences, NO-5020 Bergen, Norway
| | | | - Ralf Keidel
- Chair for Scientific Computing, University of Kaiserslautern-Landau, D-67663 Kaiserslautern, Germany
- Center for Technology and Transfer (ZTT), University of Applied Sciences Worms, Worms, Germany
| | - Chinorat Kobdaj
- Institute of Science, Suranaree University of Technology, Nakhon Ratchasima, Thailand
| | - Tobias Kortus
- Center for Technology and Transfer (ZTT), University of Applied Sciences Worms, Worms, Germany
| | - Lisa Kusch
- Chair for Scientific Computing, University of Kaiserslautern-Landau, D-67663 Kaiserslautern, Germany
| | - Viktor Leonhardt
- Scientific Visualization Lab, University of Kaiserslautern-Landau, D-67663 Kaiserslautern, Germany
| | - Shruti Mehendale
- Department of Physics and Technology, University of Bergen, NO-5007 Bergen, Norway
| | - Raju Ningappa Mulawade
- Center for Technology and Transfer (ZTT), University of Applied Sciences Worms, Worms, Germany
| | - Odd Harald Odland
- Department of Physics and Technology, University of Bergen, NO-5007 Bergen, Norway
- Department of Oncology and Medical Physics, Haukeland University Hospital, NO-5021 Bergen, Norway
| | - George O'Neill
- Department of Physics and Technology, University of Bergen, NO-5007 Bergen, Norway
| | - Gábor Papp
- Institute for Physics, Eötvös Loránd University, 1/A Pázmány P. Sétány, H-1117 Budapest, Hungary
| | - Thomas Peitzmann
- Institute for Subatomic Physics, Utrecht University/Nikhef, Utrecht, Netherlands
| | | | - Pierluigi Piersimoni
- Department of Physics and Technology, University of Bergen, NO-5007 Bergen, Norway
- FSN Department, ENEA, Frascati Research Center, I-00044, Frascati, Italy
| | - Rohit Pochampalli
- Chair for Scientific Computing, University of Kaiserslautern-Landau, D-67663 Kaiserslautern, Germany
| | - Maksym Protsenko
- Research and Production Enterprise 'LTU' (RPE LTU), Kharkiv, Ukraine
| | - Max Rauch
- Department of Physics and Technology, University of Bergen, NO-5007 Bergen, Norway
| | - Attiq Ur Rehman
- Department of Physics and Technology, University of Bergen, NO-5007 Bergen, Norway
| | | | - Dieter Röhrich
- Department of Physics and Technology, University of Bergen, NO-5007 Bergen, Norway
| | - Max Sagebaum
- Chair for Scientific Computing, University of Kaiserslautern-Landau, D-67663 Kaiserslautern, Germany
| | - Joshua Santana
- Center for Technology and Transfer (ZTT), University of Applied Sciences Worms, Worms, Germany
| | - Alexander Schilling
- Center for Technology and Transfer (ZTT), University of Applied Sciences Worms, Worms, Germany
| | - Joao Seco
- Department of Biomedical Physics in Radiation Oncology, DKFZGerman Cancer Research Center, Heidelberg, Germany
- Department of Physics and Astronomy, Heidelberg University, Heidelberg, Germany
| | - Arnon Songmoolnak
- Department of Physics and Technology, University of Bergen, NO-5007 Bergen, Norway
- Institute of Science, Suranaree University of Technology, Nakhon Ratchasima, Thailand
| | - Ákos Sudár
- Wigner Research Centre for Physics, Budapest, Hungary
| | - Ganesh Tambave
- Center for Medical and Radiation Physics (CMRP), National Institute of Science Education and Research (NISER), Bhubaneswar, India
| | - Ihor Tymchuk
- Research and Production Enterprise 'LTU' (RPE LTU), Kharkiv, Ukraine
| | - Kjetil Ullaland
- Department of Physics and Technology, University of Bergen, NO-5007 Bergen, Norway
| | | | - Lennart Volz
- Biophysics, GSI Helmholtz Center for Heavy Ion Research GmbH, Darmstadt, Germany
- Department of Medical Physics and Biomedical Engineering, University College London, London, United Kingdom
| | - Boris Wagner
- Department of Physics and Technology, University of Bergen, NO-5007 Bergen, Norway
| | - Steffen Wendzel
- Center for Technology and Transfer (ZTT), University of Applied Sciences Worms, Worms, Germany
| | - Alexander Wiebel
- Center for Technology and Transfer (ZTT), University of Applied Sciences Worms, Worms, Germany
| | - RenZheng Xiao
- Department of Physics and Technology, University of Bergen, NO-5007 Bergen, Norway
- College of Mechanical & Power Engineering, China Three Gorges University, Yichang, People's Republic of China
| | - Shiming Yang
- Department of Physics and Technology, University of Bergen, NO-5007 Bergen, Norway
| | - Sebastian Zillien
- Center for Technology and Transfer (ZTT), University of Applied Sciences Worms, Worms, Germany
| |
Collapse
|
9
|
Shimomura T, Fujiwara D, Inoue Y, Takeya A, Ohta T, Nozawa Y, Imae T, Nawa K, Nakagawa K, Haga A. Virtual cone-beam computed tomography simulator with human phantom library and its application to the elemental material decomposition. Phys Med 2023; 113:102648. [PMID: 37672845 DOI: 10.1016/j.ejmp.2023.102648] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 04/18/2023] [Revised: 06/19/2023] [Accepted: 07/29/2023] [Indexed: 09/08/2023] Open
Abstract
PURPOSE The purpose of this study is to develop a virtual CBCT simulator with a head and neck (HN) human phantom library and to demonstrate the feasibility of elemental material decomposition (EMD) for quantitative CBCT imaging using this virtual simulator. METHODS The library of 36 HN human phantoms were developed by extending the ICRP 110 adult phantoms based on human age, height, and weight statistics. To create the CBCT database for the library, a virtual CBCT simulator that simulated the direct and scattered X-ray on a flat panel detector using ray-tracing and deep-learning (DL) models was used. Gaussian distributed noise was also included on the flat panel detector, which was evaluated using a real CBCT system. The usefulness of the virtual CBCT system was demonstrated through the application of the developed DL-based EMD model for case involving virtual phantom and real patient. RESULTS The virtual simulator could generate various virtual CBCT images based on the human phantom library, and the prediction of the EMD could be successfully performed by preparing the CBCT database from the proposed virtual system, even for a real patient. The CBCT image degradation owing to the scattered X-ray and the statistical noise affected the prediction accuracy, although these effects were minimal. Furthermore, the elemental distribution using the real CBCT image was also predictable. CONCLUSIONS This study demonstrated the potential of using computer vision for medical data preparation and analysis, which could have important implications for improving patient outcomes, especially in adaptive radiation therapy.
Collapse
Affiliation(s)
- Taisei Shimomura
- Graduate School of Biomedical Sciences, Tokushima University, Tokushima 770-8503, Japan; Department of Radiology, The University of Tokyo Hospital, Bunkyo, Tokyo 113-8655, Japan
| | - Daiyu Fujiwara
- Graduate School of Biomedical Sciences, Tokushima University, Tokushima 770-8503, Japan
| | - Yuki Inoue
- Graduate School of Biomedical Sciences, Tokushima University, Tokushima 770-8503, Japan
| | - Atsushi Takeya
- Graduate School of Biomedical Sciences, Tokushima University, Tokushima 770-8503, Japan
| | - Takeshi Ohta
- Department of Radiology, The University of Tokyo Hospital, Bunkyo, Tokyo 113-8655, Japan
| | - Yuki Nozawa
- Department of Radiology, The University of Tokyo Hospital, Bunkyo, Tokyo 113-8655, Japan
| | - Toshikazu Imae
- Department of Radiology, The University of Tokyo Hospital, Bunkyo, Tokyo 113-8655, Japan
| | - Kanabu Nawa
- Department of Radiology, The University of Tokyo Hospital, Bunkyo, Tokyo 113-8655, Japan
| | - Keiichi Nakagawa
- Department of Radiology, The University of Tokyo Hospital, Bunkyo, Tokyo 113-8655, Japan
| | - Akihiro Haga
- Graduate School of Biomedical Sciences, Tokushima University, Tokushima 770-8503, Japan.
| |
Collapse
|
10
|
Peters N, Trier Taasti V, Ackermann B, Bolsi A, Vallhagen Dahlgren C, Ellerbrock M, Fracchiolla F, Gomà C, Góra J, Cambraia Lopes P, Rinaldi I, Salvo K, Sojat Tarp I, Vai A, Bortfeld T, Lomax A, Richter C, Wohlfahrt P. Consensus guide on CT-based prediction of stopping-power ratio using a Hounsfield look-up table for proton therapy. Radiother Oncol 2023; 184:109675. [PMID: 37084884 PMCID: PMC10351362 DOI: 10.1016/j.radonc.2023.109675] [Citation(s) in RCA: 11] [Impact Index Per Article: 11.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/22/2022] [Revised: 03/08/2023] [Accepted: 04/10/2023] [Indexed: 04/23/2023]
Abstract
BACKGROUND AND PURPOSE Studies have shown large variations in stopping-power ratio (SPR) prediction from computed tomography (CT) across European proton centres. To standardise this process, a step-by-step guide on specifying a Hounsfield look-up table (HLUT) is presented here. MATERIALS AND METHODS The HLUT specification process is divided into six steps: Phantom setup, CT acquisition, CT number extraction, SPR determination, HLUT specification, and HLUT validation. Appropriate CT phantoms have a head- and body-sized part, with tissue-equivalent inserts in regard to X-ray and proton interactions. CT numbers are extracted from a region-of-interest covering the inner 70% of each insert in-plane and several axial CT slices in scan direction. For optimal HLUT specification, the SPR of phantom inserts is measured in a proton beam and the SPR of tabulated human tissues is computed stoichiometrically at 100 MeV. Including both phantom inserts and tabulated human tissues increases HLUT stability. Piecewise linear regressions are performed between CT numbers and SPRs for four tissue groups (lung, adipose, soft tissue, and bone) and then connected with straight lines. Finally, a thorough but simple validation is performed. RESULTS The best practices and individual challenges are explained comprehensively for each step. A well-defined strategy for specifying the connection points between the individual line segments of the HLUT is presented. The guide was tested exemplarily on three CT scanners from different vendors, proving its feasibility. CONCLUSION The presented step-by-step guide for CT-based HLUT specification with recommendations and examples can contribute to reduce inter-centre variations in SPR prediction.
Collapse
Affiliation(s)
- Nils Peters
- OncoRay - National Center for Radiation Research in Oncology, Faculty of Medicine and University Hospital Carl Gustav Carus, Technische Universität Dresden, Helmholtz-Zentrum Dresden - Rossendorf, Germany; Massachusetts General Hospital and Harvard Medical School, Department of Radiation Oncology, Boston, MA, USA.
| | - Vicki Trier Taasti
- Department of Radiation Oncology (Maastro), GROW School for Oncology and Reproduction, Maastricht University Medical Centre+, Maastricht, The Netherlands.
| | - Benjamin Ackermann
- Heidelberg Ion-Beam Therapy Center (HIT), Department of Radiation Oncology, Heidelberg University Hospital, Heidelberg, Germany
| | - Alessandra Bolsi
- Center for Proton Therapy, Paul Scherrer Institute, Villigen, Switzerland
| | | | - Malte Ellerbrock
- Heidelberg Ion-Beam Therapy Center (HIT), Department of Radiation Oncology, Heidelberg University Hospital, Heidelberg, Germany
| | - Francesco Fracchiolla
- Azienda Provinciale per i Servizi Sanitari (APSS) Protontherapy Department, Trento, Italy
| | - Carles Gomà
- Department of Radiation Oncology, Hospital Clínic de Barcelona, Barcelona, Spain
| | - Joanna Góra
- MedAustron Ion Therapy Center, Wiener Neustadt, Austria
| | | | - Ilaria Rinaldi
- Department of Radiation Oncology (Maastro), GROW School for Oncology and Reproduction, Maastricht University Medical Centre+, Maastricht, The Netherlands
| | - Koen Salvo
- AZ Sint-Maarten, Department of Radiotherapy, Mechelen, Belgium
| | - Ivanka Sojat Tarp
- Aarhus University Hospital, Danish Center for Particle Therapy, Aarhus, Denmark
| | - Alessandro Vai
- Radiotherapy Department, Center for National Oncological Hadrontherapy (CNAO), 27100 Pavia, Italy
| | - Thomas Bortfeld
- Massachusetts General Hospital and Harvard Medical School, Department of Radiation Oncology, Boston, MA, USA
| | - Antony Lomax
- Center for Proton Therapy, Paul Scherrer Institute, Villigen, Switzerland
| | - Christian Richter
- OncoRay - National Center for Radiation Research in Oncology, Faculty of Medicine and University Hospital Carl Gustav Carus, Technische Universität Dresden, Helmholtz-Zentrum Dresden - Rossendorf, Germany; Helmholtz-Zentrum Dresden - Rossendorf, Institute of Radiooncology - OncoRay, Dresden, Germany; Department of Radiotherapy and Radiation Oncology, Faculty of Medicine and University Hospital Carl Gustav Carus, Technische Universität Dresden, Dresden, Germany; German Cancer Consortium (DKTK), partner site Dresden, and German Cancer Research Center (DKFZ), Heidelberg, Germany
| | - Patrick Wohlfahrt
- Massachusetts General Hospital and Harvard Medical School, Department of Radiation Oncology, Boston, MA, USA
| |
Collapse
|
11
|
Longarino FK, Herpel C, Tessonnier T, Mein S, Ackermann B, Debus J, Schwindling FS, Stiller W, Mairani A. Dual-energy CT-based stopping power prediction for dental materials in particle therapy. J Appl Clin Med Phys 2023:e13977. [PMID: 37032540 PMCID: PMC10402687 DOI: 10.1002/acm2.13977] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/22/2022] [Revised: 01/23/2023] [Accepted: 03/17/2023] [Indexed: 04/11/2023] Open
Abstract
Radiotherapy with protons or light ions can offer accurate and precise treatment delivery. Accurate knowledge of the stopping power ratio (SPR) distribution of the tissues in the patient is crucial for improving dose prediction in patients during planning. However, materials of uncertain stoichiometric composition such as dental implant and restoration materials can substantially impair particle therapy treatment planning due to related SPR prediction uncertainties. This study investigated the impact of using dual-energy computed tomography (DECT) imaging for characterizing and compensating for commonly used dental implant and restoration materials during particle therapy treatment planning. Radiological material parameters of ten common dental materials were determined using two different DECT techniques: sequential acquisition CT (SACT) and dual-layer spectral CT (DLCT). DECT-based direct SPR predictions of dental materials via spectral image data were compared to conventional single-energy CT (SECT)-based SPR predictions obtained via indirect CT-number-to-SPR conversion. DECT techniques were found overall to reduce uncertainty in SPR predictions in dental implant and restoration materials compared to SECT, although DECT methods showed limitations for materials containing elements of a high atomic number. To assess the influence on treatment planning, an anthropomorphic head phantom with a removable tooth containing lithium disilicate as a dental material was used. The results indicated that both DECT techniques predicted similar ranges for beams unobstructed by dental material in the head phantom. When ion beams passed through the lithium disilicate restoration, DLCT-based SPR predictions using a projection-based method showed better agreement with measured reference SPR values (range deviation: 0.2 mm) compared to SECT-based predictions. DECT-based SPR prediction may improve the management of certain non-tissue dental implant and restoration materials and subsequently increase dose prediction accuracy.
Collapse
Affiliation(s)
- Friderike K Longarino
- Clinical Cooperation Unit Radiation Oncology, German Cancer Research Center (DKFZ), Heidelberg, Germany
- Department of Radiation Oncology, Heidelberg University Hospital, Heidelberg, Germany
- Department of Physics and Astronomy, Heidelberg University, Heidelberg, Germany
| | - Christopher Herpel
- Department of Prosthodontics, Heidelberg University Hospital, Heidelberg, Germany
| | - Thomas Tessonnier
- Heidelberg Ion Beam Therapy Center (HIT), Heidelberg, Germany
- Translational Radiation Oncology, German Cancer Research Center (DKFZ), Heidelberg, Germany
- Heidelberg Institute of Radiation Oncology (HIRO), National Center for Radiation Research in Oncology (NCRO), Heidelberg, Germany
- National Center for Tumor Diseases (NCT), Heidelberg, Germany
| | - Stewart Mein
- Heidelberg Ion Beam Therapy Center (HIT), Heidelberg, Germany
- Translational Radiation Oncology, German Cancer Research Center (DKFZ), Heidelberg, Germany
- Heidelberg Institute of Radiation Oncology (HIRO), National Center for Radiation Research in Oncology (NCRO), Heidelberg, Germany
- National Center for Tumor Diseases (NCT), Heidelberg, Germany
| | | | - Jürgen Debus
- Clinical Cooperation Unit Radiation Oncology, German Cancer Research Center (DKFZ), Heidelberg, Germany
- Department of Radiation Oncology, Heidelberg University Hospital, Heidelberg, Germany
- Heidelberg Ion Beam Therapy Center (HIT), Heidelberg, Germany
- Heidelberg Institute of Radiation Oncology (HIRO), National Center for Radiation Research in Oncology (NCRO), Heidelberg, Germany
- National Center for Tumor Diseases (NCT), Heidelberg, Germany
- German Cancer Consortium (DKTK), Core Center Heidelberg, Heidelberg, Germany
| | | | - Wolfram Stiller
- Diagnostic & Interventional Radiology (DIR), Heidelberg University Hospital, Heidelberg, Germany
| | - Andrea Mairani
- Department of Radiation Oncology, Heidelberg University Hospital, Heidelberg, Germany
- Heidelberg Ion Beam Therapy Center (HIT), Heidelberg, Germany
- National Center for Tumor Diseases (NCT), Heidelberg, Germany
- Medical Physics, National Center of Oncological Hadrontherapy (CNAO), Pavia, Italy
| |
Collapse
|
12
|
Yang M, Wohlfahrt P, Shen C, Bouchard H. Dual- and multi-energy CT for particle stopping-power estimation: current state, challenges and potential. Phys Med Biol 2023; 68. [PMID: 36595276 DOI: 10.1088/1361-6560/acabfa] [Citation(s) in RCA: 7] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/27/2022] [Accepted: 12/15/2022] [Indexed: 12/23/2022]
Abstract
Range uncertainty has been a key factor preventing particle radiotherapy from reaching its full physical potential. One of the main contributing sources is the uncertainty in estimating particle stopping power (ρs) within patients. Currently, theρsdistribution in a patient is derived from a single-energy CT (SECT) scan acquired for treatment planning by converting CT number expressed in Hounsfield units (HU) of each voxel toρsusing a Hounsfield look-up table (HLUT), also known as the CT calibration curve. HU andρsshare a linear relationship with electron density but differ in their additional dependence on elemental composition through different physical properties, i.e. effective atomic number and mean excitation energy, respectively. Because of that, the HLUT approach is particularly sensitive to differences in elemental composition between real human tissues and tissue surrogates as well as tissue variations within and among individual patients. The use of dual-energy CT (DECT) forρsprediction has been shown to be effective in reducing the uncertainty inρsestimation compared to SECT. The acquisition of CT data over different x-ray spectra yields additional information on the material elemental composition. Recently, multi-energy CT (MECT) has been explored to deduct material-specific information with higher dimensionality, which has the potential to further improve the accuracy ofρsestimation. Even though various DECT and MECT methods have been proposed and evaluated over the years, these approaches are still only scarcely implemented in routine clinical practice. In this topical review, we aim at accelerating this translation process by providing: (1) a comprehensive review of the existing DECT/MECT methods forρsestimation with their respective strengths and weaknesses; (2) a general review of uncertainties associated with DECT/MECT methods; (3) a general review of different aspects related to clinical implementation of DECT/MECT methods; (4) other potential advanced DECT/MECT applications beyondρsestimation.
Collapse
Affiliation(s)
- Ming Yang
- The University of Texas MD Anderson Cancer Center, Department of Radiation Physics, 1515 Holcombe Blvd Houston, TX 77030, United States of America
| | - Patrick Wohlfahrt
- Massachusetts General Hospital and Harvard Medical School, Department of Radiation Oncology, Boston, MA 02115, United States of America
| | - Chenyang Shen
- University of Texas Southwestern Medical Center, Department of Radiation Oncology, 2280 Inwood Rd Dallas, TX 75235, United States of America
| | - Hugo Bouchard
- Département de physique, Université de Montréal, Complexe des sciences, 1375 Avenue Thérèse-Lavoie-Roux, Montréal, Québec H2V0B3, Canada.,Centre de recherche du Centre hospitalier de l'Université de Montréal, 900 Rue Saint-Denis, Montréal, Québec, H2X 0A9, Canada.,Département de radio-oncologie, Centre hospitalier de l'Université de Montréal, 1051 Rue Sanguinet, Montréal, Québec H2X 3E4, Canada
| |
Collapse
|
13
|
Sueyasu S, Takayanagi T, Miyazaki K, Kuriyama Y, Ishi Y, Uesugi T, Unlu MB, Kudo N, Chen Y, Kasamatsu K, Fujii M, Kobayashi M, Rohringer W, Matsuura T. Ionoacoustic application of an optical hydrophone to detect proton beam range in water. Med Phys 2022; 50:2438-2449. [PMID: 36565440 DOI: 10.1002/mp.16189] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/08/2022] [Revised: 11/22/2022] [Accepted: 12/14/2022] [Indexed: 12/25/2022] Open
Abstract
BACKGROUND Proton range uncertainty has been the main factor limiting the ability of proton therapy to concentrate doses to tumors to their full potential. Ionoacoustic (IA) range verification is an approach to reducing this uncertainty by detecting thermoacoustic waves emitted from an irradiated volume immediately following a pulsed proton beam delivery; however, the signal weakness has been an obstacle to its clinical application. To increase the signal-to-noise ratio (SNR) with the conventional piezoelectric hydrophone (PH), the detector-sensitive volume needs to be large, but it could narrow the range of available beam angles and disturb real-time images obtained during beam delivery. PURPOSE To prevent this issue, we investigated a millimeter-sized optical hydrophone (OH) that exploits the laser interferometric principle. For two types of IA waves [γ-wave emitted from the Bragg peak (BP) and a spherical IA wave with resonant frequency (SPIRE) emitted from the gold fiducial marker (GM)], comparisons were made with PH in terms of waveforms, SNR, range detection accuracy, and signal intensity robustness against the small detector misalignment, particularly for SPIRE. METHODS A 100-MeV proton beam with a 27 ns pulse width and 4 mm beam size was produced using a fixed-field alternating gradient accelerator and was irradiated to the water phantom. The GM was set on the beam's central axis. Acrylic plates of various thicknesses, up to 12 mm, were set in front of the phantoms to shift the proton range. OH was set distal and lateral to the beam, and the range was estimated using the time-of-flight method for γ-wave and by comparing with the calibration data (SPIRE intensity versus the distance between the GM and BP) derived from an IA wave transport simulation for SPIRE. The BP dose per pulse was 0.5-0.6 Gy. To measure the variation in SPIRE amplitude against the hydrophone misalignment, the hydrophone was shifted by ± 2 mm at a maximum in lateral directions. RESULTS Despite its small size, OH could detect γ-wave with a higher SNR than the conventional PH (diameter, 29 mm), and a single measurement was sufficient to detect the beam range with a submillimeter accuracy in water. In the SPIRE measurement, OH was far more robust against the detector misalignment than the focused PH (FPH) used in our previous study [5%/mm (OH) versus 80%/mm (FPH)], and the correlation between the measured SPIRE intensity and the distance between the GM and BP agreed well with the simulation results. However, the OH sensitivity was lower than the FPH sensitivity, and about 5.6-Gy dose was required to decrease the intensity variation among measurements to less than 10%. CONCLUSION The miniature OH was found to detect weak IA signals produced by proton beams with a BP dose used in hypofractionated regimens. The OH sensitivity improvement at the MHz regime is worth exploring as the next step.
Collapse
Affiliation(s)
- Shota Sueyasu
- Graduate School of Engineering, Hokkaido University, Sapporo, Hokkaido, Japan
| | - Taisuke Takayanagi
- Hitachi Ltd, Research and Development Group, Center for Technology Innovation-Energy, Hitachi-shi, Ibaraki, Japan
| | - Koichi Miyazaki
- Hitachi Ltd, Research and Development Group, Center for Technology Innovation-Energy, Hitachi-shi, Ibaraki, Japan
| | - Yasutoshi Kuriyama
- Institute for Integrated Radiation and Nuclear Science, Kyoto University, Osaka, Japan
| | - Yoshihiro Ishi
- Institute for Integrated Radiation and Nuclear Science, Kyoto University, Osaka, Japan
| | - Tomonori Uesugi
- Institute for Integrated Radiation and Nuclear Science, Kyoto University, Osaka, Japan
| | - Mehmet Burcin Unlu
- Faculty of Engineering, Hokkaido University, Sapporo, Hokkaido, Japan.,Department of Physics, Bogazici University, Bebek, Istanbul, Turkey
| | - Nobuki Kudo
- Faculty of Information Science and Technology, Hokkaido University, Sapporo, Hokkaido, Japan
| | - Ye Chen
- Faculty of Engineering, Hokkaido University, Sapporo, Hokkaido, Japan
| | - Koki Kasamatsu
- Graduate School of Biomedical Science and Engineering, Hokkaido University, Sapporo, Hokkaido, Japan
| | | | - Masanori Kobayashi
- Planetary Exploration Research Institute, Chiba Institute of Technology, Narashino, Chiba, Japan
| | | | - Taeko Matsuura
- Faculty of Engineering, Hokkaido University, Sapporo, Hokkaido, Japan.,Proton Beam Therapy Center, Hokkaido University Hospital, Sapporo, Hokkaido, Japan
| |
Collapse
|
14
|
Hu G, Niepel K, Risch F, Kurz C, Würl M, Kröncke T, Schwarz F, Parodi K, Landry G. Assessment of quantitative information for radiation therapy at a first-generation clinical photon-counting computed tomography scanner. Front Oncol 2022; 12:970299. [PMID: 36185297 PMCID: PMC9515409 DOI: 10.3389/fonc.2022.970299] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/15/2022] [Accepted: 08/17/2022] [Indexed: 11/28/2022] Open
Abstract
As one of the latest developments in X-ray computed tomography (CT), photon-counting technology allows spectral detection, demonstrating considerable advantages as compared to conventional CT. In this study, we investigated the use of a first-generation clinical photon-counting computed tomography (PCCT) scanner and estimated proton relative (to water) stopping power (RSP) of tissue-equivalent materials from virtual monoenergetic reconstructions provided by the scanner. A set of calibration and evaluation tissue-equivalent inserts were scanned at 120 kVp. Maps of relative electron density (RED) and effective atomic number (EAN) were estimated from the reconstructed virtual monoenergetic images (VMI) using an approach previously applied to a spectral CT scanner with dual-layer detector technology, which allows direct calculation of RSP using the Bethe-Bloch formula. The accuracy of RED, EAN, and RSP was evaluated by root-mean-square errors (RMSE) averaged over the phantom inserts. The reference RSP values were obtained experimentally using a water column in an ion beam. For RED and EAN, the reference values were calculated based on the mass density and the chemical composition of the inserts. Different combinations of low- and high-energy VMIs were investigated in this study, ranging from 40 to 190 keV. The overall lowest error was achieved using VMIs at 60 and 180 keV, with an RSP accuracy of 1.27% and 0.71% for the calibration and the evaluation phantom, respectively.
Collapse
Affiliation(s)
- Guyue Hu
- Department of Medical Physics, Faculty of Physics, Ludwig-Maximilians-Universität München (LMU), Garching bei München, Germany
- *Correspondence: Guyue Hu,
| | - Katharina Niepel
- Department of Medical Physics, Faculty of Physics, Ludwig-Maximilians-Universität München (LMU), Garching bei München, Germany
| | - Franka Risch
- Department of Diagnostic and Interventional Radiology, Universitätsklinikum Augsburg, Augsburg, Germany
| | | | - Matthias Würl
- Department of Medical Physics, Faculty of Physics, Ludwig-Maximilians-Universität München (LMU), Garching bei München, Germany
| | - Thomas Kröncke
- Department of Diagnostic and Interventional Radiology, Universitätsklinikum Augsburg, Augsburg, Germany
| | - Florian Schwarz
- Department of Diagnostic and Interventional Radiology, Universitätsklinikum Augsburg, Augsburg, Germany
- Medical Faculty, Ludwig-Maximilians-Universität München, Munich, Germany
| | - Katia Parodi
- Department of Medical Physics, Faculty of Physics, Ludwig-Maximilians-Universität München (LMU), Garching bei München, Germany
| | - Guillaume Landry
- Department of Medical Physics, Faculty of Physics, Ludwig-Maximilians-Universität München (LMU), Garching bei München, Germany
- Department of Radiation Oncology, LMU Klinikum, Munich, Germany
- German Cancer Consortium (DKTK), Partner Site Munich, Munich, Germany
| |
Collapse
|
15
|
Fujiwara D, Shimomura T, Zhao W, Li KW, Haga A, Geng LS. Virtual computed-tomography system for deep-learning-based material decomposition. Phys Med Biol 2022; 67. [DOI: 10.1088/1361-6560/ac7bcd] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/19/2022] [Accepted: 06/23/2022] [Indexed: 11/12/2022]
Abstract
Abstract
Objective. Material decomposition (MD) evaluates the elemental composition of human tissues and organs via computed tomography (CT) and is indispensable in correlating anatomical images with functional ones. A major issue in MD is inaccurate elemental information about the real human body. To overcome this problem, we developed a virtual CT system model, by which various reconstructed images can be generated based on ICRP110 human phantoms with information about six major elements (H, C, N, O, P, and Ca). Approach. We generated CT datasets labelled with accurate elemental information using the proposed generative CT model and trained a deep learning (DL)-based model to estimate the material distribution with the ICRP110 based human phantom as well as the digital Shepp–Logan phantom. The accuracy in quad-, dual-, and single-energy CT cases was investigated. The influence of beam-hardening artefacts, noise, and spectrum variations were analysed with testing datasets including elemental density and anatomical shape variations. Main results. The results indicated that this DL approach can realise precise MD, even with single-energy CT images. Moreover, noise, beam-hardening artefacts, and spectrum variations were shown to have minimal impact on the MD. Significance. Present results suggest that the difficulty to prepare a large CT database can be solved by introducing the virtual CT system and the proposed technique can be applied to clinical radiodiagnosis and radiotherapy.
Collapse
|
16
|
In vivo assessment of tissue-specific radiological parameters with intra- and inter-patient variation using dual-energy computed tomography. Radiother Oncol 2022; 175:34-41. [DOI: 10.1016/j.radonc.2022.07.021] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/28/2022] [Revised: 07/14/2022] [Accepted: 07/25/2022] [Indexed: 11/19/2022]
|
17
|
Chacko MS, Wu D, Grewal HS, Sonnad JR. Impact of beam-hardening corrections on proton relative stopping power estimates from single- and dual-energy CT. J Appl Clin Med Phys 2022; 23:e13711. [PMID: 35816460 PMCID: PMC9512361 DOI: 10.1002/acm2.13711] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/21/2022] [Revised: 05/19/2022] [Accepted: 06/14/2022] [Indexed: 12/30/2022] Open
Abstract
A major contributing factor to proton range uncertainty is the conversion of computed tomography (CT) Hounsfield units (HU) to proton relative stopping power (RSP). This uncertainty is heightened in the presence of X-ray beam-hardening artifact (BHA), which has two manifestations: cupping and streaking, especially in and near bone tissue. This uncertainty can affect the accuracy of proton RSP calculation for treatment planning in proton radiotherapy. Dual-energy CT (DECT) and iterative beam-hardening correction (iBHC) both show promise in mitigating CT BHA. This present work attempts to analyze the relative robustness of iBHC and DECT techniques on both manifestations of BHA. The stoichiometric method for HU to RSP conversion was used for single-energy CT (SECT) and DECT-based monochromatic techniques using a tissue substitute phantom. Cupping BHA was simulated by measuring the HU of a bone substitute plug in wax/3D-printed phantoms of increasing size. Streaking BHA was simulated by placing a solid water plug between two bone plugs in a wax phantom. Finally, the effect of varying calibration phantom size on RSP was calculated in an anthropomorphic head phantom. The RSP decreased -0.002 cm-1 as phantom size increased for SECT but remained largely constant when iBHC applied or with DECT techniques. The RSP varied a maximum of 2.60% in the presence of streaking BHA in SECT but was reduced to 1.40% with iBHC. For DECT techniques, the maximum difference was 2.40%, reduced to 0.6% with iBHC. Comparing calibration phantoms of 20- and 33-cm diameter, maximum voxel differences of 5 mm in the water-equivalent thickness were observed in the skull but reduced to 1.3 mm with iBHC. The DECT techniques excelled in mitigating cupping BHA, but streaking BHA still could be observed. The use of iBHC reduced RSP variation with BHA in both SECT and DECT techniques.
Collapse
Affiliation(s)
- Michael S. Chacko
- Department of Medical Physics and DosimetryOklahoma Proton CenterOklahoma CityOklahomaUSA,Department of Radiological SciencesUniversity of Oklahoma Health Sciences CenterOklahoma CityOklahomaUSA
| | - Dee Wu
- Department of Radiological SciencesUniversity of Oklahoma Health Sciences CenterOklahoma CityOklahomaUSA
| | - Hardev S. Grewal
- Department of Radiation OncologyUniversity of Florida College of MedicineGainesvilleFloridaUSA,Department of Radiation OncologyUniversity of Florida Health Proton Therapy InstituteJacksonvilleFloridaUSA
| | - Jagadeesh R. Sonnad
- Department of Radiological SciencesUniversity of Oklahoma Health Sciences CenterOklahoma CityOklahomaUSA
| |
Collapse
|
18
|
van der Heyden B, Heymans SV, Carlier B, Collado-Lara G, Sterpin E, D’hooge J. Deep learning for dose assessment in radiotherapy by the super-localization of vaporized nanodroplets in high frame rate ultrasound imaging. Phys Med Biol 2022; 67. [DOI: 10.1088/1361-6560/ac6cc3] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/21/2022] [Accepted: 05/04/2022] [Indexed: 11/12/2022]
Abstract
Abstract
Objective. External beam radiotherapy is aimed to precisely deliver a high radiation dose to malignancies, while optimally sparing surrounding healthy tissues. With the advent of increasingly complex treatment plans, the delivery should preferably be verified by quality assurance methods. Recently, online ultrasound imaging of vaporized radiosensitive nanodroplets was proposed as a promising tool for in vivo dosimetry in radiotherapy. Previously, the detection of sparse vaporization events was achieved by applying differential ultrasound (US) imaging followed by intensity thresholding using subjective parameter tuning, which is sensitive to image artifacts. Approach. A generalized deep learning solution (i.e. BubbleNet) is proposed to localize vaporized nanodroplets on differential US frames, while overcoming the aforementioned limitation. A 5-fold cross-validation was performed on a diversely composed 5747-frame training/validation dataset by manual segmentation. BubbleNet was then applied on a test dataset of 1536 differential US frames to evaluate dosimetric features. The intra-observer variability was determined by scoring the Dice similarity coefficient (DSC) on 150 frames segmented twice. Additionally, the BubbleNet generalization capability was tested on an external test dataset of 432 frames acquired by a phased array transducer at a much lower ultrasound frequency and reconstructed with unconventional pixel dimensions with respect to the training dataset. Main results. The median DSC in the 5-fold cross validation was equal to ∼0.88, which was in line with the intra-observer variability (=0.86). Next, BubbleNet was employed to detect vaporizations in differential US frames obtained during the irradiation of phantoms with a 154 MeV proton beam or a 6 MV photon beam. BubbleNet improved the bubble-count statistics by ∼30% compared to the earlier established intensity-weighted thresholding. The proton range was verified with a −0.8 mm accuracy. Significance. BubbleNet is a flexible tool to localize individual vaporized nanodroplets on experimentally acquired US images, which improves the sensitivity compared to former thresholding-weighted methods.
Collapse
|
19
|
Longarino FK, Kowalewski A, Tessonnier T, Mein S, Ackermann B, Debus J, Mairani A, Stiller W. Potential of a Second-Generation Dual-Layer Spectral CT for Dose Calculation in Particle Therapy Treatment Planning. Front Oncol 2022; 12:853495. [PMID: 35530308 PMCID: PMC9069208 DOI: 10.3389/fonc.2022.853495] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/12/2022] [Accepted: 03/14/2022] [Indexed: 11/20/2022] Open
Abstract
In particle therapy treatment planning, dose calculation is conducted using patient-specific maps of tissue ion stopping power ratio (SPR) to predict beam ranges. Improving patient-specific SPR prediction is therefore essential for accurate dose calculation. In this study, we investigated the use of the Spectral CT 7500, a second-generation dual-layer spectral computed tomography (DLCT) system, as an alternative to conventional single-energy CT (SECT) for patient-specific SPR prediction. This dual-energy CT (DECT)-based method allows for the direct prediction of SPR from quantitative measurements of relative electron density and effective atomic number using the Bethe equation, whereas the conventional SECT-based method consists of indirect image data-based prediction through the conversion of calibrated CT numbers to SPR. The performance of the Spectral CT 7500 in particle therapy treatment planning was characterized by conducting a thorough analysis of its SPR prediction accuracy for both tissue-equivalent materials and common non-tissue implant materials. In both instances, DLCT was found to reduce uncertainty in SPR predictions compared to SECT. Mean deviations of 0.7% and 1.6% from measured SPR values were found for DLCT- and SECT-based predictions, respectively, in tissue-equivalent materials. Furthermore, end-to-end analyses of DLCT-based treatment planning were performed for proton, helium, and carbon ion therapies with anthropomorphic head and pelvic phantoms. 3D gamma analysis was performed with ionization chamber array measurements as the reference. DLCT-predicted dose distributions revealed higher passing rates compared to SECT-predicted dose distributions. In the DLCT-based treatment plans, measured distal-edge evaluation layers were within 1 mm of their predicted positions, demonstrating the accuracy of DLCT-based particle range prediction. This study demonstrated that the use of the Spectral CT 7500 in particle therapy treatment planning may lead to better agreement between planned and delivered dose compared to current clinical SECT systems.
Collapse
Affiliation(s)
- Friderike K Longarino
- Clinical Cooperation Unit Radiation Oncology, German Cancer Research Center (DKFZ), Heidelberg, Germany.,Department of Radiation Oncology, Heidelberg University Hospital, Heidelberg, Germany.,Department of Physics and Astronomy, Heidelberg University, Heidelberg, Germany
| | - Antonia Kowalewski
- Department of Radiation Oncology, Heidelberg University Hospital, Heidelberg, Germany.,Translational Radiation Oncology, German Cancer Research Center (DKFZ), Heidelberg, Germany.,Department of Physics, Simon Fraser University, Burnaby, BC, Canada
| | | | - Stewart Mein
- Department of Radiation Oncology, Heidelberg University Hospital, Heidelberg, Germany.,Translational Radiation Oncology, German Cancer Research Center (DKFZ), Heidelberg, Germany.,Heidelberg Ion Beam Therapy Center (HIT), Heidelberg, Germany.,Heidelberg Institute of Radiation Oncology (HIRO), National Center for Radiation Research in Oncology (NCRO), Heidelberg, Germany.,National Center for Tumor Diseases (NCT), Heidelberg, Germany
| | | | - Jürgen Debus
- Clinical Cooperation Unit Radiation Oncology, German Cancer Research Center (DKFZ), Heidelberg, Germany.,Department of Radiation Oncology, Heidelberg University Hospital, Heidelberg, Germany.,Heidelberg Ion Beam Therapy Center (HIT), Heidelberg, Germany.,Heidelberg Institute of Radiation Oncology (HIRO), National Center for Radiation Research in Oncology (NCRO), Heidelberg, Germany.,National Center for Tumor Diseases (NCT), Heidelberg, Germany.,German Cancer Consortium (DKTK), Core Center Heidelberg, Heidelberg, Germany
| | - Andrea Mairani
- Department of Radiation Oncology, Heidelberg University Hospital, Heidelberg, Germany.,Heidelberg Ion Beam Therapy Center (HIT), Heidelberg, Germany.,National Center for Tumor Diseases (NCT), Heidelberg, Germany.,Medical Physics, National Center of Oncological Hadrontherapy (CNAO), Pavia, Italy
| | - Wolfram Stiller
- Diagnostic and Interventional Radiology (DIR), Heidelberg University Hospital, Heidelberg, Germany
| |
Collapse
|
20
|
Poon DMC, Wu S, Ho L, Cheung KY, Yu B. Proton Therapy for Prostate Cancer: Challenges and Opportunities. Cancers (Basel) 2022; 14:cancers14040925. [PMID: 35205673 PMCID: PMC8870339 DOI: 10.3390/cancers14040925] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/20/2021] [Revised: 02/10/2022] [Accepted: 02/11/2022] [Indexed: 01/02/2023] Open
Abstract
Simple Summary Reported clinical outcomes of proton therapy (PT) for localized prostate cancer are similar to photon-based external beam radiotherapy. Apparently, the dosimetric advantages of PT have yet to be translated to clinical benefits. The suboptimal clinical outcomes of PT might be attributable to inadequate dose prescription, as indicated by the ASCENDE-RT trial. Moreover, uncertainties involved in the treatment planning and delivery processes, as well as technological limitations in PT treatment systems, may lead to discrepancies between planned doses and actual doses delivered to patients. In this article, we reviewed the current status of PT for prostate cancer and discussed different clinical implementations that could potentially improve the clinical outcome of PT for prostate cancer. Various technological advancements under which uncertainties in dose calculations can be minimized, including MRI-guided PT, dual-energy photon-counting CT and high-resolution Monte Carlo-based treatment planning systems, are highlighted. Abstract The dosimetric advantages of proton therapy (PT) treatment plans are demonstrably superior to photon-based external beam radiotherapy (EBRT) for localized prostate cancer, but the reported clinical outcomes are similar. This may be due to inadequate dose prescription, especially in high-risk disease, as indicated by the ASCENDE-RT trial. Alternatively, the lack of clinical benefits with PT may be attributable to improper dose delivery, mainly due to geometric and dosimetric uncertainties during treatment planning, as well as delivery procedures that compromise the dose conformity of treatments. Advanced high-precision PT technologies, and treatment planning and beam delivery techniques are being developed to address these uncertainties. For instance, external magnetic resonance imaging (MRI)-guided patient setup rooms are being developed to improve the accuracy of patient positioning for treatment. In-room MRI-guided patient positioning systems are also being investigated to improve the geometric accuracy of PT. Soon, high-dose rate beam delivery systems will shorten beam delivery time to within one breath hold, minimizing the effects of organ motion and patient movements. Dual-energy photon-counting computed tomography and high-resolution Monte Carlo-based treatment planning systems are available to minimize uncertainties in dose planning calculations. Advanced in-room treatment verification tools such as prompt gamma detector systems will be used to verify the depth of PT. Clinical implementation of these new technologies is expected to improve the accuracy and dose conformity of PT in the treatment of localized prostate cancers, and lead to better clinical outcomes. Improvement in dose conformity may also facilitate dose escalation, improving local control and implementation of hypofractionation treatment schemes to improve patient throughput and make PT more cost effective.
Collapse
Affiliation(s)
- Darren M. C. Poon
- Comprehensive Oncology Centre, Hong Kong Sanatorium & Hospital, Hong Kong 999077, China;
| | - Stephen Wu
- Medical Physics Department, Hong Kong Sanatorium & Hospital, Hong Kong 999077, China; (L.H.); (K.Y.C.); (B.Y.)
- Correspondence: ; Tel.: +852-29171413
| | - Leon Ho
- Medical Physics Department, Hong Kong Sanatorium & Hospital, Hong Kong 999077, China; (L.H.); (K.Y.C.); (B.Y.)
| | - Kin Yin Cheung
- Medical Physics Department, Hong Kong Sanatorium & Hospital, Hong Kong 999077, China; (L.H.); (K.Y.C.); (B.Y.)
| | - Ben Yu
- Medical Physics Department, Hong Kong Sanatorium & Hospital, Hong Kong 999077, China; (L.H.); (K.Y.C.); (B.Y.)
| |
Collapse
|
21
|
Medrano M, Liu R, Zhao T, Webb T, Politte DG, Whiting BR, Liao R, Ge T, Porras-Chaverri MA, O'Sullivan JA, Williamson JF. Towards sub-percentage uncertainty proton stopping-power mapping via dual-energy CT: direct experimental validation and uncertainty analysis of a statistical iterative image reconstruction method. Med Phys 2022; 49:1599-1618. [PMID: 35029302 DOI: 10.1002/mp.15457] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/28/2021] [Revised: 10/28/2021] [Accepted: 12/22/2021] [Indexed: 11/06/2022] Open
Abstract
PURPOSE To assess the potential of a joint dual-energy CT reconstruction process (Statistical Image Reconstruction method built on a Basis Vector Model (JSIR-BVM)) implemented on a 16-slice commercial CT scanner to measure high spatial-resolution stopping-power ratio (SPR) maps with uncertainties of less than 1%. METHODS JSIR-BVM was used to reconstruct images of effective electron density and mean excitation energy from dual-energy CT (DECT) sinograms for ten high-purity samples of known density and atomic composition inserted into head and body phantoms. The measured DECT data consisted of 90 kVp and 140 kVp axial sinograms serially acquired on a Philips Brilliance Big Bore CT scanner without beam-hardening corrections. The corresponding SPRs were subsequently measured directly via ion chamber measurements on a MEVION S250 superconducting synchrocyclotron and evaluated theoretically from the known sample compositions and densities. Deviations of JSIR-BVM SPR values from their theoretically calculated and directly measured ground-truth values were evaluated for our JSIR-BVM method and for our implementation of the Hünemohr-Saito (H-S) DECT image-domain decomposition technique for SPR imaging. A thorough uncertainty analysis was then performed for 5 different scenarios (comparison of JSIR-BVM SPR/SP to International Commission on Radiation Measurements and Units (ICRU) benchmarks; comparison of JSIR-BVM SPR to measured benchmarks; and uncertainties in JSIR-BVM SPR/SP maps for patients of unknown composition) per the Joint Committee for Guides in Metrology (JCGM) and the Guide to expression of Uncertainty in Measurement (GUM), including the impact of uncertainties in measured photon spectra, sample composition and density, photon cross-section and I-value models, and random measurement uncertainty. Estimated SPR uncertainty for three main tissue groups in patients of unknown composition and the weighted proportion of each tissue type for three proton treatment sites were then used to derive a composite range uncertainty for our method. RESULTS Mean JSIR-BVM SPR estimates deviated by less than 1% from their theoretical and directly measured ground-truth values for most inserts and phantom geometries except for high density Delrin and Teflon samples with SPR error relative to proton measurements of 1.1% and -1.0% (Head Phantom) and 1.1% and -1.1% (Body Phantom). The overall RMS deviations over all samples were 0.39% and 0.52% (head phantom) and 0.43% and 0.57% (body phantom) relative to theoretical and directly measured ground-truth SPRs, respectively. The corresponding RMS (maximum) errors for the image-domain decomposition method were 2.68% and 2.73% (4.68% and 4.99%) for the head phantom and 0.71% and 0.87% (1.37% and 1.66%) for the body phantom. Compared to H-S SPR maps, JSIR-BVM yielded 30% sharper and two-fold sharper images for soft tissues and bone-like surrogates, respectively, while reducing noise by factors of 6 and 3, respectively. The uncertainty (coverage factor k = 1) of the DECT-to-benchmark values comparison ranged from 0.5% to 1.5% and is dominated by scanning-beam photon-spectra uncertainties. An analysis of the SPR uncertainty for patients of unknown composition showed a JSIR-BVM uncertainty of 0.65%, 1.21%, and 0.77% for soft-, lung-, and bony-tissue groups which led to a composite range uncertainty of 0.6%-0.9%. CONCLUSIONS Observed JSIR-BVM SPR estimation errors were all less than 50% of the estimated k = 1 total uncertainty of our benchmarking experiment, demonstrating that JSIR-BVM high spatial-resolution, low-noise SPR mapping is feasible and is robust to variations in the geometry of the scanned object. In contrast, the much larger H-S SPR estimation errors are dominated by imaging noise and residual beam-hardening artifacts. While the uncertainties characteristic of our current JSIR-BVM implementation can be as large as 1.5%, achieving <1% total uncertainty is feasible by improving the accuracy of scanner-specific scatter-profile and photon-spectrum estimates. With its robustness to beam-hardening artifact, image noise and variations in phantom size and geometry, JSIR-BVM has the potential to achieve high spatial-resolution SPR mapping with sub-percentage accuracy and estimated uncertainty in the clinical setting. This article is protected by copyright. All rights reserved.
Collapse
Affiliation(s)
- Maria Medrano
- Department of Electrical and Systems Engineering, Washington University in St. Louis, St. Louis, MO, 63130, USA
| | - Ruirui Liu
- Department of Radiation Oncology and Winship Cancer Institute, Emory University, Atlanta, GA, 30322, USA
| | - Tianyu Zhao
- Department of Radiation Oncology, Washington University in St. Louis, St. Louis, MO, 63110, USA
| | - Tyler Webb
- Department of Physics, Washington University in St. Louis, St. Louis, MO, 63130, USA
| | - David G Politte
- Mallinckrodt Institute of Radiology, St. Louis, MO, 63110, USA
| | - Bruce R Whiting
- Department of Radiology, University of Pittsburgh, Pittsburgh, PA, 15213, USA
| | - Rui Liao
- Department of Electrical and Systems Engineering, Washington University in St. Louis, St. Louis, MO, 63130, USA
| | - Tao Ge
- Department of Electrical and Systems Engineering, Washington University in St. Louis, St. Louis, MO, 63130, USA
| | - Mariela A Porras-Chaverri
- Atomic, Nuclear and Molecular Sciences Research Center (CICANUM), University of Costa Rica, San Jose, Costa Rica
| | - Joseph A O'Sullivan
- Department of Electrical and Systems Engineering, Washington University in St. Louis, St. Louis, MO, 63130, USA
| | - Jeffrey F Williamson
- Department of Radiation Oncology, Washington University in St. Louis, St. Louis, MO, 63110, USA
| |
Collapse
|
22
|
Longarino FK, Tessonnier T, Mein S, Harrabi SB, Debus J, Stiller W, Mairani A. Dual-layer spectral CT for proton, helium, and carbon ion beam therapy planning of brain tumors. J Appl Clin Med Phys 2022; 23:e13465. [PMID: 34724327 PMCID: PMC8803296 DOI: 10.1002/acm2.13465] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/23/2021] [Revised: 08/23/2021] [Accepted: 10/14/2021] [Indexed: 01/21/2023] Open
Abstract
Pretreatment computed tomography (CT) imaging is an essential component of the particle therapy treatment planning chain. Treatment planning and optimization with charged particles require accurate and precise estimations of ion beam range in tissues, characterized by the stopping power ratio (SPR). Reduction of range uncertainties arising from conventional CT-number-to-SPR conversion based on single-energy CT (SECT) imaging is of importance for improving clinical practice. Here, the application of a novel imaging and computational methodology using dual-layer spectral CT (DLCT) was performed toward refining patient-specific SPR estimates. A workflow for DLCT-based treatment planning was devised to evaluate SPR prediction for proton, helium, and carbon ion beam therapy planning in the brain. DLCT- and SECT-based SPR predictions were compared in homogeneous and heterogeneous anatomical regions. This study included eight patients scanned for diagnostic purposes with a DLCT scanner. For each patient, four different treatment plans were created, simulating tumors in different parts of the brain. For homogeneous anatomical regions, mean SPR differences of about 1% between the DLCT- and SECT-based approaches were found. In plans of heterogeneous anatomies, relative (absolute) proton range shifts of 0.6% (0.4 mm) in the mean and up to 4.4% (2.1 mm) at the distal fall-off were observed. In the investigated cohort, 12% of the evaluated organs-at-risk (OARs) presented differences in mean or maximum dose of more than 0.5 Gy (RBE) and up to 6.8 Gy (RBE) over the entire treatment. Range shifts and dose differences in OARs between DLCT and SECT in helium and carbon ion treatment plans were similar to protons. In the majority of investigated cases (75th percentile), SECT- and DLCT-based range estimations were within 0.6 mm. Nonetheless, the magnitude of patient-specific range deviations between SECT and DLCT was clinically relevant in heterogeneous anatomical sites, suggesting further study in larger, more diverse cohorts. Results indicate that patients with brain tumors may benefit from DLCT-based treatment planning.
Collapse
Affiliation(s)
- Friderike K. Longarino
- German Cancer Research Center (DKFZ)Clinical Cooperation Unit Radiation OncologyHeidelbergGermany
- Department of Radiation OncologyHeidelberg University HospitalHeidelbergGermany
- Department of Physics and AstronomyHeidelberg UniversityHeidelbergGermany
| | | | - Stewart Mein
- Department of Radiation OncologyHeidelberg University HospitalHeidelbergGermany
- Heidelberg Ion Beam Therapy Center (HIT)HeidelbergGermany
- German Cancer Research Center (DKFZ)Translational Radiation OncologyHeidelbergGermany
- National Center for Radiation Research in Oncology (NCRO)Heidelberg Institute of Radiation Oncology (HIRO)HeidelbergGermany
- National Center for Tumor Diseases (NCT)HeidelbergGermany
| | - Semi B. Harrabi
- Department of Radiation OncologyHeidelberg University HospitalHeidelbergGermany
- Heidelberg Ion Beam Therapy Center (HIT)HeidelbergGermany
- National Center for Radiation Research in Oncology (NCRO)Heidelberg Institute of Radiation Oncology (HIRO)HeidelbergGermany
- National Center for Tumor Diseases (NCT)HeidelbergGermany
| | - Jürgen Debus
- German Cancer Research Center (DKFZ)Clinical Cooperation Unit Radiation OncologyHeidelbergGermany
- Department of Radiation OncologyHeidelberg University HospitalHeidelbergGermany
- Heidelberg Ion Beam Therapy Center (HIT)HeidelbergGermany
- National Center for Radiation Research in Oncology (NCRO)Heidelberg Institute of Radiation Oncology (HIRO)HeidelbergGermany
- National Center for Tumor Diseases (NCT)HeidelbergGermany
- Partner Site HeidelbergGerman Cancer Consortium (DKTK)HeidelbergGermany
| | - Wolfram Stiller
- Diagnostic and Interventional Radiology (DIR)Heidelberg University HospitalHeidelbergGermany
| | - Andrea Mairani
- Department of Radiation OncologyHeidelberg University HospitalHeidelbergGermany
- Heidelberg Ion Beam Therapy Center (HIT)HeidelbergGermany
- National Center for Tumor Diseases (NCT)HeidelbergGermany
- Medical PhysicsNational Centre of Oncological Hadrontherapy (CNAO)PaviaItaly
| |
Collapse
|
23
|
Volz L, Collins-Fekete CA, Bär E, Brons S, Graeff C, Johnson RP, Runz A, Sarosiek C, Schulte RW, Seco J. The accuracy of helium ion CT based particle therapy range prediction: an experimental study comparing different particle and x-ray CT modalities. Phys Med Biol 2021; 66:10.1088/1361-6560/ac33ec. [PMID: 34706355 PMCID: PMC8792995 DOI: 10.1088/1361-6560/ac33ec] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/06/2021] [Accepted: 10/27/2021] [Indexed: 11/12/2022]
Abstract
This work provides a quantitative assessment of helium ion CT (HeCT) for particle therapy treatment planning. For the first time, HeCT based range prediction accuracy in a heterogeneous tissue phantom is presented and compared to single-energy x-ray CT (SECT), dual-energy x-ray CT (DECT) and proton CT (pCT). HeCT and pCT scans were acquired using the US pCT collaboration prototype particle CT scanner at the Heidelberg Ion-Beam Therapy Center. SECT and DECT scans were done with a Siemens Somatom Definition Flash and converted to RSP. A Catphan CTP404 module was used to study the RSP accuracy of HeCT. A custom phantom of 20 cm diameter containing several tissue equivalent plastic cubes was used to assess the spatial resolution of HeCT and compare it to DECT. A clinically realistic heterogeneous tissue phantom was constructed using cranial slices from a pig head placed inside a cylindrical phantom (ø150 mm). A proton beam (84.67 mm range) depth-dose measurement was acquired using a stack of GafchromicTM EBT-XD films in a central dosimetry insert in the phantom. CT scans of the phantom were acquired with each modality, and proton depth-dose estimates were simulated based on the reconstructions. The RSP accuracy of HeCT for the plastic phantom was found to be 0.3 ± 0.1%. The spatial resolution for HeCT of the cube phantom was 5.9 ± 0.4 lp cm-1for central, and 7.6 ± 0.8 lp cm-1for peripheral cubes, comparable to DECT spatial resolution (7.7 ± 0.3 lp cm-1and 7.4 ± 0.2 lp cm-1, respectively). For the pig head, HeCT, SECT, DECT and pCT predicted range accuracy was 0.25%, -1.40%, -0.45% and 0.39%, respectively. In this study, HeCT acquired with a prototype system showed potential for particle therapy treatment planning, offering RSP accuracy, spatial resolution, and range prediction accuracy comparable to that achieved with a commercial DECT scanner. Still, technical improvements of HeCT are needed to enable clinical implementation.
Collapse
Affiliation(s)
- L Volz
- Department of Biomedical Physics in Radiation Oncology, Deutsches Krebsforschungszentrum (DKFZ), Heidelberg, Germany
- Department of Physics and Astronomy, Heidelberg University, Heidelberg, Germany
| | - C-A Collins-Fekete
- Department of Medical Physics and Biomedical Engineering, University College London, London, United Kingdom
| | - E Bär
- Department of Medical Physics and Biomedical Engineering, University College London, London, United Kingdom
- Department of Radiotherapy Physics, University College London Hospitals NHS Foundation Trust, London, United Kingdom
| | - S Brons
- Heidelberg Ion-Beam Therapy Center, Universitäts Klinikum Heidelberg, Heidelberg, Germany
| | - C Graeff
- Biophysics, GSI Helmholtz Center for Heavy Ion Research GmbH, Darmstadt, Germany
| | - R P Johnson
- Department of Physics, University of California at Santa Cruz, Santa Cruz, United States of America
| | - A Runz
- Department of Medical Physics in Radiation Therapy, Deutsches Krebsforschungszentrum (DKFZ), Heidelberg, Germany
| | - C Sarosiek
- Department of Physics, Northern Illinois University, DeKalb, United States of America
| | - R W Schulte
- Department of Basic Sciences, Division of Biomedical Engineering Sciences, Loma Linda University, Loma Linda, United States of America
| | - J Seco
- Department of Biomedical Physics in Radiation Oncology, Deutsches Krebsforschungszentrum (DKFZ), Heidelberg, Germany
- Department of Physics and Astronomy, Heidelberg University, Heidelberg, Germany
| |
Collapse
|
24
|
Reduction of clinical safety margins in proton therapy enabled by the clinical implementation of dual-energy CT for direct stopping-power prediction. Radiother Oncol 2021; 166:71-78. [PMID: 34774653 DOI: 10.1016/j.radonc.2021.11.002] [Citation(s) in RCA: 43] [Impact Index Per Article: 14.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/19/2021] [Revised: 08/24/2021] [Accepted: 11/02/2021] [Indexed: 12/31/2022]
Abstract
PURPOSE To quantifiy the range uncertainty in proton treatment planning using dual-energy computed tomography (DECT) for a direct stopping-power prediction (DirectSPR) algorithm and its clinical implementation. METHODS AND MATERIALS To assess the overall uncertainty in stopping-power ratio (SPR) prediction of a DirectSPR implementation calibrated for different patient geometries, the influencing factors were categorized in imaging, modeling as well as others. The respective SPR uncertainty was quantified for lung, soft tissue and bone and translated into range uncertainty for several tumor types. The amount of healthy tissue spared was quantified for 250 patients treated with DirectSPR and the dosimetric impact was evaluated exemplarily for a representative brain-tumor patient. RESULTS For bone, soft tissue and lung, an SPR uncertainty (1σ) of 1.6%, 1.3% and 1.3% was determined for DirectSPR, respectively. This allowed for a reduction of the clinically applied range uncertainty from currently (3.5% + 2 mm) to (1.7% + 2 mm) for brain-tumor and (2.0% + 2 mm) for prostate-cancer patients. The 150 brain-tumor and 100 prostate-cancer patients treated using DirectSPR benefitted from sparing on average 2.6 mm and 4.4 mm of healthy tissue in beam direction, respectively. In the representative patient case, dose reduction in organs at risk close to the target volume was achieved, with a mean dose reduction of up to 16% in the brainstem. Patient-specific DECT-based treatment planning with reduced safety margins was successfully introduced into clinical routine. CONCLUSIONS A substantial increase in range prediction accuracy in clinical proton treatment planning was achieved by patient-specific DECT-based SPR prediction. For the first time, a relevant imaging-based reduction of range prediction uncertainty on a 2% level has been achieved.
Collapse
|
25
|
Kruis MF. Improving radiation physics, tumor visualisation, and treatment quantification in radiotherapy with spectral or dual-energy CT. J Appl Clin Med Phys 2021; 23:e13468. [PMID: 34743405 PMCID: PMC8803285 DOI: 10.1002/acm2.13468] [Citation(s) in RCA: 19] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/22/2022] [Revised: 10/13/2021] [Accepted: 10/19/2021] [Indexed: 12/11/2022] Open
Abstract
Over the past decade, spectral or dual‐energy CT has gained relevancy, especially in oncological radiology. Nonetheless, its use in the radiotherapy (RT) clinic remains limited. This review article aims to give an overview of the current state of spectral CT and to explore opportunities for applications in RT. In this article, three groups of benefits of spectral CT over conventional CT in RT are recognized. Firstly, spectral CT provides more information of physical properties of the body, which can improve dose calculation. Furthermore, it improves the visibility of tumors, for a wide variety of malignancies as well as organs‐at‐risk OARs, which could reduce treatment uncertainty. And finally, spectral CT provides quantitative physiological information, which can be used to personalize and quantify treatment.
Collapse
|
26
|
kV-kV and kV-MV DECT based estimation of proton stopping power ratio - a simulation study. Phys Med 2021; 89:182-192. [PMID: 34390901 DOI: 10.1016/j.ejmp.2021.07.038] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 04/22/2021] [Revised: 07/30/2021] [Accepted: 07/31/2021] [Indexed: 12/23/2022] Open
Abstract
PURPOSE This study aims to estimate the proton stopping power ratio (SPR) by using 80-120 kV and 120 kV-6 MV dual-energy CT (DECT) in a fully simulation-based approach for proton therapy dose calculations. METHODS Based on a virtual CT system, a two-step approach is applied to obtain the reference attenuation coefficient for image reconstruction. The effective atomic number (EAN) and electron density ratio (EDR) are estimated from two CT scans. The SPR is estimated using a calibration based on known materials to obtain a piecewise linear relationship between the EAN and the logarithm of the mean excitation energy, lnIm. The calibration phantoms are constructed from reference tissue materials taken from ICRU Report 44. Our approach is evaluated through using the ICRP110 human phantom. The respective influences of noise and beam hardening effects are studied. RESULTS With the beam hardening correction applied, the results of 120 kV-6 MV DECT are comparable to those of 80-120 kV DECT in predicting the EAN, but the former demonstrated a clear improvement in predicting the EDR and SPR. The 120 kV-6 MV DECT is able to predict the SPR within an accuracy of 10% for lung tissue and 5% for pelvis tissue, thereby outperforming the 80-120 kV DECT. CONCLUSIONS The 120 kV-6 MV DECT is less sensitive to noise but more susceptible to beam hardening effects. By applying beam hardening correction, the 120 kV-6 MV DECT can predict the SPR more accurately than the 80-120 kV DECT. To utilize our DECT approach most effectively, high-quality reconstructed images are required.
Collapse
|
27
|
Peters N, Wohlfahrt P, Dahlgren CV, de Marzi L, Ellerbrock M, Fracchiolla F, Free J, Gomà C, Góra J, Jensen MF, Kajdrowicz T, Mackay R, Molinelli S, Rinaldi I, Rompokos V, Siewert D, van der Tol P, Vermeren X, Nyström H, Lomax A, Richter C. Experimental assessment of inter-centre variation in stopping-power and range prediction in particle therapy. Radiother Oncol 2021; 163:7-13. [PMID: 34329653 DOI: 10.1016/j.radonc.2021.07.019] [Citation(s) in RCA: 17] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/15/2021] [Revised: 07/01/2021] [Accepted: 07/19/2021] [Indexed: 12/13/2022]
Abstract
PURPOSE Experimental assessment of inter-centre variation and absolute accuracy of stopping-power-ratio (SPR) prediction within 17 particle therapy centres of the European Particle Therapy Network. MATERIAL AND METHODS A head and body phantom with seventeen tissue-equivalent materials were scanned consecutively at the participating centres using their individual clinical CT scan protocol and translated into SPR with their in-house CT-number-to-SPR conversion. Inter-centre variation and absolute accuracy in SPR prediction were quantified for three tissue groups: lung, soft tissues and bones. The integral effect on range prediction for typical clinical beams traversing different tissues was determined for representative beam paths for the treatment of primary brain tumours as well as lung and prostate cancer. RESULTS An inter-centre variation in SPR prediction (2σ) of 8.7%, 6.3% and 1.5% relative to water was determined for bone, lung and soft-tissue surrogates in the head setup, respectively. Slightly smaller variations were observed in the body phantom (6.2%, 3.1%, 1.3%). This translated into inter-centre variation of integral range prediction (2σ) of 2.9%, 2.6% and 1.3% for typical beam paths of prostate-, lung- and primary brain-tumour treatments, respectively. The absolute error in range exceeded 2% in every fourth participating centre. The consideration of beam hardening and the execution of an independent HLUT validation had a positive effect, on average. CONCLUSION The large inter-centre variations in SPR and range prediction justify the currently clinically used margins accounting for range uncertainty, which are of the same magnitude as the inter-centre variation. This study underlines the necessity of higher standardisation in CT-number-to-SPR conversion.
Collapse
Affiliation(s)
- Nils Peters
- OncoRay - National Center for Radiation Research in Oncology, Faculty of Medicine and University Hospital Carl Gustav Carus, Technische Universität Dresden, Helmholtz-Zentrum Dresden - Rossendorf, Dresden, Germany; Helmholtz-Zentrum Dresden - Rossendorf, Institute of Radiooncology - OncoRay, Dresden, Germany.
| | - Patrick Wohlfahrt
- OncoRay - National Center for Radiation Research in Oncology, Faculty of Medicine and University Hospital Carl Gustav Carus, Technische Universität Dresden, Helmholtz-Zentrum Dresden - Rossendorf, Dresden, Germany; Helmholtz-Zentrum Dresden - Rossendorf, Institute of Radiooncology - OncoRay, Dresden, Germany
| | | | - Ludovic de Marzi
- Institut Curie, PSL Research University, University Paris Saclay, LITO, Orsay, France; Institut Curie, PSL Research University, Radiation Oncology Department, Proton Therapy Centre, Centre Universitaire, Orsay, France
| | - Malte Ellerbrock
- Heidelberg Ion-Beam Therapy Center (HIT), Department of Radiation Oncology, Heidelberg University Hospital, Heidelberg, Germany
| | | | - Jeffrey Free
- University of Groningen, University Medical Center Groningen, Department of Radiation Oncology, Groningen, The Netherlands
| | - Carles Gomà
- KU Leuven, Department of Oncology, Laboratory of Experimental Radiotherapy, Leuven, Belgium
| | - Joanna Góra
- EBG Medaustron GmbH, Wiener Neustadt, Austria
| | | | - Tomasz Kajdrowicz
- Institute of Nuclear Physics - Polish Academy of Sciences, Krakow, Poland
| | - Ranald Mackay
- University of Manchester - Faculty of Life Sciences, Manchester, United Kingdom
| | | | - Ilaria Rinaldi
- Department of Radiation Oncology (Maastro), GROW School for Oncology, Maastricht University Medical Centre, Maastricht, The Netherlands
| | | | | | | | - Xavier Vermeren
- Westdeutsches Protonentherapiezentrum Essen, Universitätsklinkum Essen, Germany
| | | | | | - Christian Richter
- OncoRay - National Center for Radiation Research in Oncology, Faculty of Medicine and University Hospital Carl Gustav Carus, Technische Universität Dresden, Helmholtz-Zentrum Dresden - Rossendorf, Dresden, Germany; Helmholtz-Zentrum Dresden - Rossendorf, Institute of Radiooncology - OncoRay, Dresden, Germany; Department of Radiotherapy and Radiation Oncology, Faculty of Medicine and University Hospital Carl Gustav Carus, Technische Universität Dresden, Germany; German Cancer Consortium (DKTK), partner site Dresden, and German Cancer Research Center (DKFZ), Heidelberg, Germany
| |
Collapse
|
28
|
Berthold J, Khamfongkhruea C, Petzoldt J, Thiele J, Hölscher T, Wohlfahrt P, Peters N, Jost A, Hofmann C, Janssens G, Smeets J, Richter C. First-In-Human Validation of CT-Based Proton Range Prediction Using Prompt Gamma Imaging in Prostate Cancer Treatments. Int J Radiat Oncol Biol Phys 2021; 111:1033-1043. [PMID: 34229052 DOI: 10.1016/j.ijrobp.2021.06.036] [Citation(s) in RCA: 21] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/18/2021] [Revised: 05/20/2021] [Accepted: 06/23/2021] [Indexed: 12/11/2022]
Abstract
PURPOSE Uncertainty in computed tomography (CT)-based range prediction substantially impairs the accuracy of proton therapy. Direct determination of the stopping-power ratio (SPR) from dual-energy CT (DECT) has been proposed (DirectSPR), and initial validation studies in phantoms and biological tissues have proven a high accuracy. However, a thorough validation of range prediction in patients has not yet been achieved by any means. Here, we present the first systematic validation of CT-based proton range prediction in patients using prompt gamma imaging (PGI). METHODS AND MATERIALS A PGI slit camera system with improved positioning accuracy, using a floor-based docking station, was used. Its overall uncertainty for range prediction validation was determined experimentally with both x-ray and beam measurements. The accuracy of range prediction in patients was determined from clinical PGI measurements during hypofractionated treatment of 5 patients with prostate cancer - in total 30 fractions with in-room control-CTs. For each pencil-beam-scanning spot, the range shift was obtained by comparing the PGI measurement to a control-CT-based PGI simulation. Three different SPR prediction approaches were applied in simulations: a standard CT-number-to-SPR conversion (Hounsfield look-up table [HLUT]), an adapted HLUT (DECT optimized), and DirectSPR. The spot-wise weighted mean range shift from all spots served as a measure for the accuracy of the respective range prediction approach. RESULTS A mean range prediction accuracy of 0.0% ± 0.5%, 0.3% ± 0.4%, and 1.8% ± 0.4% was obtained for DirectSPR, adapted HLUT, and standard HLUT, respectively. The overall validation uncertainty of the second-generation PGI slit camera is about 1 mm (2σ) for all approaches, which is smaller than the range prediction uncertainty for deep-seated tumors. CONCLUSIONS For the first time, range prediction accuracy was assessed in clinical routine using PGI range verification in prostate cancer treatments. Both DECT-derived range prediction approaches agree well with the measured proton range from PGI verification, whereas the standard HLUT approach differs relevantly. These results endorse the recent reduction of clinical safety margins in DirectSPR-based treatment planning in our institution.
Collapse
Affiliation(s)
- Jonathan Berthold
- OncoRay - National Center for Radiation Research in Oncology, Faculty of Medicine and University Hospital Carl Gustav Carus, Technische Universität Dresden, Helmholtz-Zentrum Dresden - Rossendorf, Dresden, Germany; Helmholtz-Zentrum Dresden - Rossendorf, Institute of Radiooncology - OncoRay, Dresden, Germany.
| | - Chirasak Khamfongkhruea
- OncoRay - National Center for Radiation Research in Oncology, Faculty of Medicine and University Hospital Carl Gustav Carus, Technische Universität Dresden, Helmholtz-Zentrum Dresden - Rossendorf, Dresden, Germany; Faculty of Medicine and Public Health, HRH Princess Chulabhorn College of Medical Science, Chulabhorn Royal Academy, Bangkok, Thailand
| | | | - Julia Thiele
- Department of Radiotherapy and Radiation Oncology, Faculty of Medicine and University Hospital Carl Gustav Carus, Technische Universität Dresden, Dresden, Germany
| | - Tobias Hölscher
- Department of Radiotherapy and Radiation Oncology, Faculty of Medicine and University Hospital Carl Gustav Carus, Technische Universität Dresden, Dresden, Germany
| | - Patrick Wohlfahrt
- OncoRay - National Center for Radiation Research in Oncology, Faculty of Medicine and University Hospital Carl Gustav Carus, Technische Universität Dresden, Helmholtz-Zentrum Dresden - Rossendorf, Dresden, Germany; Helmholtz-Zentrum Dresden - Rossendorf, Institute of Radiooncology - OncoRay, Dresden, Germany
| | - Nils Peters
- OncoRay - National Center for Radiation Research in Oncology, Faculty of Medicine and University Hospital Carl Gustav Carus, Technische Universität Dresden, Helmholtz-Zentrum Dresden - Rossendorf, Dresden, Germany; Helmholtz-Zentrum Dresden - Rossendorf, Institute of Radiooncology - OncoRay, Dresden, Germany
| | - Angelina Jost
- OncoRay - National Center for Radiation Research in Oncology, Faculty of Medicine and University Hospital Carl Gustav Carus, Technische Universität Dresden, Helmholtz-Zentrum Dresden - Rossendorf, Dresden, Germany; Beuth Hochschule für Technik, Berlin, Germany
| | | | | | | | - Christian Richter
- OncoRay - National Center for Radiation Research in Oncology, Faculty of Medicine and University Hospital Carl Gustav Carus, Technische Universität Dresden, Helmholtz-Zentrum Dresden - Rossendorf, Dresden, Germany; Helmholtz-Zentrum Dresden - Rossendorf, Institute of Radiooncology - OncoRay, Dresden, Germany; Department of Radiotherapy and Radiation Oncology, Faculty of Medicine and University Hospital Carl Gustav Carus, Technische Universität Dresden, Dresden, Germany; German Cancer Consortium (DKTK), partner site Dresden, Germany and German Cancer Research Center (DKFZ), Heidelberg, Germany
| |
Collapse
|
29
|
Dickmann J, Sarosiek C, Götz S, Pankuch M, Coutrakon G, Johnson RP, Schulte RW, Parodi K, Landry G, Dedes G. An empirical artifact correction for proton computed tomography. Phys Med 2021; 86:57-65. [PMID: 34058718 DOI: 10.1016/j.ejmp.2021.05.018] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 03/12/2021] [Revised: 04/20/2021] [Accepted: 05/12/2021] [Indexed: 12/31/2022] Open
Abstract
PURPOSE To reduce image artifacts of proton computed tomography (pCT) from a preclinical scanner, for imaging of the relative stopping power (RSP) needed for particle therapy treatment planning using a simple empirical artifact correction method. METHODS We adapted and employed a correction method previously used for beam-hardening correction in x-ray CT which makes use of a single scan of a custom-built homogeneous phantom with known RSP. Exploiting the linearity of the filtered backprojection operation, a function was found which corrects water-equivalent path lengths (RSP line integrals) in experimental scans using a prototype pCT scanner. The correction function was applied to projection values of subsequent scans of a homogeneous water phantom, a sensitometric phantom with various inserts and an anthropomorphic head phantom. Data were acquired at two different incident proton energies to test the robustness of the method. RESULTS Inaccuracies in the detection process caused an offset and known ring artifacts in the water phantom which were considerably reduced using the proposed method. The mean absolute percentage error (MAPE) of mean RSP values of all inserts of the sensitometric phantom and the water phantom was reduced from 0.87% to 0.44% and from 0.86% to 0.48% for the two incident energies respectively. In the head phantom a clear reduction of artifacts was observed. CONCLUSIONS Image artifacts of experimental pCT scans with a prototype scanner could substantially be reduced both in homogeneous, heterogeneous and anthropomorphic phantoms. RSP accuracy was also improved.
Collapse
Affiliation(s)
- Jannis Dickmann
- Department of Medical Physics, Fakultät für Physik, Ludwig-Maximilians-Universität München (LMU Munich), Am Coulombwall 1, Garching bei München, Germany.
| | - Christina Sarosiek
- Department of Physics, Northern Illinois University, 1425 W. Lincoln Highway, DeKalb, Illinois, United States.
| | - Stefanie Götz
- Department of Medical Physics, Fakultät für Physik, Ludwig-Maximilians-Universität München (LMU Munich), Am Coulombwall 1, Garching bei München, Germany.
| | - Mark Pankuch
- Northwestern Medicine Chicago Proton Center, 4455 Weaver Parkway, Warrenville, Illinois, United States.
| | - George Coutrakon
- Department of Physics, Northern Illinois University, 1425 W. Lincoln Highway, DeKalb, Illinois, United States.
| | - Robert P Johnson
- Department of Physics, U.C. Santa Cruz, 1156 High Street, Santa Cruz, California, United States.
| | - Reinhard W Schulte
- Division of Biomedical Engineering Sciences, Loma Linda University, 11175 Campus Street, Loma Linda, California, United States.
| | - Katia Parodi
- Department of Medical Physics, Fakultät für Physik, Ludwig-Maximilians-Universität München (LMU Munich), Am Coulombwall 1, Garching bei München, Germany.
| | - Guillaume Landry
- Department of Medical Physics, Fakultät für Physik, Ludwig-Maximilians-Universität München (LMU Munich), Am Coulombwall 1, Garching bei München, Germany; Department of Radiation Oncology, University Hospital, LMU Munich, Marchioninistraße 15, Munich, Germany; German Cancer Consortium (DKTK), Marchioninistraße 15, Munich, Germany.
| | - George Dedes
- Department of Medical Physics, Fakultät für Physik, Ludwig-Maximilians-Universität München (LMU Munich), Am Coulombwall 1, Garching bei München, Germany.
| |
Collapse
|
30
|
Paganetti H, Beltran C, Both S, Dong L, Flanz J, Furutani K, Grassberger C, Grosshans DR, Knopf AC, Langendijk JA, Nystrom H, Parodi K, Raaymakers BW, Richter C, Sawakuchi GO, Schippers M, Shaitelman SF, Teo BKK, Unkelbach J, Wohlfahrt P, Lomax T. Roadmap: proton therapy physics and biology. Phys Med Biol 2021; 66. [DOI: 10.1088/1361-6560/abcd16] [Citation(s) in RCA: 30] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/21/2020] [Accepted: 11/23/2020] [Indexed: 12/12/2022]
|
31
|
Hu Y, Ding X, Shen J, Bues M, Liu W, Kang Y, Leng S, Yu L. Feasibility of using megavoltage computed tomography to reduce proton range uncertainty: A simulation study. J Appl Clin Med Phys 2021; 22:131-140. [PMID: 33608992 PMCID: PMC7984486 DOI: 10.1002/acm2.13191] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/03/2019] [Revised: 01/10/2021] [Accepted: 01/11/2021] [Indexed: 12/25/2022] Open
Abstract
Purpose To demonstrate that variation in chemical composition has a negligible effect on the mapping curve from relative electron density (RED) to proton stopping power ratio (SPR), and to establish the theoretical framework of using Megavoltage (MV) computed tomography (CT), instead of kilovoltage (kV) dual energy CT, to accurately estimate proton SPR. Methods A simulation study was performed to evaluate the effect of chemical composition variation on kVCT number and proton SPR. The simulation study involved both reference and simulated human tissues. The reference human tissues, together with their physical densities and chemical compositions, came from the ICRP publication 23. The simulated human tissues were created from the reference human tissues assuming that elemental percentage weight followed a Gaussian distribution. For all tissues, kVCT number and proton SPR were obtained through (a) theoretical calculation from tissue’s physical density and chemical composition which served as the ground truth, and (b) estimation from RED using the calibration curves established from the stoichiometric method. Deviations of the estimated values from the calculated values were quantified as errors in using RED to estimate kVCT number and proton SPR. Results Given a chemical composition variation of 5% (1σ) of the nominal percentage weights, the total estimation error of using RED to estimate kVCT number was 0.34%, 0.62%, and 0.77% and the total estimation error of using RED to estimate proton SPR was 0.30%, 0.22%, and 0.16% for fat tissues, non‐fat soft tissues and bone tissues, respectively. Conclusion Chemical composition had a negligible effect on the method of using RED to determine proton SPR. RED itself is sufficient to accurately determine proton SPR. MVCT number maintains a superb linear relationship with RED because it is highly dominated by Compton scattering. Therefore, MVCT has great potential in reducing the proton range uncertainty.
Collapse
Affiliation(s)
- Yanle Hu
- Department of Radiation Oncology, Mayo Clinic, Phoenix, AZ, USA
| | - Xiaoning Ding
- Department of Radiation Oncology, Mayo Clinic, Phoenix, AZ, USA
| | - Jiajian Shen
- Department of Radiation Oncology, Mayo Clinic, Phoenix, AZ, USA
| | - Martin Bues
- Department of Radiation Oncology, Mayo Clinic, Phoenix, AZ, USA
| | - Wei Liu
- Department of Radiation Oncology, Mayo Clinic, Phoenix, AZ, USA
| | - Yixiu Kang
- Department of Radiation Oncology, Mayo Clinic, Phoenix, AZ, USA
| | - Shuai Leng
- Department of Radiology, Mayo Clinic, Rochester, MN, USA
| | - Lifeng Yu
- Department of Radiology, Mayo Clinic, Rochester, MN, USA
| |
Collapse
|
32
|
Pettersen HES, Volz L, Sølie JR, Alme J, Barnaföldi GG, Barthel R, van den Brink A, Borshchov V, Chaar M, Eikeland V, Genov G, Grøttvik O, Helstrup H, Keidel R, Kobdaj C, van der Kolk N, Mehendale S, Meric I, Harald Odland O, Papp G, Peitzmann T, Piersimoni P, Protsenko M, Ur Rehman A, Richter M, Tefre Samnøy A, Seco J, Shafiee H, Songmoolnak A, Tambave G, Tymchuk I, Ullaland K, Varga-Kofarago M, Wagner B, Xiao R, Yang S, Yokoyama H, Röhrich D. Helium radiography with a digital tracking calorimeter-a Monte Carlo study for secondary track rejection. Phys Med Biol 2021; 66:035004. [PMID: 33181502 DOI: 10.1088/1361-6560/abca03] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
Radiation therapy using protons and heavier ions is a fast-growing therapeutic option for cancer patients. A clinical system for particle imaging in particle therapy would enable online patient position verification, estimation of the dose deposition through range monitoring and a reduction of uncertainties in the calculation of the relative stopping power of the patient. Several prototype imaging modalities offer radiography and computed tomography using protons and heavy ions. A Digital Tracking Calorimeter (DTC), currently under development, has been proposed as one such detector. In the DTC 43 longitudinal layers of laterally stacked ALPIDE CMOS monolithic active pixel sensor chips are able to reconstruct a large number of simultaneously recorded proton tracks. In this study, we explored the capability of the DTC for helium imaging which offers favorable spatial resolution over proton imaging. Helium ions exhibit a larger cross section for inelastic nuclear interactions, increasing the number of produced secondaries in the imaged object and in the detector itself. To that end, a filtering process able to remove a large fraction of the secondaries was identified, and the track reconstruction process was adapted for helium ions. By filtering on the energy loss along the tracks, on the incoming angle and on the particle ranges, 97.5% of the secondaries were removed. After passing through 16 cm water, 50.0% of the primary helium ions survived; after the proposed filtering 42.4% of the primaries remained; finally after subsequent image reconstruction 31% of the primaries remained. Helium track reconstruction leads to more track matching errors compared to protons due to the increased available focus strength of the helium beam. In a head phantom radiograph, the Water Equivalent Path Length error envelope was 1.0 mm for helium and 1.1 mm for protons. This accuracy is expected to be sufficient for helium imaging for pre-treatment verification purposes.
Collapse
|
33
|
Tanaka S, Noto Y, Utsunomiya S, Yoshimura T, Matsuura T, Saito M. Proton dose calculation based on converting dual-energy CT data to stopping power ratio (DEEDZ-SPR): a beam-hardening assessment. ACTA ACUST UNITED AC 2020; 65:235046. [DOI: 10.1088/1361-6560/abae09] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
|
34
|
Permatasari FF, Eulitz J, Richter C, Wohlfahrt P, Lühr A. Material assignment for proton range prediction in Monte Carlo patient simulations using stopping-power datasets. Phys Med Biol 2020; 65:185004. [DOI: 10.1088/1361-6560/ab9702] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
|
35
|
Parodi K. Latest developments in in-vivo imaging for proton therapy. Br J Radiol 2020; 93:20190787. [PMID: 31794249 PMCID: PMC7066959 DOI: 10.1259/bjr.20190787] [Citation(s) in RCA: 26] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/13/2019] [Revised: 11/18/2019] [Accepted: 11/28/2019] [Indexed: 12/11/2022] Open
Abstract
Owing to the favorable physical and biological properties of swift ions in matter, their application to radiation therapy for highly selective cancer treatment is rapidly spreading worldwide. To date, over 90 ion therapy facilities are operational, predominantly with proton beams, and about the same amount is under construction or planning.Over the last decades, considerable developments have been achieved in accelerator technology, beam delivery and medical physics to enhance conformation of the dose delivery to complex shaped tumor volumes, with excellent sparing of surrounding normal tissue and critical organs. Nevertheless, full clinical exploitation of the ion beam advantages is still challenged, especially by uncertainties in the knowledge of the beam range in the actual patient anatomy during the fractionated course of treatment, thus calling for continued multidisciplinary research in this rapidly emerging field.This contribution will review latest developments aiming to image the patient with the same beam quality as for therapy prior to treatment, and to visualize in-vivo the treatment delivery by exploiting irradiation-induced physical emissions, with different level of maturity from proof-of-concept studies in phantoms and first in-silico studies up to clinical testing and initial clinical evaluation.
Collapse
Affiliation(s)
- Katia Parodi
- Department of Experimental Physics – Medical Physics, Ludwig-Maximilians-Universität München, Faculty of Physics, Munich, Germany
| |
Collapse
|
36
|
Held KD, Lomax AJ, Troost EGC. Proton therapy special feature: introductory editorial. Br J Radiol 2020; 93:20209004. [PMID: 32081045 DOI: 10.1259/bjr.20209004] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022] Open
Affiliation(s)
- Kathryn D Held
- Department of Radiation Oncology, Massachusetts General Hospital/Harvard Medical School, Boston, MA, USA
| | - Antony J Lomax
- Center for Proton Therapy, Paul Scherrer Institute, Villigen, Switzerland.,Department of Physics, ETH Zürich, Zürich, Switzerland
| | - Esther G C Troost
- Department of Radiotherapy and Radiation Oncology, Faculty of Medicine and University Hospital Carl Gustav Carus, Technische Universität Dresden, Dresden, Germany.,OncoRay - National Center for Radiation Research in Oncology, Dresden, Germany
| |
Collapse
|
37
|
Thorwarth D. Imaging science and development in modern high-precision radiotherapy. Phys Imaging Radiat Oncol 2019; 12:63-66. [PMID: 33458297 PMCID: PMC7807660 DOI: 10.1016/j.phro.2019.11.008] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023] Open
Affiliation(s)
- Daniela Thorwarth
- Section for Biomedical Physics, Department of Radiation Oncology, University of Tübingen, Germany
| |
Collapse
|