Evaluation of scatter compensation methods by their effects on parameter estimation from SPECT projections

Three algorithms for scatter compensation in Tc-99m brain single-photon emission computed tomography (SPECT) were optimized and compared on the basis of the accuracy and precision with which lesion and background activity could be simultaneously estimated. These performance metrics are directly related to the clinically important tasks of activity quantitation and lesion detection, in contrast to measures based solely on the fidelity of image pixel values. The scatter compensation algorithms were (a) the Compton-window (CW) method with a 20% photopeak window, a 92-126 keV scatter window, and an optimized "k-factor," (b) the triple-energy window (TEW) method, with optimized widths of the photopeak window and the abutting scatter window, and (c) a general spectral (GS) method using seventeen 4 keV windows with optimized energy weights. Each method was optimized by minimizing the sum of the mean-squared errors (MSE) of the estimates of lesion and background activity concentrations. The accuracy and precision of activity estimates were then determined for lesions of different size, location, and contrast, as well as for a more complex Bayesian estimation task in which lesion size was also estimated. For the TEW and GS methods, parameters optimized for the estimation task differed significantly from those optimized for global normalized pixel MSE. For optimal estimation, the CW bias of activity estimates was larger and varied more (-2% to 22%) with lesion location and size than that of the other methods. The magnitude of the TEW bias was less than 7% across most conditions, although its precision was worse than that of CW estimates. The GS method performed best, with bias generally less than 4% and the lowest variance; its root-mean square (rms) estimation error was within a few percent of that achievable from primary photons alone. For brain SPECT, estimation performance with an optimized, energy-based, subtractive correction may approach that of an ideal scatter-rejection procedure.     

Performance of a novel collimator for high-sensitivity brain SPECT

We assessed improvements in performance in detection and estimation tasks due to a novel brain single photon computed tomography collimator. Data were acquired on the CeraSPECT scanner using both new and standard collimators. The new variable focusing collimator SensOgrade samples the projections unequally, with central regions more heavily represented, to compensate for attenuation of counts from central brain structures. Furthermore, it utilizes more of the cylindrical crystal surface. Two phantom studies were performed. The first phantom was a 21-cm-diameter cylindrical background containing nine spheres ranging from 0.5 to 5 cm3 in volume. 99mTc sphere to background activity ratio was 10:1. Twenty-nine 10-min datasets were acquired with each collimator. The second phantom was the Radiology Support Devices (Long Beach, CA) striatal phantom with striatal-background ratios of 10:1 on the left and 5:1 on the right. Twenty-nine 4-min datasets were acquired with each collimator. Perfusion imaging using 99mTc-HMPAO was also performed in three healthy volunteers using both collimators under identical simulations. Projections were reconstructed by filtered backprojection with an unwindowed ramp filter. The nonprewhitening matched filter signal-to-noise ratio (NPW-SNR) was computed as a surrogate for human performance in detecting spherical lesions. Sphere activity concentration, radius, and location coordinates were simultaneously estimated by fitting images to an assumed model using an iterative nonlinear algorithm. Resolution recovery was implicit in the estimation procedure, as the point spread function was incorporated into the model. NPW-SNR for sphere detection was 1.5 to 2 times greater with the new collimator; for the striatal phantom the improvement in SNR was 54%. The SNR for estimating sphere activity concentration improved by 46 to 89% for spheres located more than 5 cm from the phantom center. Images acquired with the standard collimator were too noisy in the central regions to allow estimation of sphere activity. In 99mTc-HMPAO human studies, SNR was improved by 21 to 41% in the cortex, 66% in the basal ganglia, and 74% in the thalamus. The new collimator leads to substantially improved detection and estimation performance throughout the brain. The higher sensitivity will be particularly important for dynamic imaging.     

The influence of backscatter material on 99mTc and 201Tl line source responses

