Impact of event positioning algorithm on performance of a whole-body PET scanner using one-to-one coupled detectors

Resolving inter-crystal scatter in a light-sharing depth-encoding PET detector

Objective.Inter-crystal scattering (ICS) in light-sharing positron emission tomography (PET) detectors leads to ambiguity in positioning the initial interaction, which significantly degrades the contrast, quantitative accuracy, and spatial resolution of the resulting image. Here, we attempt to resolve the positioning ambiguity of ICS in a light-sharing depth-encoding detector by exploiting the confined, deterministic light-sharing enabled by the segmented light guide unique to Prism-PET.Approach.We first considered a test case of ICS between two adjacent crystals using an analytical and a neural network approach. The analytical approach used a Bayesian estimation framework constructed from a scatter absorption model-the prior-and a detector response model-the likelihood. A simple neural network was generated for the same scenario, to provide mutual validation for the findings. Finally, we generalized the solution to three-dimensional event positioning that handles all events in the photopeak using a convolutional neural network with unique architecture that separately predicts the identity and depth-of-interaction (DOI) of the crystal containing the first interaction.Main results.The analytical Bayesian method generated an estimation error of 20.5 keV in energy and 3.1 mm in DOI. Further analysis showed that the detector response model was sufficiently robust to achieve adequate performance via maximum likelihood estimation (MLE), without prior information. We then found convergent results using a simple neural network. In the generalized solution using a convolutional neural network, we found crystal identification accuracy of 83% and DOI estimation error of 3.0 mm across all events. Applying this positioning algorithm to simulated data, we demonstrated significant improvements in image quality over the baseline, centroid-based positioning approach, attaining 38.9% improvement in intrinsic spatial resolution and enhanced clarity in hot spots of diameters 0.8 to 2.5 mm.Significance.The accuracy of our findings exceeds those of previous reports in the literature. The Prism-PET light guide, mediating confined and deterministic light-sharing, plays a key role in ICS recovery, as its mathematical embodiment-the detector response model-was the essential driver of accuracy in our results.

Technical opportunities and challenges in developing total-body PET scanners for mice and rats

Positron emission tomography (PET) is the most sensitive in vivo molecular imaging technique available. Small animal PET has been widely used in studying pharmaceutical biodistribution and disease progression over time by imaging a wide range of biological processes. However, it remains true that almost all small animal PET studies using mouse or rat as preclinical models are either limited by the spatial resolution or the sensitivity (especially for dynamic studies), or both, reducing the quantitative accuracy and quantitative precision of the results. Total-body small animal PET scanners, which have axial lengths longer than the nose-to-anus length of the mouse/rat and can provide high sensitivity across the entire body of mouse/rat, can realize new opportunities for small animal PET. This article aims to discuss the technical opportunities and challenges in developing total-body small animal PET scanners for mice and rats.

Simulation study of a brain PET scanner using TOF-DOI detectors equipped with first interaction position detection

Objective. The aim of this study is to evaluate the performance characteristics of a brain positron emission tomography (PET) scanner composed of four-layer independent read-out time-of-flight depth-of-interaction (TOF-DOI) detectors capable of first interaction position (FIP) detection, using Geant4 application for tomographic emission(GATE). This includes the spatial resolution, sensitivity, count rate capability, and reconstructed image quality.Approach. The proposed TOF-DOI PET detector comprises four layers of a 50 × 50 cerium-doped lutetium-yttrium oxyorthosilicate (LYSO:Ce) scintillator array with 1 mm pitch size, coupled to a 16 × 16 multi-pixel photon counter array with 3.0 mm × 3.0 mm photosensitive segments. Along the direction distant from the center field-of-view (FOV), the scintillator thickness of the four layers is 2.5, 3, 4, and 6 mm. The four layers were simulated with a 150 ps coincidence time resolution and the independent readout make the FIP detection capable. The spatial resolution and imaging performance were compared among the true-FIP, winner-takes-all (WTA) and front-layer FIP (FL-FIP) methods (FL-FIP selects the interaction position located on the front-most interaction layer in all the interaction layers). The National Electrical Manufacturers Association NU 2-2018 procedure was referred and modified to evaluate the performance of proposed scanner.Main results. In detector evaluation, the intrinsic spatial resolutions were 0.52 and 0.76 mm full width at half-maximum (FWHM) at 0° and 30° incidentγ-rays in the first layer pair, respectively. The reconstructed spatial resolution by the filter backprojection (FBP) achieved sub-millimeter FWHM on average over the whole FOV. The maximum true count rate was 207.6 kcps at 15 kBq ml^-1and the noise equivalent count rate (NECR_2R) was 54.7 kcps at 6.0 kBq ml^-1. Total sensitivity was 45.2 cps kBq^-1and 48.4 cps kBq^-1at the center and 10 cm off-center FOV, respectively. The TOF and DOI reconstructions significantly improved the image quality in the phantom studies. Moreover, the FL-FIP outperformed the conventional WTA method in terms of the spatial resolution and image quality.Significance. The proposed brain PET scanner could achieve sub-millimeter spatial resolution and high image quality with TOF and DOI reconstruction, which is meaningful to the clinical oncology research. Meanwhile, the comparison among the three positioning methods indicated that the FL-FIP decreased the image degradation caused by Compton scatter more than WTA.

