У меня есть следующие координаты
minx, maxx, miny ,maxy = 448262.080078, 450360.750122, 6262492.020081, 6262938.950073
Я хочу создать квадратную сетку размером 1 м, используя 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
Прикреплено ли это к какой-либо конкретной ГИС-платформе или является обязательным требованием сделать это в чистом Python без какого-либо указанного выходного формата (например, шейп-файл, текстовый файл и т. Д.)
Спасибо @Dan, я хочу выступить на чистом питоне, и результат будет в формате shapefile
—
Gianni
Уровень лицензии ArcInfo ArcMap имеет инструмент Fishnet, но вы не указали, как вы собираетесь создавать шейп-файл.
Извините, я не пользуюсь коммерческим ПО. Я предпочитаю программы на чистом языке Java, Python, C ++.
—
Джанни
Но вы не против использовать такие библиотеки, как GDAL / OGR ( pypi.python.org/pypi/GDAL ) или pyshp ( pypi.python.org/pypi/pyshp )?
—
Снорфалорпаг