Giacomo Borghi

Sub-200 ps CRT in monolithic scintillator PET detectors using digital SiPM arrays and maximum likelihood interaction time estimation

Digital silicon photomultiplier (dSiPM) arrays have favorable characteristics for application in monolithic scintillator detectors for time-of-flight positron emission tomography (PET). To fully exploit these benefits, a maximum likelihood interaction time estimation (MLITE) method was developed to derive the time of interaction from the multiple time stamps obtained per scintillation event. MLITE was compared to several deterministic methods. Timing measurements were performed with monolithic scintillator detectors based on novel dSiPM arrays and LSO:Ce,0.2%Ca crystals of 16 × 16 × 10 mm(3), 16 × 16 × 20 mm(3), 24 × 24 × 10 mm(3), and 24 × 24 × 20 mm(3). The best coincidence resolving times (CRTs) for pairs of identical detectors were obtained with MLITE and measured 157 ps, 185 ps, 161 ps, and 184 ps full-width-at-half-maximum (FWHM), respectively. For comparison, a small reference detector, consisting of a 3 × 3 × 5 mm(3) LSO:Ce,0.2%Ca crystal coupled to a single pixel of a dSiPM array, was measured to have a CRT as low as 120 ps FWHM. The results of this work indicate that the influence of the optical transport of the scintillation photons on the timing performance of monolithic scintillator detectors can at least partially be corrected for by utilizing the information contained in the spatio-temporal distribution of the collection of time stamps registered per scintillation event.

A 32 mm × 32 mm × 22 mm monolithic LYSO:Ce detector with dual-sided digital photon counter readout for ultrahigh-performance TOF-PET and TOF-PET/MRI.

New applications for positron emission tomography (PET) and combined PET/magnetic resonance imaging (MRI) are currently emerging, for example in the fields of neurological, breast, and pediatric imaging. Such applications require improved image quality, reduced dose, shorter scanning times, and more precise quantification. This can be achieved by means of dedicated scanners based on ultrahigh-performance detectors, which should provide excellent spatial resolution, precise depth-of-interaction (DOI) estimation, outstanding time-of-flight (TOF) capability, and high detection efficiency. Here, we introduce such an ultrahigh-performance TOF/DOI PET detector, based on a 32 mm  ×  32 mm  ×  22 mm monolithic LYSO:Ce crystal. The 32 mm  ×  32 mm front and back faces of the crystal are coupled to a digital photon counter (DPC) array, in so-called dual-sided readout (DSR) configuration. The fully digital detector offers a spatial resolution of ~1.1 mm full width at half maximum (FWHM)/~1.2 mm mean absolute error, together with a DOI resolution of ~2.4 mm FWHM, an energy resolution of 10.2% FWHM, and a coincidence resolving time of 147 ps FWHM. The time resolution closely approaches the best results (135 ps FWHM) obtained to date with small crystals made from the same material coupled to the same DPC arrays, illustrating the excellent correction for optical and electronic transit time spreads that can be achieved in monolithic scintillators using maximum-likelihood techniques for estimating the time of interaction. The performance barely degrades for events with missing data (up to 6 out of 32 DPC dies missing), permitting the use of almost all events registered under realistic acquisition conditions. Moreover, the calibration procedures and computational methods used for position and time estimation follow recently made improvements that make them fast and practical, opening up realistic perspectives for using DSR monolithic scintillator detectors in TOF-PET and TOF-PET/MRI systems.

Towards monolithic scintillator based TOF-PET systems: practical methods for detector calibration and operation.

Gamma-ray detectors based on thick monolithic scintillator crystals can achieve spatial resolutions  <2 mm full-width-at-half-maximum (FWHM) and coincidence resolving times (CRTs) better than 200 ps FWHM. Moreover, they provide high sensitivity and depth-of-interaction (DOI) information. While these are excellent characteristics for clinical time-of-flight (TOF) positron emission tomography (PET), the application of monolithic scintillators has so far been hampered by the lengthy and complex procedures needed for position- and time-of-interaction estimation. Here, the algorithms previously developed in our group are revised to make the calibration and operation of a large number of monolithic scintillator detectors in a TOF-PET system practical. In particular, the k-nearest neighbor (k-NN) classification method for x,y-position estimation is accelerated with an algorithm that quickly preselects only the most useful reference events, reducing the computation time for position estimation by a factor of ~200 compared to the previously published k-NN 1D method. Also, the procedures for estimating the DOI and time of interaction are revised to enable full detector calibration by means of fan-beam or flood irradiations only. Moreover, a new technique is presented to allow the use of events in which some of the photosensor pixel values and/or timestamps are missing (e.g. due to dead time), so as to further increase system sensitivity. The accelerated methods were tested on a monolithic scintillator detector specifically developed for clinical PET applications, consisting of a 32 mm  ×  32 mm  ×  22 mm LYSO : Ce crystal coupled to a digital photon counter (DPC) array. This resulted in a spatial resolution of 1.7 mm FWHM, an average DOI resolution of 3.7 mm FWHM, and a CRT of 214 ps. Moreover, the possibility of using events missing the information of up to 16 out of 64 photosensor pixels is shown. This results in only a small deterioration of the detector performance.

