图像算法之一：Randon变换

1.原理： Radon变换用来计算图像矩阵在特定方向上的投影。二维函数投影是一组线积分，Radon变换计算一定方向上平行线的积分，平行线的间隔为1个像素。Radon变换可以旋转图像的中心到不同角度，来获得图像在不同方向上的投影积分。

下图是矩形图像在[0,180]度区间的Radon变换：

2、应用： 使用Radon变换检测直线

检测步骤如下：
（1）使用边缘检测函数edge函数计算二值图像；
（2）计算二值图像的Radon变换；
（3）寻找Radon变换的局部极大值，这些极大值的位置即为原始图像中直线的位置。

BW = edge(I1);
imshow(BW)
theta = 0:179;
[R,xp] = radon(BW,theta);
figure,imagesc(theta,xp,R)

3、逆Randon变换
Randon变换的逆变换，可以用于重建图像。Matlab调用语句为： I = iradon(R,theta);