Mesh modeling of system geometry and anatomy phantoms for realistic GATE simulations and their inclusion in SPECT reconstruction. (7th April 2023)

• Record Type: Journal Article
• Authors: Auer, Benjamin; Könik, Arda; Fromme, Timothy J; De Beenhouwer, Jan; Kalluri, Kesava S; Lindsay, Clifford; Furenlid, Lars R; Kuo, Philip H; King, Michael A
• Is Part Of: Physics in medicine & biology. Volume 68:Number 7(2023)
• Publication Date: 2023-04-07
• Subjects: GATE simulation -- computer aided design software -- triangulated mesh -- XCAT phantom -- simulation of complex system and phantom geometries -- brain SPECT imaging -- adaptispect-C; Biophysics -- Periodicals; Medical physics -- Periodicals; 610.153
• DOI: 10.1088/1361-6560/acbde2
• Languages: English
• ISSNs: 0031-9155

Abstract

Abstract: Objective. Monte-Carlo simulation studies have been essential for advancing various developments in single photon emission computed tomography (SPECT) imaging, such as system design and accurate image reconstruction. Among the simulation software available, Geant4 application for tomographic emission (GATE) is one of the most used simulation toolkits in nuclear medicine, which allows building systems and attenuation phantom geometries based on the combination of idealized volumes. However, these idealized volumes are inadequate for modeling free-form shape components of such geometries. Recent GATE versions alleviate these major limitations by allowing users to import triangulated surface meshes. Approach. In this study, we describe our mesh-based simulations of a next-generation multi-pinhole SPECT system dedicated to clinical brain imaging, called AdaptiSPECT-C. To simulate realistic imaging data, we incorporated in our simulation the XCAT phantom, which provides an advanced anatomical description of the human body. An additional challenge with the AdaptiSPECT-C geometry is that the default voxelized XCAT attenuation phantom was not usable in our simulation due to intersection of objects of dissimilar materials caused by overlap of the air containing regions of the XCAT beyond the surface of the phantom and the components of the imaging system. Main results. We validated our mesh-based modeling against the one constructed by idealized volumes for a simplified single vertex configuration of AdaptiSPECT-C through simulated projection data of 123 I-activity distributions. We resolved the overlap conflict by creating and incorporating a mesh-based attenuation phantom following a volume hierarchy. We then evaluated our reconstructions with attenuation and scatter correction for projections obtained from simulation consisting of mesh-based modeling of the system and the attenuation phantom for brain imaging. Our approach demonstrated similar performance as the reference scheme simulated in air for uniform and clinical-like 123 I-IMP brain perfusion source distributions. Significance. This work enables the simulation of complex SPECT acquisitions and reconstructions for emulating realistic imaging data close to those of actual patients.