Performance comparison of dual-ended readout depth-encoding PET detectors based on BGO and LYSO crystals

The effects of various penalty parameter values in Q.Clear algorithm for rectal cancer detection on 18F-FDG images using a BGO-based PET/CT scanner: a phantom and clinical study

BACKGROUND

The Q.Clear algorithm is a fully convergent iterative image reconstruction technique. We hypothesize that different PET/CT scanners with distinct crystal properties will require different optimal settings for the Q.Clear algorithm. Many studies have investigated the improvement of the Q.Clear reconstruction algorithm on PET/CT scanner with LYSO crystals and SiPM detectors. We propose an optimum penalization factor (β) for the detection of rectal cancer and its metastases using a BGO-based detector PET/CT system which obtained via accurate and comprehensive phantom and clinical studies.

METHODS

18F-FDG PET-CT scans were acquired from NEMA phantom with lesion-to-background ratio (LBR) of 2:1, 4:1, 8:1, and 15 patients with rectal cancer. Clinical lesions were classified into two size groups. OSEM and Q.Clear (β value of 100-500) reconstruction was applied. In Q.Clear, background variability (BV), contrast recovery (CR), signal-to-noise ratio (SNR), SUVmax, and signal-to-background ratio (SBR) were evaluated and compared to OSEM.

RESULTS

OSEM had 11.5-18.6% higher BV than Q.Clear using β value of 500. Conversely, RC from OSEM to Q.Clear using β value of 500 decreased by 3.3-7.7% for a sphere with a diameter of 10 mm and 2.5-5.1% for a sphere with a diameter of 37 mm. Furthermore, the increment of contrast using a β value of 500 was 5.2-8.1% in the smallest spheres compared to OSEM. When the β value was increased from 100 to 500, the SNR increased by 49.1% and 30.8% in the smallest and largest spheres at LBR 2:1, respectively. At LBR of 8:1, the relative difference of SNR between β value of 100 and 500 was 43.7% and 44.0% in the smallest and largest spheres, respectively. In the clinical study, as β increased from 100 to 500, the SUVmax decreased by 47.7% in small and 31.1% in large lesions. OSEM demonstrated the least SUVmax, SBR, and contrast. The decrement of SBR and contrast using OSEM were 13.6% and 12.9% in small and 4.2% and 3.4%, respectively, in large lesions.

CONCLUSIONS

Implementing Q.Clear enhances quantitative accuracies through a fully convergent voxel-based image approach, employing a penalization factor. In the BGO-based scanner, the optimal β value for small lesions ranges from 200 for LBR 2:1 to 300 for LBR 8:1. For large lesions, the optimal β value is between 400 for LBR 2:1 and 500 for LBR 8:1. We recommended β value of 300 for small lesions and β value of 500 for large lesions in clinical study.

Phantom and clinical evaluation of Block Sequential Regularized Expectation Maximization (BSREM) reconstruction algorithm in 68Ga-PSMA PET-CT studies

In this study, we aimed to examine the effect of varying β-values in the block sequential regularized expectation maximization (BSREM) algorithm under differing lesion sizes to determine an optimal penalty factor for clinical application. The National Electrical Manufacturers Association phantom and 15 prostate cancer patients were injected with 68Ga-PSMA and scanned using a GE Discovery IQ PET/CT scanner. Images were reconstructed using ordered subset expectation maximization (OSEM) and BSREM with different β-values. Then, the background variability (BV), contrast recovery, signal-to-noise ratio, and lung residual error were measured from the phantom data, and the signal-to-background ratio (SBR) and contrast from the clinical data. The increment of BV using a β-value of 100 was 120.0%, and the decrement of BV using a β-value of 1000 was 40.5% compared to OSEM. As β decreased from 1000 to 100, the [Formula: see text] increased by 59.0% for a sphere with a diameter of 10 mm and 26.4% for a sphere with a diameter of 37 mm. Conversely, [Formula: see text] increased by 140.5% and 29.0% in the smallest and largest spheres, respectively. Furthermore, the Δ[Formula: see text] and Δ[Formula: see text] were - 41.1% and - 36.7%, respectively. In the clinical study, OSEM exhibited the lowest SBR and contrast. When the β-value was reduced from 500 to 100, the SBR and contrast increased by 69.7% and 71.8% in small and 35.6% and 33.0%, respectively, in large lesions. Moreover, the optimal β-value decreased as lesion size decreased. In conclusion, a β-value of 400 is optimal for small lesion reconstruction, while β-values of 600 and 500 are optimal for large lesions in phantom and clinical studies, respectively.

