Получение линейных или многострочных координат в QGIS через PyQGIS

Question

Мы можем получить координаты точечного слоя с помощью следующего кода, но я хочу знать координаты линии или многострочности. Как получить эти координаты?

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

Answer 1

Для составной строки вы можете использовать следующий код (я тестировал его в 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())

Для однослойной линейки:

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

asMultiPolyline и asPolyline возвращают список точек, поэтому вы можете использовать индексы.

ОБНОВЛЕНИЕ для проверки типа геометрии слоя используйте метод слоя geometryType. Для разделения одного и нескольких типов используйте функцию isSingleType. Вот код для этого:

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

Answer 2

Если вы хотите извлечь xy вершин геометрии MultiPolyline, вы можете использовать следующий код:

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

Результат будет таким

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