积分图的每一点(x, y)的值是原图中对应位置的左上角区域的所有值得和:
而且,积分图可以只遍历一次图像即可有效的计算出来,因为积分图每一点的(x, y)值是:
一旦积分图计算完毕,对任意矩形区域的和的计算就可以在常数时间内完成。如下图中,阴影矩形区域的值:
%构建积分图的代码
clear all;
close all;
clc;
img=double(imread('lena.jpg'));
[m n]=size(img);
%计算积分图
I=zeros(m,n);
for i=1:m
for j=1:n
if i==1 && j==1 %积分图像左上角
I(i,j)=img(i,j);
elseif i==1 && j~=1 %积分图像第一行
I(i,j)=I(i,j-1)+img(i,j);
elseif i~=1 && j==1 %积分图像第一列
I(i,j)=I(i-1,j)+img(i,j);
else %积分图像其它像素
I(i,j)=img(i,j)+I(i-1,j)+I(i,j-1)-I(i-1,j-1);
end
end
end
int_img = I;
%int_img就是输出的积分图
%利用积分图求一张图片中某个矩形块的和
stap = 4;%stap是矩形块的大小4x4。
[M,N] = size(img);%image是整个图片
for i=1:stap:M
for j=1:stap:N
patch = img(i:i+stap-1,j:j+stap-1);从img中取一块
if i+stap-1==stap && j+stap-1==stap
sum =int_img(i+stap-1,j+stap-1);
elseif i+stap-1==stap && j+stap-1~=stap %第一行
sum =(int_img(i+stap-1,j+stap-1) - int_img(i+stap-1,j-1));
elseif i+stap-1~=stap && j+stap-1==stap %第一列
sum =(int_img(i+stap-1,j+stap-1) - int_img(i-1,j+stap-1));
else
sum =(int_img(i+stap-1,j+stap-1) + int_img(i-1,j-1) - int_img(i-1,j+stap-1) - int_img(i+stap-1,j-1)); %其他
end
end
end