save results after 3d matrix neighbourhood operation

how can save a series of low resolution images after operation 10x10 pixels as a single pixel (eg, mean of 10x10 pixels)

    for j=1:10  
       for k = 50:10:850;
         for m = 50:10:300;

             results(?,?,j) = mean(mean(Image(k:k+10,m:m+10,j)));

             end
    end
end

Thanks

Answer 1

[Andrei Bobrov]

result = blockproc(Image,[10 10],@(x)mean(reshape(x.data,[],1,size(Image,3))));

[Andrei Bobrov]

corrected