* Staf Pengajar Sekolah Tinggi Pertanahan Nasio- nal. Email: vinysari@yahoo.com

## UJI AKURASI  KLASIFIKASI PENGGUNAAN LAHAN DENGAN MENGGUNAKAN METODE DEFUZZIFIKASI MAXIMUM LIKELIHOOD
BERBASIS CITRA ALOS AVNIR-2
Harvini Wulansari*
Abstract: Abstract: Abstract: Abstract: Abstract: Land use information plays an important role in spatial planning and monitoring development in order to optimize landuse and to minimize land conflict. Remote sensing technology can be used to derive land cover information for land use information. The aim of this research was to study the accuracy level and the efficiency of fuzzy logic for land use classifica- tion. Defuzzification method  was implemented using  maximum  likelihood and its landuse map from classification result, using spectral data from ALOS AVNIR-2. This research used fuzzy logic approach with defuzzification method using maxi- mum likelihood. Sample plots used for modelling were chosen using stratified random sampling method, and accuracy as- sessment was obtained using stratified random sampling method also. In overall, the research process worked well, although the standpoint of accuracy and thoroughness resulted in overall kappa index was less good or less feasible; however, the results were acceptable. For the 14 landuse classes, it resulted 57% of overall accuracy and kappa index of 0.53. Time execution using maximum likelihood algorithm was approximately 10 seconds.
KKKKKey W ey W ey W ey W ey Wor or or or ords: ds: ds: ds: ds: landuse classification, ALOS AVNIR-2, fuzzy logic, maximum likelihood.
Intisari: Intisari: Intisari: Intisari: Intisari:  Informasi data penggunaan lahan sangat berperan penting diantaranya untuk melakukan perencanaan maupun pengawasan perkembangan suatu wilayah, sehingga penggunaan lahannya dapat dimanfaatkan secara optimal dan tetap menjaga kelestariannya, serta meminimalkan terjadinya konflik terhadap lahan. Teknologi penginderaan jauh merupakan salah satu cara yang dapat dimanfaatkan untuk menurunkan data penutup lahan sebelum akhirnya diterjemahkan menjadi informasi penggunaan lahan. Tujuan dari penelitian ini yaitu untuk mengetahui seberapa akurat dan efisien metode pendekatan berbasis logika samar ( fuzzy logic) dengan defuzzifikasi menggunakan algoritma maximum likelihood serta peta hasil klasifikasi penggunaan lahannya, yang melibatkan input data spektral dari Citra ALOS AVNIR-2. Metode yang digunakan dalam penelitian ini yaitu metode pendekatan berbasis logika samar (fuzzy logic) dengan defuzzifikasi menggunakan algoritma maximum likelihood. Penentuan sampel untuk training area dan uji akurasi menggunakan metode plot area sedangkan pengambilan sampel menggunakan metode stratified random sampling. Secara keseluruhan  proses penelitian berhasil dengan baik, walaupun dari sudut pandang  ketelitian  menghasilkan overall accuracy dan indeks kappa yang kurang baik atau kurang layak, namun demikian hasilnya masih dapat diterima.Untuk 14 kelas penggunaan lahan ( overall accuracy 57%, nilai indeks kappa 0,53). Hasil pencatatan waktu untuk proses eksekusi menggunakan algoritma maximum likelihood sekitar 10 detik.
Kata kunci Kata kunci
Kata kunci Kata kunci
Kata kunci: klasifikasi penggunaan lahan, Citra ALOS AVNIR-2, logika samar, maximum likelihood.

