数字图像处理图像分割matlab算法代码

本次试验主要是对图像进行分割处理，通过算法检测图像内容的边缘，并进行提取

x=imread（‘hough.bmp‘）;
subplot（221）imshow（x）title（‘原图象‘）;
x=double（x）;
[mn]=size（x）;
       %       zerocross检测图象       %
X=edge（x‘zerocross‘）;  
subplot（222）imshow（X）title（‘zerocross检测图象‘）;
       %          极坐标参数域图象    %
X=double（X）;
M=499;N=456;
H=zeros（MN）;
i1=1;
for i=1:m
    for j=1:n
        if X（ij）==1
            for t=2*pi/N:2*pi/N:2*pi
                k=i*cos（t）+j*sin（t）;
                H（round（k）+250round（t/2/pi*N））=H（round（k）+250round（t/2/pi*N））+1;
            end
        end
    end
end
c=max（max（H））;
H1=H/c*255;
H1=uint8（H1）;
subplot（223）imshow（H1）title（‘极坐标参数图象‘）;
     %      直线图象         %
d=0;
H1=double（H1）;
for i=1:499
    for j=1:256
        if H1（ij）>=120
            d=d+1;
        end
    end
end
B=zeros（2*d1）;
for i=1:499
    for j=1:256
        if H1（ij）>=120
            B（i11）=i;
            B（i1+11）=j;
            i1=i1+2;
        end
    end
end
x1=double（x）;
for l=1:2:d
    k=B（l1）-250;
    sita=B（l+11）*2*pi/256;
    x1=double（x1）;
    for i=1:m
        j=round（（k-i*cos（sita））/sin（sita））;
        if j>=1 & j<=n
            x1（ij）=255;
        end
    end
end
x1=uint8（x1）;
subplot（224）imshow（x1）title（‘直线图象‘）;

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

     文件      43558  2009-06-06 23:45  实验六\1.bmp

     文件      26464  2009-06-07 23:59  实验六\hough.bmp

     文件       1270  2009-06-08 00:28  实验六\hough.m

     文件       2865  2009-06-07 22:33  实验六\Laplacian.m

     文件       2824  2009-06-07 22:22  实验六\Sobel.m

     文件     155136  2009-06-08 21:52  实验六\实验报告.doc

     文件        128  2009-06-07 23:21  实验六\实验说明.txt

     目录          0  2009-07-01 20:01  实验六

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

               232245                    8