sigma = 5;
sz = 25;    % length of gaussFilter vector
x = linspace(-sz / 2, sz / 2, sz);
g = exp(-x .^ 2 / (2 * sigma ^ 2));
g2 = g' * g;


m = load('/mnt/flashblade/david/cache/hist_fine.txt');
m = conv2(m, g2, 'same');
% m = sqrt(m);
m /= sum(sum(m));
m /= max(max(m));
size(m)

imshow(m)
colormap(jet(256))
