NEMA NU-4 performance evaluation of PETbox4, a high sensitivity dedicated PET preclinical tomograph

Timing performance evaluation of a dual-Axis rotational PET system according to NEMA NU 4-2008 standards: A simulation study

Introduction:Positron Emission Tomography (PET) imaging's diagnostic accuracy is dependent on the scanner design and image quality, which is affected by several factors including the coincidence timing window (CTW). NEMA NU 4-2008 procedures are commonly used to assess and compare PET systems performance, including dual rotation technologies like easyPET.3D, known for high-spatial resolution and reduced parallax contribution.Aim:This study aims to identify easyPET.3D's optimal performance based on NEMA standards. In addition, explores the impact of different CTWs on PET image quality by comparing simulated electronics capable of a 300 ps CTW with a 40 ns CTW.Results:When the data is filtered by a 40 ns CTW, a sub-millimetre resolution at the field-of-view (FoV) centre and a constant behaviour in the radial direction are achieved. The absolute sensitivity was 0.18% with a maximum value of 0.31%, for a 15 mm transverse FoV. The noise equivalent count rate peaked at 18 MBq with 249 cps. Recovery coefficients ranged from 17% to 90%, and spilled-over ratios were 0.32 (water) and 0.41 (air).Conclusions:A shorter 300 ps CTW primarily impacted PET dynamic range, allowing higher activity acquisitions, with no significant changes in resolution and sensitivity under NEMA test conditions. As for the image quality test, the 300 ps CTW images have less background, better SOR values, and similar RC values when comparing the 40 ns CTW.

Performance evaluation of a small-animal PET scanner with 213 mm axis using NEMA NU 4-2008

BACKGROUND

Long-axis positron emission tomography (PET) has emerged as one of the recent research directions in PET due to its ability to significantly enhance sensitivity and counting performance for low-dose imaging, rapid imaging, and whole-body dynamic imaging.

PURPOSE

The PET system presented in this study is a long-axis animal PET based on lutetium-yttrium orthosilicate and silicon photomultiplier, designed for whole-body imaging in rats. It features a diameter of 143 mm and an axial length of 213.3 mm. This study evaluated the performance of this PET system in accordance with the National Electrical Manufacturers Association (NEMA) NU 4-2008 standards.

METHODS

The performance evaluation was conducted according to the NEMA NU 4-2008 standards in terms of spatial resolution, sensitivity, counting rate performance, scatter fraction (SF) and image quality. In addition, a rat imaging study was conducted to assess the imaging capability of this PET system.

RESULTS

The average energy resolution of the PET system was 12.87%, the average coincidence timing resolution was 751 ps. The FWHM of spatial resolution reconstructed by filtered back projection and 3D-OSEM-PSF algorithm at 5 mm radial offset from the axial center were 1.65 and 0.88 mm. The peak absolute sensitivity measured by a point source at the center of the field of view was evaluated as 6.71% (361-661 keV) and 10.31% (250-750 keV). For the mouse-like phantom, the SF was 11.0% and the peak noise equivalent counting rate (NECR) was 1193 kcps at 94.2 MBq (2.54 mCi). For the rat-like phantom, the SF was 26.8% and the NECR was 682.5 kcps at 78.6 MBq (2.12 mCi).

CONCLUSIONS

The performance measurement results demonstrate that this PET system exhibits high sensitivity and count rate performance, making it potential for high-quality whole-body dynamic imaging of rats.

A depth-of-interaction rebinning method based on both geometric and activity weights

BACKGROUND AND OBJECTIVE

For positron emission tomography (PET) scanners with depth-of-interaction (DOI) measurement, the DOI rebinning method that utilizes DOI information to process the projection data is critical to image quality. Current DOI rebinning methods map coincidence events onto the rebinned sinogram based on the correlation of lines of response (LOR). This study aims to incorporate prior radioactivity distribution of the imaging object into DOI rebinning to obtain better image quality.

METHODS

A DOI rebinning method based on both geometric and activity weights was proposed to assign coincidence events to the rebinned sinogram defined by a virtual ring. The geometric weights, representing the correlation between LORs, were calculated based on the areas of intersection. The activity weights, reflecting the activity distribution of the imaging object, were derived from the previous reconstructed image.

RESULTS

Monte Carlo simulation data from four phantoms, including the image quality phantom, Derenzo phantom, and two rat-like ROBY phantoms, was used to evaluate the proposed method. The recovery coefficient (RC), contrast recovery coefficient (CRC), structural similarity index measure (SSIM), and peak signal-to-noise ratio (PSNR) were used as image quality metrics. Compared to other DOI rebinning methods, the proposed method achieved the highest RC (maximum improvement of 32%) and CRC at the same noise level and was also optimal in terms of the SSIM and PSNR. Meanwhile, incorporating the prior activity distribution into DOI rebinning also improved the image reconstruction speed.

CONCLUSIONS

This work developed a new DOI rebinning method combining the correlation of LORs with the prior activity distribution, achieving relatively optimal image quality and reconstruction speed. Furthermore, it still needs to be evaluated on the actual equipment.

A guideline proposal for mice preparation and care in 18F-FDG PET imaging

The experimental outcomes of small-animal positron emission tomography (PET) imaging with ¹⁸F-labelled fluorodeoxyglucose (¹⁸F-FDG) can be particularly compromised by animal preparation and care. Several works intend to improve research reporting and amplify the quality and reliability of published research. Though these works provide valuable information to plan and conduct animal studies, manuscripts describe different methodologies—standardization does not exist. Consequently, the variation in details reported can explain the difference in the experimental results found in the literature. Additionally, the resources and guidelines defining protocols for small-animal imaging are scarce, making it difficult for researchers to obtain and compare accurate and reproducible data. Considering the selection of suitable procedures key to ensure animal welfare and research improvement, this paper aims to prepare the way for a future guideline on mice preparation and care for PET imaging with ¹⁸F-FDG. For this purpose, a global standard protocol was created based on recommendations and good practices described in relevant literature.

