Otsu最大类间方差法图像二值化MATLAB实现

OTSU算法是由日本学者OTSU于1979年提出的一种对图像进行二值化的高效算法，可根据算法提取平坦地区的建筑物。
参考文献地址
该方法将图像分为前景和背景两部分，背景和目标之间的类间方差越大，说明构成图像的两部分的差别越大，当部分目标错分为背景或部分背景错分为目标都会导致两部分差别变小。因此，使类间方差最大的分割意味着错分概率最小。
最大类间方差法（otsu）的公式推导：
记t为前景与背景的分割阈值，前景点数占图像比例为w0，平均灰度为u0；背景点数占图像比例为w1，平均灰度为u1。
则图像的总平均灰度为：u=w0*u0+w1*u1。
前景和背景图象的方差：g=w0*(u0-u)(u0-u)+w1(u1-u)(u1-u)=w0*w1(u0-u1)*(u0-u1)，此公式为方差公式。
当方差g最大时，可以认为此时前景和背景差异最大，此时的灰度t是最佳阈值
sb = w0*w1*(u1-u0)*(u0-u1)
MATLAB实现：

function [b,Out_image, th] = f_Otsu(a)
   [m,n] =