A finely segmented semi-monolithic detector tailored for high-resolution PET.

BACKGROUND: Preclinical research and organ-dedicated applications use and require high (spatial-)resolution positron emission tomography (PET) detectors to visualize small structures (early) and understand biological processes at a finer level of detail. Researchers seeking to improve detector and image spatial resolution have explored various detector designs. Current commercial high-resolution systems often employ finely pixelated or monolithic scintillators, each with its limitations. PURPOSE: We present a semi-monolithic detector, tailored for high-resolution PET applications with a spatial resolution in the range of 1 mm or better, merging concepts of monolithic and pixelated crystals. The detector features LYSO slabs measuring (24 × 10 × 1) mm3, coupled to a 12 × 12 readout channel photosensor with 4 mm pitch. The slabs are grouped in two arrays of 44 slabs each to achieve a higher optical photon density despite the fine segmentation. METHODS: We employ a fan beam collimator for fast calibration to train machine-learning-based positioning models for all three dimensions, including slab identification and depth-of-interaction (DOI), utilizing gradient tree boosting (GTB). The data for all dimensions was acquired in less than 2 h. Energy calculation was based on a position-dependent energy calibration. Using an analytical timing calibration, time skews were corrected for coincidence timing resolution (CTR) estimation. RESULTS: Leveraging machine-learning-based calibration in all three dimensions, we achieved high detector spatial resolution: down to 1.18 mm full width at half maximum (FWHM) detector spatial resolution and 0.75 mm mean absolute error (MAE) in the planar-monolithic direction, and 2.14 mm FWHM and 1.03 mm MAE for DOI at an energy window of (435-585) keV. Correct slab interaction identification in planar-segmented direction exceeded 80%, alongside an energy resolution of 12.7% and a CTR of 450 ps FWHM. CONCLUSIONS: The introduced finely segmented, high-resolution slab detector demonstrates appealing performance characteristics suitable for high-resolution PET applications. The current benchtop-based detector calibration routine allows these detectors to be used in PET systems.

Improving the Timing Resolution of Positron Emission Tomography Detectors Using Boosted Learning-A Residual Physics Approach.

Artificial intelligence (AI) is entering medical imaging, mainly enhancing image reconstruction. Nevertheless, improvements throughout the entire processing, from signal detection to computation, potentially offer significant benefits. This work presents a novel and versatile approach to detector optimization using machine learning (ML) and residual physics. We apply the concept to positron emission tomography (PET), intending to improve the coincidence time resolution (CTR). PET visualizes metabolic processes in the body by detecting photons with scintillation detectors. Improved CTR performance offers the advantage of reducing radioactive dose exposure for patients. Modern PET detectors with sophisticated concepts and read-out topologies represent complex physical and electronic systems requiring dedicated calibration techniques. Traditional methods primarily depend on analytical formulations successfully describing the main detector characteristics. However, when accounting for higher-order effects, additional complexities arise matching theoretical models to experimental reality. Our work addresses this challenge by combining traditional calibration with AI and residual physics, presenting a highly promising approach. We present a residual physics-based strategy using gradient tree boosting and physics-guided data generation. The explainable AI framework SHapley Additive exPlanations (SHAPs) was used to identify known physical effects with learned patterns. In addition, the models were tested against basic physical laws. We were able to improve the CTR significantly (more than 20%) for clinically relevant detectors of 19 mm height, reaching CTRs of 185 ps (450-550 keV).

Near-field coded-mask technique and its potential for proton therapy monitoring.

