Le portail francophone de la géomatique

Bonjour à tous,

Je travaille actuellement pour un institut de recherche dans le domaine de l'hydrologie.
Mon problème est le suivant :

Je dispose d'une couche de points, et d'une couche (polygones) de parcelles agricoles. Mon but est de calculer (et de stocker dans la table attributaire) pour chacun de ces points la somme des surfaces agricoles situées à l'amont.

Pour cela, je pensais mettre en place les étapes suivantes dans un script Python (pour boucler sur tous les points):

_ Déterminer le bassin versant amont d'un point (avec r.water.outlet)
_ Vectoriser ce bassin versant (r.to.vect)
_ Sélectionner l'ensemble des parcelles agricoles incluses dans ce bassin versant (v.select)
_ Calculer la somme des surfaces des parcelles ainsi sélectionnées et STOCKER la valeur obtenue quelque part (c'est là que ça bloque!... la somme peut s'obtenir avec v.report par exemple)
_ Récupérer la valeur obtenue et l'associer à mon point initial (cette étape dépend de la façon dont a été stockée la somme calculée précédemment... je pensais à quelque chose comme v.what.vect par exemple)


Voilà, j'espère m'être fait à peu près comprendre!

Merci d'avance,

Pierre P