Experimental Validation of an Efficient Fan-Beam Calibration Procedure for

Monolithic scintillator detectors can achieve excellent spatial resolution and coincidence resolving time. However, their practical use for positron emission tomography (PET) and other applications in the medical imaging field is still limited due to drawbacks of the different methods used to estimate the position of interaction. Common statistical methods for example require the collection of an extensive dataset of reference events with a narrow pencil beam aimed at a fine grid of reference positions. Such procedures are time consuming and not straightforwardly implemented in systems composed of many detectors. Here, we experimentally demonstrate for the first time a new calibration procedure for k-nearest neighbor ( k-NN) position estimation that utilizes reference data acquired with a fan beam. The procedure is tested on two detectors consisting of 16 mm ×16 mm ×10 mm and 16 mm ×16 mm ×20 mm monolithic, Ca-codoped LSO:Ce crystals and digital photon counter (DPC) arrays. For both detectors, the spatial resolution and the bias obtained with the new method are found to be practically the same as those obtained with the previously used method based on pencil-beam irradiation, while the calibration time is reduced by a factor of ~ 20. Specifically, a FWHM of ~ 1.1 mm and a FWTM of ~ 2.7 mm were obtained using the fan-beam method with the 10 mm crystal, whereas a FWHM of ~ 1.5 mm and a FWTM of ~ 6 mm were achieved with the 20 mm crystal. Using a fan beam made with a ~ 4.5 MBq 22Na point-source and a tungsten slit collimator with 0.5 mm aperture, the total measurement time needed to acquire the reference dataset was ~ 3 hours for the thinner crystal and ~ 2 hours for the thicker one.

k

Positron emission tomography (PET) is a medical imaging technique for non-invasive detection of glucose metabolism in human organs, such as the brain, heart, and lungs. A positron-emitting molecule, generally Fluorodeoxyglucose-18 (18F-FDG), is injected and subsequently identified by coincident detection of two gamma photons released during the annihilation of an electron with the positron lost by 1BF-FDG. Conventionally, coincident gamma photons are detected through scintillation by two sensors placed at 180° in a ring of identical sensors. In EndoTOFPET-US, two sensors are used, one placed in a semicircular ring outside the patient and the other placed on an endoscope. The endoscopic sensor is a miniaturized array of silicon photomultipliers (SiPMs) coupled to a mini-scintillator. To exploit the Fishburn-Seifert lower bound in timing resolution, a large number of photon timestamps must be generated, as is done in a multi-channel digital SiPM (MD-SiPM).

-Nearest Neighbor Position Estimation in Monolithic Scintillator Detectors

The goal of this simulation study is the performance evaluation and comparison of six potential designs for a time-of-flight PET scanner for pediatric patients of up to about 12 years of age. It is designed to have a high sensitivity and provide high-contrast and high-resolution images. The simulated pediatric PET is a full ring scanner, consisting of 32  ×  32 mm2 monolithic LYSO:Ce crystals coupled to digital silicon photomultiplier arrays. The six considered designs differ in axial lengths (27.2 cm, 54.4 cm and 102 cm) and crystal thicknesses (22 mm and 11 mm). The simulations are based on measured detector response data. We study two possible detector arrangements: 22 mm-thick crystals with dual-sided readout and 11 mm-thick crystals with back-sided readout. The six designs are simulated by means of the GEANT4 application for tomographic emission software, using the measured spatial, energy and time response of the monolithic scintillator detectors as input. The performance of the six designs is compared on the basis of four studies: (1) spatial resolution; (2) NEMA NU2-2012 sensitivity and scatter fraction (SF) tests; (3) non-prewhitening signal-to-noise ratio observer study; and (4) receiver operating characteristics analysis. Based on the results, two designs are identified as cost-effective solutions for fast and efficient imaging of children: one with 54.4 cm axial field-of-view (FOV) and 22 mm-thick crystals, and another one with 102 cm axial FOV and 11 cm-thick crystals. The first one has a higher center point sensitivity than the second one, but requires dual-sided readout. The second design has the advantage of allowing a whole-body scan in a single bed position acquisition. Both designs have the potential to provide an excellent spatial resolution (∼2 mm) and an ultra-high sensitivity (>100 cps [Formula: see text]).

11.4 A 67,392-SPAD PVTB-compensated multi-channel digital SiPM with 432 column-parallel 48ps 17b TDCs for endoscopic time-of-flight PET

