Le portail francophone de la géomatique

Bonjour,

Me revoila avec un nouveau problème...

Je souhaite mettre en place une analyse de coûts cumulés sur le même type d’espèces, des amphibiens (tritons et crapauds), dans le style de ce qui a été fait ici:
http://georezo.net/forum/viewtopic.php?pid=263979

Je travail sur leur capacité de dispersion à partir d'un réseau de mares que j'ai cartographié. J'ai une couche raster avec mon Occ_sol avec en valeur le coefficient de friction (perméabilité) pour chaque pixel de ma zone (300km²) et un fichier de point avec toutes mes mares (425 au total). Je voudrais savoir jusqu'où peuvent aller mes tritons sachant que la valeur limite que j'utilise est de 500m autour d'un point d'eau. Donc si le milieu est favorable, avec un coef de friction faible, 1m d'effort pour 1 mètre parcouru mon triton pourra avancé de 500 m mais si le coef de friction est plus élevé, 5m d'effort pour 1 mètre parcouru par exemple alors il ne pourra faire que 100m. C'est un exemple simpliste par rapport aux valeurs de mes coef de frictions qui change souvent (et qui vont jusqu'à 10 000)


J'ai vu qu'il y avait les outils r.cost dans GRASS du coup je les ai essayé:

r.cost: demande un point d'arriver or je veux juste qu'il me calcul les coûts cumulés tout autour des points et non vers un autre point.

r.cost.full.raster et r.cost.full: le problème c'est que ça demande beaucoup de travail au PC. J'ai essayé sur toutes la zone avec tout les points, il n'a pas réussi à le terminer en 2 jours. Du coup j'ai essayé sur un seul point avec un carré de 1000 pixel autour (1pixel=1mètre). Ça a parfaitement marché. Donc après j'ai  tenté avec 7 points sur un autre carré de 1000 pixel. Lancé hier midi et ça tourne toujours, ça devrait avoir fini ce soir... Donc je me dit qu'à ce rythme la ça ne sera pas possible de le faire sur toutes ma zone avec tout mes points...
Comme la zone est grande il se retrouve à devoir cumuler des valeurs sur des zones beaucoup plus grande que mon rayon de 500m autour d'une mares et tout ça pour rien.


Est ce que quelqu'un a déjà eu à effectuer ce genre d'analyse, sur une si grande zone et avec autant de points source?
Est ce que je dois passer par python et si oui avez vous une idées de la démarche a suivre?


J'ai également publier ce post ici :
http://www.forumsig.org/showthread.php/ … post344373

Je vous tiendrais au courant de la solution si elle m'est donné.

Encore merci pour votre aide

Tristan