I have the following coordinates minx, maxx, miny ,maxy = 448262.080078, 450360.750122, 6262492.020081, 6262938.950073 I wish to create a square grid of size 1 m using python. import math minx,maxx,miny,maxy = 448262.080078, 450360.750122, 6262492.020081, 6262938.950073 size = 1 def set_bbox(minx, maxx, miny, maxy, distx, disty): nx = int(math.ceil(abs(maxx - minx)/distx)) ny = int(math.ceil(abs(maxy - miny)/disty)) new_maxx = minx + (nx*distx) new_miny = maxy - (ny*disty) return ((minx, new_maxx, new_miny, maxy),ny,nx) # shift the bottom (right - down) coord, ny, nx = set_bbox(minx,maxx,miny,maxy,size,size) # left-up origin origin = coord[0],coord[3] # number of tiles ncell = ny*nx