nav[aria-label="Primary Navigation"] { padding: 0; & ul { list-style: none; width: 100%; display: flex; flex-direction: row; justify-content: start; align-items: start; gap: 30px; padding: 0; & li { margin: 0; } & ul li { list-style: none; } } }

Stream Order (Strahler) and Flow Direction

PeterWilson

I have generated a terrain dataset for South Africa from 20m contours and SPOT shots. In order to eliminate errors in the terrain where there isn't enough data especially low lying areas. I would like to make use of river centerlines to enforce the drainage within the terrain dataset. My proposed methodology is to intersect the river centerlines with the contours and extract the Z values of the contours into the river centerlines. I would then use the polyline Z river centerlines as hard break lines within the terrain dataset.

The problem that I have is that I need to ensure that the flow direction of the existing river centerlines that I wish to use is correct and secondly need to assign a stream order (Strahler) to the river centerline segments.

I've attached a screen shot of the river centerlines to help explain the problem that I have with the currently available river centerline dataset. There are numerous segments that should be tributaries but do not extended enough to link up with the main river. I'm looking for advice in how I can make use of the available river centerlines to embed into my terrain. Do I eliminate the segments that are not connected to the main river centerlines.

What is the most efficient methodology to rectify the river centerlines to a point that they won't add additional error to my terrain dataset, which I then would like to use to generate a river centerline network to build up an hydro data model. Any assistance will be appreciated.

Find more posts tagged with

Arc Hydro

Accepted answers

All comments

PonderosaAdvisors

Great! That works for the most part. I am using HND streams and am having trouble with all of the streams being assigned a river order. I am attaching a gdb with a sample. Here is also a clip of the message that I receive saying that 44 of the segments did not receive a stream order. If you could direct me to a solution I would greatly appreciate it.

StreamOrder.gdb.zip

Message From Results.png

Tables 3 and 4 from ms2525b.pdf

Tables 3-1 thru 3 from ms2525c.pdf

ChristineDartiguenave

StreamOrder uses flow direction as input, i.e. flow can only go one way so braided network are not supported.

If you have a braided network you can use the Assign River Order tool from the Arc Hydro toolbox or Arc Hydro Python toolbox. It is located in the Attribute tool toolset. This tool takes as optional input a flow split table that allows the tool to handle braids.

Coming in with your own streams feature class you will need to run the following tools:

· Assign HydroID

· Generate From/To Node for Lines (store connectivity)

· Find NextDownstream Line (to create the flow split table and NextDownIDs)

· Assign River Order.

Also the field OmitUp (0/1) if present is used to reset the order to 1 in a secondary braid segment.

Christine Dartiguenave

Esri Water Resources Team

PonderosaAdvisors

I am having the same issue when I use my own streams which have braids and loops. If I let ArcHydro define the streams for me then it works just fine. The jpg that you attached doesn't show any braided streams, but I am wondering if your data did contain braids. I think the Stream Order cannot handle braided networks. Can anyone confirm this?

SandeepMahajan1

Hi Aslam can you please let me know how you solved your problem? because I also want to generate drainage order on a digitized drainage network?

Thanks & Regards

Sandeep Mahajan

TRINY_GRACE_LORETLICERA

Thank you for your help Sir Mark. I already figured it out

MarkBoucher

Grace,

I haven't used HEC-GeoHMS for a while, but I'm very familiar with Arc Hydro. GeoHMS uses the same logic and tools as Arc Hydro for processing the watersheds. What GeoHMS does is take the processed grid and vector data and then allows the creation of an HMS model from that data.*

Normally, when using HEC-GeoHMS (and Arc Hydro), you start with a digital elevation model (DEM) and then process it through a few steps to create a flow direction from which you derive the rivers (drainage lines). I have started the practice of not digitizing or using existing river/stream polylines until I need them for correcting the DEM where the areas is flat or the DEM is not defined enough to put river/stream in the right place. In these cases I use the agreestream polyline layer to put the river/stream in the right place and then burn them in. See this related discussion for more on my approach to burning in streams: Stream Definition threshold value.

It looks like there is also some discussion about burning in streams in this thread. Wow! It started in 2011!!

Mark

