% REPORT ON HEAT KERNEL SMOOTHING
 
signal= imread('toy-key.tif');
signal=(double(signal)-219)/36; % makes the image into real number
figure; imagesc(signal); colormap('bone'); colorbar %596 x 368

noise= normrnd(0, 2, 596, 368);
f = signal + noise;
figure; imagesc(f); colormap('bone'); caxis( [0 1]); colorbar

F = gaussblur(f,100);
figure; imagesc(F); colormap('bone'); caxis( [0 1]); colorbar