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