## A. Pendahuluan
Informasi data penggunaan lahan berperan penting diantaranya untuk melakukan perencana- an maupun pengawasan perkembangan suatu wilayah, sehingga penggunaan lahannya dapat di- manfaatkan secara optimal dengan tetap menjaga
kelestariannya, serta meminimalkan terjadinya konflik terhadap lahan. Menurut Lillesand et al. (2008) informasi penggunaan lahan tidak selalu dapat ditafsir langsung dari penutup lahannya oleh karena itu diperlukan informasi pelengkap untuk menentukan penggunaan lahan. Untuk menurun- kan informasi penggunaan lahan dari citra pengin- deraan jauh diperlukan informasi pendukung di antaranya seperti data bentuk lahannya atau bisa
Harvini Wulansari: Uji Akurasi Klasif ikasi Penggunaan Lahan ... : 98-110
dengan local knowledge mengenai wilayah tersebut. Saat ini telah berkembang metode klasifikasi multispektral yang bisa dilakukan melalui berbagai pendekatan diantaranya menurut Jensen (2005) hard classification atau soft classifaction tergantung output yang dikehendaki, sedangkan berdasarkan distribusi datanya dapat  menggunakan algoritma pendekatan parametrik misal maximum likelihood atau dengan pendekatan non parametrik misal jaringan syaraf tiruan ( artif icial neural network ). Salah satu contoh algoritma yang menggunakan pendekatan parametrik yaitu algoritma maximum likelihood. Pendekatan parametrik mengasumsikan bahwa distribusi statistik kelas digambarkan dengan data yang terdistribusi normal atau bayesian . Sedangkan pendekatan non parametrik tidak mengharuskan distribusi statistik kelas digambarkan secara normal, salah satu contoh pendekatan non parametrik yaitu jaringan syaraf tiruan ( artif icial neural network ).
Lillesand et al. (2008) menyatakan bahwa data penginderaan jauh merupakan hasil interaksi antara tenaga elektromagnetik dengan objek yang diindera yang direkam oleh sensor, dimana setiap objek mempunyai karakteristik tertentu dalam berinteraksi dengan setiap spektrum elektromag- netik. Danoedoro (2012)  juga menyebutkan bahwa pada prinsipnya setiap benda memiliki struktur partikel yang berbeda, dimana perbedaan ini akan mempengaruhi pola respon elektromagnetiknya, yang dapat dijadikan landasan pembedaan objek.
Banyak penelitian mulai dikembangkan terkait dengan karakteristik pola respon spektral, di antaranya untuk mengektraksi data penutup lahan dari citra satelit. Pendekatan yang mulai berkem- bang diantaranya adalah klasif ikasi berbasis logika samar ( fuzzy logic ). Metode fuzzy sangat mem- bantu saat proses pengambilan sampel, di mana tidak terjadi pemaksaan keanggotaan piksel untuk masuk ke dalam satu kelas saja ( soft classif ication ). Sampel yang diambil juga tidak harus piksel murni, dapat diambil pada mixed pixel. Seperti diketahui
bahwa dari sudut pandang penginderaan jauh fenomena geografis pada level skala menengah dan rendah bersifat fuzzy artinya tidak ada batas yang jelas (dalam arti tajam) antara fenomena-fenome- na geografis tersebut. Masalah mixed pixel terutama terjadi pada citra penginderaan jauh resolusi spasial menengah dan resolusi spasial rendah, dimana di dalam satu piksel citra dapat terdiri dari dua atau lebih jenis objek.
Metode pendekatan berbasis logika samar ( fuzzy logic ) secara operasional sistematik masih jarang digunakan di instansi-instansi yang bertugas membuat peta tematik seperti misalnya Badan Pertanahan Nasional/BPN RI. Selain itu penggunaan citra penginderaan jauh resolusi spasial menengah di instansi BPN RI sebagai masukan data peng- gunaan lahan juga jarang digunakan, data yang tersedia di instansi BPN RI berupa foto udara skala besar, citra dengan resolusi spasial tinggi seperti Citra Quickbird, sedangkan untuk citra dengan resolusi spasial menengah seperti Citra ALOS 1 AVNIR-2  jarang digunakan, oleh karena itu perlu dikaji tingkat akurasi dari citra tersebut. Dari per- masalahan tersebut maka peneliti melakukan sebuah penelitian untuk mengetahui seberapa akurat dan ef isien pemetaan penggunaan lahan yang dihasilkan dengan menggunakan metode pendekatan berbasis logika samar ( fuzzy logic ) dengan defuzzif ikasi menggunakan algoritma maximum likelihood . Penelitian ini dilakukan di sebagian wilayah Daerah Istimewa Yogyakarta dan sebagian wilayah Kabupaten Klaten, dengan luasan
1 Satelit ALOS ( Advanced Land Observing Satellite ) merupakan satelit sumber daya milik Jepang, yang ter- diri dari tiga sensor yaitu Panchromatic Remote Sensing Instrument for Stereo Mapping (PRISM) dengan resolusi 2,5 meter, Advanced Visible and Near Infrared Radiom- eter type-2 (AVNIR-2) resolusi 10 meter dan Phased Array type L-band Synthetic Aperture Radar (PALSAR) resolusi 10 meter dan 100 meter. Walaupun saat ini su- dah tidak beroperasi namun data yang dihasilkan masih banyak dimanfaatkan untuk berbagai macam kebu- tuhan informasi dan ekstraksi data.

## 100 Bhumi Vol. 3 No. 1 Mei 2017
1506 x 1505 piksel atau sekitar 15 x 15 km 2 (Citra ALOS AVNIR-2 resolusi spasial sekitar 10 m). Citra ALOS AVNIR-2  yang digunakan yaitu citra yang direkam tanggal 20 Juni tahun 2009 level 1B2G, dapat dilihat pada gambar 1.
Gambar 1. Lokasi Penelitian

## B.  Klasif ikasi Multispektral
Danoedoro (2012) menyebutkan bahwa klasi- fikasi multispektral merupakan suatu metode yang dirancang untuk menurunkan informasi tematik dengan cara mengelompokkan fenomena berda- sarkan kriteria tertentu, yang pada dasarnya diasumsikan bahwa tiap objek dapat dibedakan dengan objek yang lain berdasarkan nilai spektral- nya dimana tiap-tiap objek tersebut cenderung memberikan pola respon spektral yang spesifik.
Algoritma klasifikasi multispektral berdasarkan pada asumsi distribusi datanya menurut Scho- wengerdts (1983) dapat dibedakan menjadi para- metrik dan non parametrik. Algoritma parametrik mengasumsikan bahwa distribusi statistik kelas digambarkan dengan distribusi normal dan mengestimasikan parameter yang didistribusikan dengan vektor rerata dan matriks kovarian, sedang- kan algoritma non parametrik tidak menga- sumsikan distribusi kelas.
Algoritma maximum likelihood merupakan salah satu contoh pendekatan parametrik. Danoedoro (2012) menyebutkan bahwa Algoritma maximum likelihood merupakan algoritma yang secara
statistik mapan, dengan asumsi bahwa objek homogen selalu menampilkan histogram yang terdistribusi normal atau bayesian .
Menurut Shresta (1991 dalam Danoedoro 2012) algoritma maximum likelihood dikelaskan sebagai objek tertentu berdasarkan bentuk ukuran dan orientasi sampel pada feature space yang berbentuk elipsoida, dapat dilihat gambaran algoritma maxi- mum likelihood pada gambar 2. Informasi nilai statistik yang dibutuhkan pada klasif ikasi meng- gunakan algoritma maximum likelihood yaitu nilai rerata ( mean ), simpangan baku serta variansi dan kovariansinya.
Gambar 2. Fungsi kerapatan yang digambarkan dengan algoritma maximum likelihood modifikasi dari Lillesand et al, (2008, dalam Danoedoro,
2012)

