load finalimage2;
% blurho is blurred image of house
% refim is to be used only for calculation of MSEs and comparison

bhwin = blurho.*(ones([480 1])*hanning(512)');
spec = sum(abs(fft(bhwin,[],2)),1);
subplot(211); semilogy(spec); title('Average spatial spectrum')
subplot(212); plot(blurho(125,:)); title('Plot line 125')

mseblur = mean((refim(:) - blurho(:)).^2)
