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