## C. Klasifikasi dengan Logika Fuzzy
Menurut Kusumadewi (2010), logika fuzzy merupakan salah satu komponen pembentuk soft computing, dimana peranan derajat keanggotaan atau membership function sebagai penentu kebera- daan elemen dalam suatu himpunan sangat penting. Himpunan fuzzy nilai keanggotaan terletak pada interval 0 sampai dengan 1. Novamizanti et.al. (2010) dalam tulisannya menjelaskan bahwa sebuah sistem fuzzy memiliki struktur proses seba- gai berikut: 1). Fuzzif ication (fuzzif ikasi), yaitu proses memetakan crisp input ke dalam himpunan fuzzy . Hasil dari proses ini berupa fuzzy input dalam bentuk rule fuzzy 2). Evaluation  Rule , yaitu proses melakukan penalaran terhadap fuzzy input yang dihasilkan oleh proses fuzzif ication berdasarkan
Harvini Wulansari: Uji Akurasi Klasif ikasi Penggunaan Lahan ... : 98-110
aturan fuzzy yang telah dibuat. Proses ini mengha- silkan fuzzy output 3). Defuzzif ication (defuzzi- f ikasi), yaitu proses mengubah fuzzy output men- jadi crisp value.
Jensen (2005) menggambarkan bagaimana logika fuzzy diterapkan untuk proses klasifikasi nilai piksel terhadap jenis penutup/penggunaan lahan dapat dilihat pada gambar 3. Dari gambar tersebut dapat dilihat perbandingan antara klasi- f ikasi tegas dan klasifikasi berbasis logika fuzzy.
Gambar 3. Perbandingan antara kalsifikasi tegas dan klasifikasi berbasis logika fuzzy (Jensen, 2005)

## D. Tahapan Lapangan
Tahapan penelitian yang dilakukan meliputi Tahap Pra Lapangan, Tahap Kegiatan Lapangan, Tahap Pemrosesan dan Tahap Pasca Lapangan (Uji Akurasi). Tahapan-tahapan tersebut dapat dilihat pada gambar 4.

## Gambar 4. Diagram alir penelitian

## 1. Tahap Pra Lapangan

## a. Koreksi Geometrik
Sebelum citra digital diolah terlebih dahulu dilakukan koreksi citra. Tahapan ini dilakukan karena saat sensor melakukan perekaman data di permukaan bumi terdapat beberapa faktor yang mempengaruhi kualitas citra baik dari wahananya sendiri atau di luar wahana. Koreksi geometrik dila- kukan untuk melakukan perbaikan distorsi geo- metrik citra terhadap kondisi sebenarnya di lapangan. Koreksi geometrik yang dilakukan yaitu dengan transformasi berdasarkan ground control point (GCP) dengan melakukan rektifikasi citra ke peta, prinsip rektifikasi citra ke peta bahwa peta

## 102
Bhumi Vol. 3 No. 1 Mei 2017
yang digunakan sebagai acuan memiliki sistem proyeksi dan koordinat yang lebih dapat dipercaya atau lebih dianggap benar daripada citra. Peta yang digunakan acuan untuk koreksi geometrik adalah peta RBI skala 1: 25.000. Proses koreksi geometrik dilakukan pada software ENVI dengan cara regis- trasi image to map . Peta RBI skala 1: 25.000 yang digunakan acuan sudah dilakukan proses geore- ference terlebih dahulu.
Interpolasi spasial yang digunakan mengguna- kan fungsi persamaan polinomial orde dua sebanyak 10 titik GCP yang tersebar merata di wilayah penelitian, dengan menggunakan proyeksi Universal Transverse Mercator (UTM) datum WGS
– 84 zone 49 M. Penentuan 10 titik GCP tersebut sudah mencukupi karena syarat untuk polinomial orde dua adalah sebanyak 6 titik dan juga setelah di-overlay dengan peta RBI skala 1: 25.000 sudah sama antara perpotongan jalan dan percabangan sungainya. Titik-titik kontrol diambil secara merata di wilayah penelitian yang dapat dilihat pada gambar 5. Titik-titik kontrol diambil pada perpo- tongan jalan, jembatan, percabangan aliran sungai dan bangunan yang bersifat permanen/tetap.
Koreksi geometrik menghasilkan nilai Root Mean Square Error (RMSE) 0,45. Hasil tersebut sudah memenuhi syarat untuk citra tunggal mini- mal RMSE yang diperbolehkan kurang dari 1. Pada saat proses koreksi geometrik dilakukan pula proses resample nilai piksel citra dengan metode nearest neighbour yaitu menempatkan nilai piksel tetangga terdekat untuk mengisi piksel yang digeser agar nilai piksel hasil koreksi tidak jauh berbeda dengan nilai piksel citra sebelum dilakukan proses koreksi geometrik.

## Gambar 5. Lokasi Titik-titik GCP

