Effects of multiple-interaction photon events in a high-resolution PET system that uses 3-D positioning detectors

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.

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.

Initial performance studies of a wearable brain positron emission tomography camera based on autonomous thin-film digital Geiger avalanche photodiode arrays

Using analytical and Monte Carlo modeling, we explored performance of a lightweight wearable helmet-shaped brain positron emission tomography (PET), or BET camera, based on thin-film digital Geiger avalanche photodiode arrays with Lutetium-yttrium oxyorthosilicate (LYSO) or [Formula: see text] scintillators for imaging in vivo human brain function of freely moving and acting subjects. We investigated a spherical cap BET and cylindrical brain PET (CYL) geometries with 250-mm diameter. We also considered a clinical whole-body (WB) LYSO PET/CT scanner. The simulated energy resolutions were 10.8% (LYSO) and 3.3% ([Formula: see text]), and the coincidence window was set at 2 ns. The brain was simulated as a water sphere of uniform F-18 activity with a radius of 100 mm. We found that BET achieved [Formula: see text] better noise equivalent count (NEC) performance relative to the CYL and [Formula: see text] than WB. For 10-mm-thick [Formula: see text] equivalent mass systems, LYSO (7-mm thick) had [Formula: see text] higher NEC than [Formula: see text]. We found that [Formula: see text] scintillator crystals achieved [Formula: see text] full-width-half-maximum spatial resolution without parallax errors. Additionally, our simulations showed that LYSO generally outperformed [Formula: see text] for NEC unless the timing resolution for [Formula: see text] was considerably smaller than that presently used for LYSO, i.e., well below 300 ps.

A DOI Detector With Crystal Scatter Identification Capability for High Sensitivity and High Spatial Resolution PET Imaging

A new phoswich detector is being developed at the Crump Institute, aiming to provide improvements in sensitivity, and spatial resolution for PET. The detector configuration is comprised of two layers of pixelated scintillator crystal arrays, a glass light guide and a light detector. The annihilation photon entrance (top) layer is a 48 × 48 array of 1.01 × 1.01 × 7 mm³ LYSO crystals. The bottom layer is a 32 × 32 array of 1.55 × 1.55 × 9 mm³ BGO crystals. A tapered, multiple-element glass lightguide is used to couple the exit end of the BGO crystal array (52 × 52 mm²) to the photosensitive area of the Position Sensitive Photomultiplier Tube (46 × 46 mm²), allowing the creation of flat panel detectors without gaps between the detector modules. Both simulations and measurements were performed to evaluate the characteristics and benefits of the proposed design. The GATE Monte Carlo simulation indicated that the total fraction of the cross layer crystal scatter (CLCS) events in singles detection mode for this detector geometry is 13.2%. The large majority of these CLCS events (10.1% out of 13.2%) deposit most of their energy in a scintillator layer other than the layer of first interaction. Identification of those CLCS events for rejection or correction may lead to improvements in data quality and imaging performance. Physical measurements with the prototype detector showed that the LYSO, BGO and CLCS events were successfully identified using the delayed charge integration (DCI) technique, with more than 95% of the LYSO and BGO crystal elements clearly resolved. The measured peak-to-valley ratios (PVR) in the flood histograms were 3.5 for LYSO and 2.0 for BGO. For LYSO, the energy resolution ranged from 9.7% to 37.0% full width at half maximum (FWHM), with a mean of 13.4 ± 4.8%. For BGO the energy resolution ranged from 16.0% to 33.9% FWHM, with a mean of 18.6 ± 3.2%. In conclusion, these results demonstrate that the proposed detector is feasible and can potentially lead to a high spatial resolution, high sensitivity and DOI PET system.

Recovery and normalization of triple coincidences in PET

PURPOSE

Triple coincidences in positron emission tomography (PET) are events in which three γ-rays are detected simultaneously. These events, though potentially useful for enhancing the sensitivity of PET scanners, are discarded or processed without special consideration in current systems, because there is not a clear criterion for assigning them to a unique line-of-response (LOR). Methods proposed for recovering such events usually rely on the use of highly specialized detection systems, hampering general adoption, and/or are based on Compton-scatter kinematics and, consequently, are limited in accuracy by the energy resolution of standard PET detectors. In this work, the authors propose a simple and general solution for recovering triple coincidences, which does not require specialized detectors or additional energy resolution requirements.

METHODS

