Abstract
The Inveon dedicated PET (DPET) scanner is the latest generation of preclinical PET systems devoted to high-resolution and high-sensitivity murine model imaging. In this study, we report on its performance based on the National Electrical Manufacturers Association (NEMA) NU-4 standards. Methods: The Inveon DPET consists of 64 lutetium oxyorthosilicate block detectors arranged in 4 contiguous rings, with a 16.1-cm ring diameter and a 12.7-cm axial length. Each detector block consists of a 20 × 20 lutetium oxyorthosilicate crystal array of 1.51 × 1.51 × 10.0 mm elements. The scintillation light is transmitted to position-sensitive photomultiplier tubes via optical light guides. Energy resolution, spatial resolution, sensitivity, scatter fraction, and counting-rate performance were evaluated. The NEMA NU-4 image–quality phantom and a healthy mouse injected with 18F-FDG and 18F− were scanned to evaluate the imaging capability of the Inveon DPET. Results: The energy resolution at 511 keV was 14.6% on average for the entire system. In-plane radial and tangential resolutions reconstructed with Fourier rebinning and filtered backprojection algorithms were below 1.8-mm full width at half maximum (FWHM) at the center of the field of view. The radial and tangential resolution remained under 2.0 mm, and the axial resolution remained under 2.5-mm FWHM within the central 4-cm diameter of the field of view. The absolute sensitivity of the system was 9.3% for an energy window of 250–625 keV and a timing window of 3.432 ns. At a 350- to 625-keV energy window and a 3.432-ns timing window, the peak noise equivalent counting rate was 1,670 kcps at 130 MBq for the mouse-sized phantom and 590 kcps at 110 MBq for the rat-sized phantom. The scatter fractions at the same acquisition settings were 7.8% and 17.2% for the mouse- and rat-sized phantoms, respectively. The mouse image-quality phantom results demonstrate that for typical mouse acquisitions, the image quality correlates well with the measured performance parameters in terms of image uniformity, recovery coefficients, attenuation, and scatter corrections. Conclusion: The Inveon system, compared with previous generations of preclinical PET systems from the same manufacturer, shows significantly improved energy resolution, sensitivity, axial coverage, and counting-rate capabilities. The performance of the Inveon is suitable for successful murine model imaging experiments.
Tomographic systems dedicated to noninvasive, in vivo imaging of preclinical animal models have been widely used at research institutes in recent years (1,2). With the ability to longitudinally image the same subject, each individual animal can serve as its own control. Therefore, intersubject variability can be minimized. Because of the dramatic difference in size between humans and rodents, small-animal PET imposes higher performance requirements than does clinical PET, particularly on image resolution and sensitivity. The resolution and sensitivity improvements are mainly achieved by using smaller crystal sizes, smaller detector ring diameters, and longer axial coverage. With the goal of improving the image quality of laboratory small-animal studies, researchers continually develop new techniques, including new scintillation materials, detectors, electronics, and geometries, to improve the performance of preclinical tomographs.
The Inveon dedicated PET (DPET) is the latest generation of commercial tomographs from Siemens Preclinical Solutions, Inc. DPET incorporates changes in system geometry, detectors, and electronics. This work evaluates the performance of the Inveon system in all aspects, including energy and spatial resolutions, sensitivity, scatter fraction, counting-rate performance, and imaging capabilities, on the basis of the National Electrical Manufacturers Association (NEMA) NU-4 standards for performance evaluation of small-animal PET scanners (3). The NEMA NU-4 standards offer a uniform methodology for small-animal PET performance evaluation and establish a baseline of system performance in typical imaging conditions.
MATERIALS AND METHODS
System Description
The Inveon DPET is a lutetium oxyorthosilicate (LSO)–based, high-sensitivity, high-resolution preclinical PET scanner used primarily for murine model imaging. The system consists of 64 detector blocks arranged in 4 contiguous rings, with a crystal ring diameter of 16.1 cm and an axial extent of 12.7 cm. Each detector block is composed of a 20 × 20 array of LSO crystals coupled to a position-sensitive photomultiplier tube via a light guide. Each crystal is 10.0 mm long and has a cross-sectional area of 1.51 × 1.51 mm. The crystal pitch is 1.59 mm in both axial and transverse directions.
List-mode data are acquired during measurements. From the list-mode data, coincidence events can then be sorted into 3-dimensional (3D) sinograms, with different combinations of span and ring differences, or into 2-dimensional (2D) sinograms by either single-slice rebinning (SSRB) (4) or Fourier rebinning (FORE) (5). Images can be reconstructed using analytic 2D filtered backprojection (FBP), 3D reprojection (6), or iterative methods, such as ordered-subsets expectation maximization (7,8) and maximum a posteriori (MAP) (9).
The specifications between the Inveon DPET and 3 previous preclinical systems from the same manufacturer are compared in Table 1.
Comparison of Small-Animal PET Systems
Energy Resolution
An 18F point source was placed at the center of the field of view (FOV) to acquire 2D position histograms of each detector in singles mode, with the energy window wide open. A total of 100,000 counts were acquired for each detector pixel. Lookup tables were generated for individual crystal identification (10).
Energy resolution was determined for each crystal in the system and calculated as full width at half maximum (FWHM) of the 511-keV energy peak divided by the center of the photopeak value. The mean value of all crystal energy resolutions was calculated, and the maximum and minimum energy resolutions were obtained.
Spatial Resolution
Spatial resolution was measured with a 22Na point source conforming to the NEMA NU-4 standards. The 22Na point source has a nominal size (0.3 mm), is embedded in an acrylic cube (10.0-mm extent on all sides), and has a nominal activity (198,000 Bq). The energy window setting was 350–625 keV, and the timing window was 3.432 ns (the default of the 4 available timing window settings on the Inveon system). 22Na has emission energy (Eavgβ+ = 250 keV) and positron range (∼0.23 mm) similar to 18F, which is the most widely used positron-emitting isotope. Per the NEMA protocol, the measured spatial resolutions were not corrected for source size, positron range, or photon acolinearity.
The source was fixed in the tomograph and located at 2 axial positions (the center of the axial FOV and one fourth of the axial FOV [31.75 mm from the center, along the axial direction]). For each of the 2 axial positions, the source was stepped toward the edge of the transverse FOV. For the central 5 mm of the transverse FOV, the source was stepped at 1-mm increments and then at 5-mm steps up to the edge of the FOV.
The list-mode data acquired at each location were histogrammed into 3D sinograms with delayed events subtracted to correct for random coincidences. Component-based normalization was applied to compensate for the differences in detection efficiency (11). The 3D sinograms were first Fourier rebinned into 2D sinograms and then reconstructed by 2D FBP with a ramp filter cut off at the Nyquist frequency, with a zoom selected to achieve a 0.4-mm-pixel in-plane resolution. The axial plane separation was 0.796 mm. The response function was formed by summing 1-dimensional profiles that were parallel to the radial, tangential, and axial directions. A parabolic fit of the peak point and its 2 nearest neighboring points was used to determine the maximum value of the response function. Linear interpolation between adjacent pixels was used to determine the position of half and one tenth of the maximum; the FWHM and full width at tenth maximum were determined for each extracted profile. Volumetric resolution was calculated on the basis of the FWHMs of the radial, tangential, and axial directions.
Sensitivity
The sensitivity of the system was measured with the same 22Na point source used in the spatial resolution measurement and an 18F point source. In addition, the sensitivity was also measured with an 18F line source inserted in a set of concentric aluminum sleeves, which was a more traditional methodology used in the past (12). To reduce the attenuation from the imaging bed, the 18F and 22Na point sources were taped on a thin piece of cardboard and placed into the scanner FOV. The concentric aluminum tubes were suspended on both ends without other attenuation material in the FOV.
Because of the difficulty of accurately measuring the small nominal activity of a 22Na source in a standard dose calibrator, the 22Na point source was used only to determine the relative sensitivity between different energy and timing window settings. The 22Na source was positioned at the center of the FOV and scanned for 5 min at a fixed timing window of 3.432 ns. Two sets of energy windows were used: one with a fixed lower-level discriminator (LLD) that equaled 350 keV and a changing upper-level discriminator (ULD) from 550 to 700 keV and the other with a fixed ULD that equaled 625 keV and a changing LLD from 250 to 450 keV, both at 50-keV steps. The reason for using a solid 22Na source instead of liquid 18F to measure the relative sensitivity was 2-fold: to simplify the experimental protocol and to avoid the complications of decay correction between different measurements.
The absolute sensitivity was determined by an 18F point source of approximately 40 μL and 42,550 Bq placed at the tip of a small centrifuge tube and positioned at the center of the scanner FOV. The absolute sensitivity was measured for 5 min at an energy window of 350–625 keV and a timing window of 3.432 ns. The activity of the 18F source was measured in a γ-counter (Wallac 1480 Wizard 3″; PerkinElmer Life Sciences).
The LSO scintillator crystals have intrinsic radioactivity (13–15). 176Lu emits β− particles with an average energy of 420 keV, together with 3 γ-photons of 307, 202, and 88 keV, respectively (16). The β− particles and γ-photons can make a true coincidence if both fall in the preset energy window (12,17).
For each energy window, a background measurement was acquired for 5 min. The histogrammed background true counts were subtracted from the total histogrammed true counts when the 22Na or 18F point sources were placed inside the scanner. The number of net true coincidences was normalized to the scan duration, divided by the source activity, and corrected for the branching ratio (0.906 for 22Na and 0.967 for 18F). Attenuations of the 1-cm3 cube of 22Na and the centrifuge tube for the 18F source were not compensated, but a 9% sensitivity loss was estimated on the basis of the 0.095 cm−1 attenuation of 511-keV photons in water-equivalent material. The absolute sensitivity at the 350- to 625-keV energy window was determined by the 18F measurement. On the basis of the relative sensitivity determined by the 22Na source, the absolute sensitivity for the 350-keV LLD and 625-keV ULD energy window datasets was also calculated.
We used these measurements to investigate the system sensitivity as a function of the energy window. The energy window used in typical studies should be determined on the basis of the tradeoff between absolute sensitivity, scatter fraction, and system background. For our institute, an energy window of 350–625 keV was selected as a compromise for the typical studies we perform.
The sensitivity dependence on the timing window was also measured for 3 min with the same 22Na point source at the 4 available timing window settings (2.808, 3.432, 4.056, and 4.680 ns) and a fixed energy window of 350–625 keV.
The axial sensitivity profile was measured with a set of concentric aluminum tubes and a plastic tube filled with 18F solution, with both ends sealed (12). The source tubing was placed inside the smallest metal tube, suspended in the center of the transaxial FOV, and aligned with the axis of the tomograph. The other aluminum tubes were added, one at a time, and the counting rate was measured for each set of metal tubes for 120 s at the 350- to 625-keV energy window and 3.432-ns timing window. A set of measurements with different thicknesses of aluminum tubes was used to obtain the counting rate without attenuation by exponential fitting. The list-mode data were histogrammed with SSRB without randoms or scatter correction. For each slice of each acquisition and for points farther than 1 cm from each side of the peak, the values were set to zero. The counts in all lines of response in the sinograms were summed slice by slice, corrected for decay, and scaled by the acquisition time and the 96.7% positron yield for 18F. The 18F activity was measured in a Wallac γ-counter, and the counts in the central 7-cm and the whole axial FOV were added to calculate the sensitivity of mouse- and rat-sized objects, respectively.
Scatter Fraction and Counting-Rate Performance
Scatter fraction and counting-rate performance were measured using 2 different cylindric polyethylene phantoms that simulated the geometries of a mouse and rat. The design of the phantoms conformed to the NEMA NU-4 standards.
Both phantoms were made of high-density polyethylene (0.96 g/cm3). The mouselike phantom was a 70-mm-long solid cylinder with a 25-mm diameter. A cylindric hole (diameter, 3.2 mm) was drilled parallel to the central axis, at a radial distance of 10 mm. The ratlike phantom had similar geometry but larger dimensions (length, 150 mm; diameter, 50-mm). A 3.2-mm-diameter hole was drilled at a radial offset of 17.5 mm.
A 11C solution with a concentration higher than 1,500 MBq/mL was enclosed in a flexible tube with an outer diameter fitting the 3.2-mm hole. The initial activity in the FOV was measured in a dose calibrator (Atomlab 300; Biodex Medical Systems) and was higher than 500 and 600 MBq, respectively, at the start of the acquisition for the mouse- and rat-sized phantoms. The acquisition was performed at the 350- to 625-keV energy and 3.432-ns timing windows. The phantom was centered in the FOV, and data were acquired until the total activity decayed below 10,000 Bq. The random coincidences were measured by the delayed-window technique. The list-mode data were histogrammed into 2D sinogram sets lasting 5 min by SSRB, with separate prompts and delays.
For each prompt sinogram (transaxial bin size, 0.815 mm; slice thickness, 0.796 mm), all pixels located farther than 8 mm from the edge of the phantom were set to zero. The profile of each projection angle was shifted so that the peak pixels were aligned with the center pixel of the sinogram. A sum projection was then produced by adding 160 angular projections in each slice and each frame. All pixel counts outside a 14-mm centered band were assumed to be the sum of random, scatter, and intrinsic counts. A linear interpolation between the left and right border of the 14-mm band was used to estimate these nontrue counts under the profile peak. Counts above this line were regarded as true coincidences. Random coincidences were estimated from the delayed sinogram. A background acquisition was obtained at the same energy window and timing window for 16 h with the cold mouse- and rat-sized phantoms in the FOV. The intrinsic counts were estimated from the background sinogram.
The scattered counting rate was then calculated by Equation 1:Eq. 1where Rscatter, Rtotal, Rtrue, Rrandom, and Rintrinsic are the scatter, total, true, random, and intrinsic counting rates, respectively. The scatter fraction (SF) was calculated by Equation 2:
Eq. 2
The noise equivalent counting rate (NECR) of each of the 5-min frame acquisitions was determined using the following equation (18,19):Eq. 3
Imaging Studies
NEMA Phantom Study.
The NEMA NU-4 mouse image-quality phantom (3) is composed of the following 3 regions: a main fillable uniform region chamber (diameter, 30 mm; length, 30 mm); a lid that attaches to the main fillable region, containing 2 smaller cold region chambers—one that was filled with nonradioactive water and the other with air (both chambers were composed of hollow cylinders—length, 15 mm; outer diameter [OD], 10 mm; and wall thickness, 1 mm); and a solid acrylic glass region (diameter, 30 mm; length, 20 mm), with 5 fillable rods drilled through (at 7 mm from the center) with diameters of 1, 2, 3, 4, and 5 mm, respectively. The image-quality phantom was filled with 18F-FDG solution (5.1 MBq) and acquired at 350- to 625-keV energy and 3.432-ns timing windows for 20 min. The phantom was placed on a mouse imaging chamber to simulate actual mouse imaging. The activity in the phantom was measured with a Wallac γ-counter.
A CT transmission scan of the image-quality phantom and the imaging chamber was obtained using a MicroCAT II tomograph (Siemens Preclinical Solutions). The angular sampling was 1° per projection for a full 360° scan. The radiographic source was operated at 500 μA and 70 kVp. Images were reconstructed using a Feldkamp cone-beam algorithm, with a ramp filter cut off at the Nyquist frequency. The reconstructed CT image was registered with the PET emission image to create an attenuation sinogram. The whole 3D PET sinogram data were first Fourier rebinned and then reconstructed by 2D FBP, with a ramp filter cut off at the Nyquist frequency. Normalization, dead time, random, attenuation (20), and scatter corrections (21,22) were applied.
A 22.5-mm-diameter and a 10-mm-high cylindric volume of interest (VOI) was drawn over the center of the uniform region of the image-quality phantom. The average concentration, maximum and minimum values in this VOI, and percentage SD were measured to estimate the noise performance as a measure of uniformity.
The image slices covering the central 10-mm length of the rods were averaged to obtain a single image slice of lower noise. Circular regions of interest (ROIs) were drawn in this image around each rod, with diameters twice the physical diameters of the rods (the exact size of the ROIs is not critical for this measurement). The maximum values in each of these ROIs were measured and were divided by the mean value obtained in the uniformity test to obtain the recovery coefficient (RC) for each rod size.
The transverse image pixel coordinates of the locations with the maximum ROI values were recorded and used to create 10-mm-long line profiles along the rods in the axial direction. The SD of the pixel values measured along each line profile was calculated. The SD of the RC was calculated as follows:Eq. 4
To test for the accuracy of scatter correction, VOIs were defined in the water- and air-filled cylindric inserts. The diameter of the VOI was 4 mm and encompassed the central 7.5 mm in length in the axial direction. The ratio of the mean in each cold region to the mean of the hot uniform area was reported as spillover ratio.
Mouse Study.
One healthy mouse was injected with 18F-FDG (8.9 MBq), and another was injected with 18F− (9.4 MBq). One and a half hours later, the mice were centered inside the gantry and imaged for 1 h. CT transmission scans of the mice were obtained using a MicroCAT II tomograph with the same protocol as the CT scan for the image-quality phantom. The PET data were reconstructed with FORE + 2D FBP (with a ramp filter cut off at the Nyquist frequency) and MAP reconstructions (β = 0.01), with all available corrections applied, including attenuation and scatter.
RESULTS
Energy Resolution
Energy resolution of the 511-keV photopeak was 14.6% on average for 25,600 LSO crystals, with 26.9% and 8.2%, respectively, as the worst and best energy resolutions. The energy resolution of the Inveon DPET was significantly improved, compared with the previously reported 18.5% (23), 18.3% (24), and 26% (25) energy resolutions for the Focus220, Focus120, and microPET P4 (Siemens Preclinical Solutions, Inc.), respectively.
Spatial Resolution
Figure 1 shows the radial, tangential, and axial components of the FORE- and FBP-reconstructed point source images. At the center of the FOV, the image resolutions in the transverse planes were below 1.8-mm FWHM and remained under 2.0-mm FWHM within the central 4-cm-diameter FOV.
FORE + 2D FBP reconstructed image resolution of Inveon DPET system as function of radial offset from center of the FOV. FWHM (A) and FWTM (B) of radial, tangential, and axial image resolutions and volumetric resolutions for point sources located at axial center (⧫) and 31.75 mm from axial center (⋄) (C). FWTM = full width at tenth maximum.
Sensitivity
Table 2 summarizes the two measurements of the absolute sensitivity, one at a fixed LLD equalling 350 keV and the other at a fixed ULD equalling 625 keV.
Absolute Sensitivity of Different Energy Window Settings
At the center of the axial and transaxial FOV, the absolute peak sensitivity measured using an energy window of 350–625 keV and a timing window of 3.432 ns was 6.72%. When the LLD was lowered to 250 keV, the absolute sensitivity was 9.32%. This measurement did not include a correction for the self-attenuation of the source. With consideration of self-attenuation, the absolute sensitivity is expected to be higher than 10% at the 250- to 625-keV energy window.
The relative sensitivity measured at different timing windows was compared with the sensitivity measured with a 4.680-ns timing window (the widest timing window available on the Inveon). The measurement shows that the sensitivity does not depend significantly on the timing window (within 1.5%).
The average sensitivity for a mouse-sized object (7 cm) and a rat-sized object (12.7 cm) measured at a 350- to 625-keV energy window and a 3.432-ns timing window was 4.0% and 2.8%, respectively. The axial sensitivity profile is shown in Figure 2.
Sensitivity profile over axial positions.
Scatter Fraction and Counting-Rate Performance
With an energy window of 350–625 keV and a timing window of 3.432 ns, the peak NECR is 1,670 kcps (achieved at 130 MBq) for the mouse-sized phantom and 590 kcps (achieved at 110 MBq) for the rat-sized phantom. The scatter fraction at this acquisition setting is 7.8% and 17.2% for the mouse- and rat-sized phantoms, respectively. The NECR as a function of activity is plotted in Figure 3 for the 2 phantoms.
NECR as function of total activity for mouse- and rat-sized phantoms.
The count loss due to dead time was also investigated on the basis of the mouse-sized phantom data (Fig. 4). The count loss is 25% at an activity of 50 MBq and 50% at an activity of 110 MBq. The dead-time correction works well up to an activity of 50 MBq (within 1% accuracy, compared with the expected counts).
Raw counts, decay-corrected counts, and dead time–corrected counts as function of total activity in FOV based on mouse-sized phantom data.
Imaging Studies
NEMA Phantom Study.
Figure 5 shows the images of a transverse plane with 5 rods (Fig. 5A), a coronal plane (Fig. 5B), a transverse plane of the uniform region (Fig. 5C), and a profile across the uniform area (Fig. 5D) of the NEMA image-quality phantom.
Images of NEMA NU-4 image-quality phantom scanned for 20 min with 18F-FDG (5.1 MBq): transverse plane of 5-rod region (A), coronal view (B), transverse plane of uniform region (C), and profile across uniform area (D).
With FORE + 2D FBP reconstruction, the percentage SD in the uniform region was 5.29 with all corrections applied.
The RCs reconstructed with FBP for 5 different rod sizes from 1- to 5-mm diameter are shown in Figure 6. The RC for the smallest 1-mm rod is 0.17, and for the largest 5-mm rod the RC is 0.93.
RCs for 5 rods of different sizes reconstructed with FBP.
The spillover ratios measured in the air- and water-filled chambers of the NEMA phantom after all the corrections had been applied were −0.57% and 1.65%, respectively. Without attenuation and scatter correction, the spillover ratios for air- and water-filled chambers were 13.95% and −0.21%, respectively. The residual activity is a measurement of correction accuracy. With spillover ratios below 2% after corrections for both air- and water-filled chambers, the corrections work reasonably well.
Mouse Study.
Figure 7 shows the coronal and sagittal images of mice with 18F-FDG and 18F− uptake, respectively. The images were reconstructed with FBP and MAP (β = 0.01).
Coronal plane of healthy mouse scan with injection of 18F-FDG (8.9 MBq) reconstructed with FORE + 2D FBP (A) and MAP (B). Sagittal plane of mouse bone scan with injection of 18F− (9.4 MBq) reconstructed with FORE + 2D FBP (C) and MAP (D).
DISCUSSION
The NEMA standards are not meant to produce absolute performance measurements but rather define measurements that can be easily performed, analyzed, and interpreted by the user community. They serve as tools for comparison of imaging instruments under specific operating conditions. These specifications represent a subset of measurements that define the performance of preclinical PET systems for specific imaging tasks.
The energy and spatial resolutions were primarily determined by the scintillator type and individual crystal size of the system. The LSO-based Inveon DPET system has improved energy resolution and spatial resolution similar to the Focus220, which is also based on an LSO scintillator and has the same crystal element size as the Inveon (23). The increase in energy resolution is mainly due to the use of improved light guides, which reduce light loss, in the Inveon. Because of large solid-angle coverage, the sensitivity of the Inveon is significantly higher than that of the previous systems. In addition, because of the increased sensitivity of the Inveon, the peak NECRs for both mouse- and rat-sized phantoms are also significantly improved, compared with previous preclinical systems. The scatter fractions cannot be directly compared with previous measurements for other systems (23,24) because a different size of phantom and a different location of the line source insert were used here, per the new NEMA NU-4 standard. In this study, we also have the addition of the NEMA NU-4 image-quality phantom that puts together several aspects of the system performance in a single measurement that is relevant to a rather common acquisition protocol.
The measured RCs for different rod sizes are possibly underestimated because of scatter. The RC was calculated by comparing the rod value to the uniform filled background intensity. The rods are surrounded by cold acrylic glass; therefore, there are fewer scattered coincidences in the rod region than in the uniform background area. Although we performed scatter correction during image reconstruction, the scatter correction is possibly less effective for large uniform regions. To investigate this, an increased LLD from 350 to 450 keV was used to image the same image-quality phantom. At that energy window, we saw an increase in the RC values for all rods, demonstrating that the RCs are underestimated at least in part because of scatter. The NEMA mouse image-quality phantom is rather new, and there is limited experience in its use with other systems; therefore, other tomographs need to be investigated to fully evaluate the significance of these findings.
The LSO crystals have an intrinsic background emission. With a wide open energy window, the intrinsic radioactivity will produce a uniform background and affect the ability of detecting low activities in the FOV. The total 176Lu equivalent activity contained in the scintillation crystals is calculated to be 166,500 Bq based on the percentage of radioactive 176Lu. Because of the relative compact geometry and large solid angle of the Inveon DPET, compared with other LSO-based scanners such as the Focus220, the background of the Inveon might be more significant. This LSO background can be reduced by increasing the LLD at the expense of sensitivity. The absolute system sensitivity penalty is about 30% when the LLD increases from 250 to 350 keV.
In the design of the Inveon DPET system, a small detector ring diameter and large axial extent were used to increase the solid-angle coverage and therefore to achieve higher system sensitivity. Because the Inveon DPET has a long axial FOV, it can cover the whole mouse, including the tail, in one scan. The long axial FOV provides the possibility to estimate the total injected activity from the reconstructed image itself. However, because of the relatively small detector ring diameter and long axial extension, compared with other small-animal PET systems, the Inveon is more prone to radial and axial resolution degradation due to the photon penetration at large radial offsets and large ring differences if reconstructed with analytic FBP reconstruction. Statistical reconstruction methods, such as MAP (9,26), should be the preferred image-reconstruction algorithms to obtain high-resolution images. Photon penetration and mispositioning of the coincidence events can be modeled with Monte Carlo simulation and incorporated into a system-response matrix to recover the radial and axial resolution losses (27).
CONCLUSION
This study evaluated the performance of the small-animal PET Inveon system based on the new NEMA NU-4 protocol. Although direct comparison with previous results acquired with nonstandard sources might not be straightforward, the transverse spatial resolution of the Inveon was similar to that of previous generations of small-animal PET scanners. Absolute sensitivity of the system was 9.3% for an energy window of 250–625 keV and a timing window of 3.432 ns, which is significantly improved with respect to previous tomographs. The peak NECR of the Inveon system was also greatly increased. The image-quality phantom test demonstrated that the image uniformity and RCs were good, and the scatter correction worked reasonably well. Further improvements in the overall system performance are expected to be realized with fully 3D iterative reconstruction algorithms that incorporate the estimated system response.
Acknowledgments
We thank the NEMA NU-4 Standards Committee and Stephan Siegel, Charles Landen, and Anne Smith from Siemens Preclinical Solutions, Inc., for valuable discussions and technical support. We also thank Sanghee Cho from the Signal and Image Processing Laboratory of the University of Southern California for providing the MAP-reconstructed images, Judy Edwards and Waldemar Ladno at the small-animal imaging facility of the Crump Institute for Molecular Irmaging, and the cyclotron team at University of California, Los Angeles. This work was supported by a grant from the National Institutes of Health (SAIRP NIH-NCI 2U24 CA092865).
Footnotes
-
Guest Editor: David Townsend, University of Tennessee Medical Center
-
COPYRIGHT © 2009 by the Society of Nuclear Medicine, Inc.
References
- Received for publication July 26, 2008.
- Accepted for publication November 25, 2008.