Optimization of scatter correction method in samarium-153 single-photon emission computed tomography using triple-energy window: A monte carlo simulation study :Hicham Asmi, Farida Bentayeb, Youssef Bouzekraoui, Faustino Bonutti, Indian Journal of Nuclear Medicine

ORIGINAL ARTICLE
Year : 2020 | Volume : 35 | Issue : 3 | Page : 210--215

Optimization of scatter correction method in samarium-153 single-photon emission computed tomography using triple-energy window: A monte carlo simulation study

Hicham Asmi¹, Farida Bentayeb¹, Youssef Bouzekraoui¹, Faustino Bonutti²,
¹ Department of Physics, LPHE, Modeling and Simulations, Faculty of Science, Mohammed V University, Rabat, Morocco
² Department of Medical Physics, Academic Hospital of Udine, Udine, Italy

Correspondence Address:
Dr. Youssef Bouzekraoui
Department of Physics, LPHE, Modeling and Simulations, Faculty of Science, Mohammed V University, Rabat
Morocco

Abstract

Purpose: In single-photon emission computed tomography imaging, the presence of scatter degrades image quality. The goal of this study is to optimize the main- and sub-energy windows for triple-energy window (TEW) method using Monte Carlo SImulating Medical Imaging Nuclear Detectors (SIMIND) code for samarium-153 (Sm-153) imaging. Materials and Methods: The comparison is based on the Monte Carlo simulation data with the results estimated using TEW method. Siemens Symbia gamma-camera equipped with low-energy high-resolution collimator was simulated for Sm-153 point source located in seven positions in water cylindrical phantom. Three different main-energy window widths (10%, 15%, and 20%) and three different sub-energy window widths (2, 4, and 6 keV) were evaluated. We compared the true scatter fraction determined by SIMIND and scatter fraction estimated using TEW scatter correction method at each position. In order to evaluate the image quality, we used the full width at half maximum (FWHM) computed on the PSF and image contrast using Jaszczak phantom. Results: The scatter fraction using TEW method is similar to the true scatter fraction for 20% of the main-energy window and 6 keV sub-energy windows. For these windows, the results show that the resolution and contrast were improved. Conclusion: TEW method could be a useful scatter correction method to remove the scatter event in the image for Sm-153 imaging.

How to cite this article:
Asmi H, Bentayeb F, Bouzekraoui Y, Bonutti F. Optimization of scatter correction method in samarium-153 single-photon emission computed tomography using triple-energy window: A monte carlo simulation study.Indian J Nucl Med 2020;35:210-215

How to cite this URL:
Asmi H, Bentayeb F, Bouzekraoui Y, Bonutti F. Optimization of scatter correction method in samarium-153 single-photon emission computed tomography using triple-energy window: A monte carlo simulation study. Indian J Nucl Med [serial online] 2020 [cited 2023 Mar 31 ];35:210-215
Available from: https://www.ijnm.in/text.asp?2020/35/3/210/288449

Full Text

Introduction

In nuclear medicine, samarium-153 (Sm-153) isotope is a potential choice because it emits both medium-energy β - particles (Eβ max = 0.80 MeV) with a short half-life (46.7 h) and gamma-photons suitable for imaging at 103 keV. The physical characteristics of Sm-153 allow to be considered as an excellent radiotherapeutic and diagnostic image agent.[1],[2]

Previous works [2],[3],[4] demonstrated that the isotope chelated to ethylenediamine tetramethylene phosphonate is an effective treatment of bone metastases by the excellent biolocalization. Moreover, Sm-153 is potentially suitable as an alternative to 90Y in liver cancer treatment with advantage of gamma-radiation for imaging.[5] Images obtained additionally permit be a rapid diagnostic of the therapeutic isotope. In Sm-153 single-photon emission computed tomography (SPECT) imaging with a gamma-camera, the presence of scatter introduces signifcant uncertainty in quantifcation of activity distribution. The scattered and primary photons cannot be determined experimentally. While, with the help of Monte Carlo simulation, it is possible to track the photons originating from the source that ultimately deposits its complete energy inside the crystal. It is increasingly used in nuclear medicine to develop new imaging parameters, scatter correction methods, and reconstruction algorithms. We have used Monte Carlo simulating medical imaging nuclear detectors (SIMIND) to accurately assess the contribution of scattered photons in the photopeak window. The quantification of gamma-camera imaging is improved after the correction of scattered radiation. Previous works [6],[7],[8],[9],[10],[11],[12],[13],[14],[15],[16],[17] were used triple-energy window (TEW) scatter correction to eliminate the detected scattered counts inside energy window. In order to quantify emission from the isotope Sm-153 using a gamma-camera accurately, however, it is important to correct the scattered photons which degrade image quality. We can reduce the counts of scattered photons in a photopeak energy window using TEW method, which is a simple method to use in clinical study for SPECT imaging. However, there are to date no study of measurements of the scatter fraction of Sm-153 as a function of the energy windows, in order to determine the optimal main- and sub-energy window for Sm-153 SPECT imaging. In this study, we assessed the fraction of scattered photons and determined the optimal main- and sub-energy windows for TEW scatter correction method for Sm-153 by means of Monte Carlo simulation.

