Abstract
Developing image reconstruction algorithms for diagnostic medical devices requires physically accurate and effective simulation tools. In this paper we present a hybrid Monte Carlo (MC) particle simulation method for Computed Tomography (CT) scanners. To meet the performance requirements, we combine several variance reduction techniques and tailor the algorithms for effective GPU execution. Variance reduction methods include main part separation, sample weighting, reuse, forced collision, next event estimation and table driven importance sampling. We show that the resulting method can deliver accurate simulations orders of magnitude faster than direct physical simulation.
Access this chapter
Tax calculation will be finalised at checkout
Purchases are for personal use only
References
Euclid, S.: Computed tomography: physical principles and recent technical advances. J. Med. Imaging Radiat. Sci. 41(2), 87–109 (2010). doi:10.1016/j.jmir.2010.04.001. http://www.sciencedirect.com/science/article/pii/S1939865410000317. ISSN 1939–8654
Légrády, D., Cserkaszky, A., Wirth A., Domonkos B.: PET image reconstruction with on the fly Monte Carlo using GPU. In: Proceedings of PHYSOR 2010, American Nuclear Society, Pittsburgh, Pennsylvania (2010)
Wirth, A., Cserkaszky, A., Kári, B., Légrády, D., Fehér, S., Czifrus, S., Domonkos, B.: Implementation of 3D Monte Carlo PET reconstruction algorithm on GPU. In: IEEE Nuclear Science Symposium Conference Record (NSS/MIC), pp. 4106–4109 (2009)
Woodcock E., Murphy T., Hemmings P., Longworth S.: Techniques used in the GEM code for Monte Carlo neutronics calculation. In: Proceedings of the Conference on Applications of Computing Methods to Reactors, ANL-7050 (1965)
Szirmay-Kalos, L., Tóth, B., Magdics, M., Légrády, D., Penzov, A.: Gamma photon transport on the GPU for PET. In: Lirkov, I., Margenov, S., Waśniewski, J. (eds.) LSSC 2009. LNCS, vol. 5910, pp. 435–442. Springer, Heidelberg (2010)
Szirmay-Kalos, L., Tóth, B., Magdics, M.: Free path sampling in high resolution inhomogeneous participating media. Comput. Graph. Forum 30(1), 85–97 (2011)
Yang, C.N.: The Klein-Nishina formula and quantum electrodynamics. In: Yang, C.N. (ed.) Nishina Memorial Lectures. Lecture Notes in Physics, pp. 393–397. Springer, Berlin (2008)
Jakab, G., Rácz, A., Nagy, K.: High quality cone-beam CT reconstruction on the GPU. In: 8th KÉPAF Conference, Budapest, Hungary (2011)
Magdics, M., Szirmay-Kalos, L., Tóth, B., Csendesi, Á., Penzov, A.: Scatter estimation for PET reconstruction. In: Dimov, I., Dimova, S., Kolkovska, N. (eds.) NMA 2010. LNCS, vol. 6046, pp. 77–86. Springer, Heidelberg (2011)
Jakab, G., Huszár T., Csébfalvi, B.: Iterative CT Reconstruction on the GPU. In: Sixth Hungarian Conference on Computer Graphics and Geometry, Budapest (2012)
Acknowledgement
This work has been supported by the OTKA K-104476 and by TÁMOP -4.2.2.B-10/1–2010-0009.
Author information
Authors and Affiliations
Corresponding author
Editor information
Editors and Affiliations
Rights and permissions
Copyright information
© 2014 Springer-Verlag Berlin Heidelberg
About this paper
Cite this paper
Jakab, G., Szirmay-Kalos, L. (2014). Hybrid Monte Carlo CT Simulation on GPU. In: Lirkov, I., Margenov, S., Waśniewski, J. (eds) Large-Scale Scientific Computing. LSSC 2013. Lecture Notes in Computer Science(), vol 8353. Springer, Berlin, Heidelberg. https://doi.org/10.1007/978-3-662-43880-0_17
Download citation
DOI: https://doi.org/10.1007/978-3-662-43880-0_17
Published:
Publisher Name: Springer, Berlin, Heidelberg
Print ISBN: 978-3-662-43879-4
Online ISBN: 978-3-662-43880-0
eBook Packages: Computer ScienceComputer Science (R0)