Performance Evaluation of a PET of 7T Bruker Micro-PET/MR Based on NEMA NU 4-2008 Standards

Purpose: This study aimed to measure the performance evaluation of the Bruker sequential micro-positron emission tomography/magnetic resonance imaging (PET/MRI) scanner by following National Electrical Manufacturers Association (NEMA) NU 4-2008 standards’ protocol. The system consists of a high-performance silicon photomultiplier (SiPM) advanced technology detector and a continuous lutetium-yttrium oxyorthosilicate (LYSO) crystal. Methods: A 22Na (sodium-22) point source was utilized to assess the spatial resolution and system sensitivity, and the Micro-PET scatter phantom measurements were conducted to measure count rate measurements and scatter fractions (SF). A mouse-like Micro-PET image quality (IQ) phantom was utilized as a model to analyze the uniformity, recovery coefficient (RC), and spillover ratio (SOR). A small animal PET/MRI imaging study was performed in a rat. Results: We calculated the spatial resolutions of filtered back-projection (FBP), and used 3D-MLEM to reconstruct PET images at the axial center and ¼ of the axial field of view (FOV) in axial, radial, and tangential directions. The best observed spatial resolutions in both reconstructed images were obtained in the tangential direction, and the values were 0.80 mm in 3D-MLEM and 0.94 mm in FBP. The peak noise equivalent count rate (NECR) in the 358–664 keV energy window was 477.30 kcps at 95.83 MBq and 774.45 kcps at 103.6 MBq for rat and mouse-sized scatter phantoms, respectively. The rat and mouse-sized phantoms scatter fractions (SF) were 14.2% and 6.9%, respectively. Conclusions: According to our results, the performance characteristics of the scanner are high sensitivity, good spatial resolution, low scatter fraction, and good IQ, indicating that it is suitable for preclinical imaging studies.

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.

H2RSPET: a 0.5 mm resolution high-sensitivity small-animal PET scanner, a simulation study

With the goal of developing a total-body small-animal PET system with a high spatial resolution of ∼0.5 mm and a high sensitivity >10% for mouse/rat studies, we simulated four scanners using the graphical processing unit-based Monte Carlo simulation package (gPET) and compared their performance in terms of spatial resolution and sensitivity. We also investigated the effect of depth-of-interaction (DOI) resolution on the spatial resolution. All the scanners are built upon 128 DOI encoding dual-ended readout detectors with lutetium yttrium oxyorthosilicate (LYSO) arrays arranged in 8 detector rings. The solid angle coverages of the four scanners are all ∼0.85 steradians. Each LYSO element has a cross-section of 0.44 × 0.44 mm2 and the pitch size of the LYSO arrays are all 0.5 mm. The four scanners can be divided into two groups: (1) H2RS110-C10 and H2RS110-C20 with 40 × 40 LYSO arrays, a ring diameter of 110 mm and axial length of 167 mm, and (2) H2RS160-C10 and H2RS160-C20 with 60 × 60 LYSO arrays, a diameter of 160 mm and axial length of 254 mm. C10 and C20 denote the crystal thickness of 10 and 20 mm, respectively. The simulation results show that all scanners have a spatial resolution better than 0.5 mm at the center of the field-of-view (FOV). The radial resolution strongly depends on the DOI resolution and radial offset, but not the axial resolution and tangential resolution. Comparing the C10 and C20 designs, the former provides better resolution, especially at positions away from the center of the FOV, whereas the latter has 2× higher sensitivity (∼10% versus ∼20%). This simulation study provides evidence that the 110 mm systems are a good choice for total-body mouse studies at a lower cost, whereas the 160 mm systems are suited for both total-body mouse and rat studies.

A digital phoswich detector using time-over-threshold for depth of interaction in PET

We present the performance of a digital phoswich positron emission tomography (PET) detector, composed by layers of pixilated scintillator arrays, read out by solid state light detectors and an application specific integrated circuit (ASIC). We investigated the use of integrated charge from the scintillation pulses along with time-over-threshold (ToT) to determine the layer of interaction (DOI) in the scintillator. Simulations were performed to assess the effectiveness of the ToT measurements for separating the scintillator events and identifying cross-layer-crystal-scatter (CLCS) events. These simulations indicate that ToT and charge integration from such a detector provide sufficient information to determine the layer of interaction. To demonstrate this in practice, we used a pair of prototype LYSO/BGO detectors. One detector consisted of a 19 × 19 array of 7 mm long LYSO crystals (1.36 mm pitch) coupled to a 16 × 16 array of 8 mm long BGO crystals (1.63 mm pitch). The other detector was similar except the LYSO crystal pitch was 1.63 mm. These detectors were coupled to an 8 × 8 multi-pixel photon counter mounted on a PETsys TOFPET2 ASIC. This high performance ASIC provided digital readout of the integrated charge and ToT from these detectors. We present a method to separate the events from the two scintillator layers using the ToT, and also investigate the performance of this detector. All the crystals within the proposed detector were clearly resolved, and the peak to valley ratio was 11.8 ± 4.0 and 10.1 ± 2.9 for the LYSO and BGO flood images. The measured energy resolution was 9.9% ± 1.3% and 28.5% ± 5.0% respectively for the LYSO and BGO crystals in the phoswich layers. The timing resolution between the LYSO-LYSO, LYSO-BGO and BGO-BGO coincidences was 468 ps, 1.33 ns and 2.14 ns respectively. Results show ToT can be used to identify the crystal layer where events occurred and also identify and reject the majority of CLCS events between layers.

