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.

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.

Synthetic CT imaging for PET monitoring in proton therapy: a simulation study.

Objective.This study addresses a fundamental limitation of in-beam positron emission tomography (IB-PET) in proton therapy: the lack of direct anatomical representation in the images it produces. We aim to overcome this shortcoming by pioneering the application of deep learning techniques to create synthetic control CT images (sCT) from combining IB-PET and planning CT scan data.Approach.We conducted simulations involving six patients who underwent irradiation with proton beams. Leveraging the architecture of a visual transformer (ViT) neural network, we developed a model to generate sCT images of these patients using the planning CT scans and the inter-fractional simulated PET activity maps during irradiation. To evaluate the model's performance, a comparison was conducted between the sCT images produced by the ViT model and the authentic control CT images-serving as the benchmark.Main results.The structural similarity index was computed at a mean value across all patients of 0.91, while the mean absolute error measured 22 Hounsfield Units (HU). Root mean squared error and peak signal-to-noise ratio values were 56 HU and 30 dB, respectively. The Dice similarity coefficient exhibited a value of 0.98. These values are comparable to or exceed those found in the literature. More than 70% of the synthetic morphological changes were found to be geometrically compatible with the ones reported in the real control CT scan.Significance.Our study presents an innovative approach to surface the hidden anatomical information of IB-PET in proton therapy. Our ViT-based model successfully generates sCT images from inter-fractional PET data and planning CT scans. Our model's performance stands on par with existing models relying on input from cone beam CT or magnetic resonance imaging, which contain more anatomical information than activity maps.

Few-seconds range verification with short-lived positron emitters in carbon ion therapy.

In-beam PET (Positron Emission Tomography) is one of the most precise techniques for in-vivo range monitoring in hadron therapy. Our objective was to demonstrate the feasibility of a short irradiation run for range verification before a carbon-ion treatment. To do so a PMMA target was irradiated with a 220 MeV/u carbon-ion beam and annihilation coincidences from short-lived positron emitters were acquired after irradiations lasting 0.6 s. The experiments were performed at the synchrotron-based facility CNAO (Italian National Center of Oncological Hadrontherapy) by using the INSIDE in-beam PET detector. The results show that, with 3·107 carbon ions, the reconstructed positron emitting nuclei distribution is in good agreement with the predictions of a detailed FLUKA Monte Carlo study. Moreover, the radio-nuclei production is sufficiently abundant to determine the average ion beam range with a σ of 1 mm with a 6 s measurement of the activity distribution. Since the data were acquired when the beam was off, the proposed rapid calibration method can be applied to hadron beams extracted from accelerators with very different time structures.

Design and performance evaluation of a slit-slat camera for 2D prompt gamma imaging in proton therapy monitoring: A Monte Carlo simulation study.

PURPOSE: We investigated the design of a prompt gamma camera for real-time dose delivery verification and the partial mitigation of range uncertainties. METHODS: A slit slat (SS) camera was optimized using the trade-off between the signal-to-noise ratio and spatial resolution. Then, using the GATE Monte Carlo package, the camera performances were estimated by means of target shifts, beam position quantification, changing the camera distance from the beam, and air cavity inserting. A homogeneous PMMA phantom and the air gaps induced PMMA phantom were used. The air gaps ranged from 5 mm to 30 mm by 5 mm increments were positioned in the middle of the beam range. To reduce the simulation time, phase space scoring was used. The batch method with five realizations was used for stochastic error calculations. RESULTS: The system's detection efficiency was 1.1 × 10 - 4 PGs Emitted PGs ( 1.8 × 10 - 5 $1.1 \times {10}^{-4}\frac{{\rm PGs}}{{\rm Emitted}\ {\rm PGs}}\ (1.8 \times {10}^{-5}$ PGs/proton) for a 10 × 20 cm2 detector (source-to-collimator distance = 15.0 cm). Axial and transaxial resolutions were 23 mm and 18 mm, respectively. The SS camera estimated the range as 69.0 ± 3.4 (relative stochastic error 1-sigma is 5%) and 67.6 ± 1.8 mm (2.6%) for the real range of 67.0 mm for 107 and 108 protons of 100 MeV, respectively. Considering 160 MeV, these values are 155.5 ± 3.1 (2%) and 152.2 ± 2.0 mm (1.3%) for the real range of 152.0 mm for 107 and 108 protons, respectively. Considering phantom shift, for a 100 MeV beam, the precision of the quantification (1-sigma) in the axial and lateral phantom shift estimation is 2.6 mm and 1 mm, respectively. Accordingly, the axial and lateral quantification precisions were 1.3 mm and 1 mm for a 160 MeV beam, respectively. Furthermore, the quantification of an air gap formulated as gap d e t = 0.98 × gap real ${{\rm gap}}_{det}=0.98 \times {{\rm gap}}_{{\rm real}}$ , where gap d e t ${{\rm gap}}_{det}$ and gapreal are the estimated and real air gap, respectively. The precision of the air gap quantification is 1.6 mm (1 sigma). Moreover, 2D PG images show the trajectory of the proton beam through the phantom. CONCLUSION: The proposed slit-slat imaging systems can potentially provide a real-time, in-vivo, and non-invasive treatment monitoring method for proton therapy.

