OGR ve Python kullanarak çokgenin tüm köşelerini alın

Python OGR API ile biraz sorun yaşıyorum. Yapmaya çalıştığım, bir çokgenin dış halkasının her bir tepe noktasının tüm koordinatlarını elde etmektir.

Şimdiye kadar sahip olduğum şey bu:

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()

Ben sadece forbölgenin üzerinde olduğunu duydum ama bu sadece çokgen halkaları döndürür, düğümleri değil.

Herkes yardım edebilir.

Biraz dosya formatınıza ve geometrisinize bağlıdır, ancak prensipte devam böyle görünebilir.

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

FYI - orijinal soruya dayalı tam bir kod örneği ve hemen kullanabileceğiniz bir komut dosyası için .. bkz. Https://github.com/spatialguru/NME/blob/master/ogr_explode.py

Ben sadece aynı problemle karşılaştım. Ogr ExportToJson işlevini kullanarak ve sonra bir sözlük içine Json dize okuma sona erdi. Verilerimi ve orijinal sorudaki gösterimi kullanarak, şuna benzer:

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'}

Sadece şekil dosyalarına bakıyorsanız, pyshp'yi de kullanabilirsiniz .

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