拖延症又犯了，上周第一篇： 机器学习系列之：怎么样数鸡 ？（1） 说到了先把图片转成灰度 （Grey Scale），接下来我们要做的就是 计算阈值 （Threshold）

当我们进行图片灰度化的时候，我们把 RGBA 图片每个相素4个字节转成了 亮度 1个字节，用了以下公式：

亮度 = R * 0.21 + G * 0.72 + B * 0.07

但是这对于后面要分类还是不太方便，我们理想只需要 0 和1 ， 1代表是鸡，0代表是空气，所以这一步计算阈值就是要计算出一个值，然后之后根据这个 阈值 (Threshold) 来二分成 0 或 1。比如 我们已经计算得到 threshold 然后可以遍历图片的每个相素，依次按 RGB 三个值的平均计算得到当前的密度(Intensity) 和这个阀值判断二分成前景和背景（鸡和空气） ：

本文的C#代码来自于 微软大神 Gary Short

大津算法 Otsu's Method

大津算法（Otsu's Method）是计算机图形学上用来把一个灰度图片退化为黑白（二值）图像。具体的算法可以看 WIKI，这里就不多说明了。

假设以下方法用于得到这个阀值。

private static int GetThreshold(Bitmap image)

参数为图片（灰度），首先我们得把图片变成字节数组（每个相素为1个字节8位）

byte[] imageBytes =  (byte[])new ImageConverter().ConvertTo(image, typeof(byte[]));

然后这个大津算法需要图片的 Histogram (直方图)，横坐标是0到255个格子，纵坐标是出现的次数。

int[] histogram = new int[256];
int ptr = 0;
while (ptr < imageBytes.Length)  {
    int h = 0xFF & imageBytes[ptr];
    histogram[h]++;
    ptr++;
}

把出现的次数计算成概率：

int totalPixels = imageBytes.Length;  // 相素点数
float sumOfIntensities = 0;
for (int t = 0; t < 256; t++) sumOfIntensities += t * histogram[t];

接下来就是大津算法得到这个阀值使得两类内方差（Variance）最小，也就是类间(Intra Class)方差最大

// 穷举阀值得到能使内间类别区别最大的值
float sumOfBackgroundIntensities = 0;
int backgroundWeight = 0;
int foregroundWeight = 0;
float maximumInterClassVariance = 0;
int threshold = 0;
for (int t = 0; t < 256; t++)
{
    backgroundWeight += histogram[t];
    if (backgroundWeight == 0) continue;

    foregroundWeight = totalPixels - backgroundWeight;
    if (foregroundWeight == 0) break;

    sumOfBackgroundIntensities += (float)(t * histogram[t]);

    float backgroundMean =
        sumOfBackgroundIntensities / backgroundWeight;

    float foregroundMean =
        (sumOfIntensities - sumOfBackgroundIntensities) / foregroundWeight;

    // 计算内间方差
    float intraClassVariance =
        (float)backgroundWeight * (float)foregroundWeight *
            (backgroundMean - foregroundMean) *
            (backgroundMean - foregroundMean);

    // 更新最大类间方差值
    if (intraClassVariance > maximumInterClassVariance)
    {
        maximumInterClassVariance = intraClassVariance;
        threshold = t;
    }
}
return threshold;

计得上周的这个鸡图么？

来，课后作业，留言告诉我，这张图的 阈值 是多少，前三名答对的每人奖励 3 SBD。

未完待续……

---

@justyy 是CN 区的点赞机器人，对优质内容进行点赞，只要代理给 @justyy 每天收利息(100 SP 每天0.04 SBD)并且能获得一次相应至少2倍的点赞，可以认为是VP 200%+ ，详细请看：

• cn区最低保障系统 上线了!
• cn区低保计划还会继续下去
• cn区低保计划（鼓励新人）真的适合你么？

中文区的大鱼 @htliao 都加入了这个计划(530 SP)表示支持，您还不来试试么？

@justyy 是 https://justyy.com 的博主 - 西半球知名的“土豪”博主。在大哥 @tumutanzi 的带领下加入了 STEEMIT。感谢您阅读 @justyy 的帖子 希望得到您的follow、upvote 或者是 reply 。

大津算法来计算阈值
机器学习系列之：怎么样数鸡 ？（1）

欢迎您发表高见，我会赞赏(Upvote)高质量的评论。

Steemit 在线工具和API接口
SteemIt Tools and APIs