Pobieranie współrzędnych linii lub wielu linii w QGIS przez PyQGIS

Dec 03 2020

Możemy uzyskać współrzędne warstwy punktowej za pomocą następującego kodu, ale chcę znać współrzędne linii lub multilinii. Jak zdobyć te współrzędne?

layer = QgsProject.instance().mapLayersByName("layerName")[0]
iter = layer.getFeatures()
for feature in iter:
  geom=feature.geometry()
  print(geom.asPoint().x(),geom.asPoint().y())

Odpowiedzi

4 Zoltan Dec 03 2020 at 14:01

W przypadku połączeń wieloczęściowych możesz użyć następującego kodu (przetestowałem go w 3.16.1)

layer = QgsProject.instance().mapLayersByName("folyo")[0]
for feature in layer.getFeatures():
    for part in feature.geometry().asMultiPolyline():
        print("part")
        for pnt in part:
            print(pnt.x(), pnt.y())

W przypadku pojedynczego elementu linestring:

layer = QgsProject.instance().mapLayersByName("Single parts")[0]
for feature in layer.getFeatures():
    for pnt in feature.geometry().asPolyline():
        print(pnt.x(), pnt.y())

asMultiPolyline i asPolyline zwraca listę Punktów, więc możesz używać indeksów.

UPDATE aby sprawdzić typ geometrii warstwy użyj metody geometryType warstwy. Aby oddzielić typy pojedyncze od wielu, użyj funkcji isSingleType. Oto kod:

layer = QgsProject.instance().mapLayersByName("your_layer")[0]
if layer.geometryType() ==  QgsWkbTypes.LineGeometry:
    for feature in layer.getFeatures():
        geom = feature.geometry()
        if QgsWkbTypes.isSingleType(geom.wkbType()):
            # single
            for pnt in geom.asPolyline():
                print(pnt.x(), pnt.y())
        else:
            # multipart
            for part in geom.asMultiPolyline():
                print("part")
                for pnt in part:
                    print(pnt.x(), pnt.y())
2 ahmadhanb Dec 03 2020 at 14:03

Jeśli chcesz wyodrębnić xy wierzchołków geometrii MultiPolyline, możesz użyć następującego kodu:

layer = QgsProject.instance().mapLayersByName("layerName")[0]
iter = layer.getFeatures()
for feature in iter:
    geom=feature.geometry()
    verts = geom.asMultiPolyline() # use geom.asPolyline() for line geometry
    for vert in verts:
        for i in range(len(vert)):
            print(vert[i].x(),vert[i].y())

Wynik będzie taki

-3487202.783144876 10010395.30478035
-3428476.505393841 10049595.900699263
-3428476.505393841 10049595.900699263
-3374593.995628879 9923727.261801671
-3428476.505393841 10049595.900699263
-3359626.0983730364 10059037.052994486