tnldstep.m

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

    function r = tnldStep( L, a, b, c )
    % Discrete numerical scheme of dL/dt for tensor diffusion
    N = size( L, 1 );
    M = size( L, 2 );
        Lpc = translateImage( L, 1, 0 );
        Lpp = translateImage( L, 1, 1 );
        Lcp = translateImage( L, 0, 1 );
        Lmp = translateImage( L, -1, 1 );
        Lmc = translateImage( L, -1, 0 );
        Lmm = translateImage( L, -1, -1 );
        Lcm = translateImage( L, 0, -1 );
        Lpm = translateImage( L, 1, -1 );
        amc = translateImage( a, -1, 0 );
        apc = translateImage( a, +1, 0 );
        bmc = translateImage( b, -1, 0 );
        bcm = translateImage( b, 0, -1 );
        bpc = translateImage( b, +1, 0 );
        bcp = translateImage( b, 0, +1 );
        ccp = translateImage( c, 0, +1 );
        ccm = translateImage( c, 0, -1 );
        r = -1/4 * (bmc+bcp) .* Lmp + ...
             1/2 * (ccp+c)   .* Lcp + ...
             1/4 * (bpc+bcp) .* Lpp + ...
             1/2 * (amc+a)   .* Lmc - ...
             1/2 * (amc+2*a+apc+ccm+2*c+ccp) .* L + ...
             1/2 * (apc+a)   .* Lpc + ...
             1/4 * (bmc+bcm) .* Lmm + ...
             1/2 * (ccm+c)   .* Lcm - ...
             1/4 * (bpc+bcm) .* Lpm;