L. Hernandez

ALBIRA: A small animal PET/SPECT/CT imaging system

PURPOSE The authors have developed a trimodal PET∕SPECT∕CT scanner for small animal imaging. The gamma ray subsystems are based on monolithic crystals coupled to multianode photomultiplier tubes (MA-PMTs), while computed tomography (CT) comprises a commercially available microfocus x-ray tube and a CsI scintillator 2D pixelated flat panel x-ray detector. In this study the authors will report on the design and performance evaluation of the multimodal system. METHODS X-ray transmission measurements are performed based on cone-beam geometry. Individual projections were acquired by rotating the x-ray tube and the 2D flat panel detector, thus making possible a transaxial field of view (FOV) of roughly 80 mm in diameter and an axial FOV of 65 mm for the CT system. The single photon emission computed tomography (SPECT) component has a dual head detector geometry mounted on a rotating gantry. The distance between the SPECT module detectors can be varied in order to optimize specific user requirements, including variable FOV. The positron emission tomography (PET) system is made up of eight compact modules forming an octagon with an axial FOV of 40 mm and a transaxial FOV of 80 mm in diameter. The main CT image quality parameters (spatial resolution and uniformity) have been determined. In the case of the SPECT, the tomographic spatial resolution and system sensitivity have been evaluated with a (99m)Tc solution using single-pinhole and multi-pinhole collimators. PET and SPECT images were reconstructed using three-dimensional (3D) maximum likelihood and ordered subset expectation maximization (MLEM and OSEM) algorithms developed by the authors, whereas the CT images were obtained using a 3D based FBP algorithm. RESULTS CT spatial resolution was 85 μm while a uniformity of 2.7% was obtained for a water filled phantom at 45 kV. The SPECT spatial resolution was better than 0.8 mm measured with a Derenzo-like phantom for a FOV of 20 mm using a 1-mm pinhole aperture collimator. The full width at half-maximum PET radial spatial resolution at the center of the field of view was 1.55 mm. The SPECT system sensitivity for a FOV of 20 mm and 15% energy window was 700 cps∕MBq (7.8 × 10(-2)%) using a multi-pinhole equipped with five apertures 1 mm in diameter, whereas the PET absolute sensitivity was 2% for a 350-650 keV energy window and a 5 ns timing window. Several animal images are also presented. CONCLUSIONS The new small animal PET∕SPECT∕CT proposed here exhibits high performance, producing high-quality images suitable for studies with small animals. Monolithic design for PET and SPECT scintillator crystals reduces cost and complexity without significant performance degradation.

A PET Design Based on SiPM and Monolithic LYSO Crystals: Performance Evaluation

A new small animal PET based on SiPM and monolithic LYSO crystals has been developed. Eight detector modules form the PET ring, each mounting an array of 12 × 12 SiPMs coupled to a readout providing the summed signals of the pixels on each of the 12 rows and 12 columns of the SiPM array. This design makes it possible to accurately determine the centroid of the scintillation light distribution with about 1.6 mm full width at half maximum (FWHM) resolution without correction for the 1 mm source size, and the photon depth of interaction (DOI) with nearly 2 mm FWHM. This single ring PET system has a homogeneous spatial resolution across the entire 80 mm transaxial field of view (FOV) of about 1 mm FWHM. The noise equivalent count rate (NECR) peak is estimated to occur at around 39.2 MBq with a rate of approximately 82.7 kcps for the mouse-like phantom and 22 kcps at 48.1 MBq for the rat-like phantom. Following the NEMA protocol, the peak absolute sensitivity in the center of the FOV is 2.8% for a 30% peak energy window. A pilot test injecting NaF to a mouse of 20 grams is also presented. Finally, the PET ring has been tested in front of a high field 15.2 T Magnetic Resonance (MR). No significant variation on energy and spatial resolution across the FOV has been observed due to the presence of the magnetic field.

Performance Study of a Wide-Area SiPM Array, ASICs Controlled

