## Iterative tomographic image reconstruction using Fourier-based forward and back- projectors (2004)

Venue: | IEEE Trans. Med. Imag |

Citations: | 22 - 3 self |

### BibTeX

@ARTICLE{Fessler04iterativetomographic,

author = {Jeffrey A. Fessler},

title = {Iterative tomographic image reconstruction using Fourier-based forward and back- projectors},

journal = {IEEE Trans. Med. Imag},

year = {2004},

volume = {23},

pages = {412}

}

### OpenURL

### Abstract

Fourier-based reprojection methods have the potential to reduce the computation time in iterative tomographic image reconstruction. Interpolation errors are a limitation of Fourier-based reprojection methods. We apply a min-max interpolation method for the nonuniform fast Fourier transform (NUFFT) to minimize the interpolation errors. Numerical results show that the min-max NUFFT approach provides substantially lower approximation errors in tomographic reprojection and backprojection than conventional interpolation methods.

### Citations

403 | The Mathematics of Computerized Tomography - Natterer - 2001 |

115 | Fourier volume rendering - Malzbender - 1993 |

95 |
The Fourier reconstruction of a head section
- Shepp, Logan
- 1974
(Show Context)
Citation Context ...se, corresponding to the adjoint of the forward operator, which is a form of backprojection. IV. RESULTS We evaluated the accuracy of the NUFFT-based reprojector using the Shepp-Logan digital phantom =-=[2, 55]-=- with N1 ×N2 = 128×128 pixels as shown in Fig. 1. We simulated a parallel-beam tomographic system with a sinogramsFessler NUFFT-Reprojection 5 size of Nr = 160 radial bins by Nθ = 192 angles over 180 ... |

83 | Nonuniform fast Fourier transforms using Min-Max interpolation
- Fessler, Sutton
(Show Context)
Citation Context ... min-max NUFFT method. Section III describes the application of the NUFFT method to reprojection. Section IV gives numerical results. II. NUFFT IN 2D The section reviews the NUFFT method presented in =-=[11]-=-. We focus on the 2D case since [11] emphasized only the 1D case. The extension to 3D is straightforward. A. Problem statement We are given a 2D finite-extent signal x[n1,n2], n1 = 0,...,N1−1, n2 =0,.... |

76 | Volume Rendering - Levoy - 1990 |

31 | The gridding method for image reconstruction by Fourier transformation - Schomberg, Timmer - 1995 |

31 | Fast reprojection of volume data - Dunne, Napel, et al. - 1990 |

30 | Rapid 3-D cone-beam reconstruction with the simultaneous algebraic reconstruction technique (SART) using 2-D texture mapping hardware - Mueller, Yagel - 2000 |

21 | Reconstruction algorithms: Transform methods - Lewitt - 1983 |

20 | Direct Fourier reconstruction in computer tomography - Stark, Woods, et al. - 1981 |

20 |
A new linogram algorithm for computerized tomography
- Potts, Steidl
(Show Context)
Citation Context ...ication of a NUFFT method to iterative tomographic image reconstruction. Potts et al. have applied the NUFFT to direct Fourier reconstruction using truncated gaussian bell as the interpolation kernel =-=[42,43]-=-. (A similar approach was described by Tabei and Ueda without the NUFFT moniker [32].) That method is undoubtably very good, but no claims of optimality can be made. Our min-max method is optimal for ... |

19 | A computational study of reconstruction algorithms for diffraction tomography: Interpolation versus filtered backpropagation - Pan, Kak - 1983 |

19 | An optimal and efficient new gridding algorithm using singular value decomposition - Rosenfeld - 1998 |

17 | On the direct Fourier method for computer tomography - Gottlieb, Gustafsson, et al. |

15 | A fast and accurate Fourier algorithm for iterative parallel-beam tomography - Delaney, Bresler - 1996 |

13 |
2 log2 N) Filtered backprojection reconstruction algorithm for tomography
- Basu, Bresler, et al.
- 2000
(Show Context)
Citation Context ...re the principal computational bottleneck in iterative reconstruction methods, particularly in 3D image reconstruction. A variety of methods for accelerating this process have been proposed, see e.g.,=-=[1]-=-. One natural approach to forward projection is to use the Fourier-slice theorem [2, p. 56]. This theorem relates the 1D FT of each projection to samples of the 2D FT of the object on a polar grid 1 .... |

