%% adquirir imagen1
I1 = imread('imagen1.jpg');
imshow(I1)
%% adquirir imagen2
I2 = imread('imagen2.jpg');
figure
imshow(I2)
%% capturar las matrices de coleres
I1r = I1(:,:,1);
I2r = I2(:,:,1);
I1g = I1(:,:,2);
I2g = I2(:,:,2);
I1b = I1(:,:,3);
I2b = I2(:,:,3);
Inewr = I1r-I2r;
figure
imshow(Inewr)
Inewg = I1g-I2g;
figure
imshow(Inewg)
Inewb = I1b-I2b;
figure
imshow(Inewb)
Inew(:,:,1) = Inewr;
Inew(:,:,2) = Inewg;
Inew(:,:,3) = Inewb;
figure
imshow(Inew)
%% Filtro mediana
[n,m]=size(I2r);
imshow('imagen2.jpg')
for k =1:3
for i = 2:n-1
for j = 2:m-1
v=uint8((I2(i-1,j-1,k)+ I2(i-1,j,k)+ I2(i-1,j+1,k)+ I2(i,j-1,k)+ I2(i,j,k)+ I2(i,j+1,k)+...
I2(i+1,j-1,k)+ I2(i+1,j,k)+ I2(i+1,j+1,k))/9);
I2(i,j,k) = (v);
end
end
end
figure
imshow(I2)
%% omparar imagen
for k = 1:3
Inew(:,:,k) = I1(:,:,k)-I2(:,:,k);
end
figure
imshow(Inew)
imshow(Inew)
