一维最大类间方差法

由Otsu（大津展之）于1978年提出的最大类间方差法，是引起较多关注的一种阈值选取方法。它是在判决分析或最小二乘原理的基础上推导出来的。参考文献:[1]OtsuN.Athresholdselectionmethodfromgray-levelhistogram.IEEETrans,1979;SMC-9;62-66下载地址算法思想:假设一幅图像有L个灰度级［1,2,...L］。灰度级为i的像素点的个数为nj,那么总的像素点个数就应该为N=ni+n2+.+nL。为了讨论方便，我们使用归一化的灰度级直方图并且视为这幅图像的概率分布：Pi=nJN.pi>0,£ =1. ⑴riI现在假设我们通过一个灰度级为k的门限将这些像素点划分为两类:C。和C/背景和目标，或者反之亦然）；C0表示灰度级为［1,...,k］的像素点，C1表示灰度级为［k+1,…丄］的像素点。那么，每一类出现的概率以及各类的平均灰度级分别由下面的式子给出：kTOC\o"1-5"\h\z=Pr（C（j）=》pf二血何 （2）i=I5=P「（C|）= £pi=1—少⑹ （3）b_iL■ ■以及(4)(5)如二工！Pr(i|Co)=£iPi/<Oo(4)(5)i“血严牛揺,“二I 1-■«i ji“血严牛揺,“二I 1-其中,讹）二EPii=ITOC\o"1-5"\h\z啲=E也 ⑺i分别为灰度级从1到k的累积出现概率和平均灰度级（一阶累积矩），而L“丁=p（L）=Eipi （8）iw1是整幅图像的平均灰度级。我们可以很容易验证，对于任意选定的k,都有:u)o+0J)=L (9)这两类的类内方差由下面的公式给出:此=E(f-Mo)2Pr(i|G)=S(<-MoFp血o (10)t i卅=£(i一“JPr =工(i-nJg/®(II)/=fc+1 i-k^I这需要二阶累积矩（second-ordercumulativemoment,统计学概念）。为了评价（灰度级k）这个门限''好〃的程度，我们需要引入判别式分析中使用的判别式标准来测量（类的分离性测量）：其中:(13)aw=COo^o+5(T；(13)=⑴o(»o—Pr)3+⑴［(“i—“f)‘=wo®】％-Po)2 (14)又根据式(9),可以得出：於=£0-^r)2Pi (⑸■*1这三个式子分别是类内方差、类间方差和灰度级的总方差。然后，我们的问题就简化为一个优化问题，即寻找一个门限k使(12)式中给出的一个目标函数取最大值。这个观点是出于这样一个猜想，一个好的阈值将会把灰度级分为两类，那么反过来说，就是如果一个门限能够在灰度级上将图像分割为最好的两类的话，那么这个门限就是最好的门限。上面给出的判别式标准是分别求取入、k和n的最大值。然而，对于k而言，它又等于另外一个，比如k=A+1；而对于入而言，又有口=入/(入+1)，因为始终存在下面的基本关系：(16)我们可以发现喘和於都是门限上的函数，但是话却与k无关二我们还泊差)，而於是基于一粧统辻(类均值)口因此，叮是判别k选取好坏的最简单於作为评价选择k作为门限的“好坏(分离性)的测量标准口屮我们使用下面的公式选择不同的k值烦序搜索，根据式(0和式(了)，或佳门限F使得I］取得最大值，或耆等价二使略达最大：心

可伙）=处伙）/於 （丽=［豐半）r呼$ （并且,最佳门限F就是"。*屮)=max必(k). (11£止＜丄效果展示:原始图像otsu分割后图像Dtsu分割后图像原始图像Otsu分割后图像算法评价：就最大类间方差算法而言，优点是算法简单，当目标与背景的面积相差不大时，能够有效地对图像进行分割。但是，当图像中的目标与背景的面积相差很大时，表现为直方图没有明显的双峰，或者两个峰的大小相差很大，分割效果不佳，或者目标与背景的灰度有较大的重叠时也不能准确的将目标与背景分