主要内容

fanbeam

Fan-beam transform

描述

example

F= fanbeam(I,D)computes the fan-beam projection data (sinogram)Ffrom the imageI. Each column ofFcontains fan-beam projection data at one rotation angle.D是the distance from the fan-beam vertex to the center of rotation.

F= fanbeam(I,D,名称,价值)uses name-value pairs to specify the rotation increment and sensor spacing.

example

[F,fan_sensor_positions,fan_rotation_angles] = fanbeam(___)returns the location of fan-beam sensors infan_sensor_positionsand the rotation angles where the fan-beam projections are calculated infan_rotation_angles.

例子

全部收缩

Set the IPT preference to make the axes visible.

iptsetpref('ImshowAxesVisible','on')

Create a sample image and display it.

pH =幻影(128);imshow(ph)

图包含一个轴对象。The axes object contains an object of type image.

Calculate the fanbeam projections and display them.

[F,Fpos,Fangles] = fanbeam(ph,250); figure imshow(F,[],“XData”,fangles,'ydata',fpos,...'InitialMagnification','fit') axisnormalXlabel('Rotation Angles (degrees)')ylabel(“传感器位置(度)”)Colormap(GCA,热),配色栏

图包含一个轴对象。The axes object contains an object of type image.

Compute fan-beam projections for 'arc' geometry.

I = ones(100); D = 200; dtheta = 45; [Farc,FposArcDeg,Fangles] = fanbeam(I,D,...'fansensorgemetry','arc',...“狂热”,dtheta);

Convert angular positions to linear distance along x-prime axis.

fposarc = d*tan(fposarcdeg*pi/180);

计算“线”几何形状的风扇束投影。

[Fline,FposLine] = fanbeam(I,D,...'fansensorgemetry','line',...“狂热”,dtheta);

Compute the corresponding Radon transform.

[R,Rpos]=radon(I,Fangles);

说play the three projections at one particular rotation angle. Note the three are very similar. Differences are due to the geometry of the sampling, and the numerical approximations used in the calculations.

图IDX = find(fangles == 45);绘图(rpos,r(:,idx),...fposarc,farc(:,idx),...fpoSline,fline(:,idx))传奇('Radon','Arc','Line')

图包含一个轴对象。轴对象包含3个类型行的对象。这些对象代表ra,弧,线。

Input Arguments

全部收缩

输入图像,指定为2-D数字或逻辑矩阵。

从风扇梁顶点到旋转中心的像素中的距离,指定为正数。旋转中心是图像的中心像素,定义为地板((尺寸(i)+1)/2).Dmust be large enough to ensure that the fan-beam vertex is outside of the image at all rotation angles. See提示有关指定指南D. The figure illustratesD关于一个风扇梁几何形状的风扇梁顶点。

Name-Value Arguments

Specify optional comma-separated pairs of名称,价值参数。Name是the argument name and价值是相应的值。Name必须出现在引号中。您可以按任何顺序指定几个名称和值对参数NAME1,Value1,...,Namen,Valuen.

例子:f = fanbeam(i,d,'fanrotation increment',5)

风扇梁旋转角度增量in degrees, specified as the comma-separated pair consisting of“狂热”和一个正标量。

Data Types:双倍的

Fan-beam sensor positioning, specified as the comma-separated pair consisting of'fansensorgemetry'and one of the following values.

价值

Meaning

图表

'arc'

Sensors are spaced at equal angles along a circular arc at distanceDfrom the center of rotation.

粉丝座上定义成角度的程度。

'line'

传感器沿着与平行的线相等的距离间隔x'轴。The closest sensor is distanceDfrom the center of rotation.

粉丝座上defines the distance between fan-beams on thex'axis, in pixels.

风扇蜂向传感器间距,指定为逗号分隔对“粉丝们”和一个正标量。

  • 如果FanSensorGeometry'arc', then粉丝座上定义成角度的程度。

  • 如果FanSensorGeometry'line', then粉丝座上defines the linear distance between fan-beams, in pixels. Linear spacing is measured on thex'轴。

Data Types:双倍的

Output Arguments

全部收缩

风扇梁投影数据,返回numsensors-经过-numangles数字矩阵.numsensors是the number of fan-beam sensors andnumangles是风扇梁旋转角的数量。每一列F以一个旋转角度包含风扇梁传感器样品。

Data Types:双倍的

风扇梁传感器的位置,返回numsensors-b-1数字向量。

  • 如果FanSensorGeometry'arc'(the default), thenfan_sensor_positionscontains the fan-beam spread angles.

  • 如果FanSensorGeometry'line', thenfan_sensor_positions包含沿着风扇梁的传感器位置x'轴。See粉丝座上for more information.

fanbeamdetermines the number of sensors by calculating how many beams are required to cover the entire image for any rotation angle. Fewer sensors are required to cover the image when the distanceDbetween the fan-beam vertex and the center of rotation is large.

Data Types:双倍的

风扇梁传感器的旋转角,作为1乘返回numangles数字向量.numangles360/狂热促进.

Data Types:双倍的

提示

作为指南,尝试制作Da few pixels larger than half the image diagonal dimension, calculated as follows.

sqrt(size(I,1)^2 + size(I,2)^2)

The values returned inF是风扇梁投影的数值近似。该算法取决于ra的变换,并插值为风扇梁几何形状。结果取决于所使用的参数。当图像更大时,您可以期望更准确的结果D是larger, and for points closer to the middle of the image, away from the edges.

参考

[1] Kak,A.C。和Slaney,M.,Principles of Computerized Tomographic Imaging,IEEE出版社,纽约,1988年,第92-93页。

在R2006a之前引入