Objective.Prompt-gamma imaging encompasses several approaches to the online monitoring of the beam range or deposited dose distribution in proton therapy. We test one of the imaging techniques - a coded mask approach - both experimentally and via simulations.Approach.Two imaging setups have been investigated experimentally. Each of them comprised a structured tungsten collimator in the form of a modified uniformly redundant array mask and a LYSO:Ce scintillation detector of fine granularity. The setups differed in detector dimensions and operation mode (1D or 2D imaging). A series of measurements with radioactive sources have been conducted, testing the performance of the setups for near-field gamma imaging. Additionally, Monte Carlo simulations of a larger setup of the same type were conducted, investigating its performance with a realistic gamma source distribution occurring during proton therapy.Main results.The images of point-like sources reconstructed from two small-scale prototypes' data using the maximum-likelihood expectation maximisation algorithm constitute the experimental proof of principle for the near-field coded-mask imaging modality, both in the 1D and the 2D mode. Their precision allowed us to calibrate out certain systematic offsets appearing due to the limited alignment accuracy of setup elements. The simulation of the full-scale setup yielded a mean distal falloff retrieval precision of 0.72 mm in the studies for beam energy range 89.5-107.9 MeV and with 1 × 108protons (a typical number for distal spots). The implemented algorithm of image reconstruction is relatively fast-a typical procedure needs several seconds.Significance.Coded-mask imaging appears a valid option for proton therapy monitoring. The results of simulations let us conclude that the proposed full-scale setup is competitive with the knife-edge-shaped and the multi-parallel slit cameras investigated by other groups.

Analysis of a convex time skew calibration for light sharing-based PET detectors.

Objective.Positron emission tomography (PET) detectors providing attractive coincidence time resolutions (CTRs) offer time-of-flight information, resulting in an improved signal-to-noise ratio of the PET image. In applications with photosensor arrays that employ timestampers for individual channels, timestamps typically are not time synchronized, introducing time skews due to different signal pathways. The scintillator topology and transportation of the scintillation light might provoke further skews. If not accounted for these effects, the achievable CTR deteriorates. We studied a convex timing calibration based on a matrix equation. In this work, we extended the calibration concept to arbitrary structures targeting different aspects of the time skews and focusing on optimizing the CTR performance for detector characterization. The radiation source distribution, the stability of the estimations, and the energy dependence of calibration data are subject to the analysis.Approach.A coincidence setup, equipped with a semi-monolithic detector comprising 8 LYSO slabs, each 3.9 mm × 31.9 mm × 19.0 mm, and a one-to-one coupled detector with 8 × 8 LYSO segments of 3.9 mm × 3.9 mm × 19.0 mm volume is used. Both scintillators utilize a dSiPM (DPC3200-22-44, Philips Digital Photon Counting) operated in first photon trigger. The calibration was also conducted with solely one-to-one coupled detectors and extrapolated for a slab-only setup.Main results.All analyzed hyperparameters show a strong influence on the calibration. Using multiple radiation positions improved the skew estimation. The statistical significance of the calibration dataset and the utilized energy window was of great importance. Compared to a one-to-one coupled detector pair achieving CTRs of 224 ps the slab detector configuration reached CTRs down to 222 ps, demonstrating that slabs can compete with a clinically used segmented detector design.Significance.This is the first work that systematically studies the influence of hyperparameters on skew estimation and proposes an extension to arbitrary calibration structures (e.g. scintillator volumes) of a known calibration technique.

A semi-monolithic detector providing intrinsic DOI-encoding and sub-200 ps CRT TOF-capabilities for clinical PET applications.

BACKGROUND: Current clinical positron emission tomography (PET) systems utilize detectors where the scintillator typically contains single elements of 3-6-mm width and about 20-mm height. While providing good time-of-flight performance, this design limits the spatial resolution and causes radial astigmatism as the depth-of-interaction (DOI) remains unknown. PURPOSE: We propose an alternative, aiming to combine the advantages of current detectors with the DOI capabilities shown for monolithic concepts, based on semi-monolithic scintillators (slabs). Here, the optical photons spread along one dimension enabling DOI-encoding with a still small readout area beneficial for timing performance. METHODS: An array of eight monolithic LYSO slabs of dimensions 3.9 × 32 × 19 mm3 was read out by a 64-channel photosensor containing digital SiPMs (DPC3200-22-44, Philips Digital Photon Counting). The position estimation in the detector's monolithic and DOI direction was based on a calibration with a fan beam collimator and the machine learning technique gradient tree boosting (GTB). RESULTS: We achieved a positioning performance in terms of mean absolute error (MAE) of 1.44 mm for the monolithic direction and 2.12 mm for DOI considering a wide energy window of 300-700 keV. The energy resolution was determined to be 11.3%, applying a positional-dependent energy calibration. We established both an analytical and machine-learning-based timing calibration approach and applied them for a first-photon trigger. The analytical timing calibration corrects for electronic and optical time skews leading to 240 ps coincidence resolving time (CRT) for a pair of slab-detectors. The CRT was significantly improved by utilizing GTB to predict the time difference based on specific training data and applied on top of the analytical calibration. We achieved 209 ps for the wide energy window and 198 ps for a narrow selection around the photopeak (411-561 keV). To maintain the detector's sensitivity, no filters were applied to the data during processing. CONCLUSION: Overall, the semi-monolithic detector provides attractive performance characteristics. Especially, a good CRT can be achieved while introducing DOI capabilities to the detector, making the concept suitable for clinical PET scanners.

