Note: This page has been translated by MathWorks. Click here to see

To view all translated materials including this page, select Country from the country navigator on the bottom of this page.

To view all translated materials including this page, select Country from the country navigator on the bottom of this page.

Inverse fan-beam transform

`I = ifanbeam(F,D)`

`I = ifanbeam(F,D,Name,Value)`

`[I,H] = ifanbeam(___)`

reconstructs the image `I`

= ifanbeam(`F`

,`D`

)`I`

from fan-beam
projection data in `F`

. Each column of
`F`

contains fan-beam projection data at
one rotation angle. The angle between sensors is assumed to be
uniform and equal to the increment between fan-beam rotation angles.
`D`

is the distance from the fan-beam
vertex to the center of rotation.

uses name-value pairs to control various aspects of the
reconstruction. Argument names can be abbreviated, and case does not
matter.`I`

= ifanbeam(`F`

,`D`

,`Name,Value`

)

`ifanbeam`

converts the fan-beam data to parallel beam
projections and then uses the filtered back projection algorithm to perform
the inverse Radon transform. The filter is designed directly in the
frequency domain and then multiplied by the FFT of the projections. The
projections are zero-padded to a power of 2 before filtering to prevent
spatial domain aliasing and to speed up the FFT.

[1] Kak, A. C., and M. Slaney, *Principles of
Computerized Tomographic Imaging,* New York, NY, IEEE
Press, 1988.