SPECT projections are contaminated by scatter, resulting in reduced image contrast and quantitative errors. When tissue is present behind the source, some of the detected photons backscatter via this tissue. Particularly in dual-isotope SPECT and in combined emission-transmission SPECT, backscatter constitutes a major part of the down-scatter contamination in lower-energy windows. In this paper, the effects of backscatter material were investigated. Planar images of 99mTc and 201Tl line sources between varying numbers of Perspex slabs were analysed using the photopeak windows and various scatter windows. In the 99mTc photopeak window no significant change in total counts due to backscatter material was measured. In the 201Tl photopeak window an increase of about 10% in total counts was observed. In the scatter windows an even more explicit influence of backscatter material was measured. For instance, at a forward depth of 10 cm, total counts of a 99mTc source detected in the 72 keV window eventually doubled with increasing backscatter material, compared with the situation without backscatter material. The backscatter contribution plateaued when more than 5-10 cm of scatter material was placed behind the source. In conclusion, backscatter should be taken into account, particularly in model-based down-scatter correction methods in dual-isotope SPECT and combined emission-transmission SPECT.     

Importance of the choice of the collimator for the detection of small lesions in scintimammography: a phantom study

99mTc methoxyisobutylisonitrile planar scintimammography (SMM) is mostly performed using low-energy high-resolution (LEHR) parallel collimators. We studied whether using a different collimator could improve the detection of small (< 1.5 cm) lesions for which SMM sensitivity is poor. Thirty four breast phantom configurations were considered, either with hot spheres simulating lesions or without any spheres. For each configuration, four planar acquisitions were performed using LEHR, low-energy ultra high-resolution (LEUHR), high-resolution fan-beam (HRFB) and ultra high-resolution fan-beam (UHRFB) collimators. Images corresponding to the 20% and 10% energy windows and to the Jaszczak subtraction were calculated. A database including 156 borderline images was derived. After training, 10 observers scored the images for the presence of a sphere. The performances in sphere detection were studied using receiver operating characteristic (ROC) analysis. For all types of image, the area under the ROC curve was highest with the UHRFB collimator and lowest with either the LEUHR or the HRFB collimator. For the 10% energy window images conventionally used in SMM, the detection sensitivities averaged 91%, 73%, 60% and 55% for the UHRFB, LEHR, HRFB and LEUHR collimators respectively, for the same specificity of 64%. We conclude that detection of small tumours in planar SMM might be significantly improved by using a UHRFB collimator instead of an LEHR collimator.     

Implementation and quantitative evaluation of analytical methods for attenuation correction in SPECT: a phantom study.

Three differing exact methods of inverting the two-dimensional (2D) exponential Radon transform were implemented and evaluated quantitatively with a phantom study. The phantom had the shape of a pie-chart divided into six cavities, each 480 ml in volume and 10 cm in height, that were symmetrically positioned in a cylinder that was 20 cm in diameter and 10 cm in height. This phantom tests for linearity between true activity concentration and measured activity concentration, and it is denoted as a linearity phantom in the present study. Each cavity contained a different concentration of a homogeneous solution of 99mTc (74, 148, 222, 296, 370 and 444 kBq ml(-1)). Data acquisition was performed with two energy windows: a 20% photopeak energy window set symmetrically over the 140 keV of 99mTc and a secondary 5% energy window set over the 122 keV peak. We optimized a triple-energy window scatter correction method for a gamma camera-collimator system to obtain accurate scatter-corrected projections. A circular ROI 3 cm in diameter was identified over each cavity region, and count density (counts per pixel) was calculated. This value was converted to activity concentration (kBq ml(-1)) using a cross-calibration coefficient between SPECT counts and the gamma well counter. The relation between true activity (x) and measured activity concentration (y) was fitted to a line using the least-squares method. Regression lines were y = 0.63 + 1.0255x (R2 = 0.9987), y = -2.62 + 1.0278x (R2 = 0.9995), and y = 0.092 + 1.0241x (R2 = 0.9989) for the Bellini, Inouye and Metz-Pan methods respectively. In another phantom study using two different types of phantoms, contrast of a cold region in the two was 96% and 101% for all three methods. Combined optimized scatter correction and analytical attenuation correction methods achieve good accuracy in quantification of activity distribution with a uniform attenuating medium.     

