Particle beam microstructure reconstruction and coincidence discrimination in PET monitoring for hadron therapy.

Positron emission tomography is one of the most mature techniques for monitoring the particles range in hadron therapy, aiming to reduce treatment uncertainties and therefore the extent of safety margins in the treatment plan. In-beam PET monitoring has been already performed using inter-spill and post-irradiation data, i.e. while the particle beam is off or paused. The full beam acquisition procedure is commonly discarded because the particle spills abruptly increase the random coincidence rates and therefore the image noise. This is because random coincidences cannot be separated by annihilation photons originating from radioactive decays and cannot be corrected with standard random coincidence techniques due to the time correlation of the beam-induced background with the ion beam microstructure. The aim of this paper is to provide a new method to recover in-spill data to improve the images obtained with full-beam PET acquisitions. This is done by estimating the temporal microstructure of the beam and thus selecting input PET events that are less likely to be random ones. The PET detector we used was the one developed within the INSIDE project and tested at the CNAO synchrotron-based facility. The data were taken on a PMMA phantom irradiated with 72 MeV proton pencil beams. The obtained results confirm the possibility of improving the acquired PET data without any external signal coming from the synchrotron or ad hoc detectors.

Monte Carlo simulation tool for online treatment monitoring in hadrontherapy with in-beam PET: A patient study.

Hadrontherapy is a method for treating cancer with very targeted dose distributions and enhanced radiobiological effects. To fully exploit these advantages, in vivo range monitoring systems are required. These devices measure, preferably during the treatment, the secondary radiation generated by the beam-tissue interactions. However, since correlation of the secondary radiation distribution with the dose is not straightforward, Monte Carlo (MC) simulations are very important for treatment quality assessment. The INSIDE project constructed an in-beam PET scanner to detect signals generated by the positron-emitting isotopes resulting from projectile-target fragmentation. In addition, a FLUKA-based simulation tool was developed to predict the corresponding reference PET images using a detailed scanner model. The INSIDE in-beam PET was used to monitor two consecutive proton treatment sessions on a patient at the Italian Center for Oncological Hadrontherapy (CNAO). The reconstructed PET images were updated every 10 s providing a near real-time quality assessment. By half-way through the treatment, the statistics of the measured PET images were already significant enough to be compared with the simulations with average differences in the activity range less than 2.5 mm along the beam direction. Without taking into account any preferential direction, differences within 1 mm were found. In this paper, the INSIDE MC simulation tool is described and the results of the first in vivo agreement evaluation are reported. These results have justified a clinical trial, in which the MC simulation tool will be used on a daily basis to study the compliance tolerances between the measured and simulated PET images.

Proton range monitoring with in-beam PET: Monte Carlo activity predictions and comparison with cyclotron data.

GOAL: Proton treatment monitoring with Positron-Emission-Tomography (PET) is based on comparing measured and Monte Carlo (MC) predicted ß(+) activity distributions. Here we present PET ß(+) activity data and MC predictions both during and after proton irradiation of homogeneous PMMA targets, where protons were extracted from a cyclotron. METHODS AND MATERIALS: PMMA phantoms were irradiated with 62 MeV protons extracted from the CATANA cyclotron. PET activity data were acquired with a 10 × 10 cm(2) planar PET system and compared with predictions from the FLUKA MC generator. We investigated which isotopes are produced and decay during irradiation, and compared them to the situation after irradiation. For various irradiation conditions we compared one-dimensional activity distributions of MC and data, focussing on Δw50%, i.e., the distance between the 50% rise and 50% fall-off position. RESULTS: The PET system is able to acquire data during and after cyclotron irradiation. For PMMA phantoms the difference between the FLUKA MC prediction and our data in Δw50% is less than 1 mm. The ratio of PET activity events during and after irradiation is about 1 in both data and FLUKA, when equal time-frames are considered. Some differences are observed in profile shape. CONCLUSION: We found a good agreement in Δw50% and in the ratio between beam-on and beam-off activity between the PET data and the FLUKA MC predictions in all irradiation conditions.

First full-beam PET acquisitions in proton therapy with a modular dual-head dedicated system.

During particle therapy irradiation, positron emitters with half-lives ranging from 2 to 20 min are generated from nuclear processes. The half-lives are such that it is possible either to detect the positron signal in the treatment room using an in-beam positron emission tomography (PET) system, right after the irradiation, or to quickly transfer the patient to a close PET/CT scanner. Since the activity distribution is spatially correlated with the dose, it is possible to use PET imaging as an indirect method to assure the quality of the dose delivery. In this work, we present a new dedicated PET system able to operate in-beam. The PET apparatus consists in two 10 cm × 10 cm detector heads. Each detector is composed of four scintillating matrices of 23 × 23 LYSO crystals. The crystal size is 1.9 mm × 1.9 mm × 16 mm. Each scintillation matrix is read out independently with a modularized acquisition system. The distance between the two opposing detector heads was set to 20 cm. The system has very low dead time per detector area and a 3 ns coincidence window, which is capable to sustain high single count rates and to keep the random counts relatively low. This allows a new full-beam monitoring modality that includes data acquisition also while the beam is on. The PET system was tested during the irradiation at the CATANA (INFN, Catania, Italy) cyclotron-based proton therapy facility. Four acquisitions with different doses and dose rates were analysed. In all cases the random to total coincidences ratio was equal or less than 25%. For each measurement we estimated the accuracy and precision of the activity range on a set of voxel lines within an irradiated PMMA phantom. Results show that the inclusion of data acquired during the irradiation, referred to as beam-on data, improves both the precision and accuracy of the range measurement with respect to data acquired only after irradiation. Beam-on data alone are enough to give precisions better than 1 mm when at least 5 Gy are delivered.

