Это поможет вам начать. Это изображение с пониженной дискретизацией, где каждое значение находится между нулем и квадратом размера блока. Вы самостоятельно расширяете это число в подматрицу.
bs = 8
a = imread('trees.tif');
[r,c] = size(a);
d = imresize(a,[round(r/bs), round(c/bs)]);
figure(1)
imshow(a)
figure(2)
imshow(d)
mv = max(d(:))
d = round(double(d)/double(mv)*bs*bs);
figure(3)
imagesc(d)