Improved activity estimation with MC-JOSEM versus TEW-JOSEM in 111In SPECT

We have previously developed a fast Monte Carlo (MC)-based joint ordered-subset expectation maximization (JOSEM) iterative reconstruction algorithm, MC-JOSEM. A phantom study was performed to compare quantitative imaging performance of MC-JOSEM with that of a triple-energy-window approach (TEW) in which estimated scatter was also included additively within JOSEM, TEW-JOSEM. We acquired high-count projections of a 5.5 cm3 sphere of 111In at different locations in the water-filled torso phantom; high-count projections were then obtained with 111In only in the liver or only in the soft-tissue background compartment, so that we could generate synthetic projections for spheres surrounded by various activity distributions. MC scatter estimates used by MC-JOSEM were computed once after five iterations of TEW-JOSEM. Images of different combinations of liver/background and sphere/background activity concentration ratios were reconstructed by both TEW-JOSEM and MC-JOSEM for 40 iterations. For activity estimation in the sphere, MC-JOSEM always produced better relative bias and relative standard deviation than TEW-JOSEM for each sphere location, iteration number, and activity combination. The average relative bias of activity estimates in the sphere for MC-JOSEM after 40 iterations was -6.9%, versus -15.8% for TEW-JOSEM, while the average relative standard deviation of the sphere activity estimates was 16.1% for MC-JOSEM, versus 27.4% for TEW-JOSEM. Additionally, the average relative bias of activity concentration estimates in the liver and the background for MC-JOSEM after 40 iterations was -3.9%, versus -12.2% for TEW-JOSEM, while the average relative standard deviation of these estimates was 2.5% for MC-JOSEM, versus 3.4% for TEW-JOSEM. MC-JOSEM is a promising approach for quantitative activity estimation in 111In SPECT.     

Signal to noise ratio based filter optimization in triple energy window scatter correction

Triple energy window (TEW) scatter correction estimates the contribution of scattered photons to the acquisition data by acquiring additional data through two narrow energy windows placed adjoined to the main (photopeak) energy window. The contribution is estimated by linear interpolation and then subtracted. Noise amplification is reduced by filtering both the photopeak scintigram and the scatter estimate. We have studied the filter settings of each filter using a physical phantom filled with a 201Tl-solution resulting in count densities comparable to clinical studies. The performance of order-8 Butterworth filters at different cut-off frequencies (CoFs) were compared based on signal to noise ratios (SNRs). The highest SNRs were obtained when the noisy scatter information was strongly filtered with the CoF less than or equal to 0.07 cycles/pixel (cpp). The best CoF for the filter of the photopeak image is object size dependent; smaller objects require a higher CoF. For objects with a size near the SPECT spatial resolution (approximately 15 mm) the optimal CoF is equal to 0.18 cpp. For larger objects (31.8 mm) the highest SNR was obtained with a CoF equal to 0.13 cpp. A CoF equal to 0.16 cpp is a good compromise for all objects with a diameter equal to the spatial resolution or larger. These results depend on the initial signal to noise ratio of the acquisition data and so on the count density.     

Detection efficiency of a high-pressure gas scintillation proportional chamber

The detection efficiency of a high-pressure, gas scintillation proportional chamber (GSPC), designed for medical imaging in the 30-150 keV energy range, has been investigated through measurement and Monte Carlo simulation. Measurements were conducted on a GSPC containing 4 atm of pure xenon separated from a hexagonal array of seven ultraviolet-sensitive photomultiplier tubes by 1.27-cm-thick fused-silica windows. Experimental measurements of the photopeak efficiency, fluorescence escape efficiency, and the energy collection efficiency were obtained. Results were also obtained for different photon energies and different values of temporal resolution. The measurements were compared with the results obtained from a Monte Carlo simulation designed specifically for investigating the imaging of low-energy photons (below 150 keV) with a gas-filled detector. The simulation was used to estimate photopeak efficiency, fluorescence escape efficiency, photopeak-to-fluorescence escape peak ratio, quantum interaction efficiency, energy collection efficiency, and local energy collection efficiency. The photopeak efficiency of the GSPC relative to that of a 3-in. (7.62-cm)-thick sodium iodide crystal was measured to be 0.284 +/- 0.001 at 60 keV and 0.057 +/- 0.001 at 140 keV. Of the 60-keV photons incident upon the detector, 70% +/- 4% interacted in the detector, with 28% +/- 1% being in the photopeak, as estimated both by experimentation and through the simulation. The maximum energy collection efficiency was found to be 65% at 60 keV, with 46% being deposited within 0.2 cm of the initial photon interaction. The information gained from this study is being used to design an optimized detector for use in specialized nuclear medicine studies.     

