Results 1  10
of
100
Practical conebeam algorithm
 J Opt Soc Am
, 1984
"... A convolutionbackprojection formula is deduced for direct reconstruction of a threedimensional density function from a set of twodimensional projections. The formula is approximate but has useful properties, including errors that are relatively small in many practical instances and a form that le ..."
Abstract

Cited by 153 (0 self)
 Add to MetaCart
A convolutionbackprojection formula is deduced for direct reconstruction of a threedimensional density function from a set of twodimensional projections. The formula is approximate but has useful properties, including errors that are relatively small in many practical instances and a form that leads to convenient computation. It reduces to the standard fanbeam formula in the plane that is perpendicular to the axis of rotation and contains the point source. The algorithm is applied to a mathematical phantom as an example of its performance. 1.
A local update strategy for iterative reconstruction from projections
 IEEE Tr. Sig. Proc
, 1993
"... Iterative methods for statisticallybased reconstruction from projections are computationally costly relative to convolution backprojection, but allow useful image reconstruction from sparse and noisy data. We present a method for Bayesian reconstruction which relies on updates of single pixel value ..."
Abstract

Cited by 121 (31 self)
 Add to MetaCart
Iterative methods for statisticallybased reconstruction from projections are computationally costly relative to convolution backprojection, but allow useful image reconstruction from sparse and noisy data. We present a method for Bayesian reconstruction which relies on updates of single pixel values, rather than the entire image, at each iteration. The technique is similar to GaussSeidel (GS) iteration for the solution of differential equations on finite grids. The computational cost per iteration of the GS approach is found to be approximately equal to that of gradient methods. For continuously valued images, GS is found to have significantly better convergence at modes representing high spatial frequencies. In addition, GS is well suited to segmentation when the image is constrained to be discretely valued. We demonstrate that Bayesian segmentation using GS iteration produces useful estimates at much lower signaltonoise ratios than required for continuously valued reconstruction. This paper includes analysis of the convergence properties of gradient ascent and GS for reconstruction from integral projections, and simulations of both maximumlikelihood and maximum a posteriori cases.
Nonuniform Fast Fourier Transforms Using MinMax Interpolation
 IEEE Trans. Signal Process
, 2003
"... The FFT is used widely in signal processing for efficient computation of the Fourier transform (FT) of finitelength signals over a set of uniformlyspaced frequency locations. However, in many applications, one requires nonuniform sampling in the frequency domain, i.e.,a nonuniform FT . Several pap ..."
Abstract

Cited by 88 (15 self)
 Add to MetaCart
The FFT is used widely in signal processing for efficient computation of the Fourier transform (FT) of finitelength signals over a set of uniformlyspaced frequency locations. However, in many applications, one requires nonuniform sampling in the frequency domain, i.e.,a nonuniform FT . Several papers have described fast approximations for the nonuniform FT based on interpolating an oversampled FFT. This paper presents an interpolation method for the nonuniform FT that is optimal in the minmax sense of minimizing the worstcase approximation error over all signals of unit norm. The proposed method easily generalizes to multidimensional signals. Numerical results show that the minmax approach provides substantially lower approximation errors than conventional interpolation methods. The minmax criterion is also useful for optimizing the parameters of interpolation kernels such as the KaiserBessel function.
Reconstruction from projections based on detection and estimation of objectsParts I and 11: Performance analysis and robustness analysis
 IEEE Trans. Acoust., Speech, Signal Processing
, 1984
"... AbstructThe problem of reconstructing a multidimensional field from noisy, limited projection measurements is approached using an objectbased stochastic field model. Objects within a cross section are characterized by a fiitedimensional set of parameters, which are estimated directly from limited ..."
Abstract

Cited by 54 (11 self)
 Add to MetaCart
