next up previous contents
Next: Colsher's inversion formula Up: 3D reconstruction formulas Previous: Grangeat's formula

Orlov's inversion formula

This formula inverts the ray transform

  equation1822

which comes up e.g. in 3D emission tomography (PET, Defrise et al. (1989)). If tex2html_wrap_inline3298 is restricted to a plane, then we have simply the Radon transform in this plane, and we can reconstruct f in that plane by any of the methods in the previous section. In practice g = Pf is measured for tex2html_wrap_inline3906 where tex2html_wrap_inline3908 . In Orlov's formula (Orlov (1976)), tex2html_wrap_inline3910 is a spherical zone around the equator, i.e.

displaymath3912

where tex2html_wrap_inline3302 , tex2html_wrap_inline3916 are the spherical coordinates of tex2html_wrap_inline3918 and tex2html_wrap_inline3920 . Then,

  eqnarray1830

where tex2html_wrap_inline3922 is the Laplacian acting on x and tex2html_wrap_inline3926 is the length of the intersection of tex2html_wrap_inline3910 with the plane spanned by tex2html_wrap_inline3930 . The first formula of (3.6) is - up to tex2html_wrap_inline3922 - a backprojection, while the second one a convolution in tex2html_wrap_inline3808 . Thus an implementation of (3.6) is again a filtered backprojection algorithm.

P can also be inverted by the Fourier transform. We have

  equation1842

where `` tex2html_wrap_inline3108 '' denotes the (n-1)-dimensional Fourier transform in tex2html_wrap_inline3808 on the left hand side and the Fourier transform in tex2html_wrap_inline3076 on the right hand side.

Assume that tex2html_wrap_inline3908 satisfies the Orlov condition: Every equatorial circle of tex2html_wrap_inline3784 meets tex2html_wrap_inline3910 . Note that the set tex2html_wrap_inline3910 - the spherical zone - we used above in Orlov's formula satisfies this condition. From (3.7) it follows that f is uniquely determined by tex2html_wrap_inline3956 for tex2html_wrap_inline3906 under the Orlov condition. Namely if tex2html_wrap_inline3120 is arbitrary, then Orlov's condition says that there exists tex2html_wrap_inline3962 , and tex2html_wrap_inline3964 is determined from (3.7).



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