Potential of Depth-of-Interaction-Based Detection Time Correction in Cherenkov Emitter Crystals for TOF-PET

Cherenkov light can improve the timing resolution of Positron Emission Tomography (PET) radiation detectors, thanks to its prompt emission. Coincidence time resolutions (CTR) of ~30 ps were recently reported when using 3.2 mm-thick Cherenkov emitters. However, sufficient detection efficiency requires thicker crystals, causing the timing resolution to be degraded by the optical propagation inside the crystal. We report on depth-of-interaction (DOI) correction to mitigate the time-jitter due to the photon time spread in Cherenkov-based radiation detectors. We simulated the Cherenkov and scintillation light generation and propagation in 3 × 3 mm2 lead fluoride, lutetium oxyorthosilicate, bismuth germanate, thallium chloride, and thallium bromide. Crystal thicknesses varied from 9 to 18 mm with a 3-mm step. A DOI-based time correction showed a 2-to-2.5-fold reduction of the photon time spread across all materials and thicknesses. Results showed that highly refractive crystals, though producing more Cherenkov photons, were limited by an experimentally obtained high-cutoff wavelength and refractive index, restricting the propagation and extraction of Cherenkov photons mainly emitted at shorter wavelengths. Correcting the detection time using DOI information shows a high potential to mitigate the photon time spread. These simulations highlight the complexity of Cherenkov-based detectors and the competing factors in improving timing resolution.

Image Reconstruction Analysis for Positron Emission Tomography With Heterostructured Scintillators

The concept of structure engineering has been proposed for exploring the next generation of radiation detectors with improved performance. A TOF-PET geometry with heterostructured scintillators with a pixel size of 3.0 × 3.1 × 15 mm3 was simulated using Monte Carlo. The heterostructures consisted of alternating layers of BGO as a dense material with high stopping power and plastic (EJ232) as a fast light emitter. The detector time resolution was calculated as a function of the deposited and shared energy in both materials on an event-by-event basis. While sensitivity was reduced to 32% for 100-μm thick plastic layers and 52% for 50 μm, the coincidence time resolution (CTR) distribution improved to 204 ± 49 and 220 ± 41 ps, respectively, compared to 276 ps that we considered for bulk BGO. The complex distribution of timing resolutions was accounted for in the reconstruction. We divided the events into three groups based on their CTR and modeled them with different Gaussian TOF kernels. On an NEMA IQ phantom, the heterostructures had better contrast recovery in early iterations. On the other hand, BGO achieved a better contrast-to-noise ratio (CNR) after the 15th iteration due to the higher sensitivity. The developed simulation and reconstruction methods constitute new tools for evaluating different detector designs with complex time responses.

A personal acquisition time regimen of 68Ga-DOTATATE total-body PET/CT in patients with neuroendocrine tumor (NET): a feasibility study

BACKGROUND

The injection activity of tracer, acquisition time, patient-specific photon attenuation, and large body mass, can influence on image quality. Fixed acquisition time and body mass related injection activity in clinical practice results in a large difference in image quality. Thus, this study proposes a patient-specific acquisition time regimen of ⁶⁸ Ga-DOTATATE total-body positron emission tomography-computed tomography (PET/CT) to counteract the influence of body mass (BM, kg) on image quality, and acquire an acceptable and constant image of patients with neuroendocrine tumors (NETs).

METHODS

The development cohort consisting of 19 consecutive patients with full activity (88.7-204.9 MBq, 2.0 ± 0.1 MBq/kg) was to establish the acquisition time regimen. The liver SNR (signal-to-noise ratio, SNR_L) was normalized (SNR_norm) by the product of injected activity (MBq) and acquisition time (min). Fitting of SNR_norm against body mass (BM, kg) in linear correlation was performed. Subjective assessment of image quality was performed using a 5-point Likert scale to determine the acceptable threshold of SNR_L, and an optimized acquisition regimen based on BM was proposed, and validated its feasibility through the validation cohort of 57 consecutive NET patients with half activity (66.9 ± 11.3 MBq, 1.0 ± 0.1 MBq/kg) and a fixed acquisition time regimen.

RESULTS