Gamma-ray detectors based on bright monolithic scintillation crystals coupled to pixelated photodetectors are currently being considered for several applications in the medical imaging field. In a typical monolithic detector, both the light intensity and the time of arrival of the earliest scintillation photons can be recorded by each of the photosensor pixels every time a gamma interaction occurs. Generally, the time stamps are used to determine the gamma interaction time while the light intensities are used to estimate the 3D position of the interaction point. In this work we show that the spatio-temporal distribution of the time stamps also carries information on the location of the gamma interaction point and thus the time stamps can be used as explanatory variables for position estimation. We present a model for the spatial resolution obtainable when the interaction position is estimated using exclusively the time stamp of the first photon detected on each of the photosensor pixels. The model is shown to be in agreement with experimental measurements on a 16 mm  ×  16 mm  ×  10 mm LSO : Ce,0.2%Ca crystal coupled to a digital photon counter (DPC) array where a spatial resolution of 3 mm (root mean squared error) is obtained. Finally we discuss the effects of the main parameters such as scintillator rise and decay time, light output and photosensor single photon time resolution and pixel size.

Optimization of an ultralow-dose high-resolution pediatric PET scanner design based on monolithic scintillators with dual-sided digital SiPM readout: a simulation study

PSMR 2015: 4th Conference on PET/MR and SPECT/MR, La Biodola, Isola d’Elba, Italy, 17-21 May 2015

Time-based position estimation in monolithic scintillator detectors

Recently, a monolithic scintillator detector for time-of-flight (TOF)/depth-of-interaction (DOI) positron emission tomography (PET) was developed. It has a detector spatial resolution of ~1.7 mm full-width-at-half-maximum (FWHM), a coincidence resolving time (CRT) of ~215 ps FWHM, and ~4.7 mm FWHM DOI resolution. Here, we demonstrate, for the first time, the imaging performance of this detector in a 70 cm diameter PET geometry. We built a tomographic setup representative of a whole-body clinical scanner, comprising two coaxially rotating arms, each carrying a detector module, and a central, rotating phantom table. The fully automated setup sequentially acquires all possible lines of response (LORs) of a complete detector ring, using a step-and-shoot acquisition approach. The modules contained 2  ×  2 detectors, each detector consisting of a 32 mm  ×  32 mm  ×  22 mm LYSO crystal and a digital silicon photomultiplier (dSiPM) array. The system spatial resolution was assessed using a Na-22 point source at different radial distances in the field-of-view (FOV). Using 2D filtered back projection (2D FBP, non-TOF), tangential and radial spatial resolutions of ~2.9 mm FWHM were obtained at the center of the FOV. The use of DOI information resulted in almost uniform spatial resolution throughout the FOV up to a radial distance of 25 cm, where the radial and tangential resolution are ~3.3 mm FWHM and ~4.7 mm FWHM, respectively, whereas without DOI the resolution deteriorates to ~9 mm FWHM. Additional measurements were performed with a Na-22 filled Derenzo-like phantom at different locations within the FOV. Images reconstructed with a TOF maximum-likelihood expectation-maximization (TOF ML-EM) algorithm show that the system is able to clearly resolve 3 mm diameter hot rods up to 25 cm radial distance. The excellent and uniform spatial resolution, combined with an energy resolution of 10.2% FWHM and a CRT of ~212 ps FWHM, indicates a great potential for monolithic scintillators as practical high-performance detectors in TOF/DOI-PET systems.

Geometric optimization of an ultralow-dose high-resolution pediatric PET scanner based on monolithic scintillators with dSiPM readout.

We have recently built and characterized the performance of a monolithic scintillator detector based on a 32 mm  ×  32 mm  ×  22 mm LYSO:Ce crystal read out by digital silicon photomultiplier (dSiPM) arrays coupled to the crystal front and back surfaces in a dual-sided readout (DSR) configuration. The detector spatial resolution appeared to be markedly better than that of a detector consisting of the same crystal with conventional back-sided readout (BSR). Here, we aim to evaluate the influence of this difference in the detector spatial response on the quality of reconstructed images, so as to quantify the potential benefit of the DSR approach for high-resolution, whole-body time-of-flight (TOF) positron emission tomography (PET) applications. We perform Monte Carlo simulations of clinical PET systems based on BSR and DSR detectors, using the results of our detector characterization experiments to model the detector spatial responses. We subsequently quantify the improvement in image quality obtained with DSR compared to BSR, using clinically relevant metrics such as the contrast recovery coefficient (CRC) and the area under the localized receiver operating characteristic curve (ALROC). Finally, we compare the results with simulated rings of pixelated detectors with DOI capability. Our results show that the DSR detector produces significantly higher CRC and increased ALROC values than the BSR detector. The comparison with pixelated systems indicates that one would need to choose a crystal size of 3.2 mm with three DOI layers to match the performance of the BSR detector, while a pixel size of 1.3 mm with three DOI layers would be required to get on par with the DSR detector.