Optimization and performance evaluation of the microPET II scanner for in vivo small-animal imaging

MicroPET II is a newly developed PET (positron emission tomography) scanner designed for high-resolution imaging of small animals. It consists of 17,640 LSO crystals each measuring 0.975 x 0.975 x 12.5 mm3, which are arranged in 42 contiguous rings, with 420 crystals per ring. The scanner has an axial field of view (FOV) of 4.9 cm and a transaxial FOV of 8.5 cm. The purpose of this study was to carefully evaluate the performance of the system and to optimize settings for in vivo mouse and rat imaging studies. The volumetric image resolution was found to depend strongly on the reconstruction algorithm employed and averaged 1.1 mm (1.4 microl) across the central 3 cm of the transaxial FOV when using a statistical reconstruction algorithm with accurate system modelling. The sensitivity, scatter fraction and noise-equivalent count (NEC) rate for mouse- and rat-sized phantoms were measured for different energy and timing windows. Mouse imaging was optimized with a wide open energy window (150-750 keV) and a 10 ns timing window, leading to a sensitivity of 3.3% at the centre of the FOV and a peak NEC rate of 235,000 cps for a total activity of 80 MBq (2.2 mCi) in the phantom. Rat imaging, due to the higher scatter fraction, and the activity that lies outside of the field of view, achieved a maximum NEC rate of 24,600 cps for a total activity of 80 MBq (2.2 mCi) in the phantom, with an energy window of 250-750 keV and a 6 ns timing window. The sensitivity at the centre of the FOV for these settings is 2.1%. This work demonstrates that different scanner settings are necessary to optimize the NEC count rate for different-sized animals and different injected doses. Finally, phantom and in vivo animal studies are presented to demonstrate the capabilities of microPET II for small-animal imaging studies.     

Study of the performance of a novel 1 mm resolution dual-panel PET camera design dedicated to breast cancer imaging using Monte Carlo simulation

We studied the performance of a dual-panel positron emission tomography (PET) camera dedicated to breast cancer imaging using Monte Carlo simulation. The PET camera under development has two 10x 15 cm(2) plates that are constructed from arrays of I X 1 X 3 mm(3) LSO crystals coupled to novel ultra-thin (<200 Am) silicon position-sensitive avalanche photodiodes (PSAPD). In this design the photodetectors are configured "edge-on" with respect to incoming photons which encounter a minimum of 2 cm thick of LSO with directly measured photon interaction depth. Simulations predict that this camera will have 10-15% photon sensitivity, for an 8-4 cm panel separation. Detector measurements show approximately 1 mm(3) intrinsic spatial resolution, <12% energy resolution, and approximately 2 ns coincidence time resolution. By performing simulated dual-panel PET studies using a phantom comprising active breast, heart, and torso tissue, count performance was studied as a function of coincident time and energy windows. We also studied visualization of hot spheres of 2.5-4.0 mm diameter and various locations within the simulated breast tissue for 1 X 1 X 3 mm(3), 2 x 2 x 10 mm(3), 3 x 3 x 30 mm(3), and 4 X 4 X 20 mm(3) LSO crystal resolutions and different panel separations. Images were reconstructed by focal plane tomography with attenuation and normalization corrections applied. Simulation results indicate that with an activity concentration ratio of tumor:breast:heart:torso of 10:1:10:1 and 30 s of acquisition time, only the dual-plate PET camera comprising 1 X 1 X 3 mm(3) crystals could resolve 2.5 mm diameter spheres with an average peak-to-valley ratio of 1.3.     

