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; } } }

Topographic Wetness Index

StefanieSteinbach

Hi all,
I am trying to calculate the topographic wetness index for an area and know the formula, but somehow get an error everytime I try. There are so many different approaches on the internet that I am not sure which to use and why it does not work for my project.
Can anyone tell me exactly how to use my layers with the formula?
Your help is very much appreciated.

Najděte více příspěvků označených jako

Accepted answers

JeffreyEvans

It may be as simple as a missing space around "*" in your syntax.
Try: Ln((???FLOWACC??? * 900) / Tan(???SLOPE???))

You are missing a few critical steps. I would recommend looking at this thread
http://forums.arcgis.com/threads/32896-from-DEM-to-CTI?highlight=CTI

I also have this model available in our gradient toolbox.
http://evansmurphy.wix.com/evansspatial#!arcgis-gradient-metrics-toolbox/crro

All comments

DanPatterson_Retired

give a code example of what you have tried

StefanieSteinbach

I used this one below. My cell size is 30 m (an ASTER DEM). Slope is in degrees, I have also tried the formula with Tan("SLOPE")*�?�/180 as seen in many examples.

Ln((�??FLOWACC�?�*900) / Tan(�??SLOPE�?�))

DanPatterson_Retired

is slope supposed to be in degrees or radians? are you getting a specific error message?
did you go through these posts?
http://forums.arcgis.com/search.php?searchid=1784175

StefanieSteinbach

I attached the error message. Slope is in degrees, but I have read it is supposed to be in radian so I used slope * 0.0175
I checked the forums but couldn't find anything useful for my problem. If it is something simple I apologize. I am relatively new to this.

Error.PNG

JeffreyEvans

StefanieSteinbach

Thanks a million, your method seems to work out fine for my project.
Why do you use *cellsize for the upper contributing area whereas I found *cellsize^2 in other threads?

JeffreyEvans

Here is a nice explanation provided by Tom Vanwalleghem to a similar thread on the GIS Stack Exchange.

"The units of upslope contributing area are m2/m or m (in fact A is As or specific contributing area, contributing area per unit contour width) So that's why you multiply by cell resolution = 30 m and not by cell area (30*30) A more correct way to calculate this would be to take into account the flow direction (if it is diagonal or perpendicular), but for most applications, the approximated calculus given here is enough."

toddsams

Here is a nice explanation provided by Tom Vanwalleghem to a similar thread on the GIS Stack Exchange.

"The units of upslope contributing area are m2/m or m (in fact A is As or specific contributing area, contributing area per unit contour width) So that's why you multiply by cell resolution = 30 m and not by cell area (30*30) A more correct way to calculate this would be to take into account the flow direction (if it is diagonal or perpendicular), but for most applications, the approximated calculus given here is enough."

I am not understanding why the fac grid would be multiplied by 30 rather than 900. Why would upslope contributing "area" be in units of m2/m or m, which are units of "length"?

Připojte se ke konverzaci

Přihlášení členové mohou přispívat, sledovat aktualizace a další. Jste tu noví? Zaregistrujte si bezplatný účet.

Rychlé odkazy

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