% Gaussian filtering
f=imread('Figurer.bmp');
%f=imread('building.tif'); 
f=im2double(f);
imshow(f,[ ]);pause
[M N]=size(f);
P=2*M; Q=2*N;
% Define filter
h=fspecial('gaussian',7,3)      % 7x7 with standard deviation 3
H=freqz2(h,P,Q);
H1=ifftshift(H);                % Shift to non-senter version
g=dftfiltTN(f,H1);
figure,imshow(abs(g),[ ]);pause;
close all


