Quantitative performance of advanced resolution recovery strategies on SPECT images: evaluation with use of digital phantom models

Deep learning-based attenuation correction method in 99mTc-GSA SPECT/CT hepatic imaging: a phantom study

This study aimed to evaluate a deep learning-based attenuation correction (AC) method to generate pseudo-computed tomography (CT) images from non-AC single-photon emission computed tomography images (SPECTNC) for AC in 99mTc-galactosyl human albumin diethylenetriamine pentaacetic acid (GSA) scintigraphy and to reduce patient dosage. A cycle-consistent generative network (CycleGAN) model was used to generate pseudo-CT images. The training datasets comprised approximately 850 liver phantom images obtained from SPECTNC and real CT images. The training datasets were then input to CycleGAN, and pseudo-CT images were output. SPECT images with real-time CT attenuation correction (SPECTCTAC) and pseudo-CT attenuation correction (SPECTGAN) were acquired. The difference in liver volume between real CT and pseudo-CT images was evaluated. Total counts and uniformity were then used to evaluate the effects of AC. Additionally, the similarity coefficients of SPECTCTAC and SPECTGAN were assessed using a structural similarity (SSIM) index. The pseudo-CT images produced a lower liver volume than the real CT images. SPECTCTAC exhibited a higher total count than SPECTNC and SPECTGAN, which were approximately 60% and 7% lower, respectively. The uniformities of SPECTCTAC and SPECTGAN were better than those of SPECTNC. The mean SSIM value for SPECTCTAC and SPECTGAN was 0.97. We proposed a deep learning-based AC approach to generate pseudo-CT images from SPECTNC images in 99mTc-GSA scintigraphy. SPECTGAN with AC using pseudo-CT images was similar to SPECTCTAC, demonstrating the possibility of SPECT/CT examination with reduced exposure to radiation.

Evaluation of attenuation correction method for head holder in brain perfusion single-photon emission computed tomography

Head holder attenuation affects brain perfusion single-photon emission computed tomography (SPECT) image quality. Here, we proposed a head holder-attenuation correction (AC) method using attenuation coefficient maps calculated by Chang's method from CT images. Then, we evaluated the effectiveness of the head holder-AC method by numerical phantom and clinical cerebral perfusion SPECT studies. In the numerical phantom, the posterior counts were 10.7% lower than the anterior counts without head holder-AC method. However, by performing head holder-AC, the posterior count recovered by approximately 6.8%, approaching the true value. In the clinical study, the normalized count ratio was significantly increased by performing the head holder-AC method in the posterior-middle cerebral artery, posterior cerebral artery and cerebellum regions. There were no significant increases in other regions. The head holder-AC method can correct the counts attenuated by the head holder.

Evaluation of Collimators in a High-Resolution, Whole-Body SPECT/CT Device with a Dual-Head Cadmium-Zinc-Telluride Detector for 123I-FP-CIT SPECT

The study aim was to evaluate the adaptation of collimators to 123I-N-fluoropropyl-2b-carbomethoxy-3b-(4-iodophenyl)nortropane (123I-FP-CIT) dopamine transporter SPECT (DAT-SPECT) by a high-resolution whole-body SPECT/CT system with a cadmium-zinc-telluride detector (C-SPECT) in terms of image quality, quantitation, diagnostic performance, and acquisition time. Methods: Using a C-SPECT device equipped with a wide-energy, high-resolution collimator and a medium-energy, high-resolution sensitivity (MEHRS) collimator, we evaluated the image quality and quantification of DAT-SPECT for an anthropomorphic striatal phantom. Ordered-subset expectation maximization iterative reconstruction with resolution recovery, scatter, and attenuation correction was used, and the optimal collimator was determined on the basis of the contrast-to-noise ratio (CNR), percentage contrast, and specific binding ratio. The acquisition time that could be reduced using the optimal collimator was determined. The optimal collimator was used to retrospectively evaluate diagnostic accuracy via receiver-operating-characteristic analysis and specific binding ratios for 41 consecutive patients who underwent DAT-SPECT. Results: When the collimators were compared in the phantom verification, the CNR and percentage contrast were significantly higher for the MEHRS collimator than for the wide-energy high-resolution collimator (P < 0.05). There was no significant difference in the CNR between 30 and 15 min of imaging time using the MEHRS collimator. In the clinical study, the areas under the curve for acquisition times of 30 and 15 min were 0.927 and 0.906, respectively, and the diagnostic accuracies of the DAT-SPECT images did not significantly differ between the 2 times. Conclusion: The MEHRS collimator provided the best results for DAT-SPECT with C-SPECT; shorter acquisition times (<15 min) may be possible with injected activity of 167-186 MBq.

Impact of bone-equivalent solution density in a thoracic spine phantom on bone single-photon emission computed tomography image quality and quantification

This study aimed to evaluate the effects of dipotassium hydrogen phosphate (K₂HPO₄) solution density on single-photon emission computed tomography (SPECT) image quality and quantification. We used a JSP phantom containing six cylinders filled with K₂HPO₄ solutions of varying densities. Computed tomography (CT) was performed, and CT values and linear attenuation coefficients were measured. Subsequently, SPECT images of an SIM² bone phantom filled with ^99mTc with/without K₂HPO₄ solution were acquired using a SPECT/CT camera. The full width at half maximum (FWHM), percentage coefficient of variation (%CV), recovery coefficient, and standardized uptake value (SUV) were evaluated to investigate the impact of the K₂HPO₄ solution density. The CT values and linear attenuation coefficients increased with the K₂HPO₄ solution density. The CT values for cancellous and cortical bones were reflected by K₂HPO₄ solution densities of 0.15-0.20 and 1.50-1.70 g/cm³, respectively. FWHM values were significantly lower with the K₂HPO₄ solution than those with water alone (18.0 ± 0.9 mm with water alone, 15.6 ± 0.2 mm with 0.15 g/cm³ K₂HPO₄, and 16.1 ± 0.3 mm with 1.49 g/cm³ K₂HPO₄). Although the %CVs showed no significant differences, the recovery coefficients obtained with water alone tended to be slightly lower than those obtained with the K₂HPO₄ solution. The SUV obtained using the standard density of the K₂HPO₄ solution differed from that obtained using the optimized density. In conclusion, SPECT image quality and quantification depends on the presence and concentration of the bone-equivalent solution. The optimal bone-equivalent solution density should be used to evaluate the bone image phantoms.

