Daha çox

SAGA GIS-də mənfi Topoqrafik Nəmlik İndeksi (TWI) dəyərlərini əldə etmək?


SAGA GIS-də mənfi topoqrafik indeks nəmlik dəyərləri alıram - hər kəs bunu alır?

Bütün mənfi dəyərlər deyil, araşdırdığımız bəzi dik su hövzələri üçün TWI aralığının minimumu mənfi olur. Aralar (təqribən) -1.5 qədər mənfi olur.

Bunun reallaşması üçün (xüsusi yığılma sahəsi / qaralma meylindən) yaranan dəyər 1-dən az olmalıdır, çünki 0-dan 1-ə qədər olan ln = mənfi bir dəyərdir.

Çox dik bir su hövzəsidir, buna görə 1.57 radiana (demək olar ki, 90 dərəcə) yaxınlaşan yamaclarımız var. Əslində sözügedən su hövzəsi üçün maksimum yamac 1.509 radian (və ya 86.459 dərəcə) dir. Minimum Xüsusi su yığma sahəsi dəyəri 2.2175m-dir. Beləliklə, daha da sadələşdirilmiş ln (2.2175 / 16.16159) və sonra ln (0.137208) olan ln (2.2175 / tan 1.509) hesablasaq, nəzəri mütləq minimumu göstərən mənfi bir dəyər (-1.8 və ya başqa bir şey) əldə edirik. Su hövzəmiz üçün TI dəyəri. Nəzəri səbəbi ən dik yamacın və ən kiçik SCA-nın birləşməməsi, lakin bəzi yaxın dəyərlər həqiqi minimumu əldə etmək üçün birləşdirir. -1.5.

Daha dəqiqləşdirmək üçün dik su hövzəmizdə DEM-in 8.86 m qətnaməsi var və SAGA tərəfindən hesablanan minimum tutma sahəsi bu hüceyrə ölçüsü kvadrat şəklindədir (78.4996 kvadrat metr). Bu, məntəqədə başqa heç bir hüceyrədən axını qəbul etməyən və buna görə də hüceyrə sahəsinə bərabər olan ümumi tutma sahəsinə sahib olan hüceyrələr olduğu üçün məntiqli olur. İndi bu hüceyrələr üçün, su hövzəsindəki lokallaşdırılmış zirvələrdə (və ya hətta kənarlarda maksimum zirvələrdə), Xüsusi Tutulma Sahəsi, bu hüceyrələr üçün enmə axın yollarının kontur genişliyinə bölünən ümumi tutma sahəsinə bərabər olmalıdır. Bu zirvələr bütün qonşu hüceyrələrə axır, bu 78.4996 (hazırda SAGA-da kodlandığı kimi) 35.4 ilə bölünür (8 * 0.5 * 8.86). Bu, Xüsusi Tutulma Sahəsini minimum 2.2175m verir.

Bu mənfi TI nəticəsi ilə bağlı əsas problemim bunun mümkün olmadığını düşünməyim və ya heç olmasa mənfi TI dəyərlərinin bir kağızda danışıldığını / qrafiki görmədiyini görməməmdir.

Hər hansı bir düşüncə kimisə formalaşdırır?


Bu yaxşı bir sual və məndən zaman zaman soruşulmağa meylli olduğum bir sual. Birincisi, qeyd etdiyiniz kimi, TWI = ln (a / qaralma (B)), harada a 'spesifik' su yığma sahəsidir (yəni kontur uzunluğu ölçüsü üçün normallaşdırılmış eniş zonası) və B qəfəs hüceyrəsindəki radianslarla yamac gradiyanıdır. Düzgün qeyd etdiyiniz kimi TWI yalnız mənfi dəyərləri harada alacaq a B) çünki 1-dən kiçik hər hansı bir dəyərin təbii logu mənfi olacaqdır. Beləliklə, mənfi TWI dəyərləri təəccüblü olmalıdır. Unutmayın, TWI yalnız bir indeksdir. Yalnız sizə bənzər TWI dəyərlərinə sahib olan iki qəfəs hüceyrəsinin, birləşdirilmiş qatqı sahələri və yerli yamac gradiyanları nəzərə alınmaqla, oxşar hövzə nəm şəraitində səthdə doymuş olacağıdır. Üstəlik, daha yüksək dəyərlərə sahib olan hüceyrələrin aşağı dəyərlərə sahib olan hüceyrələrdən əvvəl doymuş olacağı üçün, qəfəs hüceyrələrinin nə vaxt doymuş olacağına dair nisbi bir reytinq verir. Beləliklə, bu mənfi TWI dəyərləri həqiqətən hövzədəki ən quru yerlər arasındadır. Beven və Kirkby'nin TWI-nin əsl dahisi budur ki, bunu hövzənin mövcud nəm şərtlərini izləyən sadə bir hidroloji modelə (məsələn, TOPMODEL) bağlasanız, səth doymuş sahənin cari dərəcəsini bir eşiklə əlaqələndirə bilərsiniz. TWI dəyəri. Bu eşik dəyərindən böyük bir TWI-yə sahib olan bütün şəbəkə hüceyrələri quru axını və geri dönüş axını yaşanacağı proqnozlaşdırılan cari səth doymuş sahəni təyin edir.

