拉东变换

描述

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

你可以选择计算拉东变换使用GPU(需要并行计算工具箱™)。

例子

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

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

例子

全部折叠

使这个图像的轴比例可见。

iptsetpref (“ImshowAxesVisible”“上”

创建一个示例图像。

我= 0 (100100);I(25:75, 25:75) = 1;

计算Radon变换。

θ= 0:180;(R, xp) =氡(我,θ);

显示变换。

imshow (R, [],“Xdata”θ,“Ydata”xp,“InitialMagnification”“健康”)包含(‘\θ(度)) ylabel (“x””甘氨胆酸)colormap(热),colorbar

使轴的比例不可见。

iptsetpref (“ImshowAxesVisible”“关闭”

输入参数

全部折叠

灰度图像,指定为二维数字矩阵。

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

以度数表示的投影角,指定为数字标量或数字向量。

数据类型:

输出参数

全部折叠

图像的拉东变换,以下列方式返回。

  • 如果θ是标量吗R数值列向量是否包含Radon变换θ度。

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

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

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

算法

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

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

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

参考文献

罗纳德·N·布雷斯韦尔二维成像, Englewood Cliffs, NJ, Prentice Hall, 1995,第505-537页。

[2] Lim, Jae S,二维信号和图像处理, Englewood Cliffs, NJ, Prentice Hall, 1990,第42-45页。

扩展功能

之前介绍过的R2006a