BAB 5: Hitung Sekali, Berlaku untuk Semua: Kalkulasi Field

#Studi kasus: "Luas di tabel itu hasil ukur lama"

Kolom LUAS_LAMA_HA di tabel petak berasal dari ukuran lama di lapangan. Kepala Seksi curiga angkanya sudah tidak pas. Ia meminta tiga hal: luas yang dihitung ulang dari bentuk petak, umur tanaman, dan kelas umur (muda, sedang, tua). Mengetik 16 angka dengan kalkulator memang bisa. Tetapi bila petaknya 1.600, Anda butuh cara yang menghitung semuanya sekaligus.

#Konsep: kalkulasi field dalam tiga kalimat

Kalkulasi field adalah menyuruh GIS mengisi satu kolom untuk semua baris sekaligus memakai satu rumus, seperti rumus yang diseret ke bawah di lembar hitung. Rumus itu disebut ekspresi, dan bisa memakai isi kolom lain, angka, serta fungsi (perintah siap pakai seperti round). Khusus di GIS, ekspresi juga bisa membaca bentuk objek, misalnya luas dan panjangnya.

Istilah baru bab ini:

• Field Calculator: alat untuk mengisi kolom dengan ekspresi.
• Ekspresi: rumus yang dihitung untuk tiap baris.
• Fungsi: perintah siap pakai, ditulis dengan tanda kurung, misalnya round(...).
• Operator: tanda hitung atau sambung (+, -, /, ||).
• Field virtual: kolom yang dihitung setiap kali dibuka dan tidak ditulis ke berkas.

#Ekspresi yang akan Anda pakai

#Tiga hal yang sering mengejutkan

• Pembagian memberi desimal. 5 / 2 menghasilkan 2.5. Untuk bilangan bulat, pakai 5 // 2 yang menghasilkan 2.
• Kosong merusak hitungan. 1 + NULL menghasilkan kosong (NULL), bukan 1. Bungkus dengan coalesce(kolom, 0) bila perlu.
• Dua cara menghitung luas. area($geometry) menghitung datar di bidang peta (cocok untuk UTM). $area memakai elipsoid yang diatur di proyek. Pada data latihan, selisihnya kurang dari 1 m² per petak (diuji dengan elipsoid WGS 84), jadi hasil pembulatan ke tiga desimal hektare hampir selalu sama. Pilih satu cara dan konsisten.

#Angka dari data latihan

Hasil skrip 5.1 di QGIS 4.0.2:

Contoh baris: P-05 luasnya 1,070 ha, ukuran lama 1,00 ha, selisih +0,070 ha, umur 21 tahun, kelas Tua. P-15 luasnya 1,065 ha, umur 6 tahun, kelas Muda.

Jadi jawaban untuk Kepala Seksi: luas lama memang meleset, paling jauh sekitar 0,07 ha per petak, meski total tetap 16 ha.

round(area($geometry) / 10000, 3)

def tambah(lyr, nama, tipe, rumus, panjang=0, presisi=0):
    return processing.run("native:fieldcalculator", {
        "INPUT": lyr, "FIELD_NAME": nama, "FIELD_TYPE": tipe, "FIELD_LENGTH": panjang,
        "FIELD_PRECISION": presisi, "FORMULA": rumus, "OUTPUT": "memory:"})["OUTPUT"]

h = tambah(petak, "LUAS_HA", 0, "round(area($geometry) / 10000, 3)", 10, 3)

#Cek paham

1. Apa yang terjadi bila ekspresi memuat kolom yang isinya kosong (NULL)?
2. Mengapa kolom UMUR harus dibuat sebelum KELAS_UMUR?
3. Apa beda area($geometry) dan $area?

Jawaban:

1. Hasil hitungan biasanya ikut kosong. Gunakan coalesce(kolom, 0) untuk memberi nilai cadangan.
2. Karena rumus KELAS_UMUR membaca isi UMUR. Kolom yang dibaca harus sudah ada dan terisi.
3. area($geometry) menghitung datar di bidang peta. $area mengikuti elipsoid proyek. Pada data latihan, selisihnya di bawah 1 m² per petak.

#Kesalahan umum

• Salah tanda kutip. "UMUR" adalah nama kolom, 'Muda' adalah teks. Tertukar, ekspresi galat.
• Tipe kolom tidak cocok. Menaruh teks ke kolom angka gagal atau menghasilkan kosong. Pilih tipe sesuai isi.
• Presisi terlalu kecil. Kolom dengan panjang atau presisi pas-pasan memotong angka. Beri ruang lebih.
• Menghitung luas di CRS derajat. Luasnya bermakna ganjil, seperti di bab 5 seri B1. Hitung di CRS meter (UTM).
• Menimpa kolom penting. Gunakan kolom baru lebih dulu, bandingkan, baru ganti bila yakin.

#Ringkasan dan latihan

Ringkasan: Field Calculator mengisi satu kolom untuk semua baris dengan satu ekspresi. Ekspresi memakai kolom, fungsi, dan operator, dan bisa membaca bentuk objek. Pakai kolom baru, periksa pratinjau, lalu simpan.

Latihan: buat kolom KATEGORI dengan ekspresi yang menghasilkan Jati bila jenisnya Jati, dan Non-jati untuk yang lain. Berapa petak berkategori Jati? (Petunjuk: if("JENIS" = 'Jati', 'Jati', 'Non-jati'). Jawaban: 5 petak.)

#Tabel perbandingan: kalkulasi field