Study of Image Contrast Recovery Coefficient (CRC) of a Large Ring PET Scanner ()
1. Introduction
Positron emission tomography (PET) is a nuclear imaging technique which has been widely used over the last few decades in many clinical applications such as staging, evaluating treatment response, and predicting prognosis in malignant diseases [1] [2] .
A large ring PET “macroPET” system was constructed as a prototype by reconfiguring components from an original ECAT 951 system [3] . The original system had 32 detector modules (“buckets” or “packs”) mounted in two rings with an inner diameter of 100 cm and a 10.8 cm axial FOV. But in the macroPET system [4] , the 32 modules are mounted in a single ring with an inner diameter of 2.34 m. As in original system, the modified design has also 128 BGO blocks, 512 PMTs, 8192 individual detectors, 8 crystal rings and 1024 detectors per crystal ring. It was not practicable to construct dedicated septa [5] for macroPET, but in an attempt to investigate the value of septa some measurements were performed using the septa from the original ECAT 951 and the results compared with those obtained without any septa. Septa are lead or tungsten materials which are extended the scanner allow only data for direct and adjacent planes to be collected. Septa reduce most of the scatter events and also reduce photon flux from outside the field of view (FOV), but block many true events, and hence limit the scanner sensitivity.
Apart from these septa, macroPET has no shielding against out of field activity, which can give rise to high backgrounds of random and scattered events. Ultimately if a large ring scanner were to be used in real studies it would be necessary to add side shielding, but for the present trials these problems have been partly avoided by using phantoms. A large ring PET “macroPET” system was constructed as a prototype to demonstrate the feasibility of performing PET scans on a large scale. After design, setup and calibration of the macro PET system, initial results of its performance characteristics were experimentally studied [6] [7] . In this study, an effort has been made to investigate the influence of filter and cut-off frequency on image contrast with septa and without septa acquisition.
2. Methods
A 35.7 cm2 square phantom with warm background embedded by 3 hot and 2 cold lesions (five small cylinders) was scanned in septa and no-septa mode for the contrast recovery (CR) study. The cylinders containing high activity (hot lesions) were 4, 2, and 1 cm diameter and all have the same uptake value in each acquisition mode. The diameters of cylinders containing inactive water (cold lesions) were 3 and 1 cm. The phantom used in the study was filled with F-18 radioisotope to a depth of about 5 cm to avoid any effect of out-of-field of activity. The concentration of F-18 radioactivity ratio of the hot lesions to the warm background was approximately 7.2:1 for no-septa and 8.5:1 for septa mode.
Images were reconstructed using Hann, Hamming and Parsen filters with different cut-off frequencies.
The Hann function, named after the Austrian meteorologist Julius von Hann, is a discrete window function given by
and sometimes referred to as Hanning, presumably due to its linguistic and formulaic similarities to the Hamming window.
In each case, attenuation correction was performed using attenuation factors for each measured LOR which had previously been calculated using the known geometry and an attenuation coefficient of 0.095 cm−1 for 511 keV gamma-rays in water. This attenuation correction is quite sensitive to the orientation assumed for the phantom. Trial and error was used to find the position and orientation of the phantom in the attenuation calculations which gave minimum artefacts in the reconstructed images.
The hot lesion contrast recovery coefficient (CRChot) is calculated as [8] :
(1)
where Chot and Cwarm are the average of the counts measured using the PET system in the hot sphere ROI and the average of the counts in all warm ROIs respectively. ahot and awarm are the activities in the hot and warm background regions respectively. The cold sphere CRC (CRCcold) is calculated as:
(2)
where Ccold is the average of the counts measured in the cold sphere ROI. Counts were taken from the same place by drawing a 1 cm diameter ROI on the reconstructed PET image for each lesion.
3. Results
Figure 1 shows the image slice for the square phantom generated with two acquisition times. The images were reconstructed by a Hamming filter with cut-off frequency of 0.4. In the two counting statistics, three hot and two cold lesions are visible in each case. Image slice with 100 M events are more clearly visible than the slice with 50 M events as the more counts reduce the image noise, and hence improves the image appearance.
The square phantom images reconstructed using Hann, Hamming and Parsen filters and cut-off frequencies of 0.4 and 0.8 are shown in Figure 2. The images reconstructed by acquiring approximately 100 M coincidence events (13.3%
(a) (b)
Figure 1. Central slice for the 35.7 cm square image quality phantom generated two acquisition times in no-septa mode. (a) 50 M counts (20.74% randoms); (b) 100 M counts (13.3% randoms).
No-septa results
randoms) in no-septa configuration. Three hot and two cold lesions are observable in all filter cases. There is no significant difference between the images generated from three filters in case of lesion detection as expected. In our previous study [9] , it was very hard to make conclusion to choose the better filter from the reconstructed images. But the image generated by the Hamming filter with cut-off frequency of 0.4 and 0.8 may be slightly better visible than the other images.
Figure 3 shows the with-septa images generated using three filters (Hann, Hamming and Parsen) and its cut-off frequencies of 0.4 and 0.8. The images reconstruction performed by acquiring approximately 169 M events (12% randoms)
Septa-mode results
with septa mode. Lesions (3 hot and 2 cold) are visible in case of Hamming and Parsen filter. The 1 cm cold lesion is very difficult to detect but 1 cm hot lesion is very clearly visible in all cases. No significant changes are observable among the images generated by different filters.
The contrast recovery coefficient (CRC) as a function of the cut-off frequency for the 4 cm hot lesion with no-septa is shown in Figure 4. Figure 5 shows the corresponding results for the configuration with septa. The CRC results for all filter cases (Hamming, Hann and Parsen) are very similar, and the values are close to unity (slightly higher for the with septa mode) above a cut-off frequency of 0.4. The CRC curve is almost flat from cut-off frequency of 0.4 to 1.0 for both modes. For the large lesion, there is a slight loss (around 10%) in CRC at cut-off frequency 0.2.
Figure 6 shows the CRC as a function of the cut-off frequency for the 2 cm hot lesion with no-septa and Figure 7 for the septa configuration. The CRC values are similar for both Hamming and Hann filters but slightly lower using the Parsen filter for both modes. Without septa, the values are approximately 0.85 above a cut-off frequency of 0.4. Adding septa improves these values by approximately 5% - 12%. There is a significant reduction in contrast recovery coefficient at a cut-off frequency of 0.2.
The CRC as a function of the cut-off frequency for the 1 cm hot lesion with no-septa is shown in Figure 8 and for the with septa configuration is shown in Figure 9. Again, both the Hamming and Hann filters give very similar results but the Parsen filter gives a lower CRC in both modes. The CRC for this small lesion is around 0.3, and addition of septa does not appreciably improve the contrast at cut-off frequencies above 0.6.
Figure 4. Contrast recovery coefficient for 4 cm hot lesion as a function of cut-off frequency for three filters with no-septa mode (slice 4).
Figure 5. Contrast recovery coefficient for 4 cm hot lesion as a function of cut-off frequency for three filters with septa mode (slice 4).
Figure 6. Contrast recovery coefficient for 2 cm hot lesion as a function of cut-off frequency for three filters with no-septa mode (slice 4).
Figure 7. Contrast recovery coefficient for 2 cm hot lesion as a function of cut-off frequency for three filters with septa mode (slice 4).
Figure 8. Contrast recovery coefficient for 1 cm hot lesion as a function of cut-off frequency for three filters with no-septa mode (slice 4).
Figure 9. Contrast recovery coefficient for 1 cm hot lesion as a function of cut-off frequency for three filters with septa mode (slice 4).
4. Discussions
For the large (4 cm) lesion good values of CRC are obtained for most values of cut-off frequency, and addition of septa improves the CRC by around 10%, showing that the presence of scatter background is affecting these values. For the small (1 cm) lesion, contrast is limited by the spatial resolution of the scanner, and addition of septa has little effect. For the intermediate (2 cm) lesion, the results depend both on cut-off frequency and on the presence of septa. A general comparison of CRC results for this lesion is shown in Figure 10. The average of the CRC results for cut-off frequencies 0.4 to 1.0 is 0.83 without septa and 0.90 with septa, so that the septa increase the contrast by approximately 8% for this lesion.
The CRC for cold lesions as a function of cut-off frequency for Hamming window function in both acquisition modes is shown Figure 11. For the 3 cm cold lesion, good values of CRC are obtained for all cut-off frequencies above 0.4, and addition of septa appears to make little difference. For the 1 cm cold lesion, CRC results of up to 0.4 are obtained without septa, but addition of septa significantly worsens the CRC. The reason for this is not understood.
Figure 12 shows a horizontal profile running through the middle of the three hot lesions of the image plane 4 reconstructed by Hamming and Parsen filter with cut-off frequencies of 0.4 and 0.8. The peaks for three hot lesions (4 cm, 2 cm and 1 cm) are clearly distinguishable in the profile. The peak to warm background ratio for Hamming filter is slightly higher for 2 cm and 1 cm hot lesions than Parsen filter, and the profile result is consistent with Figures 4-9 and also with Figure 13 for profile running along the two cold lesions.
In the contrast phantom, CRC values close to unity can be achieved for the 2 cm and 4 cm hot lesions and for the 3 cm cold lesion (Figure 4). In For both the 1 cm hot and cold lesions the CRC is around 0.3. Contrast recovery is reduced when filter cut-off frequencies less than 0.4 are used. There is surprisingly little
Figure 10. Comparison of CRC for the 2 cm hot lesion as a function of cut-off frequency for Hamming filters in septa and no-septa configurations (slice 4).
Figure 11. CRC for cold lesions as a function of cut-off frequency for Hamming window function in both acquisition modes (using slice 4).
(a) (b)
Figure 12. Horizontal profile running along the middle of three hot lesions (4 cm, 2 cm and 1 cm) of slice 4 reconstructed by Hamming and Parsen filter with cut-off frequencies of 0.4 and 0.8: (a) no-septa mode; (b) septa mode.
(a) (b)
Figure 13. Horizontal profiles running along the middle of two cold lesions (3 cm and 1 cm): (a) in the image slice 4 (no-septa mode); (b) slice 10 (septa mode) reconstructed by Hamming and Parsen filter with cut-off frequency of 0.6.
difference between the results obtained with and without septa. For hot lesions, septa apparently improve CRC slightly, but for cold lesions CRC is slightly poorer using septa.
5. Conclusion
Using cut-off frequencies > 0.4, contrast recovery for both hot and cold lesions is excellent for diameters ≥ 3 cm. For a 2 cm hot lesion CRC is around 0.8 to 0.9. CRC for 1 cm hot and cold lesions is ~0.3, as expected. Using cut-off frequencies < 0.4 reduces contrast. There is surprisingly little difference between results with and without septa. For hot lesions, septa appear to improve CRC slightly, but for cold lesions CRC is slightly poorer using septa. Overall, the scanner performs as expected, taking account of its lower sensitivity and slightly poorer spatial resolution as compared to a clinical scanner. The main limitation in the study was the lack of shielding against out-of-field activity. In most cases measurements were performed with little out-of-field activity. The presence of such activity would add significantly to the background of random and scattered events. The septa used have only a limited effect in blocking out-of-field activity.