*Footnote: The first generation GeoHMS tools relied on Arc Hydro tools to process the DEM before using the GeoHMS tools. Now, it appears the GeoHMS toolbar has the Arc Hydro tools needed to get to GeoHMS datasets that can be then used to build the HMS model. Likely, you would only need the GeoHMS toolbar if your goal was to create an HMS model from the GIS data.

TRINY_GRACE_LORETLICERA

Hi. Were you able to solve the problem? I need some help here. I digitized a river network out from georeferenced topographic map. Meaning the polyline river has no data. I'll be using the river for HEC-GeoHMS. But upon processing when I tried the command adjointcatchment it says that my river polyline has no GRID CODE/ID? How can I assigned GRID ID to my river polyline? I need to produce a Strahler stream order polyline? Kindly show me how to do it.

My email address is tglicera@gmail.com

Thanks,

Grace

SeoyoungKim1

Hi Nicole,

Did you solve the problem about producing the Strahler stream order using the stream shapefiles?

I also have same problems. Would you show me your solution if you solved?

My email is

sykim2012@knights.ucf.edu

Seoyoung

AslamHanief

Got it. Thanks.

AslamHanief

Hi,

I'm unable to locate the Assign River Order toolset under the Attribute Tools. Has it been removed in ArcHydro 2?

AslamHanief

Hi Christine,

I've seen you referring to the Assign River Order toolset; however, I'm unable to locate the Assign River Order toolset under the Attribute Tools. Has it been removed in ArcHydro 2? Please see image.

Thanks,

Aslam.

Attribute Tools.png

NicoleCarlozo

Thank you for trying - I appreciate it!

When I burned the stream into the DEM, I made sure the stream value was -1000 with a buffer of -100. I used the stream raster of value 1 as the input into the stream order tool.

I think the stream raster sometimes looks like there are 2 cells because the polyline curves or runs along the edge of the raster cell. I'm not sure how to get around that.

StreamOrderTest2.jpg

CE _to_UE426.mp4

MarkBoucher

I tried using the SA Hydrology tools and Raster Calculator to do what you are doing. When I converted the stream polylines to raster I got a single cell pattern for the streams raster. I was not able to get all the tools to give me results I was expecting and then someone else doing "real work" needed the SA license and I had to stop investigating.

For your stream raster, you may need a value of greater than 1 to make the burn work. Use the Con() function to create another stream raster where 1 is changed to -1000, then add the two rasters. This will definitely burn the streams deep.

NicoleCarlozo

Thanks for the clarification. I tried what you suggested, but the result did not change. When creating the stream raster, I made sure the cell size matched the DEM. I also tried using the default cell size and a very small cell size. The same problem occurred. In environmental settings, I tried snapping the stream raster to the DEM during processing, as well as not snapping. Still no change.

My workflow:

1) Polyline to Raster conversion (1 on a background of NoData)
2) Burn the polylines into the DEM.
2) Fill the burned DEM
3) Flow Direction
4) Use the stream raster and flow direction outputs in the Stream Order tool.

Could this problem be due to a high resolution DEM? Since I am using a buffer (-100) around the stream lines, I'd think that the flow path/direction would not be a problem.

MarkBoucher

If you use the Stream Link tool you should get a raster with just a stream defined by single grid cells (one cell wide). See example attached. Your image looks like the stream is defined by several cells wide. Your image is not an example of what I would call "parallel streams".

In my image, if you look at the confluence to the right of the center of the image, you will see parallel streams defined for a short distance (mine is a 10'x10' grid). They are only one cell apart. If I were to redo this one, I would add agreestream features to burn the DEM so the parallel streams merge (that is if I knew they should not be parallel and should come together sooner).

strlnk.jpg

DataCheckReport OLD.py.zip

NicoleCarlozo

Thanks Mark. The SA Process you described makes sense, but when I followed it, the output did not make sense. Within a branch, the stream order values fluctuated (see attached for an example/screen shot).

Is this what you meant by parallel streams? I did use a small buffer as you suggested, but I could make it wider. Within the stream order tool, I used a raster version of the stream polyline and the flow direction output.

Is there an easier way to perform this analysis within ArcHydro?

StreamOrderTest.jpg

SelectMachine1.png

MarkBoucher