## a. Penentuan Skema Klasifikasi Penggunaan Lahan
Citra ALOS AVNIR-2 yang memiliki resolusi spasial 10 meter menurut Doyle (1984 dalam Sutanto, 2010) tingkat akurasi interpretasi atau akurasi klasifikasi sesuai dengan rumus a rule of thumb yaitu sama dengan tingkat akurasi pada skala peta 1: 100.000. Oleh karena itu, dalam penelitian ini dipilih skema klasifikasi penggunaan lahan untuk skala 1: 50.000 yang dibuat oleh BPN RI yang dianggap masih mendekati skala 1:100.000. Dalam penentuan peng- gunaan lahannya dilakukan modifikasi yang berbasis pada penutup lahannya  dilihat pada pada tabel 1.
Tabel 1. Skema klasifikasi penggunaan lahan yang digunakan dalam penelitian (modif ikasi atas klasifikasi penggunaan lahan BPN RI)
Singkatan
Harvini Wulansari: Uji Akurasi Klasif ikasi Penggunaan Lahan ... : 98-110

## 1. Tahap Kegiatan Lapangan
Kegiatan lapangan meliputi kerja lapangan ditujukan untuk memperoleh data lapangan yang digunakan sebagai penentuan daerah contoh atau training area sekaligus mencari informasi apakah terjadi adanya perubahan penggunaan lahan antara citra yang digunakan dengan kondisi saat ini  serta cek lapangan untuk sampel uji akurasi.
Citra yang digunakan dalam penelitian ini yaitu Citra ALOS AVNIR-2 yang direkam pada tahun 2009, sedangkan penelitian dilakukan pada tahun 2013 sehingga dimungkinkan telah terjadi peru- bahan penggunaan lahan di lapangan, sehingga perlu dilakukan  kerja lapangan untuk melakukan wawancara dengan penduduk sekitar.
Pendekatan yang digunakan dalam survei dan pemetaan penggunaan lahan yaitu photo key ap- proach yang merupakan pendekatan yang bersifat fotomorfik dimana kenampakan pada foto men- jadi kunci pengenalan objek, dengan menggunakan citra resolusi spasial tinggi berupa Citra Quickbird . Hal ini dilakukan untuk membantu mencocokkan sampel, yang sebelumnya sudah ditentukan sebe- lum berangkat ke lapangan dengan kenyataan yang sebenarnya dilapangan sebagai contoh sebelum ke lapangan hasil interpretasi bahwa objek tersebut adalah sawah irigasi setelah di lapangan apakah benar objek yang dimaksud memang sawah irigasi atau berupa tegalan/ladang, seperti yang dapat dilihat pada gambar 6.
(a). Citra ALOS AVNIR-2     (b). Citra Quickbird Gambar 6. Pendekatan photo key approach yang menunjukkan obyek sawah irigasi pada (a). Citra ALOS AVNIR-2 dan (b). Citra Quickbird
Penentuan titik-titik sampel yang digunakan untuk daerah contoh menggunakan pendekatan stratified random sampling , titik sampel ditentukan secara acak pada setiap strata dan diusahakan ter- distribusi secara merata di seluruh daerah pene- litian. Sebagaimana dikemukakan oleh Lillesand et.al. (2008) untuk dapat menghitung vektor mean dan matriks kovarian  (seperti metode maximum likelihood ) dibutuhkan minimum 10 n to 100 n piksel dimana n merupakan jumlah band, yang berarti jika dalam penelitian ini menggunakan Citra ALOS yang memiliki 4 band maka minimum 40-400 piksel diambil untuk tiap kategori dan diambil pada beberapa lokasi yang berbeda. Dalam penelitian ini plot area sampel daerah contoh dan sampel uji minimal yaitu 100 piksel, jika Citra ALOS AVNIR dengan resolusi spasial 10 m maka berkisar 10.000 m² atau 1 hektar untuk satu kelas penggunaan la- han, yang diambil pada tempat yang berbeda.
Sesuai dengan acuan skema penggunaan lahan yang digunakan dan dengan melihat keadaan di lapangan, maka dalam penelitian ini ada 14 kelas penggunaan lahan yang digunakan yaitu tegalan/ ladang, semak, permukiman, tanah jasa, tempat olahraga, sawah irigasi, kebun campuran, peter- nakan, industri, pemanfaatan berbasis air terge- nang, pemanfaatan berbasis sungai, sawah non irigasi, hutan sejenis dan penggunaan lain.
Penggunaan lahan pemanfaatan berbasis air tergenang yang dimaksudkan adalah seperti kolam ikan atau empang, sedangkan penggunaan tanah jasa seperti jalan aspal dan bandar udara (bandara). Penggunaan lahan yang dimasukkan dalam kelas penggunaan lain yaitu bangunan candi dan lahan kosong.

## 2. Tahap Pemrosesan

## a. Hasil Uji Separabilitas
Uji separabilitas objek dihitung dengan menggunakan algoritma transformed divergence (TD). Uji separabilitas dilakukan untuk mengeta- hui tingkat keterpisahan nilai spektral antar objek.