High throughput software-based gradient tree boosting positioning for PET systems.

The supervised machine learning technique Gradient Tree Boosting (GTB) has shown good accuracy for position estimation of gamma interaction in PET crystals for bench-top experiments while its computational requirements can easily be adjusted. Transitioning to preclinical and clinical applications requires near real-time processing in the scale of full PET systems. In this work, a high throughput GTB-based singles positioning C++ implementation is proposed and a series of optimizations are evaluated regarding their effect on the achievable processing throughput. Moreover, the crucial feature and parameter selection for GTB is investigated for the segmented detectors of the Hyperion IIDPET insert with two main models and a range of GTB hyperparameters. The proposed framework achieves singles positioning throughputs of more than 9.5 GB/s for smaller models and of 240 MB/s for more complex models on a recent Intel Skylake server. Detailed throughput analysis reveals the key performance limiting factors, and an empirical throughput model is derived to guide the GTB model selection process and scanner design decisions. The throughput model allows for throughput estimations with a mean absolute error (MAE) of 175.78 MB/s.

Evaluation of the radiofrequency performance of a wide-bore 1.5 T positron emission tomography/magnetic resonance imaging body coil for radiotherapy planning.

BACKGROUND AND PURPOSE: The restricted bore diameter of current simultaneous positron emission tomography/magnetic resonance imaging (PET/MRI) systems can be an impediment to achieving similar patient positioning during PET/MRI planning and radiotherapy. Our goal was to evaluate the B1 transmit (B1 +) uniformity, B1 + efficiency, and specific absorption rate (SAR) of a novel radiofrequency (RF) body coil design, in which RF shielded PET detectors were integrated with the specific aim of enabling a wide-bore PET/MRI system. MATERIALS AND METHODS: We designed and constructed a wide-bore PET/MRI RF body coil to be integrated with a clinical MRI system. To increase its inner bore diameter, the PET detectors were positioned between the conductors and the RF shield of the RF body coil. Simulations and experiments with phantoms and human volunteers were performed to compare the B1 + uniformity, B1 + efficiency, and SAR between our design and the clinical body coil. RESULTS: In the simulations, our design achieved nearly the same B1 + field uniformity as the clinical body coil and an almost identical SAR distribution. The uniformity findings were confirmed by the physical experiments. The B1 + efficiency was 38% lower compared to the clinical body coil. CONCLUSIONS: To achieve wide-bore PET/MRI, it is possible to integrate shielding for PET detectors between the body coil conductors and the RF shield without compromising MRI performance. Reduced B1 + efficiency may be compensated by adding a second RF amplifier. This finding may facilitate the application of simultaneous whole-body PET/MRI in radiotherapy planning.

RF shielding materials for highly-integrated PET/MRI systems.