Inter-crystal scattering recovery of light-sharing PET detectors using convolutional neural networks

Inter-crystal scattering (ICS) is a type of Compton scattering of photons from one crystal to adjacent crystals and causes inaccurate assignment of the annihilation photon interaction position in positron emission tomography (PET). Because ICS frequently occurs in highly light-shared PET detectors, its recovery is crucial for the spatial resolution improvement. In this study, we propose two different convolutional neural networks (CNNs) for ICS recovery, exploiting the good pattern recognition ability of CNN techniques. Using the signal distribution of a photosensor array as input, one network estimates the energy deposition in each crystal (ICS-eNet) and another network chooses the first-interacted crystal (ICS-cNet). We performed GATE Monte Carlo simulations with optical photon tracking to test PET detectors comprising different crystal arrays (8 × 8 to 21 × 21) with lengths of 20 mm and the same photosensor array (3 mm 8 × 8 array) covering an area of 25.8 × 25.8 mm². For each detector design, we trained ICS-eNet and ICS-cNet and evaluated their respective performance. ICS-eNet accurately identified whether the events were ICS (accuracy > 90%) and selected interacted crystals (accuracy > 60%) with appropriate energy estimation performance (R² > 0.7) in the 8 × 8, 12 × 12, and 16 × 16 arrays. ICS-cNet also exhibited satisfactory performance, which was less dependent on the crystal-to-sensor ratio, with an accuracy enhancement that exceeds 10% in selecting the first-interacted crystal and a reduction in error distances compared when no recovery was applied. Both ICS-eNet and ICS-cNet exhibited consistent performances under various optical property settings of the crystals. For spatial resolution measurements in PET rings, both networks achieved significant enhancements particularly for highly pixelated arrays. We also discuss approaches for training the networks in an actual experimental setup. This proof-of-concept study demonstrated the feasibility of CNNs for ICS recovery in various light-sharing designs to efficiently improve the spatial resolution of PET in various applications.

Design study of a PET detector with 0.5 mm crystal pitch for high-resolution preclinical imaging

Preclinical positron emission tomography (PET) is a sensitive and quantitative molecule imaging modality widely used in characterizing the biological processes and diseases in small animals. The purpose of this study is to investigate the methods to optimize a PET detector for high-resolution preclinical imaging. The PET detector proposed in this study consists of a 28 × 28 array of LYSO crystals 0.5 × 0.5 × 6.25 mm³in size, a wedged lightguide, and a 6 × 6 array of SiPMs 3 × 3 mm²in size. The simulation results showed that the most uniform flood map was achieved when the thickness of the lightguide was 2.35 mm. The quality of the flood map was significantly improved by suppressing the electronics noises using the simple threshold method with a best threshold. The peak-to-valley ratio of flood map improved 25.4% when the algorithm of ICS rejection was applied. An energy resolution (12.96% ± 1.03%) was measured on the prototype scanner constructed with 12 proposed detectors. Lastly, a prototype preclinic PET imager was constructed with 12 optimized detectors. The point source experiment was performed and an excellent spatial resolution (axial: 0.56 mm, tangential: 0.46 mm, radial: 0.42 mm) was achieved with the proposed high-performance PET detectors.

SiPM signal readout for inter-crystal scatter event identification in PET detectors

