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.

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).

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.