I normally use Arc Hydro that has several handy tools. You can use the Spatial Analyst (SA) tools but the process is not as well set up for someone doing it just once. This is a good link to review if you are using the SA tools. http://help.arcgis.com/en/arcgisdesktop/10.0/help/index.html#/Understanding_drainage_systems/009z0000005m000000/

This link http://help.arcgis.com/en/arcgisdesktop/10.0/help/index.html#/Deriving_runoff_characteristics/009z0000005p000000/ is a sub-section of the one above and shows a good flow chart of the process.

You can force the drainage lines to be at our polylines by "burning" them into the DEM (raster). You can do this by converting the polylines to a raster with a very low number (say -1000). Then add the DEM to the stream raster. The result is a DEM with the streams obviously below the surrounding terrain. Then you can run the fill tool and then the flow direction tool and this will set you up to create the stream features to which you can add the stream order because our flow direction grid will match the streams.

Arc Hydro adds a "buffer" the stream burning to ensure that the cells near the burned stream are lowered and you don't get parallel streams right next to the stream defined by your stream polyline based raster. This can happen if the polyline does not line up with the stream in the DEM. I'm sure you could figure out a way to lower the cells around the stream. Possibly you could create a small buffer around your stream polylines and then convert that to a raster with a smaller low number (-100). Then add the buffer based raster to the DEM and THEN the polyline based raster to that resulting DEM. This would take care of some parallel streams.

Some of this process will be influenced by the size of your DEM raster cells. If they are large, I don't think you will have many parallel streams. If the grid is small, you could have more.

Hope this isn't to confusing.

Best,

NicoleCarlozo

Thanks Mark,

The lidar is in a raster format. I am having problems with the lidar results due to the ditching/artifical paths that are abundant in my study site. The Hydrology Toolbox first output "streams" in the middle of ag fields. I changed the CON value to produce an output more in line with my polyline file, but I still had the same problem. Now I am circling back to assigning stream order values to the polyline file I already have. If my study area was at a smaller scale, I think the lidar method could be used. But I am looking at county-wide data.

Does any initial processing need to be performed before running the "Create Drainage Line Structures" or "Assign River Order" tools? I feel like I'm missing something, but I haven't located any documentation that steps me through the process.

MarkBoucher

Should I be using the lidar data as the raw DEM instead of the flow direction or flow accumulation grids? Do I need to specify file type?

The LiDAR data (points, I assume) must be converted to raster data. The rawdem is a raster file not a point file.

My experience is most polyline data fails to always follow the creek/river/stream (stream) idenditified via raster DEMs. This is because they are likely created by tracing over old maps that inaccuratly depict the streams or over ariel photos where the stream is hidden by trees. So, I use the raster with Arc Hydro to determine where the creeks are w/o the streams layer. Then I correct that where there stream is vauge in the DEM and where there are bridges/culverts, etc. that "mask" the creeks. This is where I burn in stream with the DEM reconditioning.

I suppose if you wanted to keep your stream polyline layer, you could do the Arc Hydro process create the stream order field inthe drainage layer and then transfer the stream order to your stream polyline layer.

NicoleCarlozo

After attempting to create stream order attributes via Lidar processing using the Hydrology Toolbox, I have circled back to this. Due to low coastal elevations and landuse changes, the lidar output isn't lining up with the stream polyline data that I already have.

I have a stream polyline, but I need stream order attributes. I have lidar, which I have used to create flow direction and flow accumulation grids. I attempted to use the "Create Drainage Line Structures" tool, however I received an error due to missing HydroIDs. I used the Assign HydroID tool within the Attribute Toolbox and tried again. It is still failing to execute. Am I missing something? Should I be using the lidar data as the raw DEM instead of the flow direction or flow accumulation grids? Do I need to specify file type?

If anyone has any experience with identifying stream order using known streams and lidar, I would really appreciate the help.

Thanks,
Nicole

D_Hornby_ddh_uni

There is a more efficient Strahler ordering algorithm created by Gleyzer et al 2004 which is implemented in an ESRI 10.1 AddIn tool called RivEX, worth a look?

ZichuanYe

the Create Drainage Line Structures function needs a DEM for it's cellsize and as a snapgrid.
In that sense, any DEM that provide correct cellsize and coordinate snapping will work.

NicoleCarlozo

Mark - I see the SA tool, but it requires a flow direction raster and a stream raster. I don't see a "Flow Direction with Streams" tool. There is a "Flow Direction Tool" that requires a hydro DEM, but I only have a stream polyline/shapefile.