In positron emission tomography (PET) with pixelated detectors, a significant number of annihilation photons interact with scintillation crystals through single or multiple Compton scattering events. When these partial energy depositions occur across multiple crystal elements, we call them inter-crystal scatter (ICS) events. ICS events lead to incorrect localization of the annihilation photons, thereby degrading the PET image contrast, spatial resolution, and lesion detectability. The accurate identification of ICS events is the first essential step to improve the quality of PET images by rejecting ICS events or recovering ICS events without losing PET sensitivity. In this study, we propose a novel silicon photomultiplier (SiPM) readout method to identify ICS events in one-to-one coupled PET detectors with a reduced number of data acquisition channels. For concept verification, we assembled a PET detector that consists of a 16-channel SiPM array and 4 [Formula: see text] 4 lutetium oxyorthosilicate (LSO) array with a 3.2 mm crystal pitch. The proposed SiPM readout scheme serializes the 16 SiPM anode signals into four pulse train outputs encoded with four increasing time-delays in steps of 250 ns intervals. A Sum signal of the 16 SiPM anodes provides the timing information for time-of-flight measurement and a trigger signal for coincidence detection. A time-over-threshold (TOT) method was applied for obtaining the energy information followed by a subsequent TOT-to-energy calibration. We successfully identified the ICS events and determined their interacted positions and deposited energies by analyzing the digital pulses from the four pulse train output channels. The occurrence rate of ICS events was 10.85% for the 4 × 4 PET detector module with 3.2 mm-pitch LSO crystals. The PET detector yielded an energy resolution of 10.9 [Formula: see text] 0.6% and coincidence timing resolution of 285 [Formula: see text] 12 ps FWHM. We expect that the proposed method can be a useful solution for alleviating the readout burden of SiPM-based PET scanners with ICS event identification capability.

Recovery of inter-detector and inter-crystal scattering in brain PET based on LSO and GAGG crystals

Gadolinium aluminum gallium garnet (GAGG) is a promising scintillator crystal for positron emission tomography (PET) detectors owing to its advantages of energy resolution, light yield, and absence of intrinsic radiation. However, a large portion of the incident photons undergoes Compton scattering within GAGG crystal because of its low stopping power compared to that of lutetium-based crystals such as Lu₂SiO₅ (LSO). Inter-detector scattering (IDS) and inter-crystal scattering (ICS) result in loss of sensitivity and image quality of PET, respectively. We performed a Monte Carlo simulation study to evaluate IDS recovery in our currently developing brain-dedicated PET, and extended the idea to ICS recovery. We also compared the impact of the recoveries on LSO- and GAGG-based PET scanners. We measured the sensitivity and spatial resolution of the brain PET, and analyzed the image quality using a lesion phantom, a hot-rod phantom, and a 2D Hoffman phantom with applying IDS or ICS recovery. IDS recovery increased the PET sensitivity and improved the noise level of the reconstructed images. ICS recovery enhanced the spatial resolution and the contrast of the images was improved. As the occurrence rates of IDS and ICS were higher in GAGG than in LSO, the overall impact of IDS or ICS recovery was significant in GAGG. In conclusion, we showed that the proportional method would be suitable for IDS and ICS recoveries of PET, and emphasized the importance of ICS and IDS recoveries for PET using crystals with low stopping power.

Performance assessment of the 2 γpositronium imaging with the total-body PET scanners

Purpose

In living organisms, the positron-electron annihilation (occurring during the PET imaging) proceeds in about 30% via creation of a metastable ortho-positronium atom. In the tissue, due to the pick-off and conversion processes, over 98% of ortho-positronia annihilate into two 511 keV photons. In this article, we assess the feasibility for reconstruction of the mean ortho-positronium lifetime image based on annihilations into two photons. The main objectives of this work include the (i) estimation of the sensitivity of the total-body PET scanners for the ortho-positronium mean lifetime imaging using 2γ annihilations and (ii) estimation of the spatial and time resolution of the ortho-positronium image as a function of the coincidence resolving time (CRT) of the scanner.

Methods

Simulations are conducted assuming that radiopharmaceutical is labeled with ⁴⁴Sc isotope emitting one positron and one prompt gamma. The image is reconstructed on the basis of triple coincidence events. The ortho-positronium lifetime spectrum is determined for each voxel of the image. Calculations were performed for cases of total-body detectors build of (i) LYSO scintillators as used in the EXPLORER PET and (ii) plastic scintillators as anticipated for the cost-effective total-body J-PET scanner. To assess the spatial and time resolution, the four cases were considered assuming that CRT is equal to 500 ps, 140 ps, 50 ps, and 10 ps.

Results

The estimated total-body PET sensitivity for the registration and selection of image forming triple coincidences (2γ+γ_prompt) is larger by a factor of 13.5 (for LYSO PET) and by factor of 5.2 (for plastic PET) with respect to the sensitivity for the standard 2γ imaging by LYSO PET scanners with AFOV = 20 cm. The spatial resolution of the ortho-positronium image is comparable with the resolution achievable when using TOF-FBP algorithms already for CRT = 50 ps. For the 20-min scan, the resolution better than 20 ps is expected for the mean ortho-positronium lifetime image determination.

Conclusions

Ortho-positronium mean lifetime imaging based on the annihilations into two photons and prompt gamma is shown to be feasible with the advent of the high sensitivity total-body PET systems and time resolution of the order of tens of picoseconds.

