% Call lpfilter
f=imread('TestPattern.tif');
f=im2double(f);
imshow(f,[ ]);pause
[M N]=size(f);
P=2*M; Q=2*N;
% Define filter
D0=0.03*P;
[H D]=lpfilter('gaussian',P,Q,D0);
g=dftfiltTN(f,H);
figure,imshow(abs(g),[ ]);pause;
close all


