jiaquanfilter.m

来自「均值滤波,加权滤波,中值滤波的MATLAB实例,用于图象及信号的处理」· M 代码 · 共 22 行

function G = jiaquanFilter(g)
J=zeros(g,g);
%中心值
x0=ceil(g/2);
y0=ceil(g/2);
J(x0,y0)=2^x0;

%计算周围值
for i=1:g
    for j=1:g
        juli=(abs(i-x0)+abs(j-y0));
        if(juli>=x0)
           J(i,j)=1;
        else
           J(i,j)=J(x0,y0)/2^juli;
        end
    end
end
zonghe=sum(J(:));
G=J./zonghe;