Abstract
The aim of this work was to evaluate the performance characteristics of a newly developed dedicated breast PET scanner, according to National Electrical Manufacturers Association (NEMA) NU 4-2008 standards. Methods: The dedicated breast PET scanner consists of 4 layers of a 32 × 32 lutetium oxyorthosilicate–based crystal array, a light guide, and a 64-channel position-sensitive photomultiplier tube. The size of a crystal element is 1.44 × 1.44 × 4.5 mm. The detector ring has a large solid angle with a 185-mm aperture and an axial coverage of 155.5 mm. The energy windows at depth of interaction for the first and second layers are 400–800 keV, and those at the third and fourth layers are 100–800 keV. A fixed timing window of 4.5 ns was used for all acquisitions. Spatial resolution, sensitivity, counting rate capabilities, and image quality were evaluated in accordance with NEMA NU 4-2008 standards. Human imaging was performed in addition to the evaluation. Results: Radial, tangential, and axial spatial resolution measured as minimal full width at half maximum approached 1.6, 1.7, and 2.0 mm, respectively, for filtered backprojection reconstruction and 0.8, 0.8, and 0.8 mm, respectively, for dynamic row-action maximum-likelihood algorithm reconstruction. The peak absolute sensitivity of the system was 11.2%. Scatter fraction at the same acquisition settings was 30.1% for the rat-sized phantom. Peak noise-equivalent counting rate and peak true rate for the ratlike phantom was 374 kcps at 25 MBq and 603 kcps at 31 MBq, respectively. In the image-quality phantom study, recovery coefficients and uniformity were 0.04–0.82 and 1.9%, respectively, for standard reconstruction mode and 0.09–0.97 and 4.5%, respectively, for enhanced-resolution mode. Human imaging provided high-contrast images with restricted background noise for standard reconstruction mode and high-resolution images for enhanced-resolution mode. Conclusion: The dedicated breast PET scanner has excellent spatial resolution and high sensitivity. The performance of the dedicated breast PET scanner is considered to be reasonable enough to support its use in breast cancer imaging.
Breast cancer is one of the leading causes of cancer death among women worldwide, accounting for 425,000 deaths in 187 countries in 2010 (1). A strategy for reducing breast cancer mortality is early detection for successful treatment. PET systems dedicated to the breast have been under development since 1994, when Thompson et al. reported the feasibility of the first of these systems (2). The dedicated PET systems are intended to have higher photon sensitivity and improved spatial resolution by bringing the detectors close to the breast and by using smaller detector elements than those for whole-body tomography. Most of the dedicated PET systems perform positron emission mammography with 2 planar or curved detector heads that can compress the breast mildly, providing planar and limited-angle tomography images (3–8). Fully tomographic dedicated PET systems have also been introduced. This type of scanner acquires fully tomographic images by rotating 2 or more planar heads (9–11) or by completely encircling the breast with detectors (12–14). Previous simulation studies have shown that, with depth-of-interaction (DOI) measurement capability, a rectangular tomographic scanner designed to encircle the breast outperforms dual planar designs (with or without DOI) in terms of system sensitivity and image quality (15,16).
The dedicated breast PET scanner developed by Shimadzu Co. is a fully tomographic system designed to provide high-resolution and high-sensitivity images because of its 4-layer DOI measurement capability and large field of view (FOV) in the axial direction. In this work, the performance of this dedicated breast PET camera was evaluated to determine its intrinsic features among dedicated PET systems. The performance assessment was based on the NU 4-2008 small-animal PET standards of the National Electrical Manufacturers Association (NEMA) (17), as camera-testing standards for dedicated breast PET systems have yet to be developed (3). The performance parameters evaluated in this study included spatial resolution, scatter fraction, count losses, random coincidences, sensitivity, and image-quality characteristics.
MATERIALS AND METHODS
System Description
The system consists of 36 detector blocks arranged in 3 contiguous rings, with a crystal ring diameter of 185 mm and an axial extent of 155.5 mm (Table 1; Fig. 1). Each detector block is composed of a 32 × 32 array of lutetium oxyorthosilicate (LGSO) crystals coupled to a 64-channel position-sensitive photomultiplier tube via a light guide. Each crystal is 4.5 mm long and has a cross-sectional area of 1.44 × 1.44 mm (18). The DOI information of 4 layers is extracted from one 2-dimensional position histogram made by Anger-type calculation, which is achieved by insertion of the reflector between crystal elements to control the behavior of scintillation photons (19). A previous study evaluating a prototype DOI detector showed that the spatial resolution of each DOI layer obtained by γ-ray beam scanning in depth along the side face of the detector is 4.44–4.60 mm in full width at half maximum (FWHM) on average (20). List-mode data are acquired during measurements. The energy window was set to 400–800 keV for the first and second crystal layers (the first layer is the patient side) and 100–800 keV for the third and fourth crystal layers. The coincidence timing window was 4.5 ns. A fixed energy window and a timing window were used in this study. The distance between the upper edge of the scanner and the bed surface was 13.65 mm, resulting in an approximately 18-mm-thick invisible area together with a compressed mat (about 3-mm thickness) and a substantially limited-sensitivity region at the edge of the axial FOV (2 image slices = 1.56 mm). Breast outside FOV could be larger, depending on the positioning of the body and the shape of the individual chest wall.
Main Specifications and Characteristics of Dedicated Breast PET Scanner
Dedicated breast PET scanner designed for scanning patients in prone position. Right lower image shows cylinder-shaped scanner.
Spatial Resolution
Spatial resolution was measured with a 22Na point source conforming to NEMA NU 4 standards. The 22Na point source, of a nominal size (0.3 mm), is embedded in an acrylic cube (10.0-mm extent on all sides).
Measurements were acquired with the source at the axial center of the FOV and at one fourth of the axial FOV from the axial center at the following radial distances: 5, 10, 15, 25, 50, and 75 mm. The list-mode data acquired on each measurement were graphed into 3-dimensional (3D) sinograms with delayed-events subtraction for randoms correction, followed by normalized correction. The 3D sinograms were rebinned into 2-dimensional sinograms by Fourier rebinning using John’s equation (21) and then reconstructed by 2-dimensional filtered backprojection with a ramp filter cutoff at the Nyquist frequency. The pixel size in the transverse plane was 0.2 × 0.2 mm, and the axial plane separation was 0.78 mm. The list mode was also reconstructed with a dynamic row-action maximum-likelihood algorithm (DRAMA) (22) using 128 subsets, 1 iteration, relaxation control parameter β = 100, and a voxel size of 0.39 × 0.39 × 0.39 mm. DRAMA is an iterative algorithm similar to the row-action maximum-likelihood algorithm, but the relaxation parameter is controlled in such a way that the noise propagation from projection data to the reconstructed image is independent of the access order of the input data (subsets) in each cycle of subiterations. On the reconstructed point-source images, 1-dimensional response functions were drawn in the radial, tangential, and axial directions, and their FWHM was measured.
Scatter Fraction, Count Losses, and Random Coincidence Measurements
On the basis of a precedent study (3), we used a ratlike phantom made of a solid, right and circular cylinder composed of high-density polyethylene (density, 0.96 g/cm3) 150 mm long and 50 mm in diameter. A cylindric hole (3.2-mm diameter) was drilled parallel to the central axis at a radial distance of 17.5 mm. The line-source insert was made 5 mm shorter on both ends and filled with a well-mixed 18F-FDG water solution that had a starting activity of 40.1 MBq as measured in a dose calibrator. The phantom was placed at the center of the FOV. Data were acquired until true-event losses were less than 1.0% or random-event rates were less than 1% of true rates.
Sensitivity
The same point source and measurements used for the spatial resolution study were used here. A background counting rate was acquired for 15 h and subtracted from the point-source counting rate. The number of true coincidences was normalized to the scan duration, divided by the source activity, and corrected for the branching ratio of 0.906.
The list-mode data were graphed with single-slice rebinning. 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 sensitivity (counts/s/Bq) for each slice i (Si) was calculated as
where Acal is the point-source activity, Ri is the counting rate for slice i, and RB,i is the background counting rate for slice i.
The total sensitivity of the scanner was calculated by summing Si for all slices. Because the branching ratio of 22Na is 0.906, the absolute sensitivity SA,i was specified as a percentage value,
Image Quality, Accuracy of Attenuation, and Scatter Correction
The NEMA NU 4 image-quality phantom made of polymethacrylate has internal dimensions of 50 mm in length and 30 mm in diameter. It consists of 3 parts: the first is 5 fillable rods used to measure noise and recovery coefficients as a function of rod diameter; the second is a large, uniform region connected to the rods, allowing uniformity to be measured; and the third is 2 cold-region chambers (filled with water and air) that are used to quantify spillover ratio. The recovery coefficients, percentage SD, and spillover ratio were obtained as described in the NEMA document. The image-quality phantom was filled with 18F-FDG solution (3.7 MBq), and images were acquired for 20 min. For this study, image reconstruction was performed with DRAMA using 2 different modes: standard reconstruction mode with postfiltering of 1.5 mm FWHM (β = 5) and enhanced-resolution reconstruction mode without postfiltering (β = 100), using 128 subsets, 1 iteration, and a voxel size of 0.78 × 0.78 × 0.78 mm. Attenuation corrections using a uniform attenuation map with object boundaries obtained from emission data and scatter corrections were applied for all images. The scatter correction method used was the convolution subtraction method (23) with kernels obtained by background tail-fitting.
Human Imaging
A patient with histologically proven breast cancer underwent both whole-body PET/CT and dedicated breast PET scanning. First, whole-body PET/CT scanning was performed with 3 min per bed position approximately 60 min after an injection of 188.7 MBq of 18F-FDG, followed by additional scanning for the breast in the prone position approximately 80 min after the injection using a whole-body PET/CT scanner (Discovery ST Elite; GE Healthcare). Then, the patient was brought to the dedicated breast PET imaging room and lay prone on the bed while placing one breast within the scanner. The scan of the entire unilateral breast started approximately 100 min after the injection, with a scanning time of 5 min per breast.
Whole-body PET images were attenuation-corrected using CT data and were reconstructed with a 3D ordered-subsets expectation-maximization algorithm called VUE Point Plus (14 subsets, 2 iterations, a matrix size of 128 × 128, a voxel size of 4.69 × 4.69 × 3.27 mm, and postfiltering at 5.14 mm FWHM). Dedicated breast PET images were reconstructed with the same 2 modes used for the evaluation of the NEMA NU 4 image-quality phantom, in addition to enhanced-resolution mode without DOI information.
Our institutional review board approved this study, and written informed consent was obtained from the subject before participation in this study.
RESULTS
Detector Characterization
For filtered backprojection reconstruction, the radial FWHM at 2 axial positions varied from 1.6 to 2.7 mm and tended to be greater at larger radial offsets (Fig. 2A). Tangential and axial resolution was rather constant over the whole transaxial FOV, with the tangential FWHM of 1.7–2.2 mm and the axial FWHM of 2.0–2.2 mm for all radial offsets. With DRAMA reconstruction, radial, tangential, and axial FWHM was measured as 0.8–1.3, 0.8–1.0, and 0.8–1.0 mm, respectively (Fig. 2B).
Radial, tangential, and axial spatial resolution (FWHM) as function of radial offset reported for filtered backprojection reconstruction (A) and for DRAMA reconstruction (B). Data were measured on transaxial planes at 2 axial positions: center of axial FOV and one fourth from center of axial FOV.
For the ratlike phantom, the peak noise-equivalent counting rate (NECR) at this acquisition setting was 373.8 kcps (achieved at 24.9 MBq) (Fig. 3). The peak true rate was 602.9 kcps (achieved at 30.6 MBq). The scatter fraction was 30.1% for the ratlike phantom.
Counting rate performance plots as function of total activity for ratlike phantom.
The system peak absolute sensitivity was 11.2% (101.4 cps/kBq) at the center of the axial and transaxial FOVs. The axial sensitivity profile obtained by plotting the absolute sensitivity for each slice number (SA,i) is shown in Figure 4.
Axial absolute sensitivity profile along z-axis of dedicated breast PET camera using maximum ring difference of 99.
NEMA Image-Quality Phantom Study
With 3D list-mode DRAMA reconstruction, the percentage SD in the uniform region was 1.9 for standard mode and 4.5 for enhanced-resolution mode. The smallest rod, 1-mm, was not clear with the standard mode (Fig. 5A) but was identifiable with the enhanced-resolution mode (Fig. 5B).
Images of NEMA NU 4 image-quality phantom scanned for 20 min with 18F-FDG (3.7 MBq) and reconstructed with standard mode (A) and enhanced-resolution mode (B). Transverse planes corresponding to the 5 rods are shown. With enhanced-resolution mode, faint uptake for 1-mm rod is identifiable (B; arrow).
The recovery coefficients for rods 1–5 mm in diameter are shown in Figure 6. The recovery coefficients were 0.04 and 0.09 for the 1-mm rod and 0.82 and 0.79 for the 5-mm rod in the standard and enhanced-resolution modes, respectively. The enhanced-resolution mode provided the highest recovery coefficient level, 0.97 for a 4-mm rod.
Recovery coefficients for 5 rods of different sizes reconstructed with 3D list-mode DRAMA. Data are shown for standard mode and enhanced-resolution mode.
The spillover ratios measured in the air- and water-filled chambers of the NEMA phantom were 18% and 21%, respectively, for standard mode, and 7% and 11%, respectively, for enhanced-resolution mode.
Human Imaging
Figure 7 shows images of an invasive breast cancer obtained with the whole-body PET/CT scanner while the patient was prone and images obtained with the dedicated breast PET scanner. Although an exact comparison between the 2 scans is impossible because of a substantial difference in the uptake period, the 5-min dedicated breast PET scan clearly visualized a small tumor 7 mm in diameter with intense uptake in both reconstruction modes, whereas the 3-min whole-body scan showed only mild uptake (Fig. 7A). The standard reconstruction mode provided high-contrast images with less background noise, demonstrating the practicability of the dedicated breast PET camera for cancer detection (Figs. 7B and 7C). In enhanced-resolution mode, although the background noise increased, the intratumoral distribution of 18F-FDG was clearly demonstrated, suggesting the potential utility of this camera for evaluating the detailed intratumoral morphologic and metabolic characteristics of the target lesion (Fig. 7D). When DOI measurement was not adopted, enhanced-resolution images became blurred because of decreased spatial resolution (Fig. 7E).
A 74-y-old woman with invasive breast cancer 7 mm in diameter. (A) On axial image of whole-body PET/CT study in prone position, tumor shows faint uptake providing maximum standardized uptake value of 1.0 (arrow), almost same activity level as right atrial blood pool (asterisk). (B and C) In dedicated breast PET study, maximum-intensity-projection image (B) and transaxial image (C) reconstructed with standard mode with DOI clearly depict tumor as intense focal uptake (arrow). (D) Transaxial image reconstructed with enhanced resolution with DOI shows detailed 18F-FDG distribution within tumor (arrow). (E) Contrastingly, transaxial image reconstructed with enhanced resolution without DOI is blurred compared with that reconstructed with DOI, resulting in obscuring of intratumoral 18F-FDG distribution. (F) Low-power photomicrograph (hematoxylin-eosin stain) shows that tumor cells are arranged in nests separated by fibrous septae.
DISCUSSION
The radial spatial resolution near the center of transaxial planes obtained at 2 axial positions approached 1.6 mm under the filtered backprojection reconstruction, a value that is superior to that reported for a commercially available breast positron emission mammography camera (3) and as good as a ring-type dedicated PET system for small animals (24). Spatial resolution at greater radial offsets from the center was degraded, as observed for other dedicated breast PET systems or small-animal PET systems (3,24–26), but remained at less than 2.7 mm within the central 15-cm-diameter FOV. Tangential and axial resolutions ranged from 1.7 to 2.3 mm and from 2.0 to 2.2 mm, respectively. When reconstructed with DRAMA, dedicated breast PET images achieved high and relatively invariant spatial resolution, with the radial FWHM being 0.8–1.3 mm, the tangential FWHM being 0.8–1.0 mm, and the axial FWHM being 0.8–1.0 mm. In a similar way to that described by Luo et al. (3), we also evaluated spatial resolutions with a uniform background (Supplemental Fig. 1, available at http://jnm.snmjournals.org) and found that differences in FWHM with and without the background activity were within 0.01 mm (Supplemental Table 1). This high resolution was realized by reducing the LGSO crystal size to 1.44 × 1.44 × 4.5 mm, and the high light output of the LGSO enables discrimination of input signals from different crystals. The DOI information was effective at decreasing degradation in radial resolution due to parallax error (27). The effectiveness of DOI was also demonstrated in the human imaging.
The scatter fraction was 30%. This was somewhat higher than the scatter fractions of currently available preclinical PET scanners (3,28,29), possibly because of the large solid angle of this dedicated breast PET camera and the low-energy threshold for the third and the fourth crystal layers. The energy window of the third and fourth layers was set to 100–800 keV for maintenance of sensitivity. Because the scintillation decay of LGSO is 40 ns, far shorter than that of bismuth germanium oxide, the coincidence time window was set to 4.5 ns to reduce the random rate and to maintain a high NECR. The measured peak NECR (373.8 kcps) was comparable to those of currently available preclinical PET scanners. Conversely, the measured random coincidence was also a low counting rate. In the typical clinical situation in which the normal breast-tissue activity concentration is assumed to be 2 kBq/mL, which was calculated as 5.3 kBq/mL (soft-tissue activity) (30) × 0.33 (average SUV of normal bilateral breast) (31) without consideration of activity outside the FOV, the NECR is approximately 20 kcps, a value that is sufficiently higher than that for other dedicated breast cameras (3,14). However, the optimum time window and energy window for clinical positron emission mammography may need further investigation, and peak NECR may be increased.
The measured sensitivities (11.2%; 101.4 cps/kBq) were comparable to or more than twice those of currently available breast PET cameras (3,6,10). This was achieved by the large solid angle of the detector crystals, with a high packing fraction, and the fourth-layer DOI detector (18 mm = 4.5 × 4). However, the axial sensitivity profile showed a linear drop of sensitivity from the center to the edge of the axial FOV (Fig. 4).
The results obtained with the image-quality phantom revealed good uniformity in uniform regions and high contrast in hot regions. Among several cameras that have been evaluated with NEMA NU 4 standards (3,24,28), the present dedicated breast camera provided comparable or higher visual contrast with acceptable uniformity in enhanced-resolution mode, and comparable or slightly lower contrast with excellent uniformity in standard mode. Such a property may increase the detectability of tumors, contribute to the evaluation of intratumoral morphologic and metabolic characteristics, and improve the quantification of 18F-FDG uptake. With enhanced-resolution mode, the recovery coefficient for the 4-mm rod was higher than that for the largest 5-mm rod. This effect could possibly be attributed to an overshoot at the edges, commonly referred to as the Gibbs effect (29,32). The measured activities in cold regions, that is, the water- and air-filled inserts, for enhanced-resolution mode without postfiltering processing were comparable to data obtained with another scanner (28), suggesting that corrections for attenuation and scatter were done appropriately. In general, preclinical PET attenuation correction is done with CT or with a 57Co point source by singles detection. The present dedicated breast PET camera used a uniform attenuation map. Because most of the breast consists of adipose and fibroglandular tissue, it is considered acceptable to presume that the breast has a uniform attenuation coefficient. Furthermore, the dedicated breast PET camera is primarily designed as a breast cancer screening test and could be used after whole-body PET/CT screening. The present dedicated breast PET camera has the great advantage of achieving safe screening without additional radiation exposure for the purpose of attenuation correction. However, a possible drawback of this method is the uncertain applicability to PET agents other than 18F-FDG, especially specific tracers, since it remains unknown whether these agents accumulate in the normal tissue enough to delineate the breast and enable an attenuation map to be composed. Further studies are needed to clarify the utility of this attenuation correction method for various PET agents.
In the human imaging study, the dedicated breast PET camera clearly better depicted a small tumor than did whole-body PET/CT. The enhanced-resolution mode demonstrated remarkably high-resolution images capable of depicting detailed 18F-FDG distribution within the small tumor. However, the noise level determined by visual inspection was somewhat high with the enhanced-resolution mode but significantly low with the standard reconstruction mode. There was a trade-off between spatial resolution and noise. The dynamic relaxation parameter of 3D list-mode DRAMA is adjusted by β, and it is possible to control propagation of noise from projection data to the final image. A cancer screening test requires low false-positive rates in addition to high sensitivity. Detailed morphologic and metabolic information may help further characterize the target lesion. It seems important to optimize the parameters of DRAMA based on the purpose of the study.
CONCLUSION
According to the present assessment based on the NEMA NU 4-2008 standards, the dedicated breast PET scanner has excellent spatial resolution and high sensitivity. Although NEMA NU 4 may not fully serve as an evaluation standard for dedicated breast PET systems—for example, because of the lack of consideration about activity outside the FOV—the performance of the present dedicated breast PET scanner is considered to be reasonable enough to support its use in breast cancer imaging.
DISCLOSURE
The costs of publication of this article were defrayed in part by the payment of page charges. Therefore, and solely to indicate this fact, this article is hereby marked “advertisement” in accordance with 18 USC section 1734. This work was supported by Shimadzu Co., Kyoto, Japan, and by a Grant-in-Aid for Scientific Research from the Ministry of Education, Culture, Sports, Science and Technology of Japan (C:22591329). Authors of this paper include employees of the company. Other authors have control of any data that might present a conflict of interest for employees. No other potential conflict of interest relevant to this article was reported.
Footnotes
Published online May 8, 2014.
- © 2014 by the Society of Nuclear Medicine and Molecular Imaging, Inc.
REFERENCES
- Received for publication September 3, 2013.
- Accepted for publication March 17, 2014.