Le portail francophone de la géomatique

Est ce que les noeuds des  lignes correspondent strictement aux points?
Auquel cas, la couche de point ne sert à rien...

Voila une fonction qui permet de faire ce que tu demandes :

import processing
from PyQt4.QtCore import QVariant

def calculate_length(layerSource):
    '''explose les lignes sources en bi-point, ajoute un champ 'temp' et le remplis avec la longueur de la ligne'''
    #exploser les lignes en bi-point et recuperer le layer resultat
    layer = processing.runalg("qgis:explodelines", layerSource.source(), None)
    print layer
    layer = QgsVectorLayer(layer.get('OUTPUT'), layerSource.name()+'_with_length','ogr')
    print layer
    layer.startEditing()
    #creer un champ temp pour y stocker la longueur et recuperer son id
    length_field = QgsField('temp', QVariant.Double)
    layer.dataProvider().addAttributes([length_field])
    layer.commitChanges()
    layer.startEditing()
    idx_layer = layer.fieldNameIndex('temp')
    #calculer la longueur de chaque ligne et remplir le champ temp avec cette valeur
    for f in layer.getFeatures():
        xy = f.geometry()
        length  = xy.length()
        layer.changeAttributeValue(f.id(),idx_layer, length)
    layer.commitChanges()
    #ajouter la nouvelle couche au canevas
    QgsMapLayerRegistry.instance().addMapLayer(layer)
    
calculate_length(iface.activeLayer())

Sur la capture d'écran, on a l'impression que c'est la longueur entre chaque noeud qui est affichée... Or, il semblerait que ce ne soit pas ce que vous recherchez.
Avez vous la possibilité de mettre une capture d'écran plus parlante?
Est ce que c'est la longueur entre chaque point en parcourant le réseau qu'il vous faut? Auquel cas il faut faire appel à des outils de routing.

Dans cette image on distingue bien un tronçon du réseau.
Les points correspondent à des points ou sommets du réseau.
C'est la distance entre ces points qu'il me faut

l'algorithme de dijkstra que j'ai déjà utlisé me permet de trouver le chemin le plus court dans mon réseau en donnant un point de départ et d'arrivée cela fonctionne bien.
Ce chemin court généré par l'algorithme de dijkstra passe par des points et ma problématique est de trouver à présent la distance entre tous ces points à travers lequel mon linéaire passe. Ces points constituent des sommets entre autre de mon graphe.
Je vous mets ci-joint le code complet de mon traitement

from PyQt4.QtCore import *
from PyQt4.QtGui import *
from qgis.core import *
from qgis.gui import *
from qgis.networkanalysis import *
from shapely.geometry import MultiLineString
        


vl = qgis.utils.iface.mapCanvas().currentLayer()
director = QgsLineVectorLayerDirector(vl, -1, '', '', '', 3)
properter = QgsDistanceArcProperter()
director.addProperter(properter)
crs = qgis.utils.iface.mapCanvas().mapRenderer().destinationCrs()
builder = QgsGraphBuilder(crs)

pStart = QgsPoint(794035.724801,2096061.91691)
pStop = QgsPoint(794222.102796,2095908.42915)
tiedPoints = director.makeGraph(builder, [pStart, pStop])
graph = builder.graph()

tStart = tiedPoints[0]
tStop = tiedPoints[1]

idStart = graph.findVertex(tStart)
idStop = graph.findVertex(tStop)

(tree, cost) = QgsGraphAnalyzer.dijkstra(graph, idStart, 0)


if tree[idStop] == -1:
    print "pas de cheminement"
else:
    p = []
    curPos = idStop
while curPos != idStart:
    pnt = graph.vertex(graph.arc(tree[curPos]).inVertex()).point()

    p.append(pnt)
    curPos = graph.arc(tree[curPos]).outVertex()
p.append(tStart)
geom = QgsGeometry.fromPolyline(p)
print p


v_layer = None 
for layer in QgsMapLayerRegistry.instance().mapLayers().values(): 
    if layer.name() == 'cable': 
        v_layer = layer 
        pr = v_layer.dataProvider() 
        v_layer.startEditing() 
        seg = QgsFeature() 
        seg.setGeometry(geom) 
        pr.addFeatures( [seg] ) 
        v_layer.commitChanges().

Mon feature créé dans le layer cable est ma ligne.
Le but à présent est de calculer la longueur des segments formés par les sommets du graphe que cette ligne traverse.

Cordialement

Dernière modification par hotegbe (Mon 07 August 2017 12:11)

Ce que je recherche ce n'est pas la longueur totale de la ligne créée par mon code mais plutôt la longueur de chaque tronçon (segment) contenu dans ma ligne.

Dernière modification par hotegbe (Mon 07 August 2017 15:01)

Bonjour, l'application du code proposé à la ligne créee par mon code fourni le meme réultat.
Votre code revient à utiliser explode line de qgis. Il me donne la longueur de tous les segments formés par les noeuds de ma ligne.
mais moi je recherche plutôt la longueur des segments formés par les sommets de mon graphe par le quel ma ligne passe.
C'est pas la même chose

Comme vous pouvez le voir dans l'image ci-jointe.

La même ligne est affichée avec et sans nœuds.
Au départ on a une ligne simple qui traverse les sommets du graphe.
En appliquant l'outil de nœud sur cette ligne dans Qgis, on remarque que entre les sommets du graphe il existe des nœuds qui ne s'affichent que quand on applique l'outil de nœud à la ligne.
Votre code renvoie la longueur des segments contenus entre tous les nœuds y compris ceux qui ne s'affichent que quand on applique l'outil de nœud.
Le but est d'avoir uniquement la longueur compris entre les sommets du graphe sans prendre en compte les longueurs comprises entre les nœuds cachés..
J'espère que c'est plus clair maintenant.

Merci

Dernière modification par hotegbe (Tue 08 August 2017 12:56)

Comment boucler sur une paire de points? avez vous une fonction svp?
Merci