Topografik Islaklık İndeksi şu şekilde ifade edilebilir:
Ln(a/tanB) based on the idea of Beven and Kirkby (1979)
nerede
a is the specific catchment area (a=A/L, catchment area (A)divided by contour length(L))
ve
tanB is the slope
Buradaki temel fikir basittir, ancak hem a hem de tanB'yi hesaplamanın birden fazla yolu olduğundan , bir TWI'nin sonuçları büyük ölçüde değişebilir (Qin ve ark. 2011).
Akış birikimi ve toplama alanı, örneğin şu şekilde hesaplanabilir:
D8 (O'Callaghan, J.F. / Mark, D.M. (1984))
D-infinity (Tarboton, D.G. (1997)
Triangular Multiple flow direction (Seibert, J. / McGlynn, B. (2007)
ve diğer birçok algoritma da mevcuttur.
Eğim genellikle pikselin etrafında yerel eğim olarak hesaplanır (Sorensen ve ark. 2005). Lokal eğim, piksel etrafında minimum, ortalama ve maksimum eğim olarak da hesaplanabilir. Eğimi hesaplamanın başka bir yolu Hjerdt ve ark. 2004'te eğim hücre merkezinin d metre altındaki bir noktaya hesaplanmıştır .
Eğim, CBS yazılımlarının çoğunda temel bir araçtır, ancak hesaplama farklılık gösterebilir. Birkaç örnek: ESRI: http://webhelp.esri.com/arcgisdesktop/9.2/index.cfm?TopicName=Calculating_slope SAGA: http://sourceforge.net/apps/trac/saga-gis/wiki/Terrain% 20Analysis% 20-% 20Morphometry% 20module% 20library
Gördüğünüz gibi, hem a hem de tanB'yi hesaplamak için birçok seçenek vardır. Yani, soru, pratikte, bu algoritmaları farklı kullanarak TWI'yi hesaplamanın uygun (en iyi) yolu nedir? Yoksa var mı?
Şahsen SAGA'da çalışmayı seviyorum, çünkü çok çeşitli açık kaynaklı hidroloji araçları var.
Ps Saga CBS'de havza eğiminin tam olarak nasıl hesaplandığını ve tam olarak burada ne anlama geldiğini öğrenmek için zor zamanlar yaşıyorum . (Arazi analizi -hidroloji: su toplama alanı paralel).
DÜZENLENEN: SAGA Forumları'ndan Volker Wichmann tarafından yanıtlanmıştır: "Yakalama Alanı (Paralel) modülünün havza eğim çıktı ızgarası şu şekilde hesaplanır: her hücre için yerel eğim Zevenbergen & Thorne yaklaşımı kullanılarak hesaplanır. Son olarak, her hücre için biriken eğim değerleri hücrenin türetilmiş tutma alanına bölünür. Izgara birimi radyanlardır. "
"Topografik Islaklık Endeksi (TWI) modülü, giriş olarak normal bir eğim ızgarası gerektirir."
Referanslar:
Beven ve Kirkby 1979. Havza hidrolojisinin fiziksel temelli değişken katkıda bulunan alan modeli. Hidrolojik Bilimler Bülteni, 24, s.43-69.
Hjerdt ve diğ. 2004. Yerel drenajdaki aşağı yönlü kontrolleri ölçmek için yeni bir topografik endeks. Su Kaynakları Araştırması , 40, W05602, doi: 10.1029 / 2004WR003130.
O'Callaghan, JF ve Mark, DM 1984. Dijital yükseklik verilerinden drenaj ağlarının çıkarılması. Bilgisayarla Görü, Grafik ve Görüntü İşleme , 28: 323-344
Qin ve diğ. 2011. Maksimum iniş eğimi gradyanına dayalı topografik ıslaklık indeksinin hesaplanmasına bir yaklaşım. Hassas Tarım 12: 32-43.
Seibert, J. ve McGlynn, B. 2007. Izgaralı dijital yükseklik modellerinden yukarı eğimli alanların hesaplanması için yeni bir üçgen çoklu akış yönü algoritması, Water Ressources Research , Vol. 43, W04501
Sorensen ve diğ. 2005. Topografik ıslaklık indeksinin hesaplanmasında: saha gözlemlerine göre farklı yöntemlerin değerlendirilmesi. Hidrolojik. Dünya Sys. Sci. Tartışır. , 2, 1807-1834
Tarboton, DG 1997. Izgara dijital yükseklik modellerinde akış yönleri ve yukarı eğim alanlarının belirlenmesi için yeni bir yöntem, Water Ressources Research , Cilt.33, No.2, s.309-319