Design and performance of SIAT aPET: a uniform high-resolution small animal PET scanner using dual-ended readout detectors

In this work, a small animal PET scanner named SIAT aPET was developed using dual-ended readout depth encoding detectors to simultaneously achieve high spatial resolution and high sensitivity. The scanner consists of four detector rings with 12 detector modules per ring; the ring diameter is 111 mm and the axial field of view (FOV) is 105.6 mm. The images are reconstructed using an ordered subset expectation maximization (OSEM) algorithm. The spatial resolution of the scanner was measured by using a ²²Na point source at the center axial FOV with different radial offsets. The sensitivity of the scanner was measured at center axis of the scanner with different axial positions. The count rate performance of the system was evaluated by scanning mouse-sized and rat-sized phantoms. An ultra-micro hot-rods phantom and two mice injected with ¹⁸F-NaF and ¹⁸F-FDG were scanned on the scanner. An average depth of interaction (DOI) resolution of 1.96 mm, energy resolution of 19.1% and timing resolution of 1.20 ns were obtained for the detector. Average spatial resolutions of 0.82 mm and 1.16 mm were obtained up to a distance of 30 mm radially from the center of the FOV when reconstructing a point source in 1% and 10% warm backgrounds, respectively, using OSEM reconstruction with 16 subsets and 10 iterations. Sensitivities of 16.0% and 11.9% were achieved at center of the scanner for energy windows of 250-750 keV and 350-750 keV respectively. Peak noise equivalent count rates (NECRs) of 324 kcps and 144 kcps were obtained at an activity of 26.4 MBq for the mouse-sized and rat-sized phantoms. Rods of 1.0 mm diameter can be visually resolved from the image of the ultra-micro hot-rods phantom. The capability of the scanner was demonstrated by high quality in-vivo mouse images.

Advances in Preclinical PET Instrumentation

In the light of ever-increasing demands for PET scanner with better resolvability, higher sensitivity and wide accessibility for noninvasive screening of small structures and physiological processes in laboratory rodents, several dedicated PET scanners were developed and evaluated. Understanding conceptual design constraints pros and cons of different configurations and impact of the major components will be helpful to further establish the crucial role of these miniaturized systems in a broad spectrum of modern research. Hence, a comprehensive review of preclinical PET scanners developed till early 2020 with particular emphasis on innovations in instrumentation and geometrical designs is provided.

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.

A mini-panel PET scanner-based microfluidic radiobioassay system allowing high-throughput imaging of real-time cellular pharmacokinetics

On-chip radiometric detection of biological samples using radiotracers has become an emerging research field known as microfluidic radiobioassays. Performing parallel radiobioassays is highly desirable for saving time/effort, reducing experimental variation between assays, and minimizing the cost of the radioisotope. Continuously infused microfluidic radioassay (CIMR) is one of the useful tools for investigating cellular pharmacokinetics and assessing the binding and uptakes of radiopharmaceuticals. However, existing CIMR systems can only measure the dynamics of one sample at a time due to the limited field of view (FOV) of the positron detector. To increase the throughput, we propose a new CIMR system with a custom-built miniaturized panel-based positron-emission tomography (PET) scanner and a parallel infusion setup/method, capable of imaging the cellular pharmacokinetics of three samples in one measurement. With this system, the pharmacokinetics of parallel or comparison samples can be imaged simultaneously. The increased throughput is attributed to two innovations: 1) the large 3D FOV of the mini-panel PET scanner, enabling more samples to be imaged in the microfluidic chip; and 2) a parallel infusion method, in which only one reference chamber is needed for indicating the dynamic input of the infused radiotracer medium, thus saving the total reference chambers needed compared to the current sequential infusion method. Combining the CIMR technique and the mini-panel PET scanner, this study also firstly demonstrated the feasibility of using PET, as an imaging modality, for microfluidic radiobioassays. Besides the increased throughput, the 3D imaging of PET also provides possibilities for further applications such as organoid/3D culturing systems, non-planar microfluidics, and organs-on-chips. The system is more practical for a broader range of applications in nuclear medicine, molecular imaging, and lab-on-a-chip studies.

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.

NEMA NU-4 2008 performance evaluation of Xtrim-PET: A prototype SiPM-based preclinical scanner

PURPOSE

Xtrim-PET is a newly designed Silicon Photomultipliers (SiPMs)-based prototype PET scanner dedicated for small laboratory animal imaging. We present the performance evaluation of the Xtrim-PET scanner following NEMA NU-4 2008 standards to help optimizing scanning protocols which can be achieved through standard and reliable system performance characterization.

METHODS

The performance assessment was conducted according to the National Electrical Manufacturers Association (NEMA) NU-4 2008 standards in terms of spatial resolution, sensitivity, counting rate performance, scatter fraction and image quality. The in vivo imaging capability of the scanner is also showcased through scanning a normal mouse injected with ¹⁸ F-FDG. Furthermore, the performance characteristics of the developed scanner are compared with commercially available systems and current prototypes.

RESULTS

The volumetric spatial resolution at 5 mm radial offset from the central axis of the scanner is 6.81 µl, whereas a peak absolute sensitivity of 2.99% was achieved using a 250-650 keV energy window and a 10 ns timing window. The peak noise-equivalent count rate (NECR) using a mouse-like phantom is 113.18 kcps at 0.34 KBq/cc with 12.5% scatter fraction, whereas the NECR peaked at 82.76 kcps for an activity concentration level of 0.048 KBq/cc with a scatter fraction of 25.8% for rat-like phantom. An excellent uniformity (3.8%) was obtained using NEMA image quality phantom. Recovery coefficients of 90%, 86%, 68%, 40% and 12% were calculated for rod diameters of 5, 4, 3, 2 and 1 mm, respectively. Spill-over ratios for air-filled and water-filled chambers were 35% and 25% without applying any correction for attenuation and Compton scattering effects.

