OpenCV阈值分割_二值化
Datawhale 计算机视觉基础-图像处理（上）-Task05 图像分割/二值化
5.1 简介
该部分的学习内容是对经典的阈值分割算法进行回顾，图像阈值化分割是一种传统的最常用的图像分割方法，因其实现简单、计算量小、性能较稳定而成为图像分割中最基本和应用最广泛的分割技术。它特别适用于目标和背景占据不同灰度级范围的图像。它不仅可以极大的压缩数据量，而且也大大简化了分析和处理步骤，因此在很多情况下，是进行图像分析、特征提取与模式识别之前的必要的图像预处理过程。图像阈值化的目的是要按照灰度级，对像素集合进行一个划分，得到的每个子集形成一个与现实景物相对应的区域，各个区域内部具有一致的属性，而相邻区域不具有这种一致属性。这样的划分可以通过从灰度级出发选取一个或多个阈值来实现。
5.2 学习目标
-
了解阈值分割基本概念
-
理解最大类间方差法（大津法）、自适应阈值分割的原理
-
掌握OpenCV框架下上述阈值分割算法API的使用
5.3 内容介绍
1、最大类间方差法、自适应阈值分割的原理
2、OpenCV代码实践
3、动手实践并打卡（读者完成）
5.4 算法理论介绍
5.4.1 最大类间方差法（大津阈值法）
大津法（OTSU）是一种确定图像二值化分割阈值的算法，由日本学者大津于1979年提出。从大津法的原理上来讲，该方法又称作最大类间方差法，因为按照大津法求得的阈值进行图像二值化分割后，前景与背景图像的类间方差最大。
它被认为是图像分割中阈值选取的最佳算法，计算简单，不受图像亮度和对比度的影响，因此在数字图像处理上得到了广泛的应用。它是按图像的灰度特性，将图像分成背景和前景两部分。因方差是灰度分布均匀性的一种度量,背景和前景之间的类间方差越大,说明构成图像的两部分的差别越大,当部分前景错分为背景或部分背景错分为前景都会导致两部分差别变小。因此,使类间方差最大的分割意味着错分概率最小。
其实可以这么理解，比如一张图的像素直方图如下：
OTSU就是找到一个阈值k，使左右两个分布的类间方差最小
应用： 是求图像全局阈值的最佳方法，应用不言而喻，适用于大部分需要求图像全局阈值的场合。
优点： 计算简单快速，不受图像亮度和对比度的影响。
缺点： 对图像噪声敏感；只能针对单一目标分割；当目标和背景大小比例悬殊、类间方差函数可能呈现双峰或者多峰，这个时候效果不好。
原理非常简单，涉及的知识点就是均值、方差等概念和一些公式推导。为了便于理解，我们从目的入手，反推一下这著名的OTSU算法。
求类间方差：
OTSU算法的假设是存在阈值TH将图像所有像素分为两类C1(小于TH)和C2(大于TH)，则这两类像素各自的均值就为m1、m2，图像全局均值为mG。同时像素被分为C1和C2类的概率分别为p1、p2，pi为i这个像素值的个数除以总像素值的个数。因此就有：
根据原文，式（4）还可以进一步变形：
分割：
这个分割就是二值化，OpenCV给了以下几种方式，很简单，可以参考：
!
5.4.2 三角法
三角法求阈值最早见于Zack的论文《Automatic measurement of sister chromatid exchange frequency》主要是用于染色体的研究，该方法是使用直方图数据，基于纯几何方法来寻找最佳阈值，它的成立条件是假设直方图最大波峰在靠近最亮的一侧，然后通过三角形求得最大直线距离，根据最大直线距离对应的直方图灰度等级即为分割阈值，图示如下：
对上图的详细解释：
在直方图上从最高峰处bmx到最暗对应直方图bmin(p=0)%构造一条直线，从bmin处开始计算每个对应的直方图b到直线的垂直距离，知道bmax为止，其中最大距离对应的直方图位置即为图像二值化对应的阈值T。
扩展情况：
有时候最大波峰对应位置不在直方图最亮一侧，而在暗的一侧，这样就需要翻转直方图，翻转之后求得值，用255减去即得到为阈值T。扩展情况的直方图表示如下：
三角阈值化处理细胞之类的图片会有比较好的效果
5.4.3 自适应阈值
前面介绍了OTSU算法，但这算法属于全局阈值法，所以对于某些光照不均的图像，这种全局阈值分割的方法会显得苍白无力，如下图：
显然，这样的阈值处理结果不是我们想要的，那么就需要一种方法来应对这样的情况。
这种办法就是自适应阈值法(adaptiveThreshold)，它的思想不是计算全局图像的阈值，而是根据图像不同区域亮度分布，计算其局部阈值，所以对于图像不同区域，能够自适应计算不同的阈值，因此被称为自适应阈值法。(其实就是局部阈值法)
如何确定局部阈值呢？可以计算某个邻域(局部)的均值、中值、高斯加权平均(高斯滤波)来确定阈值。值得说明的是：如果用局部的均值作为局部的阈值，就是常说的移动平均法。
例如采用方法 CV_ADAPTIVE_THRESH_MEAN_C，阈值类型：CV_THRESH_BINARY, 阈值的象素邻域大小 block_size 选取3，参数param1 取3时，opencv输出结果：
下面我们来手算一下：
( 218 + 219 + . . . + 180 + 180 ) 9 − 3 = 193 \frac{(218+219+...+180+180)}{9}-3=193 9(218+219+...+180+180)−3=193
因为193大于190，所以取0
( 191 + 190 + . . . + 183 + 144 ) 9 − 3 = 175 \frac{(191+190+...+183+144)}{9}-3=175 9(191+190+...+183+144)−3=175
因为175小于180，所以取255
高斯的处理方法也和均值相同，只不过核换成高斯核
5.5 基于OpenCV的实现
- 工具：OpenCV3.1.0+VS2013
- 平台：WIN10
函数原型（c++）
1.最大类间方差法
double cv::threshold ( InputArray src,
OutputArray dst,
double thresh,
double maxval,
int type
)
参数：
- src — input array (single-channel, 8-bit or 32-bit floating point).
- dst — output array of the same size and type as src.
- thresh — threshold value.
- maxval — maximum value to use with the THRESH_BINARY and THRESH_BINARY_INV thresholding types.
- type — thresholding type 参考：thresholdType
1.自适应阈值
void adaptiveThreshold(InputArray src, OutputArray dst,
double maxValue,
int adaptiveMethod,
int thresholdType,
int blockSize, double C)
参数：
Parameters
- src — Source 8-bit single-channel image.
- dst — Destination image of the same size and the same type as src.
- maxValue — Non-zero value assigned to the pixels for which the condition is satisfied
- adaptiveMethod — Adaptive thresholding algorithm to use,参考：cv::AdaptiveThresholdTypes
- thresholdType — Thresholding type that must be either THRESH_BINARY or THRESH_BINARY_INV, 可参考：thresholdType
blockSize Size of a pixel neighborhood that is used to calculate a threshold value for the pixel: 3, 5, 7, and so on. - C — Constant subtracted from the mean or weighted mean (see the details below). Normally, it is positive but may be zero or negative as well.
实现示例（c++)
- 1、大津阈值
#include - 2、自适应阈值
#include 进阶实现(根据原理自己实现)
实现示例（c++)
- 1、大津阈值
#include - 2、自适应阈值
#include 效果
- 1、大津阈值
- 2、自适应阈值
相关技术文档、博客、教材、项目推荐
opencv文档: https://docs.opencv.org/3.1.0/d7/d1b/group__imgproc__misc.html#gae8a4a146d1ca78c626a53577199e9c57
博客：https://blog.csdn.net/weixin_40647819/article/details/90179953
https://blog.csdn.net/weixin_40647819/article/details/90213858
python版本：https://www.kancloud.cn/aollo/aolloopencv/267591 http://www.woshicver.com/FifthSection/4_3_%E5%9B%BE%E5%83%8F%E9%98%88%E5%80%BC/
5.6 总结
该部分对两种经典阈值分割方法进行了介绍，读者可根据提供的资料进行学习，然后参考示例代码自行实现。Otsu的二值化有一些优化方法，读者可以尝试学习并实现。
Task05 阈值分割/二值化END.
— By: 小武
博客：https://blog.csdn.net/weixin_40647819
关于Datawhale：
Datawhale是一个专注于数据科学与AI领域的开源组织，汇集了众多领域院校和知名企业的优秀学习者，聚合了一群有开源精神和探索精神的团队成员。Datawhale以“for the learner，和学习者一起成长”为愿景，鼓励真实地展现自我、开放包容、互信互助、敢于试错和勇于担当。同时Datawhale 用开源的理念去探索开源内容、开源学习和开源方案，赋能人才培养，助力人才成长，建立起人与人，人与知识，人与企业和人与未来的联结。