Derivation and implementation of a cone-beam reconstruction algorithm for nonplanar orbits

B.D. Smith (ibid., vol.MI-4, p.15-25, 1985; Opt. Eng., vol.29, p.524-34, 1990) and P. Grangeat (These de doctorat, 1987; Lecture Notes in Mathematics 1497, p.66-97, 1991) derived a cone-beam inversion formula that can be applied when a nonplanar orbit satisfying the completeness condition is used. Although Grangeat's inversion formula is mathematically different from Smith's one, they have similar overall structures to each other. The contribution of the present paper is two-fold. First, based on the derivation of Smith, the authors point out that Grangeat's inversion formula and Smith's one can be conveniently described using a single formula (the Smith-Grangeat inversion formula) that is in the form of space-variant filtering followed by cone-beam back projection. Furthermore, the resulting formula is reformulated for data acquisition systems with a planar detector to obtain a new reconstruction algorithm. Second, the authors make two significant modifications to the new algorithm to reduce artifacts and numerical errors encountered in direct implementation of the new algorithm. As for exactness of the new algorithm, the following fact can be stated. The algorithm based on Grangeat's intermediate function is exact for any complete orbit, whereas that based on Smith's intermediate function should be considered as an approximate inverse excepting the special case where almost every plane in 3D space meets the orbit. The validity of the new algorithm is demonstrated by simulation studies.