主要内容

拉东变换

请注意

有关沿单一源辐射路径(称为扇束投影)的线积分创建投影数据的信息,请参见扇形波束投影.若要将平行光束投影数据转换为扇形光束投影数据,请使用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变换

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

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

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

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

计算0°和45°角度下图像的Radon变换。

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

画出0°的变换。

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

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

绘制45度的变换。

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

图中包含一个轴对象。标题为Radon Transform of a Square Function at 45°的轴对象包含一个类型为line的对象。

视Radon变换为图像

大量角度的拉东变换通常显示为图像。在这个例子中,对正方形图像的Radon变换计算角度从0°到180°,以1°增量。

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

使用180投影的Radon变换