function?H?=?adMedFilt（f）
????[M?N]?=?size（f）;
????flag?=?zeros（MN）;?%?判斷像素點是否處理，0未處理，1處理
????for?k?=?3:2:min（[M?N]）?%?自動選取模板大小
????????fmin?=?ordfilt2（f?1?ones（k?k）?‘symmetric‘）;?%?對稱擴展邊界，輸出當前窗口的灰度最小值
????????fmax?=?ordfilt2（f?k*k?ones（k?k）?‘symmetric‘）;?%?輸出當前窗口的灰度最大值
????????fmed?=?ordfilt2（f?（k*k+1）/2?ones（k?k）?‘symmetric‘）;?%?輸出當前窗口的灰度中值
????????for?i?=?1:M????????%?遍歷圖像
????????????for?j?=?1:N
????????????????if?fmed（ij）>fmin（ij）?&&?fmed（ij）