Determination of the Optimum Filterfunction for SPECT Imaging
Total Page:16
File Type:pdf, Size:1020Kb
BASIC SCIENCES Determination of the Optimum FilterFunction for SPECT Imaging DavidR. Gilland,BenjaminM. W.Tsui, WilliamH. McCartney,J. RandolphPerry, and Jan Berg Department ofRadiology and Curriculum in Biomedical Engineering, The University of North Carolina at Chapel Hill; and General Electric Company, Medical Systems Group, Milwaukee, Wisconsin An observer study was performed in order to evaluate several filters used in SPECT imaging. The filters were applied to the simulated projection data of a uniform activity density cylinder which contained a cold, spherical lesion, 2 cm in diameter. The data incorporated the effects of the detectorand scatterresponsefunctions,photonattenuation,and noise.Reconstructed transaxial images were used in 2AFC and ROC observer studies testing lesion detectability. In the 2AFC experiment, the Hanning filter scored lowest and did not show a optimum cutoff frequency. The Butterworth filter performed better and showed a well-defined optimum cutoff frequency at 0.15 cycles/pixel. The Metz filter performed as well as the optimum Butterworth but did not show an optimum power factor. In the ROC study, a high power Metz filter demonstrated an ROC curve of lower A@index and different shape from a lower power Metz filter and the optimum Butterworth filter. J Nucl Med 29:643—650, 1988 he problem of noise in single photon emission the optimum filter function and cutoff frequency for computed tomography (SPECT) is usually handled with simulated SPECT images of a simple phantom. the application of low-pass, digital filters designed to The criterion which we have used in judging the suppress high spatial frequencies where no signal is performance of a filter function is in terms of its effect present or where noise dominates the signal. Most on lesion detectability. We have performed observer SPECT filter functions allow the user to control the performance studies using computer simulated images degree of high frequency suppression by choosing a in order to quantify lesion detectability and to deter cutoff frequency, or similar filter parameter, which mine the optimum parameters for three filter types. determines where the filter rolls off to zero gain. The Subjects are shown images on a video monitor and are location of this cutoff frequency determines how the rated on their ability to detect a lesion. ifiter will affect both image noise and resolution. Low In order to make the simulated images resemble cutoff frequencies provide good noise suppression, but acquired clinical images as closely as possible, we have they can blur the image. Higher cutoff frequencies can modeled and incorporated into the simulation many of preserve resolution, but often suppress noise insuffi the physical phenomena ofan actual clinical liver study. ciently. There should exist an optimum cutoff fre The simulation incorporates the effects ofphoton atten quency for a particular filter function which compro uation, and the detector and scatter response functions. mises the trade-off between noise suppression and spa Also, noise levels are similar to those found in clinical tial resolution degradation. This optimum will depend liver studies. on factors such as the detector response function, the The results of this study can have important clinical spatial frequencies of the object, and the count density applications. SPECT filters can greatly affect the quality ofthe image. The purpose ofthis study was to determine of clinical images by their degree of smoothing. Deter mining the best filter and the proper degree of smooth ing can help to ensure the most accurate diagnosis. These results can also help speed image processing time ReceivedMar. 19, 1987;revision accepted Nov. 24, 1987. For reprints contact: David R. Gilland, Curriculum in Biomed since a proper filter function is often chosen clinically ical Engineering, Univ. of North Carolina at Chapel Hill, 152 by the tedious and time-consuming process of trial and MacNider Hall 202H, Chapel Hill, NC 27514. error. Volume29•Number5 •May1988 643 BACKGROUND select the cutoff frequency which determines the location of the filter rolloff. An example ofthis filter type is the Hanning Filter Operation filter. The Hanning filter is defined in the frequency domain The filtering ofSPECT images has been performed at three as, different stages of the reconstruction algorithm. Originally, @ H (is) = 0.5 + 0.5 cos (i-v/vs) O@ I @‘I (1) the smoothing filter was combined with the ramp filter which was applied one-dimensionallyin the transversedirection of = 0 otherwise. the projection data. However, when sagittal, coronal, or The cutoff frequency, v@,determines when the function oblique images are reconstructed from multiple slice, trans reaches zero gain. The frequency response of the Hanning axial data which has been smoothed in this way, streaking can filter is illustrated in Figure lA. occur in a direction parallel to the transaxial plane. These The Butterworth filter is another example of a smoothing streaks are the result ofpreferential smoothing in the transaxial filter and it is defined in the frequency domain as, direction. In order to overcome this problem, algorithms have been B(v)= n½ (2) developed which perform separate smoothing in the axial and I fl―' 1+1—I I transaxial directions. A one-dimensional smoothing filter is \Vc/ j applied first in the axial direction of the projection data and then in the transaxial direction along with the ramp filter In contrast to the Hanning, the Butterworth filter cutoff during the reconstruction. While this technique of separable frequency is defined as the point where the gain is down to smoothing in the axial and transaxial directions can eliminate 0.707. The shape ofthe Butterworth differs from the Hanning streaking, it has been shown that it does not result in a three in that this filter maintains a value close to one at low dimensional isotropic point response (1). Therefore, this tech frequencies and has a steeper rolloff. The parameter n deter nique is unsuitable for oblique tomography. mines the steepness of the rolloff. The frequency response of The solution to this problem has come about with the the Butterworth is illustrated in Figure lB. application ofa two-dimensional, radially symmetric smooth Another class of filterswhich have been used in SPECT ing function to the projection data before reconstruction. imaging is the enhancement and restoration ifiters (2—4). When this filter is combined with the ramp in the transaxial These filters attempt to recover the resolution lost in the direction, the resultant point response is isotropic (1). There detection process by exceeding unity gain over a desired fore, this algorithm is well-suited for the sagittal, coronal, and frequency band. Because of the domination of noise at high oblique reconstructions. In addition, this algorithm is superior frequencies,enhancement filters at some point must also in noise handling to the algorithms previously discussed (1). rolloffto zero gain. An example ofa enhancement filter is the Metz filter. This filter is defined in the frequency domain as, SPECr Filters — 1 [1 MTF2 (@)]X There are basically two types of filters that are used in M (v) = _____________ (3) SPECT imaging; the smoothing and enhancement filters. The MTF (v) smoothing filters are low-pass filters and allow the user to The Metz filter is a combination of deconvolution and Power Foctor —.— 2 ---3 I.00 —4 uc I.00 I.00 (cycles/pixel) ——.15 0 —.— 2 --- 5 4 .75 -- - 35 — .30 0 z —.5 4 w z z .50 4 .0 .2 3 .4 .5 .1 .2 .3 .4 .5 .0 .1 .2 .3 .4 .5 A CYCLES/PIXEL B CYCLES/PIXEL C CYCLES/PIXEL FiGURE 1 A:Frequencyresponseof HanningFilter.TheplotsshowtheHanningFilteratthreecutofffrequencies(va).Thecutoff frequency is the point at which the filter reaches zero gain. B: Frequency response of Butterworth Filter. The plots show the Butterworth Filter at two cutoff frequencies (vs)and two power factors (n). The cutoff frequency is the point at which the gain is down to 0.707. The power factor controls the steepness of the rolloff.C: Frequency response of Metz Filter. The plots show the Metz filter at three power factors. The power factor determines when the filter changes fromtheinverseMTFto a rolloffto zerogain.TheMTFwascalculatedbasedona sourcelocated25cmfromthe collimator face and at a depth of 12.5 cm into the scattering medium. 644 Gilland,Tsui,McCartneyetal The Journal of Nuclear Medicine smoothing filters. The parameter X determines the extent to represents the ratio of scattered to unscattered photons de which the filter follows the inverse MiT before rolling-off to tected and is a function of the depth of the source in the zero gain. The frequency response of the Metz filter is illus scattering medium. Monte Carlo studies (6) have demon trated in Figure 1C. strated that the scatter and the geometric components are . It has been reported previously that visually more pleasing related by a multiplicative function, images have resulted when a generalized exponential of the form, DRF, (v) = K (v) * DRF6 (v), (5) H (v) = exp (—i/S), (4) where K(v) is a function of the scattering conditions, i.e., energy window settings, detector energy resolution, source is used in place of the MTF in the Metz filter (2). The primary energy, and depth ofsource in the scattering medium. parameters P and S can be adjusted according the count An alternate representationof the detector response is the density and power spectrum of the image. This form of the use ofthe series-equivalent scatter response function (7) where Metz filter also produced images with a lower mean-squared the total response function is: error than the MTF form. DRFT (v) DRF0 (v) * DRF@' (v). (6) METhODS Monte Carlo studies(6) and experimental data (8) have shown that DRFs'@'(v)is essentially independent ofcollimation and is a function of only the scattering conditions.