Evaluation of event position reconstruction in monolithic crystals that are optically coupled

Evaluation of monolithic crystal detector with dual-ended readout utilizing multiplexing method

Objective.Monolithic crystal detectors are increasingly being applied in positron emission tomography (PET) devices owing to their excellent depth-of-interaction (DOI) resolution capabilities and high detection efficiency. In this study, we constructed and evaluated a dual-ended readout monolithic crystal detector based on a multiplexing method.Approach.We employed two 12 × 12 silicon photomultiplier (SiPM) arrays for readout, and the signals from the 12 × 12 array were merged into 12 X and 12 Y channels using channel multiplexing. In 2D reconstruction, three methods based on the centre of gravity (COG) were compared, and the concept of thresholds was introduced. Furthermore, a light convolutional neural network (CNN) was employed for testing. To enhance depth localization resolution, we proposed a method by utilizing the mutual information from both ends of the SiPMs. The source width and collimation effect were simulated using GEANT4, and the intrinsic spatial resolution was separated from the measured values.Main results.At an operational voltage of 29 V for the SiPM, an energy resolution of approximately 12.5 % was achieved. By subtracting a 0.8 % threshold from the total energy in every channel, a 2D spatial resolution of approximately 0.90 mm full width at half maximum (FWHM) can be obtained. Furthermore, a higher level of resolution, approximately 0.80 mm FWHM, was achieved using a CNN, with some alleviation of edge effects. With the proposed DOI method, a significant 1.36 mm FWHM average DOI resolution can be achieved. Additionally, it was found that polishing and black coating on the crystal surface yielded smaller edge effects compared to a rough surface with a black coating.Significance.The introduction of a threshold in COG method and a dual-ended readout scheme can lead to excellent spatial resolution for monolithic crystal detectors, which can help to develop PET systems with both high sensitivity and high spatial resolution.

Innovations in Small-Animal PET Instrumentation

Biomedical research has long relied on small-animal studies to elucidate disease process and develop new medical treatments. The introduction of in vivo functional imaging technology, such as PET, has allowed investigators to peer inside their subjects and follow disease progression longitudinally as well as improve understanding of normal biological processes. Recent developments in CRISPR, immuno-PET, and high-resolution in vivo imaging have only increased the importance of small-animal, or preclinical, PET imaging. Other drivers of preclinical PET innovation include new combinations of imaging technologies, such as PET/MR imaging, which require changes to PET hardware.

Timing, Energy, and 3-D Spatial Resolution of the BING PET Detector Module

We evaluated the 3D spatial, energy, and timing resolution of the Brain (or Breast)-Initiative Next-Generation (BING) PET detector. The BING detector is an array of 1-mm-thick slats of LYSO scintillator with lapped specular-reflective faces (15-mm by 52-mm) that are stacked together and oriented with their long-narrow edges normal to the imaging field of view. Interaction positions are determined from the signals of silicon-photomultiplier (SiPM) arrays placed on the entrance (top) and exit (bottom) faces. The SiPM arrays are offset to determine the slat of interaction (SOI) without requiring any optical light sharing between slats. Maximum likelihood 2D location within the SOI is determined using the sensor signals. Interaction time is determined with a modified first-optical-photon pickoff method. Performance of the BING detector was measured as a function of position using a sideways coincidence-collimated beam. Slats were accurately identified, with an effective tangential detector resolution of 1 mm. Average resolutions (and ranges) are: 0.96 mm (0.85 mm to 1.11 mm) for lateral (axial) detector resolution, 1.6 mm (1.0 mm to 2.1 mm) for depth resolution, 13.6% (12.7% to 16.0%) for energy resolution, and 317 ps (241 ps to 404 ps) for coincidence timing resolution. Initial spatial and timing resolution results demonstrated that the BING detector can be effective in a small field-of view (e.g., brain or breast) PET system.

Experimental validation of a rodent PET scanner prototype based on a single LYSO crystal tube