Tc-99m attenuation coefficients in water-filled phantoms determined with gamma cameras

Quantitative imaging with gamma cameras requires compensation for attenuation of source photons. Some methods of compensation make use of a constant or average estimated attenuation coefficient mu. A value for mu of 0.15 cm-1 for 140.5-keV photons in water or tissue is commonly used. This value, however, neglects scattered photons which are detected within the energy window in gamma camera imaging. Values for mu of 0.12 cm-1 used in attenuation compensation of Tc-99m single-photon emission computed tomography scans of uniform cylindrical sources have been shown to give improved results compared with use of mu = 0.15 cm-1. In this study, gamma cameras and a multichannel pulse-height analyzer were used to determine effective values of mu for photons in water as a function of energy window. Two cylindrical water-filled phantoms, circular and elliptical, were used with a point source of Tc-99m at depths up to 18 cm. Energy data were integrated over the top half of the photopeak, and over 10%, 20%, and 30% windows centered on the photopeak. Attenuation curves were exponential for all photopeak windows with values of mu of 0.12 +/- 0.014 cm-1 for all windows up to 20% and 0.1 cm-1 for a 30% window. This study suggests that a value of mu of 0.11-0.12 cm-1 is, in fact, appropriate for use in attenuation compensations where an average is required.     

Hybrid scatter correction applied to quantitative holmium-166 SPECT

Ho-166 is a combined beta-gamma emitter of which the betas can be used therapeutically. From the 81 keV gammas of Ho-166, SPECT images can be obtained, which give opportunities to guide Ho-166 therapy. Accurate reconstruction of Ho-166 images is currently hampered by photopeak-scatter in the patient, down-scatter in the detector, collimator and patient caused by the 1.4 MeV photons and by bremsstrahlung. We developed and validated a method for quantitative SPECT of Ho-166 that involves correction for both types of scatter plus non-uniform attenuation correction using attenuation maps. Photopeak-scatter (S) is compensated for by a rapid 3D Monte Carlo (MC) method that is incorporated in ordered subset (OS) reconstruction of the emission data, together with simultaneous correction for attenuation (A) and detector response (D); this method is referred to as OS-ADS. Additionally, for correction of down-scatter, we use a 14 keV wide energy window centred at 118 keV (OS-ADSS). Due to a limited number of available energy windows, the same 118 keV energy window was used for down-scatter correction of the simultaneously acquired Gd-153 transmission data. Validations were performed using physical phantom experiments carried out on a dual-head SPECT system; Gd-153 transmission line sources were used for acquiring attenuation maps. For quantitative comparison of OS-ADS and OS-ADSS, bottles filled with Ho-166 were placed in both a cylindrical phantom and an anthropomorphic thorax phantom. Both OS-ADS and OS-ADSS were compared with an ordered subset reconstruction without any scatter correction (OS-AD). Underestimations of about 20% in the attenuation map were reduced to a few per cent after down-scatter correction. The average deviation from the true activity contained in the bottles was +72% with OS-AD. Using OS-ADS, this average overestimation was reduced to +28% and with OS-ADSS the deviation was further reduced to 16%. With OS-AD and OS-ADS, these numbers were more sensitive to the choice of volumes of interest than with OS-ADSS. For the reconstructed activity distributions, erroneous background activity found with OS-AD was reduced by a factor of approximately 2 by applying OS-ADS and reduced by a factor of approximately 4 by applying OS-ADSS. The combined attenuation, photopeak-scatter and down-scatter correction framework proposed here greatly enhanced the quantitative accuracy of Ho-166 imaging, which is of the uppermost importance for image-guided therapies. It is expected that the method, with adapted window settings, also can be applied to other isotopes with high energy peaks that contaminate the photopeak data, such as I-131 or In-111.     

