The gridding method for image reconstruction by Fourier transformation
- Philips Research Lab., Hamburg (Germany)
- Philips Medical Systems, Best (Netherlands)
This paper explores a computational method for reconstructing an n-dimensional signal f from a sampled version of its Fourier transform {cflx f}. The method involves a window function {cflx w} and proceeds in three steps. First, the convolution {cflx g} = {cflx w} * {cflx f} is computed numerically on a Cartesian grid, using the available samples of {cflx f}. Then, g = wf is computed via the inverse discrete Fourier transform, and finally f is obtained as g/w. Due to the smoothing effect of the convolution, evaluating {cflx w} * {cflx f} is much less error prone than merely interpolating {cflx f}. The method was originally devised for image reconstruction in radio astronomy, but is actually applicable to a broad range of reconstructive imaging methods, including magnetic resonance imaging and computed tomography. In particular, it provides a fast and accurate alternative to the filtered backprojection. The basic method has several variants with other applications, such as the equidistant resampling of arbitrarily sampled signals or the fast computation of the Radon (Hough) transform.
- OSTI ID:
- 128801
- Journal Information:
- IEEE Transactions on Medical Imaging, Journal Name: IEEE Transactions on Medical Imaging Journal Issue: 3 Vol. 14; ISSN 0278-0062; ISSN ITMID4
- Country of Publication:
- United States
- Language:
- English
Similar Records
Rotational Millimeter-Wave Shoe Scanner Using the Discrete Fourier Transform for Backprojection-Based Image Reconstruction
Fourier transform resampling: Theory and application