hough-circle 读入图像

读入图像，对图像进行预处理，随后检测边缘，对得到的图像进行二值化，然后进行hough变换，获得圆边界，将此图与原图叠加，增强边缘

clear all;
%1.读入图像并显示
S=imread（‘E:\Sweden-coins.bmp‘）;
figure;imshow（S）;title（‘原图像‘）;

%2.引入均值为零，方差为a的高斯噪声对图像进行平滑，然后进行拉普拉斯变换
[mn]=size（S）;
a=1;
h=zeros（mn）;
h=double（h）;
for i=1:m
    for j=1:n
        t=i^2+j^2;
        h（ij）=（t-2*a）*exp（-（t/（2*a）））/a^2;
    end
end
H=fftshift（fft2（h））;
G=fftshift（fft2（S））;
W=H.*G;
w=ifft2（ifftshift（W））;
figure;imshow（w）;title（‘平滑后进行拉普拉斯变换后的图像‘）;

%3.对水平，垂直，+45度，-45度方向进行边缘检测，采用Prewitt梯度算子
IMG=double（w）;
H1=zeros（mn）;
H2=zeros（mn）;
G1=zeros（mn）;
G2=zeros（mn）;
for i=2:m-1
    for j=2:n-1
        H1（ij）=-IMG（i-1j-1）-IMG（i-1j）-IMG（i-1j+1）+IMG（i+1j-1）+IMG（i+1j）+IMG（i+1j+1）;
        H2（ij）=-IMG（i-1j-1）-IMG（ij-1）-IMG（i+1j-1）+IMG（i-1j+1）+IMG（ij+1）+IMG（i+1j+1）;
        G1（ij）=IMG（i-1j）+IMG（i-1j+1）+IMG（ij+1）-

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

     文件       3810  2012-06-05 11:42  王小龙（1101120432）+张芝英（1101120434）+刘炼（1101120435）大作业\myself.m

     文件     718060  2012-06-07 16:18  王小龙（1101120432）+张芝英（1101120434）+刘炼（1101120435）大作业\大作业报告.doc

     文件    1503176  2012-06-05 11:37  王小龙（1101120432）+张芝英（1101120434）+刘炼（1101120435）大作业\实验结果\+45度边缘图像.fig

     文件    1499582  2012-06-05 11:38  王小龙（1101120432）+张芝英（1101120434）+刘炼（1101120435）大作业\实验结果\-45度边缘图像.fig

     文件       7700  2012-06-05 11:34  王小龙（1101120432）+张芝英（1101120434）+刘炼（1101120435）大作业\实验结果\hough变化后的图像.fig

     文件      11402  2012-06-05 11:35  王小龙（1101120432）+张芝英（1101120434）+刘炼（1101120435）大作业\实验结果\二值化后的图像.fig

     文件     132512  2012-06-05 11:41  王小龙（1101120432）+张芝英（1101120434）+刘炼（1101120435）大作业\实验结果\原图像.fig

     文件    1500067  2012-06-05 11:39  王小龙（1101120432）+张芝英（1101120434）+刘炼（1101120435）大作业\实验结果\垂直边缘图像.fig

     文件    1467194  2012-06-05 11:41  王小龙（1101120432）+张芝英（1101120434）+刘炼（1101120435）大作业\实验结果\平滑后进行拉普拉斯变换后的图像.fig

     文件    1514450  2012-06-05 11:36  王小龙（1101120432）+张芝英（1101120434）+刘炼（1101120435）大作业\实验结果\提取的边缘图像.fig

     文件     131751  2012-06-05 11:31  王小龙（1101120432）+张芝英（1101120434）+刘炼（1101120435）大作业\实验结果\最终结果.fig

     文件    1503084  2012-06-05 11:40  王小龙（1101120432）+张芝英（1101120434）+刘炼（1101120435）大作业\实验结果\水平边缘图像.fig

     目录          0  2012-06-05 11:41  王小龙（1101120432）+张芝英（1101120434）+刘炼（1101120435）大作业\实验结果

     目录          0  2012-06-07 16:19  王小龙（1101120432）+张芝英（1101120434）+刘炼（1101120435）大作业

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

              9992788                    14