Integration of magnetic resonance imaging (MRI) and positron emission tomography (PET) into a simultaneous device calls for adaptations of the radio frequency (RF) shielding concept. Conventional PET module housings fully encase the entire PET detector to reduce mutual interference. Excluding passive components, i.e. scintillators, from the housings, offers integration advantages, e.g. by reducing the overall housing volume or utilizing bigger scintillators. However, locating the scintillator outside of the RF shielding requires an optically transparent RF shielding interface between the scintillators and the photodetector to close the aperture. Therefore, a careful evaluation and selection of RF materials is essential to ensure an excellent PET/MRI system performance. To this end, we examined 10 materials (coated glasses, coated foils, meshes). The shielding effectiveness (SE) was evaluated at 100 and 300 MHz. PET performance was tested for single event registration and coincident events by integrating the material into the PET detector stack between the digital silicon photomultiplier photodetector array and one-to-one coupled scintillator. We determined photon attenuation (PA), energy resolution (dE/E), and coincidence resolving time (CRT) and compared to reference measurements for each material group. MRI compatibility was assessed by analyzing the material influence on the main magnetic field (B0) homogeneity. The coated glasses and foils exhibited SEs of up to 25 dB at 300 MHz. Both had a PA < 23% with dE/E and CRT comparable to the reference measurements, and no measurable impact on the B0field was registered. The meshes exhibited higher SEs up to 56 dB, but also a PA > 58% with a higher impact on dE/E and CRT. Only one mesh affected B0homogeneity. Overall, we recommend the coated foil HS 9400 for integration concepts as it exhibited a good performance with SE = 25 dB, PA = 22%, resulting in a PET performance of dE/E = 12% and CRT = 274 ps.

Evaluation of the PETsys TOFPET2 ASIC in multi-channel coincidence experiments.

BACKGROUND: Aiming to measure the difference in arrival times of two coincident γ-photons with an accuracy in the order of 200ps, time-of-flight positron emission tomography systems commonly employ silicon photomultipliers (SiPMs) and high-resolution digitization electronics, application specific integrated circuits (ASICs). This work evaluates the performance of the TOFPET2 ASIC, released by PETsys Electronics S.A. in 2017, dependent on its configuration parameters in multi-channel coincidence measurements. METHODS: SiPM arrays fabricated by different vendors (KETEK, SensL, Hamamatsu, Broadcom) were tested in combination with the ASIC. Scintillator arrays featuring different reflector designs and different configurations of the TOFPET2 ASIC software parameters were evaluated. The benchtop setup used is provided with the TOFPET2 ASIC evaluation kit by PETsys Electronics S.A. RESULTS: Compared to existing studies featuring the TOFPET2 ASIC, multi-channel performance results dependent on a larger set of ASIC configuration parameters were obtained that have not been reported to this extend so far. The ASIC shows promising CRTs down to 219.9 ps in combination with two Hamamatsu S14161-3050-HS-08 SiPM arrays (128 channels read out, energy resolution 13.08%) and 216.1 ps in combination with two Broadcom AFBR-S4N44P643S SiPM arrays (32 channels read out, energy resolution 9.46%). The length of the trigger delay of the dark count suppression scheme has an impact on the ASIC performance and can be configured to further improve the coincidence resolution time. The integrator gain configuration has been investigated and allows an absolute improvement of the energy resolution by up to 1% at the cost of the linearity of the energy spectrum. CONCLUSION: Measuring up to the time-of-flight performance of state-of-the-art positron emission tomography (ToF-PET) systems while providing a uniform and stable readout for multiple channels at the same time, the TOFPET2 ASIC is treated as promising candidate for the integration in future ToF-PET systems.

Comments on the NEMA NU 4-2008 Standard on Performance Measurement of Small Animal Positron Emission Tomographs.

The National Electrical Manufacturers Association's (NEMA) NU 4-2008 standard specifies methodology for evaluating the performance of small-animal PET scanners. The standard's goal is to enable comparison of different PET scanners over a wide range of technologies and geometries used. In this work, we discuss if the NEMA standard meets these goals and we point out potential flaws and improvements to the standard.For the evaluation of spatial resolution, the NEMA standard mandates the use of filtered backprojection reconstruction. This reconstruction method can introduce star-like artifacts for detectors with an anisotropic spatial resolution, usually caused by parallax error. These artifacts can then cause a strong dependence of the resulting spatial resolution on the size of the projection window in image space, whose size is not fully specified in the NEMA standard. If the PET ring has detectors which are perpendicular to a Cartesian axis, then the resolution along this axis will typically improve with larger projection windows.We show that the standard's equations for the estimation of the random rate for PET systems with intrinsic radioactivity are circular and not satisfiable. However, a modified version can still be used to determine an approximation of the random rates under the assumption of negligible random rates for small activities and a constant scatter fraction. We compare the resulting estimated random rates to random rates obtained using a delayed coincidence window and two methods based on the singles rates. While these methods give similar estimates, the estimation method based on the NEMA equations overestimates the random rates.In the NEMA standard's protocol for the evaluation of the sensitivity, the standard specifies to axially step a point source through the scanner and to take a different scan for each source position. Later, in the data analysis section, the standard does not specify clearly how the different scans have to be incorporated into the analysis, which can lead to unclear interpretations of publicized results.The standard's definition of the recovery coefficients in the image quality phantom includes the maximum activity in a region of interest, which causes a positive correlation of noise and recovery coefficients. This leads to an unintended trade-off between desired uniformity, which is negatively correlated with variance (i.e., noise), and recovery.With this work, we want to start a discussion on possible improvements in a next version of the NEMA NU-4 standard.

