Mengekstraksi jaringan sungai dan batas Daerah Aliran Sungai (DAS) dari Digital Elevation Model (DEM) adalah tugas klasik dalam hidrologi. Namun, melakukan langkah Fill, Flow Direction, Flow Accumulation, hingga Watershed Delineation secara manual melalui antarmuka (GUI) akan sangat menyita waktu jika Anda harus memproses puluhan Sub-DAS.

Meninggalkan Klik Manual dengan Python

Sebagai Spatial Data Engineer, prinsip utama yang harus dipegang adalah: "Jika Anda mengklik alat yang sama lebih dari tiga kali, tulislah skrip untuk itu."

# Contoh skrip Python (ArcPy) untuk Delineasi DAS import arcpy from arcpy.sa import * arcpy.env.workspace = "C:/Hydrology/Data" outFill = Fill("dem_srtm.tif") outFlowDir = FlowDirection(outFill, "NORMAL") outFlowAcc = FlowAccumulation(outFlowDir) # Membuat batas sungai berdasarkan threshold stream_network = Con(outFlowAcc > 1000, 1, "") stream_network.save("streams.tif")

Mengapa Ini Penting?

Otomatisasi tidak hanya menghemat waktu, tetapi juga menjamin reproduktibilitas. Jika bulan depan tersedia data DEM resolusi lebih tinggi (misalnya LiDAR), kita cukup mengganti input pada skrip di atas dan menjalankannya kembali, alih-alih mengulang pekerjaan berjam-jam dari awal.