i, j = np.mgrid[-half_size : half_size + 1, -half_size : half_size + 1].astype(