Polar sampling in k‐space: Reconstruction effects

Magnetic resonance images are most commonly computed by taking the inverse Fourier transform of the k‐space data. This transformation can potentially create artifacts in the im age, depending on the reconstruction algorithm used. For equally spaced radial and azimuthal k‐space polar sampling, both gridding and convolution backprojection are applicable. However, these algorithms potentially can yield different res olution, signal‐to‐noise ratio, and aliasing characteristics in the reconstructed image. Here, these effects are analyzed and their tradeoffs are discussed. It is shown that, provided the modulation transfer function and the signal‐to‐noise ratio are considered together, these algorithms perform similarly. In contrast, their aliasing behavior is different, since their re spective point spread functions (PSF) differ. In gridding, the PSF is composed of the mainlobe and ringlobes that lead to aliasing. Conversely, there are no ringlobes in the convolution backprojection PSF, thus radial aliasing effects are mini mized. Also, a hybrid gridding and convolution backprojection reconstruction is presented for radially nonequidistant k‐space polar sampling.

[1]  Bob S. Hu,et al.  Fast Spiral Coronary Artery Imaging , 1992, Magnetic resonance in medicine.

[2]  J. D. O'Sullivan,et al.  A Fast Sinc Function Gridding Algorithm for Fourier Inversion in Computer Tomography , 1985, IEEE Transactions on Medical Imaging.

[3]  W. Edelstein,et al.  The intrinsic signal‐to‐noise ratio in NMR imaging , 1986, Magnetic resonance in medicine.

[4]  J. Tukey,et al.  An algorithm for the machine calculation of complex Fourier series , 1965 .

[5]  S. Rowland,et al.  Computer implementation of image reconstruction formulas , 1979 .

[6]  Terry M. Peters,et al.  Image reconstruction from finite numbers of projections , 1973 .

[7]  B K Rutt,et al.  Effects of polar sampling in k‐space , 1996, Magnetic resonance in medicine.

[8]  G. N. Ramachandran,et al.  Three-dimensional reconstruction from radiographs and electron micrographs: application of convolutions instead of Fourier transforms. , 1971, Proceedings of the National Academy of Sciences of the United States of America.

[9]  A. Macovski Noise in MRI , 1996, Magnetic resonance in medicine.

[10]  Jeffrey L. Duerk,et al.  Analytical resolution and noise characteristics of linearly reconstructed magnetic resonance data with arbitraryk‐space sampling , 1995, Magnetic resonance in medicine.

[11]  A. C. Riddle,et al.  Inversion of Fan-Beam Scans in Radio Astronomy , 1967 .

[12]  R. Bracewell Strip Integration in Radio Astronomy , 1956 .

[13]  Ronald N. Bracewell,et al.  The Main Beam and Ring Lobes of an East-West Rotation-Synthesis Array , 1973 .

[14]  A. Macovski,et al.  Selection of a convolution function for Fourier inversion using gridding [computerised tomography application]. , 1991, IEEE transactions on medical imaging.

[15]  C E Metz,et al.  Transfer function analysis of radiographic imaging systems. , 1979, Physics in medicine and biology.

[16]  R. Bracewell The Fourier Transform and Its Applications , 1966 .