自适应中值滤波器

时间:2013-11-06 11:47:32

标签: image-processing filtering matlab

我正在构建自适应中值滤波器的代码。当我执行它时,它在第12行给出了错误。没有足够的参数。在第28行。意外的MATLAB表达式。

 function f = adpmedian(g, Smax)


    %ADPMEDIAN Perform adaptive median filtering.

    %   F = ADPMEDIAN(G, SMAX) performs adaptive median filtering of

    %   image G.  The median filter starts at size 3-by-3 and iterates up   

    %   to size SMAX-by-SMAX. SMAX must be an odd integer greater than 1.

    % SMAX must be an odd, positive integer greater than 1.

      **12>>**if (Smax <= 1) || (Smax/2 == round(Smax/2)) || (Smax ~= round(Smax))
           error('SMAX must be an odd integer > 1.')
        end
        [M, N] = size(g);
        % 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');
       `28>>`    processUsingLevelB = (zmed > zmin) & (zmax > zmed) & ...
               ~alreadyProcessed; 
           zB = (g > zmin) & (zmax > g);
           outputZxy  = processUsingLevelB & zB;
           outputZmed = processUsingLevelB & ~zB;
           f(outputZxy) = g(outputZxy);
           f(outputZmed) = zmed(outputZmed);
           alreadyProcessed = alreadyProcessed | processUsingLevelB;

           if all(alreadyProcessed(:))
              break;
           end

        end

    % Output zmed for any remaining unprocessed pixels. Note that this
    % zmed was computed using a window of size Smax-by-Smax, which is
    % the final value of k in the loop.
    f(~alreadyProcessed) = zmed(~alreadyProcessed);

0 个答案:

没有答案