Improving sensitivity and spatial resolution in small animal Positron Emission Tomography imaging instrumentation constitutes one of the main goals of nuclear imaging research. These parameters are degraded by the presence of gaps between the detectors. The present manuscript experimentally validates our prototype of an edge-less pre-clinical PET system based on a single LYSO:Ce annulus with an inner diameter of 62 mm and 10 outer facets of 26 × 52 mm2. Scintillation light is read out by arrays of 8 × 8 SiPMs coupled to the facets, using a projection readout of the rows and columns signals. The readout provides accurate Depth of Interaction (DOI). We have implemented a calibration that mitigates the DOI-dependency of the transaxial and axial impact coordinates, and the energy photopeak gain. An energy resolution of 23.4 ± 1.8% was determined. Average spatial resolution of 1.4 ± 0.2 and 1.3 ± 0.4 mm FWHM were achieved for the radial and axial directions, respectively. We found a peak sensitivity of 3.8% at the system center, and a maximum NECR at 40.6 kcps for 0.27 mCi. The image quality was evaluated using reconstructed images of an array of sources and the NEMA image quality phantom was also studied.

Performance evaluation of side-by-side optically coupled monolithic LYSO crystals

Background

Significant interest has been recently shown for using monolithic scintillation crystals in molecular imaging systems, such as positron emission tomography (PET) scanners. Monolithic‐based PET scanners result in a lower cost and higher sensitivity, in contrast to systems based on the more conventional pixellated configuration. The monolithic design allows one to retrieve depth‐of‐interaction information of the impinging 511 keV photons without the need for additional hardware materials or complex positioning algorithms. However, the so‐called edge‐effect inherent to monolithic‐based approaches worsens the detector performance toward the crystal borders due to the truncation of the light distribution, thus decreasing positioning accuracy.

Purpose

The main goal of this work is to experimentally demonstrate the detector performance improvement when machine‐learning artificial neural‐network (NN) techniques are applied for positioning estimation in multiple monolithic scintillators optically coupled side‐by‐side.

Methods

In this work, we show the performance evaluation of two LYSO crystals of 33 × 25.4 × 10 mm³ optically coupled by means of a high refractive index adhesive compound (Meltmount, refractive index n = 1.70). A 12 × 12 silicon photomultiplier array has been used as photosensor. For comparison, the same detector configuration was tested for two additional coupling cases: (1) optical grease (n = 1.46) in between crystals, and (2) isolated crystals using black paint with an air gap at the interface (named standard configuration). Regarding 2D photon positioning (XY plane), we have tested two different methods: (1) a machine‐learning artificial NN algorithm and (2) a squared‐charge (SC) centroid technique.

Results

At the interface region of the detector, the SC method achieved spatial resolutions of 1.7 ± 0.3, 2.4 ± 0.3, and 2.6 ± 0.4 mm full‐width at half‐maximum (FWHM) for the Meltmount, grease, and standard configurations, respectively. These values improve to 1.0 ± 0.2, 1.2 ± 0.2, and 1.2 ± 0.3 mm FWHM when the NN algorithm was employed. Regarding energy performance, resolutions of 18 ± 2%, 20 ± 2%, and 23 ± 3% were obtained at the interface region of the detector for Meltmount, grease, and standard configurations, respectively.

Conclusions

The results suggest that optically coupling together scintillators with a high refractive index adhesive, in combination with an NN algorithm, reduces edge‐effects and makes it possible to build scanners with almost no gaps in between detectors.

Efficient Strike Artifact Reduction Based on 3D-Morphological Structure Operators from Filtered Back-Projection PET Images

Positron emission tomography (PET) can provide functional images and identify abnormal metabolic regions of the whole-body to effectively detect tumor presence and distribution. The filtered back-projection (FBP) algorithm is one of the most common images reconstruction methods. However, it will generate strike artifacts on the reconstructed image and affect the clinical diagnosis of lesions. Past studies have shown reduction in strike artifacts and improvement in quality of images by two-dimensional morphological structure operators (2D-MSO). The morphological structure method merely processes the noise distribution of 2D space and never considers the noise distribution of 3D space. This study was designed to develop three-dimensional-morphological structure operators (3D MSO) for nuclear medicine imaging and effectively eliminating strike artifacts without reducing image quality. A parallel operation was also used to calculate the minimum background standard deviation of the images for three-dimensional morphological structure operators with the optimal response curve (3D-MSO/ORC). As a result of Jaszczak phantom and rat verification, 3D-MSO/ORC showed better denoising performance and image quality than the 2D-MSO method. Thus, 3D MSO/ORC with a 3 × 3 × 3 mask can reduce noise efficiently and provide stability in FBP images.

