A CNN-based four-layer DOI encoding detector using LYSO and BGO scintillators for small animal PET imaging.

Objective. We propose a novel four-layer depth-of-interaction (DOI) encoding phoswich detector using lutetium-yttrium oxyothosilicate (LYSO) and bismuth germanate (BGO) scintillator crystal arrays for high sensitivity and high spatial resolution small animal PET imaging.Approach. The detector was comprised of a stack of four alternating LYSO and BGO scintillator crystal arrays coupled to an 8 × 8 multi-pixel photon counter (MPPC) array and read out by a PETsys TOFPET2 application specific integrated circuit. The four layers from the top (gamma ray entrance) to the bottom (facing the MPPC) consisted of a 24 × 24 array of 0.99 × 0.99 × 6 mm3LYSO crystals, a 24 × 24 array of 0.99 × 0.99 × 6 mm3BGO crystals, a 16 × 16 array of 1.53 × 1.53 × 6 mm3LYSO crystals and a 16 × 16 array of 1.53 × 1.53 × 6 mm3BGO crystals.Main results. Events that occurred in the LYSO and BGO layers were first separated by measuring the pulse energy (integrated charge) and duration (time over threshold (ToT)) from the scintillation pulses. Convolutional neural networks (CNNs) were then used to distinguish between the top and lower LYSO layers and between the upper and bottom BGO layers. Measurements with the prototype detector showed that our proposed method successfully identified events from all four layers. The CNN models achieved a classification accuracy of 91% for distinguishing the two LYSO layers and 81% for distinguishing the two BGO layers. The measured average energy resolution was 13.1% ± 1.7% for the top LYSO layer, 34.0% ± 6.3% for the upper BGO layer, 12.3% ± 1.3% for the lower LYSO layer, and 33.9% ± 6.9% for the bottom BGO layer. The timing resolution between each individual layer (from the top to the bottom) and a single crystal reference detector was 350 ps, 2.8 ns, 328 ps, and 2.1 ns respectively.Significance. In conclusion, the proposed four-layer DOI encoding detector achieved high performance and is an attractive choice for next-generation high sensitivity and high spatial resolution small animal positron emission tomography systems.

A digital phoswich detector using time-over-threshold for depth of interaction in PET.

We present the performance of a digital phoswich positron emission tomography (PET) detector, composed by layers of pixilated scintillator arrays, read out by solid state light detectors and an application specific integrated circuit (ASIC). We investigated the use of integrated charge from the scintillation pulses along with time-over-threshold (ToT) to determine the layer of interaction (DOI) in the scintillator. Simulations were performed to assess the effectiveness of the ToT measurements for separating the scintillator events and identifying cross-layer-crystal-scatter (CLCS) events. These simulations indicate that ToT and charge integration from such a detector provide sufficient information to determine the layer of interaction. To demonstrate this in practice, we used a pair of prototype LYSO/BGO detectors. One detector consisted of a 19 × 19 array of 7 mm long LYSO crystals (1.36 mm pitch) coupled to a 16 × 16 array of 8 mm long BGO crystals (1.63 mm pitch). The other detector was similar except the LYSO crystal pitch was 1.63 mm. These detectors were coupled to an 8 × 8 multi-pixel photon counter mounted on a PETsys TOFPET2 ASIC. This high performance ASIC provided digital readout of the integrated charge and ToT from these detectors. We present a method to separate the events from the two scintillator layers using the ToT, and also investigate the performance of this detector. All the crystals within the proposed detector were clearly resolved, and the peak to valley ratio was 11.8 ± 4.0 and 10.1 ± 2.9 for the LYSO and BGO flood images. The measured energy resolution was 9.9% ± 1.3% and 28.5% ± 5.0% respectively for the LYSO and BGO crystals in the phoswich layers. The timing resolution between the LYSO-LYSO, LYSO-BGO and BGO-BGO coincidences was 468 ps, 1.33 ns and 2.14 ns respectively. Results show ToT can be used to identify the crystal layer where events occurred and also identify and reject the majority of CLCS events between layers.

Performance evaluation of HiPET, a high sensitivity and high resolution preclinical PET tomograph.

HiPET is a recently developed prototype preclinical PET scanner dedicated to high sensitivity and high resolution molecular imaging. The HiPET system employs a phoswich depth of interaction (DOI) detector design, which also allows identification of the large majority of the cross layer crystal scatter (CLCS) events. This work evaluates its performance characteristics following the National Electrical Manufacturers Association (NEMA) NU4-2008 protocol. The HiPET consists of twenty flat panel type detectors arranged in two rings. The inner diameter is 160 mm and the axial field of view (FOV) is 104 mm. Each detector is comprised of two layers of phoswich scintillator crystal arrays, a tapered, pixelated glass lightguide and a multi anode photomultiplier tube (MAPMT). The front (gamma ray entrance) layer is a 48 × 48 pixelated cerium doped lutetium yttrium orthosilicate (LYSO) scintillator array with individual crystals measuring 1.01 × 1.01 × 6.1 mm. The back (towards the PMT) layer is a 32 × 32 pixelated bismuth germanate (BGO) scintillator array with individual crystals measuring 1.55 × 1.55 × 8.9 mm. For energy windows of 250-650 keV and 350-650 keV, the peak absolute sensitivity at the center of the FOV was 13.5% and 10.4% including CLCS events, and 11.8% and 8.9% excluding CLCS events, respectively. The average detector energy resolution derived by averaging the individual crystal spectra was 11.7% ± 1.4% for LYSO and 17.0% ± 1.4% for BGO. The 3D ordered-subsets expectation maximization (OSEM) reconstructed image of a point source in air, ranged from 0.73 mm to 1.19 mm, with an average value of 0.93 ± 0.09 mm at all measured locations. The peak noise equivalent count rate (NECR) and scatter fraction were 179 kcps at 12.4 MBq and 6.9% for the mouse-sized phantom, and 63 kcps at 11.3 MBq and 18.3% for the rat-sized phantom. For the NEMA image quality phantom, the uniformity was 5.8%, and the spillover ratios measured in the water- and air-filled cold region chambers were 0.047 and 0.044, respectively. The recovery coefficients (RC) ranged from 0.31 to 0.92. These results and in vivo evaluation demonstrate that the HiPET can achieve high quality molecular imaging for biomedical applications.

