Results 1  10
of
28
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 83 (13 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.
Toeplitzbased iterative image reconstruction for MRI with correction for magnetic field inhomogeneity
 IEEE Trans. Signal Process
, 2005
"... Abstract — In some types of magnetic resonance (MR) imaging, particularly functional brain scans, the conventional Fourier model for the measurements is inaccurate. Magnetic field inhomogeneities, caused by imperfect main fields and by magnetic susceptibility variations, induce distortions in images ..."
Abstract

Cited by 18 (3 self)
 Add to MetaCart
Abstract — In some types of magnetic resonance (MR) imaging, particularly functional brain scans, the conventional Fourier model for the measurements is inaccurate. Magnetic field inhomogeneities, caused by imperfect main fields and by magnetic susceptibility variations, induce distortions in images that are reconstructed by conventional Fourier methods. These artifacts hamper the use of functional MR imaging (fMRI) in brain regions near air/tissue interfaces. Recently, iterative methods that combine the conjugate gradient (CG) algorithm with nonuniform FFT (NUFFT) operations have been shown to provide considerably improved image quality relative to the conjugatephase method. However, for nonCartesian kspace trajectories, each CGNUFFT iteration requires numerous kspace interpolations, operations that are computationally expensive and poorly suited to fast hardware implementations. This paper proposes a faster iterative approach to fieldcorrected MR image reconstruction based on the CG algorithm and certain Toeplitz matrices. This CGToeplitz approach requires kspace interpolations only for the initial iteration; thereafter only FFTs are required. Simulation results show that the proposed CGToeplitz approach produces equivalent image quality as the CGNUFFT method with significantly reduced computation time. Index Terms — fMRI imaging, spiral trajectory, magnetic susceptibility, nonCartesian sampling I.
How GPUs can improve the quality of magnetic resonance imaging
 In The First Workshop on General Purpose Processing on Graphics Processing Units
, 2007
"... Abstract — In magnetic resonance imaging (MRI), nonCartesian scan trajectories are advantageous in a wide variety of emerging applications. Advanced reconstruction algorithms that operate directly on nonCartesian scan data using optimality criteria such as leastsquares (LS) can produce significan ..."
Abstract

Cited by 9 (3 self)
 Add to MetaCart
Abstract — In magnetic resonance imaging (MRI), nonCartesian scan trajectories are advantageous in a wide variety of emerging applications. Advanced reconstruction algorithms that operate directly on nonCartesian scan data using optimality criteria such as leastsquares (LS) can produce significantly better images than conventional algorithms that apply a fast Fourier transform (FFT) after interpolating the scan data onto a Cartesian grid. However, advanced LS reconstructions require significantly more computation than conventional reconstructions based on the FFT. For example, one LS algorithm requires nearly six hours to reconstruct a single threedimensional image on a modern CPU. Our work demonstrates that this advanced reconstruction can be performed quickly and efficiently on a modern GPU, with the reconstruction of a 64 3 3D image requiring just three minutes, an acceptable latency for key applications. This paper describes how the reconstruction algorithm leverages the resources of the GeForce 8800 GTX (G80) to achieve over 150 GFLOPS in performance. We find that the combination of tiling the data and storing the data in the G80’s constant memory dramatically reduces the algorithm’s required bandwidth to offchip memory. The G80’s special functional units provide substantial acceleration for the trigonometric computations in the algorithm’s inner loops. Finally, experimentdriven code transformations increase the reconstruction’s performance by as much as 60 % to 80%. I.
Iterative image reconstruction in MRI with separate magnitude and phase regularization
 In Proc. IEEE Intl. Symp. Biomed. Imag
, 2004
"... Iterative methods for image reconstruction in MRI are useful in several applications, including reconstruction from nonCartesian kspace samples, compensation for magnetic field inhomogeneities, and imaging with multiple receive coils. Existing iterative MR image reconstruction methods are either u ..."
Abstract

Cited by 9 (2 self)
 Add to MetaCart
Iterative methods for image reconstruction in MRI are useful in several applications, including reconstruction from nonCartesian kspace samples, compensation for magnetic field inhomogeneities, and imaging with multiple receive coils. Existing iterative MR image reconstruction methods are either unregularized, and therefore sensitive to noise, or have used regularization methods that smooth the complex valued image. These existing methods regularize the real and imaginary components of the image equally. In many MRI applications, including T ∗ 2weighted imaging as used in fMRI BOLD imaging, one expects most of the signal information of interest to be contained in the magnitude of the voxel value, whereas the phase values are expected to vary smoothly spatially. This paper proposes separate regularization of the magnitude and phase components, preserving the spatial resolution of the magnitude component while strongly regularizing the phase component. This leads to a nonconvex regularized leastsquares cost function. We describe a new iterative algorithm that monotonically decreases this cost function. The resulting images have reduced noise relative to conventional regularization methods. 1.
NoiseAdaptive Nonlinear Diffusion Filtering of MR Images With Spatially Varying Noise Levels
 Magnetic Resonance in Medicine
, 2004
"... Anisotropic diffusion filtering is widely used for MR image enhancement. However, the anisotropic filter is nonoptimal for MR images with spatially varying noise levels, such as images reconstructed from sensitivityencoded data and intensity inhomogeneitycorrected images. In this work, a new metho ..."
Abstract

Cited by 8 (0 self)
 Add to MetaCart
Anisotropic diffusion filtering is widely used for MR image enhancement. However, the anisotropic filter is nonoptimal for MR images with spatially varying noise levels, such as images reconstructed from sensitivityencoded data and intensity inhomogeneitycorrected images. In this work, a new method for filtering MR images with spatially varying noise levels is presented. In the new method, a priori information regarding the image noise level spatial distribution is utilized for the local adjustment of the anisotropic diffusion filter. Our new method was validated and compared with the standard filter on simulated and real MRI data. The noiseadaptive method was demonstrated to outperform the standard anisotropic diffusion filter in both image error reduction and image signaltonoise ratio (SNR) improvement. The method was also applied to inhomogeneitycorrected and sensitivity encoding (SENSE) images.
Fast, iterative image reconstruction for MRI in the presence of field inhomogeneities
 IEEE Trans. Med. Imaging
, 2003
"... Abstract—In magnetic resonance imaging, magnetic field inhomogeneities cause distortions in images that are reconstructed by conventional fast Fourier trasform (FFT) methods. Several noniterative image reconstruction methods are used currently to compensate for field inhomogeneities, but these metho ..."
Abstract

Cited by 2 (1 self)
 Add to MetaCart
Abstract—In magnetic resonance imaging, magnetic field inhomogeneities cause distortions in images that are reconstructed by conventional fast Fourier trasform (FFT) methods. Several noniterative image reconstruction methods are used currently to compensate for field inhomogeneities, but these methods assume that the field map that characterizes the offresonance frequencies is spatially smooth. Recently, iterative methods have been proposed that can circumvent this assumption and provide improved compensation for offresonance effects. However, straightforward implementations of such iterative methods suffer from inconveniently long computation times. This paper describes a tool for accelerating iterative reconstruction of fieldcorrected MR images: a novel timesegmented approximation to the MR signal equation. We use a min–max formulation to derive the temporal interpolator. Speedups of around 60 were achieved by combining this temporal interpolator with a nonuniform fast Fourier transform with normalized root mean squared approximation errors of 0.07%. The proposed method provides fast, accurate, fieldcorrected image reconstruction even when the field map is not smooth. Index Terms—Field inhomogeneity correction, image reconstruction, iterative methods, magnetic resonance imaging, temporal interpolation, time segmentation. I.
High Performance Computers
, 1991
"... ACM, (2008). This is the author’s version of the work. It is posted here by permission of ACM for your personal use. Not for ..."
Abstract

Cited by 2 (0 self)
 Add to MetaCart
ACM, (2008). This is the author’s version of the work. It is posted here by permission of ACM for your personal use. Not for
Improvements in magnetic resonance imaging using information redundancy
, 2005
"... This thesis describes a number of algorithms related to the acquisition, reconstruction and postprocessing of Magnetic Resonance data. The basic theme underlying each of these algorithms is the use of a unified systems approach to exploit information redundancy available in MR imaging. There are ..."
Abstract

Cited by 1 (1 self)
 Add to MetaCart
This thesis describes a number of algorithms related to the acquisition, reconstruction and postprocessing of Magnetic Resonance data. The basic theme underlying each of these algorithms is the use of a unified systems approach to exploit information redundancy available in MR imaging. There are three basic contributions. The first concerns the development of a new motion correction algorithm for TimeResolved MR Angiography. Motion artifacts in angiography data are very difficult to remove without affecting vascular evolution. Our algorithm uses successive POCS iterations to remove unwanted artifacts without degrading quality. Doubleblind testing has indicated significant improvement over angiograms created manually by experienced radiologist. In summary, our method seeks to exploit temporal redundancy to remove motion artifacts. The second contribution is our recent work on Parallel MR imaging in presence of sensitivity errors using a Maximum Likelihood technique. It can be shown that standard phased array reconstruction using popular parallel imaging methods is inappropriate in presence of errors in measuring sensitivity maps of coils. Since
FAST REGULARIZED RECONSTRUCTION OF NONUNIFORMLY SUBSAMPLED PARALLEL MRI DATA
"... Parallel MR imaging is an effective approach to reduce MR image acquisition time. Nonuniform subsampling allows one to tailor the subsampling scheme for improved image quality at high acceleration factors. However, nonuniform subsampling precludes fast reconstruction schemes such as SENSE, and is ..."
Abstract

Cited by 1 (1 self)
 Add to MetaCart
Parallel MR imaging is an effective approach to reduce MR image acquisition time. Nonuniform subsampling allows one to tailor the subsampling scheme for improved image quality at high acceleration factors. However, nonuniform subsampling precludes fast reconstruction schemes such as SENSE, and is more likely to require a regularized solution than reconstruction of uniformly subsampled data demands. This means that one needs to choose a good regularization parameter, typically requiring multiple expensive system solves. Here, we present an efficient LSQRHybrid algorithm which simultaneously addresses the need for rapid regularization parameter selection and fast reconstruction. This algorithm can reconstruct nonuniformly subsampled parallel MRI data, with automatic regularization and good image quality, in a time competitive with Cartesian SENSE. 1.