Improvement of the fluorine-18 fluorodeoxyglucose images in simultaneous F-18 FDG/Tc-99m collimated spect imaging.

Collimated F-18 FDG SPECT imaging has been shown to be an acceptable alternative to F-18 FDG PET imaging for evaluating injured but viable myocardium. Ultra-energy (UHE) imaging is usually performed in simultaneous F-18 FDG/Tc-99m MIBI studies. The main limitations of this technique are degradation of the Tc-99m MIBI images due to F-18 downscatter to the Tc-99m window, and loss of resolution in Tc-99m images caused by using a UHE rather than a low-energy collimator. The quality of F-18 images has not been addressed. In our clinical and phantom studies we have found that F-18 images are inferior to simultaneously acquired Tc-99m MIBI images. This paper compares two correction methods for F-18 FDG images in a realistic cardiac phantom study. One approach is subtractive scatter correction, which employs a third 410 keV energy window image to estimate scatter. The other approach is based on restoration. The phantom acquisition was performed with 7.2 MBq of F-18 and 22.2 MBq of Tc-99m injected into the left ventricular (LV) wall. Three inserts, 3 cm, 2 cm, and 1 cm in diameter, were placed in the LV wall to simulate infarcts. Circumferential profiles were drawn from three successive short-axis slices and compared with true phantom data. The differences were calculated as root-mean-square error (rmse). Scatter correction improved rmse only 4.5 +/- 0.3%, while restoration improved rmse 16.1 +/- 0.4%, when compared with raw data. The same differences, measured as rmse, were 9.5 +/- 0.5, 6.8 +/- 0.4, and 5.1 +/- 0.5 for raw, scatter corrected, and restored F-18 data, respectively, when compared with Tc-99m window 140 keV data. The amount of noise, measured as root-mean-square % (rms%) was 5.3 +/- 0.5% for the Tc-99m image, 4.9 +/- 0.7% for the F-18 restored image, 6.2 +/- 0.6% for the raw F-18 image, and 6.5 +/- 0.9% for the scatter corrected F-18 image. The contrast measured for 2 cm and 3 cm inserts was 0.17 +/- 0.07 and 0.26 +/- 0.06 for F-18 raw data, 0.19 +/- 0.08 and 0.29 +/- 0.06 for the scatter corrected F-18 image, and 0.28 +/- 0.06 and 0.43 +/- 0.07 for the restored F-18 image. The contrast was 0.20 +/- 0.07 and 0.46 +/- 0.05 for the Tc-99m 140 keV window image. The restoration approach provided F-18 images of better contrast and detectibility than uncorrected or scatter corrected F-18 images. Restored F-18 images match better with the simultaneously acquired Tc-99m images.     

Dual-energy approach to contrast-enhanced mammography using the balanced filter method: spectral optimization and preliminary phantom measurement

Dual-energy contrast agent-enhanced mammography is a technique of demonstrating breast cancers obscured by a cluttered background resulting from the contrast between soft tissues in the breast. The technique has usually been implemented by exploiting two exposures to different x-ray tube voltages. In this article, another dual-energy approach using the balanced filter method without switching the tube voltages is described. For the spectral optimization of dual-energy mammography using the balanced filters, we applied a theoretical framework reported by Lemacks et al. [Med. Phys. 29, 1739-1751 (2002)] to calculate the signal-to-noise ratio (SNR) in an iodinated contrast agent subtraction image. This permits the selection of beam parameters such as tube voltage and balanced filter material, and the optimization of the latter's thickness with respect to some critical quantity-in this case, mean glandular dose. For an imaging system with a 0.1 mm thick CsI:T1 scintillator, we predict that the optimal tube voltage would be 45 kVp for a tungsten anode using zirconium, iodine, and neodymium balanced filters. A mean glandular dose of 1.0 mGy is required to obtain an SNR of 5 in order to detect 1.0 mg/cm2 iodine in the resulting clutter-free image of a 5 cm thick breast composed of 50% adipose and 50% glandular tissue. In addition to spectral optimization, we carried out phantom measurements to demonstrate the present dual-energy approach for obtaining a clutter-free image, which preferentially shows iodine, of a breast phantom comprising three major components-acrylic spheres, olive oil, and an iodinated contrast agent. The detection of iodine details on the cluttered background originating from the contrast between acrylic spheres and olive oil is analogous to the task of distinguishing contrast agents in a mixture of glandular and adipose tissues.     

