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

12

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

mossplix
источник
Пожалуйста, уточните: у вас есть шейп- файл или файл Excel ?
whuber
1
Название действительно не отражает ваш вопрос, пожалуйста, подумайте о его редактировании.
DavidF
у меня есть шейп-
файл

Ответы:

25

Используя модуль 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
DavidF
источник
2
... и бесплатный геопространственный питон в этом; /
DavidF
4

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

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

  2. Модули \ Фигуры \ Инструменты \ Точки \ Добавить координаты к точкам Это даст вам таблицу, которая содержит координаты x и y 4 углов ограничительной рамки вашего файла многоугольника.

johanvdw
источник
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
gotchula
источник
4

Вот версия 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,]))    
}
mdsumner
источник
1

Я использую Фиона и стройные для таких задач:

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]))
Матей
источник
Это НЕ обеспечивает Lat / Lon.
Harvpan
Выходные данные выглядят следующим образом: координаты широты / долготы: -124,72583900000001,45,544321, -116,915989,49,002494 -82,626182,37.202467, -77,71951899999999,40,638801 -111.056888,40,996345999999996, -105,9996, -5, 92998671, 67, 5991, 87, 5991, 67, 599, 87, 87, ПБ
Matej