AbstructThe problem of reconstructing a multidimensional field from noisy, limited projection measurements is approached using an objectbased stochastic field model. Objects within a cross section are characterized by a fiitedimensional set of parameters, which are estimated directly from limited, noisy projection measurements using maximum likelihood estimation. In Part I, the computational structure and performance of the ML estimation procedure are investigated for the problem of locating a single object in a deterministic background; simulations are also presented. In Part 11, the issue of robustness to modeling errors is addressed. PART I PERFORMANCE ANALYSIS
Roxels: Responsibility weighted 3d volume reconstruction
 In Proc Int Conf on Computer Vision
, 1999
"... This paper examines the problem of reconstructing a voxelized representation of 3D space from a series of images. An iterative algorithm is used to find the scene model which jointly explains all the observed images by determining which region of space is responsible for each of the observations. Th ..."
Abstract

Cited by 54 (2 self)
 Add to MetaCart
This paper examines the problem of reconstructing a voxelized representation of 3D space from a series of images. An iterative algorithm is used to find the scene model which jointly explains all the observed images by determining which region of space is responsible for each of the observations. The current approach formulates the problem as one of optimization over estimates of these responsibilities. The process converges to a distribution of responsibility which accurately reflects the constraints provided by the observations, the positions and shape of both solid and transparent objects, and the uncertainty which remains. Reconstruction is robust, and gracefully represents regions of space in which there is little certainty about the exact structure due to limited, nonexistent, or contradicting data. Rendered images of voxel spaces recovered from synthetic and real observation images are shown. 1
Rapid 3D ConeBeam Reconstruction with the Simultaneous Algebraic Reconstruction Technique (SART) Using 2D Texture Mapping Hardware
, 2000
"... Algebraic reconstruction methods, such as the algebraic reconstruction technique (ART) and the related simultaneous ART (SART), reconstruct a twodimensional (2D) or threedimensional (3D) object from its Xray projections. The algebraic methods have, in certain scenarios, many advantages over t ..."
Abstract

Cited by 33 (9 self)
 Add to MetaCart
Algebraic reconstruction methods, such as the algebraic reconstruction technique (ART) and the related simultaneous ART (SART), reconstruct a twodimensional (2D) or threedimensional (3D) object from its Xray projections. The algebraic methods have, in certain scenarios, many advantages over the more popular Filtered Backprojection approaches and have also recently been shown to perform well for 3D conebeam reconstruction. However, so far the slow speed of these iterative methods have prohibited their routine use in clinical applications. In this paper, we address this shortcoming and investigate the utility of widely available 2D texture mapping graphics hardware for the purpose of accelerating the 3D algebraic reconstruction. We find that this hardware allows 3D conebeam reconstructions to be obtained at almost interactive speeds, with speedups of over 50 with respect to implementations that only use generalpurpose CPUs. However, we also find that the reconstruction quality is rather sensitive to the resolution of the framebuffer, and to address this critical issue we propose a scheme that extends the precision of a given framebuffer by 4 bits, using the color channels. With this extension, a 12bit framebuffer delivers useful reconstructions for 0.5% tissue contrast, while an 8bit framebuffer requires 4%. Since graphics hardware generates an entire image for each volume projection, it is most appropriately used with an algebraic reconstruction method that performs volume correction at that granularity as well, such as SART or SIRT. We chose SART for its faster convergence properties. Index TermsAlgebraic reconstruction technique (ART), computed tomography, cone beam reconstruction, hardware acceleration, simultaneous algebraic reconstruction techni...
The Ordered Subsets Mirror Descent Optimization Method with Applications to Tomography
 SIAM J. Optim
, 2001
"... Abstract. We describe an optimization problem arising in reconstructing 3D medical images from Positron Emission Tomography (PET). A mathematical model of the problem, based on the Maximum Likelihood principle is posed as a problem of minimizing a convex function of several millions variables over t ..."
Abstract

Cited by 32 (5 self)
 Add to MetaCart
