pmc.m

来自「超声图像实时消噪实验程序」· M 代码 · 共 22 行

    function g = pmc( f, k, stepsize, nosteps, verbose )
    % Perona and Malik Diffusion
        if verbose
            figure(verbose);
            subplot(1,2,1); imshow(f); title('Original Image'); drawnow;
        end
        g = f;
        for i=1:nosteps
                    gx = gD( g, 1, 1, 0 );
                    gy = gD( g, 1, 0, 1 );
                    grad2 = gx.*gx + gy.*gy;
                    c = exp( -grad2 / (k^2) );

            g = g + stepsize * snldStep( g, c );

                if verbose
                    figure(verbose);
                    subplot(1,2,2); imshow(g);
                    title('Perona and Malik Diffusion'); drawnow;
                end
        end