The linear correlation (R² = 0.63) between SNR_norm and BM (kg) was SNR_norm = -0.01*BM + 1.50. The threshold SNR_L of acceptable image quality was 11.2. The patient-specific variable acquisition time regimen was determined as: t (min) = 125.4/(injective activity)*(-0.01*BM + 1.50)². Based on that proposed regimen, the average acquisition time for acceptable image quality in the validation cohort was 2.99 ± 0.91 min, ranging from 2.18 to 6.35 min, which was reduced by 36.50% ~ 78.20% compared with the fixed acquisition time of 10 min. Subjective evaluation showed that acceptable image quality could be obtained at 3.00 min in the validation group, with an average subjective score of 3.44 ± 0.53 (kappa = 0.97, 95% CI: 0.96 ~ 0.98). Bland-Altman analysis revealed good agreement between the proposed regimen and the fixed acquisition time cohort.

CONCLUSION

A patient-specific acquisition time regimen was proposed in NET patients in development cohort and validated its feasibility in patients with NETs in validation cohort by ⁶⁸ Ga-DOTATATE total-body PET/CT imaging. Based on the proposed regimen, the homogenous image quality with optimal acquisition time was available independent of body mass.

Performance of Dual-Ended Readout PET Detectors Based on BGO Arrays and BaSO₄ Reflector

In this paper, the performance of two dual-ended readout PET detectors based on 15 × 15 BGO arrays were compared. The crystal elements of one BGO array have polished lateral surfaces, while the crystal elements of the other BGO array have unpolished lateral surfaces. The two ends of the BGO elements are polished. The two BGO arrays both have a pitch size of 1.6 mm and thickness of 20 mm, and BaSO4 with a thickness of 80 μm was used as the reflector. Hamamatsu S14161-0305-08 SiPM arrays were used as photodetectors. All the measurements were performed at a bias voltage of 41.0 V and a temperature of 23.5 °C. The flood histograms show that all the crystal elements in the two BGO arrays were clearly resolved. The detector based on the BGO array with polished lateral surfaces provides an energy resolution of 16.9 ± 1.3%, timing resolution of 3.2 ± 0.2 ns, and DOI resolution of 18.4 ± 2.2 mm. In comparison, the detector based on the BGO array with unpolished lateral surfaces provides an energy resolution of 17.7 ± 2.0%, timing resolution of 3.5 ± 0.3 ns, and DOI resolution of 3.2 ± 0.2 mm.

Design and modeling of a high resolution and high sensitivity PET brain scanner with double-ended readout

In the wake of recent advancements in scintillator, photodetector, and low-noise fast electronics technologies, as well as in fast reconstruction software, positron emission tomography (PET) scanners have seen considerable improvements in spatial resolution, time resolution, and absolute sensitivity. To continue this trend, we present a helmet type PET brain scanner design that combines high solid angle coverage and double-ended readout of 30 mm-thick scintillator crystals to achieve excellent absolute sensitivity, depth of interaction resolution, and time resolution. This scanner comprises 598 detector arrays, each with 8×8 Lu_1.8Y_0.2SiO₅:Ce (LYSO:Ce) crystals with dimensions 3.005×3.005×30 mm³one-to-one coupled on either end to silicon photomultipliers (SiPMs). Our Monte Carlo simulations based in the platform Geant4 predict that this scanner would attain an absolute sensitivity to a 35 cm line source placed at the center of the radial field of view of (17.1±0.1) %, a depth of interaction resolution of (3.99±0.05) mm, and a coincidence time resolution of (198±5) ps. Our simulations also predict radial, tangential, and axial spatial resolutions at the center of the field of view of 3.3 mm, 3.1 mm, and 3.3 mm, respectively. As this set of simultaneous parameters compares favorably to today's most advanced clinical PET scanners and other proposed designs, this scanner has a good chance of becoming a preferred tool for high quality brain imaging.

Reducing Calibration Time in PET Systems Based on Monolithic Crystals

In the past years, the gamma-ray detector designs based on the monolithic crystals have demonstrated to be excellent candidates for the design of high-performance PET systems. The monolithic crystals allow to achieve the intrinsic detector resolutions well below state-of-the-art; to increase packing fraction thus, increasing the system sensitivity; and to improve lesion detectability at the edges of the scanner field of view (FOV) because of their intrinsic depth of interaction (DOI) capabilities. The bottleneck to translate to the clinical PET systems based on a large number of monolithic detectors is eventually the requirement of mechanically complex and time-consuming calibration processes. To mitigate this drawback, several methods have been already proposed, such as using non-physically collimated radioactive sources or implementing the neuronal networks (NN) algorithms trained with simulated data. In this work, we aimed to simplify and fasten a calibration process of the monolithic based systems. The Normal procedure consists of individually acquiring a 11 × 11 ²²Na source array for all the detectors composing the PET system and obtaining the calibration map for each module using a method based on the Voronoi diagrams. Two reducing time methodologies are presented: (i) TEST1, where the calibration map of one detector is estimated and shared among all others, and (ii) TEST2, where the calibration map is slightly modified for each module as a function of their detector uniformity map. The experimental data from a dedicated prostate PET system was used to compare the standard calibration procedure with both the proposed methods. A greater similarity was exhibited between the TEST2 methodology and the Normal procedure; obtaining spatial resolution variances within 0.1 mm error bars and count rate deviations as small as 0.2%. Moreover, the negligible reconstructed image differences (13% deviation at most in the contrast-to-noise ratio) and almost identical contrast values were reported. Therefore, this proposed method allows us to calibrate the PET systems based on the monolithic crystals reducing the calibration time by approximately 80% compared with the Normal procedure.