Abstract. We describe an optimization problem arising in reconstructing 3D medical images from Positron Emission Tomography (PET). A mathematical model of the problem, based on the Maximum Likelihood principle is posed as a problem of minimizing a convex function of several millions variables over the standard simplex. To solve a problem of these characteristics, we develop and implement a new algorithm, Ordered Subsets Mirror Descent, and demonstrate, theoretically and computationally, that it is well suited for solving the PET reconstruction problem. Key words: positron emission tomography, maximum likelihood, image reconstruction, convex optimization, mirror descent. 1
Iterative tomographic image reconstruction using Fourierbased forward and back projectors
 IEEE Trans. Med. Imag
, 2004
"... Fourierbased reprojection methods have the potential to reduce the computation time in iterative tomographic image reconstruction. Interpolation errors are a limitation of Fourierbased reprojection methods. We apply a minmax interpolation method for the nonuniform fast Fourier transform (NUFFT) t ..."
Abstract

Cited by 24 (4 self)
 Add to MetaCart
Fourierbased reprojection methods have the potential to reduce the computation time in iterative tomographic image reconstruction. Interpolation errors are a limitation of Fourierbased reprojection methods. We apply a minmax interpolation method for the nonuniform fast Fourier transform (NUFFT) to minimize the interpolation errors. Numerical results show that the minmax NUFFT approach provides substantially lower approximation errors in tomographic reprojection and backprojection than conventional interpolation methods.
Wavelet Methods For The Inversion Of Certain Homogeneous Linear Operators In The Presence Of Noisy Data
, 1994
"... In this dissertation we explore the use of wavelets in certain linear inverse problems with discrete, noisy data. We observe discrete samples of a process y(u) = (Kf)(u)+ z(u), where K is a linear operator, z is a noise process, and f is a function we wish to recover from the data. In the problems ..."
Abstract

Cited by 22 (1 self)
 Add to MetaCart
In this dissertation we explore the use of wavelets in certain linear inverse problems with discrete, noisy data. We observe discrete samples of a process y(u) = (Kf)(u)+ z(u), where K is a linear operator, z is a noise process, and f is a function we wish to recover from the data. In the problems that we consider, the inverse of K, K \Gamma1 , either does not exist or is poorly behaved. Such problems are termed illposed i.e., ones in which small changes in the data may lead to large changes in the recovered version of f . Our methods are most effective for problems where the operator K is homogeneous with respect to dilations, such as integration, fractional integration, convolution, and the Radon transform. The theoretical framework in which we work is that of Donoho's (1992) WaveletVaguelette Decomposition (WVD). The WVD uses wavelets and vaguelettes (almost wavelets) to decompose the operator K. Although this formally resembles the Singular Value Decomposition (SVD), the use o...
Overview of methods for image reconstruction from projections in emission computed tomography
 PROC. IEEE
, 2003
"... Emission computed tomography (ECT) is a technology for medical imaging whose importance is increasing rapidly. There is a growing appreciation for the value of the functional (as opposed to anatomical) information that is provided by ECT and there are significant advancements taking place, both in t ..."
Abstract

Cited by 18 (1 self)
 Add to MetaCart
Emission computed tomography (ECT) is a technology for medical imaging whose importance is increasing rapidly. There is a growing appreciation for the value of the functional (as opposed to anatomical) information that is provided by ECT and there are significant advancements taking place, both in the instrumentation for data collection, and in the computer methods for generating images from the measured data. These computer methods are designed to solve the inverse problem known as “image reconstruction from projections.” This paper uses the various models of the data collection process as the framework for presenting an overview of the wide variety of methods that have been developed for image reconstruction in the major subfields of ECT, which are positron emission tomography (PET) and singlephoton emission computed tomography (SPECT). The overall sequence of the major sections in the paper, and the presentation within each major section, both proceed from the more realistic and general models to those that are idealized and application specific. For most of the topics, the description proceeds from the threedimensional case to the twodimensional case. The paper presents a broad overview of algorithms for PET and SPECT, giving references to the literature where these algorithms and their applications are described in more detail.