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);
log=[-2 -4 -4 -4 -2;...
    -4 0 8 0 -4;...
    -4 8 24 8 -4;...
    -4 0 8 0 -4;...
    -2 -4 -4 -4 -2];
 A = I;
 for i = 3:m-2
     for j = 3:n-2
         S=sum(sum(A(i-2:i+2,j-2:j+2).*log));
         B(i,j) = S/56;
     end
 end
 figure(2)
 imshow(B,[])
 title('LOG算法处理后')