Annonce
#1 Fri 07 August 2009 18:10
- dominique.lys
- Participant assidu
- Date d'inscription: 5 Oct 2006
- Messages: 473
- Site web
QGIS/Grass : calcul d'isochrones
Hello,
J'essaye de faire une analyse de réseau avec GRASS via QGIS pour calculer les isolignes (cost distance).
Voilà, l'enchainement que je fait (pour le moment je ne cherche pas à utiliser des données attributaires) :
-v.in.ogr : import couche réseau et couche ponctuelle
-v.patch : assembler le réseau et les points dans la même couche
-v.distance : créer des lignes reliant les points au réseau
-v.patch : assembler les lignes précédemment crées à la couche réseau + points
-v.clean.break et v.clean.snap : nettoyer le réseau
-v.net.iso : créer les isolignes
Le problème c'est que je ne comprend pas trop ce que devrait me produire v.net.iso, j'imaginais au départ obtenir une nouvelle donnée attributaire m'indiquant à quelle classe de coût appartient chaque polyligne, mais je n'ai rien d'exploitable.
Voici les éléments d'aide que j'ai pu trouvé :
http://www.gdf-hannover.de/lit_html/gra … ode51.html
http://www.hpcc.nectec.or.th/grass/gras … twork.html
http://www.gdf-hannover.de/lit_html/gra … ode57.html
http://grass.itc.it/grass64/manuals/htm … t.iso.html
Dernière modification par dominique.lys (Fri 07 August 2009 18:20)
Hors ligne
#2 Tue 11 August 2009 07:21
- Sylvain_M
- Membre
- Lieu: Lyon
- Date d'inscription: 3 Aug 2009
- Messages: 298
Re: QGIS/Grass : calcul d'isochrones
d'après la doc, v.net.iso génère un fichier vecteur (output=name Name for output vector map). regarde en bas de la page (http://grass.itc.it/grass64/manuals/htm … t.iso.html) pour voir l'exemple basé sur les données spearfish.
à priori, il suffit de faire "d.vect result_net_iso" pour afficher la couche vecteur sans affichage particulier...
Hors ligne
#3 Mon 17 August 2009 10:10
- dominique.lys
- Participant assidu
- Date d'inscription: 5 Oct 2006
- Messages: 473
- Site web
Re: QGIS/Grass : calcul d'isochrones
Finalement j'ai trouver la solution,
le problème était que je n'obtenais en sortie qu'un seul objet ligne et aucune donnée attributaire exploitable.
En fait après avoir assemblé les lignes reliant les points au réseau avec la couche réseau + points on obtient deux 'sous couches' de type ligne car les nouveaux objets intégrés dans la couche n'ont pas d'identifiant. Il faut donc faire 'v.category.add' afin de lier toute les objets linéaire dans la même couche. Les opérations de nettoyage porterons alors sur l'ensemble des lignes (auparavant les nœuds n'étaient pas correctement reliés car les objets sans 'category' étaient exclus du traitement).
En sortie, la colonne 'cat' de la couche vectorielle obtenue indique la classe de cout.
Hors ligne
