Lebih

Output aneh dari indeks basah topografi RSAGA?

Output aneh dari indeks basah topografi RSAGA?


Saya mencoba membuat raster indeks kebasahan topografi dengan paket RSAGA di R. Output yang dapat saya hasilkan tidak masuk akal bagi saya, karena memiliki banyak sel NA dan saya bertanya-tanya apakah saya melakukan sesuatu yang salah. Saya telah mengumpulkan contoh dummy sederhana untuk mengilustrasikan apa yang saya maksud.

Info sesi: Saya menjalankan R 3.2.3 di Ubuntu 15.10. Saya telah menginstal SAGA v2.2.0, serta RSAGA v0.94-5.

library(RSAGA) library(raster) library(maptools) #get DEM untuk Swiss alt <- getData('alt', country="CHE") # tulis ke format SAGA tempInfile <- 'saga_temp_twi_in' writeRaster(alt, filename=tempInfile , format="SAGA") outfile <- 'saga_temp_twi_out.sgrd' panggil <- rsaga.wetness.index(in.dem='saga_temp_twi_in.sgrd', out.wetness.index=outfile, env=rsaga.env(modules= '/usr/lib/x86_64-linux-gnu/saga/')) #membaca file yang dihasilkan kembali di res <- raster('saga_temp_twi_out.sdat') par(mfrow=c(1,2)) plot(alt) plot (res)


Saya pikir saya melihat masalahnya.

Saya dapat mereproduksi ini, baik menggunakan skrip Anda, dan SAGA itu sendiri. Keluaran dari SAGA cocok dengan keluaran dari skrip Anda. Jadi bukan RSAGA.

Tapi lihatlah lereng output dari SAGA… area yang luas hampir vertikal, dengan kemiringan sekitar 1,57 radian (hampir 90 derajat) dengan beberapa area datar. Anda dapat mengkonfirmasi ini dengan melihat histogram kemiringan di SAGA…

Anda hanya melihat air menumpuk di daerah datar… di mana kemiringannya sekitar 0. Kemiringannya bertindak seperti topeng.

Saya pikir masalahnya adalah file tidak memiliki georeferensi, jadi SAGA berpikir bahwa Swiss adalah beberapa meter (ukuran raster Anda, yang tentu saja dalam derajat). Ini mengasumsikan jarak vertikal dalam satuan yang sama dengan horizontal… tidak menyadari bahwa satuannya adalah derajat. Tapi Swiss masih tingginya beberapa ribu meter! Hasil akhirnya adalah berlebihan vertikal besar-besaran.

Sebagai pemeriksaan kewarasan cepat, saya berlari Isi Wastafel, kemudian digunakan Kalkulator kotak untuk membagi ketinggian dengan 110000 (kira-kira konversi meter ke derajat pada garis lintang 45N) Menjalankan SAGA Wetness Index pada ini, memberikan hasil yang jauh lebih masuk akal :)

Tidak yakin bagaimana melakukan ini di RSAGA. Taruhan terbaik Anda adalah menggunakan raster dengan CRS yang ditetapkan padanya, atau menetapkan CRS di SAGA itu sendiri. Trik pembagian saya tidak akan memberikan hasil seakurat yang Anda dapatkan dengan raster yang diproyeksikan dengan benar.


Tonton videonya: Jannes Münchow: R-GIS bridges: the examples of RQGIS and RSAGA