Characterization methods for comprehensive evaluations of shielding materials used in an MRI.

PURPOSE: In order to integrate electronic devices into a magnetic resonance imaging (MRI) scanner, shielding of the electronics with respect to the radio frequency (RF) transmit and receive system of the MRI scanner is required. Furthermore, MRI uses time-varying low-frequency magnetic fields for spatial encoding, i.e., the gradient magnetic fields. Time-varying magnetic fields induce eddy currents in all conductive elements. The eddy currents result in opposing magnetic fields, which can cause distortions of the magnetic resonance (MR) image. As shielding of lower frequencies is not feasible in this respect, an ideal shielding element should be transparent for gradient magnetic fields while providing a high RF shielding effectiveness. Furthermore, it should offer a low susceptibility to prevent distortion of the main magnetic field of the MRI. In this work, we characterize the aforesaid shielding parameters of different shielding samples. METHODS: We developed a nuclear magnetic resonance (NMR) probe to measure the magnetic fields to quantify the field distortions time-resolvedly. The relative distortion was introduced as a proportionality constant relating the eddy-current-inducing field changes and the field distortions. The relative distortion was measured in the frequency range from 0 to 10 kHz for all shielding samples using the NMR probe. We characterized the shielding effectiveness of the samples in the frequency range from 1 to 150 MHz using a network analyzer. We conducted all measurements with three different materials, two carbon fiber composites and copper, each in various thicknesses. RESULTS: The relative distortion of the magnetic fields induced by the carbon fiber composites samples was at least a factor of seven lower than the copper sample. A linear dependency on the sample thickness was measured for the main field distortion, the relative distortion and the shielding effectiveness. The relative distortion was roughly independent of the gradient frequency contrary to the shielding effectiveness, highly depending on the RF frequency. CONCLUSIONS: We presented a very sensitive method to characterize the distortion of the main field distortion and the gradient transparency using an NMR probe. We analyzed different shielding materials regarding the main field distortion, the gradient transparency, and the shielding effectiveness. From the tested materials, we identified a carbon fiber composite with the lowest distortion on the MRI.

PET performance evaluation of the small-animal Hyperion IID PET/MRI insert based on the NEMA NU-4 standard.

The Hyperion IID PET insert is the first scanner using fully digital silicon photomultipliers for simultaneous PET/MR imaging of small animals up to rabbit size. In this work, we evaluate the PET performance based on the National Eletrical Manufacturers Association (NEMA) NU 4-2008 standard, whose standardized measurement protocols allow comparison of different small-animal PET scanners. The Hyperion IID small-animal PET/MR insert comprises three rings of 20 detector stacks with pixelated scintillator arrays with a crystal pitch of 1 mm, read out with digital silicon photomultipliers. The scanner has a large ring diameter of 209.6 mm and an axial field of view of 96.7 mm. We evaluated the spatial resolution, energy resolution, time resolution and sensitivity by scanning a 22Na point source. The count rates and scatter fractions were measured for a wide range of 18F activity inside a mouse-sized scatter phantom. We evaluated the image quality using the mouse-sized image quality phantom specified in the NEMA NU4 standard, filled with 18F. Additionally, we verified the in-vivo imaging capabilities by performing a simultaneous PET/MRI scan of a mouse injected with 18F-FDG. We processed all measurement data with an energy window of 250 keV to 625 keV and a coincidence time window of 2 ns. The filtered-backprojection reconstruction of the point source has a full width at half maximum (FWHM) of 1.7 mm near the isocenter and degrades to 2.5 mm at a radial distance of 50 mm. The scanner's average energy resolution is 12.7% (ΔE/E FWHM) and the coincidence resolution time is 609 ps. The peak absolute sensitivity is 4.0% and the true and noise-equivalent count rates reach their peak at an activity of 46 MBq with 483 kcps and 407 kcps, respectively, with a scatter fraction of 13%. The iterative reconstruction of the image quality phantom has a uniformity of 3.7%, and recovery coefficients from 0.29, 0.91 and 0.94 for rod diameters of 1 mm, 3 mm and 5 mm, respectively. After application of scatter and attenuation corrections, the air- and water-filled cold regions have spill-over ratios of 6.3% and 5.4%, respectively. The Hyperion IID PET/MR insert provides state-of-the-art PET performance while enabling simultaneous PET/MRI acquisition of small animals up to rabbit size.