Performance evaluation of dual-ended readout PET detectors based on BGO arrays with different reflector arrangements

OBJECTIVE

Dual-ended readout depth-encoding detectors based on bismuth germanate (BGO) scintillation crystal arrays are good candidates for high-sensitivity small animal positron emission tomography used for very-low-dose imaging. In this paper, the performance of three dual-ended readout detectors based on 15 × 15 BGO arrays with three different reflector arrangements and 8 × 8 silicon photomultiplier arrays were evaluated and compared.

APPROACH

The three BGO arrays, denoted wo-ILG (without internal light guide), wp-ILG (with partial internal light guide), and wf-ILG (with full internal light guide), share a pitch size of 1.6 mm and thickness of 20 mm. Toray E60 with a thickness of 50μm was used as inter-crystal reflector. All reflector lengths in the wo-ILG and wf-ILG BGO arrays were 20 and 18 mm, respectively; the reflectors in the wp-ILG BGO array were 18 mm at the central region of the array and 20 mm at the edge. By using 18 mm reflectors, part of the crystals in the wp-ILG and wf-ILG BGO arrays worked as internal light guides.

MAIN RESULTS

The results showed that the detector based on the wo-ILG BGO array provided the best flood histogram. The energy, timing and DOI resolutions of the three detectors were similar. The energy resolutions full width at half maximum (FWHM value) based on the wo-ILG, wp-ILG and wf-ILG BGO arrays were 27.2 ± 3.9%, 28.7 ± 4.6%, and 29.5 ± 4.7%, respectively. The timing resolutions (FWHM value) were 4.7 ± 0.5 ns, 4.9 ± 0.5 ns, and 5.0 ± 0.6 ns, respectively. The DOI resolution (FWHM value) were 3.0 ± 0.2 mm, 2.9 ± 0.2 mm, and 3.0 ± 0.2 mm, respectively. Over all, the wo-ILG detector provided the best performance.

Development of dual-ended depth-of-interaction detectors using laser-induced crystals for small animal PET systems

Sensitivity and spatial resolution of positron emission tomography (PET) scanners can be improved by using thicker scintillation crystals with depth-of-interaction (DOI) encoding. Subsurface laser engraving (SSLE) can be used to segment crystals of a scintillation detector in order to fabricate a DOI detector. We previously applied SSLE to crystal bars of 3 × 3 × 20 mm³and 1.5 × 1.5 × 20 mm³and developed two dual-ended detectors with DOI segments of 3 mm and 1.5 mm, respectively. To further improve the DOI resolution, our SSLE detector design can be used with smaller pitch crystal bars, making them excellent detector candidates for small animal PET scanners with submillimetre resolution. In the present study, three small crystal bars of 1 × 1 × 20 mm³, 2 × 1 × 20 mm³, and 2 × 1 × 40 mm³were laser engraved to 12, 20 and 40 segments, respectively, by applying SSLE in their height directions. The segmented crystal bars were characterised in three prototype detector arrangements. First, the 1 × 1 × 20 mm³crystal bars were characterised in an 8 × 8 crystal array designed for DOI encoding along crystal height in a conventional small animal PET design. Second, a 4 × 8 crystal array of 2 × 1 × 20 mm³crystal bars was characterised for using the DOI information for crystal interaction positioning along the axial axis of a small animal PET scanner. Finally, the third part of the study was performed on a single 2 × 1 × 40 mm³crystal bar with 40 segments to investigate the feasibility of DOI estimation in longer crystals for application in a system with extended axial length. We evaluated the capability of segment identification and energy resolution of theses detectors. The 3D position maps of the detectors were obtained using the Anger-type calculation and the crystal identification performance was evaluated for each detector. Clear segment separation was obtained for the crystal arrays with 12 (segment pitch of 1.67 mm) and 20 (segment pitch of 1 mm) segments. Mean energy resolutions of 8.8% ± 0.4% and 9.6% ± 0.8% at 511 keV were obtained for the segments in the central regions of the 8 × 8 array with 12 segments and the 4 × 8 array with 20 segments, respectively. Clear segment identification was found to be difficult for the detector with 40 segments, especially for the segments at the middle of the crystal. Energy and interaction positioning characterisation results suggest that both prototype detectors with 12 and 20 segments are well suited for small animal PET scanners with high spatial resolution.