To recover triple coincidences, the authors' method distributes such events among their possible LORs using the relative proportions of double coincidences in these LORs. The authors show analytically that this assignment scheme represents the maximum-likelihood solution for the triple-coincidence distribution problem. The PET component of a preclinical PET/CT scanner was adapted to enable the acquisition and processing of triple coincidences. Since the efficiencies for detecting double and triple events were found to be different throughout the scanner field-of-view, a normalization procedure specific for triple coincidences was also developed. The effect of including triple coincidences using their method was compared against the cases of equally weighting the triples among their possible LORs and discarding all the triple events. The authors used as figures of merit for this comparison sensitivity, noise-equivalent count (NEC) rates and image quality calculated as described in the NEMA NU-4 protocol for the assessment of preclinical PET scanners.

RESULTS

The addition of triple-coincidence events with the authors' method increased peak NEC rates of the scanner by 26.6% and 32% for mouse- and rat-sized objects, respectively. This increase in NEC-rate performance was also reflected in the image-quality metrics. Images reconstructed using double and triple coincidences recovered using their method had better signal-to-noise ratio than those obtained using only double coincidences, while preserving spatial resolution and contrast. Distribution of triple coincidences using an equal-weighting scheme increased apparent system sensitivity but degraded image quality. The performance boost provided by the inclusion of triple coincidences using their method allowed to reduce the acquisition time of standard imaging procedures by up to ∼25%.

CONCLUSIONS

Recovering triple coincidences with the proposed method can effectively increase the sensitivity of current clinical and preclinical PET systems without compromising other parameters like spatial resolution or contrast.

Study of electrode pattern design for a CZT-based PET detector

We are developing a 1 mm resolution small animal positron emission tomography (PET) system using 3D positioning cadmium zinc telluride photon detectors comprising 40 mm × 40 mm × 5 mm crystals metalized with a cross-strip electrode pattern with a 1 mm anode strip pitch. We optimized the electrode pattern design for intrinsic sensitivity and spatial, energy and time resolution performance using a test detector comprising cathode and steering electrode strips of varying dimensions. The study found 3 and 5 mm width cathode strips locate charge-shared photon interactions near cathode strip boundaries with equal precision. 3 mm width cathode strips exhibited large time resolution variability as a function of photon interaction location between the anode and cathode planes (~26 to ~127.5 ns full width at half maximum (FWHM) for 0.5 mm and 4.2 mm depths, respectively). 5 mm width cathode strips by contrast exhibited more stable time resolution for the same interaction locations (~34 to ~83 ns FWHM), provided more linear spatial positioning in the direction orthogonal to the electrode planes, and as much as 68.4% improvement in photon sensitivity over the 3 mm wide cathode strips. The results were understood by analyzing the cathode strips' weighting functions, which indicated a stronger 'small pixel' effect in the 3 mm wide cathode strips. Photon sensitivity and anode energy resolution were seen to improve with decreasing steering electrode bias from 0 to -80 V w.r.t. the anode potential. A slight improvement in energy resolution was seen for wider steering electrode strips (400 versus 100 µm) for charge-shared photon interactions. Although this study successfully focused on electrode pattern features for PET performance, the results are generally applicable to semiconductor photon detectors employing cross-trip electrode patterns.

Promising new photon detection concepts for high-resolution clinical and preclinical PET

The ability of PET to visualize and quantify regions of low concentration of PET tracer representing subtle cellular and molecular signatures of disease depends on relatively complex biochemical, biologic, and physiologic factors that are challenging to control, as well as on instrumentation performance parameters that are, in principle, still possible to improve on. Thus, advances to the latter can somewhat offset barriers of the former. PET system performance parameters such as spatial resolution, contrast resolution, and photon sensitivity contribute significantly to PET's ability to visualize and quantify lower concentrations of signal in the presence of background. In this report we present some technology innovations under investigation toward improving these PET system performance parameters. We focus particularly on a promising advance known as 3-dimensional position-sensitive detectors, which are detectors capable of distinguishing and measuring the position, energy, and arrival time of individual interactions of multi-interaction photon events in 3 dimensions. If successful, these new strategies enable enhancements such as the detection of fewer diseased cells in tissue or the ability to characterize lower-abundance molecular targets within cells. Translating these advanced capabilities to the clinic might allow expansion of PET's roles in disease management, perhaps to earlier stages of disease. In preclinical research, such enhancements enable more sensitive and accurate studies of disease biology in living subjects.

Online detector response calculations for high-resolution PET image reconstruction