Maximum likelihood positioning algorithm for high-resolution PET scanners.

PURPOSE: In high-resolution positron emission tomography (PET), lightsharing elements are incorporated into typical detector stacks to read out scintillator arrays in which one scintillator element (crystal) is smaller than the size of the readout channel. In order to identify the hit crystal by means of the measured light distribution, a positioning algorithm is required. One commonly applied positioning algorithm uses the center of gravity (COG) of the measured light distribution. The COG algorithm is limited in spatial resolution by noise and intercrystal Compton scatter. The purpose of this work is to develop a positioning algorithm which overcomes this limitation. METHODS: The authors present a maximum likelihood (ML) algorithm which compares a set of expected light distributions given by probability density functions (PDFs) with the measured light distribution. Instead of modeling the PDFs by using an analytical model, the PDFs of the proposed ML algorithm are generated assuming a single-gamma-interaction model from measured data. The algorithm was evaluated with a hot-rod phantom measurement acquired with the preclinical hyperion II (D) PET scanner. In order to assess the performance with respect to sensitivity, energy resolution, and image quality, the ML algorithm was compared to a COG algorithm which calculates the COG from a restricted set of channels. The authors studied the energy resolution of the ML and the COG algorithm regarding incomplete light distributions (missing channel information caused by detector dead time). Furthermore, the authors investigated the effects of using a filter based on the likelihood values on sensitivity, energy resolution, and image quality. RESULTS: A sensitivity gain of up to 19% was demonstrated in comparison to the COG algorithm for the selected operation parameters. Energy resolution and image quality were on a similar level for both algorithms. Additionally, the authors demonstrated that the performance of the ML algorithm is less prone to missing channel information. A likelihood filter visually improved the image quality, i.e., the peak-to-valley increased up to a factor of 3 for 2-mm-diameter phantom rods by rejecting 87% of the coincidences. A relative improvement of the energy resolution of up to 12.8% was also measured rejecting 91% of the coincidences. CONCLUSIONS: The developed ML algorithm increases the sensitivity by correctly handling missing channel information without influencing energy resolution or image quality. Furthermore, the authors showed that energy resolution and image quality can be improved substantially by rejecting events that do not comply well with the single-gamma-interaction model, such as Compton-scattered events.

Initial PET performance evaluation of a preclinical insert for PET/MRI with digital SiPM technology.

