|
|
||||||||
CLINICAL INVESTIGATIONS |
Division of Nuclear Medicine, University Hospital, Free University of Brussels, Brussels, Belgium
| ABSTRACT |
|---|
|
|
|---|
Key Words: blood-pool SPECT automatic algorithm left ventricular ejection fraction
| INTRODUCTION |
|---|
|
|
|---|
In contrast to the widespread use of gated SPECT in myocardial perfusion scintigraphy (1), GBPT has not yet become routine for cardiac gated blood-pool studies. The small dimensions of the myocardial wall compared with the volume of left ventricular blood causes partial-volume effects in gated myocardial perfusion tomography; this error is less of a problem in GBPT. Consequently, cardiac volumes can be estimated more accurately with GBPT than with gated myocardial perfusion tomography. Moreover, the number of cardiac counts collected per millicurie of injected radioactivity is higher for GBPT (2), and the separation between left and right ventricles conceptually allows blood-pool studies to measure right ventricular function.
Despite these advantages, GBPT has not yet become widely used because of the lack of automatic, fast, and reliable quantitative algorithms for measuring ventricular ejection fraction from GBPT. Threshold methods, count-based methods, and local gradient methods have been proposed (316). In the threshold method, voxels higher than a predetermined threshold are included as part of the cardiac volume. This threshold is usually determined from phantom studies. In the count-based method, the total volume of the cavity is derived from the number of counts measured within a small volume of interest of known dimensions. In the local gradient method, the edge of the cavity volume is detected using the first and second derivatives along count density profiles. Most of these methods require a degree of manual operation for masking, setting threshold values, and setting limits for a boundary search.
The aim of this study was to develop and validate a new algorithm to calculate left ventricular ejection fraction (LVEF) from GBPT. In this algorithm, an automatic three-dimensional segmentation technique was applied, combining the threshold and the local gradient methods, to delineate the left ventricular cavity. The results were compared with those of PRNA.
| MATERIALS AND METHODS |
|---|
|
|
|---|
PRNA
PRNA was performed using 16 electrocardiographic gated frames, forwardbackward framing, and a beat acceptance window at 20% of the average R-R interval obtained just before starting acquisition of a left anterior oblique projection (best septal view). The projection was acquired using one head of a triple-head gamma camera (MultiSPECT3; Siemens, Hoffman Estates, IL) equipped with a low-energy, high-resolution collimator. Data were acquired in 64 x 64 format with a zoom of 2.29 (pixel size, 3.1 mm) until a total density of 4 million counts was reached.
For determination of LVEF, the images were transferred onto an NXT (P) computer (Sopha Medical Vision, Paris, France) and processed with a commercially available fully automatic program, Gated Equilibrium Bloodpool Processing, version 2.0.
GBPT
The GBPT studies were acquired with a MultiSPECT3 gamma camera equipped with low-energy, high-resolution collimators. Parameters of acquisition were as follows: 360° stop-and-shoot rotation, 32 stops per head (96 projections), 30 s per stop, 64 x 64 matrix, zoom of 1.23 (5.79-mm pixel size), eight time bins using 75% forwardbackward framing, and a beat acceptance window at 20% of the average R-R interval calculated just before starting the acquisition.
The projection data were reconstructed by filtered backprojection applying a Butterworth filter with a cutoff frequency of 0.35 cycles per centimeter and an order of 5. The reconstructed slices were reoriented according to the left ventricular long axis. The gated horizontal long-axis slices were used as input for the new algorithm.
The New Algorithm
The new algorithm uses geometric methods to delineate the left ventricular cavity, is fully automatic, and works in three-dimensional space. The execution time is approximately 45 s on a Macintosh G3 computer (Apple Computer, Inc., Cupertino, CA). It comprises the following steps: detection of the valve plane, detection of the interventricular septum, delineation of the left ventricular cavity at end-diastole, delineation of the left ventricular cavity for the other time bins, and calculation of left ventricular volumes and ejection fraction.
Detection of Valve Plane.
Automatic detection of the valve plane is important for separating the ventricles from the atria. The three-dimensional image at end-systole was subtracted from the end-diastolic image to detect the valve plane. End-diastole was defined as the first time bin. End-systole was defined as the time bin in which the absolute difference with respect to the end-diastolic frame was maximal. Only voxels with a count density greater than 50% of the maximum count density found at end-diastole were considered. In the subtraction image, the voxels at the borders of the ventricles become positive, whereas the voxels close to the atria turn out to be negative (Fig. 1). The valve plane was defined as the collection of voxels that satisfy the following condition:
![]() | (Eq. 1) |
|
![]() | (Eq. 2) |
This factor F was calculated using all voxels identified by their coordinates x, y, and z. pxyz+ and pxyz- represent the number of voxels with positive values positioned in front of or behind the valve plane, respectively. In contrast, nxyz- and nxyz+ are the number of voxels with negative values located behind and in front of the valve plane, respectively. Once F was maximized, voxels in front of the valve plane were identified as ventricular, whereas voxels behind the valve plane were recognized as auricular. Voxels behind the valve plane were not considered for further processing.
Detection of Interventricular Septum.
The septum was detected so that the right ventricle could be separated from the left ventricle. All short-axis slices at end-diastole containing voxels that were on the ventricular side of the valve plane were summed to obtain a two-dimensional best-septal-view image of the left and right ventricles. Again, only voxels with a count density greater than 50% of the maximum count density found at end-diastole were considered.
By summing, in this two-dimensional image, all pixels along the x-axis, a one-dimensional profile was generated (Fig. 2). The position of the local maximum (point A) in this profile was selected as the most likely location for the center of both ventricles. A profile along this position was created, and the local minimum was selected as the position of the septum (point B). The local maximum positioned on the left cardiac side of this local minimum was selected as the center of the left ventricle (point C).
|
Delineation of Left Ventricular Cavity at End-Diastole.
The segmentation started by applying a 30% threshold to the end-diastolic time bin. Only voxels having count densities higher than this threshold were considered.
First, all local maximums in the image were detected, that is, all voxels that had higher count densities than their 26 neighbors in the three-dimensional image space. These voxels were called roots. All roots were labeled with a unique index. Voxels that were not labeled as roots were linked to a root by reaching it using the steepest slope path. The steepest slope path was defined as the collection of voxels having the largest difference from one of its 26 neighbors. All voxels reaching the same root were labeled with the index of that root. Figures 3A and B show, respectively, an original and a segmented end-diastolic horizontal long-axis slice. Finally, segments with roots on the ventricular side of the valve plane and on the left side of the detected septum were identified as the left ventricle (Fig. 3C).
|
![]() | (Eq. 3) |
Delineation of Left Ventricular Cavity for Other Time Bins.
The delineation of the left ventricular cavity in the remaining time bins was done using the segmentation procedure and applying the threshold optimized for the end-diastolic delineation and corrected for a possible partial-volume effect. Accurate measurement of smaller ventricular volumes during the cardiac cycle can become difficult because of the limited resolution of the gamma camera, which leads to a partial-volume effect and results in an overestimation of the smaller volumes (8). To overcome this problem, the threshold initially applied was increased to compensate for the influence on the smaller volumes. The following formula was used to relate the loss of peak count density to the threshold required:
![]() | (Eq. 4) |
TED and TCorr represent the thresholds applied at end-diastole and after correction for a possible partial-volume effect, respectively. PED is the loss of peak count density, at a specific time bin, proportional to end-diastole and expressed as a percentage. This expression was derived after measuring a series of 99mTc-filled spheres of known volumes. After the segmentation of a specific time bin, segments with roots within the left ventricle delineated at end-diastole were identified as the left ventricle for that time bin.
Calculation of Left Ventricular Volumes and Ejection Fraction.
Once the left ventricle was delineated for each time bin, the volumes of the cavity were calculated by summing all voxels. To provide continuity to the measured volumes during the cardiac cycle, a temporal fit based on a Fourier series, keeping only the first two harmonics, was applied to the calculated volumes. End-diastolic volume (EDV) and end-systolic volume (ESV) were defined as the largest and smallest volumes, respectively, on the fitted volumetime curve. Ejection fraction was computed as (EDV - ESV)/EDV and expressed as a percentage.
Inter- and Intraobserver Reproducibility
The inter- and intraobserver reproducibility of the new algorithm, used to calculate LVEF from GBPT studies, was determined as usual. Interobserver reproducibility was determined by processing the GBPT studies two times by two observers. Intraobserver reproducibility was determined by processing the GBPT studies two times by the same observer at a 3-mo interval.
Effect of Reconstruction Filter
The effect of the reconstruction filter on the LVEF measured from GBPT was evaluated on 20 randomly selected studies reconstructed using a Butterworth filter with cutoff frequencies of 0.25, 0.30, 0.35, 0.40, and 0.45 cycles per centimeter. Each dataset was then reoriented with the same angles to eliminate variability caused by reorientation.
Statistical Analysis
Linear regression analysis was used to compare the calculated LVEFs. Slope, intercept, and correlation coefficient were calculated. The Bland-Altman method (17) was used to calculate the systematic and random errors between GBPT and PRNA. Differences between the repeated measurements were evaluated by the paired Student t test, and P < 0.05 was considered significant.
| RESULTS |
|---|
|
|
|---|
GBPT Versus PRNA
Figure 4A shows the relationship between LVEF measured with GBPT and LVEF measured with PRNA for the 50 patients with adequate delineation of the left ventricular cavity. The correlation coefficient between the two methods was 0.78 (slope = 0.94; intercept = 6.33). The mean difference (systematic error) between the GBPT and PRNA was +2.8%, with an SD (random error) of 8.8%. The t statistics indicated a P of 0.03, with mean values for the ejection fraction of 58% ± 14% and 55% ± 11% for GBPT and PRNA, respectively.
|
Intraobserver Reproducibility
Figure 5 shows the intraobserver reproducibility of the measured LVEF from GBPT. An excellent correlation between both measurements was observed (r = 0.94; slope = 0.91; intercept = 6.41). The systematic error was +1.1%, and the random error was 4.6%. No significant difference was observed between either series of measurements.
|
|
|
| DISCUSSION |
|---|
|
|
|---|
In this article, a new, fully automatic algorithm to compute LVEF from GBPT is presented. The new algorithm applies a three-dimensional segmentation technique combining threshold and local gradient methods to delineate the left ventricular cavity. The new algorithm is fast, has a high success rate, agrees well with conventional PRNA, and is reproducible.
The success rate of the algorithm was 94% among the 53 consecutive patients included in this study. The algorithm failed to correctly detect the left ventricular valve plane in two GBPT studies with poor count statistics, resulting in noisy images. The noise was further amplified when the end-diastolic and end-systolic images were subtracted. Left ventricular count density measured on the best-septal-view projection after summing all time frames was approximately 80 counts per pixel for these two patients, compared with 239 ± 62 counts per pixel for patients with adequate delineation of the valve plane. We therefore recommend that gated blood-pool studies contain at least 100 counts per pixel in the left ventricle. This goal can be achieved by increasing the dose in obese patients or, preferably, by adapting the time per stop after performing a short static acquisition in the best septal view just before starting the GBPT study.
Errors in locating the interventricular septum were observed for one patient with a very large left ventricle but a normal right ventricle. In this patient, the two ventricles merged into one and the algorithm was unable to separate them. Caution is therefore recommended when using this algorithm in patients with severe dilated cardiomyopathies involving either the left or the right ventricle.
LVEFs calculated on GBPT, although providing slightly higher values on average, agreed well with LVEFs calculated on PRNA. Although PRNA is probably not the ideal gold standard for measuring LVEF, this technique is still the most objective, reproducible, and widely used to measure LVEF in humans. That the correlation coefficient between the two techniques was only fair was not surprising, because PRNA is a planar imaging technique whereas GBPT is tomographic. The overlap between the left ventricle and the surrounding cardiac chambers (principally the left atrium) is eliminated with GBPT, and background correction is not necessary. Interestingly, the two methods compared better in patients with an LVEF less than 50% than in patients with an LVEF greater than 50%. No significant difference was seen between LVEF measured with GBPT and LVEF measured with PRNA in patients with an LVEF less than 50%. In contrast, LVEF measured with GBPT was significantly higher in patients with an LVEF greater than 50%. Higher LVEFs by GBPT than by PRNA have also been reported recently for studies by Mariano-Goulart et al. (15) and Van Kriekinge et al. (16).
In gated studies, counts have to be distributed over several time frames. For a given injected dose and acquisition time, the shorter the number of intervals the better the statistical quality of the individual time frames. Therefore, eight frames per cardiac cycle are usually preferred in gated SPECT. LVEFs measured with eight frames are, on average, 4% lower than 16-frame gated SPECT (1). This decrease is constant over a wide range of LVEFs. In this study, the predictable underestimation of LVEF was reduced by temporally fitting the calculated left ventricular volume curve during the cardiac cycle using the first two harmonics of the Fourier transform. LVEF was calculated using the maximum (EDV) and minimum (ESV) of the fitted curve.
The new algorithm is robust and reproducible. Because the algorithm is fully automatic, intra- and interobserver variability in measuring LVEF from GBPT occurred because of small variations in the reorientation of the reconstructed transverse tomograms along the long axis of the left ventricular cavity. However, the three-dimensional nature of the new algorithm minimizes this potential source of variability, as was confirmed by the excellent intra- and interobserver reproducibility. Importantly, the algorithm was also found to be little influenced by the reconstruction filter.
| CONCLUSION |
|---|
|
|
|---|
| FOOTNOTES |
|---|
For correspondence or reprints contact: Christian Vanhove, MSc, Division of Nuclear Medicine, University Hospital, Free University of Brussels (AZ VUB), 101 Laarbeeklaan, B-1090 Brussels, Belgium.
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
D. Mariano-Goulart, L. Dechaux, F. Rouzet, E. Barbotte, C. Caderas de Kerleau, M. Rossi, and D. Le Guludec Diagnosis of Diffuse and Localized Arrhythmogenic Right Ventricular Dysplasia by Gated Blood-Pool SPECT J. Nucl. Med., September 1, 2007; 48(9): 1416 - 1423. [Abstract] [Full Text] [PDF] |
||||
![]() |
P. De Bondt, T. Claessens, B. Rys, O. De Winter, S. Vandenberghe, P. Segers, P. Verdonck, and R. A. Dierckx Accuracy of 4 Different Algorithms for the Analysis of Tomographic Radionuclide Ventriculography Using a Physical, Dynamic 4-Chamber Cardiac Phantom J. Nucl. Med., January 1, 2005; 46(1): 165 - 171. [Abstract] [Full Text] [PDF] |
||||
![]() |
P. De Bondt, K. Nichols, S. Vandenberghe, P. Segers, O. De Winter, C. Van de Wiele, P. Verdonck, A. Shazad, A. H. Shoyeb, and J. De Sutter Validation of Gated Blood-Pool SPECT Cardiac Measurements Tested Using a Biventricular Dynamic Physical Phantom J. Nucl. Med., June 1, 2003; 44(6): 967 - 972. [Abstract] [Full Text] [PDF] |
||||
![]() |
G. A. Wright, S. Thackray, S. Howey, and J. G. Cleland Left Ventricular Ejection Fraction and Volumes from Gated Blood-Pool SPECT: Comparison with Planar Gated Blood-Pool Imaging and Assessment of Repeatability in Patients with Heart Failure J. Nucl. Med., April 1, 2003; 44(4): 494 - 498. [Abstract] [Full Text] [PDF] |
||||
![]() |
H. Okazawa, M. Takahashi, T. Hata, K. Sugimoto, Y. Kishibe, and T. Tsuji Quantitative Evaluation of Myocardial Blood Flow and Ejection Fraction with a Single Dose of 13NH3 and Gated PET J. Nucl. Med., August 1, 2002; 43(8): 999 - 1005. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| JOURNAL OF NUCLEAR MEDICINE TECHNOLOGY | THE JOURNAL OF NUCLEAR MEDICINE |