Materials and Methods

Detection system description

In this study, we simulated Siemens Medical System Symbia equipped with low energy high resolution [Table 1]. The images were acquired by a single-head SPECT system (Symbia) based on 103 keV peak. The dimension of detector surface was 59.1 cm × 44.5 cm and having 2.54 cm NaI (Tl) crystal thickness. A water-filled cylinder phantom (diameter: 22 cm, length: 32 cm) was placed at 12 cm from the detector surface. We used the SIMIND Monte Carlo program to acquired data from Sm-153 point sources of 0.05-cm diameter located in different seven positions at the center of the cylinder phantom and offset by ± 5 cm in the X, Y, and Z directions relative to the center. Moreover, the Jaszczak phantom consists of six spheres with different diameters (31.8, 25.4, 15.9, 19.1, 12.7, and 9.5 mm) which are used to evaluate the image contrast. [Figure 1] illustrates the SIMIND simulation for each point source.{Table 1}{Figure 1}

We collected 50 million photons during the simulation. The images had a 0.34 cm pixel size and 128 × 128 matrix size. [Figure 2] shows photon energies and intensities of Sm-153 radionuclide decay.{Figure 2}

Monte Carlo simulation

The Monte Carlo simulation SIMIND code describes a SPECT camera.[17] The SIMIND program has two main programs: CHANGE, which defines the system parameters, and SIMIND, which executes the simulation. Moreover, using CHANGE program, we can introduce the desired parameters of system. SIMIND accurately simulates all interactions of photons inside the collimator.[18] At the end of simulation, SIMIND provides the value of geometric, penetration, scatter, and X-ray component and image in separate files. We imported binary images created by SIMIND in ImageJ software (Bethesda, Maryland, USA).[19] The validation of the SIMIND program for SPECT imaging has been verified for different gamma-camera according to the previous studies.[20],[21],[22],[23],[24]

Triple-energy window method

The TEW method estimates the counts of the scattered photons in the main photopeak window from the counts acquired in two sub-windows on both sides of this window. At each pixel in planer image, the counts of scattered photons are subtracted from the total counts in photopeak window to obtain the count of primary photons. If the measured count is not enough, we can enlarge the width of the sub-windows.

In this study, we used the main-energy window widths (10%, 15%, and 20%) centered on 103 keV and sub-energy window widths (2, 4, and 6 keV) [Table 2].{Table 2}

The true scattered photon fractions for the main-energy window calculated by SIMIND were compared with the scattered photon fraction assumed from TEW scatter correction method. The counts of scattered photons and counts of primary photons were calculated by relationships below [14.15]:

[INLINE:1]

[INLINE:2]

where Cleft: counts in lower sub-energy window; Cright: counts in upper sub-energy window; Ws: width of sub-energy window; Wm: width of main window; Ctot: counts in main window; Csca: scatter counts; and Cp: primary counts. The scatter-to-total ratio (scatter fraction) was calculated as Eq. (3):[6]

[INLINE:3]

The quality of the Sm-153 SPECT image was quantitatively evaluated using the image contrast assessment. Contrast was calculated by the following formula using Eq. (4):

[INLINE:4]

where Ms and Mb are the mean pixel values of the activity of spheres and the activity of background as noise, respectively.

Results

The simulated energy spectrum is shown in [Figure 3]. The counts acquired in different windows are shown in [Table 3].{Figure 3}{Table 3}