Hyperion-II(D) is a positron emission tomography (PET) insert which allows simultaneous operation in a clinical magnetic resonance imaging (MRI) scanner. To read out the scintillation light of the employed lutetium yttrium orthosilicate crystal arrays with a pitch of 1 mm and 12 mm in height, digital silicon photomultipliers (DPC 3200-22, Philips Digital Photon Counting) (DPC) are used. The basic PET performance in terms of energy resolution, coincidence resolution time (CRT) and sensitivity as a function of the operating parameters, such as the operating temperature, the applied overvoltage, activity and configuration parameters of the DPCs, has been evaluated at system level. The measured energy resolution did not show a large dependency on the selected parameters and is in the range of 12.4%-12.9% for low activity, degrading to ∼13.6% at an activity of ∼100 MBq. The CRT strongly depends on the selected trigger scheme (trig) of the DPCs, and we measured approximately 260 ps, 440 ps, 550 ps and 1300 ps for trig 1-4, respectively. The trues sensitivity for a NEMA NU 4 mouse-sized scatter phantom with a 70 mm long tube of activity was dependent on the operating parameters and was determined to be 0.4%-1.4% at low activity. The random fraction stayed below 5% at activity up to 100 MBq and the scatter fraction was evaluated as ∼6% for an energy window of 411 keV-561 keV and ∼16% for 250 keV-625 keV. Furthermore, we performed imaging experiments using a mouse-sized hot-rod phantom and a large rabbit-sized phantom. In 2D slices of the reconstructed mouse-sized hot-rod phantom (∅ = 28 mm), the rods were distinguishable from each other down to a rod size of 0.8 mm. There was no benefit from the better CRT of trig 1 over trig 3, where in the larger rabbit-sized phantom (∅ = 114 mm) we were able to show a clear improvement in image quality using the time-of-flight information. The findings will allow system architects-aiming at a similar detector design using DPCs-to make predictions about the design requirements and the performance that can be expected.

Software-Based Real-Time Acquisition and Processing of PET Detector Raw Data.

In modern positron emission tomography (PET) readout architectures, the position and energy estimation of scintillation events (singles) and the detection of coincident events (coincidences) are typically carried out on highly integrated, programmable printed circuit boards. The implementation of advanced singles and coincidence processing (SCP) algorithms for these architectures is often limited by the strict constraints of hardware-based data processing. In this paper, we present a software-based data acquisition and processing architecture (DAPA) that offers a high degree of flexibility for advanced SCP algorithms through relaxed real-time constraints and an easily extendible data processing framework. The DAPA is designed to acquire detector raw data from independent (but synchronized) detector modules and process the data for singles and coincidences in real-time using a center-of-gravity (COG)-based, a least-squares (LS)-based, or a maximum-likelihood (ML)-based crystal position and energy estimation approach (CPEEA). To test the DAPA, we adapted it to a preclinical PET detector that outputs detector raw data from 60 independent digital silicon photomultiplier (dSiPM)-based detector stacks and evaluated it with a [(18)F]-fluorodeoxyglucose-filled hot-rod phantom. The DAPA is highly reliable with less than 0.1% of all detector raw data lost or corrupted. For high validation thresholds (37.1 ± 12.8 photons per pixel) of the dSiPM detector tiles, the DAPA is real time capable up to 55 MBq for the COG-based CPEEA, up to 31 MBq for the LS-based CPEEA, and up to 28 MBq for the ML-based CPEEA. Compared to the COG-based CPEEA, the rods in the image reconstruction of the hot-rod phantom are only slightly better separable and less blurred for the LS- and ML-based CPEEA. While the coincidence time resolution (∼ 500 ps) and energy resolution (∼12.3%) are comparable for all three CPEEA, the system sensitivity is up to 2.5 × higher for the LS- and ML-based CPEEA.

Development of an MRI-compatible digital SiPM detector stack for simultaneous PET/MRI.

Advances in solid-state photon detectors paved the way to combine positron emission tomography (PET) and magnetic resonance imaging (MRI) into highly integrated, truly simultaneous, hybrid imaging systems. Based on the most recent digital SiPM technology, we developed an MRI-compatible PET detector stack, intended as a building block for next generation simultaneous PET/MRI systems. Our detector stack comprises an array of 8 × 8 digital SiPM channels with 4 mm pitch using Philips Digital Photon Counting DPC 3200-22 devices, an FPGA for data acquisition, a supply voltage control system and a cooling infrastructure. This is the first detector design that allows the operation of digital SiPMs simultaneously inside an MRI system. We tested and optimized the MRI-compatibility of our detector stack on a laboratory test bench as well as in combination with a Philips Achieva 3 T MRI system. Our design clearly reduces distortions of the static magnetic field compared to a conventional design. The MRI static magnetic field causes weak and directional drift effects on voltage regulators, but has no direct impact on detector performance. MRI gradient switching initially degraded energy and timing resolution. Both distortions could be ascribed to voltage variations induced on the bias and the FPGA core voltage supply respectively. Based on these findings, we improved our detector design and our final design shows virtually no energy or timing degradations, even during heavy and continuous MRI gradient switching. In particular, we found no evidence that the performance of the DPC 3200-22 digital SiPM itself is degraded by the MRI system.