Analysis of image sharpness reproducibility on a novel engineered micro-CT scanner with variable geometry and embedded recalibration software.

This study investigates the reproducibility of the reconstructed image sharpness, after modifications of the geometry setup, for a variable magnification micro-CT (µCT) scanner. All the measurements were performed on a novel engineered µCT scanner for in vivo imaging of small animals (Xalt), which has been recently built at the Institute of Clinical Physiology of the National Research Council (IFC-CNR, Pisa, Italy), in partnership with the University of Pisa. The Xalt scanner is equipped with an integrated software for on-line geometric recalibration, which will be used throughout the experiments. In order to evaluate the losses of image quality due to modifications of the geometry setup, we have made 22 consecutive acquisitions by changing alternatively the system geometry between two different setups (Large FoV - LF, and High Resolution - HR). For each acquisition, the tomographic images have been reconstructed before and after the on-line geometric recalibration. For each reconstruction, the image sharpness was evaluated using two different figures of merit: (i) the percentage contrast on a small bar pattern of fixed frequency (f = 5.5 lp/mm for the LF setup and f = 10 lp/mm for the HR setup) and (ii) the image entropy. We have found that, due to the small-scale mechanical uncertainty (in the order of the voxel size), a recalibration is necessary for each geometric setup after repositioning of the system's components; the resolution losses due to the lack of recalibration are worse for the HR setup (voxel size = 18.4 µm). The integrated on-line recalibration algorithm of the Xalt scanner allowed to perform the recalibration quickly, by restoring the spatial resolution of the system to the reference resolution obtained after the initial (off-line) calibration.

Comparison of two dedicated 'in beam' PET systems via simultaneous imaging of (12)C-induced beta(+)-activity.

The selective energy deposition of hadrontherapy has led to a growing interest in quality assurance techniques such as 'in-beam' PET. Due to the current lack of commercial solutions, dedicated detectors need to be developed. In this paper, we compare the performances of two different 'in-beam' PET systems which were simultaneously operated during and after low energy carbon ion irradiation of PMMA phantoms at GSI Darmstadt. The results highlight advantages and drawbacks of a novel in-beam PET prototype against a long-term clinically operated tomograph for ion therapy monitoring.

An optimization-based method for geometrical calibration in cone-beam CT without dedicated phantoms.

In this paper we present a new method for the determination of geometrical misalignments in cone-beam CT scanners, from the analysis of the projection data of a generic object. No a priori knowledge of the object shape and positioning is required. We show that a cost function, which depends on the misalignment parameters, can be defined using the projection data and that such a cost function has a local minimum in correspondence to the actual parameters of the system. Hence, the calibration of the scanner can be carried out by minimizing the cost function using standard optimization techniques. The method is developed for a particular class of 3D object functions, for which the redundancy of the fan beam sinogram in the transaxial midplane can be extended to cone-beam projection data, even at wide cone angles. The method has an approximated validity for objects which do not belong to that class; in that case, a suitable subset of the projection data can be selected in order to compute the cost function. We show by numerical simulations that our method is capable to determine with high accuracy the most critical misalignment parameters of the scanner, i.e., the transversal shift and the skew of the detector. Additionally, the detector slant can be determined. Other parameters such as the detector tilt, the longitudinal shift and the error in the source-detector distance cannot be determined with our method, as the proposed cost function has a very weak dependence on them. However, due to the negligible influence of these latter parameters in the reconstructed image quality, they can be kept fixed at estimated values in both calibration and reconstruction processes without compromising the final result. A trade-off between computational cost and calibration accuracy must be considered when choosing the data subset used for the computation of the cost function. Results on real data of a mouse femur as obtained with a small animal micro-CT are shown as well, proving the capability of the proposed calibration method. In principle, the method can be adapted to other cone-beam imaging modalities (e.g., single photon emission computed tomography).

Experimental validation of the filtering approach for dose monitoring in proton therapy at low energy.

The higher physical selectivity of proton therapy demands higher accuracy in monitoring of the delivered dose, especially when the target volume is located next to critical organs and a fractionated therapy is applied. A method to verify a treatment plan and to ensure the high quality of the hadrontherapy is to use Positron Emission Tomography (PET), which takes advantage of the nuclear reactions between protons and nuclei in the tissue during irradiation producing beta(+)-emitting isotopes. Unfortunately, the PET image is not directly proportional to the delivered radiation dose distribution; this is the reason why, at the present time, the verification of depth dose profiles with PET techniques is limited to a comparison between the measured activity and the one predicted for the planned treatment by a Monte Carlo model. In this paper we test the feasibility of a different scheme, which permits to reconstruct the expected PET signal from the planned radiation dose distribution along beam direction in a simpler and more direct way. The considered filter model, based on the description of the PET image as a convolution of the dose distribution with a filter function, has already demonstrated its potential applicability to beam energies above 70 MeV. Our experimental investigation provides support to the possibility of extending the same approach to the lower energy range ([40, 70] MeV), in the perspective of its clinical application in eye proton therapy.

