|
|
|||||||||
Clinical Investigations |
1 Department of Nuclear Medicine, University Hospital, Nancy, France
2 UHP-INSERM U684, Nancy, France
3 UHP-INSERM ERI13, Nancy, France
4 Department of Nuclear Medicine, Arizona-VUB, Brussels, Belgium
| ABSTRACT |
|---|
|
|
|---|
Key Words: gated SPECT MRI ejection fraction sestamibi 201Tl
| INTRODUCTION |
|---|
|
|
|---|
This study aimed to determine whether the association of OSEM with depth-dependant resolution recovery prefiltering and temporal Fourier filtering, which is likely to improve both signal-to-noise ratio and spatial resolution, further enhances the results provided by 201Tl and sestamibi 16-interval gated SPECT. Cardiac MRI was the reference method, and a comparison was planned with a conventional FBP method previously optimized by screening various filter frequencies and various temporal smoothing levels.
| MATERIALS AND METHODS |
|---|
|
|
|---|
201Tl and Sestamibi Gated SPECT
Gated SPECT Acquisitions.
As already described (79), 1.5 MBq of 201Tl per kilogram of body weight were injected during exercise, and 0.5 MBq was injected 10 min before rest-reinjection imaging. Acquisitions were obtained on a double-head camera (DST-XL; SMV-GE Healthcare) equipped with high-resolution parallel-hole collimators, 32 projections being recorded on a 180° circular orbit from the 45° left posterior oblique to the 45° right anterior oblique orientations. Additional parameters were as follows: 16 bins, 50 s per projection, 50% acceptance window on the averaged length of cardiac cycles, 1.33 zoom, 43 x 43 cm field of view, 64 x 64 matrix, and a dual-energy window (6090 and 150184 keV).
When 201Tl acquisitions had been completed, 11 MBq per kilogram of body weight of 99mTc-sestamibi were injected at rest, and an additional gated SPECT acquisition was started 5060 min later using the same parameters except for energy window (126154 keV) and time per projection (30 s).
Gated SPECT Reconstructions.
For each of the 2 tracers, 25 different methods of reconstruction were analyzed (Fig. 1), including FBP and Butterworth prefilters with an order of 5 and various cutoff frequencies (0.20, 0.25, 0.30, 0.35, and 0.40 cycles per pixel) (10), and a depth-dependant resolution recovery prefilter followed by OSEM reconstruction with 2 iterations and 8 subsets (6). All these reconstructions were also obtained after an additional 2-, 3-, or 4-frame temporal smoothing of raw projections, each of the 16 bins from each projection being averaged with, respectively, the 1, 2, or 3 following bins. Finally, OSEM and depth-dependant resolution recovery prefiltering were also combined with temporal Fourier filtering applied after each iteration (5). Reconstruction time was around 1 min longer for depth-dependant resolution recovery with temporal Fourier filtering, when compared with FBP.
|
Temporal Fourier filtering was applied to a Fourier series transform of the timeactivity curves of the voxels and involves the suppression of all harmonics higher than the second one (those with a higher noise level) followed by an inverse Fourier transform (5,14).
Gated SPECT Analyses.
LV volume and ejection fraction were determined automatically using Quantitative Gated SPECT software (Cedars-Sinai Medical Center (1,15)), constrained limits of contour detection being applied only in case of evident errors with the fully automatic process. Regional contractility was assessed visually by an experienced observer with a 17-segment LV division (16), the 2 proximal septal segments being excluded. According to wall motion and to the systolic increase in myocardial count, each segment was classified as akinetic/dyskinetic, hypokinetic, or normal.
To estimate signal-to-noise ratio, we selected a midventricular end-diastolic short-axis slice on each reconstruction. Mean myocardial counts (MMy) and corresponding SD (SDMy) were determined on a ring-shaped myocardial region of interest. Mean noise counts (MNo) and corresponding SD (SDNo) were determined on a half-moonshaped region of interest, close to the left side of the lateral wall. Signal-to-noise ratio was estimated with the following formula: (MMy MNo)/
.
Acquisition and Analysis of Cardiac MRI
A breath-holding electrocardiography-triggered gradient-echo sequence (17,18) was used on a 1.5-T magnet (Infinity; GE Healthcare), with a phased-array surface coil and the following parameters: 8-mm slice thickness, 45-ms image frame duration, 30 x 30 cm field of view, and 256 x 128 matrix. Contiguous short-axis planes covering the entire LV volume were recorded, along with midventricular horizontal and horizontal slices. End-diastolic and end-systolic LV volumes, as well as LV ejection fraction, were measured on the contiguous short-axis planes using mass analysis software (Leiden University and MEDIS Medical Imaging Systems). Regional contractility was assessed visually by an experienced observer, according to the systolic increase in myocardial thickness and to wall motion and by using the same LV division and segment classification (1, akinetic or dyskinetic; 2, hypokinetic; 3, normal) as was used for gated SPECT.
Statistical Analysis
Continuous variables were expressed as mean ± SD. In a first step, the LV ejection fractions and end-diastolic volumes provided by each method of gated SPECT reconstruction and by each tracer were compared with corresponding MRI values by using linear regression analyses and paired comparisons of absolute values of relative errors (Wilcoxon tests). Relative error was computed as difference divided by mean of gated SPECT and MRI values.
In a second step, reconstruction methods providing the best determinations of LV ejection fraction and volume were selected for multiple paired comparisons of signal-to-noise ratio (Wilcoxon tests) and of exact agreement rates in segmental contractility analysis (McNemar tests). An exact agreement between gated SPECT and MRI required that the 2 techniques provide a similar contractility classification for a given segment but also that the quality of gated SPECT images be high enough for visual analysis of contractility. Contractility analyses were also assessed by determining the
-values provided by visually analyzable reconstructed gated SPECT. A P value of <0.05 was considered statistically significant.
| RESULTS |
|---|
|
|
|---|
Determinations of LV Ejection Fraction and Volume
On the FBP reconstructions provided by both tracers, the cutoff frequency of the Butterworth prefilter, but not temporal smoothing, had a clear impact on the quantifications of LV ejection fraction (Fig. 2) and volume (Fig. 3). For sestamibi, the best results were achieved at a cutoff of at least 0.30 cycles per pixel, whatever the level of temporal smoothing. Lower frequencies were associated with significant increases in the relative error for volume determination (Fig. 3). For 201Tl, the best results were achieved at a cutoff of around 0.30 cycles per pixel. Lower or higher frequencies were associated with a decrease in correlation coefficient or with an increase in relative error for ejection fraction (Fig. 2), end-diastolic volume (Fig. 3), or both. Finally, whatever the considered tracer, 0.30 cycles per pixel was the lower cutoff frequency and, thus, the cutoff frequency leading to the higher image quality among those providing the best FBP determinations of LV ejection fraction and volume.
|
|
LV volumes were lower than those from MRI, whatever the reconstruction method (all P < 0.001), and these volumes were not better determined with depth-dependant resolution recovery than with an optimal FBP (Butterworth filter of 0.30 cycles per pixel) (Fig. 3).
Compared with 201Tl, sestamibi provided higher correlation coefficients and lower relative errors in LV volume determination whatever the method of reconstruction. For LV ejection fraction, results from the 2 tracers were more comparable, even though sestamibi provided higher correlation coefficients in most FBP reconstructions.
Signal-to-Noise Ratio and Regional Contractility Analysis
Factors involved in the visual analysis of gated SPECT images were compared between a reference FBP reconstruction (filter frequency, 0.30 cycle per pixel, no temporal smoothing) and the 3 reconstruction methods providing the best determinations of LV ejection fraction: the association of OSEM with depth-dependant resolution recovery prefiltering when used alone or with the addition of either 2-frame temporal smoothing or temporal Fourier filtering.
As detailed in Figure 4, whatever the considered tracer, gradual increases in signal-to-noise ratio were documented between FBP and the association of OSEM with depth-dependant resolution recovery prefiltering and temporal Fourier filtering. In addition, the increases in signal-to-noise ratios that were observed between these 2 reconstruction methods were clearly higher for 201Tl (+80% ± 46%) than for sestamibi (+54% ± 40%, P < 0.001). Finally, whatever the method of SPECT reconstruction, signal-to-noise ratio was always significantly higher with sestamibi than with 201Tl.
|
Results of segmental contractility analysis are detailed in Figure 5. The association of OSEM with depth-dependant resolution recovery prefiltering and temporal Fourier filtering provided the higher
-values (sestamibi, 0.68; 201Tl, 0.59), and the higher exact agreement rates (sestamibi, 79%; 201Tl, 73%), and these rates were significantly better than those provided by the 3 other reconstruction methods. In addition, whatever the method of SPECT reconstruction,
-values and exact agreement rates were higher with sestamibi than with 201Tl. As illustrated by an example in Figure 6, the association of OSEM with depth-dependant resolution recovery prefiltering and temporal Fourier filtering consistently enhanced image quality and facilitated visual interpretation of gated SPECT images.
|
|
| DISCUSSION |
|---|
|
|
|---|
Ejection fractions were determined with correlation coefficients higher than 0.9, without any underestimation, and this finding agrees with previous observations in 16- but not 8-interval gated SPECT studies (1,2). In addition, this determination was clearly better than that from FBP, with lower relative errors and higher correlation coefficients with MRI, even when temporal smoothing was used to lower the noise of FBP images and despite the fact that optimal frequencies of FBP had previously been selected.
The association of OSEM with depth-dependant resolution recovery prefiltering and temporal Fourier filtering also enhanced the visual analysis of segmental contractility: All studies became of sufficient quality to be analyzed visually, and higher
-values were obtained from visually analyzable reconstructions (Fig. 5). Exact agreement rates with MRI were high for sestamibi (79%) and a little lower for 201Tl (73%). Such high rates had already been reported with tetrofosmin or sestamibi (20,21) but not with 201Tl (21,22).
By correcting collimator response, the depth-dependant resolution recovery prefilter improves and homogenizes spatial resolution on the overall myocardial volume (6,1113), and the OSEM iterative method presents numerous advantages over FBP, including better noise properties (3,4) and the possible incorporation of temporal Fourier filtering at each iteration (5). Temporal Fourier filtering consists of a truncating Fourier series of the voxel count variations between interval frames. In this way, SPECT images are obtained at each interval using the activity recorded at all intervals. This process is, therefore, definitely different from that of usual reconstructions, which are performed at each interval on only a small fraction of the recorded activity, especially when 16-interval gating is used.
On the low-count images provided by 201Tl gated SPECT, the association of OSEM with a depth-dependant resolution recovery prefilter has already been shown to substantially improve LV volume determination (6). In addition, applying temporal Fourier filtering after OSEM iterations has been found to be associated with introduction of a negligible error but also with a dramatic reduction in noise, corresponding to a more than 2-fold increase in time acquisition (5). In agreement with this observation, our study shows that, when associated with OSEM and depth-dependent recovery prefiltering, temporal Fourier filtering improves the signal-to-noise ratio and hence, visual analysis of regional contractility, without any deterioration in the determination of LV volume and ejection fraction. By contrast, determination of LV ejection fraction was affected when strong 3- or 4-frame temporal smoothing was added to depth-dependant resolution recovery reconstructions.
A limitation of the study is that the specific impacts of OSEM, temporal Fourier filtering, and depth-dependant resolution recovery filtering were not analyzed separately. These impacts will require further investigations. However, the intent of this study was to investigate only the likelihood that associationssuch as that of OSEM with temporal Fourier filtering and depth-dependant resolution recovery filteringwould reduce the noise of 16-interval gated SPECT images without affecting image resolution. A decrease in image resolution is, indeed, a common disadvantage of usual noise reduction methods (decrease in filter frequency, increase in voxel size, temporal smoothing).
The improvement that was observed between FBP and the association of OSEM with temporal Fourier filtering and depth-dependant resolution recovery filtering was more marked for low image counts and especially for 201Tl image counts. Therefore, this improvement could be less for higher image counts such as those recorded with higher injected activities and longer scanning times. Such longer times, however, increase patients discomfort and the risk of motion during acquisition, a crucial factor in the reliability of gated SPECT results (7). Moreover, further improvements in the signal-to-noise ratio of high-count images might offer the opportunity to use lower voxel sizes at acquisition (zooming, increase in matrix size). According to a previous study (23), this option might prevent the LV volume underestimation that has been documented consistently in previous gated SPECT studies (5,20,24) and in the present study.
| CONCLUSION |
|---|
|
|
|---|
| ACKNOWLEDGMENTS |
|---|
| FOOTNOTES |
|---|
For correspondence or reprints contact: Pierre-Yves Marie, MD, Service de Médecine Nucléaire, CHU-Nancy, Hôpital de Brabois, 54511 Vandoeuvre Cedex, France.
E-mail: py.marie{at}chu-nancy.fr
| REFERENCES |
|---|
|
|
|---|
Related articles in JNM:
This article has been cited by other articles:
![]() |
L. Ceriani, T. Ruberto, A. B. Delaloye, J. O. Prior, and L. Giovanella Three-Dimensional Ordered-Subset Expectation Maximization Iterative Protocol for Evaluation of Left Ventricular Volumes and Function by Quantitative Gated SPECT: A Dynamic Phantom Study J. Nucl. Med. Technol., March 1, 2010; 38(1): 18 - 23. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | RSS | TABLE OF CONTENTS |
| JOURNAL OF NUCLEAR MEDICINE TECHNOLOGY | THE JOURNAL OF NUCLEAR MEDICINE |