数字图像处理matlab函数spfilt代码、spfilt函数代码、tofloat函数代码

spfilt函数执行线性和非线性的空间滤波器

function f=spfilt(g,type,varargin)
%spfilt执行线性和非线性的空间滤波器，g为原图像，type为滤波器类型，
[m,n,Q,d]=processInputs(varargin{:});
%开始执行滤波
switch type
case 'amean' %算数平均滤波
    w=fspecial('average',[m n]);
    f=imfilter(g,w,'replicate');
case 'gmean' %几何平均滤波
    f=gmean(g,m,n);
case 'hmean'%调和平均滤波
    f=harmean(g,m,n);
case 'chmean'%反调和平均滤波，Q的默认值是1.5
    f=chmean(g,m,n,Q);
case 'median'%中值滤波
    f=medfilt2(g,[m n],'symmetric');
case 'max'%最大值滤波
    f=ordfilt2(g,m*n,ones(m,n),'symmetric');
case 'min'%最小值滤波
    f=ordfilt2(g,1,ones(m,n),'symmetric');
case 'midpoint'%中值滤波
    f1=ordfilt2(g,1,ones(m,n),'symmetric');
    f2=ordfilt2(g,mn,ones(m,n),'symmetric');
    f=imlincomb(0.5,f1,0.5,f2);
case 'atrimmed'%顺序平均值滤波
    f=alphatrim(g,m,n,d);
otherwise
    error('Unkown filter type.')
end
end
% ------------------------------------------------- %
function f=gmean(g,m,n)
[g,revertClass]=tofloat(g); %tofloat在另一个.m函数
f=exp(imfilter(log(g),ones(m,n),'replicate')).^(1/m/n);
f=revertClass(f);
end
% ------------------------------------------------- %
function f=harmean(g,m,n)
[g,revertClass]=tofloat(g);
f=m*n./imfilter(1./(g+eps),ones(m,n),'replicate');
f=revertClass(f);
end
% ------------------------------------------------- %
function f=chmean(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);
end
% ------------------------------------------------- %
function f=alphatrim(g,m,n,d)
if(d<=0)||(d/2~=round(d/2))
    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=(mn-d/2+1):mn
    f=f-ordfilt2(g,k,ones(m,n),'symmetric');
end
f=f/(m*n-d);
f=revertClass(f);
end
% ------------------------------------------------- %
function [m,n,Q,d]=processInputs(varargin)
m=3;
n=3;
Q=1.5;
d=2;
if nargin>0
    m=varargin{1};
end
if nargin>1
    n=varargin{2};
end
if nargin>2
    Q=varargin{3};
    d=varargin{3};
end
end
% ------------------------------------------------- %
function [out,revertclass]=tofloat(in)
%tofloat convert image to floating point
identity = @(x)x;
tosingle = @im2single;
table={'uint8',tosingle,@im2uint8
    'uint16',tosingle,@im2uint16
    'int16',tosingle,@im2int16
    'logical',tosingle,@logical
    'double',identity,identity
    'single',identity,identity};
classIndex=find(strcmp(class(in),table(:,1)));
if isempty(classIndex)
error('Unsupported inut image class.')
end

out=table{classIndex,2}(in);
revertclass = table{classIndex,3};
end