Le portail francophone de la géomatique

Bonjour,
Merci pour votre aide.

Cependant en rajoutant les deux lignes de code à mon script, je n'ai pas d'objet géographique créé. je n'ai pas non plus d'erreur.
J'ai juste rajouté les deux lignes que vous m'avez communiquées à la fin du script.

A votre avis pourquoi ça coince encore?

Je vous remercie pour votre réponse.  après quelques modifications j'obtiens un résultat qui m'a permis de bien avancer dans mon travail. Cependant,
J'ai une erreur python qui persiste "qgis dumpé" après exécution de mon code plusieurs fois.
J'ai désinstallé et changé la version de mon QGIS mais rien n'y fait; j'ai même changé d'ordinateur.

Comment régler ce souci svp?
Merci

J'ai pu avancer j'obtiens bien une couche dans laquelle une entité est créée.
Cependant le chemin le plus court recherché ne passe plus par tous les points intermédiaires entre le point de départ et d'arrivée.

Mon code est le suivant









geom = rb .asGeometry()
polyline = geom.asGeometryCollection()[0].asPolyline()


v_layer = QgsVectorLayer("LineString", "cable", "memory")
pr = v_layer.dataProvider()

seg = QgsFeature()

seg.setGeometry(QgsGeometry.fromPolyline([polyline]))

pr.addFeatures( [seg] )

v_layer.updateExtents()

QgsMapLayerRegistry.instance().addMapLayers([v_layer]).




Comment puis je imposer la géométrie de polyline à l'entité que je créee dans ma table??

Merci

Pour ça il n'y a aucun souci. c'est plutôt au niveau de la création du feature (entité) dans ma couche.
Mon souci est au niveau de la géométrie du feature que je vais insérer dans la couche que je crée. En renseignant directement tStart et tStop dans seg.setGeometry j'obtiens un tracé de ligne droite entre ces deux points et une insertion de ce tracé comme entité dans ma couche.
Du coup je n'arrive pas à imposer la géométrie de mon polyline à cette nouvelle que je crée comme ci bas.


seg = QgsFeature()
seg.setGeometry(QgsGeometry.fromPolyline([tStart ,tStop]))


Mon but est de garder la géométrie de polyline dans le seg.setGeometry.


Je ne sais pas si je suis plus clair.
Merci


voici le code plus détaillé











  for pnt in p:
    rb.addPoint(pnt)
   
polyline = QgsGeometry( rb.asGeometry() )


v_layer = QgsVectorLayer("LineString", "cable", "memory")
pr = v_layer.dataProvider()
seg = QgsFeature()
seg.setGeometry(QgsGeometry.fromPolyline([tStart ,tStop]))
pr.addFeatures( [seg] )
v_layer.updateExtents()
QgsMapLayerRegistry.instance().addMapLayers([v_layer])

Voici l'erreur que je récupère en faisant  seg.setGeometry(QgsGeometry.fromPolyline([polyline ]))



Traceback (most recent call last):
  File "<input>", line 1, in <module>
  File "c:/users/shount~1/appdata/local/temp/tmplczshv.py", line 53, in <module>
    seg.setGeometry(QgsGeometry.fromPolyline([polyline ]))
TypeError: QgsGeometry.fromPolyline(list-of-QgsPoint): argument 1 has unexpected type 'list'

je récupère toujours la meme erreur en essayant de passer polyline en paramètre comme ceci :


pr.addFeatures([polyline])
seg.setGeometry(QgsGeometry.fromPolyline([polyline]))

C'est pareil polyline n'est pas accepté en tant qu'argument ni par addfeatures() ni par setGeometry comme écrit plus haut.
même erreur

Il faut alors gérer en même temps un objet rubberband et un objet liste de points qgis

else:
  p = []
  coordinatePairs = []
  curPos = idStop
  while curPos != idStart:
    pnt = graph.vertex(graph.arc(tree[curPos]).inVertex()).point()
    print pnt
    p.append(pnt)
   coordinatePairs.append(QgsPoint(pnt))
    curPos = graph.arc(tree[curPos]).outVertex()

....
seg.setGeometry(QgsGeometry.fromPolyline([coordinatePairs]))

par contre il faudrait voir ce que donne le print pnt pour être certain qu'il renvoie bien un couple X,Y

Dernière modification par Sylvain PIERRE (Thu 29 June 2017 16:56)

Heu là il faudrait calmement repartir de l'ensemble du code...

Bonsoir,

vous pouvez passer par la construction de la polyligne
ou alors vous supprimez les ligne entre --- et vous décommentez les lignes en commentaires.

#  rb = QgsRubberBand(qgis.utils.iface.mapCanvas())
#  rb.setColor(Qt.red)
#-----
geom = QgsGeometry.fromPolyline(p)
#-----

#  for pnt in p:
#    rb.addPoint(pnt)

v_layer = QgsVectorLayer("LineString", "cable", "memory")
QgsMapLayerRegistry.instance().addMapLayers([v_layer])
pr = v_layer.dataProvider()
v_layer.startEditing()
seg = QgsFeature()
#geom = rb.asGeometry()
seg.setGeometry(geom)
pr.addFeatures( [seg] )
v_layer.commitChanges()

Cordialement,

Dernière modification par lejedi76 (Thu 29 June 2017 21:34)

Script + interaction avec l'interface = plugin

Je crois que vous allez devoir plonger un peu plus profondément dans le développement avec l'API PyQgis.
Cela ouvre des possibilités vraiment intéressantes

Sylvain

par rapport au résultat, j'obtiens à chaque fois la création d'un nouveau layer avec une nouvelle entité que je crée ce qui n'est pas mon objectif.
Il s'agit plutôt de mettre toutes les entités créées dans le même layer avec les champs qui sont déjà prédéfinis plitôt que de recréer à chaque fois un nouveau layer.

Je mets ci joint le code

 from PyQt4.QtCore import *
from PyQt4.QtGui import *
from qgis.core import *
from qgis.gui import *
from qgis.networkanalysis import *

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(4.83382,45.83818)
pStop = QgsPoint(4.841678,45.834846)
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 "Path not found"
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)

  rb = QgsRubberBand(qgis.utils.iface.mapCanvas())
  rb.setColor(Qt.red)
  
  
geom = QgsGeometry.fromPolyline(p)

for pnt in p:
    rb.addPoint(pnt)
    
    
v_layer = QgsVectorLayer("LineString", "cable", "memory")
QgsMapLayerRegistry.instance().addMapLayers([v_layer])
pr = v_layer.dataProvider()
v_layer.startEditing()
v_layer.addAttribute(QgsField("id", QVariant.String))
v_layer.addAttribute(QgsField("longueur", QVariant.String))
v_layer.addAttribute(QgsField("type", QVariant.String))
v_layer.addAttribute(QgsField("diametre", QVariant.String))
seg = QgsFeature()
geom = rb.asGeometry()
seg.setGeometry(geom)
pr.addFeatures( [seg] )
v_layer.commitChanges()