Deep residual-convolutional neural networks for event positioning in a monolithic annular PET scanner

PET scanners based on monolithic pieces of scintillator can potentially produce superior performance characteristics (high spatial resolution and detection sensitivity, for example) compared to conventional PET scanners. Consequently, we initiated development of a preclinical PET system based on a single 7.2 cm long annulus of LYSO, called AnnPET. While this system could facilitate creation of high-quality images, its unique geometry results in optics that can complicate estimation of event positioning in the detector. To address this challenge, we evaluated deep-residual convolutional neural networks (DR-CNN) to estimate the three-dimensional position of annihilation photon interactions. Monte Carlo simulations of the AnnPET scanner were used to replicate the physics, including optics, of the scanner. It was determined that a ten-layer-DR-CNN was most suited to application with AnnPET. The errors between known event positions, and those estimated by this network and those calculated with the commonly used center-of-mass algorithm (COM) were used to assess performance. The mean absolute errors (MAE) for the ten-layer-DR-CNN-based event positions were 0.54 mm, 0.42 mm and 0.45 mm along thex(axial)-,y(transaxial)- andz- (depth-of-interaction) axes, respectively. For COM estimates, the MAEs were 1.22 mm, 1.04 mm and 2.79 mm in thex-,y- andz-directions, respectively. Reconstruction of the network-estimated data with the 3D-FBP algorithm (5 mm source offset) yielded spatial resolutions (full-width-at-half-maximum (FWHM)) of 0.8 mm (radial), 0.7 mm (tangential) and 0.71 mm (axial). Reconstruction of the COM-derived data yielded spatial resolutions (FWHM) of 1.15 mm (radial), 0.96 mm (tangential) and 1.14 mm (axial). These findings demonstrated that use of a ten-layer-DR-CNN with a PET scanner based on a monolithic annulus of scintillator has the potential to produce excellent performance compared to standard analytical methods.

A thick semi-monolithic scintillator detector for clinical PET scanners

Both monolithic and semi-monolithic scintillator positron emission tomography (PET) detectors can measure the depth of interaction with single-ended readout. Usually scintillators with a thickness of 10 mm or less are used since the position resolutions of the detectors degrade as the scintillator thickness increases. In this work, the performance of a 20 mm thick long rectangular semi-monolithic scintillator PET detector was measured by using both single-ended and dual-ended readouts with silicon photomultiplier (SiPM) arrays to provide a high detection efficiency. The semi-monolithic scintillator detector consists of nine lutetium-yttrium oxyorthosilicate slices measuring 1.37 × 51.2 × 20 mm³ with erythrocyte sedimentation rate foils of 0.065 mm thickness in between the slices. The SiPM array at each end of the scintillator detector consists of 16 × 4 SiPMs with a pixel size of 3.0 × 3.0 mm² and a pitch of 3.2 mm. The 64 signals of each SiPM array are processed by using the TOFPET2 application-specific integrated circuit individually. All but the edge slices can be clearly resolved for the detectors with both single-ended and dual-ended readouts. The single-ended readout detector provides an average full width at half maximum (FWHM) Y (continuous direction) position resolution of 2.43 mm, Z (depth direction) position resolution of 4.77 mm, energy resolution of 25.7% and timing resolution of 779 ps. The dual-ended readout detector significantly improves the Y and Z position resolutions, slightly improves the energy and timing resolution at the cost of two photodetectors required for one detector module and provides an average FWHM Y position resolution of 1.97 mm, Z position resolution of 2.60 mm, energy resolution of 21.7% and timing resolution of 718 ps. The energy and timing resolution of the semi-monolithic scintillator detector in this work are worse than those of the segmented scintillator array detector and need to be further improved. The semi-monolithic scintillator detector described in this work reduces costs as compared to the traditional segmented scintillator array detector and reduces the edge effect as compared to the monolithic scintillator detector.

