clear;load imagefile.mat;sigma=5;L=25;
L1=(L+1)/2;for m=1:L;for n=1:L;
G(m,n)=exp(-((m-L1)^2+(n-L1)^2)/2/sigma^2)/sqrt(2*pi*sigma^2);
end;end;G=G/sum(sum(G));
Y=X+X-conv2(X,G,'same');
figure,imagesc(X),colormap(gray),axis off
figure,imagesc(Y),colormap(gray),axis off
