next up previous contents
Next: Conical tilt geometry Up: 3D reconstruction formulas Previous: Orlov's inversion formula

Colsher's inversion formula

Assume again that g = Pf is known for tex2html_wrap_inline3968 .

We want to derive an inversion procedure similar to the one in section 2.1. With the backprojection

displaymath3970

we have again

displaymath3972

provided that tex2html_wrap_inline3974 . Again the convolutions on the left hand right hand side have different meanings. Explicitly this reads

  equation1853

which corresponds to (2.2). As in (2.2) we express the relationship tex2html_wrap_inline3974 in Fourier space, obtaining

displaymath3978

see Colsher (1980). In order to get an inversion formula for P we have to determine v such that tex2html_wrap_inline3984 or tex2html_wrap_inline3986 , i.e.

  equation1868

A solution tex2html_wrap_inline3988 independent of tex2html_wrap_inline3298 is

displaymath3992

where tex2html_wrap_inline3994 is the length of tex2html_wrap_inline3996 . For the spherical zone tex2html_wrap_inline3910 from section 3.4 with tex2html_wrap_inline4000 , tex2html_wrap_inline4002 , tex2html_wrap_inline4004 a constant with tex2html_wrap_inline4006 , Colsher computed tex2html_wrap_inline4008 explicitly. With tex2html_wrap_inline4010 be obtained

  equation1881

Filters such as the Colsher filter (3.10) do not have small support. This means that g in (3.8) has to be known in all of tex2html_wrap_inline3808 . Often g is only available in part of tex2html_wrap_inline3808 (truncated projections). Let us choose

  equation1896

where tex2html_wrap_inline4020 from section (3.4) and tex2html_wrap_inline4022 is a horizontal unit vector. Since (3.11) is constant in the vertical direction, v is a tex2html_wrap_inline3098 -function in the vertical variable. Hence the integral on the right hand side of (3.8) reduces to an integral over horizontal lines in tex2html_wrap_inline3808 , making it possible to handle truncated projections. Unfortunately, (3.11) does not quite satisfy tex2html_wrap_inline3986 , i.e. it does not provide an exact inversion. Instead we only have

displaymath4032

where tex2html_wrap_inline4034 , tex2html_wrap_inline4036 . This is close to tex2html_wrap_inline3986 if tex2html_wrap_inline4004 is small. In this case reconstruction from truncated projections is possible, at least approximately.



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