主要内容

拉东变换

请注意

有关从沿着从单个源辐射到的路径的线路积分的创建投影数据的信息,称为扇形光束投影,请参阅扇形波束投影.要将并行波束投影数据转换为扇形投影数据,请使用para2fan函数。

函数计算预测图像矩阵沿指定方向的。

二维函数的投影f (x, y)是线积分的集合。的函数计算沿平行路径的多个源的线积分,在某个方向上。光束之间以1个像素单位间隔。要表示图像,则函数通过围绕图像中心旋转源,从不同角度对图像进行多次平行光束投影。下图显示了在指定旋转角度下的单个投影。

旋转角度的平行光束投影

例如,线积分f (x, y)在垂直方向上的投影是f (x, y)x设在;直线积分在水平方向的投影f (x, y)y设在。下图显示了一个简单的二维函数的水平和垂直投影。

一个简单函数的水平和垂直投影

可以沿任何角度计算投影θ(θ)。一般来说,的拉东变换f (x, y)线积分是f平行于y´设在

R θ. x f x 因为 θ. y θ. x θ. + y 因为 θ. d y

在哪里

x y 因为 θ. θ. θ. 因为 θ. x y

下图说明了拉东变换的几何形状。

Radon变换的几何

绘制图像的氡变换

此示例显示如何计算图像的Radon变换,,对于一组特定的角度,θ,使用函数。函数返回,R,其中列包含每个角的拉东变换θ.该函数还返回向量,xp,其中包含沿x轴的相应坐标。中心像素定义为地板((大小(I) + 1) / 2),即x轴上对应的像素x = 0

为本示例创建一个小示例图像,该图像包含一个方形对象并显示它。

我= 0 (100100);我(25:75 25:75)= 1;imshow(我)

图中包含一个坐标轴。轴包含一个image类型的对象。

计算图像的氡变换为角度0度和45度。

[R,xp] =氡(I,[0 45]);

画出0°的变换。

图绘制(xp, R (: 1));标题(' 0度平方函数的Radon变换'

图中包含一个坐标轴。标题为“0度平方函数拉东变换”的轴包含一个类型为line的对象。

绘制变换45度。

图绘制(xp, R (:, 2));标题(“45度平方函数的Radon变换”

图中包含一个坐标轴。标题为“矩形函数的拉东变换45度”的轴包含一个类型为line的对象。

视Radon变换为图像

大量角度的氡变换通常显示为图像。在该示例中,方形图像的氡变换以0°至180°以1°增量计算的。

θ= 0:180;(R, xp) =氡(我,θ);显示亮度图像(θ,xp, R);标题(“R_{\θ}(X \ ') ');包含(“\θ(度)”);ylabel (X \ ');集(gca, XTick, 0:20:180);colormap(热);colorbar

使用180投影的Radon变换