Kurva Hipsometrik , HI, dan Indeks R1 Bagian-1

Hipsometrik dihitung dengan DAS. Ada beberapa cara untuk melakukannya. Anda dapat membuat kelas histogram (area dengan kelas-kelas tertentu, 100 kelas) dengan alat Reklasifikasi. Anda dapat membuat kurva hipsometrik menggunakan Jerry Davis ‘Hypsometry Toolbox (area kumulatif menurut nampan, n = 100). Anda dapat melakukannya, sebagian, menggunakan skrip R oleh Veit Hoefler. Anda dapat membuat peta nilai HI menggunakan metode berbasis grid yang disediakan oleh Shantamoy Guha.

Gambaran umum konsep hypsometry disediakan di bawah ini Langkah # 1-3.

1a.) Mengklasifikasi DEM Menjadi 100 kelas elevasi dengan R Script
– Jika Anda menggunakan program R, skrip berikut, di bawah ini, oleh Veit Hoefler mengklasifikasi DEM menjadi kelas elevasi. Jika bukan pengguna R, lewati ke langkah berikutnya (“Mulailah dengan proyeksi DEM …”).

WESSO<-function(DGM,Bin){
min<-trunc(minValue(DGM))
max<-ceiling(maxValue(DGM))
DIFF<-max-min
class<-DIFF/Bin
break_v<-matrix(seq(min+class,max,class))
low<-cbind(seq(min,max-class,class))
numb<-seq(1,Bin,1)
rclmat<-cbind(low,break_v,numb)
rc<- reclassify(DGM, rclmat,include.lowest=T)
list<-list(a=rclmat,b=rc)
raster<-list$b
rclmat<-list$a
count<-freq(raster,useNA='no')
count_2<-hist(raster,breaks=Bin)
count_2<-c(count_2$counts)
count_2[1]
DIFFA<-as.numeric(count_2[1]-count[1,2])
laenge<-as.numeric(length(count[,2]))
count_3<-count[3:laenge,2]
DIFFA_2<-as.numeric(count[1,2])
as<-rbind(DIFFA_2,DIFFA)
nn<-matrix(rbind(count_2[2:length(count_2)]))
end<-rbind(as,nn)
end<-matrix(end)
info<-cbind(rclmat,end)
return(list(a=raster,b=info))}</code>

1b.) Mengklasifikasi DEM Menjadi kelas elevasi 100 dengan ArcGIS & Spatial Analyst
– Mulailah dengan DEM yang diproyeksikan terpotong ke batas DAS studi Anda

  • Konversi DEM Anda dari Float -> format Integer menggunakan Spatial Analyst> Math> Int tool.
  • Gunakan alat Reklasifikasi untuk membuat kelas elevasi Anda (n = 100) * dari DEM integer baru …

Spatial Analyst > Reclass > Reclassify tool
Input raster = watershed DEM
Reclass field = VALUE
Method = Equal interval
Number of Classes = 100

  • Catatan: Jika Anda menentukan kelas 100 di alat, tetapi hanya 10 kelas yang dibuat, maka Anda perlu menambah jumlah kelas maksimum yang diizinkan menggunakan “ArcMap Advanced Settings Utility.exe”. Mudah digunakan. Jumlah kelas maksimum adalah 200. Saya menemukan utilitas di sini: C: \ Program Files (x86) \ ArcGIS \ Desktop10.1 \ Utilities. Ubah pengaturan di bawah tab “Raster”. Anda harus memulai ulang ArcGIS setelah membuat perubahan. Periksa tabel atribut setelah memulai ulang dan menjalankan alat; itu harus berisi 100 kelas.
  • Tabel atribut tersedia jika raster dalam format integer (lihat tool Int).
  • Buka tabel. Nilai VALUE berisi nomor kelas ID. Kolom COUNT berisi jumlah piksel untuk setiap kelas.

  • Opsional: Untuk membuat grafik data dalam Arc: Table Options > Create Graph.

  • Untuk mengekspor data ke Excel: Table Options > Export to dBase (.dbf).

  • Buka .dbf yang diekspor di Excel. Pastikan untuk membuka Excel terlebih dahulu, lalu temukan file (tunjukkan Semua File *. *). Pastikan TAMPILKAN SEMUA JENIS FILE (jangan buka file .xml).

  • 2.) Bagan Histogram terkelas di Excel

    Contoh Dataset dengan Instruksi & Rumus:
    Contoh Binaan Hypsometry Spreadsheet_Spring2015.xlsx

    • Buat Bagan Kombo dengan tipe Kolom Vertikal dan Scatter. Di Excel 8, Recommended Charts > All Charts > Combo Chart…

    Seri 1 = Bilah Kelas Area (tipe kolom, sumbu y primer)
    Seri 2 = Kurva Area Kumulatif (tipe pencar, sumbu y sekunder)
    Sumbu X untuk kedua seri = Nomor Kelas

     

     

     

    3.) Bagan Kurva Hipsometrik dengan Hypsometry Toolbox (Jerry Davis)
    – Unduh di sini: Hypsometric Tools Script oleh Jerry Davis
    – Unduh folder toolbox / skrip.
    – Buka zip folder ke Desktop Anda.
    – Buka ArcToolbox> Klik kanan pada kata “ArcToolbox” di atas bingkai> Add Toolbox > navigasi ke desktop> pilih folder yang sudah unzipped> Subfolder Hypsometry> Hypsometry.tbx. Klik OK. Toolbox akan diinstal.
    – Klik menu “Hypsometry Tools: Generic Tools”> buka skrip Hypsometry:
    Input = DAS DAS
    Output = nama raster yang akan dibuat
    – Jalankan.
    – Saat selesai (sekitar 10 detik), buka table Options > Export table as .dbf (dBase Table).
    – Jangan menambahkan tabel ke peta.
    – Buka .dbf baru di Excel (tunjukkan Semua File *. *) Dan bagan seri kurva hipsometrik (tipe pencar) di atas histogram biner (jenis kolom).
    Kolom VALUE = Nomor kelas (sumbu X).
    Kolom CUMFREQ = Area Kumulatif (Sumbu Y).
    – Format sumbu-Y Primer dan Sekunder dengan hati-hati.
    – Lihat gambar contoh.

    4.) Metode Berbasis Grid untuk Memetakan Nilai Integral Hipsometrik (Guha Shantamoy)
    Dengan metode ini, Anda mulai dengan DEM daerah aliran sungai dan diakhiri dengan peta nilai HI yang didistribusikan secara spasial.
    Dokumen Word: Gridded Hypsometric Integral menggunakan Arc Map 10_Guha.doxc

     

     

     

    http://gis4geomorphology.com/hypsometric-index-integral/

    Tinggalkan Pesan