Comparison of NaI(T1), CdTe, and HgI2 surgical probes: effect of scatter compensation on probe performance

Spatial variation in the background source distribution makes tumor detection difficult for single-detector probes. Using a single energy window that brackets the photopeak helps discriminate against background events dominated by Compton scattering. Another approach is to use the information provided by an additional window in the Compton region. The performances of NaI(T1), CdTe, and HgI2 surgical probes have been compared under realistic simulations of a tumor-staging procedure using optimal single-sided energy windows and a two-window scatter-subtraction technique. Results showed that despite the differences in energy resolution of the detectors, the performances of the probes in a variable background were similar when optimal single energy windows were used. When the background variations were large, using information provided by a second window improved probe performance.     

Comparison of correction techniques for simultaneous 201Tl/99mTc myocardial perfusion SPECT imaging: a dog study

We compared two correction methods for simultaneous 201Tl/99mTc dual-isotope single-photon emission computed tomography (SPECT). Both approaches use the information from the third energy window placed between the photopeak windows of the 201Tl and 99mTc. The first approach, described by Moore et al, corrects only for the contribution of the 99mTc to the 201Tl primary 70 keV window. We developed the three-window transformation dual isotope correction method, which is a simultaneous cross-talk correction. The two correction methods were compared in a simultaneous 201Tl/99mTc sestamibi cardiac dog study. Three separate acquisitions were performed in this dog study: two single-isotope and one dual-isotope acquisition. The 201Tl single-isotope images were used as references. The total number of counts, and the contrast between the left ventricular cavity (LVC) and the myocardium, were used in 70 keV short axis slices as parameters for evaluating the results of the dual-isotope correction methods. Three consecutive short-axis slices were used to calculate averaged contrast and the averaged total number of counts. The total number of the counts was 667000+/-500 and 414500+/-400 counts for the dual isotope (201Tl+/-99mTc) and single-isotope (201Tl-only) 70 keV images, respectively. The corrected dual-isotope images had 514700+/-700 and 368000+/-600 counts for Moore's correction and our approach, respectively. Moore's method improved contrast in the dual isotope 70 keV image to 0.14+/-0.03 from 0.11+/-0.02, which was the value in the 70 keV non-corrected dual-isotope image. Our method improved the same contrast to 0.22+/-0.03. The contrast in the 201Tl single-isotope 70 keV image was 0.28+/-0.02. Both methods improved the 70 keV dual-isotope images. However, our approach provided slightly better images than Moore's correction when compared with 201Tl-only 70 keV images.     

Count-based monitoring of Anger-camera spectra--local energy shifts due to rotation

This study reports on a spectral monitoring method in which (1) a small source fixed to the camera is used, (2) a narrow, offset window is set on the side of the photopeak, and (3) variations in count rate are measured to assess energy shifts in the vicinity of the source. For one camera model, the count rate drops from 100% to 76% over a rotation of 180 degrees, implying a local energy shift of 1.4 keV. Also looked for are local count-rate variations with rotation for (1) wide-symmetric, (2) 20%-symmetric, and (3) 10%-asymmetric windows. The last is in limited use to partially compensate for Compton scattering. The effects of background and time stability are assessed.     

Partial volume effect correction in PET using regularized iterative deconvolution with variance control based on local topology

