[QGIS 1.x] Identification vertex

[david12312]

bonjour jai une petite problème je veux identifier les points d'intersection pour réaliser la segmentation au point vertex seulement jai definie une méthode pour la réaliser mais mon probleme reste identification des point pour la donner comme entrée merciii et voila la code :::

qgis.utils.iface
def lect():
     #qgis.utils.iface.addVectorLayer(file, "fichSHP", "ogr")
     canvas = qgis.utils.iface.mapCanvas()
     aLayer = canvas.currentLayer()
     print(aLayer.name())
     return aLayer
def Points():
     aLayer=lect()
     provider = aLayer.dataProvider()
     feat = QgsFeature()
     allattr=provider.attributeIndexes()
     provider.select(allattr)
     i=0
     print("le line a un type 1",QGis.Line)
     print("nombre de ligne =",provider.featureCount())
     print("nobre de champs",provider.fieldCount())
     while provider.nextFeature(feat):
        geom=feat.geometry()
        x=geom.asPolyline()
        print("la route d'id ",i)
        print("le type de la geometrie",geom.type())
        print("***************************************************************************************************")
        print("les pts constitue  la route",x)
        i=i+1
def eclate_arc(arc,indice):
     list1=[]
     list2=[]
     i=0
     while i < len(arc) :
       x=arc[i]
       list2.append(x)
       if i==indice-1 :
         list1=list2
         list2=[]
       i=i+1
     print(i)
     print("la liste originale est",arc)
     print(" 1ere et 2eme  list ",list1)
     print("la 2eme liste",list2)

[gene]

Je ne comprends pas ce que vous voulez faire

1. ce sont des lignes
2. que vous voulez découper suivant des points que vous spécifiez ?

Il y a beaucoup plus simple comme traitement que ceux qui sont indiqués dans les tutoriels de PyQgis (vous êtes sur de ne pas faire d'erreurs
La fonction suivante permet de tout sélectionner dans une couche (géométrie et attributs)

>>> def select_tout(couche):
 ...      couche.select([])
 ...      couche.setSelectedFeatures([obj.id() for obj in couche])

ensuite

 >>> macouche = qgis.utils.iface.activeLayer()
 >>> select_tout(macouche)

fait le reste, tout est dans macouche (géométries et attributs)
Si vous voulez le nom et autres choses

>>>macouche.name()
 PyQt4.QtCore.QString(u'testligne')

ensuite dans Python il vaut mieux travailler avec des boucles plutôt qu'avec la manière indiquée dans les tutoriels et celle que vous utilisez:

 >>> for elem in macouche.selectedFeatures():
 ...       geom = elem.geometry()
 ...       wkt = geom.exportToWkt()
 ...       print wkt
 LINESTRING(116935.350012 115500.583375, 117960.618676 115652.475029, 119201.067183 116184.095818, 119947.867814 116639.770780, 120897.190651 117399.229049, 121947.774591 117981.480389, 122821.151601 118411.840075, 123808.447352 118930.803226, 124226.149400 119095.352518, 124859.031291 119250.419917, 125504.570820 119393.903375, 126122.693025 119381.291081, 126746.093024 119250.510605, 127210.875871 119267.479180, 127480.348645 119110.153903, 127755.069429 118784.225242, 127693.191033 118512.405932, 127397.984044 117940.377074, 126837.736203 117219.874727, 126864.632198 116890.211110, 127285.349964 116676.704033, 127920.077596 116613.336464, 128366.270666 116608.104013, 128694.295553 116639.770780, 130466.364848 116576.482590) etc.
etc.

ou avec ce que vous voulez faire:

>>> for elem in macouche.selectedFeatures():
 ...       geom= elem.geometry()
 ...       x=geom.asPolyline()
 ...       print x  
[(116935,115501), (117961,115652), (119201,116184), (119948,116640), (120897,117399), (121948,117981), (122821,118412), (123808,118931), (124226,119095), (124859,119250), (125505,119394), (126123,119381), (126746,119251), (127211,119267), (127480,119110), (127755,118784), (127693,118512), (127398,117940), (126838,117220), (126865,116890), (127285,116677), (127920,116613), (128366,116608), (128694,116640), (130466,116576)], etc.

si vous voulez absolument votre indice i, il faut alors utiliser:

 >>> for i, elem in enumerate(macouche.selectedFeatures()):
 ...       geom= elem.geometry()
 ...       x=geom.asPolyline()
 ...       print "indice: ", i, x
 ... 
 indice: 0 [(116935,115501), (117961,115652), (119201,116184), etc.

Pour les attributs, c'est la même chose

>>> for att in macouche.selectedFeatures():
 ...      attrs = att.attributeMap()
 ...      for (k,atr) in attrs.iteritems():
 ...          print "%d: %s" % (k, atr.toString())
0, -100, test
1, -200, test2
etc.

en résumé:

def select_tout(couche):
     couche.select([])
     couche.setSelectedFeatures([obj.id() for obj in couche])
macouche = qgis.utils.iface.activeLayer()
select_tout(macouche)

for elem in macouche.selectedFeatures():
           geom= elem.geometry()
           attrs = elem.attributeMap()
          (traitements)

effectue tout votre traitement en une seule boucle simple, hormis eclate_arc(arc,indice) où je ne pige pas ce que vous voulez faire

[david12312]

merci mais j'ai besoin d'extraire les points vertex pour faire la segmentation é après determiner le plu court chemin

[nhmidi]

bjr,j'ai besoin des point vertex (intersection des ligne) pour segmenter les lignes,comment peut extraire ces points à partir de console python sous QGIS pour passer à la méthode eclate_arc l'arc et l'indice du point.

[gene]

désolé, mais vous ne pouvez pas extraire les intersections avec PyQGIS mais vous devez utiliser un module comme Shapely: voir
http://www.forumsig.org/showthread.php?t=33857

Pourquoi ne pas utiliser l'extension pgrouting qui est faite pour ça ?
http://www.portailsig.org/content/os...-openstreetmap

[nhmidi]

merci

[david12312]

merci bcp