convsepbrd.m

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

    function g = convSepBrd( f, w1, w2 )
    % convolve along colums + along rows with repetition of the border
    N = size(f,1);
    M = size(f,2);
    K = (size(w1(:),1)-1)/2;
    L = (size(w2(:),1)-1)/2;
    iind = min(max((1:(N+2*K))-K,1),N);
    jind = min(max((1:(M+2*L))-L,1),M);
    fwb = f(iind,jind);
    g=conv2(w1,w2,fwb,'valid');