fanbeam matlab,Fan-beam transform

fanbeam

Description

F = fanbeam(I,D)

computes the fan-beam projection data (sinogram)

F from the image I.

Each column of F contains fan-beam projection

data at one rotation angle. D is the distance

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

F = fanbeam(I,D,Name,Value)

uses name-value pairs to specify the rotation increment and sensor

spacing.

fan_sensor_positions and the rotation

angles where the fan-beam projections are calculated in

fan_rotation_angles.

Examples

Compute Fan-beam Projections for Rotation Angles Over Entire Image

Set the IPT preference to make the axes visible.

iptsetpref('ImshowAxesVisible','on')

Create a sample image and display it.

ph = phantom(128);

imshow(ph)

Calculate the fanbeam projections and display them.

[F,Fpos,Fangles] = fanbeam(ph,250);

figure

imshow(F,[],'XData',Fangles,'YData',Fpos,...

'InitialMagnification','fit')

axis normal

xlabel('Rotation Angles (degrees)')

ylabel('Sensor Positions (degrees)')

colormap(gca,hot), colorbar

Compute Radon and Fan-beam Projections and Compare Results

Compute fan-beam projections for 'arc' geometry.

I = ones(100