Design feasibility of gamma camera without collimator based on specific arrangement of the detectors

Design feasibility of gamma camera without collimator based on specific arrangement of the detectors 1 2 corresp first-author Corresponding author: Email: sashafiei@rums.ac.ir Corresponding author: Email: sashafiei@rums.ac.ir sashafiei@rums.ac.ir 3 Assistant Professor of Medical Physics, School of Medicine, Rafsanjan University of Medical Sciences, Rafsanjan, Iran Assistant Professor of Medical Physics, School of Medicine, Rafsanjan University of Medical Sciences, Rafsanjan, Iran Assistant Professor of Medical Physics, School of Medicine, Qom University of Medical Sciences, Qom, Iran Assistant Professor of Medical Physics, School of Medicine, Qom University of Medical Sciences, Qom, Iran Department of Medical Physics, Mashhad University of Medical Sciences, Mashhad, Iran Department of Medical Physics, Mashhad University of Medical Sciences, Mashhad, Iran In gamma camera and single-photon emission computerized tomography (SPECT), the collimator removes most photons. Here, a gamma camera without collimator utilizes a specific arrangement of detectors. Instead of bending beams (like a lens) or directing beams (by parallel hole collimator), changes are created in detectors’ field of view (FOV), so that each detector’s FOV looks different from others. Simulation proved this theory, with 98 detectors (2cm ×1.41cm) arranged in a zigzag manner for Monte Carlo simulation. A radioactive source with energy of 140 keV was situated on the detectors’ faces. Sixty projections, each 3° (0–179°) apart, were simulated by Monte Carlo N-Particle (MCNP) 4C code, rotating detectors around a radioactive point. The band containing the radioactive source is clearly visible in each projection. Counts obtained after simulation in different projections were reconstructed, and point source location emerged correctly. Simulation of gamma camera with zigzag arrangement of detectors and MCNP-4C code demonstrated that one could string the space and determine radioactive source by image reconstruction without using collimators, solely through these special detectors’ distribution. Collimator Gamma camera SPECT Monte Carlo simulation 1 Introduction 1 To be able to photograph an object, there should be a one-to-one correspondence between an object and its image. The lens performs this function in an ordinary photographic camera. However, a lens cannot create such a correspondence for high-energy radiation, since lenses ordinarily are not able to bend high-energy radiation [ 1 1 ]. Therefore, nuclear imaging (SPECT and gamma camera) is conducted with the help of a radioactive uses aperture (pinhole) and parallel holecollimator (PHC) to direct output radiation from the body. In this way, a one-to-one correspondence is created between radiopharmaceutical distribution in the body and its image [ 2 2 ]. Many studies[ 3 3 - 6 6 3-6 ] have been conductedinthe effortto find appropriate methods for design and construction of collimators of gamma-ray imaging systems. Three such influential studies arethe Jaszczak et al. study assessing cone beam collimators [ 4 4 ]; the Tsui et al. researchon the assessment of fan beam collimators [ 5 5 ]; and that of Mueller et al., which selected high-resolution low-energy collimators in brain imaging [ 6 6 ]. One type of collimator is the pinhole collimator (PC), comprised of a small pinhole aperture embedded in a piece of lead or tungsten. In PCs, where there is a small object next to the collimator aperture, there usually is good spatial resolution and detection efficiency, but the field of view (FOV) becomes smaller [ 7 7 ]. Such collimators are used in humans to image some body parts such as the wrist and thyroid [ 7 7 , 8 8 ], and are only rarely used for heart and breast [ 9 9 , 10 10 ]. The greater the distance between the object and the PC’s aperture is, the greater the FOV is, and the less the spatial resolution and sensitivity is. On the other hand, the sens itivity will be increased by using a pinhole aperture of large diameter that leads to a decrease in resolution. The more the apertures focus on different points in the FOV, the larger the FOV. Additionally, the higher the sensitivity, the better the spatial resolution [ 11 11 ]. Many studies select the proper parameters to improve the above-mentioned quantities [ 12 12 , 13 13 ]. However, using single aperture eliminates most photons containing information. The sensitivity, FOV, and spatial resolution will rise if apertures are increased [ 14 14 - 16 16 14-16 ]. Therefore, one of the methods employed to decrease the number of removed photons containing information is to use a thick lead plate with holes (PHC), in which all beams parallel to its holes are allowed to cross [ 2 2 ]. Designing the proper PHC depends on the application, sensitivity, and necessary resolution. Hence, it is constructed in different forms and designs and the real performance of the best design will usually be assessed empirically [ 17 17 , 18 18 ]. The major problem with PHCs remains the elimination of useful information-containing beams, i.e., of every 10,000 photons, just one is allowed to cross the collimator [ 19 19 ]. This leads to an increase in imaging time and patient-injected radiopharmaceutical dose, and a decrease in image quality. Since the number of emitted photons is generally very low, collimators have the most important restrictions in gamma camera and SPECT. Consequently, a great step forward will be taken if a radiopharmaceutical spatial distribution in the body is obtained by using a gamma camera without a collimator. The aim of this study was to present the idea of a gamma camera without collimator (GCWC) based on the specific arrangement of detectors in zigzag form. Ideally, this theory will be assessed through Monte Carlo simulation by MCNP code. 2 Material and Methods 1 2.1 Theory 2 In GCWC, instead of using a series of detectors on a flat plate together with a large crystal, gamma detectors are distributed in zigzag form to create a row of detectors (ROD) ( Figure 1 Figure 1 ). Detector rows are separated with lead plates, called septa. In other words, instead of bending the beams (like a lens) or directing the beams (by collimators), some changes are created in the detectors’ FOV in such a way that the FOV of each detector will be different from those of the other detectors ( Figure 2 Figure 2 ). In this way, the observation of radiopharmaceutical distribution in space will be unique for different detectors. The image can be produced by using the difference in observed radiopharmaceutical distribution and image reconstruction through the filtered back projection (FBP) method, without the necessity of a collimator. GCWC is designed in such a manner that distribution of detectors is shaped in a zigzag manner. As a result, the FOV of each detector on the right side of each triangle in Figure 2 Figure 2 is different from the adjacent detector by merely a strip (similarly, this is true for the detectors on the left side of the triangle). In other words, detectors D 1 and D 2 observe radioactive source S´ but only detector D 1 (out of the detectors on the right side) detects source S . Therefore, if a radioactive material is laid in front of this detector distribution (each ROD), N D1 –N D2 equals the absorbed radiation resulting from radioactive materials in the shaded bar ( N D1 and N D2 are respectively the number of counts read by detectors D 1 and D 2 ). In other words, a narrow strip containing radioactive material will be separated without the necessity of eliminating beams through using hardware similar to a PHC. The imaging is staged in two perpendicular projections; point S is shaped from the intersection of two bands through (D 0 i -D 0 i +1 ) and (D 90 j -D 90 j +1 ) ( Fig. 3 Fig. 3 ). Therefore, point S can be determined through image reconstruction by FBP. D φα indicates the α th detector in the projection angle of φ° . Apparently the proposed idea can be implemented theoretically, but its accuracy and verification must be proven through simulation. Therefore, the following sections will investigate the feasibility of this idea, ideally through simulation based on Monte Carlo with MCNP code. 2.2 Computer Simulation 2 Monte Carlo method is based on random sampling to find solutions for physics or math problems according to certain probabilities and laws of statistics. In addition, it is highly accurate in predicting the outcomes of processes that are repeated many times. Monte Carlo simulation (MCS) using MCNP code is a common method in gamma camera imaging and designing proper collimator experimentation [ 18 18 , 20 20 ]. Ninety-eight detectors were arranged in a zigzag manner on a line 1m in length, which formed a ROD for MCS. Dimensions of each detector’s surface were 2cm×1.41cm. Therefore, the thickness of the imaging slice that was taken by ROD was 2cm. Technetium-99m (which has an energy of 140 keV) radiopharmaceuticals are in widespread use for imaging and functional studies (over 80% of all diagnostic procedures) of major organs such as the heart, brain, liver, kidney, and thyroid [ 21 21 ]. Therefore, a radioactive source with energy of 140 keV was situated on the face of the ROD, i.e., the slice on which the ROD and some radioactive material have been placed ( Figs. 2 Figs. 2 and 3 3 ). Among the active detectors, 49 detectors were laid on the right side of the triangle in ROD. Sixty projections at a distance of 3° (0–179°) from each other were simulated by MCNP-4C code. In other words, simulation was performed 60 separate times, with 60 different projections at a distance of 3°, and the photons passing through each detector were counted.It was verified that in each simulation period, 100,000,000 photons were released from the radioactive source and were traced by MCS. 3 Results 1 After simulation and data collection, each detector’s counts were estimated in different projections. For example, Figure 4 Figure 4 shows the first projection (angle of 0°) in which the vertical axis is the gamma ray counts and the horizontal axis is the number of detectors on the right side of the zigzag arrangement. In Figure 2 Figure 2 , the number of counts necessary for the relevant strip will greatly decrease with the angle between the incident radiation and detector surface. Therefore, Eq. (1) observes the number of counts resulting from the difference counts of two consecutive detectors. D φα indicates counts of the α th detector in the projection angle of φ° . Two consecutive detectors determine the bandwidthvalue of each band ( d φα ), calculated by Eq. (1); in addition, ρ and μ are arbitrary values for reinforcing these differences. <math display="block" id="M2"><mrow><msub><mi>d</mi><mrow><mi>φ</mi><mi>α</mi></mrow></msub><mo>=</mo><mi>ρ</mi><mo>*</mo><msup><mrow><mrow><mo>(</mo><mrow><mfrac><mrow><msub><mi>D</mi><mrow><mi>φ</mi><mi>α</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>−</mo><msub><mi>D</mi><mrow><mi>φ</mi><mi>α</mi></mrow></msub></mrow><mrow><msub><mi>D</mi><mrow><mi>φ</mi><mi>α</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>+</mo><msub><mi>D</mi><mrow><mi>φ</mi><mi>α</mi></mrow></msub></mrow></mfrac></mrow><mo>)</mo></mrow></mrow><mi>μ</mi></msup></mrow></math> Figure 5 Figure 5 is the result of the applying Eq. (1) on the collected counts obtained from different detectors in the first projection (at 0°). The band containing the radioactive source is clearly visible. The counts obtained after simulation for each detector in different projections were reconstructed after being transferred to MATLAB software by the command “iradon”, based on the FBP technique; the final image is illustrated in Figure 6 Figure 6 . An imaginary point source was also placed within a uniform matrix through using MATLAB software, and then 60 projections were provided by the command “radon”. Then, the mentioned point was reconstructed by the command “iradon” ( Figure6A Figure6A ). As can be seen in Figure 6 Figure 6 , the point image emerged as plus (+) after reconstruction from its projections due to the limitation of the command “iradon” in MATLAB software.This software uses the back projection method for image reconstruction based on projections; the main problem with this method is the star artifact that can be seen in Figure6A Figure6A . All photons received by the detectors have energy of 140 keV due to the lack of scattering in the environment. Thus, spectrum evaluation of detected photons’ energy was not possible and deemed meaningless. Sensitivity is defined as the number of counts passing through the detectors per second per megabecquerel of source (cps/MBq). Geometric sensitivity of ROD in each projection was different with respect to other projections, because radioactive source position varied in the face of ROD in each projection. Mean, maximum, and minimum sensitivity in 60 projections are presented in Table 1 Table 1 . Sensitivity was calculated before and after subtracting two consecutive detectors by Eq. (1). SD = standard deviation. 1 1   1 1 Mean±SD (cps/MBq) 1 1 Max (cps/MBq) 1 1 Min (cps/MBq) 1 1 Before subtraction 1 1 (4.1±0.72)×10 3 1 1 5.4×10 3 1 1 3.3×10 3 1 1 After subtraction 1 1 29.5.6±96.3 1 1 540 1 1 203.8 4 Discussion 1 Sixty simulations were conducted by rotation of detectors around the radioactive point in 60 projections using MCNP-4C code; then data were transferred, differences were counted with Eq. (1) by MATLAB software, and the space was threaded in each projection. Figure 6B Figure 6B is the outcome of this action. The star pattern around the point source results from the image reconstruction process [ 22 22 ], as a sharp curve appeared after subtracting ( Figure 5 Figure 5 ), and the band containing radioactive source was clearly visible. Proving this issue, a zero matrix was built in order to determine the contribution of blurring the image in running the command “iradon” in MATLAB software for the image reconstruction process. Digit 1 was placed as one of its elements, similar to radioactive source, in front of the detectors. After that, ideal projections were made by the command “radon” and the point’s image was reconstructed by the command “iradon”. Figure 6A Figure 6A is the outcome of this action. The comparison made between this image and the reconstructed image through 60 projections in MCNP-4C code ( Figs. 6B and 6C Figs. 6B and 6C ) indicates that lack of clarity in the image is related to running the command “iradon”. As a result of the comparison made between point T, resulting from the simulation with point + due to the command “radon”, the hypothesis of the elimination of collimators through detectors’ zigzag distribution ideally is verified, even though the radioactivity in the reconstructed images can be achieved without shadow by edge filters (e.g., Roberts filter). Regarding the study’s aim (investigating the possibility of omitting collimator and observing a radioactive source by a special distribution of detectors), the simulation of various sources at different distances from the camera or the distribution of radioactive material has not been dealt with and will be a future research topic. MCS and practical construction of slit-slat collimators have been investigated for resolution and sensitivity along the axis and perpendicular to the axis as a function of aperture [ 23 23 , 24 24 ]. Saad et al.also used MCS to find the best mesh in PHC for breast imaging [ 25 25 ]. Islamian et al. have applied MCS to investigate the suitable thickness of collimators [ 26 26 ]. Since the real structure of the detector was not considered in the simulations, the crossing of each photon from the small plate—considered a detector—was recorded, and the beam was assumed to have been completely absorbed. Therefore, the scattering due to photons’ interaction with the detector components was not taken into account; instead, it was assumed that each beam crossing from the detector surface was recorded by the detector, put aside, and could not penetrate adjacent detectors. In fact, penetration of beam to adjacent detectors occurs and should be considered. Penetration and scattered photons in the collimator have a major impact on nuclear medicine imaging [ 27 27 - 29 29 27-29 ]. The goal behind using collimators is to thread the space in each projection and to determine the bands including radioactive points, and then to obtain the distribution of radioactive material through intersection of stripped spaces by different projections. The main problem in using collimators is that most of the radiation containing information is eliminated for threading. However, in the above-mentioned method, the detectors receive the majority of information-containing radiations. As a result, the sensitivity of GCWC ( Table 1 Table 1 ) is higher than PHC and PCs [ 30 30 - 32 32 30-32 ].Therefore, according to Eq. (1), through the difference between the two detectors’ data, a string can be created that can be seen by one detector and not by the other detectors ( Fig. 2 Fig. 2 ). Apparently there is no difference between the presence of collimators or the above-mentioned method, and sensitivity after subtracting ( Table 1 Table 1 ) is approximately similar to conventional gamma camera [ 16 16 , 33 33 ]. For example, sensitivity was obtained with 10 cameras used for Tc-99m thyroid scintigraphy by Seret et al. [ 31 31 ]. Measured minimum-maximum values of sensitivity were 62-131 cps/MBq for six pinhole cameras and 55-85 cps/MBq for four non-pinhole cameras. Additionally, Fleming etal. [ 32 32 ] reported sensitivity of two kinds of SPECT, using 99m Tc at 120.2 and 153 cps/MBq.Nevertheless, it is worth mentioning that the strings made up of two adjacent detectors produce the maximum spatial resolution. The image quality can be improved by making strings at the distance of two detectors or more; besides, it is possible to correct the final image through combining several images with various strings’ widths. When two adjacent detectors observe the radioactive material clearly, the counts of these two detectors will rise. Therefore, their difference may be much greater than the difference between two border detectors (one discerns radioactive source but the other does not) with just a few counts. Each difference will be divided by the sum of counts of two detectors in order to reinforce the difference between the two small values relative to that of the two large values. The point will be observed vividly in this stage, although it is possible to affect border reinforcement or final image clarification by changing the coefficient or power in Eq. (1). For example, in Figure 6 Figure 6 the effect of selecting ( ρ =1, μ =2) is noted compared with ( ρ =1, μ =1). Clearly, if several RODs are separated by septa, detectors will be distributed in two dimensions ( Figure 1 Figure 1 ). Therefore, a volumetric image will be formed instead of a sectional one when the camera rotates around the patient. If two parallel heads are used instead of a one-head system, increasing image-making potential, the applied corrections and imaging speed during recording will be observed. As a result, some corrections can be imagined for GCWC, which could inspire future research. The topics could well include optimizing the number and size of detectors, suitable angle of detectors, taking account of the corrections due to the influence of the radiations from adjacent detectors, and absorption of scattered radiations. 5 Conclusion 1 Simulation of gamma camera with arrangement of zigzag detectors using MCNP-4C code in 60 projections has revealed that the space can be threaded and the radioactive source can be determined through image reconstruction without using collimator, and caninsteadbe accomplished through a special distribution of the detectors. Simulation was conducted ideally to prove the main idea; issues such as scattering of beam due to striking detector components, the beam passing through a detector and recording in another detector, etc., are related to photons’ collision with the material in a real location, which can affect image quality and have not been considered in the simulation. These issues should be taken into account for simulation of real GCWC in the future. References M.N. Wernick , J.N. Aarsvold , Emission tomography: the fundamentals of PET and SPECT . Amsterdam; Boston : Elsevier Academic Press , 2004 , 65 - 75 . Emission tomography: the fundamentals of PET and SPECT R.A. Powsner , E.R. Powsner , Essential nuclear medicine physics . Malden, Mass. : Blackwell Pub ., 2006 , 16 - 17 . Essential nuclear medicine physics C. Robert , G. Montémont , V. Rebuffel et al ., Optimization of a parallel hole collimator/CdZnTe gamma-camera architecture for scintimammography . Med. Phys ., 2011 , 38 , 1806 - 1819 . doi: 10.1118/1.3560423 http://doi.org/10.1118/1.3560423 Optimization of a parallel hole collimator/CdZnTe gamma-camera architecture for scintimammography R.J. Jaszczak , K.L. Greer , R.E. Coleman , SPECT using a specially designed cone beam collimator . J Nucl. Med ., 1988 , 29 , 1398 - 1405 . SPECT using a specially designed cone beam collimator B.M.W. Tsui , G.T. Gullberg , E.R. Edgerton , Design and clinical utility of a fan beam collimator for SPECT imaging of the head . J. Nucl. Med ., 1986 , 27 , 810 - 819 . Design and clinical utility of a fan beam collimator for SPECT imaging of the head S.P. Mueller , J.F. Polak , M.F. Kijewski et al ., Collimator selection for SPECT brain imaging: The advantage of high resolution . J. Nucl. Med ., 1986 , 27 , 1729 - 1738 . Collimator selection for SPECT brain imaging: The advantage of high resolution B. Sivertsen , S. Omvik , S. Pallesen et al ., Cognitive behavioral therapy vs zopiclone for treatment of chronic primary insomnia in older adults: a randomized controlled trial . JAMA , 2006 , 295 , 2851 - 2858 . doi: 10.1001/jama.295.24.2851 http://doi.org/10.1001/jama.295.24.2851 Cognitive behavioral therapy vs zopiclone for treatment of chronic primary insomnia in older adults: a randomized controlled trial J.P. Soeffing , K.L. Lichstein , S.D. Nau et al ., Psychological treatment of insomnia in hypnotic-dependant older adults . Sleep Med ., 2008 , 9 , 165 - 171 . doi: 10.1016/j.sleep.2007.02.009 http://doi.org/10.1016/j.sleep.2007.02.009 Psychological treatment of insomnia in hypnotic-dependant older adults P. Paredes , S. Vidal-Sicart , G. Zanón et al ., Radioguided occult lesion localisation in breast cancer using an intraoperative portable gamma camera: first results . Eur. J. Nucl. Med. Mol. Imaging , 2008 , 35 , 230 - 235 . doi: 10.1007/s00259-007-0640-x http://doi.org/10.1007/s00259-007-0640-x Radioguided occult lesion localisation in breast cancer using an intraoperative portable gamma camera: first results M.A. Ghanem , A.H. Elgazzar , M.M. Elsaid et al ., Comparison of pinhole and high-resolution parallel-hole imaging for nodular thyroid disease . Clin. Nucl. Med ., 2011 , 36 , 770 - 771 . doi: 10.1097/RLU.0b013e318219ac29 http://doi.org/10.1097/RLU.0b013e318219ac29 Comparison of pinhole and high-resolution parallel-hole imaging for nodular thyroid disease W. Branderhorst , B. Vastenhouw , F. Van Der Have et al ., Targeted multi-pinhole SPECT . Eur. J. Nucl. Med. Mol. Imaging , 2011 , 38 , 552 - 561 . Targeted multi-pinhole SPECT Q.S. Dai , Channel model of pinhole gamma camera . Chinese Sci. Bull ., 2011 , 56 , 2758 - 2763 . doi: 10.1007/s11434-011-4606-9 http://doi.org/10.1007/s11434-011-4606-9 Channel model of pinhole gamma camera J. Mejia , O.Y. Galvis-Alonso , A.A. de Castro , A clinical gamma camera-based pinhole collimated system for high resolution small animal SPECT imaging . Braz. J. Med. Biol. Res ., 2010 , 43 , 1160 - 1166 . doi: 10.1590/S0100-879X2010007500128 http://doi.org/10.1590/S0100-879X2010007500128 A clinical gamma camera-based pinhole collimated system for high resolution small animal SPECT imaging F.P. DiFilippo , Geometric characterization of multi-axis multi-pinhole SPECT . Med. Phys ., 2008 , 35 , 181 - 194 . doi: 10.1118/1.2821703 http://doi.org/10.1118/1.2821703 Geometric characterization of multi-axis multi-pinhole SPECT S. Shokouhi , S.D. Metzler , D. Wilson et al ., Multi-pinhole collimator design for small-object imaging with SiliSPECT: a high-resolution SPECT . Phys. Med. Biol ., 2009 , 54 , 207 - 225 . doi: 10.1088/0031-9155/54/2/003 http://doi.org/10.1088/0031-9155/54/2/003 Multi-pinhole collimator design for small-object imaging with SiliSPECT: a high-resolution SPECT C. Vanhove , M. Defrise , T. Lahoutte et al ., Three-pinhole collimator to improve axial spatial resolution and sensitivity in pinhole SPECT . Eur. J. Nucl. Med. Mol. Imaging , 2008 , 35 , 407 - 415 . doi: 10.1007/s00259-007-0579-y http://doi.org/10.1007/s00259-007-0579-y Three-pinhole collimator to improve axial spatial resolution and sensitivity in pinhole SPECT M. Rasouli , A. Takavar , M.R. Ay et al ., Effects of crystal pixel size and collimator geometry on the performance of a pixelated crystal γ-camera using Monte Carlo simulation . J. Nucl. Med. Technol ., 2010 , 38 , 199 - 204 . doi: 10.2967/jnmt.109.073072 http://doi.org/10.2967/jnmt.109.073072 Effects of crystal pixel size and collimator geometry on the performance of a pixelated crystal γ-camera using Monte Carlo simulation I. Buvat , S. Laffont , J. Le Cloirec et al ., Importance of the choice of the collimator for the detection of small lesions in scintimammography: a phantom study . Phys. Med. Biol ., 2001 , 46 , 1343 - 1355 . doi: 10.1088/0031-9155/46/5/301 http://doi.org/10.1088/0031-9155/46/5/301 Importance of the choice of the collimator for the detection of small lesions in scintimammography: a phantom study M.N. Wernick , J.N. Aarsvold . Emission tomography: the fundamentals of PET and SPECT . Amsterdam; Boston : Elsevier Academic Press , 2004 , 153 - 163 Emission tomography: the fundamentals of PET and SPECT M. Rafecas , G. Boning , B.J. Pichler et al ., Inter-crystal scatter in a dual layer, high resolution LSO-APD positron emission tomograph . Phys. Med. Biol ., 2003 , 48 , 821 - 848 . Inter-crystal scatter in a dual layer, high resolution LSO-APD positron emission tomograph International Atomic Energy A , Technetium-99m radiopharmaceuticals: Manufacture of kits . Vienna : International Atomic Energy Agency , 2008 , 1 - 9 . Technetium-99m radiopharmaceuticals: Manufacture of kits J.T. Bushberg , The essential physics of medical imaging . Philadelphia : Lippincott Williams & Wilkins , 2002 , 37 - 373 . The essential physics of medical imaging S.D. Metzler , R. Accorsi , J.R. Novak et al ., On-axis sensitivity and resolution of a slit-slat collimator . J. Nucl. Med ., 2006 , 47 , 1884 - 1890 . On-axis sensitivity and resolution of a slit-slat collimator A. Kamali-Asl , S. Sarkar , M. Shahriari et al ., Slit slat collimator optimization with respect to MTF . Appl. Radiat. Isot ., 2005 , 62 , 461 - 468 . doi: 10.1016/j.apradiso.2004.07.009 http://doi.org/10.1016/j.apradiso.2004.07.009 Slit slat collimator optimization with respect to MTF W.H.M. Saad , R.E. Roslan , M.A. Mahdi et al ., Monte Carlo design of optimal wire mesh collimator for breast tumor imaging process . Nucl. Instr. Meth. Phys. Res. A , 2011 , 648 , 254 - 260 . doi: 10.1016/j.nima.2011.05.064 http://doi.org/10.1016/j.nima.2011.05.064 Monte Carlo design of optimal wire mesh collimator for breast tumor imaging process J.P. Islamian , M.T. Toossi , M. Momennezhad et al ., Monte Carlo study of the effect of collimator thickness on T-99m source response in single photon emission computed tomography . World J. Nucl. Med ., 2012 , 11 , 70 - 74 . doi: 10.4103/1450-1147.103419 http://doi.org/10.4103/1450-1147.103419 Monte Carlo study of the effect of collimator thickness on T-99m source response in single photon emission computed tomography D. Autret , A. Bitar , L. Ferrer et al ., Monte Carlo modeling of gamma cameras for I-131 imaging in targeted radiotherapy . Cancer Biother. Radio ., 2005 , 20 , 77 - 84 . doi: 10.1089/cbr.2005.20.77 http://doi.org/10.1089/cbr.2005.20.77 Monte Carlo modeling of gamma cameras for I-131 imaging in targeted radiotherapy A. Larsson , M. Ljungberg , S.J. Mo et al ., Correction for scatter and septal penetration using convolution subtraction methods and model-based compensation in 123I brain SPECT imaging-a Monte Carlo study . Phys. Med. Biol ., 2006 , 51 , 5753 - 5767 . doi: 10.1088/0031-9155/51/22/003 http://doi.org/10.1088/0031-9155/51/22/003 Correction for scatter and septal penetration using convolution subtraction methods and model-based compensation in 123I brain SPECT imaging-a Monte Carlo study M.F. Smith , R.J. Jaszczak , The effect of gamma ray penetration on angle-dependent sensitivity for pinhole collimation in nuclear medicine . Med. Phys ., 1997 , 24 , 1701 - 1709 . The effect of gamma ray penetration on angle-dependent sensitivity for pinhole collimation in nuclear medicine F. Beekman , F. van der Have , The pinhole: gateway to ultra-high-resolution three-dimensional radionuclide imaging . Eur. J. Nucl. Med. Mol. Imaging , 2007 , 34 , 151 - 161 . doi: 10.1007/s00259-006-0248-6 http://doi.org/10.1007/s00259-006-0248-6 The pinhole: gateway to ultra-high-resolution three-dimensional radionuclide imaging A.A. Seret , A comparison of contrast and sensitivity in Tc-99m thyroid scintigraphy between nine nuclear medicine centres of a geographic area . Alasbimn Journal , 2006 , 32 , 32 - 33 . A comparison of contrast and sensitivity in Tc-99m thyroid scintigraphy between nine nuclear medicine centres of a geographic area J.S. Fleming , A.S. Alaamer , Influence of collimator characteristics on quantification in SPECT . J. Nucl. Med ., 1996 , 37 , 1832 - 1836 . Influence of collimator characteristics on quantification in SPECT M. Rasouli , A. Takavar , M.R. Ay et al ., Effects of crystal pixel size and collimator geometry on the performance of a pixelated crystal gamma-camera using Monte Carlo simulation . J. Nucl. Med. Technol ., 2010 , 38 , 199 - 204 . doi: 10.2967/jnmt.109.073072 http://doi.org/10.2967/jnmt.109.073072 Effects of crystal pixel size and collimator geometry on the performance of a pixelated crystal gamma-camera using Monte Carlo simulation The authors would like to thank S. Skies and S.A. Sabbagh for editing this manuscript. 10.1007/s41365-016-0088-4 This work was supported by grants from Rafsanjan University of Medical Sciences, Rafsanjan, Iran. 2015-10-28 2016-05-10 2016-05-11 2016-07-11 2016-08-20 2021-04-11 © China Science Publishing & Media Ltd. (Science Press), Shanghai Institute of Applied Physics, the Chinese Academy of Sciences, ChineseNuclear Society and Springer Nature Singapore Pte Ltd. 2016 This is a post-peer-review, pre-copyedit version of an article published in Nuclear Science and Techniques. The final authenticated version is available online at: http://dx.doi.org/10.1007/s41365-016-0088-4 http://dx.doi.org/10.1007/s41365-016-0088-4 . 2016 Nuclear Science and Techniques 04 27