自适应中值滤波MATLAB

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）