图像积分图的计算

/*
*@function IntegralImage.cpp
*@brief 计算积分图
*@author ltc
*@date 18:14 Tuesday，December 1st，2015
*/

#include<opencv2\core\core.hpp>
#include<opencv2\highgui\highgui.hpp>
#include<opencv2\imgproc\imgproc.hpp>
#include<iostream>

using namespace std;
using namespace cv;

int main(int argc,char* argv[])
{
//读取图像
Mat image;
image=imread("lena.bmp",IMREAD_COLOR);
if(image.empty())
{
cerr<<"Failure in loading image"<<endl;
return -1;
}
//转换成灰度图像
Mat grayImage;
cvtColor(image,grayImage,COLOR_BGR2GRAY);
//定义、计算积分图像，积分图像比灰度图像多一行一列
Mat integralImage=Mat::zeros(grayImage.rows+1,grayImage.cols+1,CV_32SC1);
for(int i=0;i<grayImage.rows;i++)
{
for(int j=0;j<grayImage.cols;j++)
{
integralImage.at<int>(i+1,j+1)=integralImage.at<int>(i,j+1)
+integralImage.at<int>(i+1,j)
-integralImage.at<int>(i,j)
+grayImage.at<unsigned char>(i,j);

}
}
//真正的积分图像
Mat realIntegralImage=integralImage(Range(1,integralImage.rows),Range(1,integralImage.cols));
/*cout<<realIntegralImage.row(0)<<endl;*/
//Opencv中自带的积分图像计算
Mat integralImage_integral;
integral(grayImage,integralImage_integral);
//下面两行代码用来对比自己的积分图像和Opencv自己计算的积分图像是否一样
//cout<<integralImage.row(1)<<endl;
//cout<<integralImage_integral.row(1)<<endl;
imshow("src",image);
imshow("gray",grayImage);
waitKey(0);
return 0;
}