Practical considerations for image-based PSF and blobs reconstruction in PET.

Published on May 16, 2013in Physics in Medicine and Biology3.609
路 DOI :10.1088/0031-9155/58/11/3849
Simon Stute1
Estimated H-index: 1
(CEA: French Alternative Energies and Atomic Energy Commission),
Claude Comtat27
Estimated H-index: 27
Sources
Abstract
Iterative reconstructions in positron emission tomography (PET) need a model relating the recorded data to the object/patient being imaged, called the system matrix (SM). The more realistic this model, the better the spatial resolution in the reconstructed images. However, a serious concern when using a SM that accurately models the resolution properties of the PET system is the undesirable edge artefact, visible through oscillations near sharp discontinuities in the reconstructed images. This artefact is a natural consequence of solving an ill-conditioned inverse problem, where the recorded data are band-limited. In this paper, we focus on practical aspects when considering image-based point-spread function (PSF) reconstructions. To remove the edge artefact, we propose to use a particular case of the method of sieves (Grenander 1981 Abstract Inference New York: Wiley), which simply consists in performing a standard PSF reconstruction, followed by a post-smoothing using the PSF as the convolution kernel. Using analytical simulations, we investigate the impact of different reconstruction and PSF modelling parameters on the edge artefact and its suppression, in the case of noise-free data and an exactly known PSF. Using Monte-Carlo simulations, we assess the proposed method of sieves with respect to the choice of the geometric projector and the PSF model used in the reconstruction. When the PSF model is accurately known, we show that the proposed method of sieves succeeds in completely suppressing the edge artefact, though after a number of iterations higher than typically used in practice. When applying the method to realistic data (i.e. unknown true SM and noisy data), we show that the choice of the geometric projector and the PSF model does not impact the results in terms of noise and contrast recovery, as long as the PSF has a width close to the true PSF one. Equivalent results were obtained using either blobs or voxels in the same conditions (i.e. the blob鈥檚 density function being the same as the voxel-based PSF). From a practical point-of-view, the method can be used to perform fast reconstructions based on very simple models (compared to sinogram-based PSF modelling), producing artefact-free images with a better compromise between noise and spatial resolution than images reconstructed without or with under-estimated PSF. Besides, the method inherently limits the spatial resolution in the reconstructed images to the intrinsic one of the PET system.
馃摉 Papers frequently viewed together
References44
Newest
The response of a PET system can be described by its characteristic Point Spread Function (PSF) representing the spatial degradation of a point source due to physical effects and system design. If the PSF is accounted for in the reconstruction algorithm, better image quality and spatial resolution may be achieved. Un- fortunately, a common behaviour of unregularized iterative reconstruction techniques is represented by the increase of noise as the iterations proceed, while-on the other hand-a hi...
Source
#1Simon StuteH-Index: 12
#2Johan Nuyts (Katholieke Universiteit Leuven)H-Index: 62
Last. Claude ComtatH-Index: 27
view all 7 authors...
The use of iterative image reconstruction algorithms with resolution modeling allows for reduced partial volume effect without noise increase. However, it is now recognized that EM-ML type algorithms are biased in very low counts images, in particular for cold regions. Alternative ML algorithms that allow for negative image voxels have been proposed to reduce the bias: NEG-ML of Nuyts et al and AB-ML of Byrne.
Source
#1Audrey PepinH-Index: 1
#2Simon StuteH-Index: 12
Last. Claude ComtatH-Index: 27
view all 4 authors...
The reconstruction of PET data without proper correction for non-uniform detection efficiency leads to severe artifacts in the image. This is also true for simulated data, particularly for Monte Carlo simulations that aim to realistically model the detectors. In the present study, we showed that the simulation with GATE software of the ECAT HR+, a block detector ring scanner, results in non negligible variations across the lines of response in the detection sensitivity. Like for real scanners, t...
Source
#1Johan Nuyts (Katholieke Universiteit Leuven)H-Index: 62
#2Simon StuteH-Index: 12
Last. Claude ComtatH-Index: 27
view all 6 authors...
The AB-ML algorithm [2] and the proposed NEG-ML algorithm can reduce or eliminate the bias in reconstructions from noisy data. AB-ML shifts the Poisson distribution, while NEG-ML modifies it near zero. Both algorithms have a single parameter that defines the balance between MLEM-like and unweighted-least-squares-like behaviour.
Source
#1Shan Tong (UW: University of Washington)H-Index: 7
#2Adam M. Alessio (UW: University of Washington)H-Index: 33
Last. Paul E. Kinahan (UW: University of Washington)H-Index: 70
view all 6 authors...
PSF (point spread function) based image reconstruction causes an overshoot at sharp intensity transitions (edges) of the object. This edge artifact, or ringing, has not been fully studied. In this work, we analyze the properties of edge artifacts in PSF-based reconstruction in an effort to develop mitigation methods. Our study is based on 1D and 2D simulation experiments. Two approaches are adopted to analyze the artifacts. In the system theory approach, we relate the presence of edge artifacts ...
Source
This paper proposes a 2-step image reconstruction method in which the nonnegativity constraint in the iterative maximum-likelihood expectation maximization (MLEM) algorithm is used to effectively reduce Gibbs ringing artifacts.Gibbs artifacts are difficult to control during imaging reconstruction. The proposed method uses the postprocessing strategy to suppress Gibbs artifacts. In the first step, a raw image is reconstructed from projections without correction for point spread function (PSF). Th...
Source
#1S茅bastien Jan (CEA: French Alternative Energies and Atomic Energy Commission)H-Index: 18
#2D. BenoitH-Index: 3
Last. Ir猫ne BuvatH-Index: 47
view all 20 authors...
GATE (Geant4 Application for Emission Tomography) is a Monte Carlo simulation platform developed by the OpenGATE collaboration since 2001 and first publicly released in 2004. Dedicated to the modelling of planar scintigraphy, single photon emission computed tomography (SPECT) and positron emission tomography (PET) acquisitions, this platform is widely used to assist PET and SPECT research. A recent extension of this platform, released by the OpenGATE collaboration as GATE V6, now also enables mo...
Source
#1Simon Stute (CNRS: Centre national de la recherche scientifique)H-Index: 12
#2D Benoit (CNRS: Centre national de la recherche scientifique)H-Index: 1
Last. Ir猫ne Buvat (CNRS: Centre national de la recherche scientifique)H-Index: 47
view all 5 authors...
Positron emission tomography (PET) images suffer from low spatial resolution and signal-to-noise ratio. Accurate modelling of the effects affecting resolution within iterative reconstruction algorithms can improve the trade-off between spatial resolution and signal-to-noise ratio in PET images. In this work, we present an original approach for modelling the resolution loss introduced by physical interactions between and within the crystals of the tomograph and we investigate the impact of such m...
Source
#1Sangtae Ahn (GE: General Electric)H-Index: 23
#2Evren Asma (GE: General Electric)H-Index: 16
Last. Charles W. Stearns (GE Healthcare)H-Index: 21
view all 6 authors...
Accurate system modeling is essential for improved quantitation and lesion detection. Many investigators have made efforts to accurately model detector blurring using point spread functions (PSFs) in sinogram space and to incorporate them into image reconstruction for accurate quantitation. It has been observed that incorporating detector PSF into reconstruction leads to improved contrast recovery and resolution with reduced noise but introduces edge artifacts. It is not straightforward to inves...
Source
#1Pablo Aguiar (University of Santiago de Compostela)H-Index: 19
#2Magdalena Rafecas (CSIC: Spanish National Research Council)H-Index: 22
Last. Dom猫nec Ros (University of Barcelona)H-Index: 19
view all 7 authors...
Purpose: In the present work, the authors compare geometrical and Monte Carlo projectors in detail. The geometrical projectors considered were the conventional geometrical Siddon ray-tracer S-RT and the orthogonal distance-based ray-tracer OD-RT, based on computing the orthogonal distance from the center of image voxel to the line-of-response. A comparison of these geometrical projectors was performed using different point spread function PSF models. The Monte Carlobased method under considerati...
Source
Cited By30
Newest
In this work we continue studies of the uncertainty quantification problem in emission tomographies such as PET or SPECT. In particular, we consider a scenario when additional multimodal data (e.g., anatomical MRI images) are available. To solve the aforementioned problem we adapt the recently proposed nonparametric posterior learning technique to the context of Poisson-type data in emission tomography. Using this approach we derive sampling algorithms which are trivially parallelizable, scalabl...
Noise suppression is particularly important in low count positron emission tomography (PET) imaging. Post-smoothing (PS) and regularization methods which aim to reduce noise also tend to reduce resolution and introduce bias. Alternatively, anatomical information from another modality such as magnetic resonance (MR) imaging can be used to improve image quality. Convolutional neural networks (CNNs) are particularly well suited to such joint image processing, but usually require large amounts of tr...
Source
#1Marina Filipovi膰 (University of Paris-Sud)H-Index: 1
#2Claude Comtat (University of Paris-Sud)H-Index: 27
Last. Simon Stute (University of Nantes)H-Index: 1
view all 3 authors...
The time-of-flight (TOF) feature of PET scanners has been used for a long time in PET reconstruction, but many implementational aspects are still incomplete or ambiguous in the literature. Here we formalize and present theoretical and practical implementation details for the reconstruction of clinical TOF histogram and list-mode data using ML-EM. Relevant aspects include the computation of the TOF component of the system matrix, the processing of TOF bins, the use of estimations of random and sc...
Source
#1James Bland (University of Cambridge)H-Index: 3
#2Abolfazl Mehranian (University of Cambridge)H-Index: 1
Last. Andrew J. Reader (University of Cambridge)H-Index: 32
view all 8 authors...
PURPOSE: Numerous image reconstruction methodologies for positron emission tomography (PET) have been developed that incorporate magnetic resonance (MR) imaging structural information, producing reconstructed images with improved suppression of noise and reduced partial volume effects. However, the influence of MR structural information also increases the possibility of suppression or bias of structures present only in the PET data (PET-unique regions). To address this, further developments for ...
Source
#1Yu-Fen ChangH-Index: 7
#2Kyle ThompsonH-Index: 1
Last. Zisis PapandreouH-Index: 6
view all 5 authors...
Positron emission tomography (PET) is a powerful tool for mapping the dynamic distribution of substances in visually opaque systems. A modular, scalable, transformable, and relocatable PET system is custom built to provide optimized flexible PET imaging for a wider variety of applications. In coordination with the purposely designed hardware, we develop adaptive data processing and image reconstruction pipelines. We perform Monte Carlo simulation and demonstrate the application of the developed ...
Source
#1Yanna Cruz Cavalcanti (University of Toulouse)H-Index: 3
#2Thomas Oberlin (University of Toulouse)H-Index: 12
Last. Clovis Tauber (Fran莽ois Rabelais University)H-Index: 16
view all 7 authors...
Factor analysis has proven to be a relevant tool for extracting tissue time-activity curves (TACs) in dynamic PET images, since it allows for an unsupervised analysis of the data. Reliable and interpretable results are possible only if it is considered with respect to suitable noise statistics. However, the noise in reconstructed dynamic PET images is very difficult to characterize, despite the Poissonian nature of the count rates. Rather than explicitly modeling the noise distribution, this pap...
Source
#1Yue-Meng Feng (INSA Lyon: Institut national des sciences Appliqu茅es de Lyon)H-Index: 2
#2Ane Etxebeste (INSA Lyon: Institut national des sciences Appliqu茅es de Lyon)H-Index: 1
Last. Voichita Maxim (INSA Lyon: Institut national des sciences Appliqu茅es de Lyon)H-Index: 1
view all 5 authors...
The Compton camera is a gamma ray imaging device already employed in astronomy and still in investigation for clinical domain. A key point in the imaging process is the tomographic reconstruction step. When the acquisition parameters and the a priori information are correctly accounted for, iterative algorithms are able to produce accurate images by compensating for measurement uncertainties and statistical noise. In this work we focus on the list-mode maximum likelihood expectation maximization...
Source
#1Mael Millardet (French Institute of Health and Medical Research)H-Index: 1
#2Sa茂d Moussaoui (CNRS: Centre national de la recherche scientifique)H-Index: 18
Last. Thomas Carlier (French Institute of Health and Medical Research)H-Index: 12
view all 6 authors...
Y TTRIUM-90 ( 90 Y) is an element used in nuclear medicine for the radioembolization of hepatic tumours. 90 Y is a 尾 鈭 emitter whose decay also produces positrons with a very low branching ratio of 3.2 脳 10 鈭 5 . The possibility of using the emitted positrons for a post-treatment PET imaging has been established in 2004 [1] . PET imaging using 90 Y has then gained significant interest [2] , [3] , but remains challenging because the classical reconstruction algorithms, such as OSEM 1 , introduce ...
Source
#1Yu-Fen ChangH-Index: 7
#2Steven D. SicilianoH-Index: 61
Last. Aram TeymurazyanH-Index: 2
view all 11 authors...
Source
#1James Bland (St Thomas' Hospital)H-Index: 3
#2Martin A. Belzunce (St Thomas' Hospital)H-Index: 7
Last. Andrew J. Reader (St Thomas' Hospital)H-Index: 32
view all 6 authors...
Positron emission tomography (PET) is a highly sensitive functional and molecular imaging modality which can measure picomolar concentrations of an injected radionuclide. However, the physical sensitivity of PET is limited, and reducing the injected dose leads to low count data and noisy reconstructed images. A highly effective way of reducing noise is to reparameterize the reconstruction in terms of MR-derived spatial basis functions. Spatial basis functions derived using the kernel method have...
Source
This website uses cookies.
We use cookies to improve your online experience. By continuing to use our website we assume you agree to the placement of these cookies.
To learn more, you can find in our Privacy Policy.