CONCLUSION

Our findings revealed that beyond compactness, lightweight, easy installation and good energy resolution, the Xtrim-PET prototype presents a reasonable performance making it suitable for preclinical molecular imaging-based research.

Performance evaluation of a small animal PET scanner a high level of multiplexing and charge-signal transmission

Small animal positron emission tomography (PET) is a noninvasive imaging modality that enables in vivo imaging and quantification of the biological processes of small experimental animals. We have developed a small animal PET that utilizes a high-resolution multiplexed readout and charge signal transmission (CST) method. The small animal PET was composed of six detector blocks consisting of SiPMs and LYSO arrays. Six detector blocks were mounted on a PET gantry having an inner diameter of 76 mm, outer diameter of 112 mm, and axial length of 40.8 mm. The charge signals of SiPM output were transmitted to the input of multiplexed readout using 4 m flexible flat cables. The multiplexed readout was composed of six main boards, each of which included 36 detector boards, to reduce the number of readout channels by a factor of 36, with a multiplexing ratio of 144:4. The performance of the small animal PET was evaluated using NEMA NU 4-2008 standards, and its imaging capability was demonstrated by in vivo mouse imaging studies. The average energy and time resolutions were 13.2%  ±  0.3% and 3.8 ns, respectively. The spatial resolution at the center of the transaxial FOV was 1.1 mm, and the peak sensitivity at the center of the axial FOV was 1.5%. The peak noise equivalent count (NEC) rate and scatter fraction were 21.1 kcps at 18.2 MBq and 21%, respectively. The acquired images demonstrated high quality tracer uptake patterns of small experimental animals. The results of performance evaluation and animal imaging indicate that the small animal PET developed in this study can provide high-quality small animal imaging with cost-effectiveness and compactness.

Performance Evaluation of G8, a High-Sensitivity Benchtop Preclinical PET/CT Tomograph

G8 is a benchtop integrated PET/CT scanner dedicated to high-sensitivity and high-resolution imaging of mice. This work characterizes its National Electrical Manufacturers Association NU 4-2008 performance where applicable and also assesses the basic imaging performance of the CT subsystem. Methods: The PET subsystem in G8 consists of 4 flat-panel detectors arranged in a boxlike geometry. Each panel consists of 2 modules of a 26 × 26 pixelated bismuth germanate scintillator array with individual crystals measuring 1.75 × 1.75 × 7.2 mm. The crystal arrays are coupled to multichannel photomultiplier tubes via a tapered, pixelated glass lightguide. A cone-beam CT scanner consisting of a MicroFocus x-ray source and a complementary metal oxide semiconductor detector provides anatomic information. Sensitivity, spatial resolution, energy resolution, scatter fraction, count-rate performance, and the capability of performing phantom and mouse imaging were evaluated for the PET subsystem. Noise, dose level, contrast, and resolution were evaluated for the CT subsystem. Results: With an energy window of 350-650 keV, the peak sensitivity was 9.0% near the center of the field of view. The crystal energy resolution ranged from 15.0% to 69.6% in full width at half maximum (FWHM), with a mean of 19.3% ± 3.7%. The average intrinsic spatial resolution was 1.30 and 1.38 mm FWHM in the transverse and axial directions, respectively. The maximum-likelihood expectation maximization reconstructed image of a point source in air averaged 0.81 ± 0.11 mm FWHM. The peak noise-equivalent count rate for the mouse-sized phantom was 44 kcps for a total activity of 2.9 MBq (78 μCi), and the scatter fraction was 11%. For the CT subsystem, the value of the modulation transfer function at 10% was 2.05 cycles/mm. Conclusion: The overall performance demonstrates that the G8 can produce high-quality images for molecular imaging-based biomedical research.

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

The Hyperion II^D 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 II^D 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 ²²Na point source. The count rates and scatter fractions were measured for a wide range of ¹⁸F 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 ¹⁸F. Additionally, we verified the in-vivo imaging capabilities by performing a simultaneous PET/MRI scan of a mouse injected with ¹⁸F-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 II^D PET/MR insert provides state-of-the-art PET performance while enabling simultaneous PET/MRI acquisition of small animals up to rabbit size.

PET performance evaluation of MADPET4: a small animal PET insert for a 7 T MRI scanner

MADPET4 is the first small animal PET insert with two layers of individually read out crystals in combination with silicon photomultiplier technology. It has a novel detector arrangement, in which all crystals face the center of field of view transaxially. In this work, the PET performance of MADPET4 was evaluated and compared to other preclinical PET scanners using the NEMA NU 4 measurements, followed by imaging a mouse-size hot-rod resolution phantom and two in vivo simultaneous PET/MRI scans in a 7 T MRI scanner. The insert had a peak sensitivity of 0.49%, using an energy threshold of 350 keV. A uniform transaxial resolution was obtained up to 15 mm radial offset from the axial center, using filtered back-projection with single-slice rebinning. The measured average radial and tangential resolutions (FWHM) were 1.38 mm and 1.39 mm, respectively. The 1.2 mm rods were separable in the hot-rod phantom using an iterative image reconstruction algorithm. The scatter fraction was 7.3% and peak noise equivalent count rate was 15.5 kcps at 65.1 MBq of activity. The FDG uptake in a mouse heart and brain were visible in the two in vivo simultaneous PET/MRI scans without applying image corrections. In conclusion, the insert demonstrated a good overall performance and can be used for small animal multi-modal research applications.

Open-field mouse brain PET: design optimisation and detector characterisation