## 104 Bhumi Vol. 3 No. 1 Mei 2017
Penelitian ini menggunakan software ENVI untuk menghitung tingkat keterpisahan objeknya, dimana untuk nilai maksimumnya adalah 2 atau setara dengan nilai 2000 apabila menggunakan software Idrisi Selva, yang berarti bahwa kelas yang diuji benar-benar terpisah secara baik atau maksimal sehingga dapat dibedakan antara kelas satu dengan yang lain secara baik. Nilai separabilitas dibawah 1,7 dianggap kurang baik yang berarti bahwa ada mixed piksel/piksel campuran didalamnya. Dalam penelitian ini ditetapkan ambang nilai minimum 1, sebagai batas bahwa kelas tersebut masih bisa dibedakan, apabila kurang dari 1 maka harus dija- dikan satu kelas. Indeks separabilitas 14 kelas peng- gunaan lahan dapat dilihat pada tabel 2.

## Tabel 2. Indeks separabilitas untuk 14 kelas penggunaan lahan
Sumber: Pengolahan penelitian, 2013 Keterangan: = adanya mixed piksel/piksel campuran
Dari tabel 2 dapat dilihat contoh bahwa kelas penggunaan lahan tegalan/ladang (tg) memiliki tingkat keterpisahan yang kurang baik terhadap kelas penggunaan lahan permukiman (pm), peternakan (pt) dan sawah non irigasi (sn), kemudian kelas penggunaan lahan semak yang juga memiliki keterpisahan yang kurang baik terhadap sawah irigasi (si), kebun campuran (kc) dan hutan sejenis.

## b. Hasil Penentuan Himpunan

## Keanggotaan Kelas Penggunaan Lahan
Himpunan nilai derajat keanggotaan piksel untuk masuk dalam satu kelas penggunaan lahan
tertentu ditentukan dari parameter statistik hasil perhitungan statistik sampel. Nilai mean (rerata) dan standar deviasi merupakan parameter statistik utama dalam penentuan himpunan keanggotaan piksel, dimana nilai mean (rerata) menunjukkan nilai tengah kurva keanggotaan piksel terhadap suatu kelas penggunaan lahan, sedangkan standar deviasi sebagai acuan lebar kurva.
Hasil perhitungan statistik untuk 14 kelas peng- gunaan lahan dapat dilihat pada tabel 3 sampai dengan tabel 6.
Tabel 3. Hasil perhitungan statistik band 1 untuk 14 kelas penggunaan lahan

## Sumber : Pengolahan penelitian, 2013
Tabel 4. Hasil perhitungan statistik band 2 untuk 14 kelas penggunaan lahan
Sumber : Pengolahan penelitian, 2013
Kelas Tg sm pm tj to si kc pt ti at sg sn hs pl Tg 0 Sm 1,7 0 Pm 1,2 1,9 0 Tj 1,9 1,9 1,4 0 To 1,9 1,8 1,9 1,9 0 Si 1,7 1,6 1,8 1,9 1,9 0 Kc 1,9 1,4 1,9 1,9 1,9 1,4 0 Pt 1,5 1,7 1,1 1,7 1,6 1,6 1,8 0 Ti 1,9 1,9 1,5 1,4 1,9 1,9 1,9 1,5 0 At 1,8 1,9 1,7 1,7 1,9 1,6 1,9 1,8 1,9 0 Sg 1,9 1,8 1,7 1,8 1,9 1,3 1,6 1,6 1,9 1,0 0 Sn 1,0 1,7 1,4 1,9 1,6 1,6 1,9 1,4 1,9 1,8 1,8 0 Hs 1,8 1,0 1,9 1,9 1,9 1,9 1,3 1,7 1,9 1,9 1,9 1,9 0 Pl 1,7 1,9 1,2 1,0 1,9 1,9 1,9 1,7 1,4 1,7 1,8 1,9 1,9 0 No Kelas LU Band 2
Minimum Maximum Mean Stand. Deviation b c 1 Tg 64 95 72,434 6,812 65,622 79,246 2 Sm 48 76 61,358 5,626 55,732 66,984 3 Pm 52 111 76,78 10,319 66,461 87,099 4 Tj 64 176 77,517 17,862 59,655 95,379 5 To 69 107 77,747 6,340 71,407 84,087 6 Si 57 86 63,266 4,985 58,281 68,251 7 Kc 49 63 55,528 2,513 53,015 58,041 8 Pt 58 205 88,299 32,367 55,932 120,666 9 Ti 59 253 133,92 52,459 81,462 186,38 10 At 50 77 59,472 4,809 54,663 64,281 11 Sg 49 77 63,438 5,122 58,316 68,56 12 Sn 62 99 77,898 6,856 71,042 84,754 13 Hs 46 68 52,594 3,275 49,319 55,869 14 Pl 50 102 75,596 8,806 66,79 84,402
Tabel 5. Hasil perhitungan statistik band 3 untuk 14 kelas penggunaan lahan
Sumber : Pengolahan penelitian, 2013
Tabel 6. Hasil perhitungan statistik band 4 untuk 14 kelas penggunaan lahan

