自适应均值滤波

function f = adpmedian（g Smax）
if（Smax<=1）||（Smax/2==round（Smax/2））||（Smax~=round（Smax））
   error（‘SMAX must be an odd integer > 1.‘）
end 

% Initial setup.
f = g;
f（:） = 0;
alreadyProcessed = false（size（g））;

% Begin filtering.
for k = 3:2:Smax
   zmin = ordfilt2（g 1 ones（k k） ‘symmetric‘）;
   zmax = ordfilt2（g k * k ones（k k） ‘symmetric‘）;
   zmed = medfilt2（g [k k] ‘symmetric‘）;
   
   processUsingLevelB = （zmed > zmin） & （zmax > zmed） & ...
       ~alreadyProcessed; 
   zB = （g > zmi