map1=imread('BMW.jpg');
[m,n,l]=size(map1); %行,列,深度值
I=zeros(m,n);
%I=uint8(I);
%灰度化
for i=1:m
for j=1:n
I(i,j)=0.11*map1(i,j,1)+0.59*map1(i,j,2)+0.3*map1(i,j,3);
end
end
%灰度化完成,输出灰度图
figure(1)
imshow(I,[])
title('原始灰度图像')
B = zeros(m,n);
A = I;
for i=1:m-1;
for j=1:n-1;
B(i,j)=max(abs(A(i+1,j+1)-A(i,j)),abs(A(i+1,j)-A(i,j+1)));
end
end
figure(4)
imshow(B,[])
title('Robert算法处理后')