autosmooth.m

用Matlab编写的自动多阈值分割程序

function  [imm,n]=autosmooth(im)
% im=imread('rice.tif');
% im=imnoise(a,'gaussian',0,0.005);
% a=imread('mytu.tif');
% im=rgb2gray(a);
H=imhist(im);
for i=1:255
    u(i)=H(i+1)-H(i); %求差分曲线
end
k=1;
for i=1:254         %找出差分曲线由负到正的过零点的位置（谷底位置）
    if u(i)<0
        if u(i+1)<0      
            i=i+1;
        elseif u(i+1)>0
               d(k)=i+0.5;                
               k=k+1;
        elseif u(i+1)==0
               for j=1:50
                   if (i+j<255&u(i+j)==0&u(i+j+1)>0)                      
                       d(k)=i+0.5*j+0.5;                
                       k=k+1;
                       break;
                   else 
                       break;
                   end
               end
           end
       end
  end
  for i=1:k-2
      D(i)=d(i+1)-d(i);     %求各小峰的峰宽
  end
  dd=uint8(D);
  Hh=imhist(dd);            %统计峰宽度的直方图
  DD=find(Hh);            
  kd=size(DD);
  n=round(0.7*kd(1));        %求出平滑窗口的尺寸
  imm=wiener2(im,[n,n]);    %平滑原始直方图