next up previous contents
Next: 3D reconstruction formulas Up: The filtered backprojection algorithm Previous: Linear fan beam geometry

Fast backprojection

The backprojection (i.e. Step 2 in Algorithms 1-3) is the most time consuming part of the filtered backprojection algorithm. The filtering or convolution step (Step 1 in Algorithms 1-3) requires in principle the same number of operations, but this can easily be reduced drastically either by cutting off the filter tex2html_wrap_inline3364 or by using the fast Fourier transform (FFT).

The backprojection consists is the evaluation of the sums

displaymath3694

on a tex2html_wrap_inline3696 grid. This is the simplest case of Algorithm 1, the resolution of the image being adjusted to the number of views p according to the sampling theorem. Nilsson (1997) suggested a devide and conquer strategy for doing this with O tex2html_wrap_inline3700 operations, as opposed to the O tex2html_wrap_inline3702 operations of a direct evaluation. Suppose tex2html_wrap_inline3704 .

Step 1: For tex2html_wrap_inline3706 compute

displaymath3708

Since tex2html_wrap_inline3710 is constant along the lines tex2html_wrap_inline3712 it suffices to compute tex2html_wrap_inline3714 at 2p points x. We need tex2html_wrap_inline3720 operations.

Step 2: For tex2html_wrap_inline3722 compute

displaymath3724

Since tex2html_wrap_inline3710 , tex2html_wrap_inline3728 are constant along the lines tex2html_wrap_inline3712 , tex2html_wrap_inline3732 , respectively, tex2html_wrap_inline3734 is almost constant along the lines tex2html_wrap_inline3736 where tex2html_wrap_inline3738 . Hence it suffices to compute tex2html_wrap_inline3734 only for a few, say 2, points on each of those lines. This means that we have to evaluate tex2html_wrap_inline3734 at 4p points, requiring tex2html_wrap_inline3746 operations.

Step 3: For tex2html_wrap_inline3748 compute

displaymath3750

With the same reasoning as in step 1 and 2 we find that it suffices to compute tex2html_wrap_inline3752 for only for 8p points, requiring tex2html_wrap_inline3756 operations.

Proceeding in this fashion we arrive in step m at the approximation tex2html_wrap_inline3760 to f, which has to be evaluated at tex2html_wrap_inline3764 points. Hence the number of operations in step 1 to m is

displaymath3768

and this is O tex2html_wrap_inline3700 . Of course this derivation is heuristic, and we simply ignored the necessary interpolations and approximations. But practical experience demonstrates that such an algorithm can be made work.


next up previous contents
Next: 3D reconstruction formulas Up: The filtered backprojection algorithm Previous: Linear fan beam geometry

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