Parallel data resampling and Fourier inversion by the scan-line method

Fourier inversion is an efficient method for image reconstruction in a variety of applications, for example, in computed tomography and magnetic resonance imaging. Fourier inversion normally consists of two steps, interpolation of data onto a rectilinear grid, if necessary, and inverse Fourier transformation. Here, the authors present interpolation by the scan-line method, in which the interpolation algorithm is implemented in a form consisting only of row operations and data transposes. The two-dimensional inverse Fourier transformation can also be implemented with only row operations and data transposes. Accordingly, Fourier inversion can easily be implemented on a parallel computer that supports row operations and data transposes on row distributed data. The conditions under which the scan-line implementations are algorithmically equivalent to the original serial computer implementation are described and methods for improving accuracy outside of those conditions are presented. The scan-line algorithm is implemented on the iWarp parallel computer using the Adapt language for parallel image processing. This implementation is applied to magnetic resonance data acquired along radial-lines and spiral trajectories through Fourier transform space.

[1]  Jack Dongarra,et al.  Pvm 3 user's guide and reference manual , 1993 .

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

[3]  Orhan Nalcioglu,et al.  Spectroscopic imaging by quadrature modulated echo time shifting , 1986 .

[4]  S. Ljunggren.  A simple graphical representation of fourier-based imaging methods , 1983 .

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

[6]  C. K. Yuen,et al.  Digital Filters , 1979, IEEE Transactions on Systems, Man, and Cybernetics.

[7]  Douglas C. Noll,et al.  Deblurring for non‐2D fourier transform magnetic resonance imaging , 1992, Magnetic resonance in medicine.

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

[9]  Gabor T. Herman,et al.  Image Reconstruction From Projections , 1975, Real Time Imaging.

[10]  D C Noll,et al.  MR imaging of lung parenchyma: a solution to susceptibility. , 1992, Radiology.

[11]  R. Hingorani,et al.  An Investigation of Computerized Tomography by Direct Fourier Inversion and Optimum Interpolation , 1981 .

[12]  Jon A. Webb.  Steps toward architecture-independent image processing , 1992, Computer.

[13]  Alvy Ray Smith,et al.  3-D transformations of images in scanline order , 1980, SIGGRAPH '80.

[14]  George Wolberg,et al.  Digital image warping , 1990 .

[15]  J. Pauly,et al.  A homogeneity correction method for magnetic resonance imaging with time-varying gradients. , 1991, IEEE transactions on medical imaging.

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

[17]  Terrance E. Boult,et al.  Separable image warping with spatial lookup tables , 1989, SIGGRAPH '89.

[18]  R. Lewitt.  Reconstruction algorithms: Transform methods , 1983, Proceedings of the IEEE.