matlab中图像的阈值分割,基于MATLAB的图像阈值分割技术汇总
数字图像处理课程论文
基于MATLAB的图像阈值分割技术
摘要：本文主要针对图像阈值分割做一个基于MATLAB的分析。通过双峰法，迭
代法以及OUTS法三种算法来实现图像阈值分割，并且就这三种算法做了一定的分析和比较，在加椒盐的图片上同时进行三种实验，做出比较，最终得出实践结论。
关键词：图像分割 MATLAB 阈值分割 算法 引言：图像分割是图像处理与计算机视觉领域低层次视觉中最为基础和重要的领域之一，它是对图像进行视觉分析和模式识别的基本前提．同时它也是一个经典难题，到目前为止既不存在一种通用的图像分割方法，也不存在一种判断是否分割成功的客观标准，图像阈值分割即是其中的一种方法。 阈值分割技术因其实现简单、计算量小、性能较稳定而成为图像分割中最基本和应用最广泛的分割技术，已被应用于很多的领域，在很多图像处理系统中都是必不可少的一个环节。
1、阈值分割思想和原理
若图像中目标和背景具有不同的灰度集合：目标灰度集合与背景灰度集合，且两个灰度集合可用一个灰度级阈值T进行分割。这样就可以用阈值分割灰度级的方法在图像中分割出目标区域与背景区域，这种方法称为灰度阈值分割方法。
在物体与背景有较强的对比度的图像中，此种方法应用特别有效。比如说物体内部灰度分布均匀一致，背景在另一个灰度级上也分布均匀，这时利用阈值可以将目标与背景分割得很好。如果目标和背景的差别是某些其他特征而不是灰度特征时，那么先将这些特征差别转化为灰度差别，然后再应用阈值分割方法进行处理，这样使用阈值分割技术也可能是有效的
设图像为f(x,y)，其灰度集范围是[0,L]，在0和L之间选择一个合适的灰度阈值T，则图像分割方法可由下式描述：
f
g(x,y){1(x,y)?T0f(x,y)?T
这样得到的g(x,y)是一幅二值图像。
(一)原理研究
图像阈值分割的方法有很多，在这里就其中三种方法进行研究，双峰法，迭代法，以及OUTS法。
方法一：双峰法 方法原理：
双峰法是一种简单的阈值分割方法。双峰法先将原图转为灰度图，然后将灰度图转为灰度直方图，灰度直方图就是灰度级的像素数ni与灰度i的二维关系，它反映了一幅图像上灰度分布的统计特性，在MATLAB中使用函数imhist
数字图像处理课程论文
来实现。如果得到的灰度直方图呈现明显的双峰状，则选取双峰之间的谷底所对应的灰度级作为阈值分割。
方法二：迭代法 方法原理：
开始时候选择一个阈值作为初始估计值，然后按着某种策略不断得改进这个估计值，直到满足给定的准则为止。
(1) 求出图像最大灰度值Max和最小灰度值Min，初始阈值估计值
T0.=(Max+Min)*0.5.
(2) 用T0将图像分割为目标和前景。图像分成两组像素：Z1由所有灰度值大
于或等于T0的像素组成，而Z0由所有灰度值小于T0的像素组成。 (3) 分别求出两者区域中的所有像素计算平均灰度值a1和a2。 (4) 计算新的阈值T1=(a1+a2)*0.5.
如果|Ti+1-Ti|<0.5，则退出循环，Ti+1即为所求阈值；否则，将Ti+1复制给Ti，重复(2)~(5)。
方法三：OTSU法 方法原理 OTUS又称：最大类间法，该算法是在使类间方差最大的自动确定阈值的方法，
是在判决分析最小二乘法原理的基础上推到得出的，其算法比较简单，是一种方便可行的阈值选取方法。
设原始灰度图像灰度级范围为[0,L],灰度级为i的像素点数为ni，则图像的全部像素数为：
N?no?n1?......?nL?1
?Pi?0L?1i?1
把图像中的像素按灰度值用阈值t分成两类C0和C1，由灰度值在[0,t]之间的像素组成，由灰度值在[t+1,L-1]之间的像素组成，对于灰度分布概率，整幅图的均值为
ut??iPi
i?0L?1因此，C0和C1的均值为：
μ0??ipi?otiwo?μ?t?w?t?μ1?其中
ori?t?1?ipt0L-1iμT?μ?t??w11?w?t?i
w?P(C)??p?w?t?i?0
w?P(C)??p?1?w(t)1r1i?t?1iL?1上面三式可得
数字图像处理课程论文
ut= w0u0+ w1u1 类间方差定义为：
σ
算法：
2B?w0(μ0?μT)?w1(μ0?μT)?w0w1(μ1?μo)222
让t在[0,L-1]范围一次取值，使类间方差最大的他值即为OUTS法的最佳阈值。
MATLAB工具箱提供的graythresh函数求取阈值。
1、双峰法图像阈值分割matlab code：
I=imread('cat.jpeg'); I=rgb2gray(I); imhist(I) 直方图：
说明：根据双峰法原理，观察到灰度图像直方图呈现明显的双峰状，则选取双峰之间的谷底所对应的灰度级作为阈值分割。如上图，选取230作为分割点。
I=imread('cat.jpeg'); I=rgb2gray(I); figure
subplot(1,2,1) imshow(I);
[width,height]=size(I); title('原图')
数字图像处理课程论文
for i=1:width for j=1:height if(I(i,j)<230) RC(i,j)=0; else
RC(i,j)=1; end end end
subplot(1,2,2) imshow(RC)
title('双峰法图像阈值分割处理效果图')
此图为C=230
此图为C=150