Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userApplication of Monte Carlo algorithms to cardiac imaging reconstruction
116
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
Monte Carlo algorithms have a growing impact on nuclear medicine reconstruction processes. One of the main limitations of myocardial perfusion imaging (MPI) is the effective mitigation of the scattering component, which is particularly challenging in Single Photon Emission Computed Tomography (SPECT). In SPECT, no timing information can be retrieved to locate the primary source photons. Monte Carlo methods allow an event-by-event simulation of the scattering kinematics, which can be incorporated into a model of the imaging system response. This approach was adopted since the late Nineties by several authors, and recently took advantage of the increased computational power made available by high-performance CPUs and GPUs. These recent developments enable a fast image reconstruction with an improved image quality, compared to deterministic approaches. Deterministic approaches are based on energy-windowing of the detector response, and on the cumulative estimate and subtraction of the scattering component. In this paper, we review the main strategies and algorithms to correct for the scattering effect in SPECT and focus on Monte Carlo developments, which nowadays allow the three-dimensional reconstruction of SPECT cardiac images in a few seconds.
rate research
Read More
We present a novel approach aimed at high-performance uncertainty quantification for time-dependent problems governed by partial differential equations. In particular, we consider input uncertainties described by a Karhunen-Loeeve expansion and compute statistics of high-dimensional quantities-of-interest, such as the cardiac activation potential. Our methodology relies on a close integration of multilevel Monte Carlo methods, parallel iterative solvers, and a space-time discretization. This combination allows for space-time adaptivity, time-changing domains, and to take advantage of past samples to initialize the space-time solution. The resulting sequence of problems is distributed using a multilevel parallelization strategy, allocating batches of samples having different sizes to a different number of processors. We assess the performance of the proposed framework by showing in detail its application to the solution of nonlinear equations arising from cardiac electrophysiology. Specifically, we study the effect of spatially-correlated perturbations of the heart fibers conductivities on the mean and variance of the resulting activation map. As shown by the experiments, the theoretical rates of convergence of multilevel Monte Carlo are achieved. Moreover, the total computational work for a prescribed accuracy is reduced by an order of magnitude with respect to standard Monte Carlo methods.
In Single Photon Emission Computed Tomography (SPECT), 3D images usually reconstructed by performing a set of bidimensional (2D) analytical or iterative reconstructions can also be reconstructed using an iterative reconstruction algorithm involving a 3D projector. Accurate Monte Carlo (MC) simulations modeling all the physical effects that affect the imaging process can be used to estimate this projector. However, the accuracy of the projector is affected by the stochastic nature of MC simulations. In this paper, we study the accuracy of the reconstructed images with respect to the number of simulated histories used to estimate the MC projector. Furthermore, we study the impact of applying the bootstrapping technique when estimating the projector
INTPIX4NA is an integration-type silicon-on-insulator pixel detector. This detector has a 14.1 x 8.7 mm^2 sensitive area, 425,984 (832 column x 512 row matrix) pixels and the pixel size is 17 x 17 um^2. This detector was developed for residual stress measurement using X-rays (the cos alpha method). The performance of INTPIX4NA was tested with the synchrotron beamlines of the Photon Factory (KEK), and the following results were obtained. The modulation transfer function, the index of the spatial resolution, was more than 50% at the Nyquist frequency (29.4 cycle/mm). The energy resolution analyzed from the collected charge counts is 35.3%--46.2% at 5.415 keV, 21.7%--35.6% at 8 keV, and 15.7%--19.4% at 12 keV. The X-ray signal can be separated from the noise even at a low energy of 5.415 keV at room temperature (approximately 25--27 degree Celsius). The maximum frame rate at which the signal quality can be maintained is 153 fps in the current measurement system. These results satisfy the required performance in the air and at room temperature (approximately 25--27 degree Celsius) condition that is assumed for the environment of the residual stress measurement.
Monte Carlo simulations are increasingly used in scintigraphic imaging to model imaging systems and to develop and assess tomographic reconstruction algorithms and correction methods for improved image quantitation. GATE (GEANT 4 Application for Tomographic Emission) is a new Monte Carlo simulation platform based on GEANT4 dedicated to nuclear imaging applications. This paper describes the GATE simulation of a prototype of scintillation camera dedicated to small animal imaging and consisting of a CsI(Tl) crystal array coupled to a position sensitive photomultiplier tube. The relevance of GATE to model the camera prototype was assessed by comparing simulated 99mTc point spread functions, energy spectra, sensitivities, scatter fractions and image of a capillary phantom with the corresponding experimental measurements. Results showed an excellent agreement between simulated and experimental data: experimental spatial resolutions were predicted with an error less than 100 mu m. The difference between experimental and simulated system sensitivities for different source-to-collimator distances was within 2%. Simulated and experimental scatter fractions in a [98-182 keV] energy window differed by less than 2% for sources located in water. Simulated and experimental energy spectra agreed very well between 40 and 180 keV. These results demonstrate the ability and flexibility of GATE for simulating original detector designs. The main weakness of GATE concerns the long computation time it requires: this issue is currently under investigation by the GEANT4 and the GATE collaborations
We describe the Monte Carlo (MC) simulation package of the Borexino detector and discuss the agreement of its output with data. The Borexino MC ab initio simulates the energy loss of particles in all detector components and generates the resulting scintillation photons and their propagation within the liquid scintillator volume. The simulation accounts for absorption, reemission, and scattering of the optical photons and tracks them until they either are absorbed or reach the photocathode of one of the photomultiplier tubes. Photon detection is followed by a comprehensive simulation of the readout electronics response. The algorithm proceeds with a detailed simulation of the electronics chain. The MC is tuned using data collected with radioactive calibration sources deployed inside and around the scintillator volume. The simulation reproduces the energy response of the detector, its uniformity within the fiducial scintillator volume relevant to neutrino physics, and the time distribution of detected photons to better than 1% between 100 keV and several MeV. The techniques developed to simulate the Borexino detector and their level of refinement are of possible interest to the neutrino community, especially for current and future large-volume liquid scintillator experiments such as Kamland-Zen, SNO+, and Juno.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
