我有以下坐标
minx, maxx, miny ,maxy = 448262.080078, 450360.750122, 6262492.020081, 6262938.950073
我希望使用python创建一个大小为1 m的正方形网格。
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
这是附加到任何特定的GIS平台还是以纯python的方式进行而没有任何指定的输出格式(例如shapefile,textfile等)的要求
感谢@Dan,我希望以纯Python格式执行,输出将为shapefile格式
—
Gianni
抱歉,我没有使用商业软件。我更喜欢使用Java,Python,C ++等纯语言编写的程序。
—
詹尼
但是您不介意使用诸如GDAL / OGR(pypi.python.org/pypi/GDAL)或pyshp(pypi.python.org/pypi/pyshp)之类的库吗?
—
Snorfalorpagus