六、图像的重建与复原（MATLAB）下

一、仅有噪声的复原——空间滤波

如果出现的退化仅是噪声，所选择的方法是空间滤波。
1，空间噪声滤波器
根据各种类型滤波器公式，定义M函数，实现各种滤波器

function f=gmean(g,m,n)
%几何平均滤波器
[g,revertclass]=tofloat(g);
f=exp(imfilter(log(g),ones(m,n),'replicate')).^(1/m/n);
f=revertclass(f);

function f=harmean(g,m,n)
%调和均值滤波器
[g,revertclass]=tofloat(g);
f=m*n./imfilter(1./(g+eps),ones(m,n),'replicate');
f=revertclass(f);

function f=charmean(g,m,n,q)
%逆调和均值滤波器
[g,revertclass]=tofloat(g);
f=imfilter(g.^(q+1),ones(m,n),'replicate');
f=f./(imfilter(g.^q,ones(m,n),'replicate')+eps);
f=revertclass(f);

function f=alphatrim(g,m,n,d)
%修正的alpha均值滤波器
if (d<=0)||(d/2~=round(d/2))%round 四舍五入函数
    error('d must be a positive,even integer')
end
[g,revertclass]=tofloat(g);
f=imfilter(g,ones(m,n),'symmetric');
for k=1:d/2
    f=f-ordfilt2(g,k,ones(m,n),'symmetric');
end
for k=(m*n-(d/2)+1):m*n
    f=f-ordfilt2(g,k,ones(m,n),'symmetric');
end
f=f/(m*n-d);
f=revertclass(f);

定义spfilt函数执行执行上述空间滤波

function f=spfilt(g,type,varargin)
[m,n,Q,d]=processInputs(varargin{:});
switch type
    case 'amean'
        w=fspecial('