i, j = np.mgrid[0 : image.shape[0], 0 : image.shape[1]]