13 | A new polar Fourier-transform for computeraided tomography and spotlight synthetic aperture radar - Lawton - 1988 |

10 | Macovski A. Fast spiral coronary artery imaging. Magn Reson Med - Meyer, Hu, et al. - 1992 |

9 | Direct-Fourier reconstruction in tomography and synthetic aperture radar, Int
- Choi, Jr, et al.
- 1998
(Show Context)
Citation Context ...ke no claim of optimality of our min-max interpolator for the direct Fourier reconstruction method. We conjecture that a min-max approach could be devised for the direct Fourier reconstruction method =-=[39]-=-. Section II reviews the min-max NUFFT method. Section III describes the application of the NUFFT method to reprojection. Section IV gives numerical results. II. NUFFT IN 2D The section reviews the NU... |

8 |
3D-FRP: Direct Fourier reconstruction with Fourier reprojection for fully 3-D
- Matej, Lewitt
- 2001
(Show Context)
Citation Context ...useful for quickly forming an initial image estimate for subsequent iterations by another algorithm based on statistical weighting. Very recently, Matej et al. have applied Fourier-based reprojection =-=[14]-=- for (noniterative) fully 3D PET reconstruction and [15] for calculation of attenuation correction factors in PET. They used Kaiser-Bessel windows for interpolation, which are known to be reasonably a... |

8 | On the Use of Graphics Hardware to Accelerate Algebraic Reconstruction Methods - Mueller, Yagel - 1999 |

7 | New Fourier reconstruction algorithms for computerized tomography
- Potts, Steidl
- 2000
(Show Context)
Citation Context ...ication of a NUFFT method to iterative tomographic image reconstruction. Potts et al. have applied the NUFFT to direct Fourier reconstruction using truncated gaussian bell as the interpolation kernel =-=[42,43]-=-. (A similar approach was described by Tabei and Ueda without the NUFFT moniker [32].) That method is undoubtably very good, but no claims of optimality can be made. Our min-max method is optimal for ... |

6 | Recovering Multidimensional Signals from their Projections - Mersereau - 1973 |

6 | Direct Fourier transform techniques in 3-D image reconstruction - Mersereau - 1976 |

6 | Analysis of the direct Fourier method for computer tomography - Walden - 2000 |

5 |
A min-max approach to the multidimensional nonuniform FFT: Application to tomographic image reconstruction
- Fessler, Sutton
- 2001
(Show Context)
Citation Context ...y claims we make are for the forward problem only. The prospect of using graphics hardware to assist in this projection process is enticing [50–54]. To our knowledge, the work herein and presented in =-=[8]-=- represents the first application of a NUFFT method to iterative tomographic image reconstruction. Potts et al. have applied the NUFFT to direct Fourier reconstruction using truncated gaussian bell as... |

5 | A high-speed computerized tomography image reconstruction using direct two-dimensional Fourier transform method - Niki, Mizutani, et al. - 1983 |

5 | Fourier reconstruction in tomography - Natterer - 1985 |

4 |
Threedimensional image reconstruction in the Fourier domain
- Stearns, Chesler, et al.
- 1987
(Show Context)
Citation Context ...e use of the Fourier-slice theorem as a tool for reprojection was noted in the late 1980’s by Crawford et al. [3, 4], in the context of correction of beam hardening in X-ray CT, and by Stearns et al. =-=[5, 6]-=-, in the context of filling in missing projections in (noniterative) 3D PET image reconstruction. These approaches were apparently largely abandoned thereafter due in part to unacceptable image artifa... |

4 | Direct Fourier reconstruction in fan-beam tomography - Peng, Stark - 1987 |

