Problème avec les lignes de flux dans ArcMap

Question

<\/HEAD>Bonjour à la communauté, George ici !<\/P> Dans le cadre de mon doctorat, je dois cartographier les rivières du nord du Péloponnèse, Grèce. J'ai téléchargé deux DEMs depuis Nasa Reverb, que j'ai reliés avec "mosaic to new raster". J'ai commencé le processus habituel de "fill", "flow direction" et "flow accumulation" afin de pouvoir analyser les rivières de la région. Mais malheureusement, je fais face à deux problèmes.<\/P> Le premier est que les résultats que j'obtiens avec l'outil Flow Direction sont erronés et je ne trouve pas la solution, peu importe ce que j'essaie. Comme vous pouvez le voir sur la première capture d'écran, il y a des lignes de flux horizontales qui n'ont aucun sens. J'ai essayé d'éditer les valeurs (j'ai essayé de faire en sorte qu'arcmap exclue les valeurs inférieures à 1) du raster original puis d'utiliser "fill" mais cela n'empêche pas l'outil "fill" d'ignorer ces zones (troisième capture d'écran).<\/P> Le deuxième problème (qui semble être lié au premier) est que j'ai décidé de passer à l'étape suivante et d'utiliser quand même "Flow Accumulation" car certaines zones semblaient correctes pour voir ce qui se passe.<\/P> Comme vous pouvez le voir sur la deuxième capture d'écran, les lignes d'accumulation de flux suivent la forme du littoral puis s'étendent en mer. D'après ce que je comprends, ces lignes d'accumulation de flux suivent celles de flow direction donc le problème est créé lorsque j'utilise l'outil "flow direction".<\/P> Pouvez-vous m'aider s'il vous plaît ? Je suis assez novice dans ce domaine donc excusez mon ignorance.<\/P> Si c'est un problème avec les données originales et que le problème ne peut pas être corrigé, y a-t-il une solution alternative ?<\/P> <\/P> Merci !<\/P> <\/P> <\/P> <\/P> <\/P> <\/P> <\/P> <\/P><\/BODY><\/HTML>

MarkBoucher · Answer

Georgios,

I found the area you are working in using Google Maps (see attached image). The area where you have the strange lines is in the water (Gulf of Corinth). If the elevation in your DEM is level in the water area, the flow direction would be somewhat random. Since you are dealing with a coastal area, you will have some unique issues to deal with. However, I don't see anything wrong with the results. There was some recent discussion on coastal drainages. Here is the link to that includes that discussion: Re: Catchment vs watershed.​

You are right that the flow accumulation is based on the flow direction. If the flow direction is wrong (has errors) the rest of the analysis will have errors.

Since you have creeks, rivers, and streams that have to cross under roads on their way to the sea, I suspect you will want to create an agreestream layer and "burn" the streams into the rawdem at the road crossings so the flow direction will be correct at those locations. It looks like the topography is steep enough that for the most part you won't need to burn in streams anywhere else except in very flat place or where man-made drainage systems control the flow direction instead of the topography.

So, basically you are not doing anything wrong to get the "horizontal flow lines". A possible way to get rid of them would be to use the Con() function to make every raster grid below a certain elevation (i.e. in the water) "Null". I'd have to do a bit or research myself to give you direction on how to do that. I'm sure there is a tool in Spatial Analyst that would do work. This link might show you how. ArcGIS Desktop Help 9.3 - Setting values to NoData: SetNull.​ Once the raster has null (or "no data") values in the water area, it is possible that the flow direction grid will not be defined there and so you won't get flow accumulation there either.

The flow direction grid is only a foundational layer and will not show in your final analysis where I assume you will define watersheds. You might simply move on and establish your analysis process for a small area until you settle on a process you are happy with. Then perform that process on the whole dataset.  I use Model Builder in this way. I allows me to build and edit the process and they reapply it on different datasets.

Best,

Mark

GabrielUpchurch1 · Answer

Hi George,

It seems to be something related to the data but that doesn't necessarily mean the data is bad.  It may just take some preprocessing to get it to work properly.  I have never observed this behavior before but let's see what we can figure out.

For starters, what specifically do you mean by:  I tried editing the values (tried to make arcmap exclude values lower than 1) of the original raster  My first thought would be to eliminate the "sea" area from the data so I am curious what method you used when attempting this.