文档

拉东变换

语法

氡(I)
R =氡(I,theta)
[R,xp] =氡(___
[gpuarrayR,gpuarrayXp] = radon(gpuarrayI,theta)

描述

R=氡(返回Radon变换R灰度图像的对于[0,179]度范围内的角度。Radon变换是图像强度沿某一特定角度的径向线的投影。

例子

R=氡(θ返回指定角度的拉东变换θ

Rxp] =氡(___返回一个向量xp包含图像的每一行对应的径向坐标。

例子

gpuarrayRgpuarrayXp] =氡(gpuarrayIθ在GPU上执行Radon变换。此语法需要并行计算工具箱™。

例子

全部折叠

使此图像的轴刻度可见。

iptsetpref (“ImshowAxesVisible”“上”

创建一个示例映像。

I = 0 (100,100);I(25:75, 25:75) = 1;

计算拉东变换。

Theta = 0:180;[R,xp] =氡(I, θ);

显示转换。

imshow (R, [],“Xdata”θ,“Ydata”xp,“InitialMagnification”“健康”)包含(‘\θ(度)) ylabel (“x””) colormap(gca,hot), colorbar

让坐标轴的刻度不可见。

iptsetpref (“ImshowAxesVisible”“关闭”

在GPU上计算Radon变换并将其可视化。

iptsetpref (“ImshowAxesVisible”“上”) I = 0 (100,100);I(25:75, 25:75) = 1;Theta = 0:180;[R,xp] = radon(gpuArray(I),theta);imshow (R, [],“Xdata”θ,“Ydata”xp,...“InitialMagnification”“健康”)包含(‘\θ(度)) ylabel (“x””) colormap(gca,hot), colorbar (“ImshowAxesVisible”“关闭”

输入参数

全部折叠

灰度图像,指定为数值矩阵。

数据类型:||int8|int16|int32|int64|uint8|uint16|uint32|uint64|逻辑

投影角度(以度为单位),指定为数值标量或数值向量。

如果输入图像为gpuarrayI,然后θ也可以是gpuArray包含数字标量或向量的。

数据类型:

图形处理器上的灰度图像,指定为gpuArray包含数据类型的映像uint8uint16uint32int8int16,或int32

输出参数

全部折叠

像的拉东变换,返回为下列之一。

  • 如果θ是标量,R是一个数字列向量,包含θ度。

  • 如果θ是一个向量,R是一个矩阵,其中每一列都是其中一个角的拉东变换θ

对应于每一行的径向坐标R,作为数字向量返回。径向坐标是沿x-轴,它的方向是θ逆时针方向x设在。两个轴的原点都是图像的中心像素,定义为

地板((大小(I) + 1) / 2)
例如,在一个20 × 30的图像中,中心像素是(10,15)。

像的拉东变换gpuarrayI当在GPU上运行时,返回为gpuArray

对应于每一行的径向坐标gpuarrayR当在GPU上运行时,返回为gpuArray

算法

一个图像的Radon变换是每个像素的Radon变换的和。

算法首先将图像中的像素划分为4个子像素,并将每个子像素分别投影,如下图所示。

每个亚像素的贡献是按比例分成两个最近的箱子,根据投影位置和箱子中心之间的距离。如果亚像素投影击中一个容器的中心点,坐标轴上的容器将获得该亚像素的全部值,或该像素值的四分之一。如果亚像素投影击中两个箱子之间的边界,则亚像素值在箱子之间平均分配。

参考文献

罗纳德·N·布雷斯韦尔,二维成像,恩格尔伍德悬崖,新泽西州,普伦蒂斯霍尔,1995年,第505-537页。

[2]林宰S,二维信号与图像处理,恩格尔伍德悬崖,新泽西州,普伦蒂斯霍尔,1990年,第42-45页。

R2006a之前介绍

这个话题有用吗?