% 积分图算法
% 样本灰度图(double)
[m n] = size(Im);
s = zeros(m+1,n+1); % 初始化迭代矩阵
ii = zeros(m+1,n+1);
%根据公式迭代计算积分图
for i = 2:m+1
for j = 2:n+1
s(i,j) = s(i,j-1)+Im(i-1,j-1);
ii(i,j) = ii(i-1,j)+s(i,j);
end
end