OGR 및 Python을 사용하여 다각형의 모든 정점을 가져옵니다.

Python OGR API에 약간의 문제가 있습니다. 내가하려고하는 것은 다각형 바깥 고리의 각 꼭지점의 모든 좌표를 얻는 것입니다.

이것이 내가 지금까지 가진 것입니다.

import osgeo.ogr
import glob

path = "/home/woo/maps/"
out = path + 'output.txt'

file = open(out,'w')
for filename in glob.glob(path + "*.shp"):
    ds = osgeo.ogr.Open(filename)
    layer1 = ds.GetLayer(0)
    print layer1.GetExtent()    
    for feat in layer1:
        geom = feat.GetGeometryRef()
        ring = geom.GetGeometryRef(0)
        points = ring.GetPointCount()
        #Not sure what to do here


file.close()

나는 당신 for이 그 지역을 뛰어 넘을 수는 있지만 노드가 아닌 다각형의 고리 만 반환한다고 들었습니다 .

도울 수있는 사람

파일 형식과 지오메트리에 따라 조금씩 다르지만 원칙적으로 연속은 다음과 같습니다.

  for p in xrange(points):
        lon, lat, z = ring.GetPoint(p)

참고-원래 질문을 기반으로 한 완전한 코드 예제 및 즉시 사용할 수있는 스크립트 .. https://github.com/spatialguru/NME/blob/master/ogr_explode.py 참조

방금 같은 문제가 발생했습니다. ogr에서 ExportToJson 함수를 사용한 다음 Json 문자열을 사전으로 읽었습니다. 내 데이터와 원래 질문의 표기법을 사용하면 다음과 같습니다.

import json
...
ring_dict = json.loads(ring.ExportToJson())
ring_dict

{'coordinates': [[-4.94237, 55.725449],
  [-4.941922, 55.725585],
  [-4.9420024, 55.7252119],
  [-4.9422001, 55.7250997],
  [-4.9423197, 55.7251789],
  [-4.9425472, 55.7253089],
  [-4.94237, 55.725449]],
 'type': 'LineString'}

shapefile 만보고 있다면 pyshp 를 사용할 수도 있습니다 .

import shapefile
sf = shapefile.Reader("shapefiles/blockgroups")
shapes = sf.shapes()
for shape in shapes:
  for vertex in shape.points:
    #do something with the vertex