In this paper, the capabilities of a wide-area gamma ray photosensor based on a SiPM array are investigated. For this purpose, we have mounted an array of 144 (12×12) SiPMs with individual active area of 3 ×3 mm2 and a pitch of 4.2 mm, thus covering an active area of 50.2 ×50.2 mm2. The measurements were performed by coupling the SiPM array to LYSO crystal arrays of different pixel size ( 2×2 mm2, 1.5 ×1.5 mm2, and 1 ×1 mm2) and 10-12 mm thicknesses. The SiPM array was controlled by means of three ASICs, and the SiPM signals were multiplexed in order to determine the gamma ray impact position by means of implementing the Anger logic algorithm in the ASIC. The optimum bias voltage and temperature dependence of the gamma ray sensor were determined. An energy resolution as good as 8%, for individual crystal pixels, were reached at 5 V overvoltage. The ASICs design allows one to “activate” different photosensor array areas. This feature has been used to evaluate the detector performance as a function of the crystal pixel size and the photosensor dark noise contribution. In this work we also show the system capability to provide depth-of-interaction (DOI) information by means of implementing a two-layer staggered approach. We have found that accurate DOI information is obtained when the ASICs enabled an SiPM active area as high as 32×32 mm2( 8 ×8 SiPMs).

Calibration and Performance Tests of Detectors for Laser-Accelerated Protons

We present the calibration and performance tests carried out with two detectors for intense proton pulses accelerated by lasers. Most of the procedures were realized with proton beams of 0.46-5.60 MeV from a tandem accelerator. One approach made use of radiochromic films, for which we calibrated the relation between optical density and energy deposition over more than three orders of magnitude. The validity of these results and of our analysis algorithms has been confirmed by controlled irradiation of film stacks and reconstruction of the total beam charge for strongly non-uniform beam profiles. For the spectral analysis of protons from repeated laser shots, we have designed an online monitor based on a plastic scintillator. The resulting signal from a photomultiplier directly measured on a fast oscilloscope is especially useful for time-of-flight applications. Variable optical filters allow for suppression of saturation and an extension of the dynamic range. With pulsed proton beams we have tested the detector response to a wide range of beam intensities from single particles to 3 ×105 protons per 100 ns time interval.

Timing Results Using an FPGA-Based TDC with Large Arrays of 144 SiPMs

Silicon photomultipliers (SiPMs) have become an alternative to traditional tubes due to several features. However, their implementation to form large arrays is still a challenge especially due to their relatively high intrinsic noise, depending on the chosen readout. In this contribution, two modules composed of 12 ×12 SiPMs with an area of roughly 50 mm×50 mm are used in coincidence. Coincidence resolving time (CRT) results with a field-programmable gate array, in combination with a time to digital converter, are shown as a function of both the sensor bias voltage and the digitizer threshold. The dependence of the CRT on the sensor matrix temperature, the amount of SiPM active area and the crystal type is also analyzed. Measurements carried out with a crystal array of 2 mm pixel size and 10 mm height have shown time resolutions for the entire 288 SiPM two-detector set-up as good as 800 ps full width at half maximum (FWHM).

3-D photon impact determination using fitting approaches to the Light Distribution

In Positron Emission Tomography (PET) detectors based in monolithic scintillators, the spatial resolution is limited by the accuracy in the determination of the interaction coordinates from the 511 keV photons. When linear algorithms, such as Center of Gravity (CoG) are used a poor estimation of the interaction positions, specially towards the edges is the major limitation in spatial resolution. A novel PET detector block, where complete information of Light Distribution (LD) for each event is available, allows to fit each event to a theoretical model, improving the estimation of the interaction coordinates, and minimizing border effects. In this work, by means of the LD fitting approach, we were able to obtain an average spatial resolution of 1.2 mm in the entire scintillator volume and an average depth of interaction (DOI) resolution of 1.5 mm. Moreover, splitting the data in three DOI regions, we obtained an average spatial resolution of 1.0 mm at the DOI region closer to the photodetectors. Finally, it is remarkable that the implementation of the LD fitting approach is capable of processing up to 50 kcps in a octacore system.

