... и бесплатный геопространственный питон в этом; /
DavidF
4
Один из возможных способов продолжения использования SAGA GIS http://www.saga-gis.org
После открытия шейп-файла запустите следующие 3 модуля: 1. Модули \ Shapes \ Tools \ Get Shapes экстент
Модули \ Фигуры \ Инструменты \ Точки \ Точки из линий [вопреки тому, что следует из названия, ты также можешь использовать это для получения точек из многоугольника]
Модули \ Фигуры \ Инструменты \ Точки \ Добавить координаты к точкам Это даст вам таблицу, которая содержит координаты x и y 4 углов ограничительной рамки вашего файла многоугольника.
В arcgis вот код Python. результат - список minx, miny, maxx, maxy, minM, maxM, minZ, maxZ (
import arcpy
for feat in arcpy.SearchCursor(r"c:\data\f.gdb\counties"):print feat.Shape.extent
-2.6685272725154649.4265363633626-2.5284818181812149.5079454546192NaNNaNNaNNaN-10.46333636378251.4455454544593-6.0130545458304555.3799909091533NaNNaNNaNNaN-4.7777818182761454.0555454544593-4.3534727268846854.4100000000002NaNNaNNaNNaN
Вот версия R, использующая пример данных из пакета rgdal:
library(rgdal)
dsn <- system.file("vectors/ps_cant_31.MIF",package="rgdal")[1]
d <- readOGR(dsn = dsn, layer="ps_cant_31")## transform if this is not longlatif(is.projected(d)) d <- spTransform(d, CRS("+proj=longlat +ellps=WGS84"))for(i in1:nrow(d)){print(bbox(d[i,]))}
import fiona
from shapely.geometry import shape
with fiona.open(r'd:\Projects\_00_Data\_USstates\fe_2007_us_state00.shp','r')as features:for i, feat in enumerate(features):
geom = shape(feat['geometry'])
name = feat['properties']['NAME00']print','.join((name,)+ tuple([str(i)for i in geom.bounds]))
Ответы:
Используя модуль Python ogr из OSGEO, этот пример даст вам кортеж, содержащий координаты, которые определяют конверт для каждой функции.
источник
Один из возможных способов продолжения использования SAGA GIS http://www.saga-gis.org После открытия шейп-файла запустите следующие 3 модуля: 1. Модули \ Shapes \ Tools \ Get Shapes экстент
Модули \ Фигуры \ Инструменты \ Точки \ Точки из линий [вопреки тому, что следует из названия, ты также можешь использовать это для получения точек из многоугольника]
Модули \ Фигуры \ Инструменты \ Точки \ Добавить координаты к точкам Это даст вам таблицу, которая содержит координаты x и y 4 углов ограничительной рамки вашего файла многоугольника.
источник
В arcgis вот код Python. результат - список minx, miny, maxx, maxy, minM, maxM, minZ, maxZ (
источник
Вот версия R, использующая пример данных из пакета rgdal:
источник
Я использую Фиона и стройные для таких задач:
источник