Maybe
blockproc
would help?
fun = @(block_struct) ...
mode(block_struct.data(:)) * ones(size(block_struct.data));
output_matrix = blockproc(input_matrix,[<size_x> <size_y>],fun);
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!