Zye - The "Create Drainage Line Structures" tool requires a Raw DEM. Do you have any more information about the DEM specifics? Does resolution matter?

ZichuanYe

(1) to to assign river order, you can use Assign River Order tool of ArcHydroToolbox (Toolboxes\System Toolboxes\Arc Hydro Tools.tbx\Attribute Tools). This GP tool handles multiple river system at the same time, so you don't need to delete unconnected river lines.
(2) to make sure a known river line is considered in the flowdirection processing (burn-in river line), you can use the Create Drainage Line Structures tool (under Toolboxes\System Toolboxes\Arc Hydro Tools.tbx\Terrain Preprocessing\Create Drainage Line Structures)

MarkBoucher

There is one in SA:
Toolboxes\System Toolboxes\Spatial Analyst Tools.tbx\Hydrology\Stream Order

NicoleCarlozo

There is a tool in the Arc Hydro Tools toolbox > Attribute Tools toolset called Assign River Order that will calculates the stream order using 3 different methods

I am interested in assigning river order to the NHD dataset. I just downloaded Arc Hydro for Arc 10.0, but I do not see the "Assign River Order" tool. I do see a tool for assigning a HydroID. Doe this need tobe used before determining river order, and where is the River Order tool located?

Thanks,
Nicole

ChristineDartiguenave

There is a tool in the Arc Hydro Tools toolbox > Attribute Tools toolset called Assign River Order that will calculates the stream order using 3 different methods:

Strahler: if different upstream orders, then maximum of the upstream orders, if more than 1 with maximum order, then max order +1.
Shreve: add upstream orders.
PU_Order: downstream feature river order = maximum of upstream features river order + 1 (PU stands for Processing Unit).

This tool uses the attributes HydroID/NextDownID in the input feature class to determine the downstream direction. It does not support flow splits or loops at this time.

Christine Dartiguenave
Esri Water Resources Team

MarkBoucher

Spatial Analyst>Hydrology tools are able to create a stream order. They start with a stream raster AND you need a flow direction grid. So, if all you have is polylines of streams, you will need to make a raster with them. This should be an easy polyline to raster function. Once you have the stream raster, run stream order tool to get another raster that depicts the stream order. Then you need to convert the stream to a feature to get back in to vector data where you can symbolize the streams based on the stream order (grid code in resulting polyline feature class).

If you don't have a flow direction grid, or need to make sure you have one that exactly matches your stream polylines, you can burn the flow direction of the streams into a raster using Arc Hydro Toolbar>Terrain Processing>Flow Direction with Streams. I believe the direction of the polyline is key to making this work right. To check the direction of the polylines, change the symbology to show arrows on the lines. The arrows will point in a direction based on the direction the polylines are drawn. There is a tool to called Flip Line under Editing tools that allows you to flip the lines that don't point in the right direction.

Also, see this recent thread: http://forums.arcgis.com/threads/75629-Redrawing-rivers-in-proportion-to-their-discharge?highlight=stream+order

JanSlaats

Peter,
Did you develop any code to create a Strahler Stream Order attribute
for vector streams?
Just wondering.
Thanks,

PeterWilson

Hi Biniam

Firstly I must apologise for not getting back to you sooner. I'll gladly have a look at the algorithm, please can you email me a copy of the pdf as I don't have access to Wiley's site.

Regards

BiniamTes

hi peter ,

Though I did not see any reply to you post , I am sure u must have already solved your problem. I have somehow similar problem with your case.
I am using ArcGIS 10 for sewer network analysis. and i want to adopt river network ordering method(strahler order) to my sewer network(shape file). the sewer network is characterized by branched and mesh type. therefore, for ordering methods, i have to apply the ordering method for braided networks. I have found one algorithm below ....

http://onlinelibrary.wiley.com/doi/10.1111/j.1752-1688.1990.tb01432.x/pdf

that can exactly do the same for my network. can u give me little idea how to convert the algorithm to vba code.
thank you in advance for your assistance.

Join the Conversation

Signed in members can post, follow updates, and more. New here? Register a free account.

Quick Links

Community Help
Find useful guides, FAQs, and documents to help you navigate and make the most of Esri Community.