'Open-field' PET, in which an animal is free to move within an enclosed space during imaging, is a very promising advance for neuroscientific research. It provides a key advantage over conventional imaging under anesthesia by enabling functional changes in the brain to be correlated with an animal's behavioural response to environmental or pharmacologic stimuli. Previously we have demonstrated the feasibility of open-field imaging of rats using motion compensation techniques applied to a commercially available PET scanner. However, this approach of 'retro-fitting' motion compensation techniques to an existing system is limited by the inherent geometric and performance constraints of the system. The goal of this project is to develop a purpose-built PET scanner with geometry, motion tracking and imaging performance tailored and optimised for open-field imaging of the mouse brain. The design concept is a rail-based sliding tomograph which moves according to the animal's motion. Our specific aim in this work was to evaluate candidate scanner designs and characterise the performance of a depth-of-interaction detector module for the open-field system. We performed Monte Carlo simulations to estimate and compare the sensitivity and spatial resolution performance of four scanner geometries: a ring, parallel plate, and two box variants. Each system was based on a detector block consisting of a 23  ×  23 array of 0.785  ×  0.785  ×  20 mm³ LSO crystals (overall dim. 19.6  ×  19.6  ×  20 mm). We found that a DoI resolution capability of 3 mm was necessary to achieve approximately uniform sub-millimetre spatial resolution throughout the FoV for all scanners except the parallel-plate geometry. With this DoI performance, the sensitivity advantage afforded by the box geometry with overlapping panels (16% peak absolute sensitivity, a 36% improvement over the ring design) suggests this unconventional design is best suited for imaging the mouse brain. We also built and characterised the block detector modelled in the simulations, including a dual-ended readout based on 6  ×  6 arrays of through-silicon-via silicon photomultipliers (active area 84%) for DoI estimation. Identification of individual crystals in the flood map was excellent, energy resolution varied from 12.4%  ±  0.6% near the centre to 24.4%  ±  3.4% at the ends of the crystal, and the average DoI resolution was 2.8 mm  ±  0.35 mm near the central depth (10 mm) and 3.5 mm  ±  1.0 mm near the ends. Timing resolution was 1.4  ±  0.14 ns. Therefore, the DoI detector module meets the target specifications for the application and will be used as the basis for a prototype open-field mouse PET scanner.

Preclinical positron emission tomography scanner based on a monolithic annulus of scintillator: initial design study

Positron emission tomography (PET) scanners designed for imaging of small animals have transformed translational research by reducing the necessity to invasively monitor physiology and disease progression. Virtually all of these scanners are based on the use of pixelated detector modules arranged in rings. This design, while generally successful, has some limitations. Specifically, use of discrete detector modules to construct PET scanners reduces detection sensitivity and can introduce artifacts in reconstructed images, requiring the use of correction methods. To address these challenges, and facilitate measurement of photon depth-of-interaction in the detector, we investigated a small animal PET scanner (called AnnPET) based on a monolithic annulus of scintillator. The scanner was created by placing 12 flat facets around the outer surface of the scintillator to accommodate placement of silicon photomultiplier arrays. Its performance characteristics were explored using Monte Carlo simulations and sections of the NEMA NU4-2008 protocol. Results from this study revealed that AnnPET's reconstructed spatial resolution is predicted to be [Formula: see text] full width at half maximum in the radial, tangential, and axial directions. Peak detection sensitivity is predicted to be 10.1%. Images of simulated phantoms (mini-hot rod and mouse whole body) yielded promising results, indicating the potential of this system for enhancing PET imaging of small animals.

Evaluation of a BGO-Based PET System for Single-Cell Tracking Performance by Simulation and Phantom Studies

A recent method based on positron emission was reported for tracking moving point sources using the Inveon PET system. However, the effect of scanner background noise was not further explored. Here, we evaluate tracking with the Genisys4, a bismuth germanate-based PET system, which has no significant intrinsic background and may be better suited to tracking lower and/or faster activity sources. Position-dependent sensitivity of the Genisys4 was simulated in Geant4 Application for Tomographic Emission (GATE) using a static ¹⁸F point source. Trajectories of helically moving point sources with varying activity and rotation speed were reconstructed from list-mode data as described previously. Simulations showed that the Inveon’s ability to track sources within 2 mm of localization error is limited to objects with a velocity-to-activity ratio < 0.13 mm/decay, compared to < 0.29 mm/decay for the Genisys4. Tracking with the Genisys4 was then validated using a physical phantom of helically moving [¹⁸F] fluorodeoxyglucose-in-oil droplets (< 0.24 mm diameter, 139-296 Bq), yielding < 1 mm localization error under the tested conditions, with good agreement between simulated sensitivity and measured activity (Pearson correlation R = .64, P << .05 in a representative example). We have investigated the tracking performance with the Genisys4, and results suggest the feasibility of tracking low activity, point source-like objects with this system.

A New Pulse Pileup Rejection Method Based on Position Shift Identification

Pulse pileup events degrade the signal-to-noise ratio (SNR) of nuclear medicine data. When such events occur in multiplexed detectors, they cause spatial misposition, energy spectrum distortion and degraded timing resolution, which leads to image artifacts. Pulse pileup is pronounced in PETbox4, a bench top PET scanner dedicated to high sensitivity and high resolution imaging of mice. In that system, the combination of high absolute sensitivity, long scintillator decay time (BGO) and highly multiplexed electronics lead to a significant fraction of pulse pileup, reached at lower total activity than for comparable instruments. In this manuscript, a new pulse pileup rejection method named position shift rejection (PSR) is introduced. The performance of PSR is compared with a conventional leading edge rejection (LER) method and with no pileup rejection implemented (NoPR). A comprehensive digital pulse library was developed for objective evaluation and optimization of the PSR and LER, in which pulse waveforms were directly recorded from real measurements exactly representing the signals to be processed. Physical measurements including singles event acquisition, peak system sensitivity and NEMA NU-4 image quality phantom were also performed in the PETbox4 system to validate and compare the different pulse pile-up rejection methods. The evaluation of both physical measurements and model pulse trains demonstrated that the new PSR performs more accurate pileup event identification and avoids erroneous rejection of valid events. For the PETbox4 system, this improvement leads to a significant recovery of sensitivity at low count rates, amounting to about 1/4th of the expected true coincidence events, compared to the LER method. Furthermore, with the implementation of PSR, optimal image quality can be achieved near the peak noise equivalent count rate (NECR).