Impact of patient body habitus on image quality and quantitative value in bone SPECT/CT

OBJECTIVE

The first edition of guidelines for standardization of bone single photon emission computed tomography (SPECT) imaging was published in 2017, and the optimization and standardization are widely promoted. To the purpose, clarification of the factors related to image quality and quantitative values and their influence are required. The present study aimed to clarify and optimize the influence of patient body habitus on image quality and quantitative values in bone SPECT/CT.

METHODS

National Electrical Manufacturers Association body phantom (S-size) and custom-made large body phantoms (M-size and L-size) that simulate the abdomens of Japanese patients weighing 60, 80, and 100 kg, were used. Each phantom was filled with ^99mTc-solutions of 108  and 18 kBq/mL for the hot spheres and background, respectively. Dynamic SPECT acquisition was performed for 6000 s (150 s /rotation × 40 rotation). The data were divided into six projection data and reconstructed each acquisition time (150, 300, 450, 600, 750, 900 s, and single projection 6000 s). Image quality was evaluated for contrast (Q_{H, 17 mm}), background noise (N_{B, 17 mm}), contrast-to-noise ratio (CNR), maximum standardized uptake value (SUV_{max, 17 mm}), and visual assessment for a 17 mm hot sphere.

RESULTS

Image quality in the 300 s acquisition showed that values of Q_{H, 17 mm}, CNR, and SUV_{max, 17 mm} decreased (-16.7%, -11.8%, and -11.3%) for M-size and (-28.2%, -30.1%, and -21.7%) for L-size compared with S-size, respectively. No significant difference was observed in N_{B, 17 mm} values. M-size and L-size required 1.2 and 2.3 times longer acquisition, to achieve same CNR as S-size. In visual assessment, 17 mm hot sphere could not be detected only in the L-size. When the Japanese bone SPECT guidelines criteria were applied in 600 s, the sphere could be detected between all phantoms.

CONCLUSIONS

Patient body habitus significantly affects image quality and decreases the quantitative value in bone SPECT/CT. For the optimization, extend acquisition time according to the patient body habitus is effective for image quality. And for the standardization, it is important to achieve imaging conditions that meet the Japanese bone SPECT guidelines criteria to ensure adequate detectability.

Improvement of Quantitative Single-Photon Emission Computed Tomography Image Quality by the New Step-and-Shoot Scan Mode

The step-and-shoot (SS) mode and continuous mode are currently used for single-photon emission computed tomography (SPECT) scan mode, and a new scan mode that combines both modes, step-and-shoot plus continuous (SSC) mode, was developed. It is expected to allow a shorter scan time and lower injected dose because the SSC mode is more sensitive than the SS mode. We confirmed the image quality of this scan mode, including various quantitative correction methods for scatter (SC), attenuation (AC), and resolution recovery (RR) in a phantom study and clinical case study. Image quality was evaluated by the count, contrast-to-noise ratio (CNR), and percent of the coefficient of variation (%CV). Independent of the correction methods, the count, CNR, and %CV of the SSC mode were superior to those of the SS mode. The ACSCRR was the best method, with a maximum increased rate of 66.4% in counts and 57.8% in CNR for the 13-mm sphere and 19.6% in CNR for other sphere sizes. The %CV for the SSC mode was the best for AC and ACRR, which was at 15.1%. With regards to attaining short bone SPECT scan time, the combination of the SSC mode and ACRR or ACSCRR demonstrated the best physical performance.

Optimization of image reconstruction method of cerebral blood flow perfusion imaging with digital CZT SPECT

PURPOSE

The purpose of this study is to evaluate the effects of filtered back projection (FBP), ordered subset expectation maximisation (OSEM), and different filters on cadmium zinc telluride single-photon emission computed tomography [CZT single-photon emission computed tomography (SPECT)] cerebral blood perfusion image quality to optimise the image reconstruction method.

METHODS

Under routine clinical conditions, tomographic imaging was performed on the phantom and patients. Image processing included image reconstruction using FBP and OSEM, and the filtering method used Butterworth (Bw) and Gaussian (Gs) filters. Visual and semi-quantitative parameters [integral uniformity, root mean square (RMS) noise and contrast and contrast-to-noise ratio (CNR)] were used to evaluate image quality to optimise image reconstruction parameters. One-way and two-way analysis of variance were used to process phantom and clinical data.

RESULTS

In the tomographic images of the phantom, the semi-quantitative analysis showed that the integral uniformity of FBP+Bw was better than that of OSEM+Bw and OSEM+Gs (P < 0.05), and that the RMS noise of FBP+Bw was lower than that of OSEM+Bw and OSEM+Gs (P < 0.001). The contrast of FBP+Bw and OSEM+Bw in the cold area diameter ≥2 cm group was higher than that of OSEM+Gs (P < 0.001), whereas the CNR of FBP+Bw was higher than that of OSEM+Bw and OSEM+Gs (P < 0.001); the contrast of OSEM+Bw cold area diameter <2 cm was higher than that of FBP+Bw (P < 0.01). The semi-quantitative analysis results of the clinical images were consistent with the phantom's.

CONCLUSION

In CZT SPECT cerebral blood flow perfusion imaging, it is suggested that the image postprocessing method of FBP+Bw (fc = 0.40; n = 10) should be used routinely in clinical application, and if there are uncertain small lesions in the processed image, it is suggested to use the reconstruction method of OSEM+Bw (EM-equivalent iterations = 60; fc = 0.45; n = 10) instead.

