Calcular longitudes de línea con Python en QGIS

Encontré este artículo que describe la "forma de QGIS GUI" para calcular las longitudes de línea. http://www.qgistutorials.com/en/docs/calculating_line_lengths.html

¿Cómo puedo calcular la longitud de una línea (característica) en una capa vectorial con código Python?

Como ejemplo, quiero saber la longitud de la línea amarilla.

qgis python pyqgis

— gustavgans
fuente

Respuestas:

De los Documentos QGIS: Manejo de Geometría , puede usar el siguiente código para obtener la longitud de cualquier línea seleccionada:

layer = qgis.utils.iface.activeLayer()
features = layer.selectedFeatures()
for f in features:
    geom = f.geometry()
    print "Length:", geom.length()

¡Espero que esto ayude!

— José
fuente

Vencido por segundos ... :-)

— Alexandre Neto

@AlexandreNeto - ¡Jaja lo siento! Fue pura coincidencia porque también estaba usando este código exacto para mis capas hace solo unos momentos :)

— Joseph

@gustavgans - Muy bienvenido amigo :)

— Joseph

@Joseph, ¿cómo puedo ajustarlo para que lo calcule para todos los atributos en todas las capas de un proyecto en lugar de solo los seleccionados? He experimentado un poco, pero mi código no está haciendo el truco: layer = layer en QgsMapLayerRegistry.instance (). MapLayers (). Values () features = layer.allFeatures () para f en las características: geom = f.geometry ( ) imprimir "Longitud:", geom.length ()

— Hannes Ledegen

@HannesLedegen - ¡Estás cerca! Intenta for layer in QgsMapLayerRegistry.instance().mapLayers().values(): features = layer.getFeatures() for f in features:...

— Joseph

Gran código, sin embargo, esto solo funciona para la capa seleccionada y solo la imprime. Con algo de ayuda de otras publicaciones y Joseph lo convertí en un código que agrega un atributo a todas las capas en su proyecto con la longitud.

from PyQt4.QtCore import QVariant
for layer in QgsMapLayerRegistry.instance().mapLayers().values():
    features = layer.getFeatures()
    for f in features:
        geom = f.geometry()
        leng = geom.length()
        res = layer.dataProvider().addAttributes([QgsField("Length", QVariant.Int)])
        layer.updateFields()
        fieldIndex = layer.dataProvider().fieldNameIndex( "Length" )
        attrFeatMap = {}
        attrMap = { fieldIndex : leng }
        for feature in layer.getFeatures():
            attrFeatMap[ feature.id() ] = attrMap
        layer.dataProvider().changeAttributeValues( attrFeatMap )

— Hannes Ledegen
fuente