1
|
Kraan AC, Moglioni M, Battistoni G, Bersani D, Berti A, Carra P, Cerello P, Ciocca M, Ferrero V, Fiorina E, Mazzoni E, Morrocchi M, Muraro S, Orlandi E, Pennazio F, Retico A, Rosso V, Sportelli G, Vischioni B, Vitolo V, Bisogni MG. Using the gamma-index analysis for inter-fractional comparison of in-beam PET images for head-and-neck treatment monitoring in proton therapy: A Monte Carlo simulation study. Phys Med 2024; 120:103329. [PMID: 38492331 DOI: 10.1016/j.ejmp.2024.103329] [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: 05/15/2023] [Revised: 02/13/2024] [Accepted: 03/07/2024] [Indexed: 03/18/2024] Open
Abstract
GOAL In-beam Positron Emission Tomography (PET) is a technique for in-vivo non-invasive treatment monitoring for proton therapy. To detect anatomical changes in patients with PET, various analysis methods exist, but their clinical interpretation is problematic. The goal of this work is to investigate whether the gamma-index analysis, widely used for dose comparisons, is an appropriate tool for comparing in-beam PET distributions. Focusing on a head-and-neck patient, we investigate whether the gamma-index map and the passing rate are sensitive to progressive anatomical changes. METHODS/MATERIALS We simulated a treatment course of a proton therapy patient using FLUKA Monte Carlo simulations. Gradual emptying of the sinonasal cavity was modeled through a series of artificially modified CT scans. The in-beam PET activity distributions from three fields were evaluated, simulating a planar dual head geometry. We applied the 3D-gamma evaluation method to compare the PET images with a reference image without changes. Various tolerance criteria and parameters were tested, and results were compared to the CT-scans. RESULTS Based on 210 MC simulations we identified appropriate parameters for the gamma-index analysis. Tolerance values of 3 mm/3% and 2 mm/2% were suited for comparison of simulated in-beam PET distributions. The gamma passing rate decreased with increasing volume change for all fields. CONCLUSION The gamma-index analysis was found to be a useful tool for comparing simulated in-beam PET images, sensitive to sinonasal cavity emptying. Monitoring the gamma passing rate behavior over the treatment course is useful to detect anatomical changes occurring during the treatment course.
Collapse
Affiliation(s)
- Aafke Christine Kraan
- Istituto Nazionale di Fisica Nucleare, Sezione di Pisa, Largo Bruno Pontecorvo 3, Pisa, 56127, Italy
| | - Martina Moglioni
- Istituto Nazionale di Fisica Nucleare, Sezione di Pisa, Largo Bruno Pontecorvo 3, Pisa, 56127, Italy; Dipartimento di Fisica, Università di Pisa, Largo Bruno Pontecorvo 3, Pisa, 56127, Italy.
| | - Giuseppe Battistoni
- Istituto Nazionale di Fisica Nucleare, Sezione di Milano, Via Giovanni Celoria 16, Milano, 20133, Italy
| | - Davide Bersani
- Istituto Nazionale di Fisica Nucleare, Sezione di Torino, Via Pietro Giuria 1, Torino, 10125, Italy
| | - Andrea Berti
- Istituto Nazionale di Fisica Nucleare, Sezione di Pisa, Largo Bruno Pontecorvo 3, Pisa, 56127, Italy; Dipartimento di Fisica, Università di Pisa, Largo Bruno Pontecorvo 3, Pisa, 56127, Italy
| | - Pietro Carra
- Istituto Nazionale di Fisica Nucleare, Sezione di Pisa, Largo Bruno Pontecorvo 3, Pisa, 56127, Italy; Dipartimento di Fisica, Università di Pisa, Largo Bruno Pontecorvo 3, Pisa, 56127, Italy
| | - Piergiorgio Cerello
- Istituto Nazionale di Fisica Nucleare, Sezione di Torino, Via Pietro Giuria 1, Torino, 10125, Italy
| | - Mario Ciocca
- Centro Nazionale di Adroterapia Oncologica, Strada Privata Campeggi 53, Pavia, 27100, Italy
| | - Veronica Ferrero
- Istituto Nazionale di Fisica Nucleare, Sezione di Torino, Via Pietro Giuria 1, Torino, 10125, Italy
| | - Elisa Fiorina
- Istituto Nazionale di Fisica Nucleare, Sezione di Torino, Via Pietro Giuria 1, Torino, 10125, Italy
| | - Enrico Mazzoni
- Istituto Nazionale di Fisica Nucleare, Sezione di Pisa, Largo Bruno Pontecorvo 3, Pisa, 56127, Italy
| | - Matteo Morrocchi
- Istituto Nazionale di Fisica Nucleare, Sezione di Pisa, Largo Bruno Pontecorvo 3, Pisa, 56127, Italy; Dipartimento di Fisica, Università di Pisa, Largo Bruno Pontecorvo 3, Pisa, 56127, Italy
| | - Silvia Muraro
- Istituto Nazionale di Fisica Nucleare, Sezione di Milano, Via Giovanni Celoria 16, Milano, 20133, Italy
| | - Ester Orlandi
- Centro Nazionale di Adroterapia Oncologica, Strada Privata Campeggi 53, Pavia, 27100, Italy
| | - Francesco Pennazio
- Istituto Nazionale di Fisica Nucleare, Sezione di Torino, Via Pietro Giuria 1, Torino, 10125, Italy
| | - Alessandra Retico
- Istituto Nazionale di Fisica Nucleare, Sezione di Pisa, Largo Bruno Pontecorvo 3, Pisa, 56127, Italy
| | - Valeria Rosso
- Istituto Nazionale di Fisica Nucleare, Sezione di Pisa, Largo Bruno Pontecorvo 3, Pisa, 56127, Italy; Dipartimento di Fisica, Università di Pisa, Largo Bruno Pontecorvo 3, Pisa, 56127, Italy
| | - Giancarlo Sportelli
- Istituto Nazionale di Fisica Nucleare, Sezione di Pisa, Largo Bruno Pontecorvo 3, Pisa, 56127, Italy; Dipartimento di Fisica, Università di Pisa, Largo Bruno Pontecorvo 3, Pisa, 56127, Italy
| | - Barbara Vischioni
- Centro Nazionale di Adroterapia Oncologica, Strada Privata Campeggi 53, Pavia, 27100, Italy
| | - Viviana Vitolo
- Centro Nazionale di Adroterapia Oncologica, Strada Privata Campeggi 53, Pavia, 27100, Italy
| | - Maria Giuseppina Bisogni
- Istituto Nazionale di Fisica Nucleare, Sezione di Pisa, Largo Bruno Pontecorvo 3, Pisa, 56127, Italy; Dipartimento di Fisica, Università di Pisa, Largo Bruno Pontecorvo 3, Pisa, 56127, Italy
| |
Collapse
|
2
|
Yu L, Baker A, Kairn T, Livingstone A, Trapp J, Crowe SB. A structure-based gamma evaluation method for identifying clinically relevant dose differences in organs at risk. Phys Eng Sci Med 2023; 46:1033-1041. [PMID: 37219798 PMCID: PMC10480250 DOI: 10.1007/s13246-023-01270-3] [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: 08/02/2022] [Accepted: 05/01/2023] [Indexed: 05/24/2023]
Abstract
Gamma evaluation is currently the most widely used dose comparison method for patient specific quality assurance (PSQA). However, existing methods for normalising the dose difference, using either the dose at the global maximum dose point or at each local point, can respectively lead to under- and over-sensitivity to dose differences in organ-at-risk structures. This may be of concern for plan evaluation from clinical perspectives. This study has explored and proposed a new method called structural gamma, which takes structural dose tolerances into consideration while performing gamma analysis for PSQA. As a demonstration of the structural gamma method, a total of 78 retrospective plans on four treatment sites were re-calculated on an in-house Monte Carlo system and compared with doses calculated from the treatment planning system. Structural gamma evaluations were performed using both QUANTEC dose tolerances and radiation oncologist specified dose tolerances, then compared with conventional global and local gamma evaluations. Results demonstrated that structural gamma evaluation is especially sensitive to errors in structures with restrictive dose constraints. The structural gamma map provides both geometric and dosimetric information on PSQA results, allowing straightforward clinical interpretation. The proposed structure-based gamma method accounts for dose tolerances for specific anatomical structures. This method can provide a clinically useful method to assess and communicate PSQA results, offering radiation oncologists a more intuitive way of examining agreement in surrounding critical normal structures.
Collapse
Affiliation(s)
- Liting Yu
- Royal Brisbane and Women's Hospital, Herston, 4029, Australia.
- Queensland University of Technology, Brisbane, 4001, Australia.
| | - Anthony Baker
- Royal Brisbane and Women's Hospital, Herston, 4029, Australia
- Nepean Hospital, Kingswood, 2747, Australia
| | - Tanya Kairn
- Royal Brisbane and Women's Hospital, Herston, 4029, Australia
- Queensland University of Technology, Brisbane, 4001, Australia
| | | | - Jamie Trapp
- Queensland University of Technology, Brisbane, 4001, Australia
| | - Scott B Crowe
- Royal Brisbane and Women's Hospital, Herston, 4029, Australia
- Queensland University of Technology, Brisbane, 4001, Australia
| |
Collapse
|
3
|
Lu W, Li Y, Huang W, Cui H, Zhang H, Yi X. Optimizing the Region for Evaluation of Global Gamma Analysis for Nasopharyngeal Cancer (NPC) Pretreatment IMRT QA by COMPASS: A Retrospective Study. Front Oncol 2022; 12:859415. [PMID: 35774127 PMCID: PMC9238439 DOI: 10.3389/fonc.2022.859415] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/21/2022] [Accepted: 05/10/2022] [Indexed: 11/13/2022] Open
Abstract
Background The global gamma passing rate is the most commonly used metric for patient-specific pretreatment quality assurance in radiation therapy. However, the optimal region for evaluation and specific action limits (ALs) need to be explored. Therefore, this study was carried out to explore the optimal region for evaluation of the global gamma passing rate and define ALs by using the COMPASS software. Methods A total of 93 intensity-modulated radiation therapy (IMRT) plans for nasopharyngeal cancer (NPC) patients, including 61 original plans and 32 multileaf collimator (MLC) error-introduced test plans, were selected for retrospective analysis. Firstly, the dose distribution was divided into six isodose regions (“≥10%”, “≥20%”, “≥30%”, “≥40%”, “≥50%”, and “≥60%”) based on the prescribed dose and one clinically oriented region for evaluation (“whole”) to perform the three-dimensional (3D) global gamma reanalysis. Meanwhile, the percentage gamma passing rate (%GP), mean gamma index (μGI) based on 3%/2 mm criteria, and percentage dose error (%DE) of the dose–volume histogram (DVH) metrics were recorded by COMPASS application. Secondly, the Pearson’s correlation coefficient was used to analyze the correlation between %GP and %DE and between μGI and %DE in different regions. Additionally, receiver operating characteristic (ROC) methodology was applied to quantify the fraction of patient-specific plans evaluated as “fail” and “pass”. In order to improve the correlation between gamma analysis result and clinical criteria, ROC analysis was carried out in accordance with hybridization analysis criteria (%DE ≤3%, %GP ≥90% and %DE ≤3%, μGI ≤0.6). ROC was performed for two purposes: 1) to analyze the sensitivity and specificity of %GP and μGI in different regions for evaluation and 2) to define the ALs of %GP and μGI in the optimal region for evaluation. Finally, the plans introduced with MLC errors were prepared for validation. Moreover, we also compared the positive rate of ALs of both %GP and μGI in detecting MLC error-introduced plans in different regions. Results 1) In our study, a number of DVH-based metrics were found to be correlated with the evaluation parameters. The corresponding number was 4, 2, 1, 1, 1, 1, and 3 in γwhole, γ10%, γ20%, γ30%, γ40%, γ50%, and γ60%, respectively, and 5, 3, 0, 1, 1, 4, and 2 in μGIwhole, μGI10%, μGI20%, μGI30%, μGI40%, μGI50%, and μGI60%, respectively. The results by COMPASS have revealed that the %DE of specific structures has a slightly higher correlation with both %GP and μGI of the “whole” region compared with that of any other region. However, it is a moderate correlation (0.5 ≤ |r| < 0.8). 2) The areas under the curves (AUCs) of γwhole, μGIwhole, μGI40%, μGI50%, and μGI60% were >0.7 based on 3%/2 mm criteria. According to the Youden coefficient, we defined the ALs of γwhole ≥92%, μGIwhole ≤0.36, μGI40% ≤0.43, and μGI60% ≤0.40 based on 3%/2 mm criteria. 3) In the validation, for original plans, the accuracy of ALγwhole, ALγ10%, ALμGIwhole, ALμGI40%, ALμGI50%, and ALμGI60% was 23%, 9.8%, 90%, 80.3%, 9.8%, and 88.5%, respectively. For test plans with systematic MLC errors smaller than 0.8 mm, the positive rates of ALγwhole, ALγ10%, ALμGIwhole, ALμGI40%, ALμGI50%, and ALμGI60% were 25%, 58%, 92%, 92%, 42%, and 100%, respectively. For the plans with systematic MLC errors higher than 0.8 mm, the positive rates of all AL%GP and ALμGI were 100%. From the COMPASS validation results, the accuracy of γwhole, μGIwhole, μGI40%, and μGI60% was higher than that of the conventional γ10% and commonly used μGI50%. Conclusions Compared with the traditional evaluation region (i.e., the criteria with a threshold of 10% or a threshold of 50%, it was the same with the isodose regions of “≥10%”, “≥50%” based on the prescribed dose in our study), the “whole” region is more meaningful to the clinic by COMPASS. The accuracy of μGIwhole is higher than that of the conventional γ10% and the commonly used μGI50%.
Collapse
|