In-vivo range verification analysis with in-beam PET data for patients treated with proton therapy at CNAO.

Morphological changes that may arise through a treatment course are probably one of the most significant sources of range uncertainty in proton therapy. Non-invasive in-vivo treatment monitoring is useful to increase treatment quality. The INSIDE in-beam Positron Emission Tomography (PET) scanner performs in-vivo range monitoring in proton and carbon therapy treatments at the National Center of Oncological Hadrontherapy (CNAO). It is currently in a clinical trial (ID: NCT03662373) and has acquired in-beam PET data during the treatment of various patients. In this work we analyze the in-beam PET (IB-PET) data of eight patients treated with proton therapy at CNAO. The goal of the analysis is twofold. First, we assess the level of experimental fluctuations in inter-fractional range differences (sensitivity) of the INSIDE PET system by studying patients without morphological changes. Second, we use the obtained results to see whether we can observe anomalously large range variations in patients where morphological changes have occurred. The sensitivity of the INSIDE IB-PET scanner was quantified as the standard deviation of the range difference distributions observed for six patients that did not show morphological changes. Inter-fractional range variations with respect to a reference distribution were estimated using the Most-Likely-Shift (MLS) method. To establish the efficacy of this method, we made a comparison with the Beam's Eye View (BEV) method. For patients showing no morphological changes in the control CT the average range variation standard deviation was found to be 2.5 mm with the MLS method and 2.3 mm with the BEV method. On the other hand, for patients where some small anatomical changes occurred, we found larger standard deviation values. In these patients we evaluated where anomalous range differences were found and compared them with the CT. We found that the identified regions were mostly in agreement with the morphological changes seen in the CT scan.

Proton therapy monitoring: spatiotemporal emission reconstruction with prompt gamma timing and implementation with PET detectors.

Objective. In this study we introduce spatiotemporal emission reconstruction prompt gamma timing (SER-PGT), a new method to directly reconstruct the prompt photon emission in the space and time domains inside the patient in proton therapy.Approach. SER-PGT is based on the numerical optimisation of a multidimensional likelihood function, followed by a post-processing of the results. The current approach relies on a specific implementation of the maximum-likelihood expectation maximisation algorithm. The robustness of the method is guaranteed by the complete absence of any information about the target composition in the algorithm.Main results. Accurate Monte Carlo simulations indicate a range resolution of about 0.5 cm (standard deviation) when considering 107primary protons impinging on an homogeneous phantom. Preliminary results on an anthropomorphic phantom are also reported.Significance. By showing the feasibility for the reconstruction of the primary particle range using PET detectors, this study provides significant basis for the development of an hybrid in-beam PET and prompt photon device.

Localization of anatomical changes in patients during proton therapy with in-beam PET monitoring: A voxel-based morphometry approach exploiting Monte Carlo simulations.

PURPOSE: In-beam positron emission tomography (PET) is one of the modalities that can be used for in vivo noninvasive treatment monitoring in proton therapy. Although PET monitoring has been frequently applied for this purpose, there is still no straightforward method to translate the information obtained from the PET images into easy-to-interpret information for clinical personnel. The purpose of this work is to propose a statistical method for analyzing in-beam PET monitoring images that can be used to locate, quantify, and visualize regions with possible morphological changes occurring over the course of treatment. METHODS: We selected a patient treated for squamous cell carcinoma (SCC) with proton therapy, to perform multiple Monte Carlo (MC) simulations of the expected PET signal at the start of treatment, and to study how the PET signal may change along the treatment course due to morphological changes. We performed voxel-wise two-tailed statistical tests of the simulated PET images, resembling the voxel-based morphometry (VBM) method commonly used in neuroimaging data analysis, to locate regions with significant morphological changes and to quantify the change. RESULTS: The VBM resembling method has been successfully applied to the simulated in-beam PET images, despite the fact that such images suffer from image artifacts and limited statistics. Three dimensional probability maps were obtained, that allowed to identify interfractional morphological changes and to visualize them superimposed on the computed tomography (CT) scan. In particular, the characteristic color patterns resulting from the two-tailed statistical tests lend themselves to trigger alarms in case of morphological changes along the course of treatment. CONCLUSIONS: The statistical method presented in this work is a promising method to apply to PET monitoring data to reveal interfractional morphological changes in patients, occurring over the course of treatment. Based on simulated in-beam PET treatment monitoring images, we showed that with our method it was possible to correctly identify the regions that changed. Moreover we could quantify the changes, and visualize them superimposed on the CT scan. The proposed method can possibly help clinical personnel in the replanning procedure in adaptive proton therapy treatments.

