next up previous contents
Next: Orlov's inversion formula Up: 3D reconstruction formulas Previous: The FDK approximate formula

Grangeat's formula

Grangeat's formula requires sources on a curve C with the following property: Each plane meeting supp(f) contains at least one source. This condition is obviously not satisfied for C a plane circle for which the FDK approximation has been derived.

The data for Grangeat's formula is

displaymath3862

The condition on the source curve means that for each x with tex2html_wrap_inline3866 and each tex2html_wrap_inline3868 there exists a source tex2html_wrap_inline3870 such that tex2html_wrap_inline3872 .

The gist of Grangeat's inversion is a relation between g and the 3D Radon transform Rf of f. This relation reads (Grangeat (1991))

  equation1803

where tex2html_wrap_inline3880 stands for the derivative in the direction tex2html_wrap_inline3868 , acting on the second argument. For this to make sense we have to extend g to all of tex2html_wrap_inline3886 by using the above definition not only for tex2html_wrap_inline3868 , but for all of tex2html_wrap_inline3890 . This is equivalent to extending g by homogeneity of degree -1 in the second argument. With help of the 3D inversion formula

displaymath3896

for the 3D Radon transform, Grangeat's formula leads immediately to an inversion procedure for the data g. Related inversion formulas for cone beam tomography have been derived by Tuy (1983), Gelfand and Goncharov (1987). For the details of the implementation see Defrise and Clack (1995).



Frank Wuebbeling
Thu Sep 10 10:51:17 MET DST 1998