Reporte de actividades del proyecto de Investigación en ...

im=double(im);
[ancho largo]=size(im)
for x=2:1:ancho-1
for y=2:1:largo-1
a=im(x-1,y-1);b=im(x,y-1);c=im(x+1,y-1);d=im(x-1,y);e=im(x,y);
f=im(x+1,y); g=im(x-1,y+1);h=im(x,y+1);i=im(x+1,y+1);
a=double(a);b=double(b);c=double(c);d=double(d);e=double(e);
f=double(f);g=double(g);h=double(h);i=double(i);
Roberts(x,y)=abs(e-a)+abs(b-d);
SobelM1=a*(-1)+d*(-2)+g*(-1)+c+f*2+i;
SobelM2=a*(-1)+b*(-2)+c*(-1)+g+h*2+i;
valores(1:2)=[SobelM1 SobelM2];
Sobel(x,y)=uint8(max(valores));
PrewittM1=a*(-1)+b*(-1)+c*(-1)+g+h+i;
PrewittM2=c*(-1)+f*(-1)+i*(-1)+a+d+g;
valores(1:2)=[PrewittM1 PrewittM2];
Prewitt(x,y)=uint8(max(valores));
KirchM1=a*(+5)+b*(+5)+c*(+5)+d*(-3)+f*(-3)+g*(-3)+h*(-3)+i*(-3);
KirchM2=a*(-3)+b*(+5)+c*(+5)+d*(-3)+f*(+5)+g*(-3)+h*(-3)+i*(-3);
KirchM3=a*(-3)+b*(-3)+c*(+5)+d*(-3)+f*(+5)+g*(-3)+h*(-3)+i*(+5);
KirchM4=a*(-3)+b*(-3)+c*(-3)+d*(-3)+f*(+5)+g*(-3)+h*(+5)+i*(+5);
KirchM5=a*(-3)+b*(-3)+c*(-3)+d*(-3)+f*(-3)+g*(+5)+h*(+5)+i*(+5);
KirchM6=a*(-3)+b*(-3)+c*(-3)+d*(+5)+f*(-3)+g*(+5)+h*(+5)+i*(-3);
KirchM7=a*(+5)+b*(-3)+c*(-3)+d*(+5)+f*(-3)+g*(+5)+h*(-3)+i*(-3);
KirchM8=a*(+5)+b*(+5)+c*(-3)+d*(+5)+f*(-3)+g*(-3)+h*(-3)+i*(-3);
valores(1:8)=[KirchM1 KirchM2 KirchM3 KirchM4 KirchM5 KirchM6 KirchM7
KirchM8];
Kirch(x,y)=uint8(max(valores));
LaplacianoM1(x,y)=b*(-1)+d*(-1)+e*(4)+f*(-1)+h*(-1);
LaplacianoM2(x,y)=a+b*(-2)+c+d*(-2)+e*(4)+f*(-2)+g+h*(-2)+i;
LaplacianoM3(x,y)=a*(-1)+b*(-1)+c*(-1)+d*(-1)+e*(8)+f*(-1)+g*(-1)+h*(-1)+i*(-1);
end
end
figure, imshow(im) ; title('Original');
figure, imshow(Roberts);title('Roberts');
figure, imshow(Sobel);title('Sobel');
figure, imshow(Prewitt);title('Prewitt');
figure, imshow(Kirch);title('Kirch');
figure, imshow(LaplacianoM1);title('Laplaciano a');
figure, imshow(LaplacianoM2);title('Laplaciano b');
figure, imshow(LaplacianoM3);title('Laplaciano c');
Programa #7: Detección de contornos
96