Monitoring Carbon Ion Beams Transverse Position Detecting Charged Secondary Fragments: Results From Patient Treatment Performed at CNAO.

Particle therapy in which deep seated tumours are treated using 12C ions (Carbon Ions RadioTherapy or CIRT) exploits the high conformity in the dose release, the high relative biological effectiveness and low oxygen enhancement ratio of such projectiles. The advantages of CIRT are driving a rapid increase in the number of centres that are trying to implement such technique. To fully profit from the ballistic precision achievable in delivering the dose to the target volume an online range verification system would be needed, but currently missing. The 12C ions beams range could only be monitored by looking at the secondary radiation emitted by the primary beam interaction with the patient tissues and no technical solution capable of the needed precision has been adopted in the clinical centres yet. The detection of charged secondary fragments, mainly protons, emitted by the patient is a promising approach, and is currently being explored in clinical trials at CNAO. Charged particles are easy to detect and can be back-tracked to the emission point with high efficiency in an almost background-free environment. These fragments are the product of projectiles fragmentation, and are hence mainly produced along the beam path inside the patient. This experimental signature can be used to monitor the beam position in the plane orthogonal to its flight direction, providing an online feedback to the beam transverse position monitor chambers used in the clinical centres. This information could be used to cross-check, validate and calibrate, whenever needed, the information provided by the ion chambers already implemented in most clinical centres as beam control detectors. In this paper we study the feasibility of such strategy in the clinical routine, analysing the data collected during the clinical trial performed at the CNAO facility on patients treated using 12C ions and monitored using the Dose Profiler (DP) detector developed within the INSIDE project. On the basis of the data collected monitoring three patients, the technique potential and limitations will be discussed.

Technical Design Report for a Carbon-11 Treatment Facility.

Particle therapy relies on the advantageous dose deposition which permits to highly conform the dose to the target and better spare the surrounding healthy tissues and organs at risk with respect to conventional radiotherapy. In the case of treatments with heavier ions (like carbon ions already clinically used), another advantage is the enhanced radiobiological effectiveness due to high linear energy transfer radiation. These particle therapy advantages are unfortunately not thoroughly exploited due to particle range uncertainties. The possibility to monitor the compliance between the ongoing and prescribed dose distribution is a crucial step toward new optimizations in treatment planning and adaptive therapy. The Positron Emission Tomography (PET) is an established quantitative 3D imaging technique for particle treatment verification and, among the isotopes used for PET imaging, the 11C has gained more attention from the scientific and clinical communities for its application as new radioactive projectile for particle therapy. This is an interesting option clinically because of an enhanced imaging potential, without dosimetry drawbacks; technically, because the stable isotope 12C is successfully already in use in clinics. The MEDICIS-Promed network led an initiative to study the possible technical solutions for the implementation of 11C radioisotopes in an accelerator-based particle therapy center. We present here the result of this study, consisting in a Technical Design Report for a 11C Treatment Facility. The clinical usefulness is reviewed based on existing experimental data, complemented by Monte Carlo simulations using the FLUKA code. The technical analysis starts from reviewing the layout and results of the facilities which produced 11C beams in the past, for testing purposes. It then focuses on the elaboration of the feasible upgrades of an existing 12C particle therapy center, to accommodate the production of 11C beams for therapy. The analysis covers the options to produce the 11C atoms in sufficient amounts (as required for therapy), to ionize them as required by the existing accelerator layouts, to accelerate and transport them to the irradiation rooms. The results of the analysis and the identified challenges define the possible implementation scenario and timeline.

Carbon ions beam therapy monitoring with the INSIDE in-beam PET.

