function?II=adp_median（ISmax）
%自適應中值濾波

if?（Smax<=1）||（Smax/2==round（Smax/2））||（Smax~=round（Smax））%Smax約束，噪聲密度越大，Smax應當設置越高
????error（‘Smax?must?be?an?odd?integer?>1‘）;
end
%初始化
II=I;
II（:）=0;
alreadyProcessed=false（size（I））;
%迭代
for?k=3:2:Smax
zmin=ordfilt2（I1ones（kk）‘symmetric‘）;
zmax=ordfilt2（Ik*kones（kk）‘symmetric‘）;
zmed=medfilt2（I[k?k]‘symmetric‘）;
processUsingLevelB=（zmed>zmin）&（zm