积分图加速

积分图加速

常用于计算面积、加速卷积运算、避免重复计算等

1.一种简单的快速计算方法

void integral(unsigned char*inputMatrix,unsigned char*outpMatirx,int width,int height)
{
    unsigned long *columnSum=new unsigned long[width];
    for(int i=0;i<width;++i)
    {
        columnSum[i]=inputMatrix[i];
        outputMatrix[i]=inputMatrix[i];
        if(i>0)
        {
            outputMatrix[i]+=outputMatrix[i-1];
        }
    }
    for(int i=1;i<height;i++)
    {
        int offset=i*width;
        //第一列
        outputMatrix[offset]=outputMatrix[offset-width]+inputMatrix[offset];
        //其他列
        for(int j=1;j<width;++j)
        {
            outputMatrix[offset+j] = outputMatrix[offset+j-1] + outputMatrix[offset-width] + outputMatrix[offset-width-1] + inputMatrix[offset]; 
        }
    }
    return outputMatrix;
}

2.一种改进的快速积分图

void fastIntegral(unsigned char* inputMatrix, unsigned long* outputMatrix, int width, int height){
	unsigned long *columnSum = new unsigned long[width]; // sum of each column
	// calculate integral of the first line
	for(int i=0;i<width;i++){
		columnSum[i]=inputMatrix[i];
		outputMatrix[i] = inputMatrix[i];
		if(i>0){
			outputMatrix[i] += outputMatrix[i-1];
		}
	}
	for (int i=1;i<height;i++){
		int offset = i*width;
		// first column of each line
		columnSum[0] +=inputMatrix[offset];
		outputMatrix[offset] = columnSum[0];
		 // other columns 
		for(int j=1;j<width;j++){
			columnSum[j] += inputMatrix[offset+j];
			outputMatrix[offset+j] = outputMatrix[offset+j-1] + columnSum[j]; 
		}
	}
	return ;
}

参考：https://blog.csdn.net/xiaowei_cqu/article/details/17928733