## Sumber : Pengolahan penelitian, 2013
Cara penentuan derajat keanggotaan karena bersifat keahlian penganalisis maka sebetulnya harus dilakukan oleh orang yang benar-benar ahli bukan pemula.  Penelitian ini dalam pertimbangan penentuan himpunan dan derajat keanggotaan selain pengetahuan lokal ( local knowledge ) terhadap wilayah penelitian juga dibantu dengan diagram pencar atau feature space , membuat kurva kosinus secara manual dengan menggunakan software microsoft excel dan juga dengan bantuan citra fuzzy kelas penggunaan lahannya.
Contoh graf ik dapat dilihat pada gambar 7, sedangkan contoh citra fuzzy dapat dilihat pada gambar 8. Pada gambar 8, contoh citra kelas fuzzy menggambarkan bahwa semakin mendekati nilai 1 (warna merah) makin masuk kedalam kelas yang bersangkutan dalam hal ini yaitu hutan sejenis, seba- liknya semakin mendekati nilai 0 (biru tua) maka semakin tidak masuk kelas hutan sejenis.
Setelah penentuan himpunan dan derajat ke- anggotaan dibuat, maka kemudian menyusun matrik partisi fuzzy. Matrik partisi fuzzy ditentukan atas dasar nilai derajat keanggotaan setiap sampel. Matrik partizi fuzzy berupa data tabular yang terdiri dari baris dan kolom, dimana tiap baris berisikan informasi ID kelas penggunaan lahan pada train- ing site yang telah dibuat sementara kolom berisi- kan kelas penggunaan lahannya. Tabel ini diisi sesuai dengan derajat keanggotaan yang telah ditentukan, dapat dilihat pada tabel 7.
No Kelas LU
Band 3 Minimum Maximum Mean Stand. Deviation b c 1 tg 46 101 62,307 11,658 50,649 73,965 2 sm 30 56 40,140 5,264 34,876 45,404 3 pm 38 127 72,386 15,707 56,679 88,093 4 tj 50 165 65,915 19,334 46,581 85,249 5 to 41 110 59,590 12,579 47,011 72,169 6 si 36 69 43,24 6,567 36,673 49,807 7 kc 30 44 34,853 2,443 32,41 37,296 8 pt 36 197 74,39 35,209 39,181 109,599 9 ti 45 255 129,023 54,884 74,139 183,907 10 at 32 60 41,341 5,251 36,09 46,592 11 sg 31 62 42,823 4,836 37,987 47,659 12 sn 45 111 70,153 13,998 56,155 84,151 13 hs 28 58 33,738 4,249 29,489 37,987 14 pl 36 99 66,579 10,762 55,817 77,341 No Kelas LU
Band 4 Minimum Maximum Mean Stand. Deviation b
c 1 tg 40 87 60,592 9,943 50,649 70,535 2 sm 59 104 85,126 8,416 76,71 93,542 3 pm 36
85 52,166 9,498
42,668 61,664 4 tj 24
68 34,021 9,306 24,715 43,327
5 to 62 112 94,645 10,564 84,081 105,209 6 si 41 85 73,026 8,361 64,665 81,387 7 kc 51 101 74,362 9,506 64,856 83,868
8 pt 44 99 68,89 13,281
55,609 82,171 9 ti 33 105 58,023 17,464 40,559 75,487
10 at 16 84 34,093 15,827 18,266 49,92
11 sg 12 108 43,624 19,558 24,066 63,182
12 sn 35 96 66,114
12,853
53,261 78,967
13 hs 55 92 76,082 7,282 68,8 83,364 14 pl 18 70 36,421 7,000
29,421 43,421

## 106
Bhumi Vol. 3 No. 1 Mei 2017
Tabel 7. Matriks partisi untuk 14 kelas penggunaan lahan

## Sumber: Pengolahan penelitian, 2013
Dalam mengisi baris-baris tabel matrik partisi memperhatikan hubungan tingkat keterpisahan kelas yang dapat dilihat pada tabel 2. Sebagai contoh kelas penggunaan lahan tegalan/ladang (tg) memiliki keterpisahan yang kurang baik terhadap kelas penggunaan lahan permukiman (pm), peternakan (pt) dan sawah non irigasi (sn), maka baris-baris tersebut yang nantinya diisikan bobot atau atau nilai derajat keanggotaan, besarnya nilai bobot dalam penelitian ini berdasarkan pertimbangan lo- cal knowledge wilayah kajian juga dibantu dengan diagram pencar ( feature space ), grafik kosinus dan citra kelas fuzzy-nya. Caranya yaitu dengan bantuan citra kelas fuzzy dan local knowledge wilayah kajian, diambil suatu luasan kajian yang mencakup kelas penggunaan lahan yang akan ditentukan nilai derajat keanggotaannya, maka yang lebih dominan pada luasan tersebut yang diberikan bobot yang lebih besar, sedangkan yang tidak dominan diberi bobot yang lebih rendah. Jumlah total setiap baris dari matriks partisi tersebut sama dengan 1.
Matriks partisi yang sudah dibuat kemudian dihubungkan dengan sampel (training site) yang telah ditentukan sebelumnya. Training site yang telah memiliki nilai keanggotaan baru, kemudian dilakukan signature development lewat fuzsig menjadi daerah contoh baru menjadi training area fuzzy yang akan digunakan dalam eksekusi klasifikasi atau proses defuzzifikasi.
c. Proses eksekusi klasifikasi menggunakan algoritma maximum likelihood
Dari citra hasil proses klasifikasi digital yang masih berbasis fuzzy dilakukan proses defuzzifikasi dengan maksud untuk mempertegas (memperta- jam) batas-batas antar kelasnya.
Pada proses eksekusi ini training area fuzzy digunakan sebagai input, yang pembobotannya secara otomatis dihitung oleh komputer apabila kita memilih option use equal prior probabilities for each signature. Bobot yang diberikan pada pene- litian ini yaitu 0,071, nilai bobot tersebut diberikan sama nilainya untuk setiap kelas penggunaan lahan. Untuk melihat seberapa efisienkah proses ini maka diperlukan pencatatan waktu, dari hasil pencatatan waktu untuk proses eksekusi menggunakan algorit- ma maximum likelihood sekitar 10 detik. Software yang digunakan pada pemrosesan ini yaitu menggu- nakan software Idrisi Selva. Citra hasil eksekusi klasifikasi menggunakan algoritma maximum like- lihood dapat dilihat pada gambar 9. Sedangkan Peta Penggunaan Lahan hasil defuzzifikasi meng- gunakan algoritma maximum likelihood, yang telah diolah legenda dan proses layout -nya pada software ArcGIS  dapat dilihat pada gambar 10.
Gambar 9. Citra hasil eksekusi klasifikasi menggunakan algoritma maximum likelihood pada 14 kelas penggunaan lahan
ID tg sm Pm tj to Si Kc pt Ti at Sg sn hs pl tg 0,65 0 0,05 0 0 0 0 0,2 0 0 0 0,1 0 0 sm 0 0,6 0 0 0 0 0,05 0 0 0 0 0 0,35 0 pm 0 0 0,7 0 0 0 0 0,2 0 0 0 0 0 0,1 tj 0 0 0,04 0,72 0 0 0 0 0 0 0 0 0 0,26 to 0 0 0 0 0,6 0 0 0,4 0 0 0 0 0 0 si 0 0 0 0 0 0,89 0,11 0 0 0 0 0 0 0 kc 0 0 0 0 0 0,05 0,95 0 0 0 0 0 0 0 pt 0 0 0,15 0 0 0 0 0,8 0,05 0 0 0 0 0 ti 0 0 0,1 0,15 0 0 0 0,13 0,62 0 0 0 0 0 at 0 0 0 0 0 0 0 0 0 0,75 0,25 0 0 0 sg 0 0 0 0 0 0 0 0 0 0,16 0,84 0 0 0 sn 0,1 0 0,18 0 0 0 0 0,1 0 0 0 0,62 0 0 hs 0 0,18 0 0 0 0 0,3 0 0 0 0 0 0,52 0 pl 0 0 0,2 0,15 0 0 0 0 0 0 0 0 0 0,65

