Le portail francophone de la géomatique

Bonjour,

Vous pouvez également réaliser librement ce genre de calcul avec GRASS

r.basins.fill pour bassins versants:
http://grass.itc.it/grass64/manuals/htm … .fill.html

r.flow pour les ecoulements
http://grass.itc.it/grass64/manuals/htm … .flow.html

r.slope pour le calcul de pente
http://grass.itc.it/grass64/manuals/htm … spect.html

@+

Nicolab

Bonjour,
Je pense que la solution de combler les puits de façon artificielle n'est pas la meilleure.
Il y a probablement une possibilité de changer de pas de la grille.
Comptez le nombre de cases totales et comparez au nombre de puits trouvés : 2632, je pense que c'est là que se situe le problème.
La superficie traitée est grande, mais le fichier définissent le MNT est gros ?
D'ailleurs, je ne sais pas si c'est un fichier représentant un semi de points XYZ ou si c'est un raster.
Cordialement.

Bonjour,
Le nombre de cases dont je parlais est le produit du nombre de lignes par le nombre ce colonnes.
Evidemment 53MO, c'est énorme.
Vous devriez essayer de visualiser l'image pour essayer de savoir de quelle façon se disposent les puits. Ce n'est pas du tout la même chose s'il s'agit d'un lac sur lequel il n'y a aucun point, ou au contraire des puits isolés.
Vous devriez essayer aussi de découper ce MNT de façon à pouvoir le traiter, au moins en partie, par morceaux. Autrement dit, même si vous n'obtenez, pour l'instant, qu'un résultat partiel, voire très limité, vous aurez déjà quelque-chose et ça vous permettra peut-être de trouver ce qui bloque. 
Il me revient à l'esprit un principe en géographie : sur un même continent, on peut aller de n'importe quel point vers n'importe quel point sans se mouiller les pieds. Il est très probable que le traitement des bassins versants est basé sur ce principe. Si le programme trouve un ile au milieu d'un lac, il y a de fortes chances qu'il soit perdu. A mon avis, c'est ce genre d'"anomalie" qu'il faut chercher.
Bon courage.

Vous ne pouvez joindre de fichiers de 53 Mo, cf l'aide du forum. Le plus simple est de faire une visualisation du MNT et de faire une impression écran. On pourra déjà voir si ce lac est si anormal que ça.


@Pierre, un grand lac ne pose pas de problème, on peut toujours aller à n'importe quel endroit du continent sans se mouiller les pieds, ou alors, on a droit à un canal du panama.

Bonjour,
Je suis content de voir revenir le sujet.
Comme le sait Robin, je travaille avec un logiciel que je développe moi-même. J'ai bien avancé sur cette question des bassins versants qui m'intéresse beaucoup, mais je n'ai pas terminé, d'ici 8 jours, j'espère.
A Romfly : le MNT me parait parfaitement correct, c'est d'ailleurs très joli. Je joins une copie pour le plaisir.
A Jerome, C'est vrai qu'un grand lac ne pose pas de problème (sauf s'il contient une ile et que par hasard on y met les pieds).
A Robin, Si il y a un puits, c'est une faute. Soit on constate que c'est une faute et on la corrige, soit c'est une anomalie, et ce n'est pas la peine de continuer. Un bassin versant est déterminé par son périmètre, c'est à dire une ligne de partage des eaux et son exutoire, le point bas. Autrement dit, si la ligne définissant le périmètre est fermée et si le point bas est unique, on a un bassin versant. Tout ce qui se trouve à l'intérieur n'a aucune importance (hydroliquement - sauf pour l'évaluation de temps de concentration, et encore). Tout cela pour dire que le comblement des puits me parait un mauvais choix D'ailleurs, il ne me semble pas que le fichier concerné en comporte.
Je vous tiendrai au courant, mais il faut me laisser un peu de temps.
cordialement.
Pierre

Dernière modification par Pierre Dolez (Thu 12 February 2009 12:19)

Ca c'est un très gros sujet de discussion, et qui n'a rien à voir avec un fichier raster.

Je fais une différence entre un puits (sans fond) et une dépression. Si on a un puits sans fond, je considère que c'est une faute, parce que ça n'existe pas, donc je le corrige.
Si c'est une dépression, ça ne me gène pas dans la détermination des bassins versants, parce que je n'ai pas de ligne de point bas à rechercher.
Sauf erreur de ma part, on cherche à déterminer un bassin versant pour connaître son aire, l'altitude du point le plus haut, celle du point le plus bas, la longueur du chemin hydraulique et son coefficient de ruissellement.

Sur le plan hydraulique, si on a une dépression avec ou non un lac au fond, il s'agit d'un petit incident de la nature qui ne doit empêcher en aucun cas le calcul du bassin versant qui contient cette dépression mais qui possède un exutoire différent de cette dépression.
Autrement dit, ce serait une erreur que d'essayer de la corriger. Elle existe et c'est tout.
Mais on s'éloigne des SIG

Dernière modification par Pierre Dolez (Thu 12 February 2009 15:51)

Bonjour et merci à vous 2 !!
Je vois que vous etes bien calé sur le sujet et j'ai encore 2/3 questions à vous posez.
J'ai donc générer mes BV mais je constate des anomalies, et j'ai bien l'impression qu'elles concernent essentiellement les zones d'étendues d'eau.
Je vous joins dans un .zip, 3 couches :
- BV généré depuis mon MNT (avec l'outil ArcHydro)
- Zones d'étendues d'eau (tout les lacs de la zone)
- limite administrative de la zone.

Je vous joins la limite admin car mes BV ne couvrent pas toutes la zone admin. Je pense que cela provient de la valeur seuil choisie pour générer mes BV. Pouvez vous me le confirmer et peut être m'apporter une solution.
Sinon pour résoudre mes problèmes sur les étendues d'eau, je ne sais pas si il existe une méthode qui permette de renseigner les lacs pour rectifier mes BV.
Merci pour votre aide.
Bon Robin, juste pour t'annoncer que ce w-e je vais pouvoir me baigner "dans un lac azur aux plages chaudes dans des températures tropicales au milieu d'une île" car je ne vous ai peut etre pas encore dit, que je travaille en Nouvelle Calédonie (d'où les réponses décalés dù au décalage horaire de 10h !!).
Merci en tout cas,
a+

Romain