hough变换的matlab源代码

Hough变换可以检测已知形状的目标，而且受噪声和曲线间断的影响小，matlab代码。

function [hough_spacehough_circlepara] = hough_circle（BWstep_rstep_angler_minr_maxp）

% %%%%%%%%%%%%%%%%%%%%%%%%%%
% input
% BW:二值图像；
% step_r:检测的圆半径步长
% step_angle:角度步长，单位为弧度
% r_min:最小圆半径
% r_max:最大圆半径
% p:阈值，0，1之间的数

% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% output
% hough_space:参数空间，h（abr）表示圆心在（ab）半径为r的圆上的点数
% hough_circl:二值图像，检测到的圆
% para:检测到的圆的圆心、半径

[mn] = size（BW）;
size_r = round（（r_max-r_min）/step_r）+1;
size_angle = round（2*pi/step_angle）;

hough_space = zeros（mnsize_r）;

[rowscols] = find（BW）;
ecount = size（rows）;

% Hough变换
% 将图像空间（xy）对应到参数空间（abr）
% a = x-r*cos（angle）
% b = y-r*sin（angle）
for i=1:ecount
    for r=1:size_r
        for k=1:size_angle
            a = round（rows（i）-（r_min+（r-1）*step_r）*cos（k*step_angle））;
            b = round（cols（i）-（r_min+（r-1）*step_r）*sin（k*step_angle））;
            if（a>0&a<=m&b>0&b<=n）
                hough_space（abr） = hough_space（abr）+1;
            end
        end
    end
end

% 搜索超过阈值的聚集点
max_para = max（max（max（hough_space）））;
index = find（hough_space>=max_para*p）;
length = size（index）;
hough_circle = false（mn）;
for i=1:ecount
    for k=1:length
        par3 = floor（index（k）/（m*n））+1;
        par2 = floor（（index（k）-（par3-1）*（m*n））/m）+1;
        par1 = index（k）-（par3-1）*（m*n）-（par2-1）*m;
        if（（rows（i）-par1）^2+（cols（i）-par2）^2<（r_min+（par3-1）*step_r）^2+5&...
                （rows（i）-par1）^2+（cols（i）-par2）^2>（r_min+（par3-1）*step_r）^2-5）
            hough_circle（rows（i）cols（i）） = true;
        end
    end
end

% 打印检测结果
for k=1:length
    par3 = floor（index（k）/（m*n））+1;
    par2 = floor（（index（k）-（par3-1）*（m*n））/m）+1;
    par1 = index（k）-（par3-1）*（m*n）-（par2-1）*m;
    par3 = r_min+（par3-1）*step_r;
    fprintf（1‘Center %d %d radius %d\n‘par1par2par3）;
    para（:k） = [par1par2par3];
end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件        406  2005-12-17 12:39  main.m

     文件      36352  2005-12-17 13:37  Hough变换.doc

     文件      46926  2005-12-16 11:34  2.bmp

     文件       1973  2005-12-17 13:33  hough_circle.m

*    文件       4862  2006-09-29 09:08  中国图象图形网下载说明.html

----------- ---------  ---------- -----  ----

                90519                    5