|Year : 2010 | Volume
| Issue : 1 | Page : 42-47
SIMIND Monte Carlo simulation of a single photon emission CT
MT Bahreyni Toossi1, J Pirayesh Islamian2, M Momennezhad3, M Ljungberg4, SH Naseri2
1 Medical Physics Research Center, Mashhad University of Medical Sciences, Iran
2 Medical Physics Department, Faculty of Medicine, Mashhad University of Medical Sciences, Mashhad, Iran
3 Nuclear Medicine Department, Faculty of Medicine, Mashhad University of Medical Sciences, Mashhad, Iran
4 Medical Radiation Physics Department, Clinical Sciences-Lund, Lund University, Lund, Sweden
|Date of Submission||20-Mar-2009|
|Date of Decision||28-Apr-2009|
|Date of Acceptance||01-May-2009|
|Date of Web Publication||7-Jan-2010|
J Pirayesh Islamian
Medical Physics Department, Mashhad University of Medical Sciences, Post box 9177948564, Mashhad
Source of Support: None, Conflict of Interest: None
| Abstract|| |
In this study, we simulated a Siemens E.CAM SPECT system using SIMIND Monte Carlo program to acquire its experimental characterization in terms of energy resolution, sensitivity, spatial resolution and imaging of phantoms using 99m Tc. The experimental and simulation data for SPECT imaging was acquired from a point source and Jaszczak phantom . Verification of the simulation was done by comparing two sets of images and related data obtained from the actual and simulated systems. Image quality was assessed by comparing image contrast and resolution. Simulated and measured energy spectra (with or without a collimator) and spatial resolution from point sources in air were compared. The resulted energy spectra present similar peaks for the gamma energy of 99m Tc at 140 KeV. FWHM for the simulation calculated to14.01 KeV and 13.80 KeV for experimental data, corresponding to energy resolution of 10.01and 9.86% compared to defined 9.9% for both systems, respectively. Sensitivities of the real and virtual gamma cameras were calculated to 85.11 and 85.39 cps/MBq, respectively. The energy spectra of both simulated and real gamma cameras were matched. Images obtained from Jaszczak phantom, experimentally and by simulation, showed similarity in contrast and resolution. SIMIND Monte Carlo could successfully simulate the Siemens E.CAM gamma camera. The results validate the use of the simulated system for further investigation, including modification, planning, and developing a SPECT system to improve the quality of images.
Keywords: Emission tomography, Monte Carlo simulation, SIMIND, SPECT, 99m Tc imaging
|How to cite this article:|
Bahreyni Toossi M T, Islamian J P, Momennezhad M, Ljungberg M, Naseri S H. SIMIND Monte Carlo simulation of a single photon emission CT. J Med Phys 2010;35:42-7
|How to cite this URL:|
Bahreyni Toossi M T, Islamian J P, Momennezhad M, Ljungberg M, Naseri S H. SIMIND Monte Carlo simulation of a single photon emission CT. J Med Phys [serial online] 2010 [cited 2019 Nov 21];35:42-7. Available from: http://www.jmp.org.in/text.asp?2010/35/1/42/55967
| Introduction|| |
Recent developments in nuclear medicine instrumentation and multiple-processor parallel processing systems have created a need for Monte Carlo simulation opportunities in nuclear medicine imaging. Since nuclear medicine imaging deals with random phenomena such as radioactive decay, emission of radiation energy through photons and particles, and the detection of these quanta and particles in various material, Monte Carlo simulations are nowadays employed as an essential tool in nuclear medicine imaging, both in single-photon emission computed tomography (SPECT) and positron emission tomography (PET). ,, One of the aims of a medical physicist, involved in nuclear medical imaging research, is to optimize the design of imaging systems and improve qualitative and quantitative accuracy of reconstructed images. Several factors affect image quality and accuracy of the data obtained from a nuclear medicine scan. These include: physical properties of detectors, collimator and gantry design, attenuation and scatter compensation and reconstruction algorithms. , Integrated improvement in these parameters with current tracers and sensitive and specific tracers under development will provide major advantages to the general nuclear medicine clinician and research investigator. Mathematical modeling is necessary for assessment of various parameters in nuclear medical imaging systems since no analytical solution is possible when solving the transport equation describing the interaction of photons with nonuniformly attenuating body structures and complex detector geometries. One of the most frequent applications of Monte Carlo methods in nuclear medicine imaging is detector modeling with three main purposes: to study interactions within the radiation sensor for each photon and thus correct for sources of image degradation; evaluate techniques of image treatment to quantify the effects of dispersion and attenuation, and to optimize reconstruction algorithms; and patient dosimetry calculations. To improve quality, some authors have suggested improving the detection system itself, especially the collimator. ,, Monte Carlo simulation appears to be the best solution to model gamma-camera behavior, ,,,,,,, although the long calculation time implied has often limited its application. A detailed description of the general principles of the Monte Carlo method is given in a number of publications. ,,,
A number of Monte Carlo simulation codes (e.g. EGS4, MCNP4, SimSET, Geant4, Gate, and SIMIND) applicable to nuclear medicine have been developed and described in literature. ,,
The design of SPECT and PET systems using the Monte Carlo method has received considerable attention and such investigations resulted in several applications. , The study used the SIMIND Monte Carlo simulation program, which is well established for SPECT with low-energy photons, for photonic physics and other applications. This program was originally designed for the calibration of whole-body counters, but soon evolved to simulate scintillation cameras. ,, It is now available in Fortran-90 and can be freely downloaded from http://www.radfys.lu.se/simind and run on major computer platforms including PCs.  The SIMIND program actually consists of two programs, CHANGE, which defines the parameters, and SIMIND, which performs the actual simulation. The program can simulate nonuniform attenuation from voxel-based phantoms and includes several types of variance reduction techniques. Transmission SPECT images can also be simulated. For particular projects, the user can write scoring routines linked to the main code. The program has been modified to allow it to be run on parallel computers using the MPE command language. The use of SIMIND is well documented. 
In this article, we present the simulation of a Siemens E.CAM gamma camera,  used in clinical practice in our center.
| Materials and Methods|| |
The Siemens E.CAM gamma camera was modeled using SIMIND Monte Carlo program. The camera consists of a removable low energy high resolution (LEHR) collimator, a NaI (Tl) scintillation crystal, a light guide and an array of photomultiplier tubes (PMTs). The parameters of LEHR collimator, used for low energy sources such as 99m Tc, for experiment and simulation were as follows: parallel hexagonal holes with cells of 1.11 mm diameter, 2.405 cm height, and 0.16 mm septal thickness. The NaI (Tl) crystal specifications are as follows: planar, 9.5 mm in thickness, 59.1 x 44.5 cm 2 in area, light yield 40k photons /MeV, and a peak emission spectrum at 415 nm. , Generated light in the crystal is collected by a matrix composed of 59 PMTs, 53 with 7.6 cm and 6 with 5.1 cm in diameter. The photocathode is a bialkali type with quantum efficiency of approximately 30% for the wavelength of maximum NaI (Tl) emission.  A light guide ensures a good optical coupling between the scintillating crystal and PMTs. The SIMIND Monte Carlo program was utilized to simulate the aforementioned structures. When 99m Tc is used, various structures attached to back of the crystal contribute to backscattering of the emitted photons. To assess the effect of these parts, a single 6 cm slab of Pyrex  was substituted and simulated.
The scintillation process was simulated for generation of optical photons within the NaI (Tl) crystal. An isotropic source emitting 140 keV photons was assigned to 99m Tc. To obtain the output spectra of the simulated gamma camera a 99m Tc point source was simulated for the following geometries and activities:
The experimental output spectra of the real gamma camera were acquired for exactly the same conditions as applied to the simulated system. Acquired data representing the system spectra was stored in digital format. An excel program was employed to display the energy spectra.
- 3.7 MBq, without collimator, positioned at 0, 10, 15, 20, and 25cm from the detector surface.
- 37 MBq, with collimator, and at the same distances as in (a).
In this work we have studied the following properties of the gamma camera: ,,,,,, energy resolution, spatial resolution and image contrast. An energy window was centered on the 99m Tc photo peak (130 - 151 keV).  The images were reconstructed in matrices of 128Χ 128 pixels, with a pixel size of 0.39 mm.
Energy resolution of the gamma camera was measured with a 99m Tc (3.7MBq) point source, positioned at the center of the field of view (FOV), 25 cm from the crystal surface. The energy spectrum was acquired for 10 7 photons / projection. The simulated and experimental energy spectra were obtained with and without a low energy high resolution collimator for a source-detector distance of 25 cm.
Spatial resolution of the real and simulated gamma camera were determined by placing a 99m Tc (1.9MBq) source (1.5 mm in diameter) at the center of the FOV. A study of the SPECT reconstructed spatial resolution was also carried out both experimentally and by SIMIND simulation. SPECT projections of a Jaszczak Deluxe Phantom  along the axis of rotation were acquired. The phantom was uniformly filled in with 370MBq 99m Tc, and positioned 15 cm from the collimator surface. The projections were obtained in the 15% window (130 - 151 keV). Spatial resolution (in mm) was obtained by the smallest visible and recognizable rods.
Sensitivity of the gamma camera is determined by taking ratio of the detected counts per second in the selected energy window per unit activity in the source (cps/MBq). The system sensitivity was experimentally measured and estimated by simulation for the point source with the condition specified in section 2. The measured sensitivity was obtained for a 3.7 MBq source counted over a period of 60 seconds, whereas the corresponding simulation was performed for 10 million photons generated for the same source position. Decay time and background radiation were taken into consideration. Finally, the response of the E-CAM in extrinsic mode (i.e., with the collimator) was studied for a point source placed 25 cm from the detector.
The Jaszczak Deluxe phantom was also employed to evaluate the quality measurement of images obtained experimentally and by simulation. , The SPECT parameters were the same as already described in section 3, also same configuration was applied to SIMIND simulation of the phantom. The acquisition parameters were identical for experiments and simulations: 128x128 matrix, 128 views, 1.23 zoom factor, 3.9 mm pixel size, and with approximately 1 million counts / projection resulting in 128 million total counts.
| Results and Discussion|| |
[Figure 1] shows the energy spectra produced by the real and simulated gamma cameras for a 99m Tc point source, positioned 25 cm from the detector surface. Some minor differences may be observed between the simulated and the experimental energy spectra, the most striking being that the experimental spectrum presents a wider peak which may be explained by the superimposition of the energy peaks corresponding to the X-ray escape of 53 I present in the NaI (Tl) crystal (~110 keV), the 99m Tc photo peak (140 keV) and the sum of 140 keV with the x-ray energy of 99m Tc (total ~160 keV),  which cannot be separated by the detector. However, the results obtained for the spectrum peak at 140 keV is satisfactory.
The calculated FWHMs for 140 keV photo peak were equal to 14.01 keV and 13.80 keV for simulated and experimental normalized energy spectra, respectively. The corresponding relative energy resolution, with and without a LEHR collimator, obtained for simulated and experimentally acquired spectra are 9.94%, 9.61% and 10.01%, 9.86%, respectively. Contributions of Compton and photoelectric interactions in the whole spectra (simulated) and in the selected window are given in [Table 1]. There is a good agreement between the results obtained from simulations and experiments. However, there were some differences between simulated results obtained for the FWHMs that were constant over the range of source-collimator selected distances. The differences are due to the fact that the size of point source defined by the SIMIND is negligible compared to the 1.5 mm diameter of the source used in experimental measurements. Our measured and computed relative energy resolution (intrinsic) compare very well with the 9.9% value claimed by the manufacturer. 
Spatial resolution for a source to collimator distance of 10 cm was found to be 8.4 plus/minus 0.1 mm and 7.8 plus/minus 0.1mm for the actual and simulated systems, respectively. The simulated and experimental spatial resolutions were closely the value provided by manufacturer for a point source at this distance is 7.8 mm,  which is in good agreement with our experimental value.
To evaluate the spatial resolution of reconstructed SPECT of a Jaszczak phantom, the images were obtained by experiment and modeling. [Figure 2] shows corresponding images of the cold rods and spheres in the phantoms, prepared by simulation and by experiment. The best results were obtained for the simulated hot rods and spheres [Figure 3]. The qualities of the produced images were compared in terms of image contrast and spatial resolution.  From [Figure 2] and [Figure 3] it is evident that reconstructed spatial resolution of both SPECTs is nearly equal to 9.5mm. Image contrasts for the six spheres were calculated by equations (1) and (2) and presented in [Table 2].
Contrast H = N sp. / N b. Eq.1
Contrast C = 1- (M sp ./M cy. ) Eq.2
The N sp and N b are the mean pixel values of the hot spheres and background, respectively. M sp and M cy correspond to the minimum pixel value in cold spheres and the maximum pixel values in the phantom cylinder, respectively.
From [Table 2] it can be seen that contrast decreased with decreasing sphere diameter, as would be expected. The calculated contrasts for cold and hot simulated spheres with the diameter larger than 15.9 mm have a remarkable value compared to the diameter of 9.5 mm. Similar results were also seen for the experiment with an apparent similarity for sphere sizes above the diameter of 15.9 mm. However, there was no significant change in contrast for the 12.7- and 9.5-mm spheres. Therefore, the obtained contrast similarity for cold spheres, as an imaging parameter, is another evident for verification of the simulation.
Sensitivities of the actual and simulated gamma cameras were equal to 85.11 and 85.39 cps/MBq. Although these figures are not significantly different, but the small difference is due to acquisition dead time and signal overlap in the detector, which are not taken into account for the latter system. In [Figure 1], the backscattering pattern obtained for the simulated and real systems (Compton region) are very similar, confirming that the six cm Pyrex slab is a justifiable substitute for the whole structures attached to the back of NaI (Tl) crystal in the simulated system.
| Conclusion|| |
Comparison of acquired results by simulation and experiment is good evidence that SIMIND Monte Carlo simulation is capable of simulating the Siemens E.CAM gamma camera successfully. Rodrigues et al,  have also modeled a similar gamma camera with GATE.  The two codes are slightly different. Physical and geometrical parameters applied in these two studies are not exactly identical. However, taking into account the effect of these variations, the minor differences between two sets of results are negligible.
In summary, comparing the results acquired experimentally to those obtained by Monte Carlo simulation validate the simulated SPECT system. This step is a prerequisite for using the simulated system for complex studies on optimization of the performance characteristics, and finally improving the quality of SPECT produced images.
| Acknowledgment|| |
The authors are gratefully to the head and all staff members of the nuclear medicine department of Imam Reza hospital of Mashhad University of Medical Sciences for their sincere co-operation. This work was financially supported by the office of vice-president for research in Mashhad University of Medical Sciences, Iran.
| References|| |
|1.||Ljungberg M, Strand SE. A Monte Carlo program for the simulation of scintillation camera characteristics. Comp Meth Progr Biomed 1989;29:257-72. |
|2.||Ljungberg M. Simulation Techniques and Phantoms. Emission Tomography: The Fundamentals of PET and SPECT. In: Wernick M, Aarsvold J, editors. New York, USA: Academic Press Inc.; 2000. |
|3.||Zaidi H. Relevance of accurate Monte Carlo modeling in nuclear medical imaging. Med Phys 1999;26:574-608. |
|4.||Rosenthal MS, Cullom J, Hawkins W, Moore SC, Tsui BM, Yester M. Quantitative SPECT imaging: A review and recommendations by the Focus Committee of the Society of Nuclear Medicine Computer and Instrumentation Council. J Nucl Med 1995;36:1489-513. |
|5.||Zaidi H. Quantitative SPECT: Recent developments in detector response, attenuation and scatter correction techniques. Phys Med 1996;12:101-17. |
|6.||Zaidi H. Monte Carlo Modelling in Nuclear Medicine Imaging. In: Quantitative Analysis in Nuclear Medicine Imaging. Springer Sciences Business Media Inc.; 2006. p. 358-89. |
|7.||Lazaro D, Buvat I. Validation of the GATE Monte Carlo simulation platform for modelling a CsI(Tl) scintillation camera dedicated to small-animal imaging. Phys Med Biol 2004;49:271-85. |
|8.||Buvat I, Castiglioni I. Monte Carlo simulations in SPECT and PET. Q J Nucl Med 2002;46:48-61. |
|9.||Zaidi H, Koral KF. Scatter modeling and compensation in emission tomography. Eur J Nucl Med Mol Imaging 2004;31:761-82. |
|10.||Bice AN, Eary JF, Nelp WB. Quantification of iodine- 131 distribution by gamma camera imaging. Eur J Nucl Med 1991;18:142-8. |
|11.||Dewaraja YK, Ljungberg M, Koral KF. Accuracy of 131 I tumor quantification in radioimmunotherapy using SPECT imaging with an ultra-high-energy collimator: Monte Carlo study. J Nucl Med 2000;41:1760-7. |
|12.||Formiconi AR. Collimators. Q J Nucl Med 2002;46:8-15. |
|13.||Du Y, Frey EC, Wang WT, Tocharoenchai C, Baird WH, Tsui BM. Combination of MCNP and SimSET for Monte Carlo simulation of SPECT with medium- and high-energy photons. IEEE Trans Nucl Sci 2002;49:668-74. |
|14.||Dewaraja YK, Ljungberg M, Koral KF. Characterization of scatter and penetration using Monte Carlo simulation in 131 I imaging. J Nucl Med 2000;41:123-30. |
|15.||Lazaro D, Buvat I, Loudos G, Strul D, Santin G, Giokaris N, et al. Validation of the GATE Monte Carlo simulation platform for modelling a CsI(Tl) scintillation camera dedicated to small animal imaging. Phys Med Biol 2004;49:271-85. |
|16.||Staelens S, Strul D, Santin G, Staelens, Vandenberghe S, Koole M, et al. Monte Carlo simulations of a scintillation camera using GATE: Validation and application modeling. Phys Med Biol 2003;48:3021-42. |
|17.||Raeside DE. Monte Carlo principles and applications. Phys Med Biol 1976;21:181-97. |
|18.||Turner JE, Wright HA, Hamm RN. A Monte Carlo primer for health physicists. Health Phys 1985;48:717-33. |
|19.||James F. Monte Carlo theory and practice. Rep Prog Phys 1980;43:1145-89. |
|20.||Lund T. An introduction to the Monte Carlo method. Rep CERN/HS/067 CERN data handling division. Geneva: 1981. |
|21.||Bradshaw J, Burnham C, Correia J. Application of Monte Carlo methods to the design of SPECT detector systems. IEEE Trans Nucl Sci 1985;32:753-7. |
|22.||Lupton LR, Keller NA. Performance study of single-slice positron emission tomography scanners by Monte Carlo techniques. IEEE Trans Med Imaging 1983;2:154-68. |
|23.||Dewaraja YK, Ljungberg M, Koral K. Characterization of Scatter and Penetration Using Monte Carlo Simulation in 131 I Imaging. J Nucl Med 2000;1:123-30. |
|24.||Ljungberg M. The SIMIND Monte Carlo program. In: Monte Carlo calculation in nuclear medicine: Applications in diagnostic imaging. In: Ljungberg M, Strand SE, King MA, editors. Bristol and Philadelphia: IOP Publishing; 1998. p. 145-63. |
|25.||The SIMIND Monte Carlo Program. Available from: http://www.radfys.lu.se/simind. [cited in 2009]. |
|26.||Siemens Medical Solutions Home Page. Available from: http://www.medical.siemens.com. [cited in 2009]. |
|27.||Motta D, Schnert S. Optical properties of bialkali photocathodes. Nucl Instr Meth Phys Res A 2005;539:217-35. |
|28.||Data sheet of Siemens E.CAM Dual Head gamma camera. 1992. p. 1-12. |
|29.||Vittori F, Notaristefani F, Malatesta T. Crystals and light collection in Nuclear Medicine. Nucl Phys B (Proc Suppl) 1999;78:616-21. |
|30.||De Vries DJ, Moore SC, Zimmerman RE, Mueller SP, Friedland B, Lanza RC. Development and validation of a Monte Carlo simulation of photon transport in an anger camera. IEEE Trans Med Imaging 1990;9:430-8. |
|31.||Rotating Scintillation camera SPECT acceptance testing and quality control. AAPM Report No.22, American Institute of Physics, USA: 1987. |
|32.||National Electrical Manufacturers Association. Recommendations for implementing SPECT instrumentation quality control. J Nucl Med 2000;41:383-9. |
|33.||National Electrical Manufacturers Association (NEMA). Performance Measurements of Scintillation Cameras. NEMA NU 1-2001. Washington, USA; 2001. |
|34.||Quality Control of Nuclear Medicine Instruments. Vienna, Austria: IAEA; 1991. |
|35.||Quality Control of Gamma Camera Systems. IPEM 86, Institute of Physics and Engineering in Medicine, York, England: 2003. |
|36.||Nuclear Medicine Resources Manual IAEA, Vienna, Austria: 2006. |
|37.||Jarritt PH, Perkins AC, Woods SD. Audit of nuclear medicine scientific and technical standards. Nucl Med Commun 2004;25:771-5. |
|38.||Data sheet of Flangeless Deluxe Jaszczak SPECT Phantom. Available from: http://www.spect.com/pub/Flangeless_Jaszczak_Phantom.doc.pdf. [cited in 2002]. |
|39.||Groch MW, Erwin WD. Single-photon emission computed tomography in the year 2001: Instrumentation and quality control. J Nucl Med Technol 2001;29:9-15. |
|40.||Rodrigues S, Tome B, Abreu MC, Santos N, Rato Mendes P, Peralta L. Monte Carlo simulation and experimental characterization of a dual head gamma camera. 2007;14:2577. |
|41.||Holstensson M, Hindrof C, Ljungberg M, Partridge M, Flux GD. Optimization of energy-window setting for scatter correction in quantitative 111 in imaging: Comparison of measurements and Monte Carlo simulations. Can Biother Radiopharm 2007;22:136-42. |
|42.||Jan S, Santin G, Strul D, Staelens S, Assi΄e K, Autret D, et al. Gate: A simulation toolkit for PET and SPECT. Phys Med Biol 2004;49:4543-61. |
[Figure 1], [Figure 2], [Figure 3]
[Table 1], [Table 2]
|This article has been cited by|
||Attenuation correction of 111In planar images by use of dual energy, fundamental study by Monte Carlo simulation
| ||Seiji Shirakawa,Masanori Tadokoro,Hiroshi Hashimoto,Tomoya Ushiroda,Hiroshi Toyama |
| ||Radiological Physics and Technology. 2014; |
|[Pubmed] | [DOI]|
||Modulation transfer function assessment in parallel beam and fan beam collimators with square and cylindrical holes
| ||Abdollah Khorshidi,Mansour Ashoor |
| ||Annals of Nuclear Medicine. 2014; 28(4): 363 |
|[Pubmed] | [DOI]|
||Radius dependence of FP-CIT quantification: a Monte Carlo-based simulation study
| ||Walter Koch,Peter Bartenstein,Christian Fougère |
| ||Annals of Nuclear Medicine. 2013; |
|[Pubmed] | [DOI]|
||3D-OSEM and FP-CIT SPECT quantification
| ||Walter Koch,Peter Bartenstein,Christian la Fougère |
| ||Nuclear Medicine Communications. 2013; 34(10): 971 |
|[Pubmed] | [DOI]|
||Development of GATE Monte Carlo simulation for a dual-head gamma camera
| ||Mehdi Momennezhad,Ramin Sadeghi,Shahrokh Nasseri |
| ||Radiological Physics and Technology. 2012; 5(2): 222 |
|[Pubmed] | [DOI]|
||Evaluation of collimatorsæ response: Round and hexagonal holes in parallel and fan beam
| ||Khorshidi, A. and Ashoor, M. and Hamed Hosseini, S. and Rajaee, A. |
| ||Progress in Biophysics and Molecular Biology. 2012; 109(3): 59-66 |
||Development of GATE Monte Carlo simulation for a dual-head gamma camera
| ||Momennezhad, M. and Sadeghi, R. and Nasseri, S. |
| ||Radiological Physics and Technology. 2012; 5(2): 222-228 |
||Simulation of a quality control jaszczak phantom with SIMIND Monte Carlo and adding the phantom as an accessory to the program
| ||Islamian, J.P. and Toossi, M.T.B. and Momennezhad, M. and Naseri, S. and Ljungberg, M. |
| ||Iranian Journal of Medical Physics. 2012; 9(2 A): 135-140 |