[Improvement of Standardized Uptake Value Accuracy in the 99mTc Body SPECT and SPECT/CT: Optimization of the Phantom for Calculating Becquerel Calibration Factor and Correction Method]

PURPOSE

The purpose of this study was to evaluate the best phantom for calculating the becquerel calibration factor (BCF) and correction method to obtain the improvement of standardized uptake value (SUV) accuracy in both single photon emission computed tomography (SPECT) and SPECT/CT.

METHOD

A SPECT/CT scanner was used in this study. BCFs were calculated using four phantoms with different cross sections including National Electrical Manufacturers Association International Electrotechnical Commission body phantom (NEMA IEC body phantom) filled with ^99mTcO₄^-, and five correction methods were used for reconstruction. SUVs were calculated by the NEMA IEC body phantom and pediatric phantom in house with these BCFs. We then measured SUV_mean in the background region of the NEMA IEC body phantom, SUV_max and SUV_peak of the 37-mm-diameter sphere.

RESULTS

In the SPECT scanner, SUV_mean and SUV_max measured 1.04 and 4.02, respectively, in the case of BCF calculation and SUV measurement using NEMA IEC body phantoms without corrections. In the SPECT/CT scanner, SUV_mean with CT attenuation correction (AC) was in agreement with the theoretical values using each phantom. SUV_max showed the same trend.

CONCLUSION

In the SPECT scanner, it is possible to obtain a highly accurate SUV by using a phantom that matches the size of the subject for BCF calculation and without correction. In the SPECT/CT scanner, highly accurate SUVs can be obtained by using CT-based attenuation correction, and these values do not depend on the size of the BCF calculation phantom.

Comparison of Myocardial Ischemia Detection Between Semiconductor and Conventional Anger-type Three-detector SPECT

Objective: Although semiconductor single-photon emission computed tomography (D-SPECT) has been used for myocardial perfusion imaging, few studies have compared its ability to detect myocardial ischemia with that of 3-detector SPECT (GCA9300R). This study used invasive coronary angiography to determine whether the detectability of myocardial ischemia differs between D-SPECT and GCA9300R. Materials and methods: This study included 24 patients who were assessed by coronary angiography within 60 days of myocardial perfusion D-SPECT and GCA9300R. Two nuclear medicine physicians interpreted myocardial perfusion D-SPECT and GCA9300R images with five grades of confidence, then defined regions of ischemia on polar maps. The gold standard was determined by another nuclear cardiology specialist based on integrated assessment of the coronary angiography findings and other clinical information derived from medical charts. The concordance rate and the Cohen kappa (κ) between D-SPECT and GCA9300R were calculated. Results: The sensitivity, specificity, negative and positive predictive values, and the accuracy of patient-based diagnoses were 66.7%, 91.7%, 89.2%, 72.8%, and 85.5%, respectively, for GCA9300R, and 83.3%, 83.3%, 93.7%, 62.4%, and 83.3%, respectively, for D-SPECT. Interpretations of ischemia did not uncover any significant differences between D-SPECT and GCA9300R. The Cohen κ values of D-SPECT and GCA9300 agreed substantially, moderately and marginally for the left circumflex coronary artery (LCX) (0.68), right coronary artery (RCA) (0.43), and left anterior descending coronary artery (LAD) (0.39), respectively. Conclusions: The detectability of myocardial ischemia is comparable between D-SPECT and GCA9300R. Sensitivity is better for D-SPECT than GCA9300R. However, false-positive D-SPECT findings, especially in the apex and inferior wall should be interpreted with caution.

Performance of SwiftScan planar and SPECT technology using low-energy high-resolution and sensitivity collimator compared with Siemens SPECT system

PURPOSE

A new low-energy high-resolution-sensitivity (LEHRS) collimator was developed by General Electric (GE) Healthcare. SwiftScan planar and single photon emission computed tomography (SPECT) systems using LEHRS collimator were developed to achieve the low-dose and/or short-time acquisition. We demonstrated the performance of SwiftScan planar and SPECT system with LEHRS collimator using phantoms.

METHODS

Line source, cylindrical and flat plastic dish phantoms were used to evaluate the performance of planar and SPECT images for four patterns of Siemens LEHR, GE LEHR, GE LEHRS and SwiftScan using two SPECT-CT scanners. Each phantom was filled with 99mTc solution, and the spatial resolution, sensitivity and image uniformity were calculated from the planar and SPECT data.

RESULTS

The full-width at half maximum (FWHM) values as a system spatial resolution of Siemens LEHR, GE LEHR and GE LEHRS were approximately 7.4 mm. GE LEHRS showed a lower FWHM value by increasing the blend ratio in Clarity2D processing. The system sensitivity of GE LEHRS increased by approximately 30% compared with that of GE LEHR and was similar to that of Siemens LEHR. The FWHM values of SPECT with an filtered back projection (FBP) method were approximately 10.3 mm. The FWHM values of the ordered subset expectation maximization (OSEM) method were better with an increase in iteration values. The differential uniformities of Siemens LEHR, GE LEHR, GE LEHRS and GE SwiftScan using the FBP method were approximately 15.1%. The differential uniformity of OSEM method was higher with an increase in the iteration value.

CONCLUSION

The SwiftScan planar and SPECT have a high sensitivity while maintaining the spatial resolution compared with the conventional system.

Multicenter Study of Quantitative SPECT: Reproducibility of 99mTc Quantitation Using a Conjugated-Gradient Minimization Reconstruction Algorithm