TandemPET- A High Resolution, Small Animal, Virtual Pinhole-Based PET Scanner: Initial Design Study

Mice are the perhaps the most common species of rodents used in biomedical research, but many of the current generation of small animal PET scanners are non-optimal for imaging these small rodents due to their relatively low resolution. Consequently, a number of researchers have investigated the development of high-resolution scanners to address this need. In this investigation, the design of a novel, high-resolution system based on the dual-detector, virtual-pinhole PET concept was explored via Monte Carlo simulations. Specifically, this system, called TandemPET, consists of a 5 cm × 5 cm high-resolution detector made-up of a 90 × 90 array of 0.5 mm × 0.5 mm × 10 mm (pitch= 0.55 mm) LYSO detector elements in coincidence with a lower resolution detector consisting of a 68 × 68 array of 1.5 mm × 1.5 mm × 10 mm LYSO detector elements (total size= 10.5 cm × 10.5 cm). Analyses indicated that TandemPET's optimal geometry is to position the high-resolution detector 3 cm from the center-of-rotation, with the lower resolution detector positioned 9 cm from center. Measurements using modified NEMA NU4-2008-based protocols revealed that the spatial resolution of the system is ~0.5 mm FWHM, after correction of positron range effects. Peak sensitivity is 2.1%, which is comparable to current small animal PET scanners. Images from a digital mouse brain phantom demonstrated the potential of the system for identifying important neurological structures.

Performance evaluation of the small-animal PET scanner ClairvivoPET using NEMA NU 4-2008 Standards

The aim of this study was to evaluate the performance of ClairvivoPET using NEMA NU4 standards. The ClairvivoPET incorporates a LYSO dual depth-of-interaction detector system with 151 mm axial field of view (FOV). Spatial resolution, sensitivity, counting rate capabilities, and image quality were evaluated using NEMA NU4-2008 standards. Normal mouse imaging was also performed for 10 min after intravenous injection of (18)F(-)-NaF. Data were compared with 19 other preclinical PET scanners. Spatial resolution measured using full width at half maximum on FBP-ramp reconstructed images was 2.16 mm at radial offset 5 mm of the axial centre FOV. The maximum absolute sensitivity for a point source at the FOV centre was 8.72%. Peak noise equivalent counting rate (NECR) was 415 kcps at 14.6 MBq ml(-1). The uniformity with the image-quality phantom was 4.62%. Spillover ratios in the images of air and water filled chambers were 0.19 and 0.06, respectively. Our results were comparable with the 19 other preclinical PET scanners based on NEMA NU4 standards, with excellent sensitivity because of the large FOV. The ClairvivoPET with iterative reconstruction algorithm also provided sufficient visualization of the mouse spine. The high sensitivity and resolution of the ClairvivoPET scanner provided high quality images for preclinical studies.

Ultrasound-aided Multi-parametric Photoacoustic Microscopy of the Mouse Brain

High-resolution quantitative imaging of cerebral oxygen metabolism in mice is crucial for understanding brain functions and formulating new strategies to treat neurological disorders, but remains a challenge. Here, we report on our newly developed ultrasound-aided multi-parametric photoacoustic microscopy (PAM), which enables simultaneous quantification of the total concentration of hemoglobin (CHb), the oxygen saturation of hemoglobin (sO2), and cerebral blood flow (CBF) at the microscopic level and through the intact mouse skull. The three-dimensional skull and vascular anatomies delineated by the dual-contrast (i.e., ultrasonic and photoacoustic) system provide important guidance for dynamically focused contour scan and vessel orientation-dependent correction of CBF, respectively. Moreover, bi-directional raster scan allows determining the direction of blood flow in individual vessels. Capable of imaging all three hemodynamic parameters at the same spatiotemporal scale, our ultrasound-aided PAM fills a critical gap in preclinical neuroimaging and lays the foundation for high-resolution mapping of the cerebral metabolic rate of oxygen (CMRO2)-a quantitative index of cerebral oxygen metabolism. This technical innovation is expected to shed new light on the mechanism and treatment of a broad spectrum of neurological disorders, including Alzheimer's disease and ischemic stroke.

A deformable atlas of the laboratory mouse

PURPOSE

This paper presents a deformable mouse atlas of the laboratory mouse anatomy. This atlas is fully articulated and can be positioned into arbitrary body poses. The atlas can also adapt body weight by changing body length and fat amount.

PROCEDURES

A training set of 103 micro-CT images was used to construct the atlas. A cage-based deformation method was applied to realize the articulated pose change. The weight-related body deformation was learned from the training set using a linear regression method. A conditional Gaussian model and thin-plate spline mapping were used to deform the internal organs following the changes of pose and weight.

RESULTS

The atlas was deformed into different body poses and weights, and the deformation results were more realistic compared to the results achieved with other mouse atlases. The organ weights of this atlas matched well with the measurements of real mouse organ weights. This atlas can also be converted into voxelized images with labeled organs, pseudo CT images and tetrahedral mesh for phantom studies.

CONCLUSIONS

With the unique ability of articulated pose and weight changes, the deformable laboratory mouse atlas can become a valuable tool for preclinical image analysis.

