Jinyi Qi

High-resolution 3D Bayesian image reconstruction using the microPET small-animal scanner

A Bayesian method is described for reconstruction of high-resolution 3D images from the microPET small-animal scanner. Resolution recovery is achieved by explicitly modelling the depth dependent geometric sensitivity for each voxel in combination with an accurate detector response model that includes factors due to photon pair non-collinearity and inter-crystal scatter and penetration. To reduce storage and computational costs we use a factored matrix in which the detector response is modelled using a sinogram blurring kernel. Maximum a posteriori (MAP) images are reconstructed using this model in combination with a Poisson likelihood function and a Gibbs prior on the image. Reconstructions obtained from point source data using the accurate system model demonstrate a potential for near-isotropic FWHM resolution of approximately 1.2 mm at the center of the field of view compared with approximately 2 mm when using an analytic 3D reprojection (3DRP) method with a ramp filter. These results also show the ability of the accurate system model to compensate for resolution loss due to crystal penetration producing nearly constant radial FWHM resolution of 1 mm out to a 4 mm radius. Studies with a point source in a uniform cylinder indicate that as the resolution of the image is reduced to control noise propagation the resolution obtained using the accurate system model is superior to that obtained using 3DRP at matched background noise levels. Additional studies using pie phantoms with hot and cold cylinders of diameter 1-2.5 mm and 18FDG animal studies appear to confirm this observation.

Resolution and noise properties of MAP reconstruction for fully 3-D PET

Derives approximate analytical expressions for the local impulse response and covariance of images reconstructed from fully three-dimensional (3-D) positron emission tomography (PET) data using maximum a posteriori (MAP) estimation. These expressions explicitly account for the spatially variant detector response and sensitivity of a 3-D tomograph. The resulting spatially variant impulse response and covariance are computed using 3-D Fourier transforms. A truncated Gaussian distribution is used to account for the effect on the variance of the nonnegativity constraint used in MAP reconstruction. Using Monte Carlo simulations and phantom data from the microPET small animal scanner, the authors show that the approximations provide reasonably accurate estimates of contrast recovery and covariance of MAP reconstruction for priors with quadratic energy functions. They also describe how these analytical results can be used to achieve near-uniform contrast recovery throughout the reconstructed volume.

Simultaneous in vivo positron emission tomography and magnetic resonance imaging

Positron emission tomography (PET) and magnetic resonance imaging (MRI) are widely used in vivo imaging technologies with both clinical and biomedical research applications. The strengths of MRI include high-resolution, high-contrast morphologic imaging of soft tissues; the ability to image physiologic parameters such as diffusion and changes in oxygenation level resulting from neuronal stimulation; and the measurement of metabolites using chemical shift imaging. PET images the distribution of biologically targeted radiotracers with high sensitivity, but images generally lack anatomic context and are of lower spatial resolution. Integration of these technologies permits the acquisition of temporally correlated data showing the distribution of PET radiotracers and MRI contrast agents or MR-detectable metabolites, with registration to the underlying anatomy. An MRI-compatible PET scanner has been built for biomedical research applications that allows data from both modalities to be acquired simultaneously. Experiments demonstrate no effect of the MRI system on the spatial resolution of the PET system and <10% reduction in the fraction of radioactive decay events detected by the PET scanner inside the MRI. The signal-to-noise ratio and uniformity of the MR images, with the exception of one particular pulse sequence, were little affected by the presence of the PET scanner. In vivo simultaneous PET and MRI studies were performed in mice. Proof-of-principle in vivo MR spectroscopy and functional MRI experiments were also demonstrated with the combined scanner.

Iterative reconstruction techniques in emission computed tomography

In emission tomography statistically based iterative methods can improve image quality relative to analytic image reconstruction through more accurate physical and statistical modelling of high-energy photon production and detection processes. Continued exponential improvements in computing power, coupled with the development of fast algorithms, have made routine use of iterative techniques practical, resulting in their increasing popularity in both clinical and research environments. Here we review recent progress in developing statistically based iterative techniques for emission computed tomography. We describe the different formulations of the emission image reconstruction problem and their properties. We then describe the numerical algorithms that are used for optimizing these functions and illustrate their behaviour using small scale simulations.

Spatiotemporal reconstruction of list-mode PET data

We describe a method for computing a continuous time estimate of tracer density using list-mode positron emission tomography data. The rate function in each voxel is modeled as an inhomogeneous Poisson process whose rate function can be represented using a cubic B-spline basis. The rate functions are estimated by maximizing the likelihood of the arrival times of detected photon pairs over the control vertices of the spline, modified by quadratic spatial and temporal smoothness penalties and a penalty term to enforce nonnegativity. Randoms rate functions are estimated by assuming independence between the spatial and temporal randoms distributions. Similarly, scatter rate functions are estimated by assuming spatiotemporal independence and that the temporal distribution of the scatter is proportional to the temporal distribution of the trues. A quantitative evaluation was performed using simulated data and the method is also demonstrated in a human study using /sup 11/C-raclopride.

A theoretical study of the contrast recovery and variance of MAP reconstructions from PET data

