Not sure how you want to define block number but how about:
blockNumber = 1;
for row = 1 : blockSizeR : rows
for col = 1 : blockSizeC : columns
row1 = row;
row2 = row1 + blockSizeR - 1;
row2 = min(rows, row2);
col1 = col;
col2 = col1 + blockSizeC - 1;
col2 = min(columns, col2);
oneBlock = grayImage(row1:row2, col1:col2);
Mean(blockNumber) = mean2(oneBlock);
fprintf('The mean for block #%d = $f', ...
blockNumber, Mean(blockNumber));
blockNumber = blockNumber + 1;
end
end
