How can I efficiently process a numpy array in blocks similar to Matlab's blkproc (blockproc) function

Question

I\'m looking for a good approach for efficiently dividing an image into small regions, processing each region separately, and then re-assembling the results from each proces

Answer 1

Process by slices/views. Concatenation is very expensive.

for x in xrange(0, 160, 16):
    for y in xrange(0, 160, 16):
        view = A[x:x+16, y:y+16]
        view[:,:] = fun(view)