3 |
Performance of 3D RAMLA with Smooth Basis Functions on Fully 3D
- Matej, Daube-Witherspoon, et al.
- 2001
(Show Context)
Citation Context ... subsequent iterations by another algorithm based on statistical weighting. Very recently, Matej et al. have applied Fourier-based reprojection [14] for (noniterative) fully 3D PET reconstruction and =-=[15]-=- for calculation of attenuation correction factors in PET. They used Kaiser-Bessel windows for interpolation, which are known to be reasonably accurate [16,17] but without explicitly evaluating the ac... |

3 |
Backprojection by upsampled Fourier series expansion and interpolated FFT
- Tabei, Ueda
- 1992
(Show Context)
Citation Context ...ave applied the NUFFT to direct Fourier reconstruction using truncated gaussian bell as the interpolation kernel [42,43]. (A similar approach was described by Tabei and Ueda without the NUFFT moniker =-=[32]-=-.) That method is undoubtably very good, but no claims of optimality can be made. Our min-max method is optimal for the NUFFT need for reprojection, but we make no claim of optimality of our min-max i... |

3 |
On the possibility of direct Fourier reconstruction from divergent-beam projections
- Alliney, Matej, et al.
- 1993
(Show Context)
Citation Context ...not an obvious approach to forming fan-beam projections, at least not without further interpolations that conceivably may degrade spatial resolution. An exception would be when the fan angle is small =-=[56]-=-. Perhaps a minmax approach could be found for the parallel-fan interpolation, perhaps by some change of variables. Thus the applicability to X-ray CT is an open question. Cone-beam CT is an additiona... |

2 |
System for reprojecting images using transform techniques
- Crawford
- 1986
(Show Context)
Citation Context ... the 1D FT of each projection to samples of the 2D FT of the object on a polar grid 1 . The use of the Fourier-slice theorem as a tool for reprojection was noted in the late 1980’s by Crawford et al. =-=[3, 4]-=-, in the context of correction of beam hardening in X-ray CT, and by Stearns et al. [5, 6], in the context of filling in missing projections in (noniterative) 3D PET image reconstruction. These approa... |

2 |
Accelerated image reconstruction for a cylindrical positron tomograph using Fourier domain methods
- Stearns, Chesler, et al.
- 1990
(Show Context)
Citation Context ...e use of the Fourier-slice theorem as a tool for reprojection was noted in the late 1980’s by Crawford et al. [3, 4], in the context of correction of beam hardening in X-ray CT, and by Stearns et al. =-=[5, 6]-=-, in the context of filling in missing projections in (noniterative) 3D PET image reconstruction. These approaches were apparently largely abandoned thereafter due in part to unacceptable image artifa... |

2 |
Voir: A volumetric image reconstruction algorithm based on fourier techniques for inversion of the 3-d radon transform
- Dusaussoy
- 1996
(Show Context)
Citation Context ...minmax approach could be found for the parallel-fan interpolation, perhaps by some change of variables. Thus the applicability to X-ray CT is an open question. Cone-beam CT is an additional challenge =-=[36]-=-. The Fourier approach appears unsuitable for SPECT since it cannot easily accommodate nonuniform attenuation and depth-dependent detector response. Perhaps some form of the frequency distance princip... |

2 | Fast direct Fourier methods, based on one- and two-pass coordinate transformations, yield accurate reconstructions of x-ray CT clinical images,” Phys - Bellon, Lanzavecchia - 1997 |

2 | Iterative reconstruction in diffraction tomography using NUFFT - Bronstein, Bronstein, et al. - 2002 |

1 | Mersereau and A V oppenheim, “Digital reconstruction of multidimensional signals from their projections - M - 1974 |

1 | Electron tomography in conical tilt geometry. The accuracy of a direct Fourier method (DFM) and the suppression of nontomographic noise - Lanzavecchia, Bellon - 1996 |

1 |
Cox and Raoqiong Tong, “Two- and threedimensional image rotation using the FFT
- W
- 1999
(Show Context)
Citation Context ...tion and depth-dependent detector response. Perhaps some form of the frequency distance principle could overcome this limitation? Combining with fast rotation methods may be another avenue to explore =-=[57]-=-. The simplest version of the method is only suitable for shift-invariant radial blur, which is only an approximation to the PSF of PET systems. However, one could combine ideal “line integral” reproj... |