Matlab的hough圆检测-soble_hough.m

Matlab的hough圆检测-soble_hough.m
先用sobel算子对图像进行预处理，提取边缘。对于不连续的圆，在利用hough算法检测圆心坐标和半径，再将圆进行重画，得到完整的圆并显示。

4.jpg 待处理图片    

sobel边缘提取后：

sobel.jpg

圆没有完全闭合。

hough检测并重画圆：
sobel_hough.jpg

a=imread（‘1.jpg‘）;
b=rgb2gray（a）;
in=edge（b‘sobel‘）;

[hw]=size（in）;
if w>h  g=floor（h/2）;
else g=floor（w/2）;
end
B=zeros（ghw）;%创建一个三维数组

%具体算法
 for j=1:w    
    for i=1:h        
        if in（ij）==1
            for aa=1:w
                for bb=1:h
                    rr=floor（sqrt（（j-aa）^2+（i-bb）^2））;                                   
                    if rr                    end
                end