a = imread('PSWebImageBitMap0001.bmp',  'bmp');
b = double(a(:,:,2));

f = inline('fLOG(x)');
c = nlfilter(b,[5 5],f);

figure(4)
clf
imagesc(c)
colormap(gray)
title('Laplacian of gradient (LOG) mask')
