主要内容

伊拉登

逆氡变换

描述

例子

=伊拉顿(R,西塔)重建图像从中的投影数据R以投影角度捕捉西塔.

=伊拉顿(R,西塔,内页,筛选,频率标度,输出大小)指定要在Radon逆变换中使用的参数。可以指定最后四个参数的任意组合。伊拉登对省略的参数使用默认值。

[,H] = Iradon(___)还返回滤波器的频率响应,H.

例子

全部收缩

创建幻影的图像。显示图像。

p = phantom(128);imshow(p)标题('原始图像​​')

图中包含一个轴对象。标题为“原始图像”的轴对象包含一个类型为“图像”的对象。

执行图像的氡变换。

R=氡(P,0:179);

执行过滤反投影。

I1=iradon(R,0:179);

执行未经过滤的反调。

I2=iradon(R,0:179,“线性”,“没有”);

显示重建图像。

图形子图(1,2,1)IMSHOW(I1,[])标题('过滤反冲')子图(1,2,2)IMSHOW(I2,[])标题('未经过滤的背部重点')

图中包含2个轴对象。带有标题过滤反投影的轴对象1包含图像类型的对象。带有标题未过滤反投影的轴对象2包含图像类型的对象。

创建幻影的图像。

p = phantom(128);

执行图像的氡变换,然后以45度角获得对应于投影的投影载体。

R=氡(P,0:179);R45 = R(:,46);

执行该单个投影向量的逆氡变换。这伊拉登语法不允许您直接执行此操作,因为如果西塔是一个标量,它被视为增量。您可以通过传入两个投影向量副本,然后将结果除以2来完成此任务。

i = Iradon([R45 R45],[45 45])/ 2;

显示结果。

imshow(I,[])标题(“从45度反向投影”)

图包含轴对象。具有从45度的标题背部注射的轴对象包含类型图像的对象。

输入参数

全部收缩

并行波束投影数据,指定为以下之一。

  • 如果西塔是标量,然后指定R作为数值列向量,包含西塔度。

  • 如果西塔是一个向量,然后指定R作为一个二维矩阵,其中每列都是其中一个角度的Radon变换西塔.

数据类型:单一的|双倍的|int8|int16|INT32.|INT64.|uint8|uint16|uint32|uint64|符合逻辑的

投影角度(以度为单位),指定为以下之一。

价值 描述
数字矢量 投影角度。两个角之间的间距必须相等。
数字标量 投影之间的增量角度。投影在角度上拍摄m *θ., 在哪里m=0,1,2,…,尺寸(R,2)-1.
[] 自动在投影之间设置增量角度180/尺寸(R,2)

数据类型:双倍的

反投影中使用的插值类型,指定为这些值之一,按精度和计算复杂度的增加顺序列出。

价值

描述

'最近'

最近邻的插值

“线性”

线性插值

'样条曲线'

样条插值

“pchip” 保形分段三次插值
'v5cbic' 三次卷积在MATLAB中的应用®5.

数据类型:char|一串

过滤器用于频域过滤,指定为其中一个值。

价值

描述

'Ram-Lak'

裁剪Ram Lak或斜坡滤波器。该滤波器的频率响应为|F|。因为该过滤器对投影中的噪声敏感,所以下面列出的一个过滤器可能是优选的。这些过滤器将RAM-LAK过滤器乘以解除高频的窗口。

“谢普·洛根”

将RAM-LAK过滤器乘以a信诺功能

'余弦'

将RAM-LAK过滤器乘以a余辉功能

“哈明”

将Ram Lak过滤器乘以汉明窗口

“汉恩”

将Ram Lak过滤器乘以Hann窗口

'没有任何' 没有过滤。伊拉登返回未过滤的反投影数据。

数据类型:char|一串

用于重新缩放频率轴的比例因子,指定为范围(0,1)中的正数。如果频率标度小于1,则滤波器被压缩以适应频率范围[0,频率缩放],归一化频率;所有频率上面频率标度准备好了吗0.

重建图像中的行数和列数,指定为正整数。如果输出大小如果未指定,则根据以下内容根据投影长度确定尺寸:

output_size = 2 *楼层(大小(R,1)/(2*sqrt(2)))

如果您指定输出大小那么伊拉登重建图像的较小或较大部分,但不改变数据的缩放比例。如果使用功能,然后重建的图像可能与原始图像的大小不同。

输出参数

全部收缩

灰度图像,作为数字矩阵返回。如果输入投影数据R是数据类型单一的那么单一的;除此以外双倍的.

数据类型:单一的|双倍的

过滤器的频率响应,作为数字矢量返回。

数据类型:双倍的

算法

伊拉登假设旋转中心是投影的中心点,其被定义为天花板(尺寸(R,1)/2).

伊拉登使用滤波反投影算法执行逆Radon变换。滤波器直接在频域中设计,然后乘以投影的FFT。滤波前,投影被零填充到2的幂,以防止空域混叠并加快FFT。

参考

[1] Kak,A.C.和M.Slaney,计算机断层成像原理,纽约,纽约,IEEE出版社,1988年。

扩展能力

在R2006A之前介绍