iradon - Inverse Radon transform
Syntax
I = iradon(R, theta)
I = iradon(P, theta, interp,
filter, frequency_scaling, output_size)
[I,H] = iradon(...)
Description
I = iradon(R, theta) reconstructs the image I
from projection data in the two-dimensional array R. The
columns of R are parallel beam projection data.
iradon assumes that the center of rotation is the center
point of the projections, which is defined as
ceil(size(R,1)/2).
theta describes the angles (in degrees) at which the
projections were taken. It can be either a vector containing the
angles or a scalar specifying D_theta, the incremental
angle between projections. If theta is a vector, it must
contain angles with equal spacing between them. If theta
is a scalar specifying D_theta, the projections were taken
at angles theta = m*D_theta, where m =
0,1,2,...,size(R,2)-1. If the input is the empty matrix
([]), D_theta defaults to
180/size(R,2).
iradon uses the filtered back-projection algorithm to
perform the inverse Radon transform. The filter is designed
directly in the freq