function af = matmedfilt( a, m )
% MATMEDFILT - Moving median filter
% af = medfilt( x, m )

if( nargin ~=2 ), m=5;, end;
if(mod(m,2)==0),error('m must be odd'),end
m2 = fix(m/2);
[nr,nc] = size(a);
af = a;

for j=1:nc,
  for i=(m2+1):(nr-m2-1),
    af(i,j) = nanmedian( a((i-m2:i+m2),j) );
  end
end