The authors examine the spatial resolution and variance properties of PET images reconstructed using maximum a posteriori (MAP) or penalized-likelihood methods. Resolution is characterized by the contrast recovery coefficient (CRC) of the local impulse response. Simplified approximate expressions are derived for the local impulse response CRCs and variances for each voxel. Using these results the authors propose a practical scheme for selecting spatially variant smoothing parameters to optimize lesion detectability through maximization of the local CRC-to-noise ratio in the reconstructed image.

List-mode maximum-likelihood reconstruction applied to positron emission mammography (PEM) with irregular sampling

Presents a preliminary study of list-mode likelihood reconstruction of images for a rectangular positron emission tomograph (PET) specifically designed to image the human breast. The prospective device consists of small arrays of scintillation crystals for which depth of interaction is estimated. Except in very rare instances, the number of annihilation events detected is expected to be far less than the number of distinguishable events. If one were to histogram the acquired data, most histogram bins would remain vacant. Therefore, it seems natural to investigate the efficacy of processing events one at a time rather than processing the data in histogram format. From a reconstruction perspective, the new tomograph presents a challenge in that the rectangular geometry leads to irregular radial and angular sampling, and the field of view extends completely to the detector faces. Simulations are presented that indicate that the proposed tomograph can detect 8-mm-diameter spherical tumors with a tumor-to-background tracer density ratio of 3:1 using realistic image acquisition parameters. Spherical tumors of 4-mm diameter are near the limit of detectability with the image acquisition parameters used. Expressions are presented to estimate the loss of image contrast due to Compton scattering.

Resolution modeling in PET imaging: Theory, practice, benefits, and pitfalls

In this paper, the authors review the field of resolution modeling in positron emission tomography (PET) image reconstruction, also referred to as point-spread-function modeling. The review includes theoretical analysis of the resolution modeling framework as well as an overview of various approaches in the literature. It also discusses potential advantages gained via this approach, as discussed with reference to various metrics and tasks, including lesion detection observer studies. Furthermore, attention is paid to issues arising from this approach including the pervasive problem of edge artifacts, as well as explanation and potential remedies for this phenomenon. Furthermore, the authors emphasize limitations encountered in the context of quantitative PET imaging, wherein increased intervoxel correlations due to resolution modeling can lead to significant loss of precision (reproducibility) for small regions of interest, which can be a considerable pitfall depending on the task of interest.In this paper, the authors review the field of resolution modeling in positron emission tomography (PET) image reconstruction, also referred to as point-spread-function modeling. The review includes theoretical analysis of the resolution modeling framework as well as an overview of various approaches in the literature. It also discusses potential advantages gained via this approach, as discussed with reference to various metrics and tasks, including lesion detection observer studies. Furthermore, attention is paid to issues arising from this approach including the pervasive problem of edge artifacts, as well as explanation and potential remedies for this phenomenon. Furthermore, the authors emphasize limitations encountered in the context of quantitative PET imaging, wherein increased intervoxel correlations due to resolution modeling can lead to significant loss of precision (reproducibility) for small regions of interest, which can be a considerable pitfall depending on the task of interest.

Fully 3D Bayesian image reconstruction for the ECAT EXACT HR

A fully 3D Bayesian method is described for high resolution reconstruction of images from the Siemens/CTI ECAT EXACT HR+ whole body positron emission tomography (PET) scanner. To maximize resolution recovery from the system the authors model depth dependent geometric efficiency, intrinsic detector efficiency, photon pair non-colinearity, crystal penetration and inter-crystal scatter. They also explicitly model the effects of axial rebinning and angular mashing on the detection probability or system matrix. By fully exploiting sinogram symmetries and using a factored system matrix and automated indexing schemes, the authors are able to achieve substantial savings in both the storage size and time required to compute forward and backward projections. Reconstruction times are further reduced using multi-threaded programming on a four processor Unix server. Bayesian reconstructions are computed using a Huber prior and a shifted-Poisson likelihood model that accounts for the effects of randoms subtraction and scatter. Reconstructions of phantom data show that the 3D Bayesian method can achieve improved FWHM resolution and contrast recovery ratios at matched background noise levels compared to both the 3D reprojection method and an OSEM method based on the shifted-Poisson model.

Approximate maximum likelihood hyperparameter estimation for Gibbs priors

The parameters of the prior, the hyperparameters, play an important role in Bayesian image estimation. Of particular importance for the case of Gibbs priors is the global hyperparameter, beta, which multiplies the Hamiltonian. Here we consider maximum likelihood (ML) estimation of beta from incomplete data, i.e., problems in which the image, which is drawn from a Gibbs prior, is observed indirectly through some degradation or blurring process. Important applications include image restoration and image reconstruction from projections. Exact ML estimation of beta from incomplete data is intractable for most image processing. Here we present an approximate ML estimator that is computed simultaneously with a maximum a posteriori (MAP) image estimate. The algorithm is based on a mean field approximation technique through which multidimensional Gibbs distributions are approximated by a separable function equal to a product of one-dimensional (1-D) densities. We show how this approach can be used to simplify the ML estimation problem. We also show how the Gibbs-Bogoliubov-Feynman (GBF) bound can be used to optimize the approximation for a restricted class of problems. We present the results of a Monte Carlo study that examines the bias and variance of this estimator when applied to image restoration.