function?f?=adpmedian（gSmax）
%ADPMEDIAN?Perform?adap?tive?median?filtering
%首先判斷Smax是否符合要求:?Smax應(yīng)該是大于1的奇數(shù)
if?（Smax?<=?1）|（Smax/2?==?round（Smax/2））|（Smax?~=?round（Smax））
????error（‘Smax?must?be?an?odd?integer>1.‘）
end
[MN]=size（g）;
%初始化.
f=g;
f（:）=0;
alreadyProcessed=false（size（g））;
%開始濾波.
for?k=3:2:Smax
%相當(dāng)于k*k的極小值濾波
?????zmin=ordfilt2（g1ones（kk））;
%相當(dāng)于k*k的極大值濾波
?????zmax=ordfilt2（gk*kones（kk））;
%?k*k的中值濾波
?????zmed=medfilt2（g[k?k]）;
%判斷Zmed是否為一脈沖
?????processUsingLevelB=（zmed>zmin）?&?（zmax>zmed）?&?~alreadyProcessed;
%判斷Zxy是否為一脈沖
?????zB=（g>zmin）?&?（zmax>g）;
?????outputZxy=processUsingLevelB?&?zB;
?????outputZmed=processUsingLevelB?&?~zB;
%輸出一個不變的像素值Zxy來代替鄰域中值作為輸出
?????f（outputZxy）=g（outputZxy）;
%輸出中值濾波的值作為輸出
?????f（outputZmed）=g（outputZmed）;
?????a

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????9388??2014-05-28?17:20??untitled.m

?????文件?????141724??2014-05-27?22:50??1.jpg

?????文件???????1016??2014-05-28?16:46??adpmedian.m

?????文件???????5399??2014-05-28?17:18??untitled.fig

-----------?---------??----------?-----??----

???????????????157527????????????????????4