积分图像

原理

代码

  close all;
  clear all;
  clc;

  I=imread('C:\Users\luoqi\Desktop\图片\影院\tree.jpg');
  if ndims(I)==3
      I=rgb2gray(I);
  end
  [M,N]=size(I);
  img=zeros(M,N);

  for i=1:M
      for j=1:N
          if(i==1&&j==1)
              img(i,j)=I(i,j);
          elseif (i==1&&j~=1)
              img(i,j)=I(i,j)+img(i,j-1);
          elseif (i~=1&&j==1)
              img(i,j)=I(i,j)+img(i-1,j);
          else
              img(i,j)=I(i,j)+img(i-1,j)+img(i,j-1)-img(i-1,j-1);
          end
      end
  end

  figure;
  subplot(121),imshow(I);
  subplot(122),imshow(img);