It may be likely for me to make errors for the coordinate issues. The coordinate for center of (i,j)-th pixel grid should be (i+1/2, j+1/2)