Problem mit Fließlinien in ArcMap

Question

<\/HEAD>Hallo an die Community, hier ist George!<\/P> Im Rahmen meiner Promotion muss ich die Flüsse im nördlichen Peloponnes, Griechenland, kartieren. Ich habe zwei DEMs von Nasa Reverb heruntergeladen, die ich mit "mosaic to new raster" verbunden habe. Ich habe den üblichen Prozess von "fill", "flow direction" und "flow accumulation" gestartet, um die Flüsse der Region analysieren zu können. Leider stehe ich vor zwei Problemen.<\/P> Das erste Problem ist, dass die Ergebnisse des Flow Direction Tools falsch sind und ich keine Lösung finde, egal was ich versuche. Wie Sie im ersten Screenshot sehen können, gibt es einige horizontale Fließlinien, die keinen Sinn ergeben. Ich habe versucht, die Werte (versucht ArcMap dazu zu bringen, Werte unter 1 auszuschließen) des ursprünglichen Rasters zu bearbeiten und dann "fill" zu verwenden, aber das verhindert nicht, dass das "fill"-Tool diese Bereiche ignoriert (dritter Screenshot).<\/P> Das zweite Problem (das scheinbar mit dem ersten zusammenhängt) ist, dass ich beschlossen habe, trotzdem zum nächsten Schritt überzugehen und "Flow Accumulation" zu verwenden, da einige Bereiche in Ordnung schienen und um zu sehen was passiert.<\/P> Wie Sie im zweiten Screenshot sehen können folgen die Flussakkumulationslinien der Form der Küstenlinie und erstrecken sich dann ins Meer. Soweit ich verstehe folgen diese Flussakkumulationslinien den Flow Direction Linien, daher entsteht das Problem beim Einsatz des "flow direction" Tools.<\/P> Können Sie bitte helfen? Ich bin noch ziemlich neu darin und bitte um Nachsicht für meine Unwissenheit.<\/P> Falls es ein Problem mit den Originaldaten ist und das Problem nicht korrigiert werden kann – gibt es eine Möglichkeit dies zu umgehen?<\/P> <\/P> Vielen Dank!<\/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.