In vivo range monitoring techniques are necessary in order to fully take advantage of the high dose gradients deliverable in hadrontherapy treatments. Positron emission tomography (PET) scanners can be used to monitor beam-induced activation in tissues and hence measure the range. The INSIDE (Innovative Solutions for In-beam DosimEtry in Hadrontherapy) in-beam PET scanner, installed at the Italian National Center of Oncological Hadrontherapy (CNAO, Pavia, Italy) synchrotron facility, has already been successfully tested in vivo during a proton therapy treatment. We discuss here the system performance evaluation with carbon ion beams, in view of future in vivo tests. The work is focused on the analysis of activity images obtained with therapeutic treatments delivered to polymethyl methacrylate (PMMA) phantoms, as well as on the test of an innovative and robust Monte Carlo simulation technique for the production of reliable prior activity maps. Images are reconstructed using different integration intervals, so as to monitor the activity evolution during and after the treatment. Three procedures to compare activity images are presented, namely Pearson correlation coefficient, Beam's eye view and overall view. Images of repeated irradiations of the same treatments are compared to assess the integration time necessary to provide reproducible images. The range agreement between simulated and experimental images is also evaluated, so as to validate the simulation capability to provide sound prior information. The results indicate that at treatment end, or at most 20 s afterwards, the range measurement is reliable within 1-2 mm, when comparing both different experimental sessions and data with simulations. In conclusion, this work shows that the INSIDE in-beam PET scanner performance is promising towards its in vivo test with carbon ions.

Online proton therapy monitoring: clinical test of a Silicon-photodetector-based in-beam PET.

Particle therapy exploits the energy deposition pattern of hadron beams. The narrow Bragg Peak at the end of range is a major advantage but range uncertainties can cause severe damage and require online verification to maximise the effectiveness in clinics. In-beam Positron Emission Tomography (PET) is a non-invasive, promising in-vivo technique, which consists in the measurement of the ß+ activity induced by beam-tissue interactions during treatment, and presents the highest correlation of the measured activity distribution with the deposited dose, since it is not much influenced by biological washout. Here we report the first clinical results obtained with a state-of-the-art in-beam PET scanner, with on-the-fly reconstruction of the activity distribution during irradiation. An automated time-resolved quantitative analysis was tested on a lacrimal gland carcinoma case, monitored during two consecutive treatment sessions. The 3D activity map was reconstructed every 10 s, with an average delay between beam delivery and image availability of about 6 s. The correlation coefficient of 3D activity maps for the two sessions (above 0.9 after 120 s) and the range agreement (within 1 mm) prove the suitability of in-beam PET for online range verification during treatment, a crucial step towards adaptive strategies in particle therapy.

INSIDE in-beam positron emission tomography system for particle range monitoring in hadrontherapy.

The quality assurance of particle therapy treatment is a fundamental issue that can be addressed by developing reliable monitoring techniques and indicators of the treatment plan correctness. Among the available imaging techniques, positron emission tomography (PET) has long been investigated and then clinically applied to proton and carbon beams. In 2013, the Innovative Solutions for Dosimetry in Hadrontherapy (INSIDE) collaboration proposed an innovative bimodal imaging concept that combines an in-beam PET scanner with a tracking system for charged particle imaging. This paper presents the general architecture of the INSIDE project but focuses on the in-beam PET scanner that has been designed to reconstruct the particles range with millimetric resolution within a fraction of the dose delivered in a treatment of head and neck tumors. The in-beam PET scanner has been recently installed at the Italian National Center of Oncologic Hadrontherapy (CNAO) in Pavia, Italy, and the commissioning phase has just started. The results of the first beam test with clinical proton beams on phantoms clearly show the capability of the in-beam PET to operate during the irradiation delivery and to reconstruct on-line the beam-induced activity map. The accuracy in the activity distal fall-off determination is millimetric for therapeutic doses.

Predictive Models Based on Support Vector Machines: Whole-Brain versus Regional Analysis of Structural MRI in the Alzheimer's Disease.

Decision-making systems trained on structural magnetic resonance imaging data of subjects affected by the Alzheimer's disease (AD) and healthy controls (CTRL) are becoming widespread prognostic tools for subjects with mild cognitive impairment (MCI). This study compares the performances of three classification methods based on support vector machines (SVMs), using as initial sets of brain voxels (ie, features): (1) the segmented grey matter (GM); (2) regions of interest (ROIs) by voxel-wise t-test filtering; (3) parceled ROIs, according to prior knowledge. The recursive feature elimination (RFE) is applied in all cases to investigate whether feature reduction improves the classification accuracy. We analyzed more than 600 AD Neuroimaging Initiative (ADNI) subjects, training the SVMs on the AD/CTRL dataset, and evaluating them on a trial MCI dataset. The classification performance, evaluated as the area under the receiver operating characteristic curve (AUC), reaches AUC = (88.9 ± .5)% in 20-fold cross-validation on the AD/CTRL dataset, when the GM is classified as a whole. The highest discrimination accuracy between MCI converters and nonconverters is achieved when the SVM-RFE is applied to the whole GM: with AUC reaching (70.7 ± .9)%, it outperforms both ROI-based approaches in predicting the AD conversion.