## 108 Bhumi Vol. 3 No. 1 Mei 2017

## 3. Tahap Pasca Lapangan (Uji Akurasi)
Uji akurasi hasil klasifikasi dilakukan untuk menguji tingkat akurasi peta penggunaan yang dihasilkan dari proses klasif ikasi digital dengan sampel uji dari hasil kegiatan lapangan. Antara sampel yang digunakan sebagai training area dengan sampel yang digunakan untuk uji akurasi bukan sampel yang sama tetapi sampel uji akurasi diambil di tempat yang berbeda, hal ini agar lebih diterima keakuratannya.
Metode yang digunakan untuk menghitung akurasi klasif ikasi dengan menggunakan matriks kesalahan atau confusion matrix/error matrix dapat dilihat pada tabel 8 , untuk selanjutnya menurut Jensen (2005) dapat dilakukan perhitungan producer’s accuracy , user’s accuracy, overall accuracy dan nilai indeks kappa.
Tabel 8 Matriks kesalahan ( error matrix )
Keterangan: A,B,C,……..., G = jumlah piksel benar dari hasil interpretasi dan cek lapangan.
a,b,c...,a’,b’,…,p’ = jumlah piksel dalam satu kelas hasil pengujian lapangan
Producer’s accuracy merupakan akurasi yang dilihat dari sisi penghasil peta, sedangkan user’s accuracy merupakan akurasi yang dilihat dari sisi pengguna petanya. Pada tabel perhitungan juga terdapat istilah omisi kesalahan yaitu kesalahan karena adanya penghilangan, sebaliknya komisi kesalahan yaitu kesalahan karena adanya penam- bahan.
Uji akurasi dihitung dengan menggunakan rumus producer’s accuracy dan user’s accuracy yang dapat dilihat pada tabel 8. Selain itu juga akan dila- kukan perhitungan overall accuracy dengan rumus:
Tabel 8. Perhitungan akurasi penghasil dan pengguna berdasarkan matriks kesalahan
a. Hasil perhitungan producer’s accuracy dan user’s accuracy 14 kelas penggunaan lahan
Hasil perhitungan producer’s accuracy dan user’s accuracy 14 kelas penggunaan  lahan dapat dilihat pada tabel 9.
Tabel 9. Perhitungan akurasi penghasil dan pengguna untuk 14 kelas penggunaan lahan (defuzzif ikasi menggunakan maximum likelihood )
Kelas
Producer’s accuracy
User’s accuracy Akurasi (%) Omisi kesalahan (%) Akurasi (%) Komisi kesalahan (%) Tg 38,80 61,20 27,73 72,27 Sm 75,52 24,48 65,47 34,53 Pm 49,77 50,23 56,65 43,35 Tj 42,76 57,24 70,45 29,55 To 90,16 9,84 94,83 5,17 Si 17,73 82,27 15,34 84,66 Kc 29,52 70,48 24,14 75,86 Pt 4,46 95,54 4,00 96,00 Ti 57,24 42,76 78,30 21,70 At 58,19 41,81 87,29 12,71 Sg 73,54 26,46 80,77 19,23 Sn 71,16 28,84 64,02 35,98 Hs 79,55 20,45 75,92 24,08 Pl 33,33 66,67 30,95 69,05
Harvini Wulansari: Uji Akurasi Klasif ikasi Penggunaan Lahan ... : 98-110
Dari tabel 9 dapat dilihat bahwa akurasi tertinggi untuk 14 kelas penggunaan lahan pada proses defuzzifikasi menggunakan maximum likelihood tingkat akurasi penghasil atau producer’s accuracy sebesar 90,16 % yaitu kelas penggunaan lahan tempat olahraga (to) dengan omisi kesalahan sebesar 9,84% sedangkan akurasi terendah sebesar 4,46 % yaitu kelas penggunaan lahan peternakan dengan omisi kesalahan sebesar 95,54%, Pada tingkat akurasi pengguna atau user’s accuracy akurasi tertinggi sebesar 94,83 % yaitu kelas peng- gunaan lahan tempat olahraga dengan komisi kesa- lahan sebesar 5,17%, sedangkan akurasi terendah sebesar 4% yaitu kelas penggunaan lahan peter- nakan dengan komisi kesalahan sebesar 96%.
a. Hasil Perhitungan Overall Accuracy dan Indeks Kappa
Hasil perhitungan overall accuracy dan indeks kappa untuk 14 kelas penggunaan lahan proses defuzzif ikasi maximum likelihood dapat dilihat pada tabel 10.
Tabel 10. Perhitungan overall accuracy dan indeks kappa dengan proses defuzzif ikasi menggunakan algoritma maximum likelihood
Sumber: Pengolahan penelitian, 2013

