Obtendo as coordenadas de linha ou MultiLine no QGIS por meio do PyQGIS

Question

Podemos obter coordenadas de camada de ponto por meio do código a seguir, mas eu quero saber as coordenadas de linha ou multilinha. Como obter essas coordenadas?

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

Para cadeias de linha multipartes, você pode usar o seguinte código (eu testei em 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())

Para cadeia de linha de parte única:

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

asMultiPolyline e asPolyline retornam uma lista de pontos, para que você possa usar índices.

ATUALIZE para verificar o tipo de geometria da camada, use o método geometryType da camada. Para separar tipos únicos e múltiplos, use a função isSingleType. Aqui está o código para isso:

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

Se você deseja extrair o xy dos vértices de uma geometria MultiPolyline, você pode usar o seguinte código:

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

A saída será assim

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