Performance Evaluation of G8, a High-Sensitivity Benchtop Preclinical PET/CT Tomograph.

G8 is a benchtop integrated PET/CT scanner dedicated to high-sensitivity and high-resolution imaging of mice. This work characterizes its National Electrical Manufacturers Association NU 4-2008 performance where applicable and also assesses the basic imaging performance of the CT subsystem. Methods: The PET subsystem in G8 consists of 4 flat-panel detectors arranged in a boxlike geometry. Each panel consists of 2 modules of a 26 × 26 pixelated bismuth germanate scintillator array with individual crystals measuring 1.75 × 1.75 × 7.2 mm. The crystal arrays are coupled to multichannel photomultiplier tubes via a tapered, pixelated glass lightguide. A cone-beam CT scanner consisting of a MicroFocus x-ray source and a complementary metal oxide semiconductor detector provides anatomic information. Sensitivity, spatial resolution, energy resolution, scatter fraction, count-rate performance, and the capability of performing phantom and mouse imaging were evaluated for the PET subsystem. Noise, dose level, contrast, and resolution were evaluated for the CT subsystem. Results: With an energy window of 350-650 keV, the peak sensitivity was 9.0% near the center of the field of view. The crystal energy resolution ranged from 15.0% to 69.6% in full width at half maximum (FWHM), with a mean of 19.3% ± 3.7%. The average intrinsic spatial resolution was 1.30 and 1.38 mm FWHM in the transverse and axial directions, respectively. The maximum-likelihood expectation maximization reconstructed image of a point source in air averaged 0.81 ± 0.11 mm FWHM. The peak noise-equivalent count rate for the mouse-sized phantom was 44 kcps for a total activity of 2.9 MBq (78 µCi), and the scatter fraction was 11%. For the CT subsystem, the value of the modulation transfer function at 10% was 2.05 cycles/mm. Conclusion: The overall performance demonstrates that the G8 can produce high-quality images for molecular imaging-based biomedical research.

MARS: a mouse atlas registration system based on a planar x-ray projector and an optical camera.

This paper introduces a mouse atlas registration system (MARS), composed of a stationary top-view x-ray projector and a side-view optical camera, coupled to a mouse atlas registration algorithm. This system uses the x-ray and optical images to guide a fully automatic co-registration of a mouse atlas with each subject, in order to provide anatomical reference for small animal molecular imaging systems such as positron emission tomography (PET). To facilitate the registration, a statistical atlas that accounts for inter-subject anatomical variations was constructed based on 83 organ-labeled mouse micro-computed tomography (CT) images. The statistical shape model and conditional Gaussian model techniques were used to register the atlas with the x-ray image and optical photo. The accuracy of the atlas registration was evaluated by comparing the registered atlas with the organ-labeled micro-CT images of the test subjects. The results showed excellent registration accuracy of the whole-body region, and good accuracy for the brain, liver, heart, lungs and kidneys. In its implementation, the MARS was integrated with a preclinical PET scanner to deliver combined PET/MARS imaging, and to facilitate atlas-assisted analysis of the preclinical PET images.

Investigation of OPET Performance Using GATE, a Geant4-Based Simulation Software.

A combined optical positron emission tomography (OPET) system is capable of both optical and PET imaging in the same setting, and it can provide information/interpretation not possible in single-mode imaging. The scintillator array here serves the dual function of coupling the optical signal from bioluminescence/fluorescence to the photodetector and also of channeling optical scintillations from the gamma rays. We report simulation results of the PET part of OPET using GATE, a Geant4 simulation package. The purpose of this investigation is the definition of the geometric parameters of the OPET tomograph. OPET is composed of six detector blocks arranged in a hexagonal ring-shaped pattern with an inner radius of 15.6 mm. Each detector consists of a two-dimensional array of 8 × 8 scintillator crystals each measuring 2 × 2 × 10 mm(3). Monte Carlo simulations were performed using the GATE software to measure absolute sensitivity, depth of interaction, and spatial resolution for two ring configurations, with and without gantry rotations, two crystal materials, and several crystal lengths. Images were reconstructed with filtered backprojection after angular interleaving and transverse one-dimensional interpolation of the sinogram. We report absolute sensitivities nearly seven times that of the prototype microPET at the center of field of view and 2.0 mm tangential and 2.3 mm radial resolutions with gantry rotations up to an 8.0 mm radial offset. These performance parameters indicate that the imaging spatial resolution and sensitivity of the OPET system will be suitable for high-resolution and high-sensitivity small-animal PET imaging.