This multicenter study aimed to determine the reproducibility of quantitative SPECT images reconstructed using a commercially available method of ordered-subset conjugate-gradient minimization. Methods: A common cylindric phantom containing a 100 kBq/mL concentration of ^99mTc-pertechnetate solution in a volume of 7 L was scanned under standard imaging conditions at 6 institutions using the local clinical protocol of each. Interinstitutional variation among the quantitative SPECT images was evaluated using the coefficient of variation. Dose calibrator accuracy was also investigated by measuring the same lot of commercially available ^99mTc vials at each institution. Results: The respective radioactivity concentrations under standard and clinical conditions ranged from 95.71 ± 0.60 (mean ± SD) to 108.35 ± 0.36 kBq/mL and from 96.78 ± 0.64 to 108.49 ± 0.11 kBq/mL, respectively. Interinstitutional variation in radioactivity concentration was 4.20%. The bias in the radioactivity concentrations in SPECT images was associated with the accuracy of the dose calibrator at each institution. Conclusion: The reproducibility of the commercially available quantitative SPECT reconstruction method is high and comparable to that of PET, for comparatively large (∼7 L), homogeneous objects.

Experimental evaluation of the GE NM/CT 870 CZT clinical SPECT system equipped with WEHR and MEHRS collimator

Purpose

A high‐energy‐resolution whole‐body SPECT‐CT device (NM/CT 870 CZT; C‐SPECT) equipped with a CZT detector has been developed and is being used clinically. A MEHRS collimator has also been developed recently, with an expected improvement in imaging accuracy using medium‐energy radionuclides. The objective of this study was to compare and analyze the accuracies of the following devices: a WEHR collimator and the MEHRS collimator installed on a C‐SPECT, and a NaI scintillation detector‐equipped Anger‐type SPECT (A‐SPECT) scanner, with a LEHR and LMEGP.

Methods

A line phantom was used to measure the energy resolutions including collimator characteristics in the planar acquisition of each device using ^99mTc and ¹²³I. We also measured the system's sensitivity and high‐contrast resolution using a lead bar phantom. We evaluated SPECT spatial resolution, high‐contrast resolution, radioactivity concentration linearity, and homogeneity, using a basic performance evaluation phantom. In addition, the effect of scatter correction was evaluated by varying the sub window (SW) employed for scattering correction.

Results

The energy resolution with ^99mTc was 5.6% in C‐SPECT with WEHR and 9.9% in A‐SPECT with LEHR. Using 123I, the results were 9.1% in C‐SPECT with WEHR, 5.5% in C‐SPECT with MEHRS, and 10.4% in A‐SPECT with LMEGP. The planar spatial resolution was similar under all conditions, but C‐SPECT performed better in SPECT acquisition. High‐contrast resolution was improved in C‐SPECT under planar condition and SPECT. The sensitivity and homogeneity were improved by setting the SW for scattering correction to 3% of the main peak in C‐SPECT.

Conclusion

C‐SPECT demonstrates excellent energy resolution and improved high‐contrast resolution for each radionuclide. In addition, when using 123I, careful attention should be paid to SW for scatter correction. By setting the appropriate SW, C‐SPECT with MEHRS has an excellent scattered ray removal effect, and highly homogenous imaging is possible while maintaining the high‐contrast resolution.

Optimization of Number of Iterations as a Reconstruction Parameter in Bone SPECT Imaging Using a Novel Thoracic Spine Phantom

The aim of this study was to optimize the number of iterations in bone SPECT imaging using a novel thoracic spine phantom (ISMM phantom). Methods: The quality and quantitative accuracy of bone SPECT images were evaluated by changing the number of iterations and the size of the hot spot in the phantom. True SUVs in the vertebra, tumor, and background parts were 9.8, 52.2, and 1.0, respectively. The phantom image was reconstructed using the ordered-subset expectation-maximization algorithm with CT-based attenuation correction, scatter correction, and resolution recovery; the number of ordered-subset expectation-maximization subsets was fixed at 10, with iterations ranging from 1 to 40. Full width at half maximum, percentage coefficient of variation, contrast ratio for the sphere and background (contrast), and recovery coefficient were evaluated as a function of the number of iterations for a given number of subsets (10) using the reconstructed images. In addition, SUV_max, SUV_peak, and SUV_mean were calculated with various numbers of iterations for each sphere (13, 17, 22, and 28 mm) simulating a tumor. Results: Full width at half maximum decreased as the number of iterations was increased, and full width at half maximum converged uniformly when the number of iterations exceeded 10. The percentage coefficient of variation increased as the number of iterations was increased. Recovery coefficient decreased with decreasing sphere size. Contrast and all SUVs increased as the number of iterations was increased, and contrast and all SUVs converged uniformly when the number of iterations exceeded 5 and 10, respectively, for all sphere sizes. When the SUV was defined as the converged value for 10 iterations in the 28-mm sphere, the converged values of SUV_max, SUV_peak, and SUV_mean were 75.1, 66.5, and 55.6, respectively. The relative error in the converged values for SUV_max, SUV_peak, and SUV_mean were 43.8%, 27.3%, and 7.2% of the true value (52.2); all SUVs were overestimated. Conclusion: Using a thoracic spine phantom to evaluate the optimal reconstruction parameters in bone SPECT imaging, we determined the optimal number of iterations for 10 subsets to be 10.

Imaging technology for myocardial perfusion single-photon emission computed tomography 2018 in Japan

AIM

Recently, nuclear cardiology has dramatically advanced by a new technology development such as the device, short-term acquisition system, image reconstruction algorithm and image analysis. Although these innovations have been gradually employed in routine examinations, we did not investigate the current use of image acquisition, image reconstruction, and image analysis with myocardial perfusion single-photon emission computed tomography (MPS). We investigated the current status of MPS imaging technology in Japan.

METHODS

We carried out a survey using a Web-based questionnaire system, the opening of which was announced via e-mail, and it was available on a website for 3 months. We collected data on the current use of MPS with ²⁰¹Tl and/or ^99mTc agents with respect to routine protocols, image acquisition, image reconstruction, and image analysis.

RESULTS