The effects of inter-crystal scattering events on the performance of PET detectors

The probability of inter-crystal scattering (ICS) events for 511 keV gamma rays in all current scintillation crystals is high and the ICS events degrade the spatial resolution of PET scanners. In this work, Monte Carlo simulations were performed to study the effects of ICS events on the sensitivity and spatial resolution of PET detectors. LaBr₃, LYSO, and PWO that represent scintillation crystals of low, medium and high density, respectively, were used. For a point source placed in the middle of two scintillation detectors of 50  ×  50  ×  20 mm³ and a lower energy threshold (LET) of 350 keV, the probabilities that at least one gamma ray undergoes ICS are 94%, 84% and 76% for LaBr₃, LYSO, and PWO, respectively. The ICS events still provide useful spatial information. The full width at half maximum (FWHM), the full width at tenth maximum (FWTM) and the mean absolute error (MAE) of the curve of the mispositioning of a point source caused by ICS events are 0.45, 3.0 and 0.9 mm if the most popular PET scintillator LYSO is used. The MAE is smaller than the spatial resolution of most current PET scanners. The effect of ICS increases as the detector LET increases, scintillator density decreases, and crystal size decreases. The intrinsic spatial resolutions of a pair of LYSO detectors were calculated using curves of the coincidence counts between one column of the crystals in the two detectors and the sum of the coincidence counts between two opposite crystals of the columns in the two detectors that are in line with the point source changing with the source positions. The latter method removes almost all of ICS events. The FWHM (FWTM) intrinsic spatial resolutions obtained by the two methods are 0.40 (2.0) mm and 0.33 (0.8) mm if the crystal size is 0.5 mm, and are 0.8 (3.0) and 0.68 (1.5) mm if the crystal size is 1.0 mm. ICS events have much bigger contributions to the FWTM rather than the FWHM of the intrinsic spatial resolution of PET detectors. The spatial resolution of a PET scanner can still be improved by decreasing the crystal size to as small as 0.5 mm.

Intercrystal scatter studies for a 1 mm3 resolution clinical PET system prototype

Positron emission tomography (PET) systems designed with multiplexed readout do not usually have the capability to resolve individual intercrystal scatter (ICS) interactions, leading to interaction mispositioning that degrades spatial resolution and contrast. A 3D position sensitive scintillation detector capable of individual ICS readout has been designed and incorporated into a 1 mm³ resolution clinical PET system used for locoregional imaging. Incorporating ICS events increases photon sensitivity by 51.5% compared to using only photoelectric events. A Compton scatter angle error minimization algorithm is used to estimate the first ICS interaction location for accurate line-of-response pairing of coincident photons. An optimal scatter angle error threshold of 15 degrees is used to discard ICS events with a high mismatch between energy-derived and position-derived intercrystal scatter angles. Finally, positioning rather than rejecting ICS events boosts peak contrast to noise ratio by 8.1%, and allows for an equivalent dose reduction of 12% while maintaining equivalent image quality.

Highly Integrated FPGA-Only Signal Digitization Method Using Single-Ended Memory Interface Input Receivers for Time-of-Flight PET Detectors

We propose a new highly integrated field-programm-able gate array (FPGA) only signal digitization method for individual signal digitization of time-of-flight positron emission tomography (TOF PET). We configured I/O port of the FPGA with a single-ended memory interface (SeMI) input receiver. The SeMI is a single-ended voltage-referenced interface that has a common reference voltage per I/O Bank, such that each SeMI input receiver can serve as a voltage comparator. The FPGA-only digitizer that uses the single-ended input receivers does not require a separate digitizing integrated chip, and can obtain twice as many signals as that using LVDS input receivers. We implemented a highly integrated digitizer consisting of 82 energy and 82 timing channels using a 28-nm FPGA. The energy and arrival time were measured using a 625-ps binary counter, and a 10-ps time-to-digital converter (TDC), respectively. We first measured the intrinsic characteristics of the proposed FPGA-only digitizer. The SeMI input receiver functioned as the voltage comparator without undesirable offset voltage. The standard deviation value of the time difference measured using two SeMI input receivers with respective TDCs was less than 14.6 ps RMS. In addition, we fed signals from the TOF PET detectors to the SeMI input receivers directly and collected data. The TOF PET detector consisted of a 3 × 3 × 20 mm³ LYSO crystal coupled with a silicon photomultiplier. The energy resolutions were 7.7% and 7.1% for two TOF PET detectors. The coincidence resolving time was 204 ps full width at half maximum. The SeMI digitizer with a high-performance signal digitizer, processor, and high-speed transceivers provides a compact all-in-one data acquisition system.