光照归一化算法——DoG滤波，自商图

DoG code：
function pic=DoG(I)
if size(I,3)==3
    I=rgb2gray(I);
end
I=double(I);
h1=fspecial('gaussian',15,1);
h2=fspecial('gaussian',15,2);
m=filter2(h2,I)-filter2(h1,I);
R=m;
mi=min(min(R));
ma=max(max(R));
newR=(R-mi)*255/(ma-mi);
pic=uint8(newR);

处理效果：

理论：

Given a m-channels, n-dimensional image