Correcting positron emission tomography (PET) images for the partial volume effect (PVE) due to the limited resolution of PET has been a long-standing challenge. Various approaches including incorporation of the system response function in the reconstruction have been previously tested. We present a post-reconstruction PVE correction based on iterative deconvolution using a 3D maximum likelihood expectation-maximization (MLEM) algorithm. To achieve convergence we used a one step late (OSL) regularization procedure based on the assumption of local monotonic behavior of the PET signal following Alenius et al. This technique was further modified to selectively control variance depending on the local topology of the PET image. No prior 'anatomic' information is needed in this approach. An estimate of the noise properties of the image is used instead. The procedure was tested for symmetric and isotropic deconvolution functions with Gaussian shape and full width at half-maximum (FWHM) ranging from 6.31 mm to infinity. The method was applied to simulated and experimental scans of the NEMA NU 2 image quality phantom with the GE Discovery LS PET/CT scanner. The phantom contained uniform activity spheres with diameters ranging from 1 cm to 3.7 cm within uniform background. The optimal sphere activity to variance ratio was obtained when the deconvolution function was replaced by a step function few voxels wide. In this case, the deconvolution method converged in approximately 3-5 iterations for most points on both the simulated and experimental images. For the 1 cm diameter sphere, the contrast recovery improved from 12% to 36% in the simulated and from 21% to 55% in the experimental data. Recovery coefficients between 80% and 120% were obtained for all larger spheres, except for the 13 mm diameter sphere in the simulated scan (68%). No increase in variance was observed except for a few voxels neighboring strong activity gradients and inside the largest spheres. Testing the method for patient images increased the visibility of small lesions in non-uniform background and preserved the overall image quality. Regularized iterative deconvolution with variance control based on the local properties of the PET image and on estimated image noise is a promising approach for partial volume effect corrections in PET.     

Dual-window scatter correction and energy window setting in cerebral blood flow SPECT: a Monte Carlo study

The image quality in SPECT studies of the regional cerebral blood flow (rCBF) performed with 99mTc-HMPAO is degraded by scattered photons. The finite energy resolution of the gamma camera makes the detection of scattered photons unavoidable, and this is observed in the image as an impaired contrast between grey and white matter structures. In this work, a Monte Carlo simulated SPECT study of a realistic voxel-based brain phantom was used to evaluate the resulting contrast-to-noise ratio for a number of energy window settings, with and without the dual-window scatter correction. Values of the scaling factor k, used to obtain the fraction of scattered photons in the photopeak window, were estimated for each energy window. The use of a narrower, asymmetric, energy discrimination window improved the contrast, with a subsequent increase in statistical noise due to the lower number of counts. The photopeak-window setting giving the best contrast-to-noise ratio was found to be the same whether or not scatter correction was applied. Its value was 17% centred at 142 keV. At the optimum photopeak-window setting, the contrast was improved by using scatter correction, but the contrast-to-noise ratio was made worse.     

A hybrid scatter correction for 3D PET based on an estimation of the distribution of unscattered coincidences: implementation on the ECAT EXACT HR+

We implemented a hybrid scatter-correction method for 3D PET that combines two scatter-correction methods in a complementary way. The implemented scheme uses a method based on the discrimination of the energy of events (the estimation of trues method (ETM)) and an auxiliary method (the single scatter simulation method (SSSI) or the convolution-subtraction method (CONV)) in an attempt to increase the accuracy of the correction over a wider range of acquisitions. The ETM takes into account the scatter from outside the field-of-view (FOV), which is not estimated with the auxiliary method. On the other hand, the auxiliary method accounts for events that have scattered with small angles, which have an energy that cannot be discriminated from that of unscattered events using the ETM. The ETM uses the data acquired in an upper energy window above the photopeak (550-650 keV) to obtain a noisy estimate of the unscattered events in the standard window (350-650 keV). Our implementation uses the auxiliary method to correct the residual scatter in the upper window. After appropriate scaling, the upper window data are subtracted from the total coincidences acquired in the standard window, resulting in the final scatter estimate, after smoothing. In this work we compare the hybrid method with the corrections used by default in the 2D and 3D modes of the ECAT EXACT HR+ using phantom measurements. Generally, the contrast was better with the hybrid method, although the relative errors of quantification were similar. We conclude that hybrid techniques such as the one implemented in this work can provide an accurate, general-purpose and practical way to correct the scatter in 3D PET, taking into account the scatter from outside the FOV.