matlablunkuotiqu.txt

图象轮廓提取 比较简单 是本人在MATLAB环境下编译的

img=imread('lena.bmp');
imshow(img);
[m,n]=size(img);
img=double(img);
for i=2:m
  for  j=2:n
   img_pingyi1(i-1,j)=img(i,j) ; img_pingyi2(i+1,j)=img(i,j) ;  img_pingyi3(i,j-1)=img(i,j)  ; img_pingyi4(i,j+1)=img(i,j) ;
   img_pingyi5(i-1,j-1)=img(i,j) ; img_pingyi6(i+1,j+1)=img(i,j) ;  img_pingyi7(i-1,j+1)=img(i,j) ; img_pingyi8(i+1,j-1)=img(i,j) ;
end
end
[img_pingyi1_h,img_pingyi1_k]=size(img_pingyi1);
[img_pingyi2_h,img_pingyi2_k]=size(img_pingyi2);
[img_pingyi3_h,img_pingyi3_k]=size(img_pingyi3);
[img_pingyi4_h,img_pingyi4_k]=size(img_pingyi4);
[img_pingyi5_h,img_pingyi5_k]=size(img_pingyi5);
[img_pingyi6_h,img_pingyi6_k]=size(img_pingyi6);
[img_pingyi7_h,img_pingyi7_k]=size(img_pingyi7);
[img_pingyi8_h,img_pingyi8_k]=size(img_pingyi8);
h1=min(img_pingyi2_h,img_pingyi1_h);k1=min(img_pingyi2_k,img_pingyi1_k);
h2=min(img_pingyi4_h,img_pingyi3_h);k2=min(img_pingyi4_k,img_pingyi3_k);
h3=min(img_pingyi5_h,img_pingyi6_h);k3=min(img_pingyi5_k,img_pingyi6_k);
h4=min(img_pingyi7_h,img_pingyi8_h);k4=min(img_pingyi7_k,img_pingyi8_k);
for i=1:h1
  for  j=1:k1
  imgj1(i,j)=img_pingyi2(i,j)-img_pingyi1(i,j);imgi1(i,j)=img_pingyi1(i,j)-img_pingyi2(i,j);
end
end
for i=1:h2
  for  j=1:k2
  imgj2(i,j)=img_pingyi4(i,j)-img_pingyi3(i,j); imgi2(i,j)=img_pingyi3(i,j)-img_pingyi4(i,j);
end
end
for i=1:h3
  for  j=1:k3
 imgj3(i,j)=img_pingyi5(i,j)-img_pingyi6(i,j);imgi3(i,j)=img_pingyi6(i,j)-img_pingyi5(i,j);
end
end
for i=1:h4
  for  j=1:k4
  imgj4(i,j)=img_pingyi7(i,j)-img_pingyi8(i,j);imgi4(i,j)=img_pingyi8(i,j)-img_pingyi7(i,j);
end
end

figure;
imshow(imgj1,[0,255]);
figure;
imshow(imgj2,[0,255]);
figure;
imshow(imgj3,[0,255]);
figure;
imshow(imgj4,[0,255]);
figure;
imshow(imgi1,[0,255]);
figure;
imshow(imgi2,[0,255]);
figure;
imshow(imgi3,[0,255]);
figure;
imshow(imgi4,[0,255]);
[a,b]=size(imgj1);[c,d]=size(imgj2);[e,f]=size(imgj3);[g,l]=size(imgj4);[a1,b1]=size(imgi1);[c1,d1]=size(imgi2);[e1,f1]=size(imgi3);[g1,l1]=size(imgi4);
a2=min(a,c);e2=min(e,g);a3=min(a1,c1);e3=min(e1,g1);
a4=min(a2,e2);a5=min(a3,e3);
u=min(a4,a5);
b2=min(b,d);f2=min(f,l);b3=min(b1,d1);f3=min(f1,l1);
b4=min(b2,f2);b5=min(b3,f3);
v=min(b4,b5);

%取阈值
[m1,n1]=size(img_pingjun2);
for i=1:m1
  for  j=1:n1
    if   img_pingjun2(i,j)<=200
      img_pingjun3(i,j)=0;
    else
      img_pingjun3(i,j)=255;
end
end
end
 figure;
imshow(img_pingjun3,[0,255]);