As shown in [Figure 4], It is clear that the geometric component is large and remains constant with increase in photopeak window. [Table 4] shows the comparison of true scatter fraction (%) and scatter fraction estimated by TEW scatter correction method at each position. The scatter fraction depended on the source position in cylindrical phantom. The results showed that, for 20% of the main-energy window and 6 keV sub-energy window of Sm-153, the scatter fraction estimated by TEW is similar to the true scatter fraction determined by SIMIND.{Figure 4}{Table 4}

[Figure 5] demonstrates the images for Sm-153 point source at seven positions. The calculated vertical and horizontal full width at half maximum (FWHM) on the images is shown in [Table 5]. It shows that, for each position, the TEW method is decreasing the FWHM.{Figure 5}{Table 5}

[Figure 6] shows the reconstructed images of simulated Jaszczak phantom SPECT before and after scatter correction. It noted that the contrast of the six spheres was improved with 20% of the main-energy window and 6 keV sub-energy windows for TEW, as shown in [Figure 7]. Therefore, 6 keV sub-window with a 20% main-energy was possible energy windows setting for the TEW method in Sm-153.{Figure 6}{Figure 7}

Discussion

A previous study [12] has assessed the TEW method for Sm-153 SPECT data. But until now, no study has compared the fraction of scattered photons of Sm-153 with varying main- and sub-energy windows for the implementation of the TEW method. The scattered photons are a major factor degrade resolution image and image contrast.[25],[26],[27] The resolution is mostly expressed as the FWHM of the Point Spread Function (PSF). A smaller FWHM implies the better quality of the image. The FWHM and contrast were utilized to assess the effect of the processing method on image quality. The choice of optimal energy window in TEW method has a key role which appears the lowest scatter fraction. We have determined the optimal windows that it displays the similar scatter fraction calculated by means of simulation and TEW method. Different scatter fractions were obtained for each source location. There was a smaller difference with 20% main-energy window and 6 keV secondary energy window. By looking at [Table 4], it is obvious that the calculated spatial resolution was improved with TEW method. As shown in [Figure 5], the 20% main-energy window with 6 keV sub-energy windows produces a better contrast when using the TEW correction. The simplicity of this method makes it feasible in a clinical study. However, the estimation of scattered photons cannot be determined experimentally; with the help of Monte Carlo simulation, accurate assessment of the scattered photon fractions inside photopeak window can be made.

Conclusion

In this study, we used the Monte Carlo SIMIND code with TEW scatter correction method to correct the scatter events detected in photopeak window for Sm-153. The results indicate that it is better to use a 20% main photopeak window with 6 keV sub-windows when TEW method is applied.

Financial support and sponsorship

Nil.

Conflicts of interest

There are no conflicts of interest.

References

