1

My main goal here is to end up with a shapefile that simply defines the ridgeline in a given area. I followed the directions at How To: Identify ridgelines from a DEM and am having some trouble. Firstly I went through every step exactly as it is shown on that page, and ended up with a shapefile that defined almost every cell in the raster as "ridgeline." I went back through in order to diagnose where I went wrong. What I found was that in the step where you invert your filled DEM, the resulting raster has all negative values for elevation. I modified the equation in the raster calculator to be:

(("Rasterlayer" - Max_value) * -1) + Min_value

This made it so that my inverted and filled DEM had positive values. I threw the two DEMs (inverted and not inverted) into arcscene just to check my sanity and yes, the new inverted DEM showed valleys where previously there were ridges. All made sense.

However, when I then move on to the flow direction tool, I get weird results. I specify the new inverted DEM as the input, and choose to save the output as the same name plus "flow_dir" somewhere in the name. I checked the box to make the edge cells flow outward, and I left it at D8, which is what I wanted. In theory, this should give me a raster where the cells are defined as one of 8 discrete values, 1, 2, 4, 8, 16, 32, 64, & 128. In practice, however, I get a raster that has not 8 discrete cell values but 208 discrete values, from 1-208. What have done wrong here and how do I fix it?enter image description here

It should be of note that when I run the original filled DEM (non-inverted) through the flow direction tool with the exact same parameters, it spits out a raster that has the expected 8 discrete values. The screen grab here shows the raster that I get with 208 values (Flow_dir_7) which is not what I want but was resultant of using the tool on the inverted DEM. As you can see, right below that is the Flow_Dir_5 raster which is defined by the 8 values as it should be for D8, however this raster was produced using the non-inverted DEM as the input.

Do you know where I went wrong or know another way to define ridgelines other than inverting my DEM and pretending they are streams?

I'm running Windows 10 and ArcMap 10.7.1.

2 Answers 2

1

If you follow exact steps from ESRI, it works and will produce something like that:

enter image description here

I don't like that they are disconnected. Perhaps create sub-catchments from stream network and assume that divides between them ARE ridges:

enter image description here

Note smaller flow accumulation threshold you'll use to derive streams will result in more detailed ridges network.

Consider erasing part of them near streams.

You might try to produce same size subcatchments:

enter image description here

but this is tricky.

2
  • Hi FelixIP, thanks for your answer. I just now went back through and followed the instructions again word for word. I still came up with the same result; a "stretched" raster rather than one with unique values. I decided to follow through the rest of the steps to output a shapefile, and what I got was not a shapefile that defined ridges, but one that haphazardly define almost the entire DEM extent. Would you be able to share screengrabs of the settings you used at each step? I really am unsure at which point in the process I am going astray. Jan 23, 2020 at 4:51
  • Not sure I can make it more clear, sorry. Perhaps increase flow accumulation threshold and skip stream order. Just derive them using flow direction and raster after threshold.
    – FelixIP
    Jan 23, 2020 at 20:40
0

If you want to have exactly a flow direction raster with 8 numbers from 1 to 128, you should use the "fill" tool for your inverted DEM, then use the resulting raster as input for the "flow direction" tool. I had the same problem this morning! and just figured it out!

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge that you have read and understand our privacy policy and code of conduct.

Not the answer you're looking for? Browse other questions tagged or ask your own question.