Small animal PET: a review of what we have done and where we are going

Small animal research is an essential tool in studying both pharmaceutical biodistributions and disease progression over time. Furthermore, through the rapid development of in vivo imaging technology over the last few decades, small animal imaging (also referred to as preclinical imaging) has become a mainstay for all fields of biologic research and a center point for most preclinical cancer research. Preclinical imaging modalities include optical, MRI and MRS, microCT, small animal PET, ultrasound, and photoacoustic, each with their individual strengths. The strong points of small animal PET are its translatability to the clinic; its quantitative imaging capabilities; its whole-body imaging ability to dynamically trace functional/biochemical processes; its ability to provide useful images with only nano- to pico‑ molar concentrations of administered compounds; and its ability to study animals serially over time. This review paper gives an overview of the development and evolution of small animal PET imaging. It provides an overview of detector designs; system configurations; multimodality PET imaging systems; image reconstruction and analysis tools; and an overview of research and commercially available small animal PET systems. It concludes with a look toward developing technologies/methodologies that will further enhance the impact of small animal PET imaging on medical research in the future.

Development of depth encoding small animal PET detectors using dual-ended readout of pixelated scintillator arrays with SiPMs

PURPOSE

The performance of current small animal PET scanners is mainly limited by the detector performance and depth encoding detectors are required to develop PET scanner to simultaneously achieve high spatial resolution and high sensitivity. Among all depth encoding PET detector approaches, dual-ended readout detector has the advantage to achieve the highest depth of interaction (DOI) resolution and spatial resolution. Silicon photomultiplier (SiPM) is believed to be the photodetector of the future for PET detector due to its excellent properties as compared to the traditional photodetectors such as photomultiplier tube (PMT) and avalanche photodiode (APD). The purpose of this work is to develop high resolution depth encoding small animal PET detector using dual-ended readout of finely pixelated scintillator arrays with SiPMs.

METHODS

Four lutetium-yttrium oxyorthosilicate (LYSO) arrays with 11 × 11 crystals and 11.6 × 11.6 × 20 mm³ outside dimension were made using ESR, Toray and BaSO₄ reflectors. The LYSO arrays were read out with Hamamatsu 4 × 4 SiPM arrays from both ends. The SiPM array has a pixel size of 3 × 3 mm² , 0.2 mm gap in between the pixels and a total active area of 12.6 × 12.6 mm² . The flood histograms, DOI resolution, energy resolution and timing resolution of the four detector modules were measured and compared.

RESULTS

All crystals can be clearly resolved from the measured flood histograms of all four arrays. The BaSO₄ arrays provide the best and the ESR array provides the worst flood histograms. The DOI resolution obtained from the DOI profiles of the individual crystals of the four array is from 2.1 to 2.35 mm for events with E > 350 keV. The DOI ratio variation among crystals is bigger for the BaSO₄ arrays as compared to both the ESR and Toray arrays. The BaSO₄ arrays provide worse detector based DOI resolution. The photopeak amplitude of the Toray array had the maximum change with depth, it provides the worst energy resolution of 21.3%. The photopeak amplitude of the BaSO₄ array with 80 μm reflector almost doesn't change with depth, it provides the best energy resolution of 12.9%. A maximum timing shift of 1.37 ns to 1.61 ns among the corner and the center crystals in the four arrays was obtained due to the use of resistor network readout. A crystal based timing resolution of 0.68 ns to 0.83 ns and a detector based timing resolution of 1.26 ns to 1.45 ns were obtained for the four detector modules.

CONCLUSIONS

Four high resolution depth encoding small animal PET detectors were developed using dual-ended readout of pixelated scintillator arrays with SiPMs. The performance results show that those detectors can be used to build a small animal PET scanner to simultaneously achieve uniform high spatial resolution and high sensitivity.

Performance of a SiPM based semi-monolithic scintillator PET detector

