Создайте матрицу простой выборки.

I = magic(5)

I = 5×5

    17    24     1     8    15
    23     5     7    14    16
     4     6    13    20    22
    10    12    19    21     3
    11    18    25     2     9

В интегральном изображении каждый пиксель является совокупной суммой пикселя непосредственно выше его и с его левой стороны от него. Например, в этой тривиальной матрице, пикселе в исходной матрице в строке 1, столбец 1 (значение 17) был бы неизменен в интегральном изображении, потому что вы добавляете 0s в значение. (Функция integralImage добавляет строку 0s выше и слева от исходной матрицы.) Вычисление интегрального значения изображений для пикселя (1,2) в исходной матрице, вы добавляете пиксель выше его (0) и пиксель с его левой стороны от него (17): 24 + 17 +0 = 41. Вычисляя интегральное значение изображений для Пикселя (1,3) в originl матрице, вы добавляете пиксель выше его (0) и пиксель с его левой стороны от него, которая была уже суммирована, 41. Таким образом значение на уровне пикселя (2,4) в интегральном изображении равняется 1 + 41 + 0 = 42. Этот процесс продолжается для каждого пикселя в оригинальном изображении, как вы видите в интегральном изображении исходной матрицы, вычисленной функцией integralImage.

J = integralImage(I)

J = 6×6

     0     0     0     0     0     0
     0    17    41    42    50    65
     0    40    69    77    99   130
     0    44    79   100   142   195
     0    54   101   141   204   260
     0    65   130   195   260   325

Создайте демонстрационное изображение.

I = magic(5)

I = 5×5

    17    24     1     8    15
    23     5     7    14    16
     4     6    13    20    22
    10    12    19    21     3
    11    18    25     2     9

% Define rotated rectangular region as [x, y, width, height] where x, y
% denote the indices of the top corner of the rectangle. Width and height
% are along 45 degree lines from the top corner.
[x, y, w, h] = deal(3, 1, 3, 2);

Создайте интегральное изображение.

J = integralImage(I, 'rotated');

Вычислите сумму по области с помощью интегрального изображения.

regionSum = J(y+w+h,x+w-h+1) + J(y,x+1) - J(y+h,x-h+1) - J(y+w,x+w+1);