若该文为原创文章,未经允许不得转载
原博主博客地址:https://blog.csdn.net/qq21497936
原博主博客导航:https://blog.csdn.net/qq21497936/article/details/102478062
本文章博客地址:https://blog.csdn.net/qq21497936/article/details/104833663
各位读者,知识无穷而人力有穷,要么改需求,要么找专业人士,要么自己研究
OpenCV开发专栏
《OpenCV开发笔记(〇):使用mingw530_32编译openCV3.4.1源码,搭建Qt5.9.3的openCV开发环境》
《OpenCV开发笔记(三):OpenCV图像的概念和基本操作》
《OpenCV开发笔记(四):OpenCV图片和视频数据的读取与存储》
《OpenCV开发笔记(五):OpenCV读取与操作摄像头》
《OpenCV开发笔记(六):OpenCV基础数据结构、颜色转换函数和颜色空间》
《OpenCV开发笔记(八):OpenCV常用操作之计时、缩放、旋转、镜像》
《OpenCV开发笔记(九):OpenCV区域图像(ROI)和整体、局部图像混合》
《OpenCV开发笔记(十):OpenCV图像颜色通道分离和图像颜色多通道混合》
《OpenCV开发笔记(十一):OpenCV编译支持Gpu(cuda) 加速开发之win-qt-mingw32编译》
《OpenCV开发笔记(十二):OpenCV编译支持Gpu(cuda) 加速开发之win-qt-msvc2015编译(opencv3.4.0、cuda9.0、VS2015)》
《OpenCV开发笔记(十三):OpenCV图像对比度、亮度的调整》
《OpenCV开发笔记(十四):算法基础之线性滤波-方框滤波》
《OpenCV开发笔记(十五):算法基础之线性滤波-均值滤波》
《OpenCV开发笔记(十六):算法基础之线性滤波-高斯滤波》
《OpenCV开发笔记(十七):算法基础之线性滤波对比-方框、均值、高斯滤波》
《OpenCV开发笔记(十八):算法基础之非线性滤波-中值滤波》
《OpenCV开发笔记(十九):算法基础之非线性滤波-双边滤波》
《OpenCV开发笔记(二十):算法基础之非线性滤波对比-中值、双边滤波》
《OpenCV开发笔记(二十一):算法基础之形态学滤波-膨胀》
《OpenCV开发笔记(二十二):算法基础之形态学滤波-腐蚀》
《OpenCV开发笔记(二十三):算法基础之形态学滤波-开运算》
《OpenCV开发笔记(二十四):算法基础之形态学滤波-闭运算》
《OpenCV开发笔记(二十五):算法基础之形态学滤波-形态学梯度》
《OpenCV开发笔记(二十六):算法基础之形态学滤波-顶帽(礼帽)》
《OpenCV开发笔记(二十七):算法基础之形态学滤波-黑帽》
《OpenCV开发笔记(二十八):带你学习图像识别之阈值化》
《OpenCV开发笔记(二十九):带你学习图像识别之自适应阈值》
《OpenCV开发笔记(三十):带你学习图像识别之经典OTSU算法阈值化》
《OpenCV开发笔记(三十一):带你学习图像识别之双阈值化》
《OpenCV开发笔记(三十二):带你学习图像识别之半阈值化》
持续补充中…
OpenCV开发笔记(三十):带你学习图像识别之经典OTSU算法阈值化
前言
本篇继续讲解阈值化,本篇章有点特殊是需要手码算法的OTSU阈值化。
Demo
(注意:左上为原图,右上为灰度图,左下角为threshold做otsu算法阈值化,中下为代码的otsu算法,右下为threshold做三角形算法)
OTSU阈值化
概述
OTSU是经典常用的阈值化算法,是由日本学者OTSU于1979年提出的一种对图像进行二值化的高效算法,是一种自适应的阈值确定的方法,又称大津阈值分割法,是最小二乘法意义下的最优分割是一种寻找图像阈值的最大类间方差的算法。
OTSU算法的具体步骤如下:
- 步骤一:计算灰度级中每个像素在整幅图像中的个数;
- 步骤二:计算每个像素在整幅图像的概率分布;
- 步骤三:对灰度级进行遍历搜索,计算当前灰度值下前景背景类间概率;
- 步骤四:通过目标函数计算出类内与类间方差下对应的阈值。
阈值化函数原型
double threshold( InputArray src,
OutputArray dst,
double thresh,
double maxval,
int type );
- 参数一:InputArray类型,一般是cv::Mat,且可以处理多通道,8或者32位浮点(注意:当使用THRESH_BINARY处理多通道的时候,每个通道都会进行阈值化,比如RGB三通道,那么可能R比G,B大,当阈值设置为大于G、B小于R时,则R为最大是,显示红色,其他类型的阈值形式类推)。
- 参数二;OutputArray类型,输出的目标图像,需要和原图片有一样的尺寸和类型。
- 参数三:double类型的thresh,阈值。
- 参数四:double类型的maxval,与“THRESH_BINARY”枚举和“THRESH_BINARY_INV”枚举一起使用才有效果,其他枚举忽略。
- 参数五:int类型的type,阈值类型。
Demo源码
void OpenCVManager::testOtsu()
{
QString fileName1 = "E:/qtProject/openCVDemo/openCVDemo/modules/openCVManager/images/1.jpg";
cv::Mat srcMat = cv::imread(fileName1.toStdString());
int width = 300;
int height = 200;
cv::resize(srcMat, srcMat, cv::Size(width, height));
cv::String windowName = _windowTitle.toStdString();
cvui::init(windowName);
if(!srcMat.data)
{
qDebug() << __FILE__ << __LINE__
<< "Failed to load image:" << fileName1;
return;
}
// 灰度化图像
cv::Mat grayMat;
cv::cvtColor(srcMat, grayMat, cv::COLOR_BGR2GRAY);
// 传入otsm算法获取计算的阈值
int thresh = otsuAlgorithm(grayMat);
qDebug() << __FILE__ << __LINE__ <<thresh;
cv::Mat windowMat = cv::Mat(cv::Size(srcMat.cols * 3, srcMat.rows * 2),
srcMat.type());
cv::Mat dstMat;
while(true)
{
windowMat = cv::Scalar(0, 0, 0);
cvui::printf(windowMat, width * 1 + 60, height * 0 + 100, "otsuThresh = %d", thresh);
cv::Mat grayMat2;
cv::cvtColor(grayMat, grayMat2, CV_GRAY2BGR);
// 原始图像
dstMat = windowMat(cv::Range(srcMat.rows * 0, srcMat.rows * 1),
cv::Range(srcMat.cols * 0, srcMat.cols * 1));
cv::addWeighted(dstMat, 0.0, srcMat, 1.0, 0.0, dstMat);
// 转化的灰度图像
dstMat = windowMat(cv::Range(srcMat.rows * 0, srcMat.rows * 1),
cv::Range(srcMat.cols * 2, srcMat.cols * 3));
cv::addWeighted(dstMat, 0.0, grayMat2, 1.0, 0.0, dstMat);
// 使用阈值函数计算otsu
dstMat = windowMat(cv::Range(srcMat.rows * 1, srcMat.rows * 2),
cv::Range(srcMat.cols * 0, srcMat.cols * 1));
cv::Mat grayMat3;
cv::threshold(grayMat, grayMat3, thresh, 255, CV_THRESH_OTSU);
cv::cvtColor(grayMat3, grayMat2, CV_GRAY2BGR);
cv::addWeighted(dstMat, 0.0, grayMat2, 1.0, 0.0, dstMat);
// 使用自己写的计算的otsu
dstMat = windowMat(cv::Range(srcMat.rows * 1, srcMat.rows * 2),
cv::Range(srcMat.cols * 1, srcMat.cols * 2));
cv::Mat grayMat4;
cv::threshold(grayMat, grayMat4, thresh, 255, CV_THRESH_BINARY);
cv::cvtColor(grayMat4, grayMat2, CV_GRAY2BGR);
cv::addWeighted(dstMat, 0.0, grayMat2, 1.0, 0.0, dstMat);
// 额外:三角形使用自己写的计算的otsu
dstMat = windowMat(cv::Range(srcMat.rows * 1, srcMat.rows * 2),
cv::Range(srcMat.cols * 2, srcMat.cols * 3));
cv::Mat grayMat5;
cv::threshold(grayMat, grayMat5, thresh, 255, CV_THRESH_TRIANGLE);
cv::cvtColor(grayMat5, grayMat2, CV_GRAY2BGR);
cv::addWeighted(dstMat, 0.0, grayMat2, 1.0, 0.0, dstMat);
// 更新
cvui::update();
// 显示
cv::imshow(windowName, windowMat);
// esc键退出
if(cv::waitKey(25) == 27)
{
break;
}
}
}
int OpenCVManager::otsuAlgorithm(cv::Mat srcMat)
{
int threshold = -1;
// 检查图像是否是8位数据,单通道
if(srcMat.type() != CV_8UC1)
{
qDebug() << __FILE__ << __LINE__ << "srcMat must be CV_8UC1";
return -1;
}
int grayScale[256] = {0}; // 每个灰度级所占像素的
double grayPro[256] = {0}; // 每个灰度级所占像素比例: 为 该像素出现的次数/总像素
// 步骤一:计算灰度级中每个像素在整幅图像中的个数;灰度级8位为256级别0~255
for(int row = 0; row < srcMat.rows; row++)
{
for(int col = 0; col < srcMat.cols; col++)
{
grayScale[srcMat.at<uchar>(row, col)]++;
#if 0
// 整数行设置成白色
if(row % 10 == 0)
{
srcMat.at<uchar>(row, col) = 255;
}
// 整数列设置成黑色
if(col % 10 == 0)
{
srcMat.at<uchar>(row, col) = 0;
}
#endif
}
}
// 步骤二:计算每个像素在整幅图像的概率分布;
int pixSum = srcMat.rows * srcMat.cols;
for(int index = 0; index < 256; index++)
{
grayPro[index] = (float)grayScale[index] / pixSum;
}
// 步骤三:便利所有灰度级的阈值分割条件,测试哪一个的类间方差最大
// 注意:方差公式 g = w0 * w1 * (u1 - u2)^2
double deltaMax = 0;
double w0 = 0.0f;
double w1 = 0.0f;
double u0 = 0.0f;
double u1 = 0.0f;
double u0Temp = 0.0f;
double u1Temp = 0.0f;
double deltaTemp = 0.0f;
for(int i = 0; i < 256; i++)
{
w0 = 0.0f;
w1 = 0.0f;
u0 = 0.0f;
u1 = 0.0f;
u0Temp = 0.0f;
u1Temp = 0.0f;
deltaTemp = 0.0f;
// 步骤三:对灰度级进行遍历搜索,计算当前灰度值下前景背景类间概率;
for(int j = 0; j < 256; j++)
{
// 背景
if(j <= i)
{
w0 += grayPro[j];
u0Temp += j * grayPro[j];
}else // 前景
{
w1 += grayPro[j];
u1Temp += j * grayPro[j];
}
}
u0 = u0Temp / w0;
u1 = u1Temp / w1;
// 方差公式 g = w0 * w1 * (u1 - u2)^2
deltaTemp = (float)(w0 * w1 * pow((u0 - u1), 2));
qDebug() << __FILE__ <<__LINE__
<< w0 << w1 << u0 << u1;
// 步骤四:通过目标函数计算出类内与类间方差下对应的阈值。
if(deltaTemp > deltaMax)
{
deltaMax = deltaTemp;
threshold = i;
}
}
return threshold;
}
工程模板:对应版本号v1.25.0
对应版本号v1.25.0
原博主博客地址:https://blog.csdn.net/qq21497936
原博主博客导航:https://blog.csdn.net/qq21497936/article/details/102478062
本文章博客地址:https://blog.csdn.net/qq21497936/article/details/104833663