Müəyyən bir DEM üçün həqiqətən bu mənfi TWI dəyərlərini hansı şərtlərdə tapacağınızı anlaya bilərsiniz. DEM qətnamə üçün g, mümkün olan ən aşağı a dəyər bərabərdir (və ya TWI tətbiqinin kontur uzunluğunu necə təyin etdiyinə görə təxminən bərabərdir) g (yəni g^2 / g). Bu minimum dəyəri deyə bilərik amin. Bununla şəbəkə hüceyrələri a dəyər yalnız topoqrafik silsilə xətləri boyunca baş verəcək heç bir daxil olan hüceyrə olmadığı yerdə baş verir. Maraqlıdır ki, buralar yüksək yerli yamac qradiyentləri ilə qarşılaşma ehtimalı olan yerlərdir. Beləliklə, indi TWI'nin mənfi dəyərləri almağa başladığı nöqtədə yamac həddini hesablaya bilərsiniz: Bthresh = atan (amin). 8.86 m bir DEM çözünürlüğünüz olduğundan, bununla birlikdə hər hansı bir grid hüceyrəsi var a yamacının 83.56 dərəcədən böyük olduğu mənfi TWI dəyərlərini alacaq. Bir çox insanın TWI-nin mənfi ola bilməyəcəyinə inandığını görə bilərsiniz; bu kifayət qədər yüksək bir yamacdır və mənfi TWI alacağınız şərtlər bir qədər nadirdir. O zaman belə maraqlı bir mühitdə çalışdığınız üçün şanslısınız. (DİQQƏT: TWI-nin bu günə qədər tətbiq olunduğu çox geniş mənzərələrdə istifadə edilməsi heç vaxt nəzərdə tutulmamışdır.)

Yuxarıdakı TWI rasterində, mənfi TWI olan hər bir hüceyrə hüceyrəsini vurğuladım. Gördüyünüz kimi, şəbəkədəki bütün bunlar ümumi deyil. Əslində, bu vəziyyətdə, çox dik və kiçik olduqları üçün ərazi cisimləri ilə əlaqəli (bu LiDAR DEM-dəki binalar) şəbəkə hüceyrələri ilə tamamilə məhdudlaşır. a dəyərlər. Bu konkret hal proses baxımından çox real deyil (binanın rütubət indeksi nə deməkdir?), Lakin buna baxmayaraq, fikrimi yaxşı vurğulayır.

Jeffery Evans şərhində mənfi dəyərlərin səbəbinin sıfır yamac dəyərlərinin işlənmə üsulu ilə olması ehtimalını qaldırır. Əslində, bu düz sahələrdə mənfi dəyərlərdən çox çox böyük TWI dəyərləri olacaqdır. Çox aşağı yamaclarda vadinin dibi mövqelərinə meylli olan düz ərazilər, hövzədəki ən nəm və doymuş yerlər arasındadır. İndi açıq şəkildə B sıfır, qaralma (B) eyni zamanda sıfırdır və sıfıra bölmə əldə edirsiniz, bu da əlbəttə sonsuzluqdur. Bu, mükəmməl düz bir yeriniz varsa, həmişə aşağıya doğru su tökə bilmədiyi üçün doymuş olacaq demək kimi bir şeydir. Yenə də unutmayın ki, TWI topoqrafiyaya əsaslanan nisbi bir nəmlik göstəricisi deməkdir, TWI dəyəri 2 olan bir hüceyrənin TWI 1 olan hüceyrədən iki qat daha doymuş olduğunu demirik. Sonsuzluq dəyərləri olan bir şəbəkəyə sahib deyilsiniz və DEM-in keyfiyyətinə və qaynaq məlumatlarına görə sıfır gradiyentlər çox güman ki. Beləliklə, TWI alətlərini inkişaf etdirənlərin sıfır lokal qradiyentli şəbəkə hüceyrələrini idarə etmək üçün xüsusi qaydalara sahib olmaları lazımdır. Ümumiyyətlə istifadəçiyə bu saytın kifayət qədər nəm olduğunu bildirmək üçün sadəcə "həqiqətən yüksək bir dəyər" əvəz edəcəksiniz. Bu yanaşmanın problemi budur ki, mənfi TWI dəyərlərini niyə aldığınızı soruşduğunuz kimi, istifadəçilər də ümumiyyətlə 'niyə TWI dəyərlərini 9999 qədər alıram?' Deyəcəklər. Buna görə alternativ sadəcə bu hüceyrələrə NoData dəyərini qoymaqdır.

Ümid edirəm sualları cavablandırmağa kömək edir.


Videoya baxın: Базовый топографический анализ в SAGA GIS. (Oktyabr 2021).