图像处理遇到求积分图像,自己写了个小代码,看似很简单,但调试过程中出现了两个问题。最后测试代码如下:
#include <stdio.h>
#include <iostream>
#include <vector>
using namespace std;
int main()
{
int width=9;
int height=4;
vector<double> rowsum;
vector<double> outimg;
vector<double> test;
for (int i=0; i<height+1; i++)
{
//vector<double> rowsum;
if (i!=0)
{
rowsum.clear();
}
for (int j=0; j<width+1 ;j++)
{
if (i==0)
{
outimg.push_back(0);
rowsum.push_back(0);
}
else if (j==0)
{
outimg.push_back(0);
rowsum.push_back(0);
}
else
{
rowsum.push_back(rowsum[j-1]+1);
outimg.push_back(outimg[(i-1)*(width+1)+j]+rowsum[j]);
}
test.push_back(outimg[i*(width+1)+j]);
cout<<i<<" "<<j<<": "<<test[i*(width+1)+j]<<endl;
}
}
getchar();
return 0;
}
逻辑很简单,就是一个迭代过程。rowsum代表当行该像素前像素和,outimg中存储积分图像。第一次错误发生在未加上rowsun.clear()这部分,第二个错误发生在for循环中终止条件改为i<height+1, j<width+1,但是下面outimg中行数还是height,导致最后输出结果不对,折腾了半个多小时才想起来,新矩阵行列都加1,代码中未改正,改正后正确输出。
代码rowsum.push_back(rowsum[j-1]+1); 这一行最后的1表示原来图像当前位置的像素值,我都假设为1,实际可以是srcimg[(i-1)*width+j-1]。