Fully 3D Monte Carlo image reconstruction in SPECT using functional regions

Image reconstruction in Single Photon Emission Computed Tomography (SPECT) is affected by physical effects such as photon attenuation, Compton scatter and detector response. These effects can be compensated for by modeling the corresponding spread of photons in 3D within the system matrix used for tomographic reconstruction. The fully 3D Monte Carlo (F3DMC) reconstruction technique consists in calculating this system matrix using Monte Carlo simulations. The inverse problem of tomographic reconstruction is then solved using conventional iterative algorithms such as maximum likelihood expectation maximization (MLEM). Although F3DMC has already shown promising results, its use is currently limited by two major issues: huge size of the fully 3D system matrix and long computation time required for calculating a robust and accurate system matrix. To address these two issues, we propose to calculate the F3DMC system matrix using a spatial sampling matching the functional regions to be reconstructed. In this approach, different regions of interest can be reconstructed with different spatial sampling. For instance, a single value is reconstructed for a functional region assumed to contain uniform activity. To assess the value of this approach, Monte Carlo simulations have been performed using GATE. Results suggest that F3DMC reconstruction using functional regions improves quantitative accuracy compared to the F3DMC reconstruction method proposed so far. In addition, it considerably reduces disk space requirement and duration of the simulations needed to estimate the system matrix. The concept of functional regions might therefore make F3DMC reconstruction practically feasible.Comment: 6 pages, 3 figures, 3rd International Conference on maging Technologies in Biomedical Sciences : ITBS2005, Milos Island, Greece, 25-28 september 2005, submitted to NIM

A multi-resolution image reconstruction method in X-ray computed tomography

International audienceWe propose a multiresolution X-ray imaging method designed for non-destructive testing/ evaluation (NDT/NDE) applications which can also be used for small animal imaging studies. Two sets of projections taken at different magnifications are combined and a multiresolution image is reconstructed. A geometrical relation is introduced in order to combine properly the two sets of data and the processing using wavelet transforms is described. The accuracy of the reconstruction procedure is verified through a comparison to the standard filtered backprojection (FBP) algorithm on simulated data

Reconstruction tri-dimensionnelle complete d'images en spect-ct par modelisation Monte-Carlo

présenté par Z. El Bitar, proceedings sous forme de CDEn tomographie d'émission monophotonique (SPECT), les images 3D normalement reconstruites par des algorithmes de reconstruction analytiques ou itératives bidimensionnelles (2D) pourraient aussi bien être reconstruites avec des algorithmes de reconstruction itérative (3D) qui permettent de compenser les effets physiques perturbant le processus de formation de l'image notamment l'atténuation et la diffusion. Nous avons étudié une technique de reconstruction 3D complète (F3DMC) (Lazaro et al. NIM 2004), dans laquelle le projecteur 3D impliqué dans la reconstruction est estimé par des simulations Monte-Carlo effectuées à partir de données tomodensitométriques du patient

Parallelization of monte Carlo simulations and submission to a grid environment

presentation L. Maign