function H=hpfilter(type,M,N,D0,n )
%HPFilter
[U V]=dftuvTN(M,N);
D=sqrt(U.^2+V.^2);
if nargin ==4
   n=1;
end
Hlp=lpfilter(type,M,N,D0,n);
H=1-Hlp;