PET performance and MRI compatibility evaluation of a digital, ToF-capable PET/MRI insert equipped with clinical scintillators.

We evaluate the MR compatibility of the Hyperion-II(D) positron emission tomography (PET) insert, which allows simultaneous operation in a clinical magnetic resonance imaging (MRI) scanner. In contrast to previous investigations, this work aims at the evaluation of a clinical crystal configuration. An imaging-capable demonstrator with an axial field-of-view of 32 mm and a crystal-to-crystal spacing of 217.6 mm was equipped with LYSO scintillators with a pitch of 4 mm which were read out in a one-to-one coupling scheme by sensor tiles composed of digital silicon photomultipliers from Philips Digital Photon Counting (DPC 3200-22). The PET performance degradation (energy resolution and coincidence resolution time (CRT)) was evaluated during simultaneous operation of the MRI scanner. We used clinically motivated imaging sequences as well as synthetic gradient stress test sequences. Without activity of the MRI scanner, we measured for trigger scheme 1 (first photon trigger) an energy resolution of 11.4% and a CRT of 213 ps for a narrow energy (NE) window using five (22)Na point-like sources. When applying the synthetic gradient sequences, we found worst-case relative degradations of the energy resolution by 5.1% and of the CRT by 33.9%. After identifying the origin of the degradations and implementing a fix to the read-out hardware, the same evaluation revealed no degradation of the PET performance anymore even when the most demanding gradient stress tests were applied. The PET performance of the insert was initially evaluated using the point sources, a high-activity phantom and hot-rod phantoms in order to assess the spatial resolution. Trigger schemes 2-4 delivered an energy resolution of 11.4% as well and CRTs of 279 ps, 333 ps and 557 ps for the NE window, respectively. An isocenter sensitivity of 0.41% using the NE window and 0.71% with a wide energy window was measured. Using a hot-rod phantom, a spatial resolution in the order of 2 mm was demonstrated and the benefit of time-of-flight PET was shown with a larger rabbit-sized phantom. In conclusion, the Hyperion architecture is an interesting platform for clinically driven hybrid PET/MRI systems.

A Digital Preclinical PET/MRI Insert and Initial Results.

Combining Positron Emission Tomography (PET) with Magnetic Resonance Imaging (MRI) results in a promising hybrid molecular imaging modality as it unifies the high sensitivity of PET for molecular and cellular processes with the functional and anatomical information from MRI. Digital Silicon Photomultipliers (dSiPMs) are the digital evolution in scintillation light detector technology and promise high PET SNR. DSiPMs from Philips Digital Photon Counting (PDPC) were used to develop a preclinical PET/RF gantry with 1-mm scintillation crystal pitch as an insert for clinical MRI scanners. With three exchangeable RF coils, the hybrid field of view has a maximum size of 160 mm × 96.6 mm (transaxial × axial). 0.1 ppm volume-root-mean-square B 0-homogeneity is kept within a spherical diameter of 96 mm (automatic volume shimming). Depending on the coil, MRI SNR is decreased by 13% or 5% by the PET system. PET count rates, energy resolution of 12.6% FWHM, and spatial resolution of 0.73 mm (3) (isometric volume resolution at isocenter) are not affected by applied MRI sequences. PET time resolution of 565 ps (FWHM) degraded by 6 ps during an EPI sequence. Timing-optimized settings yielded 260 ps time resolution. PET and MR images of a hot-rod phantom show no visible differences when the other modality was in operation and both resolve 0.8-mm rods. Versatility of the insert is shown by successfully combining multi-nuclei MRI ((1)H/(19)F) with simultaneously measured PET ((18)F-FDG). A longitudinal study of a tumor-bearing mouse verifies the operability, stability, and in vivo capabilities of the system. Cardiac- and respiratory-gated PET/MRI motion-capturing (CINE) images of the mouse heart demonstrate the advantage of simultaneous acquisition for temporal and spatial image registration.