We received responses to the Web-based questionnaire from 178 and 174 people for ^99mTc and ²⁰¹Tl MPS, respectively. The routine protocols of MPS of stress-rest and rest-stress MPS on 1-day protocols with ^99mTc were 41.2% and 14.5%, respectively, and the rest-only scan response rate was 23.7%, whereas that of ²⁰¹Tl MPS was 65.9% with stress-rest MPS, 19.0% with rest-only MPS, and 10.9% with stress-rest MPS adding a rest scan 24 h after injection. The filtered back projection (FBP) method is most commonly used image reconstruction method, yielding 70.5% for ^99mTc MPS and 76.8% for ²⁰¹Tl MPS, including combined FBP and ordered subset expectation maximization method. The results for no-correction (NC) images were 49.2% with ^99mTc MPS and 55.2% with ²⁰¹Tl MPS including the response of NC and combined attenuation correction (AC) and scatter correction (SC) (i.e., ACSC) images. The AC or ACSC images of ^99mTc and ²⁰¹Tl were provided by 30-40% of the institutions surveyed.

CONCLUSIONS

We investigated the current status of MPS imaging technology in Japan, and found that although the use of various technical developments has been reported, some of these technologies have not been utilized effectively. Hence, we expect that nuclear medicine technology will be used more effectively to improve diagnosis.

Impact of iterative reconstruction with resolution recovery in myocardial perfusion SPECT: phantom and clinical studies

The corrections of photon attenuation, scatter, and depth-dependent blurring improve image quality in myocardial perfusion single-photon emission computed tomography (SPECT) imaging; however, the combined corrections induce artifacts. Here, we present the single correction method of depth-dependent blurring and its impact for myocardial perfusion distribution in phantom and clinical studies. The phantom and clinical patient images were acquired with two conditions: circular and noncircular orbits of gamma cameras yielded constant and variable depth-dependent blurring, respectively. An iterative reconstruction with the correction method of depth-dependent was used to reconstruct the phantom and clinical patient images. We found that the single correction method improved the robustness of phantom images whether the images contained constant or variable depth-dependent blurring. The myocardial perfusion databases generated from 72 normal patients exhibited uniform perfusion distribution of whole myocardium. In summary, the single correction method of depth-dependent blurring with iterative reconstruction is helpful for myocardial perfusion SPECT.

Evaluation of edge-preserving and noise-reducing effects using the nonlinear diffusion method in bone single-photon emission computed tomography

OBJECTIVES

This study aims to evaluate nonlinear diffusion (NLD) processing to smoothen images while suppressing resolution degradation in single-photon emission computed tomography (SPECT) images. Phantom data were used for NLD method optimization. The resultant optimal settings were used for NLD processing of clinical images.

MATERIALS AND METHODS

Tc was used to simulate tumors and normal soft tissues. Using the data collected, images were reconstructed. Images were processed using various k values and iteration. The background region's coefficient of variation (CV) was determined, and the effects of parameters on image properties were examined. NLD-processed images with optimal parameters were compared with Butterworth (BW)-filtered and nine-point smoothing (SM)-processed images to evaluate smoothing filter properties in real and frequency space. Receiver operating characteristic curve analysis was carried out on NLD-processed and BW048-processed bone SPECT images.

RESULTS

From CVs in background, with NLD, increased k value and iteration led to a low CV, indicating enhanced smoothing effect. At k=0.9, a strong noise-reducing effect with less iteration was achieved. Contrasts and recovery coefficients of NLD were the highest. The visual score for SPECT image quality was significantly higher with NLD than with BW048, BW090, and SM. In the low-frequency and high-frequency ranges, BW048, BW090, and NLD showed similar signal strengths and NLD and BW090 showed high signal strength, respectively. SM processing reduced the signal strength at all frequency ranges. On receiver operating characteristic analysis, noise reduction using NLD processing enhanced diagnostic performance than with the use of BW processing.

CONCLUSION

NLD processing of bone SPECT images using optimized parameters enabled smoothing with less resolution degradation.

Clinical evaluation of General Electric new Swiftscan solution in bone scintigraphy on NaI-camera: A head to head comparison with Siemens Symbia

PURPOSE

The General Electric (GE) Swiftscan solution combines a new Low Energy High Resolution and Sensitivity collimator (LEHRS) with image processing (Clarity 2D) and tomographic step and shoot continuous mode. The aim of this study was to compare clinical and physical performances of this new technology in bone scintigraphy.

METHODS

Physical phantom measurements were performed using GE LEHRS, GE Low Energy High Resolution (LEHR) and Siemens LEHR collimators. These measurements were associated with a prospective clinical study. Sixty-seven patients referred for bone scintigraphy were enrolled from February to July 2018. Each patient underwent two acquisitions consecutively on GE and Siemens gamma camera, using respectively Swiftscan solution and LEHR collimator.

RESULTS

On planar acquisitions, maximum sensitivity was 100 cts/MBq for Siemens LEHR. GE SwiftScan LEHRS and GE LEHR maximum sensitivity were respectively 9% and 22% lower. Using Clarity 2D, GE Swiftscan LEHRS spatial resolution was the best with 9.2 mm versus 10.1 mm and 10.6 mm for GE LEHR and Siemens LEHR collimators. In tomographic mode, the sensitivity of GE Swiftscan solution was superior to both LEHR systems (16% and 25% respectively for Siemens and GE). There was no significant difference in spatial resolution. In clinical use, signal was higher on Siemens system and noise was lower on GE Swiftscan solution. Contrast-to-noise ratios were not significantly different between the two systems. There was a significant image quality improvement with GE SwiftScan in planar images and in whole body scan. No significant difference in image quality was observed on SPECT images.

CONCLUSION

New GE SwiftScan collimator design improved sensitivity compared to "classical" GE LEHR collimator without compromising resolution. GE SwiftScan solution enhances planar image quality with a better Clarity 2D resolution recovery and noise treatment. In SPECT mode, GE SwiftScan solution improves volumetric sensitivity without significant impact on image quality, and could lead to time or dose reduction.

Quantitation of specific binding ratio in 123I-FP-CIT SPECT: accurate processing strategy for cerebral ventricular enlargement with use of 3D-striatal digital brain phantom

