Radon data interpolation is a necessary procedure in computed tomography (CT), especially for reconstruction from divergent beam scanning. In a polar-grid representation, the Radon data of a fanbeam projection are populated on an arc, rather on a radial line. Collectively, the Radon data generated from a fanbeam CT system are unevenly populated: The population becomes sparser as the polar distance increases. In CT reconstruction, the Fourier central slice theorem requires a radial scanline full of Radon data. Therefore the vacant entries of a scanline must be filled by interpolation. In addition, interpolation is also required in polar-to-Cartesian conversion. In this paper we propose a practical interpolation technique for filling the vacant entries by local convex combination. It is a linear interpolant that generates a value for a grid point from the available data lying in its neighborhood, by a weighted average, with the weights corresponding to the inverse distances. In fact, the linear convex combination serves as a general flat-smoothing operation in filling a vacancy. Specifically, this technique realizes a variety of linear interpolations, including nearest-neighbor replication, two-point collinear, three-point triangulation, and four-point quadrilateral, and local extrapolation, in a unified framework. Algorithms and a simulation demonstration are provided.
© 2003 Optical Society of AmericaFull Article | PDF Article
L. A. Feldkamp, L. C. Davis, and J. W. Kress
J. Opt. Soc. Am. A 1(6) 612-619 (1984)
Nagaaki Ohyama, Shigeru Ohki, Soichi Inoue, Jumpei Tsujiuchi, and Toshio Honda
J. Opt. Soc. Am. A 4(2) 318-324 (1987)
Hiroyuki Kudo and Tsuneo Saito
J. Opt. Soc. Am. A 7(12) 2169-2183 (1990)