Optimization of scintillator-reflector optical interfaces for the LUT Davis model

PURPOSE

Designing and optimizing scintillator-based gamma detector using Monte Carlo simulation is of great importance in nuclear medicine and high energy physics. In scintillation detectors, understanding the light transport in the scintillator and the light collection by the photodetector plays a crucial role in achieving high performance. Thus, accurately modeling them is critical.

METHODS

In previous works, we developed a model to compute crystal reflectance from the crystal 3D surface measurement and store it in look-up tables to be used in the Monte Carlo simulation software GATE. The relative light output comparison showed excellent agreement between simulations and experiments for both polished and rough surfaces in several configurations, that is, without and with reflector. However, when comparing them at the irradiation depth closest to the photodetector face, rough crystals with a reflector overestimated the predicted light output. Investigating the cause of this overestimation, we optimized the LUT algorithm to improve the reflectance computation accuracy, especially for rough surfaces. However, optical Monte Carlo simulations carried out with these newly generated LUTs still overestimate the light output. Based on previous observations, one probable cause is the erroneous assumption of perfect couplings between the reflector and crystal and between the crystal and photodetector, which likely results in an important overestimation of the light output compared to experimental values. In practice, several factors could degrade it. Here, we investigated possible suboptimal optical experimental configurations that could lead to a degraded light collection when using Teflon or ESR reflectors coupled to the crystal with air or grease. We generated look-up tables with a mixture of air and grease and showed the effect of three possible sources of light loss: the presence of a small gap between the crystal and the reflector edges close to the photodetector face, the infiltration of grease in the crystal-reflector coupling, and the presence of inhomogeneities in the photodetector-crystal interface.

RESULTS

The strongest effect is linked to the presence of a small gap of grease between the edges of the reflector material and the crystal (light loss of 10%-12% for 0.2 mm gap). The optical grease infiltrating the crystal-reflector air coupling decreases the light output, depending on the infiltration's extent and the amount of grease infiltrated. Five percent of air in the crystal-photodetector coupling can cause a light output decrease of 2% to 4%. The individual and combined effect of these advanced models can explain the discrepancy of the relative light output obtained with ESR in simulations and experiments. With Teflon, the study indicates that the light output loss strongly depends on the reflectance deterioration caused by grease absorption.

CONCLUSIONS

Our results indicate that when studying scintillation detector performance with different finishes, performing simulations in ideal coupling conditions can lead to light output overestimation. To perform an accurate light output comparison and ultimately have a reliable detector performance estimation, all potential sources of practical limitations must be carefully considered. To broadly enable high-fidelity modeling, we developed an interface for users to compute their own LUTs, using their surface, scintillator, and reflector characteristics.

A high resolution and high detection efficiency depth-encoding detector for brain positron emission tomography based on a 0.75 mm pitch scintillator array

The quantitative accuracy and precision of brain positron emission tomography (PET) studies can be considerably improved using dedicated brain PET scanners with a uniform high resolution and a high sensitivity across the brain volume. One approach to building such a system is to construct the PET scanner using depth-of-interaction (DOI) encoding detectors with finely segmented and thick crystal arrays. In this paper, the performance of a DOI PET detector based on two 16 × 16 arrays of 2 × 2 mm2 SiPMs coupled to both ends of a 44 × 44 array of 0.69 × 0.69 × 30 mm3 polished LYSO crystals was evaluated at different temperatures (-9°C, 0°C, 10°C, and 20°C) for brain PET applications. The pitch size of the LYSO array is 0.75 mm. The flood histograms show that all the crystal elements in the LYSO array can be resolved except some edge crystals, due to the limited light sharing. The average energy resolution, average DOI resolution, and average timing resolution across crystal elements are 21.1 ± 3.0%, 3.47 ± 0.17 mm, and 1.38 ± 0.09 ns, respectively, which were obtained at a bias voltage of 56.5 V and a temperature of 0°C.