radon变换 matlab 代码,Radon变换及其Matlab代码实现

Radon变换和Hough变换类似，最初是用于检测图像中的直线(例如笔直的街道边沿、房屋的边沿、

笔直的电线等)。

关于Hough变换，可以参考OpenCV中的代码和示例(其实除了Hough Lines还有Hough Circles等等变种)，

此处不再赘述。

关于Radon变换，可以参考wiki或者百科，或者网络上的其他资料介绍。

这里做一个简单的总结。

首先准备一张灰度化的图像，及黑白图像，然后检测图像的边缘(如Canny算子，Sobel算子等等，

主要是灰度梯度方法吧)

假设我们得到一幅图像的边缘图像，然后找出这些边缘中的直线(段)，除了之前的Hough算法之外，

这里使用Radon变换。

Radon变换核心在于“线积分”，针对一幅图像，具体来说就是：

设置一个方向(水平自左至右为0度，垂直自下至上为90度，以此类推)

按照这个方向对图像进行“列向量求和”

以0度为例，就是从图像的第1列开始，计算这一列上所有像素的灰度值之和；如此一直计算到最后一列。

以30度为例，如下图，X'为积分方向，对红色垂线所在列的像素进行灰度值叠加求和，

得到积分如图中的淡蓝色小块(仅作示意，可能不准确)

其他角度类似。

如果图像中有一条亮白的直线段，那么会存在这样的情况，

沿着某个方向积分得到的积分图没有“突变”(平缓，没有特别的尖峰)

沿着与之垂直的聆听一个方向则存在“突变”

下图是一个简单的示例

至此，你应当能体会到Radon变换的原理了。

接下来用Matlab进行测试

%% Radon Transform

function RadonTest()

fileName=input('Input image file name:');

srcImage=imread(fileName);

grayImage=rgb2gray(srcImage);

cannyImage=edge(grayImage,'canny');

theta=0:180;

[R,x]=radon(cannyImage,theta);

figure,imagesc(theta,x,R);

title('R_theta X');

xlabel('theta(degree)');

ylabel('X\prime');

colormap(hot);

colorbar;

end

首先载入一张图片，然后灰度化，接着检测边缘

然后针对边缘原图像进行Radon变换

最后得到结果

示例1

原图像(有一条直线段)

Radon变换结果

其中的极亮(暗)点对应的是 X'=-45, theta=50° (大概位置)，这表明原图像中截距45处，50+90=140°方向有一条直线(段)

示例2

原图像

Radon变换结果

我正在尝试利用OpenCV或者EmguCV实现Radon变换，请关注后续更新。

本文原创，文章链接

http://blog..net/fengyhack/article/details/41873119