Accurate and efficient modeling of the detector response in small animal multi-head PET systems.

In fully three-dimensional PET imaging, iterative image reconstruction techniques usually outperform analytical algorithms in terms of image quality provided that an appropriate system model is used. In this study we concentrate on the calculation of an accurate system model for the YAP-(S)PET II small animal scanner, with the aim to obtain fully resolution- and contrast-recovered images at low levels of image roughness. For this purpose we calculate the system model by decomposing it into a product of five matrices: (1) a detector response component obtained via Monte Carlo simulations, (2) a geometric component which describes the scanner geometry and which is calculated via a multi-ray method, (3) a detector normalization component derived from the acquisition of a planar source, (4) a photon attenuation component calculated from x-ray computed tomography data, and finally, (5) a positron range component is formally included. This system model factorization allows the optimization of each component in terms of computation time, storage requirements and accuracy. The main contribution of this work is a new, efficient way to calculate the detector response component for rotating, planar detectors, that consists of a GEANT4 based simulation of a subset of lines of flight (LOFs) for a single detector head whereas the missing LOFs are obtained by using intrinsic detector symmetries. Additionally, we introduce and analyze a probability threshold for matrix elements of the detector component to optimize the trade-off between the matrix size in terms of non-zero elements and the resulting quality of the reconstructed images. In order to evaluate our proposed system model we reconstructed various images of objects, acquired according to the NEMA NU 4-2008 standard, and we compared them to the images reconstructed with two other system models: a model that does not include any detector response component and a model that approximates analytically the depth of interaction as detector response component. The comparisons confirm previous research results, showing that the usage of an accurate system model with a realistic detector response leads to reconstructed images with better resolution and contrast recovery at low levels of image roughness.

Multi-ray-based system matrix generation for 3D PET reconstruction.

Iterative image reconstruction algorithms for positron emission tomography (PET) require a sophisticated system matrix (model) of the scanner. Our aim is to set up such a model offline for the YAP-(S)PET II small animal imaging tomograph in order to use it subsequently with standard ML-EM (maximum-likelihood expectation maximization) and OSEM (ordered subset expectation maximization) for fully three-dimensional image reconstruction. In general, the system model can be obtained analytically, via measurements or via Monte Carlo simulations. In this paper, we present the multi-ray method, which can be considered as a hybrid method to set up the system model offline. It incorporates accurate analytical (geometric) considerations as well as crystal depth and crystal scatter effects. At the same time, it has the potential to model seamlessly other physical aspects such as the positron range. The proposed method is based on multiple rays which are traced from/to the detector crystals through the image volume. Such a ray-tracing approach itself is not new; however, we derive a novel mathematical formulation of the approach and investigate the positioning of the integration (ray-end) points. First, we study single system matrix entries and show that the positioning and weighting of the ray-end points according to Gaussian integration give better results compared to equally spaced integration points (trapezoidal integration), especially if only a small number of integration points (rays) are used. Additionally, we show that, for a given variance of the single matrix entries, the number of rays (events) required to calculate the whole matrix is a factor of 20 larger when using a pure Monte-Carlo-based method. Finally, we analyse the quality of the model by reconstructing phantom data from the YAP-(S)PET II scanner.

Advanced radiation measurement techniques in diagnostic radiology and molecular imaging.

This paper reports some technological advances recently achieved in the fields of micro-CT and small animal PET instrumentation. It highlights a balance between image-quality improvement and dose reduction. Most of the recent accomplishments in these fields are due to the use of novel imaging sensors such as CMOS-based X-ray detectors and silicon photomultipliers (SiPM). Some of the research projects carried out at the University of Pisa for the development of such advanced radiation imaging technology are also described.

A detector head design for small-animal PET with silicon photomultipliers (SiPM).

Small-animal PET systems are now striving for sub-millimetre resolution. Current systems based upon PSPMTs and finely pixellated scintillators can be pushed to higher resolution, but at the expense of other performance parameters and a rapidly escalating cost. Moreover, depth of interaction (DOI) information is usually difficult to assess in such systems, even though this information is highly desirable to reduce the parallax error, which is often the dominant error for such high-resolution systems. In this study we propose a high-resolution detector head for a small-animal PET imaging system with intrinsic DOI information. Instead of a pixellated scintillator, our design is based upon the classic Anger camera principle, i.e. the head is constructed of modular layers each consisting of a continuous slab of scintillator, viewed by a new type of compact silicon photodetector. The photodetector is the recently developed silicon photomultiplier (SiPM) that as well as being very compact has many other attractive properties: high gain at low bias voltage, excellent single-photoelectron resolution and fast timing. A detector head of about 4 x 4 cm2 in area is proposed, constructed from three modular layers of the type described above. We perform a simulation study, using the Monte Carlo simulation package Geant4. The simulation results are used to optimize the geometry of the detector head and characterize its performance. Additionally, hit estimation algorithms are studied to determine the interaction position of annihilation photons correctly over the whole detector surface. The resulting detector has a nearly uniform efficiency for 511 keV photons of approximately 70% and an intrinsic spatial resolution of less than approximately 0.4 mm full width at half maximum (fwhm).