Improving Depth, Energy and Timing Estimation in PET Detectors with Deconvolution and Maximum Likelihood Pulse Shape Discrimination

In a scintillation detector, the light generated in the scintillator by a gamma interaction is converted to photoelectrons by a photodetector and produces a time-dependent waveform, the shape of which depends on the scintillator properties and the photodetector response. Several depth-of-interaction (DOI) encoding strategies have been developed that manipulate the scintillator's temporal response along the crystal length and therefore require pulse shape discrimination techniques to differentiate waveform shapes. In this work, we demonstrate how maximum likelihood (ML) estimation methods can be applied to pulse shape discrimination to better estimate deposited energy, DOI and interaction time (for time-of-flight (TOF) PET) of a gamma ray in a scintillation detector. We developed likelihood models based on either the estimated detection times of individual photoelectrons or the number of photoelectrons in discrete time bins, and applied to two phosphor-coated crystals (LFS and LYSO) used in a previously developed TOF-DOI detector concept. Compared with conventional analytical methods, ML pulse shape discrimination improved DOI encoding by 27% for both crystals. Using the ML DOI estimate, we were able to counter depth-dependent changes in light collection inherent to long scintillator crystals and recover the energy resolution measured with fixed depth irradiation (~11.5% for both crystals). Lastly, we demonstrated how the Richardson-Lucy algorithm, an iterative, ML-based deconvolution technique, can be applied to the digitized waveforms to deconvolve the photodetector's single photoelectron response and produce waveforms with a faster rising edge. After deconvolution and applying DOI and time-walk corrections, we demonstrated a 13% improvement in coincidence timing resolution (from 290 to 254 ps) with the LFS crystal and an 8% improvement (323 to 297 ps) with the LYSO crystal.

Predicting the timing properties of phosphor-coated scintillators using Monte Carlo light transport simulation

Current developments in positron emission tomography focus on improving timing performance for scanners with time-of-flight (TOF) capability, and incorporating depth-of-interaction (DOI) information. Recent studies have shown that incorporating DOI correction in TOF detectors can improve timing resolution, and that DOI also becomes more important in long axial field-of-view scanners. We have previously reported the development of DOI-encoding detectors using phosphor-coated scintillation crystals; here we study the timing properties of those crystals to assess the feasibility of providing some level of DOI information without significantly degrading the timing performance. We used Monte Carlo simulations to provide a detailed understanding of light transport in phosphor-coated crystals which cannot be fully characterized experimentally. Our simulations used a custom reflectance model based on 3D crystal surface measurements. Lutetium oxyorthosilicate crystals were simulated with a phosphor coating in contact with the scintillator surfaces and an external diffuse reflector (teflon). Light output, energy resolution, and pulse shape showed excellent agreement with experimental data obtained on 3 × 3 × 10 mm³ crystals coupled to a photomultiplier tube. Scintillator intrinsic timing resolution was simulated with head-on and side-on configurations, confirming the trends observed experimentally. These results indicate that the model may be used to predict timing properties in phosphor-coated crystals and guide the coating for optimal DOI resolution/timing performance trade-off for a given crystal geometry. Simulation data suggested that a time stamp generated from early photoelectrons minimizes degradation of the timing resolution, thus making this method potentially more useful for TOF-DOI detectors than our initial experiments suggested. Finally, this approach could easily be extended to the study of timing properties in other scintillation crystals, with a range of treatments and materials attached to the surface.

Recent development in PET instrumentation

Positron emission tomography (PET) is used in the clinic and in vivo small animal research to study molecular processes associated with diseases such as cancer, heart disease, and neurological disorders, and to guide the discovery and development of new treatments. This paper reviews current challenges of advancing PET technology and some of newly developed PET detectors and systems. The paper focuses on four aspects of PET instrumentation: high photon detection sensitivity; improved spatial resolution; depth-of-interaction (DOI) resolution and time-of-flight (TOF). Improved system geometry, novel non-scintillator based detectors, and tapered scintillation crystal arrays are able to enhance the photon detection sensitivity of a PET system. Several challenges for achieving high resolution with standard scintillator-based PET detectors are discussed. Novel detectors with 3-D positioning capability have great potential to be deployed in PET for achieving spatial resolution better than 1 mm, such as cadmium-zinc-telluride (CZT) and position-sensitive avalanche photodiodes (PSAPDs). DOI capability enables a PET system to mitigate parallax error and achieve uniform spatial resolution across the field-of-view (FOV). Six common DOI designs, as well as advantages and limitations of each design, are discussed. The availability of fast scintillation crystals such as LaBr(3), and the silicon photomultiplier (SiPM) greatly advances TOF-PET development. Recent instrumentation and initial results of clinical trials are briefly presented. If successful, these technology advances, together with new probe molecules, will substantially enhance the molecular sensitivity of PET and thus increase its role in preclinical and clinical research as well as evaluating and managing disease in the clinic.