[图像处理]快速计算积分图.docVIP

?  HYPERLINK /xiaowei_cqu/article/details【图像处理】快速计算积分图 分类：? HYPERLINK /xiaowei_cqu/article/category/1280539 【算法分析】? HYPERLINK /xiaowei_cqu/article/category/1845851 【机器视觉】? HYPERLINK /xiaowei_cqu/article/category/923660 【图像处理】2014-01-06 21:39?4398人阅读? HYPERLINK /xiaowei_cqu/article/details\l comments 评论(14)? HYPERLINK javascript:void(0); \o 收藏 收藏? HYPERLINK /xiaowei_cqu/article/details\l report \o 举报 举报 目录 HYPERLINK /xiaowei_cqu/article/details\o 系统根据文章中H1到H6标签自动生成文章目录 (?) HYPERLINK /xiaowei_cqu/article/details\o 展开 [+] 积分图是图像中十分常用的方法，最初是在Haar特征的快速计算中学到（参考博文： HYPERLINK /xiaowei_cqu/article/details/8219324 \t _blank 利用积分图像法快速计算Haar特征），后来发现在均值滤波，二值化等图像处理方法中也十分常见。 积分图的简要介绍可以参考博文： HYPERLINK /xiaowei_cqu/article/details/8219324 \t _blank 利用积分图像法快速计算Haar特征，这里不再重复了。本篇主要是小记一下积分图的计算方法。 由于积分图中每个单元存储的信息是原图中此位置左上角所有像素之和，所以对一张W*H的图像直接求取积分图，需要： (1+2+...+W-1)*H+(1+2+...+W)*(H=1) = (w-2)*H/(W-1)+2(W+1)*(H-1)/W? 次加法。 一种简单的快速计算方法 最直接的快速计算方法是利用以计算积分求当前位置的积分，其思想正如快速计算Haar特征的方法。 即，Integral(i,j) = Integral(i,j-1) +?Integral(i-1,j) - Integral(i-1,j-1) + Image(i,j); 于是，对一张W*H的图像直接求取积分图，需要： (W-1)+(H-1)+3*(W-1)*(H-1) 次加法。代码如下： [cpp]? HYPERLINK /xiaowei_cqu/article/details\o view plain view plain HYPERLINK /xiaowei_cqu/article/details\o copy copy void?integral(unsigned?char*?inputMatrix,?unsigned?long*?outputMatrix,?int?width,?int?height){?? ????//?calculate?integral?of?the?first?line?? ????for(int?i=0;iwidth;i++){?? ????????outputMatrix[i]?=?inputMatrix[i];?? ????????if(i0){?? ????????????outputMatrix[i]?+=?outputMatrix[i-1]; ????????}?? ????}?? ????for?(int?i=1;iheight;i++){?? ????????int?offset?=?i*width;?? ????????//?first?column?of?each?line?? ????????outputMatrix[offset]?=?outputMatrix[offset-width]+inputMatrix[offset];?? ????????//?other?columns??? ????????????outputMatrix[offset+j]?=?outputMatrix[offset+j-1]?+?outputMatrix[offset-width]?+?outputMatrix[offset-width-1]?+?inputMatr