## C. Penutup
Proses penelitian berhasil dengan baik, meng- hasilkan peta yang informatif dengan klasifikasi yang bervariasi, walaupun dari sudut pandang ketelitian menghasilkan overall accuracy dan indeks kappa yang kurang baik atau kurang layak apabila menggunakan acuan akurasi yang dianggap layak untuk penggunaan citra satelit 80% - 85% (Sutan- to,2013). Hasil perhitungan overall accuracy untuk 14 kelas penggunaan lahan proses defuzzif ikasi maximum likelihood yaitu 57% sedangkan indeks kappa 0,53 namun demikian hasilnya masih dapat
diterima. Sedangkan hasil pencatatan waktu untuk proses eksekusi menggunakan algoritma maximum likelihood sekitar 10 detik.
Hasil akurasi yang kurang baik dalam penelitian ini bisa disebabkan oleh beberapa kemungkinan diantaranya yaitu dalam interpretasi sampel di citra dengan kondisi sebenarnya di lapangan ada perbedaan disebabkan karena citra yang digunakan merupakan hasil perekaman tahun 2009 sedang- kan penelitian dilakukan tahun 2013 tentu saja ada kemungkinan kondisi lapangan sudah berubah meskipun wawancara yang dilakukan dengan pen- duduk sekitar sudah dilakukan tetapi belum tentu informasi tersebut akurat karena rentang waktu yang cukup lama.
Selain itu kemungkinan dalam hal penentuan derajat keanggotaan himpunan fuzzy dibutuhkan analisis yang baik oleh orang yang benar benar ahli- nya, karena dalam penelitian ini peneliti masih pemula dalam melakukan analisis/bukan ahlinya maka bisa saja terjadi kesalahan dalam pemberian bobot-bobot ( membership grade ) keanggotaan.

## Daftar Pustaka
Danoedoro, P., 2012, Pengantar Penginderaan Jauh Digital , Penerbit Andi, Yogyakarta.
Denaswidhi, E., 2008, Pendekatan Logika Samar (Fuzzy Logic) untuk Pemetaan Penggunaan Lahan Sebagian Kota Semarang dari Citra Landsat ETM+ Multitemporal , Skripsi, Fakultas Geograf i Universitas Gadjah Mada Yogyakarta.
Direktorat Pemetaan Tematik, 2012, Norma, Standar, Prosedur dan Kriteria , Deputi Bidang Survei, Pengukuran dan Pemetaan, BPN RI. Jensen, J.R., 2005, Introductory Digital Image Pro- cessing : A Remote Sensing Perspective, Third Edition , Pearson Education, Inc., United States of America.
Kusumadewi, S., dan Purnomo, H., 2010, Aplikasi Logika Fuzzy Untuk Pendukung Keputusan, Edisi 2, Penerbit Graha Ilmu, Yogyakarta.
Jumlah kelas overall accuracy indeks kappa 14 kelas 57% 0,53

## 110 Bhumi Vol. 3 No. 1 Mei 2017
Lillesand, T.M., Kiefer, R.W.,  and Chipman, J., 2008, Remote Sensing and Image Interpretation , Sixth Edition , John Wiley and Sons, Inc.
Novamizanti et.al., 2010, Optimasi Logika Fuzzy Menggunakan Algoritma Genetika Pada Iden- tif ikasi Pola Tanda Tangan, Jurnal, Konferensi Nasional Sistem dan Informatika Bali, hal: 153. Schowengerdt, R.A., 1983, Techniques for Image Pro- cessing and Classif ication in Remote Sensing ,
Academic Press, London.
Sutanto, 2010, Remote Sensing Research: A User’s Perspective , Indonesian Journal of Geography,
Faculty of Geography Gadjah Mada Univ. and The Indonesian Geographers Association, hal: 131.
Sutanto, 2013, Metode Penelitian Penginderaan Jauh , Penerbit Ombak, Yogyakarta.
http://www.ga.gov.au/webtemp/image_cache / GA10286.pdf ý, diakses tanggal 21 September 2013