myfilter.m

来自「滤波器与伪彩色变换」· M 代码 · 共 18 行

function result = myfilter (f , H , D0);

%根据输入图象f求输出图象result
F = fft2(f);
F1 = movetocenter(F);
[x,y] = size(F1);
%G = F.*H;

G =zeros(x,y);
for i = 1:x
    for j = 1:y
        G(i,j) = F1(i,j) * H(i,j);
    end
end

result = ifft2(G);

result = uint8( abs(result) );