import numpy as np

def rebin(a, n):
   s = a.shape
   return np.tile(np.transpose(np.tile(a, (1,n)).reshape(s[0],n,s[1]), axes=(0,2,1)), (n,1)).reshape(s[0],n,n*s[1]).reshape(n*s[0],n*s[1])