Evaluation of Matrix9 silicon photomultiplier array for small-animal PET

PURPOSE

The MatrixSL-9-30035-OEM (Matrix9) from SensL is a large-area silicon photomultiplier (SiPM) photodetector module consisting of a 3 × 3 array of 4 × 4 element SiPM arrays (total of 144 SiPM pixels) and incorporates SensL's front-end electronics board and coincidence board. Each SiPM pixel measures 3.16 × 3.16 mm(2) and the total size of the detector head is 47.8 × 46.3 mm(2). Using 8 × 8 polished LSO/LYSO arrays (pitch 1.5 mm) the performance of this detector system (SiPM array and readout electronics) was evaluated with a view for its eventual use in small-animal positron emission tomography (PET).

METHODS

Measurements of noise, signal, signal-to-noise ratio, energy resolution, flood histogram quality, timing resolution, and array trigger error were obtained at different bias voltages (28.0-32.5 V in 0.5 V intervals) and at different temperatures (5 °C-25 °C in 5 °C degree steps) to find the optimal operating conditions.

RESULTS

The best measured signal-to-noise ratio and flood histogram quality for 511 keV gamma photons were obtained at a bias voltage of 30.0 V and a temperature of 5 °C. The energy resolution and timing resolution under these conditions were 14.2% ± 0.1% and 4.2 ± 0.1 ns, respectively. The flood histograms show that all the crystals in the 1.5 mm pitch LSO array can be clearly identified and that smaller crystal pitches can also be resolved. Flood histogram quality was also calculated using different center of gravity based positioning algorithms. Improved and more robust results were achieved using the local 9 pixels for positioning along with an energy offset calibration. To evaluate the front-end detector readout, and multiplexing efficiency, an array trigger error metric is introduced and measured at different lower energy thresholds. Using a lower energy threshold greater than 150 keV effectively eliminates any mispositioning between SiPM arrays.

CONCLUSIONS

In summary, the Matrix9 detector system can resolve high-resolution scintillator arrays common in small-animal PET with adequate energy resolution and timing resolution over a large detector area. The modular design of the Matrix9 detector allows it to be used as a building block for simple, low channel-count, yet high performance, small animal PET or PET/MRI systems.

A DOI Detector With Crystal Scatter Identification Capability for High Sensitivity and High Spatial Resolution PET Imaging

A new phoswich detector is being developed at the Crump Institute, aiming to provide improvements in sensitivity, and spatial resolution for PET. The detector configuration is comprised of two layers of pixelated scintillator crystal arrays, a glass light guide and a light detector. The annihilation photon entrance (top) layer is a 48 × 48 array of 1.01 × 1.01 × 7 mm³ LYSO crystals. The bottom layer is a 32 × 32 array of 1.55 × 1.55 × 9 mm³ BGO crystals. A tapered, multiple-element glass lightguide is used to couple the exit end of the BGO crystal array (52 × 52 mm²) to the photosensitive area of the Position Sensitive Photomultiplier Tube (46 × 46 mm²), allowing the creation of flat panel detectors without gaps between the detector modules. Both simulations and measurements were performed to evaluate the characteristics and benefits of the proposed design. The GATE Monte Carlo simulation indicated that the total fraction of the cross layer crystal scatter (CLCS) events in singles detection mode for this detector geometry is 13.2%. The large majority of these CLCS events (10.1% out of 13.2%) deposit most of their energy in a scintillator layer other than the layer of first interaction. Identification of those CLCS events for rejection or correction may lead to improvements in data quality and imaging performance. Physical measurements with the prototype detector showed that the LYSO, BGO and CLCS events were successfully identified using the delayed charge integration (DCI) technique, with more than 95% of the LYSO and BGO crystal elements clearly resolved. The measured peak-to-valley ratios (PVR) in the flood histograms were 3.5 for LYSO and 2.0 for BGO. For LYSO, the energy resolution ranged from 9.7% to 37.0% full width at half maximum (FWHM), with a mean of 13.4 ± 4.8%. For BGO the energy resolution ranged from 16.0% to 33.9% FWHM, with a mean of 18.6 ± 3.2%. In conclusion, these results demonstrate that the proposed detector is feasible and can potentially lead to a high spatial resolution, high sensitivity and DOI PET system.

DigiPET: sub-millimeter spatial resolution small-animal PET imaging using thin monolithic scintillators

A new preclinical PET system based on dSiPMs, called DigiPET, is presented. The system is based on thin monolithic scintillation crystals and exhibits superior spatial resolution at low-cost compared to systems based on pixelated crystals. Current dedicated small-rodent PET scanners have a spatial resolution in the order of 1 mm. Most of them have a large footprint, requiring considerable laboratory space. For rodent brain imaging, a PET scanner with sub-millimeter resolution is desired. To achieve this, crystals with a pixel pitch down to 0.5 mm have been used. However, fine pixels are difficult to produce and will render systems expensive. In this work, we present the first results with a high-resolution preclinical PET scanner based on thin monolithic scintillators and a large solid angle. The design is dedicated to rat-brain imaging and therefore has a very compact geometry. Four detectors were placed in a square arrangement with a distance of 34.5 mm between two opposing detector modules, defining a field of view (FOV) of 32 × 32 × 32 mm(3). Each detector consists of a thin monolithic LYSO crystal of 32 × 32 × 2 mm(3) optically coupled to a digital silicon photomultiplier (dSiPM). Event positioning within each detector was obtained using the maximum likelihood estimation (MLE) method. To evaluate the system performance, we measured the energy resolution, coincidence resolving time (CRT), sensitivity and spatial resolution. The image quality was evaluated by acquiring a hot-rod phantom filled with (18)F-FDG and a rat head one hour after an (18)F-FDG injection. The MLE yielded an average intrinsic spatial resolution on the detector of 0.54 mm FWHM. We obtained a CRT of 680 ps and an energy resolution of 18% FWHM at 511 keV. The sensitivity and spatial resolution obtained at the center of the FOV were 6.0 cps kBq(-1) and 0.7 mm, respectively. In the reconstructed images of the hot-rod phantom, hot rods down to 0.7 mm can be discriminated. In conclusion, a compact PET scanner was built using dSiPM technology and thin monolithic LYSO crystals. Excellent spatial resolution and acceptable sensitivity were demonstrated. Promising results were also obtained in a hot-rod phantom and in rat-brain imaging.

