Le portail francophone de la géomatique

Bonjour spl33n,
Alors comme personne ne t'a répondu, je me lance. Je ne suis pas du tout expert en Qgis, et je ne l'utilise que pour visualiser rapidement des choses réaliser avec GRASS. Donc je n'en parlerais pas. De mon côté et avec les indications que tu nous as données. J'aurai effectivement commencé avec la création d'un MNT. Il faut que tu regardes de plus près les méthodes d'interpolation suivant le résultat que tu cherches a obtenir. Personnellement j'ai travaillé sur les micros reliefs, et je n'ai pas pu utiliser r.surf.rast parce que la résolution du SHP couplé a l'interpolation, je perdais mes détailles (peut être que je m'y prenais comme un pied).
Ensuite à partir du MNT je reconstruirais les basin versant et les court d'eau (très bon tuto ici http://www.ing.unitn.it/~grass/docs/tut … /dtm4.html ils sont bon dans le Trento ;-) ) et je comparerais avec mes données existantes ça permettra de valider ton MNT.
Pour les ondes de de charge des cours d'eau j'aurai utilisé un module qui permet de faire le calcule des isochrones (j'ai vu un truc ici mais j'avais lu ça en diagonal http://www.forumsig.org/showthread.php?t=25912)
Ou alors j'aurais directement utilisé un script avec une boucle et r.mapcalc ....
Voilà, j'espère, que je n'ai pas dit trop de bêtises et que ça te donnera des Pise.
Bons courages et tiens-nous au courant
Etienne

PS je suis sur qu'il y a des fondus qui ont développé des modules pour faire ce qui t'intéresse dans GRASS il y en a de partout!!

Merci pour vos réponses, je vais potasser tout ça!

Bonjour Spl33n,
je rentre juste et je vois que personne ne t'a répondu. Alors pour intégrer ton chenal je pense que c'est possible. Il doit falloir un peu bidouiller avec r.mapcalc une fois que tu l'as rastérisé. Je ne sais pas si tu as suivi le tuto de l'université de trento jusqu'au bout ou jusqu'a la création des cours d'eau. Mais  si c'est le cas, tu peux véritablement ajouter les raster avec quelque chose du genre :

r.mapcalc 'riverfinal=rivertuto+riverPolygone'

tu peux procéder de la même manière pour l'ajouter au raster... en fait le secret pour moi c'est d'arriver a désincarné les objets et de réfléchir sur des matrices de pixels qui ont des valeurs et ce sont les évolutions de ces valeurs qui nous intéresse.

Pour l'histoire du script je pense aussi que c'est une bonne idée pour voir l'évolution dans le temps avec une boucle for ou while et ton compteur qui devient véritablement un chronomètre.

Voilà en espérant que tu ne sois pas resté bloqué et que ma réponse ne serve à rien

Etienne

Bonjour Fab,
ne te lance pas dans R ni dans Python tout de suite, commence par bash :-)
pour cela il faut procéder en plusieurs étapes :
1 définir le nombre pas dans ton phénomène de crue (on est dans une sorte d'abstraction ou tu découpes les temps en fonction des étapes que tu veux connaître)
2 tu fais une boucle (ou pas selon le nombre d'opérations que tu prévois)
3 tu lance le script dans grass (je suis sous Linux donc je fais sudo chmod 777 monfichier.sh et ensuite ./monfichier.sh)
et là .... tu corriges les erreurs ou ton ordinateur mouline pour avoir les résultats.

tu as un exemple de script Grass ici : http://www.grassbook.org/scripts2nd/rsun_crop.sh.txt ce n’est pas pour ton problème, mais il y a une boucle pour t'inspirer.

Bon courage

Etienne

pourquoi un "script" ne serait pas indiqué ? "script" veut simplement dire un fichier contenant des commandes qui sont interprétées et non pas compilées avant exécution ... je ne voit pas trop ce qu'on script en python aurait de moins par rapport à du code C, à
GRASS fonctionne très bien avec Bash ou python, il est très facile d'automatiser des tâches répétitives. Il faut simplement prévoir tous les calculs à effectuer (par exemple avec r.mapcalc), et ensuite coder la boucle ...

la principale difficulté n'est pas de faire une boucle dans un script, mais bien de savoir à chaque étape comment combiner les informations pour effectuer les bons calculs !