A depth encoding PET detector module using semi-monolithic scintillation crystal single-ended readout by a SiPM array was built and its performance was measured. The semi-monolithic scintillator detector consists of 11 polished LYSO slices measuring 1  ×  11.6  ×  10 mm³. The slices are glued together with enhanced specular reflector (ESR) in between and outside of the slices. The bottom surface of the slices is coupled to a 4  ×  4 SiPM array with a 1 mm light guide and silicon grease between them. No reflector is used on the top surface and two sides of the slices to reduce the scintillation photon reflection. The signals of the 4  ×  4 SiPM array are grouped along rows and columns separately into eight signals. Four SiPM column signals are used to identify the slices according to the center of the gravity of the scintillation photon distribution in the pixelated direction. Four SiPM row signals are used to estimate the y (monolithic direction) and z (depth of interaction) positions according to the center of the gravity and the width of the scintillation photon distribution in the monolithic direction, respectively. The detector was measured with 1 mm sampling interval in both the y and z directions with electronic collimation by using a 0.25 mm diameter ²²Na point source and a 1  ×  1  ×  20 mm³ LYSO crystal detector. An average slice based energy resolution of 14.9% was obtained. All slices of 1 mm thick were clearly resolved and a detector with even thinner slices could be used. The y positions calculated with the center of gravity method are different for interactions happening at the same y, but different z positions due to depth dependent edge effects. The least-square minimization and the maximum likelihood positioning algorithms were developed and both methods improved the spatial resolution at the edges of the detector as compared with the center of gravity method. A mean absolute error (MAE) which is defined as the probability-weighted mean of the absolute value of the positioning error is used to evaluate the spatial resolution. An average MAE spatial resolution of ~1.15 mm was obtained in both y and z directions without rejection of the multiple scattering events. The average MAE spatial resolution was ~0.7 mm in both y and z directions after the multiple scattering events were rejected. The timing resolution of the detector is 575 ps. In the next step, long rectangle detector will be built to reduce edge effects and improve the spatial resolution of the semi-monolithic detector. Thick detector up to 20 mm will be explored and the positioning algorithms will be further optimized.

Depth of interaction determination in monolithic scintillator with double side SiPM readout

Background

Monolithic scintillators read out by arrays of photodetectors represent a promising solution to obtain high spatial resolution and the depth of interaction (DOI) of the annihilation photon. We have recently investigated a detector geometry composed of a monolithic scintillator readout on two sides by silicon photomultiplier (SiPM) arrays, and we have proposed two parameters for the DOI determination: the difference in the number of triggered SiPMs on the two sides of the detector and the difference in the maximum collected signal on a single SiPM on each side. This work is focused on the DOI calibration and on the determination of the capability of our detector. For the DOI calibration, we studied a method which can be implemented also in detectors mounted in a full PET scanner. We used a PET detector module composed of a monolithic 20 × 20 × 10 mm³ LYSO scintillator crystal coupled on two opposite faces to two arrays of SiPMs. On each side, the scintillator was coupled to 6 × 6 SiPMs. In this paper, the two parameters previously proposed for the DOI determination were calibrated with two different methods. The first used a lateral scan of the detector with a collimated 511 keV pencil beam at steps of 0.5 mm to study the detector DOI capability, while the second used the background radiation of the ¹⁷⁶Lu in the scintillator. The DOI determination capability was tested on different regions of the detector using each parameter and the combination of the two.

Results

With both parameters for the DOI determination, in the lateral scan, the bias between the mean reconstructed DOI and the real beam position was lower than 0.3 mm, and the DOI distribution had a standard deviation of about 1.5 mm. When using the calibration with the radioactivity of the LYSO, the mean bias increased of about 0.2 mm but with no degradation of the standard deviation of the DOI distribution.

Conclusions

The two parameters allow to achieve a DOI resolution comparable with the state of the art, giving a continuous information about the three-dimensional interaction position of the scintillation. These results were obtained by using simple estimators and a detector scalable to a whole PET system. The DOI calibration obtained using lutetium natural radioactivity gives results comparable to the other standard method but appears more readily applicable to detectors mounted in a full PET scanner.