Optimization of the Energy Window for PETbox4, a Preclinical PET Tomograph With a Small Inner Diameter

Small animal positron emission tomography (PET) systems are often designed by employing close geometry configurations. Due to the different characteristics caused by geometrical factors, these tomographs require data acquisition protocols that differ from those optimized for conventional large diameter ring systems. In this work we optimized the energy window for data acquisitions with PETbox4, a 50 mm detector separation (box-like geometry) pre-clinical PET scanner, using the Geant4 Application for Tomographic Emission (GATE). The fractions of different types of events were estimated using a voxelized phantom including a mouse as well as its supporting chamber, mimicking a realistic mouse imaging environment. Separate code was developed to extract additional information about the gamma interactions for more accurate event type classification. Three types of detector backscatter events were identified in addition to the trues, phantom scatters and randoms. The energy window was optimized based on the noise equivalent count rate (NECR) and scatter fraction (SF) with lower-level discriminators (LLD) corresponding to energies from 150 keV to 450 keV. The results were validated based on the calculated image uniformity, spillover ratio (SOR) and recovery coefficient (RC) from physical measurements using the National Electrical Manufacturers Association (NEMA) NU-4 image quality phantom. These results indicate that when PETbox4 is operated with a more narrow energy window (350-650 keV), detector backscatter rejection is unnecessary. For the NEMA NU-4 image quality phantom, the SOR for the water chamber decreases by about 45% from 15.1% to 8.3%, and the SOR for the air chamber decreases by 31% from 12.0% to 8.3% at the LLDs of 150 and 350 keV, without obvious change in uniformity, further supporting the simulation based optimization. The optimization described in this work is not limited to PETbox4, but also applicable or helpful to other small inner diameter geometry scanners.

Cold wall effect eliminating method to determine the contrast recovery coefficient for small animal PET scanners using the NEMA NU-4 image quality phantom

The contrast recovery coefficients (CRC) were evaluated for five different small animal PET scanners: GE Explore Vista, Genisys4, MiniPET-2, nanoScan PC and Siemens Inveon. The NEMA NU-4 2008 performance test with the suggested image quality phantom (NU4IQ) does not allow the determination of the CRC values for the hot regions in the phantom. This drawback of NU4IQ phantom motivated us to develop a new method for this purpose. The method includes special acquisition and reconstruction protocols using the original phantom, and results in an artificially merged image enabling the evaluation of CRC values. An advantageous feature of this method is that it stops the cold wall effect from distorting the CRC calculation. Our suggested protocol results in a set of CRC values contributing to the characterization of small animal PET scanners. GATE simulations were also performed to validate the new method and verify the evaluated CRC values. We also demonstrated that the numerical values of this parameter depend on the actual object contrast of the hot region(s) and this mainly comes from the spillover effect. This effect was also studied while analysing the background activity level around the hot rods. We revealed that the calculated background mean values depended on the target contrast in a scanner specific manner. Performing the artificially merged imaging procedure and additional simulations using the micro hollow sphere (MHS) phantom geometry, we also proved that the inactive wall around the hot spheres can have a remarkable impact on the calculated CRC. In conclusion, we have shown that the proposed artificial merging procedure and the commonly used NU4IQ phantom prescribed by the NEMA NU-4 can easily deliver reliable CRC data otherwise unavailable for the NU4IQ phantom in the conventional protocol or the MHS phantom.

Accurate and efficient modeling of the detector response in small animal multi-head PET systems

In fully three-dimensional PET imaging, iterative image reconstruction techniques usually outperform analytical algorithms in terms of image quality provided that an appropriate system model is used. In this study we concentrate on the calculation of an accurate system model for the YAP-(S)PET II small animal scanner, with the aim to obtain fully resolution- and contrast-recovered images at low levels of image roughness. For this purpose we calculate the system model by decomposing it into a product of five matrices: (1) a detector response component obtained via Monte Carlo simulations, (2) a geometric component which describes the scanner geometry and which is calculated via a multi-ray method, (3) a detector normalization component derived from the acquisition of a planar source, (4) a photon attenuation component calculated from x-ray computed tomography data, and finally, (5) a positron range component is formally included. This system model factorization allows the optimization of each component in terms of computation time, storage requirements and accuracy. The main contribution of this work is a new, efficient way to calculate the detector response component for rotating, planar detectors, that consists of a GEANT4 based simulation of a subset of lines of flight (LOFs) for a single detector head whereas the missing LOFs are obtained by using intrinsic detector symmetries. Additionally, we introduce and analyze a probability threshold for matrix elements of the detector component to optimize the trade-off between the matrix size in terms of non-zero elements and the resulting quality of the reconstructed images. In order to evaluate our proposed system model we reconstructed various images of objects, acquired according to the NEMA NU 4-2008 standard, and we compared them to the images reconstructed with two other system models: a model that does not include any detector response component and a model that approximates analytically the depth of interaction as detector response component. The comparisons confirm previous research results, showing that the usage of an accurate system model with a realistic detector response leads to reconstructed images with better resolution and contrast recovery at low levels of image roughness.