Fast 3D-EM reconstruction using Planograms for stationary planar positron emission mammography camera.

At the University of Pisa we are building a PEM prototype, the YAP-PEM camera, consisting of two opposite 6 x 6 x 3 cm3 detector heads of 30 x 30 YAP:Ce finger crystals, 2 x 2 x 30 mm3 each. The camera will be equipped with breast compressors. The acquisition will be stationary. Compared with a whole body PET scanner, a planar Positron Emission Mammography (PEM) camera allows a better, easier and more flexible positioning around the breast in the vicinity of the tumor: this increases the sensitivity and solid angle coverage, and reduces cost. To avoid software rejection of data during the reconstruction, resulting in a reduced sensitivity, we adopted a 3D-EM reconstruction which uses all of the collected Lines Of Response (LORs). This skips the PSF distortion given by data rebinning procedures and/or Fourier methods. The traditional 3D-EM reconstruction requires several times the computation of the LOR-voxel correlation matrix, or probability matrix {p(ij)}; therefore is highly time-consuming. We use the sparse and symmetry properties of the matrix {p(ij)} to perform fast 3D-EM reconstruction. Geometrically, a 3D grid of cubic voxels (FOV) is crossed by several divergent 3D line sets (LORs). The symmetries occur when tracing different LORs produces the same p(ij) value. Parallel LORs of different sets cross the FOV in the same way, and the repetition of p(ij) values depends on the ratio between the tube and voxel sizes. By optimizing this ratio, the occurrence of symmetries is increased. We identify a nucleus of symmetry of LORs: for each set of symmetrical LORs we choose just one LOR to be put in the nucleus, while the others lie outside. All of the possible p(ij) values are obtainable by tracking only the LORs of this nucleus. The coordinates of the voxels of all of the other LORs are given by means of simple translation rules. Before making the reconstruction, we trace the LORs of the nucleus to find the intersecting voxels, whose p(ij) values are computed and stored with their voxel coordinates on a hard disk. Only the non-zero p(ij) are considered and their computation is performed just once. During the reconstruction, the stored values are loaded and are available in the random access memory for all of the operations of normalization, backprojection and projection: these are now performed rapidly, because the application of the translation rules is much faster than the probability computations. We tested the algorithm on Monte Carlo data fully simulating the typical YAP-PEM clinical condition. The adopted algorithm gives an excellent positioning capability for hot spots in the camera FOV. To use all of the possible skew LORs in the FOV avoids the software rejection of collected data. Reconstructed images indicate that a 5mm diameter tumor of 37 kBq/cm3, in an active breast with a 10:1 Tissue to Background ratio (T/B), with a 10 min acquisition, for a head distance of 5 cm, can be detected by the YAP-PEM with a SNR of 8.7+/-1.0. The obtained SNR values depend linearly on the tumor volume. The algorithm allows one to discriminate between two hot sources of 5.0 mm diameter if they do not lie on the same axis. The YAP-PEM is now in the assembly stage.

Advances in animal PET scanners.

Nowadays, a growing number of research groups shows a great interest for the application of PET and SPECT techniques to the development of new drugs. Preliminary studies on small animals require high performance dedicated scanners with a higher spatial resolution and sensitivity than those of clinical systems. In this paper the potential applications of such innovative instruments are shown together with a brief review of the dedicated PET and SPECT tomographs developed worldwide. Most of the scanners have been built as research prototypes. Only two are commercially available: micro-PET(R), designed and developed at UCLA, Los Angeles as a research prototype, and now produced and distributed by Concorde Microsystems Inc. (USA) and HIDAC PET produced by Oxford Positron Systems Ltd. (UK). Also in Italy, a high performance tomograph, YAP-(S) PET able to perform both PET and SPECT studies, has been developed at the University of Ferrara. The technical characteristics and performance of this scanner are described. Tomographs with combined imaging techniques, such as PET/CT or SPECT/CT, are now under study in various international research centers. The advantages of this new generation of animal scanners will be briefly outlined.

Coulomb energy differences in t = 1 mirror rotational bands in (50)Fe and (50)Cr.

Gamma rays from the N = Z-2 nucleus (50)Fe have been observed, establishing the rotational ground state band up to the state J(pi) = 11+ at 6.994 MeV excitation energy. The experimental Coulomb energy differences, obtained by comparison with the isobaric analog states in its mirror (50)Cr, confirm the qualitative interpretation of the backbending patterns in terms of successive alignments of proton and neutron pairs. A quantitative agreement with experiment has been achieved by exact shell model calculations, incorporating the differences in radii along the yrast bands, and properly renormalizing the Coulomb matrix elements in the pf model space.