Как я могу преобразовать шейп-файл в широту и долготу?

у меня есть шейп-файл страновых подразделений, и я хотел бы извлечь массив ограничивающих лат и долг для каждого подразделения ... это можно сделать?

Используя модуль Python ogr из OSGEO, этот пример даст вам кортеж, содержащий координаты, которые определяют конверт для каждой функции.

from osgeo import ogr

ds = ogr.Open("mn_counties.shp")
lyr = ds.GetLayerByName("mn_counties")

lyr.ResetReading()

for feat in lyr:
    # get bounding coords in minx, maxx, miny, maxy format
    env = feat.GetGeometryRef().GetEnvelope()
    # get bounding coords in minx, miny, maxx, maxy format
    bbox = [env[0], env[2], env[1], env[3]]
    print env
    print bbox
    print

Один из возможных способов продолжения использования SAGA GIS http://www.saga-gis.org После открытия шейп-файла запустите следующие 3 модуля: 1. Модули \ Shapes \ Tools \ Get Shapes экстент

1. Модули \ Фигуры \ Инструменты \ Точки \ Точки из линий [вопреки тому, что следует из названия, ты также можешь использовать это для получения точек из многоугольника]

2. Модули \ Фигуры \ Инструменты \ Точки \ Добавить координаты к точкам Это даст вам таблицу, которая содержит координаты 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.66852727251546 49.4265363633626 -2.52848181818121 49.5079454546192 NaN NaN NaN NaN
-10.463336363782 51.4455454544593 -6.01305454583045 55.3799909091533 NaN NaN NaN NaN
-4.77778181827614 54.0555454544593 -4.35347272688468 54.4100000000002 NaN NaN NaN NaN

Вот версия 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 longlat
if (is.projected(d)) d <- spTransform(d, CRS("+proj=longlat +ellps=WGS84"))

for (i in 1: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]))