MATLAB数字图像处理 图像分割技术 迭代法

本文为博主原创，转载请注明出处:

http://blog.sina.com.cn/s/blog_156532ba70102wmjf.html

迭代法首先选择一个近似阈值T，将图像分割成两部分R1和R2，计算区域R1和R2的均值μ1和μ2，选择新的分割阈值T = （μ1 + μ2）/2 ，重复上述步骤知道μ1和μ2（或T）不再变化为止。这里我已经将迭代法写成了函数，直接调用就可以了。函数源代码如下：

function [ TK ] = diedai( I )  % 图像分割之迭代法获取阈值  % 形参I为原始图像灰度图，返回值为阈值    I = imread('rice.jpg');  ZMAX = max(I);%取出最大灰度值  ZMIN = min(I);%取出最小灰度值  TK = (ZMAX+ZMIN)/2;%初始阈值  BCal = 1;  iSize = size(I);%图像大小  while (BCal)      %定义前景和背景数      iForeground = 0;      iBackground = 0;      %定义前景和背景灰度总和      ForegroundSUM = 0;      BackgroundSUM = 0;      for i = 1:iSize(1)          for j = 1:iSize(2)              tmp = I(i,j);              if tmp >= TK                  %前景灰度值                  iForeground = iForeground + 1;                  ForegroundSUM = ForegroundSUM + double(tmp);              else                  iBackground = iBackground + 1;                  BackgroundSUM = BackgroundSUM + double(tmp);              end          end      end      %计算前景和背景的平均值      ZO = ForegroundSUM/iForeground;      ZB = BackgroundSUM/iBackground;      TKtemp = uint8((ZO + ZB)/2);%新的阈值      if TKtemp == TK          BCal = 0;      else         TK = TKtemp;      end      %当阈值不再变化时，迭代结束      end  disp(strcat('迭代后的阈值：',num2str(TK)));  end  灰度图。

下面简单说下怎么调用，只有一点需要注意，迭代法操作对象为灰度矩阵，不是灰度图的要先转为灰度图。代码如下：

I1 = imread('rice.jpg');  I = rgb2gray(I1);  TK = diedai(I);  newI = imbinarize(I,double(TK/255));  

我这里的‘rice.jpg’图片不是灰度图，所以就要先转为灰度图。
效果图我就不放了。图像分割的方法有很多，所以并不是所有的图都适合用迭代法来分割，有的图甚至效果还会非常差，所以想要获得好的分割效果只能不断尝试好的方法。所以，加油吧！我也是刚学，写的不好请见谅！