1	Yeong CH, Abdullah BJ, Ng KH, Chung LY, Goh KL, Sarji SA, et al. Production and first use of 153SmCl3-ion exchange resin capsule formulation for assessing gastrointestinal motility. Appl Radiat Isot 2012;70:450-5.
2	Iagaru AH, Mittra E, Colletti PM, Jadvar H. Bone-targeted imaging and radionuclide therapy in prostate cancer. J Nucl Med 2016;57:19S-24S.
3	Naseri Z, Jalilian AR, Kharat AN, Bahrami-Samani A, Ghannadi-Maragheh M. Production, quality control and biological evaluation of 153Sm-TTHMP as a possible bone palliation agent. Iran J Nucl Med. 2011;19:60-8.
4	Kathiresan RK, Begum B, Rangarajan. Comparison of Tc-99m MDP and Sm-153 EDTMP bone scan. Indian J Nuc Med 2011;26:163-4.
5	Hashikin NA, Yeong CH, Abdullah BJ, Ng KH, Chung LY, Dahalan R et al. Neutron activated samarium-153 microparticles for transarterial radioembolization of liver tumour with post procedure imaging capabilities. PLoS One 2015;10:e0138106.
6	Lee YS, Kim JS, Kim KM, Lim SM, Kim HJ. Determination of energy windows for the triple energy window scatter correction method in I-131 on a Siemens SYMBIA gamma camera: A GATE simulation study. J Inst 2015;15:10:1-8.
7	Takayama T, Ichihara T, Motomura N, Ogawa K. Determination of energy window width and position for scintigraphic imaging using different energy resolution detection with the 'triple energy window (TEW) Scatter Compensation Method IEEE. Kaku Igaku 1998;35:51-9.
8	Changizi V, Takavar A, Babakhani A, Sohrabi M. Scatter correction for heart SPECT images using TEW method. J Appl Clin Med Phys 2008;9:136-40.
9	Dewaraja YK, Li J, Koral K. Quantitative 131I SPECT with triple energy window Compton scatter correction. IEEE Trans Nucl Sci 1998;45:3109-14.
10	Fujioka H, Inoue T, Ishimaru Y, Akamune A, Murase K, Tanada S, et al. Compton scatter correction using the triple energy window (TEW) method in conventional single photon emission computed tomography without TEW acquisition hardware. Ka Iga1997;34:251-8.
11	Ogawa K. Simulation study of triple-energy-window scatter correction in combined TI-201. Tc-99m SPECT. Ann Nucl Med 1994;8:277-81.
12	Asgari A, Ashoor M, Sohrabpour M, Shokrani P and Rezaei A. Evaluation of various energy windows at different radionuclides for scatter and attenuation correction in nuclear medicine. Ann Nucl Med 2015;29:375-83.
13	Kim KM, Varrone A, Watabe H, Shidahara M, Fujita M, Innis RB, et al. Contribution of scatter and attenuation compensation to SPECT images of nonuniformly distributed brain activities. J Nucl Med 2003;44:512-9.
14	Ogawa K, Harata Y, Ichihara T, Kubo A, Hashimoto S. A practical method for position-dependent Compton-scatter correction in single photon emission CT. IEEE Trans Med Imaging 1991;10:408-12.
15	Bong JK, Son HK, Lee JD, Kim HJ. Improved scatter correction for SPECT images: A Monte Carlo study. IEEE Trans Nucl Sci 2005;52:1263-70.
16	Bouzekraoui Y, Bentayeb F, Asmi H, Bonutti F. Determination of the energy windows for the triple energy window scatter correction method in gadolinium-159 single photon emission computed tomography using Monte Carlo simulation. Iranian J Med Phys 2019.
17	Ljungberg M. The SIMIND Monte Carlo Program Home Page. Available from: https://www.msf.lu.se/research/the-simind-monte-carlo-program. [Last accessed on 2019 Nov 05].
18	Pandey AK, Sharma SK, Karunanithi S, Kumar P, Bal C, Kumar R. Characterization of parallel–hole collimator using Monte Carlo simulation Indian J Nucl Med 2015;30:128-34.
19	ImageJ Program. Available from: https://imagej.nih.gov/ij/download.html. [Last accessed on 2019 Nov 05].
20	Ejeh JE, van Staden JA, du Raan H. Validation of SIMIND Monte Carlo simulation software for modelling a siemens symbia T SPECT scintillation camera. IFMBE Proceedings 2018;573-6.
21	Toossi BM, Islamian PJ, Momennezhad M, Ljungberg M, Naseri SH. SIMIND Monte Carlo simulation of a single photon emission CT. J Med Phys 2010;35:42-7.
22	Azarm A, Islamian JP, Mahmoudian B, Gharepapagh E. The effect of parallel-hole collimator material on image and functional parameters in SPECT imaging: A SIMIND Monte Carlo study. World J Nucl Med 2015;14:160-4.
23	Vicente EM, Lodge MA, Rowe SP, Wahl RL, Frey EC. Simplifying volumes-of-interest (VOIs) definition in quantitative SPECT: Beyond manual definition of 3D whole-organ VOIs. Med Phys 2017;44:1707-17.
24	Rafati M, Rouhani H, Bitarafan-Rajabi A, Noori-Asl M, et al. Assessment of the scatter correction procedures in single photon emission computed tomography imaging using simulation and clinical study. J Cancer Res Ther 2017;13:936-942.
25	ALehyani SH. Application of single photon emission computed tomography (SPECT) parameters for bone scintigraphy. J King Saud Univ – Sci 2009;21:109-17.
26	Asl MN, Sadremomtaz A, Bitarafan-Rajabi A. Evaluation of six scatter correction methods based on spectral analysis in (99m) Tc SPECT imaging using SIMIND Monte Carlo simulation J Med Phys 2013;38:189-97.
27	Dewaraja YK. Ljungberg M and Koral KF. Characterization of scatter and penetration using Monte Carlo simulation in 131I imaging. J Nucl Med 2000;41:123-30.