This study aimed to evaluate the effect of ventricular enlargement on the specific binding ratio (SBR) and to validate the cerebrospinal fluid (CSF)-Mask algorithm for quantitative SBR assessment of ¹²³I-FP-CIT single-photon emission computed tomography (SPECT) images with the use of a 3D-striatum digital brain (SDB) phantom. Ventricular enlargement was simulated by three-dimensional extensions in a 3D-SDB phantom comprising segments representing the striatum, ventricle, brain parenchyma, and skull bone. The Evans Index (EI) was measured in 3D-SDB phantom images of an enlarged ventricle. Projection data sets were generated from the 3D-SDB phantoms with blurring, scatter, and attenuation. Images were reconstructed using the ordered subset expectation maximization (OSEM) algorithm and corrected for attenuation, scatter, and resolution recovery. We bundled DaTView (Southampton method) with the CSF-Mask processing software for SBR. We assessed SBR with the use of various coefficients (f factor) of the CSF-Mask. Specific binding ratios of 1, 2, 3, 4, and 5 corresponded to SDB phantom simulations with true values. Measured SBRs > 50% that were underestimated with EI increased compared with the true SBR and this trend was outstanding at low SBR. The CSF-Mask improved 20% underestimates and brought the measured SBR closer to the true values at an f factor of 1.0 despite an increase in EI. We connected the linear regression function (y = - 3.53x + 1.95; r = 0.95) with the EI and f factor using root-mean-square error. Processing with CSF-Mask generates accurate quantitative SBR from dopamine transporter SPECT images of patients with ventricular enlargement.

Validation of a calibration method using the cross-calibration factor and system planar sensitivity in quantitative single-photon emission computed tomography imaging

The present study aimed to validate the absolute quantitative accuracy of a calibration method for single-photon emission computed tomography (SPECT) using cross-calibration factor (CCF)- and system sensitivity-based calibration methods. The CCF obtained with different reconstruction parameters was evaluated using a cylindrical phantom (diameter 20 cm, height 20 cm). SPECT images were acquired with a positron emission tomography/computed tomography (CT) phantom. Subsequently, they were reconstructed by using ordered subset expectation maximization with resolution recovery, scatter, and CT-based attenuation correction. All reconstructed SPECT counts were converted to activity concentrations based on the CCF and system planar sensitivity. We placed 12 circular regions of interest, 37 mm in diameter, on the phantom background, and the converted activity concentration and relative measurement error were assessed. The CCF obtained using a cylindrical phantom was affected by the iterative update number and post-smoothing filter function. The activity concentration calibrated using the CCF showed over- and underestimation. However, the activity concentration obtained from the system planar sensitivity was similar to that gained using the phantom. The values obtained using the system planar sensitivity were within 10% of the activity concentrations obtained with the phantom. These findings demonstrated that the calibration method using system planar sensitivity provides accurate quantification within 10% of the true activity concentration. Further clinical examination is required to validate the present results.

Effect of resolution recovery using graph plots on regional cerebral blood flow in healthy volunteers

PURPOSE

We evaluated the effect of resolution recovery (RR) using graph plots on regional cerebral blood flow (rCBF) in brain perfusion single-photon emission computed tomography (SPECT) images derived from healthy volunteers and patients diagnosed with probable Alzheimer's disease.

METHOD

We acquired brain perfusion SPECT images with scatter correction (SC), computed tomography-based attenuation correction (CTAC), and RR from a three-dimensional brain phantom and from healthy volunteers. We then compared contrast-to-noise ratio, count density ratios, increase maps, and rCBF using statistical parametric mapping 8.

RESULTS

Regional brain counts were significantly increased from 20-24% with SC, CTAC, and RR compared with SC and CTAC. Mean CBF in healthy volunteers was 42.5 ± 5.4 mL/100 g/min. Average rCBF determined using SC, CTAC and RR increased 7.5, 2.0, and 3.7% at the thalamus, posterior cingulate, and whole brain, respectively, compared with SC and CTAC.

CONCLUSION

Resolution recovery caused variations in normal rCBF because counts increased in cerebral regions.

Does applying resolution recovery to normal databases confer an advantage over conventional 3D-stereotactic surface projection techniques?

We evaluated a novel normal database (NDB) generated using single photon emission computed tomography (SPECT) data obtained from healthy brains by using a SPECT/CT system, analyzed using a resolution recovery (RR) technique applied to the three-dimensional stereotactic surface projection (3D-SSP) technique. We used a three-dimensional ordered subset expectation maximization method (3D-OSEM) with applied scatter correction (SC), attenuation correction, and RR to reconstruct the data. We verified the accuracy of the novel NDB's values (Z, extent, and error scores), and compared the novel NDB to the 3D-SSP technique by using simulated misery perfusion-related patient data from a conventional NDB. In addition, Z, extent, and error scores at the precuneus, cuneus, and posterior cingulate were compared under different reconstruction conditions by using the patient data. In the simulation, Z scores decreased when using the novel NDB corrected using computed tomography-based attenuation correction (CTAC), SC, and RR. The extent scores of the posterior cingulate increased using the novel NDB, relative to the other NDBs. The error score with the novel NDB without RR decreased by 15% compared to that of the conventional NDB. Z scores generated from patient data decreased in the novel NDB with RR. The extent scores tended to decrease in the novel NDB with RR. The extent scores in the novel NDB with RR improved at the posterior cingulate, compared to the scores with the other NDBs. However, applying RR to the novel NDB conferred no advantage because the cut-off of the current Z score must be reconsidered when using the additive RR technique.

Evaluation of Resolution Recovery for Each Collimator in Brain Perfusion Image

PURPOSE

This study aimed to verify the resolution recovery for each collimator in the brain perfusion image.

METHOD