Positron emission tomography systems are best described by a linear shift-varying model. However, image reconstruction often assumes simplified shift-invariant models to the detriment of image quality and quantitative accuracy. We investigated a shift-varying model of the geometrical system response based on an analytical formulation. The model was incorporated within a list-mode, fully 3D iterative reconstruction process in which the system response coefficients are calculated online on a graphics processing unit (GPU). The implementation requires less than 512 Mb of GPU memory and can process two million events per minute (forward and backprojection). For small detector volume elements, the analytical model compared well to reference calculations. Images reconstructed with the shift-varying model achieved higher quality and quantitative accuracy than those that used a simpler shift-invariant model. For an 8 mm sphere in a warm background, the contrast recovery was 95.8% for the shift-varying model versus 85.9% for the shift-invariant model. In addition, the spatial resolution was more uniform across the field-of-view: for an array of 1.75 mm hot spheres in air, the variation in reconstructed sphere size was 0.5 mm RMS for the shift-invariant model, compared to 0.07 mm RMS for the shift-varying model.

Study of a high-resolution, 3D positioning cadmium zinc telluride detector for PET

This paper investigates the performance of 1 mm resolution cadmium zinc telluride (CZT) detectors for positron emission tomography (PET) capable of positioning the 3D coordinates of individual 511 keV photon interactions. The detectors comprise 40 mm × 40 mm × 5 mm monolithic CZT crystals that employ a novel cross-strip readout with interspersed steering electrodes to obtain high spatial and energy resolution. The study found a single anode FWHM energy resolution of 3.06 ± 0.39% at 511 keV throughout most of the detector volume. Improved resolution is expected with properly shielded front-end electronics. Measurements made using a collimated beam established the efficacy of the steering electrodes in facilitating enhanced charge collection across anodes, as well as a spatial resolution of 0.44 ± 0.07 mm in the direction orthogonal to the electrode planes. Finally, measurements based on coincidence electronic collimation yielded a point spread function with 0.78 ± 0.10 mm FWHM, demonstrating 1 mm spatial resolution capability transverse to the anodes-as expected from the 1 mm anode pitch. These findings indicate that the CZT-based detector concept has excellent performance and shows great promise for a high-resolution PET system.

Physical effects of mechanical design parameters on photon sensitivity and spatial resolution performance of a breast-dedicated PET system

PURPOSE

This study aims to address design considerations of a high resolution, high sensitivity positron emission tomography scanner dedicated to breast imaging.

METHODS

The methodology uses a detailed Monte Carlo model of the system structures to obtain a quantitative evaluation of several performance parameters. Special focus was given to the effect of dense mechanical structures designed to provide mechanical robustness and thermal regulation to the minuscule and temperature sensitive detectors.

RESULTS

For the energies of interest around the photopeak (450-700 keV energy window), the simulation results predict a 6.5% reduction in the single photon detection efficiency and a 12.5% reduction in the coincidence photon detection efficiency in the case that the mechanical structures are interspersed between the detectors. However for lower energies, a substantial increase in the number of detected events (approximately 14% and 7% for singles at a 100-200 keV energy window and coincidences at a lower energy threshold of 100 keV, respectively) was observed with the presence of these structures due to backscatter. The number of photon events that involve multiple interactions in various crystal elements is also affected by the presence of the structures. For photon events involving multiple interactions among various crystal elements, the coincidence photon sensitivity is reduced by as much as 20% for a point source at the center of the field of view. There is no observable effect on the intrinsic and the reconstructed spatial resolution and spatial resolution uniformity.

CONCLUSIONS

Mechanical structures can have a considerable effect on system sensitivity, especially for systems processing multi-interaction photon events. This effect, however, does not impact the spatial resolution. Various mechanical structure designs are currently under evaluation in order to achieve optimum trade-off between temperature stability, accurate detector positioning, and minimum influence on system performance.

Analog signal multiplexing for PSAPD-based PET detectors: simulation and experimental validation

A 1 mm(3) resolution clinical positron emission tomography (PET) system employing 4608 position-sensitive avalanche photodiodes (PSAPDs) is under development. This paper describes a detector multiplexing technique that simplifies the readout electronics and reduces the density of the circuit board design. The multiplexing scheme was validated using a simulation framework that models the PSAPDs and front-end multiplexing circuits to predict the signal-to-noise ratio and flood histogram performance. Two independent experimental setups measured the energy resolution, time resolution, crystal identification ability and count rate both with and without multiplexing. With multiplexing, there was no significant degradation in energy resolution, time resolution and count rate. There was a relative 6.9 ± 1.0% and 9.4 ± 1.0% degradation in the figure of merit that characterizes the crystal identification ability observed in the measured and simulated ceramic-mounted PSAPD module flood histograms, respectively.