function w=Ahp(A) [M N K]=size(A); for m=1:M for n=1:N w(m,n,:)=highpass(A(m,n,:)); end end end