Expectation maximization (EM) algorithms using polar symmetries for computed tomography (CT) image reconstruction

We suggest a symmetric-polar pixellation scheme which makes possible a reduction of the computational cost for expectation maximization (EM) iterative algorithms. The proposed symmetric-polar pixellation allows us to deal with 3D images as a whole problem without dividing the 3D problem into 2D slices approach. Performance evaluation of each approach in terms of stability and image quality is presented. Exhaustive comparisons between all approaches were conducted in a 2D based image reconstruction model. From these 2D approaches, that showing the best performances were finally implemented and evaluated in a 3D based image reconstruction model. Comparison to 3D images reconstructed with FBP is also presented. Although the algorithm is presented in the context of computed tomography (CT) image reconstruction, it can be applied to any other tomographic technique as well, due to the fact that the only requirement is a scanning geometry involving measurements of an object under different projection angles. Real data have been acquired with a small animal (CT) scanner to verify the proposed mathematical description of the CT system.

EM tomographic image reconstruction using polar voxels

The splitting of the field of view (FOV) in polar voxels is pro posed in this work in order to obtain an efficient description of a cone-beam compu ted tomography (CT) scanner. The proposed symmetric-polar pixelation makes it possible to deal with the 3D iterative reconstruction considering a number of projections and voxel sizes typical in CT preclinical imaging. The performance comparison, between the filtered backproje ction (FBP) and 3D maximum likelihood expectation maximization (MLEM) reconstruction algorithm for CT, is presented. It is feasible to achieve the hardware spatial resolution limi t with the considered pixelation. The image quality achieved with MLEM and FBP have been analyzed. The results obtained with both algorithms in clinical images have been compared too. Although the polar-symmetric pixelation is presented in the context of CT imaging, it can be applied to any other tomographic technique as long as the scan comprises the measurement of an object under several projection angles.

Simulation Study of Resistor Networks Applied to an Array of 256 SiPMs

In this work we describe a procedure to reduce the number of signals detected by an array of 256 Silicon Photomultipliers (SiPMs) using a resistor network to divide the signal charge into few readout channels. Several configurations were modeled, and the pulsed signal at the readout contacts were simulated. These simulation results were experimentally tested on a specifically designed and manufactured set of printed circuit boards. Three network configurations were modeled. The modeling provided encouraging results for all three configurations. The measurements on the prototypes constructed for this study, however, provided useful position-sensitivity for only one of the network configurations. The lack of input signal amplification into the networks, the SiPM dark current, as well as the complexity of an eight layers board with parasitic capacitances, could have caused the degradation of resolving the impact photon position. This is hard to overcome with external printed circuit boards and components.

First results of an ASIC controlled γ-detector based on a SiPM-array and a monolithic LYSO

We present the first results of an MR compatible PET system based on continuous crystals and SiPM detectors, along with a dedicated application-specific integrated circuit (ASIC) for the readout stage. The use of continuous scintillation crystals preserves the spatial distribution of scintillation light for each γ--ray impact, which can be reconstructed with a small number of statistical moments Therefore, it is possible to reduce the number of analogue-to-digital conversion channels. The current ASIC can measure up to 8 moments of the light distribution, providing information about the X and Y photon impact coordinates, photon energy, depth of interaction and other higher moments that can improve the position accuracy of the photon impact. SiPMs exhibit their best performance for reduce active areas, where the dark counts are minimized. An innovate coupling design between crystal and SiPM-array has been made in order to keep the compromise between reduced active areas and efficient collection of the scintillation light. SiPM detectors are also suitable to work in the presence of magnetic fields. Moreover, they could be integrated in a RF coil of the MRI system due to their reduced dimensions [1]. Experimental results show very promising possibilities for the system.