To verify the effect of the resolution recovery for each collimator, we evaluated via the three-dimensional brain phantom (phantom) and the normal brain perfusion single photon emission computed tomography (SPECT) data. These data were reconstructed using the three-dimensional ordered subset expectation maximization method (3D-OSEM) (Evolution for bone^TM) that was performed with scatter correction, attenuation correction, and resolution recovery (RR). The performance of resolution recovery was evaluated in the two collimator systems (ELEGP and MEGP) reconstruction condition via the contrast value, mean counts, normalized mean square error (NMSE), and regional brain activity.

RESULT

In the "with resolution recovery (+RR)", the NMSE indicated minimum value with SI (subset×iteration) = 100, cut-off frequency (Fc) = 0.50 cycles/cm. The contrast value in the "+RR" increased 20% for the cortical region and decreased 28% and 6% at ELEGP collimator and MEGP collimator for the central region, as compared to the "without resolution recovery (-RR)". In the phantom study, the error of the brain activity using MEGP collimator at the temporal lobe and sub-lobar decreased 15%, compared with ELEGP collimator in the + RR. In the clinical study, the error of the regional brain activity using MEGP collimator in the "+RR" increased from 3% to 8%, compared with "-RR".

DISCUSSION

The accurate resolution recovery was obtained at SI = 100 and Fc = 0.50 cycles/cm. The contrast value and regional brain activity at the central region decreased due to incomplete resolution recovery by use of ELEGP collimator.

[Accuracy of Resolution Recovery in PSF-based Fully-3D PET Image Reconstruction: Simulation and Phantom Study in Multicenter Trial]

PURPOSE

Recently, the quality of positron emission tomography (PET) images has rapidly improved using resolution recovery algorithm with point spread function (PSF). The aim of this study was to investigate the accuracy of the resolution recovery algorithm using three different PET systems.

METHODS

Three PET scanner models, the GE Discovery 600 M (D600M), SIEMENS Biograph mCT (mCT), and SHIMADZU SET-3000GCT/X (3000GCT) were used in this study. The radial dependences of spatial resolution (full width at half maximum: FWHM) were obtained by point source measurements (0.9 mmφ). All PET data were acquired in three-dimensional (3D) mode and reconstructed using the filtered back projection (FBP) , 3D-ordered subsets expectation maximization (3D-OSEM or dynamic row-action maximum likelihood algorithm) , and 3D-OSEM+PSF (PSF) algorithms. Two indicators, aspect ratio (ASR) and resolution recovery ratio (RRR), were calculated from measured FWHMs and compared among the three PET scanners.

RESULTS

In D600 and 3000GCT, distortions of the radial direction were slightly increased at circumference of field of view (FOV). On the other hand, random distortions were occurred in both radial and tangential direction in mCT. ASRs calculated from 3D-OSEM images at circumference of FOV were 2.06, 1.22, and 2.04 on D600M, mCT, and 3000GCT, respectively. ASR improved with PSF in all PET scanners. On the other hand, RRR with PSF were calculated 57.6%, 61.4%, and 31.6%, respectively.

CONCLUSION

Our results suggest that the spatial resolutions of PET images could be improved with PSF algorithm in all PET systems; however, effect of PSF was different depending on PET systems. Furthermore, PSF algorithm could not completely improve spatial resolutions in circumference of FOV.

Influence of reconstruction algorithms on image quality in SPECT myocardial perfusion imaging

INTRODUCTION

We investigated if image- and diagnostic quality in SPECT MPI could be maintained despite a reduced acquisition time adding Depth Dependent Resolution Recovery (DDRR) for image reconstruction. Images were compared with filtered back projection (FBP) and iterative reconstruction using Ordered Subsets Expectation Maximization with (IRAC) and without (IRNC) attenuation correction (AC).

MATERIALS AND METHODS

Stress- and rest imaging for 15 min was performed on 21 subjects with a dual head gamma camera (Infinia Hawkeye; GE Healthcare), ECG-gating with 8 frames/cardiac cycle and a low-dose CT-scan. A 9 min acquisition was generated using five instead of eight gated frames and was reconstructed with DDRR, with (IRACRR) and without AC (IRNCRR) as well as with FBP. Three experienced nuclear medicine specialists visually assessed anonymized images according to eight criteria on a four point scale, three related to image quality and five to diagnostic confidence. Statistical analysis was performed using Visual Grading Regression (VGR).

RESULTS

Observer confidence in statements on image quality was highest for the images that were reconstructed using DDRR (P<0·01 compared to FBP). Iterative reconstruction without DDRR was not superior to FBP. Interobserver variability was significant for statements on image quality (P<0·05) but lower in the diagnostic statements on ischemia and scar. The confidence in assessing ischemia and scar was not different between the reconstruction techniques (P = n.s.).

CONCLUSION

SPECT MPI collected in 9 min, reconstructed with DDRR and AC, produced better image quality than the standard procedure. The observers expressed the highest diagnostic confidence in the DDRR reconstruction.

[Evaluation of Resolution Correction in Single Photon Emission Computed Tomography Reconstruction Method Using a Body Phantom: Study of Three Different Models]

PURPOSE

The aim of this study was to evaluate the resolution recovery techniques of Flash3D, Astonish, and Evolution in single photon emission computed tomography(SPECT) using a body phantom.

METHODS

We scanned a National Electrical Manufactures Association body phantom filled with 99mTc. The body of the phantom with radioactive sphere and background was filled with either water or radioactive solution. We investigated image quality using profile curves, recovery coefficient, and image contrast.

RESULTS

The profile curve at the edge of the hot sphere showed artifact due to Gibbs oscillation for all techniques, and also over estimation of recovery coefficient was seen in the hot sphere, as had been previously reported in a simulation study. These phenomena were more remarkable than Evolution in the Flash3D and Astonish techniques. For the contrast between hot sphere and background, the contrast recover was enough for the <17-mm hot spheres. These results showed that the effect of contrast correction was less as the radius of rotation diameter became large.

CONCLUSION

In the present study using the body phantom, overestimated counts and edge artifacts due to Gibbs oscillation were shown. These phenomena were different by each resolution correction algorithms. Also, there were limitation regarding image quality improvement by resolution correction depending on sphere size and length of radius of rotation.

