/**************************************************************************************************/
// 函数名称:OnMenu020503()
// 函数功能:“对角积分图”菜单。
// 函数参数:
// 输入参数: 无
// 输出参数: 无
// 返 回 值:void
// 创建作者:(QQ:370711753)
// 修改日期:2017/11/08 16:25:58
/**************************************************************************************************/
void COpenCVDlg::OnMenu020503()
{
try // 错误处理
{
Mat lv_MatImageIntegralSum = Mat();
Mat lv_MatImageIntegralSqSum = Mat();
Mat lv_MatImageIntegralTilted = Mat();
Mat lv_MatImageIntegralNorm = Mat();
// 计算积分图像
cv::integral(m_MatImageRead, lv_MatImageIntegralSum,
lv_MatImageIntegralSqSum, lv_MatImageIntegralTilted, CV_64F, CV_64F);
// 图像数据归一化
cv::normalize(lv_MatImageIntegralTilted, lv_MatImageIntegralTilted, 0, 255, CV_MINMAX);
convertScaleAbs(lv_MatImageIntegralTilted,lv_MatImageIntegralNorm); // 精度转换为8位INT整型
// 显示结果
HV_ImageShow(lv_MatImageIntegralNorm, m_strNameWindow);
}
catch (Exception& lv_Exception) // 抛出异常
{
CString lv_strException = _T("");
lv_strException.Format(_T(":%s!"), lv_Exception.what());
MessageBox(lv_strException, _T("错误!"), MB_ICONERROR);
return;
}
return;
}
//************************************************
积分图算法由Crow在1984年首次提出,是为了在多尺度透视投影中提高渲染速度。积分图算法是一种快速计算图像区域和以及图像区域平方和的算法。它的核心思想就是对每一个图像建立起自己的积分图查找表,在图像处理的阶段就可以根据预先建立积分图查找表直接查找从而实现对均值卷积的线性时间计算。做到了卷积执行的时间与窗口大小无关。这种算法被应用到基于NCC的快速匹配、对象检测和SURF变换、边缘检测、基于统计学的快速滤波器等方面。
积分图算法的缺点:1、需要较大的内存来存储积分图。2、如果图片较大,会导致数据溢出。
积分图由原图像计算而来,假设原图大小为W*H,则积分图大小为(W+1)*(H+1)。在积分图(integral image)上任意坐标(x,y)处的值,表示原图中坐标为(x,y)的点的左上角所有像素点像素值的和(平方和表中的就是平方和)。和表和平方和表建立公式如下:
和表:
平方和表:
如下图所示,假设输入图像大小为2x2,则积分图大小为3x3
如上图所示,如果想求输入图像中蓝色区域内的像素值之和(3+2+5+4=14),只要根据和表积分图进行两次减法和一次减法即可:46+10-22-20=14。也就是右下角+左上角-右上角-左下角。
opencv中计算积分图的API介绍:
//API一
void integral( InputArray src, //输入图像
OutputArray sum, //和表
int sdepth = -1,); //和表深度
//API二
void integral( InputArray src, //输入图像
OutputArray sum, //和表
OutputArray sqsum, //平方和表
int sdepth = -1, //和表深度,一般为CV_32S
int sqdepth = -1 ); //平方和表深度,一般为CV_32F
计算积分图并显示:
//
//opencv4.1.0
//
#include
using namespace std;
using namespace cv;
int main() {
Mat src, sum, sqrsum;
src = imread("1.png", 0);
integral(src, sum, sqrsum, CV_32S, CV_32F);
normalize(sum, sum, 0, 255, NORM_MINMAX, CV_8UC1, Mat());
normalize(sqrsum, sqrsum, 0, 255, NORM_MINMAX, CV_8UC1, Mat());
imshow("原图", src);
imshow("和表积分图", sum);
imshow("平方和表积分图", sqrsum);
waitKey(0);
return 0;
}
//
//opencv4.1.0
//
#include
using namespace std;
using namespace cv;
int get_block_sum(Mat &sum, int x1, int y1, int x2, int y2);
int main() {
Mat src, sum, sqrsum;
src = imread("1.png", 0);
//获取ROI,便于观察(使用ImageWatch插件)
int LR = 4;
int LC = 4;
int width = 6;
int height = 6;
Mat roi = src(Rect(LR, LC, width, height)).clone();
//计算积分图
integral(src, sum, sqrsum, CV_32S, CV_32F);
//利用积分图计算ROI内的像素值总和
int value = get_block_sum(sum, LR, LC, LR + height, LC + width);
cout << value << endl;
imshow("原图", src);
waitKey(0);
return 0;
}
int get_block_sum(Mat &sum, int x1, int y1, int x2, int y2) {
int BottomRight = sum.at
int TopLeft = sum.at
int TopRight = sum.at
int BottomLeft = sum.at
int sum_value = (BottomRight + TopLeft - TopRight - BottomLeft);
return sum_value;
}
//********************************************
Paul Viola和Michael Jones在2001年首次将积分图应用在图像特征提取上,在他们的论文“Rapid Object Detection using a Boosted Cascade of Simple Features”中,积分图被当作一种新的图像特征表征方式,可以把检测的Haar特征非常高效的计算出来,用于实时人脸检测系统。
积分图是一种能够描述全局信息的矩阵表示方法,其构造方式是积分图像上位置(i,j)处的值ii(i,j)是原图像(i,j)左上角方向所有像素的和。
利用积分图可以可以快速的计算图像上某一区域内的像素和,如下图:
要计算区域D内的像素和,只需要获取到积分图上1、2、3/4点各自的像素值,分别表述为ii(1)、ii(2)、ii(3)、ii(4),则区域D的像素和=ii(4)+ii(1)-ii(2)-ii(3)。
传统的计算像素和的方式需要遍历区域D内所有的像素,再执行累加,计算量随着区域D面积的增大而增大,而对积分图方式来说,只需要在计算积分图之后,通过简单几次加减运算就可以得到某一区域内“像素和”这一特征,计算速度非常快,并且这种速度的提升效果随着区域面积的增大和计算次数的增多表现的更为明显。
Opencv中使用integral函数计算积分图。
void integral( InputArray src, OutputArray sum, int sdepth=-1 );
第一个参数src,可以使灰度图或RGB彩色图,单通道和三通道均可作为输入,但每个通道的精度必须是8位int或32位、64位浮点型;
第二个参数sum,积分图,若输入src是灰度图,则积分图也是灰度图,若输入src是RGB三通道图,则积分图sum也是RGB三通道彩色。sum的图像深度是32位整型或32位、64位浮点型,这取决于第三个参数sdepth的定义;
第三个参数sdepth,定义积分图的深度(depth),32位整型或者32位、64位浮点型。注意图像的深度跟图像的通道数是无关的,相关概念可以参看这里:Opencv Mat矩阵中data、size、depth、elemSize、step等属性的理解 ;
所以在使用sum之前,声明即可,可以不事先定义Mat矩阵的大小和数据类型。
#include "highgui/highgui.hpp"
#include "imgproc/imgproc.hpp"
using namespace cv;
int main(int argc,char *argv[])
{
Mat image=imread(argv[1]);
//cvtColor(image,image,CV_RGB2GRAY); //原图像是三通道,积分图也是三通道
Mat imageIntegral;
integral(image,imageIntegral,CV_32F); //计算积分图
normalize(imageIntegral,imageIntegral,0,255,CV_MINMAX); //归一化,方便显示
Mat imageIntegralNorm;
convertScaleAbs(imageIntegral,imageIntegralNorm); //精度转换为8位int整型
imshow("Source Image",image);
imshow("Integral Image",imageIntegralNorm);
waitKey();
}
原图:
积分图:
显示的积分图是经过归一化后的。可以看到从左上角到右下角,图像是越来越亮的,也就是说积分图上像素值是越来越大的。
积分图的图像通道数跟原始图像保存一致,当计算积分图的原始图像是彩色图像时,积分图也是彩色图像,这时候积分图计算的不是灰度,而是颜色:
//************************************