[The Optimal Reconstruction Parameters by Scatter and Attenuation Corrections Using Multi-focus Collimator System in Thallium-201 Myocardial Perfusion SPECT Study]

The aim of this study was to reveal the optimal reconstruction parameters of ordered subset conjugates gradient minimizer (OSCGM) by no correction (NC), attenuation correction (AC), and AC+scatter correction (ACSC) using IQ-single photon emission computed tomography (SPECT) system in thallium-201 myocardial perfusion SPECT. Myocardial phantom acquired two patterns, with or without defect. Myocardial images were performed 5-point scale visual score and quantitative evaluations using contrast, uptake, and uniformity about the subset and update (subset×iteration) of OSCGM and the full width at half maximum (FWHM) of Gaussian filter by three corrections. We decided on optimal reconstruction parameters of OSCGM by three corrections. The number of subsets to create suitable images were 3 or 5 for NC and AC, 2 or 3 for ACSC. The updates to create suitable images were 30 or 40 for NC, 40 or 60 for AC, and 30 for ACSC. Furthermore, the FWHM of Gaussian filters were 9.6 mm or 12 mm for NC and ACSC, 7.2 mm or 9.6 mm for AC. In conclusion, the following optimal reconstruction parameters of OSCGM were decided; NC: subset 5, iteration 8 and FWHM 9.6 mm, AC: subset 5, iteration 8 and FWHM 7.2 mm, ACSC: subset 3, iteration 10 and FWHM 9.6 mm.

[Validation of noise reduction in iterative reconstruction: a simulation phantom study]

PURPOSE

In the iterative reconstruction method, image noise tends to increase in proportion to falling available photon count and increasing update number. Image filtering is an important factor in single photon emission computed tomography (SPECT) image reconstruction, but it is frequently treated in a subjective way. The aim of this study was to evaluate the effects of pre-reconstruction filtering and post-reconstruction filtering on the iterative reconstruction process.

METHODS

Using simulation phantoms, projection data were reconstructed using ordered subsets expectation maximization (OSEM) with or without compensation for resolution recovery. Pre- and post-reconstruction filtering was performed using a Butterworth filter (BW) (range: 0.3-1.3 cycles/cm) and a Gaussian filter (GA) (range: 0.3-1.3 mm) with various parameters. We evaluated the variances of full width at half maximum (FWHM), coefficients of variation (CV), image contrast and normalized mean squared error (NMSE) values.

RESULTS

The FWHM values for pre-reconstruction filtering tended to be lower than those observed for post-reconstruction filtering. These values were 5.1 mm (pre-reconstruction) and 6.7 mm (post-reconstruction). The CV on pre- and post-reconstruction filtering was 7.5% and 11.6%. Pre-reconstruction filtering reduced image noise more effectively than post-reconstruction filtering. The contrast for pre-reconstruction filtering was similar to that observed after post-reconstruction filtering. However, contrast after filtering with a GA slowly decreased as compared to the BW. NMSE values obtained by pre-reconstruction filtering tended to be lower than those observed for post-reconstruction filtering.

CONCLUSIONS

Pre-reconstruction filtering provided SPECT image quality comparable to that from post-reconstruction filtering, especially when using the BW. Our results suggest that pre-reconstruction filtering is a beneficial method when applied to the iterative reconstruction method with or without compensation for resolution recovery.

Optimization of iterative reconstruction parameters with attenuation correction, scatter correction and resolution recovery in myocardial perfusion SPECT/CT

OBJECTIVE

The aim of this study was to characterize the optimal reconstruction parameters for ordered-subset expectation maximization (OSEM) with attenuation correction, scatter correction, and depth-dependent resolution recovery (OSEMACSCRR). We assessed the optimal parameters for OSEMACSCRR in an anthropomorphic torso phantom study, and evaluated the validity of the reconstruction parameters in the groups of normal volunteers and patients with abnormal perfusion.

METHODS

Images of the anthropomorphic torso phantom, 9 normal volunteers and 7 patients undergoing myocardial perfusion SPECT were acquired with a SPECT/CT scanner. SPECT data comprised a 64×64 matrix with an acquisition pixel size of 6.6 mm. A normalized mean square error (NMSE) of the phantom image was calculated to determine both optimal OSEM update and a full width at half maximum (FWHM) of Gaussian filter. We validated the myocardial count, contrast and noise characteristic for clinical subjects derived from OSEMACSCRR processing. OSEM with depth-dependent resolution recovery (OSEMRR) and filtered back projection (FBP) were simultaneously performed to compare OSEMACSCRR.

RESULTS

The combination of OSEMACSCRR with 90-120 OSEM updates and Gaussian filter with 13.2-14.85 mm FWHM yielded low NMSE value in the phantom study. When we used OSEMACSCRR with 120 updates and Gaussian filter with 13.2 mm FWHM in the normal volunteers, myocardial contrast showed significantly higher value than that derived from 120 updates and 14.85 mm FWHM. OSEMACSCRR with the combination of 90-120 OSEM updates and 14.85 mm FWHM produced lowest % root mean square (RMS) noise. Regarding the defect contrast of patients with abnormal perfusion, OSEMACSCRR with the combination of 90-120 OSEM updates and 13.2 mm FWHM produced significantly higher value than that derived from 90-120 OSEM updates and 14.85 mm FWHM. OSEMACSCRR was superior to FBP for the % RMS noise (8.52±1.08 vs. 9.55±1.71, p=0.02) and defect contrast (0.368±0.061 vs. 0.327±0.052, p=0.01), respectively.

CONCLUSIONS

Clinically optimized the number of OSEM updates and FWHM of Gaussian filter were (1) 120 updates and 13.2 mm, and (2) 90-120 updates and 14.85 mm on the OSEMACSCRR processing, respectively. Further assessment may be required to determine the optimal iterative reconstruction parameters in a larger patient population.