PENYUSUNAN MODEL PENDUGAAN DAN PEMETAAN BIOMASSA PERMUKAAN PADA TEGAKAN JATI (Tectona grandis Linn.F) MENGGUNAKAN CITRA ALOS PALSAR RESOLUSI 50 M DAN 12,5 M (Studi Kasus : KPH Kebonharjo Perhutani Unit I Jawa Tengah)
RISA DESIANA SYARIF
DEPARTEMEN MANAJEMEN HUTAN FAKULTAS KEHUTANAN INSTITUT PERTANIAN BOGOR 2011
RINGKASAN RISA DESIANA SYARIF. Penyusunan Model Pendugaan dan Pemetaan Biomassa Permukaan Pada Tegakan Jati (Tectona grandis Linn. F) Menggunakan Citra Alos Palsar Resolusi 50 M dan 12,5 M (Studi Kasus : KPH Kebonharjo Perhutani Unit I Jawa Tengah). Dibimbing oleh NINING PUSPANINGSIH dan M. BUCE SALEH. Hutan sebagai sistem dinamis disusun oleh berbagai komponen salah satunya adalah biomassa. Biomassa memiliki peranan penting dalam pengelolaan ekosistem hutan dan menjadi salah satu parameter yang digunakan dalam mengetahui perubahan struktur hutan, karena jumlah stok biomassa bergantung pada ada atau tidaknya permudaan alam, terganggu atau tidaknya hutan, dan peruntukan hutan (IPCC 1995). Penelitian yang menghasilkan pendekatan untuk mengestimasi biomassa telah banyak dilakukan, salah satu diantaranya adalah penggunaan informasi dijital pada data SAR (Synthetic Aperture Radar). Pada tahun 2006, diluncurkan sensor baru dengan sistem SAR pada satelit ALOS milik pemerintah Jepang, yaitu sensor PALSAR. Penelitian ini bertujuan untuk menyusun model pendugaan biomassa pada hutan tanaman jati (Tectona grandis) di KPH Kebonharjo perhutani unit I Jawa Tengah, dengan menggunakan citra ALOS PALSAR resolusi spasial 50 m dan 12,5 m serta pembuatan peta sebaran dengan model terpilih. Pada citra ALOS PALSAR dilakukan ekstraksi nilai dijital yang kemudian dikalibrasi menggunakan persamaan NRCS (Normalize Radar Cross Section, dB) yang dikembangkan oleh Shimada et al. (2009) untuk menghasilkan nilai backscatter. Nilai backscatter citra ALOS PALSAR kemudian diregresikan mengguanakan kaidah non-linier dengan biomassa dimana perhitungan biomassa dilakukan menggunakan model alometrik untuk tegakan jati dan koefisien BEF (Biomass Expansion Factor). Model terpilih diperoleh dengan memperhatikan nilai R2adj dan Root mean square error (RMSE). Verifikasi pada model terpilih dilakukan menggunakan uji T berpasangan untuk mengetahui kelayakan model. Pemetaan biomassa dibuat berdasarkan model-model yang telah diverifikasi. Pada peta sebaran biomassa dilakukan perhitungan overall accuracy dan kappa accuracy. Nilai biomassa minimum pada tegakan jati yang dihitung dengan menggunakan alometrik tegakan jati adalah sebesar 7,52 ton/ha pada KU I dan nilai biomassa maksimum sebesar 328,69 ton/ha pada KU XI. Perhitungan dengan menggunakan BEF cenderung overestimate pada kelas umur tua sehingga pada penelitian ini pembuatan model pendugaan dengan BEF tidak dilakukan. Berdasarkan hasil regresi antara biomassa dan backscatter didapatkan 16 model penduga biomassa. Dari 16 model yang telah dibuat, dipilih 4 model terbaik. Verifikasi dilakukan pada 4 model terpilih menggunakan uji T berpasangan. Verifikasi menunjukan bahwa ke 4 model tersebut memiliki nilai T hitung < T tabel dan signifikansi > 0,05 sehingga model layak digunakan untuk menduga biomassa di daerah penelitian. Model yang menghasilkan peta sebaran biomassa dengan akurasi terbaik pada citra ALOS PALSAR resolusi 50 m adalah model 2 sedangkan pada citra ALOS PALSAR resolusi 12,5 adalah model 4. Kata kunci : biomassa, backscatter, ALOS PALSAR, pemetaan.
SUMMARY RISA DESIANA SYARIF. Above-Ground Biomass Estimation and Mapping in Teak (Tectona grandis Linn. F) Plantation Area using ALOS PALSAR Image with 12,5 m and 50 m Spatial Resolution (Case Study: KPH Kebonharjo, Perhutani Unit I, Central Java). Supervised by NINING PUSPANINGSIH and M. BUCE SALEH. Forest as a dynamic system is consisted by many components, one of its components is biomass. Biomass can be used to determine the forest structure and condition, since biomass is depend on natural regeneration of the forest, disturbed state of the forest, and the forest use (IPCC 1995). Researchs than estimate biomass has been done using various methods, one of these methods is using SAR data in estimating biomass. In year 2006 Japanese government launched a SARbased satellite, ALOS satellite. PALSAR sensor in ALOS satellite can produce SAR-based image that can be used in estimating biomass. This research aimed to create biomass estimation model and biomass distribution map for teak plantation area in KPH Kebonharjo Unit I Central Java, using ALOS PALSAR image with 12,5 m and 50 m spatial resolution. Backscatter value is generated using digital number from ALOS PALSAR image and calibrated using Normalized Radar Cross Section (NRCS), which developed by Shimada et al. (2009). Therefore the backscatter value is then regressed using non-linear rulings with biomass that is estimated from teak allometric equations and Biomass Expansion Factor (BEF) coefficient. Certain model is selected by ranking the R2adj and Root Mean Square Error (RMSE). Chosen model is then verified using Paired T-test to analyze the appropriateness of the models. Biomass distribution map is constructed based on verified models, and then assessed using two accuracy assessment, Overall Accuracy and Kappa Accuracy. Minimum estimated biomass value in teak stands area using allometric equation is numbered at 7,52 ton/ha in Age Class (KU) I and maximum estimated biomass value in teak stands area using allometric equation is numbered at 328,69 ton/ha in KU XI. The estimation of biomass that calculated using BEF coefficient are overestimate. This research focused on estimating biomass using Hendri allometric equation. 16 models are generated based on biomass and backscatter regression. From 16 models, 4 models are chosen. Verification using Paired T-test on the 4 models shows that all the models having Tscore < Ttable and significance value > 0,05. Model that used in construction of biomass distribution map which has the best accuration for 50 m spatial resolution is model 2, meanwhile best accuration for 12,5 spatial resolution is model 4.
Keyword: biomass, backscatter, ALOS PALSAR, mapping.
PENYUSUNAN MODEL PENDUGAAN DAN PEMETAAN BIOMASSA PERMUKAAN PADA TEGAKAN JATI (Tectona grandis Linn.F) MENGGUNAKAN CITRA ALOS PALSAR RESOLUSI 50 M DAN 12,5 M (Studi Kasus : KPH Kebonharjo Perhutani Unit I Jawa Tengah)
RISA DESIANA SYARIF
Skripsi Sebagai salah satu syarat untuk memperoleh gelar Sarjana Kehutanan pada Fakultas Kehutanan Institut Pertanian Bogor
DEPARTEMEN MANAJEMEN HUTAN FAKULTAS KEHUTANAN INSTITUT PERTANIAN BOGOR 2011
PERNYATAAN Dengan ini saya menyatakan bahwa skripsi berjudul Penyusunan Model Pendugaan dan Pemetaan Biomassa Permukaan Pada Tegakan Jati (Tectona grandis Linn. F) Menggunakan Citra ALOS PALSAR Resolusi 50 M dan 12,5 M (Studi Kasus : KPH Kebonharjo Perhutani Unit I Jawa Tengah) adalah benarbenar hasil karya saya sendiri dengan bimbingan dosen pembimbing dan belum pernah digunakan sebagai karya ilmiah pada perguruan tingga manapun. Sumber informasi yang berasal atau dikutip dari karya yang diterbitkan maupun tidak diterbitkan dari penulis lain telah disebutkan dalam teks dan dicantumkan dalam Daftar Pustaka di bagian akhir skripsi ini.
Bogor, Juni 2011
Risa Desiana Syarif NRP E14052331
Judul Penelitian
: Penyusunan Model Pendugaan dan Pemetaan Biomassa Permukaan Pada Tegakan Jati (Tectona grandis Linn. F) Menggunakan Citra Alos Palsar Resolusi 50 M dan 12,5 M (Studi Kasus : KPH Kebonharjo Perhutani Unit I Jawa Tengah)
Nama Mahasiswa : Risa Desiana Syarif Nomor Pokok
: E14052331
Menyetujui Komisi Pembimbing
Ketua,
Anggota,
Dr. Nining Puspaningsih, M.Si. NIP. 19630612 199003 2 014
Dr. M. Buce Saleh, MS. NIP. 19571005 198303 1 002
Mengetahui, Ketua Departemen Manajemen Hutan Fakultas Kehutanan IPB
Dr. Ir. Didik Suharjito, MS. NIP. 19630401 199403 1 001 Tanggal Lulus:
KATA PENGANTAR
Segala puji dan syukur kehadirat Allah SWT yang telah memberikan kemudahan dan jalan kepada penulis dalam menyelesaikan skripsi yang berjudul Penyusunan Model Pendugaan Biomassa Permukaan pada Tegakan Jati (Tectona grandis Linn. F) menggunakan citra ALOS PALSAR resolusi 50 m dan 12,5 m (Studi Kasus : KPH Kebonharjo Perhutani Unit I Jawa Tengah). Skripsi ini disusun sebagai salah satu syarat untuk memperoleh gelar Sarjana Kehutanan di Fakultas Kehutanan IPB. Skripsi yang telah disusun ini menggambarkan pola hubungan antara backscatter citra ALOS PALSAR pada resolusi yang berbeda, yaitu citra ALOS PALSAR resolusi 50 m dan 12,5 m dengan biomassa tegakan jati yang dihitung berdasarkan model alometriknya. Dimana kemudian dibuatlah peta sebaran biomassa pada daerah penelitian berdasarkan pola hubungan tersebut. Penulis menyadari dalam penulisan skripsi ini masih belum sempurna. Oleh karena itu saran, kritik, dan masukan guna memperbaiki tulisan ini sangat diharapkan. Semoga tulisan ini dapat bermanfaat bagi dunia pendidikan, khususnya dunia kehutanan.
Bogor, Juli 2011
Penulis
DAFTAR ISI Halaman KATA PENGANTAR ................................................................................... i DAFTAR ISI .................................................................................................. iii DAFTAR GAMBAR ..................................................................................... v DAFTAR TABEL .......................................................................................... vii BAB I PENDAHULUAN .............................................................................. 1 1.1 Latar Belakang .......................................................................... 1 1.2 Tujuan Penelitian ...................................................................... 2 BAB II TINJAUAN PUSTAKA................................................................... 3 2.1 Teknologi Penginderaan Jauh ................................................... 3 2.2 Citra Satelit Sistem RADAR .................................................... 6 2.2.1 Citra Satelit ALOS PALSAR .......................................... 10 2.3 Sistem Informasi Geografi ........................................................ 13 2.4 Jati (Tectona grandis Linn F.) .................................................. 15 2.5 Biomassa ................................................................................... 16 BAB III METODOLOGI .............................................................................. 21 3.1 Lokasi dan Waktu Penelitian .................................................... 21 3.2 Bahan dan Alat.......................................................................... 21 3.3 Tahapan Pelaksanaan ................................................................ 22 3.3.1 Persiapan dan Pengambilan Data Lapangan ................... 22 3.3.2 Bentuk dan Ukuran Plot Contoh ..................................... 22 3.3.3 Pengambilan Data Lapangan ........................................... 24 3.3.4 Pengolahan Data .............................................................. 25 3.3.5 Pembuatan Citra Pendugaan Biomassa dan Peta Kelas Biomassa ......................................................................... 28 BAB IV KONDISI UMUM LOKASI PENELITIAN .................................. 31 4.1 Risalah Umum KPH Kebonharjo ............................................. 31 4.1.1 Letak ................................................................................ 31 4.1.2 Batas Wilayah.................................................................. 31 4.1.3 Topografi ......................................................................... 33 4.1.4 Enclave ............................................................................ 33 BAB V HASIL DAN PEMBAHASAN ......................................................... 35 5.1 Rata-rata Biomassa di KPH Kebonharjo .................................. 35 5.2 Hasil Pengolahan Data Citra ..................................................... 35 5.3 Penyusunan Model Pendugaan Biomassa................................. 36 5.3.1 Biomassa Hendri dan Biomassa BEF .............................. 36 5.3.2 Pola Hubungan Backscatter Citra ALOS PALSAR dengan Biomassa ............................................................. 38 5.3.3 Penyusunan Model Hubungan Backscatter dengan Biomassa Alometrik Hendri ............................................ 39 5.3.4 Verifikasi Model.............................................................. 42 5.4 Pembuatan Peta Sebaran Biomassa dan Analisis Akurasi ........ 45 5.4.1 Peta Sebaran Biomassa .................................................... 46 5.4.2 Uji Akurasi Peta Sebaran Biomassa ................................ 47
iv
BAB VI KESIMPULAN DAN SARAN .................................................. .... 6.1 Kesimpulan............................................................................... DAFTAR PUSTAKA .................................................................................... LAMPIRAN ..................................................................................................
55 55 57 61
DAFTAR GAMBAR No. 1. 2. 3. 4. 5. 6. 7. 8. 9. 10. 11. 12. 13. 14. 15. 16.
Halaman Spektrum elektromagnetik (Purwadhi 2001)......................................... 5 Transmisi dan reflektansi pada radar (NASA 1996).............................. 7 Wahana dan arah tembak sensor (NASA 1996)..................................... 8 Mekanisme hamburan balik pada radar di setiap jenis permukaan (NASA 1996)....................................................................................................... 9 Mekanisme hamburan balik pada tegakan hutan (Dobson et al. 1990).. 10 Instrumen PALSAR (JAXA 2010)......................................................... 12 Prinsip geometri dari PALSAR (JAXA 2010)....................................... 12 Lokasi penelitian.................................................................................... 21 Diagram Alir Penelitian......................................................................... 23 Sketsa plot pengukuran.......................................................................... 24 Perbandingan nilai biomassa alometrik Hendri dan biomassa BEF...... 37 Scatterplot hubungan antara biomassa Hendri dengan backscatter ALOS PALSAR................................................................................................ 38 Peta sebaran biomassa model 1............................................................. 48 Peta sebaran biomassa model 2............................................................. 49 Peta sebaran biomassa model 3............................................................. 51 Peta sebaran biomassa model 4............................................................. 53
DAFTAR TABEL No. 1. 2. 3. 4. 5. 6.
Halaman Spektrum elektromagnet. ........................................................................ 6 Keterangan umum ALOS ....................................................................... 11 Karakteristik PALSAR ........................................................................... 13 Beberapa metode untuk menduga biomassa ........................................... 17 Model yang digunakan untuk pendugaan biomassa ............................... 26 Rata-rata biomassa di KPH Kebonharjo per KU dan jumlah plot pengamatan per KU berdasarkan alometrik tegakan jati dan BEF .................................................................................................. 35 7. Model pendugaan biomassa berdasarkan hubungan biomassa Hendri dengan backscatter HH citra ALOS PALSAR resolusi 50 m ........................................................................................... 40 8. Model pendugaan biomassa berdasarkan hubungan biomassa Hendri dengan backscatter HV citra ALOS PALSAR resolusi 50 m ........................................................................................... 40 9. Model pendugaan biomassa berdasarkan hubungan biomassa Hendri dengan backscatter HH citra ALOS PALSAR resolusi 12,5 m. ....................................................................................... 41 10. Model pendugaan biomassa berdasarkan hubungan biomassa Hendri dengan backscatter HV citra ALOS PALSAR resolusi 12,5 m. ....................................................................................... 41 11 . Hasil verifikasi model pendugaan biomassa pada citra ALOS PALSAR resolusi 50 m .......................................................................... 43 12. Hasil verifikasi model pendugaan biomassa pada citra ALOS PALSAR resolusi 12,5 m ........................................................................ 43 13. Model pendugaan biomassa terbaik ALOS PALSAR resolusi 50 m dan ALOS PALSAR resolusi 12,5 m ..................................................... 46 14 . Selang kelas biomassa yang digunakan untuk membuat peta sebaran biomassa .................................................................................... 46 15. Rekapitulasi hasil perhitungan overall accuracy dan kappa accuracy pada peta sebaran dengan citra ALOS PALSAR resolusi 50 m ........................................................................................... 47 16. Rekapitulasi hasil perhitungan overall accuracy dan kappa accuracy pada citra ALOS PALSAR resolusi 12,5 m ............................ 50
BAB I PENDAHULUAN 1.1 Latar Belakang Hutan sebagai sistem yang dinamis memiliki peranan yang sangat besar terhadap terjadinya perubahan lingkungan. Biomassa sebagai salah satu komponen yang terdapat dalam sistem dinamis hutan memiliki peranan penting dalam pengelolaan ekosistem hutan. Di negara-negara maju biomassa tegakan antara lain digunakan sebagai dasar pertimbangan bagi kegiatan pengelolaan hutan lestari, karena jumlah stok biomassa tergantung pada terganggu atau tidaknya hutan, ada atau tidaknya permudaan alam, dan peruntukan hutan (IPCC 1995). Informasi keadaan hutan diketahui melalui monitoring yang selalu dilakukan pada luasan hutan tertentu dimana monitoring yang dilakukan tidaklah mudah, monitoring dilakukan untuk memahami perubahan-perubahan yang terjadi pada komponen hutan seperti volume batang dan biomassa dimana keduanya berhubungan dan bersifat esensial serta kompleks (Kumar 2009). Di masa lalu informasi mengenai biomassa didapatkan secara manual melalui ground survey (survey lapangan), tetapi saat ini ketika informasi mengenai komponen hutan beserta perubahannya semakin banyak dibutuhkan maka pengambilan informasi dengan ground survey dianggap kurang efisien untuk dilakukan pada daerah pengelolaan yang relatif luas. Pengambilan informasi skala besar dan pada daerah pengelolaan hutan yang relatif luas dengan menggunakan teknologi penginderaan jauh telah dilakukan dan menunjang kebutuhan informasi yang berkembang dari waktu ke waktu. Teknologi penginderaan jauh telah digunakan di banyak bidang seperti identifikasi tutupan lahan, pemetaan, pembagian wilayah dan tipe hutan, serta pendugaan parameter tegakan hutan. Saat ini biomassa sebagai parameter yang penting dalam mengukur perubahan struktur hutan dapat diduga dengan menggunakan teknologi penginderaan jauh (Schreuder et al. (1997) dalam Husch et al. 2003). Banyak studi telah memberikan pendekatan yang berguna untuk
2
mengestimasi biomassa, salah satu diantaranya adalah penggunaan informasi dijital pada data SAR (Synthetic Aperture Radar). Dalam perkembangan teknologi penginderaan jauh banyak dilakukan penelitian terhadap data SAR dalam mengestimasi biomassa. Sensitifitas SAR terhadap struktur tegakan hutan serta hubungan antara backscatter pada radar dengan biomassa dan parameter tegakan yang lain telah banyak diteliti di berbagai daerah tropis maupun subtropis. Diprediksikan, dengan pengembanganpengembangan yang terus dilakukan pada data SAR memungkinkan pendugaan biomassa akan memiliki keterwakilan yang lebih baik (Sarker dan Nichol 2010). Teknologi
penginderaan
jauh
merupakan
teknologi
yang
dapat
dikembangkan dari waktu ke waktu dalam mempermudah pendugaan biomassa, walaupun demikian keberadaan teknologi ini tidak dapat berdiri sendiri, pendekatan-pendekatan
alometrik
yang
telah
dikembangkan
sebelumnya
kemudian mendukung efektifitas penggunaan teknologi penginderaan jauh dalam melakukan pendugaan biomassa (Brown et al. 1989). Hutan tanaman Jati merupakan salah satu pemanfaatan kawasan hutan yang produknya memiliki sumbangsih dalam memajukan perekonomian negara. Selain fungsi ekonomi, Hutan tanaman Jati memiliki peran penting dalam menjaga kestabilan iklim mikro didaerah setempat dan sekaligus memiliki fungsi penyimpanan karbon. Berkaitan dengan upaya pencegahan pemanasan global maka dilakukan studi mengenai kandungan karbon pada tegakan di hutan tropis salah satunya adalah hutan tanaman Jati. Untuk mengetahui kandungan karbon pada suatu tegakan Jati, dilakukan dengan menggunakan pendekatan biomassa, dimana pendekatan biomassa merupakan ukuran yang paling sering digunakan untuk menggambarkan dan mempelajari pertumbuhan tanaman. 1.2 Tujuan Penelitian Penelitian ini bertujuan untuk menyusun model pendugaan biomassa pada hutan tanaman Jati (Tectona grandis) di KPH Kebonharjo Perhutani Unit I, Jawa Tengah, menggunakan citra ALOS PALSAR resolusi spasial 50 m dan 12,5 m serta pembuatan peta sebaran biomassa dari model yang terpilih.
BAB II TINJAUAN PUSTAKA 2.1 Teknologi Penginderaan Jauh Penginderaan jauh adalah ilmu dan seni untuk memperoleh informasi tentang suatu objek, daerah, atau fenomena melalui analisis data yang diperoleh dengan suatu alat tanpa kontak langsung dengan objek, daerah, atau fenomena yang dikaji (Lillesand & Kiefer 1979). Fussel et al. (1986) mendefinisikan penginderaan jauh sebagai pengumpulan dan pencatatan informasi tanpa kontak langsung pada julat elektromagnetik ultraviolet, tampak, inframerah, dan mikro dengan menggunakan peralatan seperti penyiam (scanner) dan kamera yang ditempatkan pada wahana bergerak seperti, pesawat udara atau pesawat angkasa, dan menganalisis informasi yang diterima dengan teknik interpretasi foto, citra dan pengolahan citra. Teknik penginderaan jauh berkembang sangat pesat sejak diluncurkannya satelit penginderaan jauh ERTS (Earth Resources Technology Satellite) pada 1972. Hal itu memungkinkan pengumpulan data permukaan bumi dalam jumlah besar. Perkembangan sistem penginderaan jauh khususnya dalam penggunaan sensor
dan
cara
perekaman
data
penginderaan,
telah
diikuti
dengan
pengembangan dalam cara pengolahan dan analisis data penginderaan jauh. Pengembangan penginderaan jauh yang terjadi tahun 1960-an terbatas pada penelitian serta analisis foto udara multispectral scanner (Perangkat pengindera yang terdiri atas kurang dari sepuluh spectrum elektromagnetik yang berbeda) dan dijitalisasi dari foto udara yang diambil dengan sensor kamera. Semenjak ERTS-1 atau Landsat-1 (Land Satellite) diluncurkan tahun 1972, data citra dijital rekaman satelit menjadi semakin penting, khususnya untuk penggunaan analisis permukaan bumi (Purwadhi 2001).
4
Samsuri (2004) menyatakan bahwa dalam penginderaan jauh terdapat beberapa proses melibatkan interaksi antara radiasi dan target yang dituju mencakup tujuh elemen penting yakni: 1. Sumber Energi atau Illumination, merupakan elemen pertama dalam menyediakan energi elektromagnetik ke target interes 2. Radiasi dan Atmosfer, adalah perjalanan energi dari sumber ke targetnya dan sebaliknya. Energi akan mengalami kontak dengan target dan berinteraksi dengan atmosfer yang dilewatinya. 3. Interaksi dengan Target. 4. Perekaman Energi oleh Sensor setelah energi dipancarkan atau dilepaskan dari target, elemen penting yang dibutuhkan adalah sensor
untuk
mengumpulkan
dan
merekam
radiasi
elektromagnetik. 5. Transmisi, penerimaan dan pemrosesan energi yang terekam oleh sensor
harus
ditransmisikan
untuk
diterima
oleh
stasiun
pengolahan, dimana data diolah menjadi citra (hardcopy ataupun digital) 6. Interpretasi dan Analisis, merupakan pengolahan image dengan interpretasi secara visual atau digital untuk mengekstrak informasi tentang target. 7. Aplikasi, elemen terakhir adalah pengaplikasian informasi tentang target untuk memperoleh pengertian yang lebih baik, menerima beberapa informasi baru, dan membantu pemecahan masalah. Sensor adalah alat perekam objek bumi. sensor dipasang pada wahana (platform) dan letaknya jauh dari objek yang diindera, maka diperlukan tenaga elektromagnetik yang dipancarkan atau dipantulkan oleh objek tersebut. sensor elektronik membangkitkan sinyal elektrik yang sesuai dengan variasi tenaga elektromagnetik. Tenaga radiasi yang dipantulkan oleh objek dan ditangkap oleh sensor, dapat menghasilkan citra yang sesuai dengan wujud aslinya. setiap sensor mempunyai kepekaan spektral yang terbatas, batas kemampuan sensor untuk memisahkan objek disebut resolusi. Resolusi suatu citra merupakan indikator tentang kemampuan sensor dan atau kualitas sensor dalam merekam objek. Empat
5
resolusi yang biasa digunakan sebagai parameter kemampuan sensor, yaitu resolusi spasial, resolusi spektral, resolusi radiometrik, dan resolusi temporal (Purwadhi, 2001). Data penginderaan jauh digital (citra digital) direkam menggunakan sensor non-kamera, antara lain scanner, radiometer, spectometer. Detektor yang digunakan dalam sensor penginderaan jauh adalah detektor elektronik dengan menggunakan tenaga elektromagnetik yang luas, yaitu spektrum tampak, ultraviolet, inframerah dekat, inframerah termal, dan gelombang mikro (Gambar 1). Citra digital dibentuk dari elemen-elemen gambar atau pixel (picture element) yang menyatakan tingkat keabuan pada gambar. Citra bersifat digital dapat secara langsung disimpan pada suatu pita magnetik. Citra digital dapat direkam dalam beberapa spektrum secara sekaligus sehingga disebut citra multispektral.
Gambar 1 Spektrum elektromagnetik (Purwadhi 2001). Sumber data citra penginderaan jauh digital saat ini mulai berkembang dan bervariasi, yaitu data dapat diperoleh dari rekaman satelit penginderaan jauh komersial, rekaman data satelit meteorologi, rekaman scanner dari pesawat udara, data digital dari kamera foto udara, citra digital dari rekaman sensor mikro densitometer, dan citra digital dari hasil rekaman beberapa sistem dijitasi dengan resolusi tinggi (Purwadhi 2001). Komponen dasar pengambilan data penginderaan jauh sistem pasif meliputi sumber tenaga, atmosfer, interaksi tenaga dengan objek dipermukaan bumi, sensor, sistem pengolahan data, dan berbagai penggunaan data. Penginderaan jauh sistem aktif menggunakan tenaga elektromagnetik yang dibangkitkan oleh sensor Radar (Radio Detection And Ranging).
6
Tabel 1 Spektrum elektromagnet. Spektrum/saluran Panjang gelombang Gamma
≤ 0,03 nm
X
0.03 – 3 nm
Ultraviolet (UV) UV Fotografik
0,03 - 0,4 μm 0,3 – 0,4 μm
Tampak Biru Hijau Merah Inframerah (IM) NIR SWIR MWIR LWIR FIR Gelombang Mikro
0,4 – 0,7 μm 0,4 – 0,5 μm 0,5 – 0,6 μm 0,6 – 0,7 μm 0,7 – 300 μm 0,7 – 1,5 μm 1,5 – 3 μm 3 – 8 μm 8 – 15 μm ≥ 15 μm 8 – 14 μm
Radar
0,75 – 300 cm Ka K Ku X C S L P
Keterangan Diserap oleh atmosfer, benda radioaktif dapat diindera dari pesawat terbang rendah. Diserap oleh atmosfer, sinar buatan digunakan dalam kedokteran. Radiasi UV diserap oleh ozon. UV Fotografik dapat direkam oleh film fotografi, diperlukan lensa kuarsa dalam kamera.
Jendela atmosfer terpisah oleh saluran absorpsi. Film khusus dapat merekam hingga panjang gelombang hampir 1,2 μm. MWIR dan LWIR dikenal dengan sebagai thermal infrared (TIR). Gelombang panjang yang mampu menembus awan, citra dapat dibuat dengan cara aktif. Penginderaan jauh sistem aktif. Yang paling sering digunakan.
0,75 – 1,1 cm 1,1 – 1,7 cm 1,7 – 2,4 cm 2,4 – 3,8 cm 3,8 – 7,5 cm 7,5 – 15 cm 15 – 30 cm 30 – 100 cm
Radio Sumber : Prahasta (2009)
Tidak digunakan dalam penginderaan jauh.
2.2 Citra Satelit Sistem RADAR Radar (Radio Detecting and Ranging) dikembangkan sebagai suatu cara untuk mendeteksi adanya objek dan menentukan posisi objek tersebut dengan menggunakan radio. Karena penginderaan jauh sistem radar merupakan penginderaan jauh sistem aktif, tenaga elektromagnetik yang digunakan di dalam penginderaan jauh dibangkitkan pada sensor. Tenaga ini berupa pulsa bertenaga tinggi yang dipancarkan dalam waktu yang sangat pendek yaitu sekitar 10 -6 detik (Purwadhi 2001). Antena pada radar mentransmisi dan menerima gelombang (pulsa) pada panjang gelombang dan polarisasi tertentu. Energi gelombang radar menyebar ke seluruh bagian permukaan bumi, dengan sebagian energi yang dikenal sebagai backscatter atau hamburan balik.
7
Gambar 2 Transmisi dan reflektansi pada radar (NASA 1996). Hamburan balik ini dipantulkan kembali pada radar sebagai pantulan gelombang radar yang lemah dan diterima oleh antena pada bentuk polarisasi tertentu (horizontal atau vertikal, tidak selalu sama dengan yang ditransmisikan). Pantulan gelombang tersebut dikonversikan menjadi data dijital dan dikirim ke perekaman data kemudian ditampilkan menjadi image (citra satelit). Biasanya lama waktu sebuah gelombang sampai pada objek digunakan sebagai penghitung jarak ke objek (bandwidth). Semakin besar bandwidth semakin baik resolusi yang dihasilkan pada dimensi objek tersebut. Panjang antena radar menentukan resolusi pada image searah azimuth, semakin panjang antena semakin baik resolusi yang di hasilkan. Synthetic Aperture Radar (SAR) menunjuk pada sebuah teknik yang digunakan untuk mensintetis antena yang sangat panjang dengan mengombinasikan sinyal yang diterima radar yang bergerak pada jalur terbangnya. Aperture berarti pembukaan yang terjadi dalam proses penyerapan refleksi energi yang digunakan dalam pembuatan gambar, sebagai contoh, pada kamera pembukaan ini berarti pembukaan lensa kamera, sedangkan pada radar adalah pembukaan antena. Sebuah aperture sintetis dibangun oleh pergerakan antena pada berbagai posisi di sepanjang jalur penerbangan. Pada radar, antena dipasang dibagian bawah pesawat dan diarahkan kesamping, sistem ini dikenal sebagai Side Looking Aperture Radar (SLAR). SAR merupakan teknik yang digunakan untuk menghasilkan radar image, dan menyediakan kemampuan yang unik sebagai alat pencitraan. SAR dapat menghasilkan penerangan sendiri (pulsa radar) tidak bergantung pada penerangan matahari, sehingga dapat melakukan peliputan baik di siang hari maupun di malam hari. Dan, karena panjang gelombang radar lebih besar dari sinar tampak
8
maupun Infra merah, SAR dapat menembus awan maupun debu dimana kondisi ini tidak memungkinkan untuk sistem optik (NASA 1996). Sinyal radar dapat disaring sedemikian rupa sehingga getaran gelombang elektrik dibatasi hanya pada bidang datar yang tegak lurus arah perjalanan gelombang.
Gambar 3 Wahana dan arah tembak sensor (NASA 1996). Satu sinyal radar dapat ditransmisikan pada bidang horizontal (H) ataupun vertikal (V), demikian pula dapat diterima pada bidang mendatar maupun tegak sehingga ada empat kombinasi sinyal transmisi dan penerimaan yang berbeda, yaitu dikirim H diterima H (HH), dikirim H diterima V (HV), dikirim V diterima H (VH), dan dikirim V diterima V (VV). Karena berbagai objek mengubah polarisasi tenaga yang dipantulkan dalam berbagai tingkatan maka bentuk polarisasi sinyal mempengaruhi kenampakan objek pada citra yang dihasilkan. Banyak sifat khas medan yang bekerja bersama panjang gelombang dan polarisasi sinyal radar untuk menentukan intensitas hasil balik radar dari objek. Akan tetapi faktor utama yang mempengaruhi intensitas hasil balik sinyal objek adalah ukuran (geometris) dan sifat khas elektrik objek. Efek geometri sensor/objek dari intensitas backscatter radar terpadu dengan efek kekasaran permukaan. Permukaan yang kasar bertindak sebagai pemantul baur dan memencar tenaga datang ke semua arah dan hanya mengembalikan sebagian kecil ke antena. Suatu permukaan halus pada umumnya memantulkan sebagian besar
9
tenaga menjauhi sensor dan mengakibatkan sinyal hasil balik yang rendah. Meskipun demikian orientasi objek terhadap sensor harus dipikirkan juga karena permukaan halus yang mengarah ke sensor akan menghasilkan sinyal balik yang sangat kuat (Lillesand & Kiefer 1979).
Gambar 4 Mekanisme hamburan balik pada radar di setiap jenis permukaan (NASA 1996). NASA (1996) mengategorikan nilai hamburan balik pada radar ke dalam beberapa kelas, yaitu nilai backscatter sangat tinggi (berkisar -5dB ke atas) biasanya terjadi pada objek lereng menghadap sensor, incident angle kecil, permukaan objek yang sangat kasar, hutan yang tergenang, dan objek buatan. Pada kelas nilai backscatter tinggi (berkisar 0 sampai -10dB) bisanya terjadi pada objek dengan permukaan yang kasar dan vegetasi rapat. Hamburan balik pada radar merupakan ukuran kuantitatif dari intensitas energi yang kembali ke antena. Nilai hamburan balik yang dihasilkan pada sebuah sensor radar dipengaruhi beberapa faktor antara lain kedalaman penetrasi dari gelombang radar, kekasaran permukaan objek dan sifat-sifat dielektrik volume objek. Michigan
Microwave
Canopy
Scattering
Model
(MIMICS)
telah
dikembangkan untuk memberikan pemahaman terhadap hamburan balik (backscatter) radar pada vegetasi. Beberapa bentuk hamburan yang dapat dikalkulasi adalah hamburan pada permukaan dan volume tajuk, hamburan langsung pada permukaan tanah, hamburan langsung pada batang, hamburan dari permukaan tanah ke batang, dan hamburan dari permukaan tanah ke tajuk (Dobson et al. 1990). Mekanisma hamburan balik ini digambarkan sebagai berikut.
10
Gambar 5 Mekanisme hamburan balik pada tegakan hutan (Dobson et al. 1990). Gelombang radar yang lebih panjang menghasilkan nilai backscatter yang tinggi pada penetrasi batang, percabangan, permukaan tanah dan tajuk. Sedangkan gelombang yang lebih pendek menghasilkan nilai backscatter yang tinggi hanya pada tajuk saja. Kemampuan gelombang panjang untuk mempenetrasikan kanopi hutan dengan lebih baik menjadi dasar kemampuan dari sistem SAR untuk secara langsung menduga kuantiti dari struktur tegakan dalam hal ini yang berkaitan dengan biomassa dimana sebagian besar biomassa berada pada batang dan percabangan (ranting-ranting besar). Hamburan balik yang kuat dari vegetasi akan dihasilkan oleh tipe vegetasi rapat. Sistem radar L-band bekerja pada gelombang maksimum untuk citra radar yang tersedia. L-band memiliki kemampuan besar untuk menembus daun-daunan hingga ke pokok batang yang paling bawah. Banyak studi yang telah dilakukan dan menemukan hubungan yang kuat antara biomassa dan hamburan balik pada SAR (Mitchard et al. 2009; Sarker & Nichol 2010). Studi-studi tersebut selain menemukan hubungan yang kuat antara biomassa dan backscatter juga menemukan bahwa backscatter SAR meningkat seiring peningkatan biomassa sampai mencapai nilai saturasi tertentu yang mana nilai pendugaan tertinggi akan berada pada frekuensi yang rendah atau memiliki nilai backscatter rendah (Bergen & Dobson 1999). 2.2.1 Citra Satelit ALOS PALSAR ALOS (Advanced land Observing Sattelite) merupakan satelit yang diluncurkan oleh Badan Luar Angkasa Jepang pada bulan Januari 2006. Satelit ALOS ini membawa tiga jenis sensor yaitu PALSAR (Phased Array L-band
11
Synthetic Aperture Radar), PRISM (Panchromatic Remote-sensing Instrument for Stereo Mapping), dan AVNIR-2 (Advanced Visible and Near Infrared Radiometer type-2). PALSAR merupakan sensor gelombang mikro aktif yang bekerja pada frekuensi band L. Sensor PALSAR mempunyai kemampuan untuk menembus awan, sehingga informasi permukaan bumi dapat diperoleh setiap saat, baik malam ataupun siang hari. Data PALSAR ini dapat digunakan untuk pembuatan DEM, Interferometri untuk mendapatkan pergeseran tanah, maupun kandungan biomassa, monitoring kehutanan, pertanian, tumpahan minyak (oil spill), soil moisture, mineral, dan lain-lain (Rosenqvist et al. 2004). Untuk dapat bekerja dengan ketiga instrumen di atas, ALOS dilengkapi dengan dua teknologi yang lebih maju. Pertama teknologi yang mampu mengerjakan data dalam kapasitas yang sangat besar dengan kecepatan tinggi, dan selanjutnya kapasitas untuk menentukan posisi satelit dengan ketinggian yang lebih tepat. Keterangan umum tentang ALOS disajikan pada Tabel 2. Tabel 2 Keterangan umum ALOS Alat Peluncuran Tempat Peluncuran Berat Satelit Power Waktu Operasional Orbit
Roket H-IIA Pusat Ruang Angkasa Tanagashima 4000 Kg 7000 W 3-5 Tahun Sun-Synchronous Sub-Recurr Orbit Recurrent Period: 46 Hari Sub Cycle 2 hari Tinggi Lintasan: 692 km diatas Ekuator Inklinasi: 98,2°
Sumber: Jaxa 2010
Secara ringkas terdapat lima misi dari satelit ALOS (JAXA 2010), yaitu: 1.
Kartografi
: Untuk menyediakan peta wilayah Jepang dan wilayah Asia Pasifik
2.
Pemantauan Regional
: Melakukan pemantauan regional untuk pengembangan berkelanjutan kesediaan
pembangunan dan sumber
harmonisasi daya
pengembangan pembangunan 3.
Monitoring Bencana
: Melakukan monitoring bencana alam
4.
Survei Sumberdaya
: Untuk survei sumber daya alam
yang antara alam
12
5.
Pengembangan Teknologi
: Mengembangkan teknologi penginderaan jauh yang tepat untuk masa sekarang dan akan datang.
PALSAR merupakan salah satu instrumen ALOS dengan sensor aktif untuk pengamatan cuaca dan permukaan daratan pada siang dan malam hari dengan sistem yang lebih maju dari JERS-1 SAR. Sensor PALSAR mempunyai sorotan yang dapat disetir dalam elevasi, disamping mode ScanSAR. Bentuk dari instrumen PALSAR dan prinsip pengambilan obyeknya disajikan pada Gambar 6 dan Gambar 7. Sedangkan karakter teknik sensor PALSAR disajikan pada Tabel 3.
Gambar 6 Instrumen PALSAR (JAXA 2010).
Gambar 7 Prinsip geometri dari PALSAR (JAXA 2010).
13
Sensor ini merupakan sensor gelombang mikro aktif yang dapat melakukan observasi siang dan malam tanpa terpengaruh pada kondisi cuaca. Melalui salah satu observasinya, yaitu ScanSAR sensor ini memungkinkan untuk melakukan pengamatan permukaan bumi dengan cakupan area yang cukup luas yaitu 250 hingga 350 km. ScanSAR mempunyai kemudi berkas cahaya (yang dapat diatur) pada elevasi (ketinggian) dan didesain untuk memperoleh cakupan yang lebih lebar daripada SAR konvensional. Bentuk dari instrument PALSAR dan prinsip pengambilan objeknya disajikan pada Tabel 3. Tabel 3 Karakteristik PALSAR Mode Fine
ScanSAR
Frekuensi Lebar Kanal Polarisasi
1.270 MHz (L-Band) 28/114 MHz HH atau VV
HH+HV+VH+VV
100 m (multi look) 250-350 km 18-43 derajat <-25 dB
30 m 30 km 8-30 derajat <-29 dB
5 bit AZ:8.9 m x EL :2.9 m
3 bit atau 5 bit
Resolusi Spasial Lebar cakupan Incidence Angle NE Sigma 0 Panjang Bit Ukuran
HH/VV/HH+HV VV+VH
atau
10 m (2 look)/20 m (4 look) 70 km 8-60 derajat <-23 dB (70 km) <-25 dB (60 km) 3 bit atau 5 bit
Polarimetric (Experiment Mode)
Sumber : Jaxa (2010)
2.3 Sistem Informasi Geografi Johnson (1992) mendefinisikan SIG sebagai sebuah sistem yang berguna untuk menangani dan menganalisis data geografi untuk banyak pemakai dan aplikasi. Data yang digunakan dalam SIG adalah data geografis yang terdiri dari data geometrik dan data deskriptif. Data geometrik berhubungan dengan lokasi, bentuk dan hubungan antar kenampakan, misal peta-peta atau data dari penginderaan jauh. Sementara itu, data deskriptif berhubungan dengan sifat-sifat dari kenampakan, misal tabel, grafis dan keterangan lainnya. Data-data tersebut dipakai sebagai visualisasi dan menerangkan keadaan dunia yang sesungguhnya. SIG adalah informasi yang dibuat untuk berbagai data yang dikumpulkan dengan keruangan atau koordinat geografi. Dengan kata lain, SIG adalah sebuah sistem database dengan kemampuan spesifik untuk data keruangan dan juga sebuah perangkat operasi untuk bekerja dengan data.
14
Menurut Paryono (1994) SIG memerlukan data masukan agar berfungsi dan memberikan informasi hasil analisisnya. Data masukan tersebut dapat diperoleh dari tiga sumber, yaitu : (a) Data lapangan, data ini diperoleh langsung dari pengukuran lapangan secara langsung, seperti misalnya pH tanah, salinitas air, curah hujan, jenis tanah, dan sebagainya (b) Data peta, informasi yang lebih terekam pada peta kertas atau film, dikonversikan ke dalam bentuk digital. Misalnya, peta geologi, peta tanah dan sebagainya. Apabila data sudah terekam dalam bentuk peta, tidak lagi diperlukan data lapangan, kecuali untuk mengecek kebenarannya. (c) Data citra pengideraan jauh, citra penginderaan jauh yang berupa foto udara atau radar dapat diinterpretasi terlebih dahulu sebelum dikonversi ke dalam bentuk digital. Sementara itu, citra yang diperoleh dari satelit yang sudah dalam bentuk digital dapat langsung digunakan setelah
diadakan
koreksi seperlunya. Lebih lanjut dikatakan bahwa ketiga sumber tersebut saling mendukung satu terhadap yang lain. Data lapangan dapat digunakan untuk membuat peta fisis, sedangkan data penginderaan jauh juga memerlukan data lapangan untuk lebih memastikan kebenaran data tersebut. Jadi ketiga sumber data saling berkaitan, melengkapi dan mendukung, sehingga tidak boleh ada yang diabaikan. Menurut Jaya (2002) pada bidang kehutanan, SIG sangat diperlukan guna mendukung pengambilan keputusan untuk memecahkan masalah keruangan (spasial) mulai dari tahap perencanaan, pengelolaan sampai dengan pengawasan. SIG sangat membantu memecahkan permasalahan yang menyangkut luasan (polygon), batas (line atau Arc) dan lokasi (point). Data spasial (peta) yang umum digunakan di bidang kehutanan antara lain adalah: 1. Peta Rencana Tata Ruang, 2. Peta Rencana Tata Guna Hutan, 3. Peta Rupa Bumi (Kontur), 4. Peta Jaringan Jalan, 5. Peta Jaringan Sungai, 6. Peta Tata Batas, 7. Peta Batas Unit Pengelolaan Hutan, 8. Peta Batas Administrasi Kehutanan,
15
9. Peta Tanah, 10. Peta Iklim, 11. Peta Geologi, 12. Peta Vegetasi, 13. Peta Potensi Sumberdaya Hutan. 2.4 Jati (Tectona grandis Linn F.) Jati (Tectona grandis Linn F.) hingga saat ini masih menjadi komoditas mewah yang banyak diminati masyarakat. Tanaman jati diklasifikasikan ke dalam famili verbenaceae, genus tectona, dengan nama species terbanyak di Indonesia adalah Tectona grandis Linn F, dimana jenis ini merupakan jenis terbaik dibandingkan dengn jenis Jati lainnya. Kayu jati memiliki berat jenis antara 0.62 – 0.75. Saat ini budidaya jati banyak dilakukan di daerah-daerah tropis, di Panama lahan tidur yang lama tidak dimanfaatkan kemudia dilakukan penanaman jati dan telah dilakukan reforestasi dari 11.000 hingga 34.600 ha. Jati merupakan kayu terkuat ke tiga di dunia, dan budidaya jati telah di lakukan di banyak daerah tropis untuk menghasilkan kayu kualitas tinggi. Pertumbuhan jati terbaik biasanya pada daerah-daerah rendah dan panas, di pulau Jawa terutama pada dataran rendah dan berbukit-bukit, kurang air, yang terdiri dari formasi tua kapur dan megalit. Tanaman jati juga tumbuh pada daerah yang memiliki musim kering 3-5 bulan kering, curah hujan rata-rata 1.250 – 2.500 mm/tahun dengan ketinggian kurang dari 700 mdpl dan temperatur rata-rata 22 26oC (Khaerudin 1994). Tajuk pepohonan dalam hutan jati akan menyerap dan menguraikan zat-zat pencemar (polutan) dan cahaya berlebihan. Seperti halnya daun lain, daun jati juga melakukan fotosintesis yang menyerap karbondioksida dari udara dan kemudian melepaskannya kembali sebagai oksigen dan uap air. Hal ini membantu menjaga kestabilan iklim di dalam dan di sekitar hutan. Sejak abad ke-9 tanaman jati yang merupakan tanaman tropika dan subtropika telah dikenal sebagai pohon yang memiliki kayu kualitas tinggi dan bernilai jual tinggi. Jati digolongkan sebagai kayu mewah dan memiliki kelas awet tinggi yang tahan terhadap gangguan rayap serta jamur dan mampu bertahan sampai 500 tahun.
16
Secara alamiah tanaman jati banyak ditemukan di negara-negara Asia Tenggara dan Asia Selatan. Pada abad ke-19 jati juga mulai dibudidayakan di Amerika tropik. Keberhasilan permudaan sejak akhir abad 19 telah dapat mengembangkan luas kawasan hutan jati di pulau Jawa. Pada akhir abad 19 luas hutan jati di pulau jawa seluruhnya diperkirakan berkisar 650.000 ha. Dan hingga tahun 1985 luas hutan jati terus bertambah hingga mencapai 1.069.712 ha. Hutan jati yang sebagian besar berada di Pulau Jawa dikelola oleh PT. Perhutani. Dimana luas daerah pengelolaannya mencapai 2.6 juta ha, terdiri dari 54 KPH (Kesatuan Pemangkuan Hutan). Produksi hutan jati yang dikelola oleh Perhutani rata-rata 800.000 m3/tahun, dimana sebagian besar produksi hutan jati dijual dalam bentuk log. 2.5 Biomassa Brown (1997) mendefinisikan biomassa hutan sebagai bobot total materi organisma hidup setiap pohon di atas permukaan tanah dinyatakan dalam bobot kering ton per unit area. Biomassa dapat pula didefinisikan sebagai bobot dari material tumbuhan hidup per unit area. Total biomassa yang terdapat pada areal hutan dibagi ke dalam dua bagian yaitu biomassa di atas permukaan dan biomassa di bawah permukaan tanah. Biomassa pada dasarnya terdiri dari bobot organisma hidup di atas permukaan dan di bawah permukaan tanah, seperti pohon, semak belukar, tumbuhan menjalar, akar dan berat organisma mati dan sampah kasar yang terasosiasi dengan tanah. Karena terdapat kesulitan pada pengumpulan data lapangan biomassa di bawah permukaan (Below-Ground Biomass, BGB), penelitian estimasi biomassa yang telah banyak dilakukan sebelumnya terfokus pada biomassa di atas permukaan (Above-Ground Biomass, AGB) (Lu 2006). Tabel 4 menyajikan rangkuman dari beberapa teknik pendugaan biomassa yang berbeda berdasarkan (1) pengukuran lapangan, (2) remote sensing, dan (3) GIS. Biomassa tegakan hutan dipengaruhi oleh umur tegakan hutan, sejarah perkembangan vegetasi, komposisi dan struktur tegakan (Lugo & Snedaker 1974). Pengukuran biomassa pada dasarnya mengacu pada empat teknik pengukuran (Lu 2006): (a) teknik pemetaan pemanenan atau teknik sampling destruktif (b) teknik sampling non-destruktif (c) pengukuran berdasarkan data remote sensing yang
17
dihasilkan oleh sistem airborne/spaceborne, and (d) estimasi menggunakan model. Tabel 4 Beberapa metode untuk menduga biomassa Ketegori Metode Data Karakteristik digunakan Metode DasarPengukuran Lapangan
Metode DasarPenginderaan Jauh
Metode DasarSIG
Sumber: Lu 2005
Referensi
Penarikan contoh Destruktif
Pohon contoh
Individu pohon
Klinge et al. (1975)
Allometric equations (persamaan allometric)
Pohon contoh
Individu pohon
Konversi dari volume ke biomassa
Volume dari pohon contoh atau tegakan
Individu pohon atau tegakan vegetasi
Metode berdasarkan fine spatialresolution data Metode berdasarkan medium spatialresolution data Metode berdasarkan coarse spatialresolution data Metode berdasarkan data Radar
Aerial photographs, IKONOS
Per-pixel level
Overman et al.(1994), Honzak et al. (1996), Nelson et al. (1999) Brown dan Lugo (1984), Brown et al. (1989), Brown dan Lugo (1992), Gillespie et al. (1992), Segura dan Kanninen (2005) Tiwari dan Singh (1984), Thenkabail et al. (2004)
Landsat TM/ETM+, SPOT
Per-pixel level
IRS-1C WiFS, AVHRR
Per-pixel level
Radar, lidar
Per-pixel level
Harrel et al. (1997), Lefsky et al. (1999b), Santos et al. (2002, 2003)
Metode berdasarkan Ancillary data
Elevasi, kemiringan, tanah, presipitasi, dll.
Per-pixel level atau per-field level
Brown et al. (1994), Iverson et al. (1994), Brown dan Gaston (1995)
Roy dan Ravan (1996), Nelson et al. (2000a), Steininger (2000), Foody et al. (2003), Lu (2005) Barbosa et al. (1999), Wylie et al. (2002), Dong et al. (2003)
18
Metode Pendugaan Biomassa dengan Non-Destruktif Sampling Pendekatan destruktif untuk menduga biomassa memberikan hasil yang paling akurat, tetapi penerapan teknik ini tidak dapat dilakukan pada seluruh areal hutan karena kerusakan yang diakibatkan cukup besar. Selain kerusakan yang diakibatkan,
mahalnya
biaya
dan
banyaknya
waktu
yang
dibutuhkan
dibandingkan dengan teknik pendugaan biomassa yang lain menjadi bahan pertimbangan dalam penggunaan teknik ini. Teknik pendugaan biomasa yang banyak dilakukan pada saat ini adalah pendekatan non-destruktif sampling dimana pada pendekatan ini tidak diperlukan pemanenan pohon. Pendekatan non-destruktif sampling memilik persamaan regresi yang berbeda-beda, dengan parameter penyusunya seperti tinggi pohon, diameter pohon, volume batang, dan basal area untuk menduga biomassa. Berikut merupakan persamaan regresi yang banyak digunakan untuk menduga biomassa menurut Husch et al. (2003);
B = C0 + C 1 S B = C0 + C1S + C2S2 B = exp (C1S) B = exp (C2SC1) B = C0 + C1logS Menurut Chave et al. (2005) berdasarkan pengujian yang telah dilakukan didapatkan bahwa diameter pohon merupakan salah satu variabel yang penting bagi pendugaan biomassa selain berat jenis pohon dan tipe hutan. Teknik estimasi biomassa non-destruktif memiliki efisiensi yang baik jika dibandingkan dengan teknik sampling destruktif. Pada kasus area hutan tropis campuran yang heterogen, survey lapangan tidak mungkin dilakukan untuk mengambil sample pada plot dengan aksesibilitas yang rendah. Untuk memonitoring area hutan tropis campuran yang heterogen maupun hutan tropis homogen, penginderaan jauh menyediakan alat yang paling sesuai dan efektivitas waktu serta biaya jauh lebih baik dibandingkan pengukuran in-situ. Beberapa tahun terakhir, teknik penginderaan jauh (remote sensing) telah menjadi hal umum dalam menduga AGB (Lu 2005). Penginderaan jauh sistem
19
optik pada dasarnya merespon pada struktur kimia daun seperti Normilize Difference Vegetation Index (NDVI) (Dong et al. 2003). Berdasarkan resolusi spasial, Lu (2006) mengategorikan data penginderaan jauh (citra satelit) untuk estimasi AGB ke dalam tiga kategori yaitu fine spatial-resolution data (resolusi spasial kurang dari 5 m), medium spatial-resolution data (resolusi spasial pada kisaran antara 10 m hingga 100 m), dan coarse spatial-resolution data (resolusi spasial lebih dari 100 m). Penggunaan coarse spatial-resolution data memiliki keterbatasan dikarenakan ukuran pixelnya yang besar, dimana berbagai jenis pohon dari bermacam-macam area hutan terdapat di dalam pixel tersebut dan perbedaan yang amat besar antara ukuran pixel dengan plot yang dibuat untuk pengukuran lapangan. Pendugaan Biomassa Hutan Jati Berdasarkan Model Alometrik Pada pendugaan nilai biomassa tegakan jati di lokasi peneltian digunakan model alometrik biomassa terukur yang telah dikembangkan oleh Hendri (2001) di daerah KPH Cepu yang merupakan daerah Hutan Tanaman Jati yang berseberangan langsung dengan lokasi peneltian. Dimana R2 pada persamaan alometik dengan destructive sampling ini sebesar 94,1%. Selain itu kondisi umum lapangan baik kondisi tegakan maupun kondisi iklim daerah KPH Cepu memiliki kesamaan dengan kondisi umum lapangan KPH Kebonharjo. Berikut merupakan persamaan alometrik yang dikembangkan oleh Hendri (2001) (Hendri dalam Tiryana 2011). W = 0.2759D2.2227 (R2 = 0.941)
Pendugaan Biomassa Hutan Jati Berdasarkan Nilai Biomass Expansion Factor (BEF) Biomass Expansion Factor (BEF) didefinisikan sebagai rasio total bobot biomassa kering tanur di atas permukaan tanah pada diameter minimum (dbh) 10 cm atau lebih dengan bobot biomassa kering tanur pada volume yang diinventarisasi atau rasio antara AGB total dengan biomassa batang yang dapat dimanfaatkan. BEF pada setiap tegakan hutan di berbagai negara sangat bervariasi,
penghitungan
BEF
yang
banyak
dilakukan
adalah
dengan
menggunakan dua metode, yaitu perhitungan BEF dari diameter pohon secara
20
langsung dan BEF yang dihitung dengan menggunakan Double Sampling. Expansion Factor berdasarkan diameter saja secara konsisten memiliki nilai yang lebih baik dari Expansion Factor yang dicari berdasarkan Double sampling. Hipotesis yang diuji sebelumnya menunjukkan bahwa perbedaan nilai Expansion Factor pada daerah-daerah yang berbeda disebabkan oleh perbedaan rata-rata ukuran pohon yang diinventarisasi. Biomass Expansion Factor (BEF) pada tegakan jati yang dikembangkan di daerah tropis Panama, di hitung dengan membagi total prororsi biomassa dengan biomassa batang sehingga menghasilkan nilai BEF sebesar 1,53186 (Kraenzel et al. 2003).
21
BAB III METODOLOGI 3.1 Lokasi dan Waktu Penelitian Penelitian ini dilakukan di KPH Kebonharjo Perum Perhutani Unit I, Jawa Tengah. Meliputi Bagian Hutan (BH) Tuder dan Balo, pada Kelas Perusahaan Jati. Penelitian di lakukan pada bulan November 2010 sampai dengan Maret 2011. Lokasi penelitian disajikan pada Gambar 8.
Gambar 8 Lokasi penelitian. 3.2 Bahan dan Alat Alat yang digunakan pada penelitian ini berupa seperangkat komputer dengan kelengkapan aplikasi penunjang pengolahan data seperti MS Office, Erdas Imagine, ArcView GIS Ver 3.3, SPSS Ver 17. Sedangkan alat yang menunjang pengambilan data di lapangan antara lain GPS, Phi-Band, galah bantu, meteran, kamera dijital dan tallysheet. Data yang digunakan dalam penelitian terdiri dari: 1. Citra ALOS PALSAR resolusi spasial 50 m tahun perekaman/peliputan 2009 dan citra ALOS PALSAR resolusi spasial 12.5 m tahun
22
perekaman/peliputan 2009 daerah Kebonharjo, Kabupaten Rembang, Jawa Tengah. 2. Data hasil pengukuran dimensi tegakan (Diameter setinggi dada/D dan Tinggi total/Tt) keadaan tapak, jarak antar pohon dan azimut per plot pengamatan. 3. Data informasi petak dan anak petak KPH Kebonharjo. 3.3 Tahapan Pelaksanaan Tahapan dalam kegiatan penelitian ini secara umum terdiri dari persiapan dan pengambilan data lapangan, pengolahan data lapangan, pengolahan data citra, pembuatan model penduga biomassa, dan pemetaan biomassa (Gambar 10). 3.3.1 Persiapan dan Pengambilan Data Lapangan Pengambilan data lapangan direncanakan di atas peta kerja dan peta administrasi KPH Kebonharjo, Perhutani Unit II Jawa Tengah. Pemilihan titik plot pengukuran lapangan dilakukan berdasarkan sebaran kelas umur di lokasi penelitian dan kenampakan Citra ALOS PALSAR. Penentuan titik dilakukan secara purposive sampling. Titik pengambilan sampel di tentukan sebanyak 64 titik yang menyebar di seluruh lokasi penelitian, berada di dua bagian hutan yaitu Bagian Hutan Balo dan Bagian Hutan Tuder. 3.3.2 Bentuk dan Ukuran Plot Contoh Pengambilan data lapangan dilakukan di hutan tanaman jati, KPH Kebonharjo Perhutani Unit II Jawa Tengah. Plot contoh yang digunakan berbentuk lingkaran dengan luasan plot disesuaikan pada titik dengan KU (kelas umur) tertentu yang akan dikunjungi di lapangan.
23
Persiapan dan pengumpulan data
Mulai
Hasil Inventarisasi Tegakan & Dimensi Tegakan
Citra ALOS PALSAR 2009 Citra Backscatter Ekstraksi Nilai Dijital Setiap plot
Perhitungan Biomassa
Biomassa Alometrik Hendri
Biomassa BEF
Konversi Nilai Dijital ke Backscatter
Nilai Backscatter
Analisis Statistik & Penyusunan Model
Peta Petak KPH Kebonharjo
Verifikasi Model & Pemilihan Model Terbaik
Model Terbaik Pembuatan Peta Sebaran Biomassa
Penghitungan Akurasi Kelas Biomassa
Selesai Gambar 9 Diagram Alir Penelitian
Peta Sebaran Biomassa
24
Pada KU I – II dilakukan pengukuran dengan plot lingkaran seluas 0.02 ha, pada KU III – IV dilakukan pengukuran dengan plot lingkaran seluas 0.04 ha, sedangkan untuk KU ≥ V dilakukan pengukuran dengan plot lingkaran seluas 0.1 ha. Unit contoh yang digunakan merupakan hasil klasifikasi visual pada citra ALOS PALSAR resolusi spasial 50 m dan 12,5 m berada dalam areal kerja KPH Kebonharjo, Bagian Hutan Tuder dan Bagian Hutan Balo, Perhutani Unit I Jawa Tengah.
Gambar 10 Sketsa plot pengukuran. 3.3.3 Pengambilan Data Lapangan 3.3.3.1 Pemilihan Titik Pengukuran Lapangan Pemilihan titik pengukuran lapangan dilakukan berdasarkan peta petak dan anak petak KPH Kebonharjo Perhutani Unit I Jawa Tengah. Pemilihan titik dilakukan dengan metode purposive sampling, dimana kriteria pengambilan titik didasarkan pada sebaran kelas umur pada peta kerja KPH Kebonharjo Perhutani Unit I Jawa Tengah. Pengukuran koordinat titik pengamatan dapat diukur dengan GPS atau menggunakan koordinat peta yang ada. Kemudian, dilakukan perekaman posisi area contoh menggunakan GPS. 3.3.3.2 Pengukuran Parameter Tegakan Parameter tegakan yang diukur berupa diameter pohon setinggi dada (D), tinggi total (Tt), tinggi bebas cabang (Tbc) dan jenis pohon. Tingkat vegetasi yang diukur berdasarkan kelas umur (KU).
25
3.3.4
Pengolahan Data
3.3.4.1 Model Pendugaan Biomassa Lapangan. Model alometrik pendugaan biomassa yang digunakan pada penelitian ini adalah model alometrik biomassa Hendri (2001) yang diformulasikan kembali oleh Tiryana (2011), dimana parameter penduga biomasa (B) yang diukur adalah diameter setinggi dada (D). Model tersebut secara matematis ditentukan dengan persamaan sebagai berikut :
Biomassa yang diukur dalam penelitian ini merupakan biomassa di atas permukaan tanah (above-ground biomass) dari tegakan jati pada tingkat usia tertentu (KU). Selain menggunakan persamaan alometrik pada daerah penelitian, nilai biomassa dihitung dengan menggunakan BEF (Biomass Expansion Factor), dimana BEF didefinisikan sebagai rasio total bobot biomassa kering tanur diatas permukaan tanah pada diameter (diameter setinggi dada/D) minimum 10 cm atau lebih dengan bobot biomassa kering tanur pada volume yang dapat dimanfaatkan/BEF pada biomassa kering tanur pada volume batang. BEF merupakan faktor koreksi bagi nilai biomassa yang diduga melalui volume (BVakt, ton/ha). Perhitungan biomassa di atas permukaan dengan menggunakan BEF dilakukan dengan menggunakan rumus : Volume (V, m3) dihitung dengan menggunakan formulasi :
Bbef = Biomassa diduga menggunakan BEF (ton/ha) Vb = Volume untuk Bagian Hutan Balo (m3) Vt = Volume untuk Bagian Hutan Tuder (m3) ρ = Berat jenis rata-rata pohon jati sebesar 0.67 ton/m3 BEF = Biomass Expansion Factor dengan nilai koefisien 1,53186 untuk Jati pada hutan tropis (Kraenzel et al. 2003).
26
3.3.4.2 Ekstraksi Nilai Dijital pada Citra ALOS PALSAR Dengan menggunakan ekstensi square buffer pada ArcView 3.3 dibuat buffer pada titik pengamatan dilapangan dengan ukuran buffer 5 piksel × 5 piksel atau setara dengan 250 m × 250 m yang ditentukan berdasarkan pertimbangan error GPS dan pergeseran citra. Square buffer yang dihasilkan kemudian digunakan sebagai AOI (Area of Interest) sehingga didapat nilai dijital rata-rata pada buffer titik pengamatan. Nilai dijital yang dihasilkan kemudian di konversi menjadi nilai hamburan balik (backscatter) dengan menggunakan formulasi sebagai berikut (Shimada et al. 2009) :
Dimana : BS = Backscatter (dB) dN = Nilai dijital (degree) CF = Calibration factor dari Citra ALOS PALSAR peliputan tahun 2009 sebesar -83 (JAXA Publication) 3.3.4.3 Penyusunan dan Pemilihan Model Penyusunan model hubungan antara biomassa di atas permukaan tanah dengan nilai backscatter pada citra ALOS PALSAR menggunakan beberapa model matematika sebagai berikut : Tabel 5 Model yang digunakan untuk pendugaan biomassa Jenis Model Bentuk Model Model linier Model eksponensial Inverse polynomial Schumacher
Y = a + bX Y = a+(1/b)X Y = Exp(a + bX) Y = a*bX Y = X/(a+bX) Y = a*(Exp(b/X))
Keterangan : Y = Biomassa; X = Backscatter ALOS PALSAR
Pemilihan model dilakukan dengan memperhatikan koefisien determinasi terkoreksi R2adj dan Root Mean Square Error (RMSE) yang dihasilkan oleh masing-masing persamaan. Koefisien determinasi terkoreksi adalah koefisien determinasi yang telah terkoreksi dari derajat bebas sisa dan derajat bebas totalnya. Dimana koefisien determinasi terkoreksi dihitung menggunakan formulasi sebagai berikut:
27
Dimana : JKS = Jumlah kuadrat sisa JKT = Jumlah kuadrat total (n - p) = Derajat bebas sisa (n - 1) = Derajat bebas total Sedangkan akar kuadrat eror dihitung berdasarkan formula : MSE = RMSE = Dimana : MSE = Kuadrat tengah sisa RMSE = Akar kuadrat tengah sisa yi = Biomassa ke-i = Rata-rata biomassa ke-i n = Jumlah plot sampel p = Jumlah parameter yang digunakan 3.3.4.4 Verifikasi Model Verifikasi model dilakukan dengan membandingkan antara hasil kandungan biomassa di atas permukaan tanah dengan menggunakan model terpilih dan hasil pengukuran di lapangan menggunakan persamaan alometrik yang diasumsikan sebagai biomassa aktual. Verifikasi model dilakukan dengan cara mengambil piksel secara purposive pada cira sebanyak 26 plot kemudian dilakukan pengukuran di lapangan sesuai dengan prosedur yang dilakukan pada penyusunan model. Untuk membandingkan hasil pendugaan biomassa pada model terbaik yang dihasilkan dalam penelitian ini dengan hasil pengukuran biomassa di lapangan menggunakan alometrik digunakan uji t-student berpasangan (Mattjik & Sumertajaya 2000).
28
Dengan menggunakan hipotesis uji sebagai berikut : H0 : µ1 - µ2 = 0 (Biomassa aktual = biomassa model) H1 : µ1 - µ2 ≠ 0 (Biomassa aktual ≠ biomassa model) Model yang dianggap mewakili data dan layak digunakan didasarkan pada thitung dengan kriteria apabila thitung
0,05 (taraf nyata 5%), maka model pendugaannya layak digunakan dan sebaliknya jika thitung > t(α/2) atau nilai signifikansi <0,05 (taraf nyata 5%), maka model penduganya kurang layak digunakan. 3.3.5 Pembuatan Citra Pendugaan Biomassa dan Peta Kelas Biomassa Citra pendugaan biomassa diturunkan dari citra backscatter yang diproses dengan menggunakan Modeler (
) pada perangkat lunak ERDAS Imagine 9.1
dengan menggunakan model pendugaan terpilih yang telah terverifikasi. Dari citra pendugaan biomassa dibuat peta sebaran kelas biomassa dengan menggunakan perangkat lunak ArcView 3.3, kemudian untuk mengetahui tingkat keterwakilan dan akurasi pembuatan peta sebaran dilakukan penghitungan kappa accuracy (K) dan overall accuracy (OA) pada peta sebaran kelas biomassa yang telah dibuat.
Dimana : Xii
= nilai diagonal dari matrik kontingensi baris ke-i dan kolom ke-i
Xi+
= jumlah piksel dalam kolom ke-i
X+i
= jumlah piksel dalam baris ke-i
N
= banyaknya titik contoh
29
Dalam matrik kontingensi dapat pula dihitung besarnya akurasi pembuat (Producer Accuracy/PA) dan akurasi pengguna (User Accuracy/UA) dari setiap kelas.
31
BAB IV KONDISI UMUM LOKASI PENELITIAN
4.1 Risalah Umum KPH Kebonharjo 4.1.1
Letak Secara geografis terletak pada 111°20′ - 111°30′ BT dan 6°30′ - 6°60′ LS.
KPH Kebonharjo memiliki luas wilayah ±17.801,30 ha. Berdasarkan pembagian wilayah administratif
KPH Kebonharjo berada pada wilayah Propinsi Jawa
Tengah dan Jawa Timur. Meliputi wilayah Kabupaten Rembang Jawa Tengah dengan luas ± 11.946,7 ha yang terdiri dari 9 kecamatan yaitu Kecamatan Gunem, Sale, Sarang, Sedan, Pamotan, Pancur, Kragan, Sluke, dan Kecamatan Lasem (terdapat sekitar 42 desa sekitar hutan). Wilayah Kabupaten Blora Jawa Tengah, seluas ± 2.140,4 ha yang terdiri dari 2 kecamatan yaitu Kecamatan Bogorejo dan Kecamatan Jiken (terdapat sekitar 7 desa sekitar hutan). Wilayah pemerintahan Kabupaten Tuban Jawa Timur, seluas ± 3.714,2 ha yang terdiri dari 2 kecamatan yaitu Kecamatan Kenduruan dan Kecamatan Jatirogo (terdapat sekitar 11 desa sekitar hutan). Dalam pembagian wilayah kerja Perum Perhutani, KPH Kebonharjo terletak di wilayah Perum Perhutani Unit I Jawa Tengah wilayah Seksi Perencanaan Hutan IV Rembang. 4.1.2
Batas Wilayah Batas wilayah KPH Kebonharjo meliputi :
1) Sebelah Utara : KPH Kebonharjo berbatasan langsung dengan Laut Jawa bermula di Kali Lasem, menyusuri pantai ke arah Timur sampai di Desa Bulu berpotongan tegak lurus dengan perpanjangan jalan Jatirogo – Bulu. 2) Sebelah Timur : Berbatasan dengan batas KPH Jatirogo bagian Barat, mulai dari pantai Laut Jawa di desa Bulu ke Selatan mengikuti jalan besar Bulu – Jatirogo melalui desa Ngepon / Ketangi sampai Pal DK No. 1 TPK Jatirogo ke Timur mengikuti jalan raya sampai dipertigaan jalan besar dari Blora. Dari sini terus ke Selatan
32
mengikuti jalan besar jurusan Blora sampai dukuh Bogel, terus menyusuri kali Kedung Gede sampai desa Sidomukti melintang dengan jalan DK lanjutan Alur D terus ke Selatan berpapasan dengan kali Gunung Wangon Kendik hingga dukuh Gunung Wangon Alur HV. 3) Sebelah Selatan : Berbatasan dengan wilayah KPH Cepu bagian Utara mulai dukuh Gunung Wangon dari Pal Batas Hutan B 509 (petak 163) Bagian Hutan Merah ke Barat mengikuti Alur HV hingga Pal Batas Hutan B 150 (petak 162) Bagian Hutan Merah terus mengikuti jalan DK sampai Pal Batas Hutan B 526/DK A1 (petak 169) mengikuti batas hutan sampai Pal Batas Hutan B 577 (petak 156) Bagian Hutan Merah terus mengikuti Alur HK sampai pertemuan dengan Alur GL, terus membelok ke Barat mengikuti Alur GL sampai Pal Batas Hutan B 578 (petak 132) Bagian Hutan Merah ( B 86 Bagian Hutan Nanas) ke Barat melalui desa Gempol Ketringan, hingga persimpangan jalan Ketringan dengan Ngaglik terus mengikuti Alur FO sampai Pal Batas Hutan B 1199 ( petak 122) Bagian Hutan Merah terus mengikuti jalan desa Gempol, mengikuti jalan desa sampai Bogorejo, selanjutnya mengikuti jalan besar Jatirogo – Blora sampai di pertigaan pasar desa Karang. 4) Sebelah Barat : Berbatasan dengan batas KPH Mantingan bagian Timur mulai dari pertigaan pasar desa Karang ke Utara melewati desa Poncosari, desa Juwet, desa Jurangjero sampai Batas Hutan B 1566 (petak 18) Bagian Hutan Merah, mengikuti Batas Hutan sampai di Pal B 10 (petak 16) Bagian Hutan Merah, belok ke Utara mengikuti Kecamatan Pamotan sampai Pal B 1008 (petak 47) Bagian Hutan Tuder, langsung mengikuti Batas Hutan sampai dengan Pal B 1050 dan terus ke Utara menyusuri batas Kecamatan Pamotan sampai desa Ukir terus mengikuti jalan jurusan desa Bamban. Dari belokan jalan desa Bamban menyeberang ke Utara sampai bertemu jalan jurusan dari desa Bamban ke desa Bedog (Bangunrejo) terus ke Utara sampai bertemu pada pertigaan jalan besar jurusan Kebonharjo – Lasem.
33
Dari pertigaan jalan terus ke Barat menyusuri jalan sampai dipertigaan Lasem antara jalan jurusan Bulu – Rembang, belok ke Barat sampai di sungai Lasem sampai Laut Jawa. Adapun letak kantor KPH Kebonharjo, berkedudukan di Kebonharjo, Kecamatan Jatirogo Kabupaten Tuban. Sedangkan KPH – KPH yang berbatasan langsung dengan KPH Kebonharjo :
4.1.3
1. Sebelah Timur
: KPH Jatirogo
2. Sebelah Selatan
: KPH Cepu
3. Sebelah Barat
: KPH Mantingan
Topografi Secara garis besar KPH Kebonharjo dapat digambarkan terletak pada
lereng Gunung Kendeng Utara dengan keadaan topografi landai di bagian utara (kecuali pada kelompok – kelompok Gunung Lasem), dimana semakin menuju arah selatan kelerengan semakin meningkat dan semakin bergelombang. Hal ini disebabkan adanya perbukitan – perbukitan dan lembah – lembah. Perbukitan disebelah utara adalah komplek Gunung Lasem dan disebelah selatan adalah Gunung Butak, sedangkan lembah – lembah terjadi karena adanya aliran Kali Kuning. 4.1.4
Enclave Enclave hanya terdapat di Bagian Hutan Tuder yang terletak di desa Karas
Kecamatan Sedan dengan luas ± 8,0 ha. Tanaman pertanian yang banyak ditemukan adalah kelapa, polowijo (palawija) dan persawahan (padi). Sebagian besar penduduknya bekerja sebagai petani yang mengandalkan pada musim (petani musim), terutama musim hujan. Pada musim kering masyarakat sulit bercocok tanam, dan untuk mencukupi kebutuhan hidupnya, mereka banyak yang mengembara ke daerah lain sebagai pekerja bangunan musiman. Sedang sebagian yang lain menetap di desa dengan sesekali mencari rencek di hutan.
BAB V HASIL DAN PEMBAHASAN 5.1 Rata-rata Biomassa di KPH Kebonharjo Berdasarkan Tabel 1 rata-rata biomassa yang ditemukan di lapangan dengan menggunakan model alometrik Hendri tidak lebih besar dari rata-rata biomassa lapangan yang dihitung berdasarkan Biomass Expansion Factor (BEF). Nilai biomassa rata-rata pada kelas umur II, III, V, VI, VIII, IX, dan XI pada biomassa Hendri memiliki nilai biomassa yang lebih rendah dari biomassa BEF. Nilai biomassa tertinggi baik pada perhitungan biomassa dengan menggunakan alometrik Hendri maupun dengan menggunakan BEF, yaitu pada kelas umur XI dengan nilai biomassa 328,60 Ton/ha dan 461,17 Ton/ha. Tabel 6 Rata-rata biomassa di KPH Kebonharjo per KU dan jumlah plot pengamatan per KU berdasarkan alometrik tegakan jati dan BEF Kelas Umur
Σ Plot Per KU
I II III IV V VI VII VIII IX XI
17 13 9 5 5 5 3 4 1 1 Total
63
Rata-rata Biomassa Per KU (Ton/ha) Alometrik Hendri BEF 42.72 35.88 111.76 124.31 118.42 126.22 184.55 155.91 125.78 274.43 140.54 144.99 191.74 176.44 211.22 245.83 199.07 288.63 328.69 461.17 1654.49 2033.87
5.2 Hasil Pengolahan Data Citra Nilai hamburan balik (backscatter) pada citra ALOS PALSAR didapat dengan mengekstraksi nilai dijital pada citra. Dilakukan pembuatan square buffer pada titik pengamatan ukuran 5 × 5 piksel setara dengan 250 m × 250 m ukuran sesungguhnya dilapangan pada citra ALOS PALSAR resolusi spasial 50 m dan square buffer ukuran 20 × 20 piksel setara dengan 250 m × 250 m ukuran sesungguhnya di lapangan pada citra ALOS PALSAR resolusi 12,5 m. Buffer titik ini dibuat dengan tujuan untuk mempermudah ekstraksi nilai dijital pada citra. Nilai dijital yang diperoleh merupakan rata-rata nilai dijital dari 25 (dua puluh lima) piksel terdekat pada citra ALOS PALSAR resolusi 50 m dan 400 (empat
36
ratus) piksel terdekat pada citra ALOS PALSAR resolusi 12,5 m dengan pusat plot, dimana nilai yang didapatkan adalah nilai tengah dari ekstraksi nilai dijital. Kemudian nilai dijital yang terekstraksi dikonversi menjadi nilai hamburan balik dengan menggunakan persamaan BS (dB) = 10*Log10(DN2) + CF (Shimada et al. 2009). Nilai hamburan balik (backscatter) pada masing-masing plot pengamatan disajikan pada Lampiran 3. Berdasarkan hasil perhitungan dapat dilihat bahwa nilai backscatter pada polarisasi HH memiliki nilai yang lebih besar dari backscatter pada polarisasi HV. Ukuran 5 × 5 piksel dan 20 × 20 piksel yang digunakan untuk mendapatkan nilai backscatter ini berfungsi sebagai pereduksi efek dari speckle dan error rektifikasi pada nilai koefisien backscatter (Austin et al. 2003). Selain itu penentuan jumlah piksel yang digunakan juga didasari oleh pertimbangan error GPS di lapangan pada saat ground check. 5.3 Penyusunan Model Pendugaan Biomassa 5.3.1 Biomassa Hendri dan Biomassa BEF Berdasarkan tabel pada Lampiran 1 hasil perhitungan biomassa dengan menggunakan model alometrik Hendri memiliki nilai biomassa terendah sebesar 7.52 ton/ha pada KU I dan nilai biomassa tertinggi sebesar 328.69 ton/ha pada KU XI. Berdasarkan tabel pada Lampiran 2 hasil perhitungan biomassa dengan menggunakan Biomass Expansion Factor memiliki nilai biomassa yang lebih besar dari biomassa alometrik Hendri. Tetapi pada kelas umur I (KU I) dan II (KU II) biomassa hasil perhitungan BEF lebih kecil dari biomassa alometrik Hendri. Nilai biomassa minimum sebesar 4.95 ton/ha pada KU I dan nilai biomassa tertinggi sebesar 461.71 ton/ha pada KU XI. Gambar 11 pada daerah dengan KU yang tinggi menunjukkan perhitungan biomassa dengan menggunakan BEF menghasilkan nilai biomassa yang lebih tinggi jika dibandingkan dengan nilai biomassa alometrik Hendri.
37
Gambar 11 Perbandingan nilai biomassa alometrik Hendri dan biomassa BEF. Besarnya nilai pendugaan dengan menggunakan BEF dipengaruhi beberapa faktor, diantaranya adalah BEF yang digunakan untuk menghitung biomassa tidak dihasilkan dari data pada daerah penelitian. Nilai BEF pada tegakan jati ini dikembangkan oleh Kraenzel et al. (2003) berdasarkan data perhitungan biomassa tegakan jati secara destruktif di daerah Panama. Selain itu penghitungan BEF ini dikonsentrasikan pada pohon jati berusia 0 - 20 tahun (KU I - KU II), sedangkan pada daerah penelitian pohon jati yang diambil sebagai sampel pengukuran memiliki umur yang bervariasi. Penggunaan koefisien BEF dalam pembuatan model dimaksudkan untuk menghasilkan model pendugaan yang bersifat general sehingga model pendugaan diharapkan dapat digunakan untuk menduga biomassa jati di daerah manapun. Tetapi karena koefisien BEF yang digunakan untuk menghitung biomassa dihasilkan dari data yang tidak mewakili kondisi tegakan jati di KPH Kebonharjo sehingga pada penelitian ini penggunaan koefisien BEF tidak disarankan. Hasil pembuatan model pendugaan dengan koefisien BEF tetap dilampirkan. Berdasarkan analisis terhadap nilai pendugaan biomassa dengan menggunakan BEF maka pembuatan model pendugaan biomassa dikonsentrasikan pada hubungan biomassa alometrik Hendri dengan backscatter citra ALOS PALSAR.
38
5.3.2 Pola Hubungan Backscatter Citra ALOS PALSAR dengan Biomassa Penyusunan model pendugaan biomassa di atas permukaan tanah pada tegakan jati dilakukan dengan menganalisis sebaran data backscatter yang diekstraksi dari citra satelit ALOS PALSAR dengan nilai biomassa lapangan yang dihitung berdasarkan model alometrik Hendri.
a
b
c
d
Gambar 12 Scatterplot hubungan antara biomassa Hendri dengan backscatter ALOS PALSAR : a) backscatter HV resolusi 50 m, b) backscatter HH resolusi 50 m, c) backscatter HV resolusi 12,5 m, d) backscatter HH resolusi 12,5 m. Gambar 12 menunjukan pola sebaran pada hubungan backscatter citra ALOS PALSAR resolusi 50 m dan 12,5 m dengan biomassa Hendri. Berdasarkan sebaran data pada Gambar 12 dapat dilihat bahwa data memiliki pola yang tidak linier sehingga pembuatan model dilakukan dengan memperhatikan kaidah regresi nonlinier yaitu membuat jumlah kuadrat galat sekecil mungkin. Lucas et al. (2006) menemukan pola hubungan backscatter pada citra radar dengan nilai biomassa di atas permukaan membentuk pola yang nonlinier. Sehingga model terbaik ditentukan oleh nilai RMSE terkecil, Fhitung terbesar, koefisien determinasi disesuaikan (R2adj) terbesar, dan nilai signifikansi > 0.05.
39
Dengan memperhatikan bentuk kurva pada sebaran data (Baker 2006) dimana hubungan variabel respon dan prediktornya mengikuti bentuk fungsional tertentu (Smyth 2002). Pada persamaan yang dibangun dengan kaidah nonlinier, nilai koefisien determinasi disesuaikan dan RMSE menjelaskan keeratan hubungan nonlinier antara biomassa dan nilai tengah backscatter. Pada plot-plot dengan kelas umur muda, kondisi vegetasi di lapangan yang rapat, pendugaan biomassa dengan menggunakan backscatter cenderung memiliki nilai yang lebih besar, atau sebaliknya pada plot-plot dengan kelas umur tua dan memiliki kondisi vegetasi di lapangan yang kurang rapat, pendugaan biomassa dengan menggunakan backscatter akan lebih kecil. Hal ini dikarenakan hubungan antara biomassa dan backscatter dipengaruhi oleh banyak faktor, diantaranya adalah heterogenitas atau homogenitas hutan, topografi, tutupan tajuk, dan salah satunya adalah kerapatan tegakan. 5.3.3 Penyusunan Model Hubungan Backscatter dengan Biomassa Alometrik Hendri Hasil analisis hubungan antara nilai biomassa Alometrik Hendri dan backscatter citra ALOS PALSAR resolusi 50 m ditunjukkan pada Tabel 7 dan Tabel 8. Baik backscatter HH maupun HV menunjukkan terdapat hubungan dengan biomassa walaupun hubungan yang dihasilkan dari regresi antara biomassa dengan backscatter HH maupun HV tidak terlalu besar. Pada Tabel 8 dapat dilihat bahwa hubungan backscatter HV dengan biomassa menunjukkan nilai statistik lebih baik dibandingkan dengan hubungan backscatter HH dengan biomassa yang disajikan pada Tabel 7. Pada studi-studi yang telah dilakukan sebelumnya didapatkan bahwa polarisasi HV lebih sensitif dalam menduga nilai biomassa di atas permukaan pada kondisi permukaan yang datar dibandingkan di tempat yang bergelombang (Wijaya 2009). Dimana pada lokasi penelitian, daerah-daerah pengambilan sampel relatif memiliki permukaan yang datar. Nilai backscatter HV akan menghasilkan nilai pendugaan biomassa yang lebih baik karena memiliki mekanisma hamburan lebih dominan terhadap hamburan balik tajuk (Franson & Israelsson 1999 dalam Wijaya 2009).
40
Berdasarkan Tabel 7 model untuk menduga biomassa terbaik dengan menggunakan backscatter HH adalah model 4 dengan nilai RMSE pada model = 38,25; Fhitung = 173,12 dan R2adj = 60,7% pada bentuk persamaan Y = 318,289+(1/0,03)BS_HH. Pada model pendugaan 1 sampai dengan model 4, nilai koefisien determinasi disesuaikan (R2adj ) adalah 30 - 60%. Hal ini berarti variasi dari variabel tak bebas (biomassa) dapat dijelaskan oleh variabel bebas (backscatter HH) sebesar 30 – 60% menggunakan biomassa alometrik Hendri. Tabel 7 Model pendugaan biomassa berdasarkan hubungan biomassa Hendri dengan backscatter HH citra ALOS PALSAR resolusi 50 m Model 1 2 3 4
Bentuk persamaan Y=Exp(6,45+0,287BS_HH) Y=BS_HH/(0,066+0,021BS_HH) Y=23,214*(Exp(-9,145/BS_HH)) Y=318,289+(1/0,03)BS_HH
R2adj 53.40% 36.50% 43.70% 60.70%
RMSE 41.62 48.44 45.77 38.25
Fhit 143.59 101.51 115.76 173.12
p-value 0 0 0 0
Ftabel 3.276 3.276 3.276 3.276
Y = Biomassa (ton/ha); BS_HH = nilai backscatter HH pada Citra ALOS PALSAR 50 m
Berdasarkan Tabel 8 model pendugaan biomassa terbaik dengan menggunakan backscatter HV adalah model 1 dimana hasil regresi biomassa Hendri dan backscatter HV pada model 1 menunjukan nilai RMSE terkecil sebesar 27,70; Fhitung = 345,49 dan R2adj = 79,4% pada bentuk persamaan Y = Exp(9,291+0,38BS_HV). Pada model 1 sampai model 4 (Tabel 8), nilai R2adj sebesar 60 - 80%, hal ini menunjukkan bahwa variasi dari variabel tak bebas (biomassa) dapat dijelaskan oleh variabel bebas (backscatter HV) sebesar 60 – 80%. Tabel 8 Model pendugaan biomassa berdasarkan hubungan biomassa Hendri dengan backscatter HV citra ALOS PALSAR resolusi 50 m Model 1 2 3 4
Bentuk Persamaan Y=Exp(9,291+0,38BS_HV) Y=BS_HV/(0,431+0,045BS_HV) Y=1,089*(Exp(-55,31/BS_HV)) Y=479,331+(1/0,034)BS_HV
R2adj 79.40% 69.90% 77.00% 75.30%
RMSE 27.7 33.67 29.27 30.28
Fhit 345.49 228.37 307.71 286.23
p-value 0 0 0 0
Ftabel 3.276 3.276 3.276 3.276
Y = Biomassa (ton/ha); BS_HV = nilai backscatter HV pada Citra ALOS PALSAR 50 m
Model pendugaan biomassa yang diformulasikan berdasarkan hubungan biomassa Hendri dengan backscatter citra ALOS PALSAR resolusi 12,5 m dapat dilihat pada Tabel 9 dan Tabel 10. Berdasarkan Tabel 9 maka model pendugaan biomassa terbaik dengan menggunakan backscatter HH citra ALOS PALSAR resolusi 12,5 m terdapat pada model 1, dengan nilai RMSE = 37,29 Fhitung =183,03
41
dan R2adj = 62,20% pada bentuk persamaan Y=Exp(6,676+0,274BS_HH125). Tabel 9 memperlihatkan nilai koefisien determinasi disesuaikan (R2adj) pada setiap model (model 1 - 4) berkisar antara 40 sampai 60%, hal ini berarti bahwa variasi dari biomassa dapat dijelaskan 40 sampai 60 % oleh backscatter HH pada citra ALOS PALSAR resolusi 12,5 m. Tabel 9 Model pendugaan biomassa berdasarkan hubungan biomassa Hendri dengan backscatter HH citra ALOS PALSAR resolusi 12,5 m. Model 1 2 3 4
Bentuk Persamaan Y=Exp(6,676+0,274BS_HH125) Y=BS_HH125/(0,081+0,021BS_HH125) Y=19,443*(Exp(-12,123/BS_HH125)) Y=272,917+(1/0,047)BS_HH125
R2adj 62.20% 43.60% 53.20% 59.60%
RMSE 37.29 45.78 41.74 38.74
Fhit 183.03 115.67 142.64 168.30
p-value 0 0 0 0
Ftabel 3.276 3.276 3.276 3.276
Y = Biomassa (ton/ha); BS_HH125 = nilai backscatter HH pada Citra ALOS PALSAR 12,5 m
Berdasarkan Tabel 10 maka model pendugaan biomassa terbaik dengan menggunakan backscatter HV citra alos palsar resolusi 12,5 m terdapat pada model 1, dengan nilai RMSE = 32,43 Fhitung = 247,43 dan R2adj = 71,70% pada bentuk persamaan Y=Exp(8,811+0,302BS_HV125). Tabel 10 memperlihatkan nilai koefisien determinasi disesuaikan (R2adj) pada setiap model (model 1 - 4) berkisar antara 50 sampai 70%, hal ini berarti bahwa variasi dari biomassa dapat dijelaskan 50 sampai 70 % oleh backscatter HV pada citra ALOS PALSAR resolusi 12,5 m. Tabel 10 Model pendugaan biomassa berdasarkan hubungan biomassa Hendri dengan backscatter HV citra ALOS PALSAR resolusi 12,5 m. Model 1 2 3 4
Bentuk Persamaan Y=Exp(8,811+0,302BS_HV125) Y=BS_HV125/(0,387+0,038BS_HV125) Y=1,933*(Exp(-54,587/BS_HV125)) Y=374,285+(1/0,054)BS_HV125
R2adj 71.70% 58.60% 68.50% 63.30%
RMSE 32.43 39.21 34.23 36.96
Fhit 247.43 163.84 220.42 186.54
p-value 0 0 0 0
Ftabel 3.276 3.276 3.276 3.276
Y = Biomassa (ton/ha); BS_HV125 = nilai backscatter HV pada Citra ALOS PALSAR 12,5 m
Pada model pendugaan biomassa terbaik dengan backscatter HH dan backscatter HV baik dengan menggunakan citra ALOS PALSAR resolusi 50 m, maupun citra ALOS PALSAR 12,5 m didapatkan nilai F hitung > F tabel, hal ini menunjukkan bahwa variabel backscatter HH maupun backscatter HV mempengaruhi nilai biomassa, didukung dengan nilai koefisien determinasi disesuaikan yang cukup besar maka kedua variabel bebas baik backscatter HH maupun backscatter HV berpengaruh nyata terhadap nilai pendugaan biomassa.
42
Berdasarkan nilai RMSE dan koefisien determinasi disesuaikan pada modelmodel yang telah dibuat dapat dikatakan bahwa pendugaan biomassa dengan mengguanakan backscatter HV akan menghasilkan nilai pendugaan yang lebih baik. Hal tersebut dapat dilihat dari nilai statistik pada model dugaan biomassa di atas permukaan yang telah dibuat, dimana backscatter HV menunjukan hubungan yang lebih baik dibandingkan hubungan biomassa dengan backscatter HH. Seluruh model pendugaan dengan menggunakan backscatter HV sebagai variabel bebas menunjukkan pola hubungan yang lebih baik jika dibandingkan dengan hubungan biomassa dengan backscatter HH, hal ini dikarenakan polarisasi HV memiliki sensitifitas lebih baik terhadap komponen penyusun biomassa yaitu batang dan tutupan tajuk. 5.3.4 Verifikasi Model Verifikasi model dilakukan untuk mengetahui apakah nilai dugaan biomassa yang dihasilkan oleh model terpilih tidak berbeda dengan nilai biomassa di lapangan. Hasil verifikasi model dengan menggunakan analisis uji T berpasangan disajikan pada Tabel 10 dan Tabel 11. Apabila Thit > Ttabel maka tolak Ho atau nilai signifikasnsi < 0,05. Apabila Thit < Ttabel maka terima Ho atau nilai signifikansi > 0,05. Dimana Hipotesis uji yang diberlakukan adalah sebagai berikut : H0 : µ1 - µ2 = 0 (Biomassa aktual = biomassa model) H1 : µ1 - µ2 ≠ 0 (Biomassa aktual ≠ biomassa model) Berdasarkan hasil uji T berpasangan yang dilakukan pada model pendugaan dengan variabel bebas backscatter ALOS PALSAR resolusi 50 m, terdapat 4 (empat) model terbaik yang dapat digunakan dalam menduga biomassa, dua model pertama merupakan model pendugaan biomassa berdasarkan backscatter ALOS PALSAR 50 m yang dihasilkan dari regresi biomassa Hendri dengan backscatter HH maupun HV. Sedangkan dua model yang lain merupakan model pendugaan biomassa berdasarkan backscatter ALOS PALSAR 12,5 m yang dihasilkan dari regresi biomassa alometrik Hendri dengan backscatter HH maupun HV.
43
Pada Tabel 11 dapat dilihat bahwa model 1 sampai dengan model 4 memiliki nilai T hitung berturut-turut sebesar 0.182, dan 0.444, nilai ini lebih kecil dari T tabel (2.056) dan nilai signifikansi lebih besar dari 0,05 sehingga dapat diartikan bahwa nilai pendugaan biomassa tidak berbeda nyata dengan nilai biomassa aktual di lapangan. Lampiran 5 menyajikan hasil perhitungan dengan menggunakan model hubungan terpilih pada citra ALOS PALSAR resolusi 50 m. Tabel 11 Hasil verifikasi model pendugaan biomassa pada citra ALOS PALSAR resolusi 50 m No.
Model Terpilih
R2adj
RMSE
T hit
T tabel
sig.
1 2
Alometrik Hendri – Backscatter Y = 318,289 + 1/0.030×BS_HH Y = Exp (9,291 + 0.38×BS_HV)
60,7% 79,4%
38,25 27.70
0,182 0,444
2,056 2,056
0,857 0,661
Y = Biomassa, BS = Backscatter citra ALOS PALSAR
Model pendugaan biomassa di atas permukaan dengan menggunakan polarisasi HH (Lampiran 5) pada citra ALOS PALSAR 50 m menunjukkan hasil yang tidak berbeda jauh dengan hasil pendugaan biomassa dengan polarisasi HV. Sehingga baik pada model pendugaan biomassa yang dihasilkan dari hubungan biomassa alometrik Hendri dengan backscatter HH citra ALOS PALSAR 50 m maupun model pendugaan biomassa yang dihasilkan dari hubungan biomassa alometrik Hendri dengan backscatter HV citra ALOS PALSAR 50 m dapat digunakan dalam pendugaan biomassa di atas permukaan. Hasil verifikasi model hubungan biomassa alometrik Hendri dan backscatter citra ALOS PALSAR 12,5 m disajikan pada Tabel 12. Nilai T hitung pada model 3 dan model 4 sebesar 0.647 dan 0.745 menunjukkan bahwa T hitung setiap model bernilai lebih kecil dari T tabel (2.056) dan nilai signifikansi lebih besar dari 0,05 sehingga dapat diartikan bahwa nilai pendugaan biomassa tidak berbeda nyata dengan nilai biomassa aktual di lapangan baik pada model hubungan biomassa model alometrik Hendri dengan backscatter HH, maupun model hubungan biomassa alometrik Hendri dengan backscatter HV. Tabel 12 Hasil verifikasi model pendugaan biomassa pada citra ALOS PALSAR resolusi 12,5 m No.
Model Terpilih
R2adj
RMSE
T hit
T tabel
sig.
1 2
Alometrik Hendri – Backscatter Y = Exp(6.676 + 0.274×BS_HH125) Y = Exp(8.811 + 0.302×BS_HV125)
62.20% 71.70%
37.28 32.43
0.647 0.745
2,056 2,056
0.523 0.463
Y = Biomassa, BS = Backscatter citra ALOS PALSAR
44
Model pendugaan biomassa di atas permukaan dengan menggunakan polarisasi HH menunjukkan hasil yang tidak berbeda jauh dengan hasil pendugaan biomassa dengan polarisasi HV (Lampiran 7). Sehingga baik pada model pendugaan biomassa yang dihasilkan dari hubungan biomassa alometrik hendri dengan backscatter HH citra ALOS PALSAR 12,5 m maupun model pendugaan biomassa yang dihasilkan dari hubungan biomassa alometrik Hendri dengan backscatter HV citra ALOS PALSAR 12,5 m dapat digunakan dalam pendugaan biomassa di atas permukaan. Pendugaan biomassa dengan menggunakan polarisasi HV lebih baik dari pendugaan biomassa dengan menggunakan polarisasi HH. Hal ini dibuktikan dengan nilai RMSE, R2adj dan nilai signifikansi yang lebih baik dari pendugaan biomassa dengan menggunakan polarisasi HH (Tabel 11 dan Tabel 12). Bukti ini didukung oleh hasil penelitian sebelumnya dengan menggunakan Citra ALOS PALSAR, dimana hubungan antara polarisasi HV dengan biomassa lebih baik dibandingkan polarisasi HH dengan biomassa (Sarker & Nichol 2010). Selanjutnya Sarker dan Nichol (2010) menjelaskan bahwa polarisasi dua arah (HV atau VH) lebih sensitif terhadap volume maka pendugaan biomassa dengan menggunakan polarisasi dua arah menjadi lebih baik dibandingkan pendugaan biomassa dengan menggunakan polarisasi satu arah (HH atau VV). Pada pembuatan model pendugaan biomassa dengan menggunakan polarisasi HV selalu memberikan nilai statistik dan nilai verifikasi yang lebih baik bila dibandingkan dengan polarisasi HH. Hamburan balik (backscatter) pada radar L-band menembus sebagian daun dan ranting dimana hamburan terbesar dihasilkan oleh batang dan cabang, pada kelas umur tua yang memiliki jarak antar pohon semakin jauh diakibatkan penjarangan, nilai biomassa yang dapat di duga oleh backscatter radar cenderung kecil. Ditemukan pada beberapa plot yang diukur di lapangan, pada daerah yang memiliki jarak tanam yang teratur sehingga memungkinkan mekanisme backscatter dikembalikan dengan langsung tanpa ada gap yang terlalu besar antar pohon, membuat nilai pendugaan biomassa cukup besar. Polarisasi searah menuju sensor dihasilkan oleh pantulan energi gelombang mikro dari kanopi, batang, atau cabang. Backscatter yang dihasilkan cenderung
45
sangat kuat dan direkam sebagai objek yang memiliki rona cerah pada citra radar dengan polarisasi searah (HH atau VV). Kecerahan rona ini menandakan nilai dijital yang tinggi sehingga backscatter yang dihasilkan pun tinggi. Sebaliknya pada backscatter dengan polarisasi dua arah (HV atau VH) dihasilkan dari pancaran energi gelombang mikro yang menembus permukaan kanopi suatu pohon, pada daun, cabang dan batangnya. Hamburan balik (backscatter) dua arah ini seringkali disebut sebagai hamburan balik Volume. Nilai dijital yang relatif rendah ditandakan dengan rona yang tidak terlalu cerah, sehingga nilai backscatter pada citra dengan polarisasi ini cenderung rendah. Variasi yang terjadi pada nilai pendugaan biomassa dipengaruhi oleh jarak antar individu pohon (Sarker & Nichol 2010). Reflektansi dari biomassa (baik biomassa total maupun komponen biomassa) pada SAR backscatter akan lebih baik pada hutan yang homogen, baik yang memiliki kesamaan bentuk, tinggi, maupun diameter, dimana pada hutan yang homogen alokasi biomassa pada setiap komponen tegakan yang berbeda akan relatif lebih konsisten (Lucas et al. 2004). 5.4 Pembuatan Peta Sebaran Biomassa dan Analisis Akurasi Peta sebaran biomassa dibuat berdasarkan empat model terbaik yang disajikan pada Tabel 13 Pembuatan peta sebaran biomassa dilakukan dengan dua cara, yaitu dengan menggunakan image (citra) hasil pengolahan dimana informasi nilai dijital pada image ini telah dirubah menjadi nilai biomassa dan klasifikasi pada image dilakukan dengan menganalisis rata-rata setiap nilai 5×5 piksel untuk citra 50 m dan 20×20 piksel untuk citra 12,5 m. Cara yag kedua adalah klasifikasi dilakukan dengan menganalisis nilai rata-rata biomassa pada setiap piksel dalam petak areal kerja KPH Kebonharjo. Peta sebaran biomassa dibuat ke dalam tiga kelas biomassa. Kelas biomassa dibuat berdasarkan perubahan bentuk kurva biomassa terendah hingga tertinggi. Selang kelas biomassa yang dihasilkan disajikan pada Tabel 14.
46
Tabel 13 Model pendugaan biomassa terbaik ALOS PALSAR resolusi 50 m dan ALOS PALSAR resolusi 12,5 m Model Bentuk Model Keterangan 1
Y = 318,289 + 1/0.030×BS_HH
2
Y = Exp (9,291 + 0.38×BS_HV)
3
Y = Exp(6.676 + 0.274×BS_HH125)
4
Y = Exp(8.811 + 0.302×BS_HV125)
BS_HH = Backscatter HH Citra ALOS PALSAR resolusi 50 m, BS_HV = Backscatter HV Citra ALOS PALSAR resolusi 50 m, BS_HH125 = Backscatter HH Citra ALOS PALSAR resolusi 12,5 m, BS_HV125 = Backscatter HV Citra ALOS PALSAR resolusi 12,5 m.
Selang kelas biomassa pada Tabel 14 digunakan sebagai acuan klasifikasi dalam pembuatan peta sebaran biomassa berdasarkan model hubungan Biomassa alometrik Hendri dengan nilai backscatter pada Citra ALOS PALSAR baik resolusi 50 m maupun resolusi 12,5 m. Pada klasifikasi biomassa yang telah dilakukan kemudian dihitung akurasi pengkelasan yang telah dibuat tersebut. Tabel 14 Selang kelas biomassa yang digunakan untuk membuat peta sebaran biomassa Kelas Biomassa Kisaran Nilai Biomassa Hendri (ton/ha) 1 2 3
0 ~ 80.00 80.00 ~ 160.00 160.00 ~ up
5.4.1 Peta Sebaran Biomassa Berdasarkan model pendugaan biomassa pada Tabel 13 dibuat peta sebaran biomassa pada citra ALOS PALSAR resolusi 50 m dan citra ALOS PALSAR resolusi 12,5 m. Model 1 dan model 2 merupakan model hubungan antara backscatter citra ALOS PALSAR resolusi 50 m dengan biomassa alometrik Hendri yang digunakan untuk membuat peta sebaran biomassa pada citra ALOS PALSAR resolusi 50 m. Model 3 dan model 4 merupakan model hubungan antara backscatter citra ALOS PALSAR resolusi 12,5 m dengan biomassa alometrik Hendri yang digunakan untuk membuat peta sebaran biomassa pada citra ALOS PALSAR resolusi 12,5 m. Dengan menggunakan model-model terpilih dihasilkan 8 (delapan) peta sebaran biomassa pada daerah penelitian (Gambar 13 sampai Gambar 16). Rekapitulasi hasil perhitungan akurasi pada setiap peta sebaran biomassa yang telah dibuat disajikan pada Tabel 15 dan Tabel 16.
47
5.4.2
Uji Akurasi Peta Sebaran Biomassa
Tabel 15 Rekapitulasi hasil perhitungan overall accuracy dan kappa accuracy pada peta sebaran dengan citra ALOS PALSAR resolusi 50 m Akurasi (%) Model Bentuk Persamaan Piksel Petak OA KA OA KA 1 2
Y = 318,289 + 1/0.030×BS_HH Y = Exp (9,291 + 0.38×BS_HV)
50.79 53.97
20.93 27.3
55.56 55.56
28.64 30.06
Berdasarkan hasil rekapitulasi perhitungan akurasi pada Tabel 15, model 1 menghasilkan overall accuracy sebesar 50,79% dan kappa accuracy sebesar 20,93% pada analisis akurasi klasifikasi dengan menggunakan rata-rata nilai biomassa dalam 5×5 piksel (analisis akurasi peta per piksel) pada Gambar 13a. Peta sebaran biomassa dibuat berdasarkan nilai setiap piksel pada citra ALOS PALSAR resolusi 50 m polarisasi HH ditunjukkan dengan Gambar 13a, sebaran kelas biomassa dominan berada pada kelas biomassa 1 dan kelas biomassa 2. Titik pengamatan menyebar pada kelas biomassa 2, terdapat kesalahan interpretasi pada pengkelasan dengan model 1 sebanyak 29 titik. Kesalahan interpretasi terbanyak ditemukan pada interpretasi kelas biomassa 3, dimana kelas biomassa 3 pada titik diklasifikasikan sebagai kelas biomassa 2 sebanyak 17 titik dan kelas biomassa 1 diklasifikasikan sebagai kelas biomassa 2 sebanyak 7 titik. Banyakya kesalahan interpretasi menyebabkan kecilnya nilai kappa accuracy. Nilai overall accuracy 50,79% menjelaskan kesalahan interpretasi 29 titik dari 63 titik sebaran biomassa. Pada Gambar 13b ditunjukkan sebaran biomassa pada areal kerja KPH Kebonharjo berada pada kelas biomassa 1 dan kelas biomassa 2. Dengan menggunakan peta areal kerja KPH Kebonharjo, penghitungan akurasi peta sebaran menghasilkan overall accuracy sebesar 55,56% dan kappa accuracy sebesar 28,64% (Tabel 15). Titik pengamatan menyebar pada kelas biomassa 2, dengan kesalahan interpretasi sebanyak 29 titik dari 63 titik sebaran biomassa. Hasil interpretasi pada peta sebaran biomassa areal kerja KPH Kebonharjo (Gambar 13b) menunjukkan kesalahan interpretasi kelas biomassa 3 menjadi kelas biomassa 2 sebanyak 16 titik dan kelas biomassa 1 menjadi kelas biomassa 2 sebanyak 7 titik. Berdasarkan Gambar 13b terlihat bahwa kelas biomassa 3 tidak terlalu banyak menyebar pada peta sebaran menggunakan peta areal kerja KPH Kebonharjo.
48
a)
b) Gambar 13 Peta sebaran biomassa modelb 1, a) peta sebaran biomassa per piksel, ) b) peta sebaran biomassa berdasarkan peta areal kerja KPH Kebonharjo. Gambar 14 menunjukan peta sebaran yang dibuat berdasarkan model 2, yaitu hubungan backscatter HV dengan biomassa alometrik Hendri. Pada peta sebaran biomassa Gambar 14a, kelas biomassa menyebar pada kelas biomassa 1 serta kelas biomassa 2, dimana terdapat sekumpulan kecil kelas biomassa 3 yang jika dibandingkan dengan peta sebaran model 1 (Gambar 13) sebaran kelas biomassa 3 terlihat lebih banyak. Titik pengamatan menyebar pada kelas biomassa 2, terdapat kesalahan interpretasi pada pengkelasan dengan model 2 sebanyak 29 titik.
49
Kesalahan interpretasi terbanyak ditemukan pada interpretasi kelas biomassa 3, dimana kelas biomassa 3 diklasifikasikan sebagai kelas biomassa 2 sebanyak 14 titik dan kelas biomassa diklasifikasikan sebagai kelas biomassa 2 sebanyak 6 titik. Kappa accuracy yang dihasilkan sebesar 27,30% dengan nilai overall accuracy 53,97%. Terjadi kenaikan kappa accuracy pada klasifikasi yang dilakukan dengan model 2 jika dibandingkan dengan model 1, hal ini karena jumlah titik kesalahan interpretasi pada pengklasifikasian kelas biomassa 3 pada model 2 lebih sedikit jika dibandingkan dengan model 1.
a)
b) Gambar 14 Peta sebaran biomassa model 2, a) peta sebaran biomassa per piksel, b) peta sebaran biomassa berdasarkan peta areal kerja KPH Kebonharjo.
50
Dengan menggunakan model 2, peta sebaran biomassa pada areal kerja KPH Kebonharjo pada Gambar 14b menunjukkan bahwa kelas biomassa 1 dan kelas biomassa 2 mendominasi sebaran biomassa pada daerah penelitian. Titik pengamatan menyebar pada kelas biomassa 2, dimana terdapat kesalahan interpretasi sebanyak 28 titik. Kesalahan interpretasi dalam pengklasifikasian kelas biomassa 3 sebagai kelas biomassa 2 sebanyak 14 titik dan kesalahan interpretasi dalam pengklasifikasian kelas biomassa 1 sebagai kelas biomassa 2 sebanyak 5 titik, kappa accuracy yang dihasilkan sebesar 30,06% dan overall accuracy yang dihasilkan sebesar 55,56%. Pada pembuatan peta sebaran dengan menggunakan model 2, pembuatan peta sebaran dengan peta areal kerja KPH Kebonharjo menunjukan akurasi yang lebih besar dari peta sebaran biomassa per piksel. Hasil perhitungan akurasi pada peta sebaran biomassa dengan menggukana citra ALOS PALSAR 12,5 m ditunjukkan pada Tabel 16. Pada Tabel 16 dapat dilihat bahwa akurasi terbesar pada pembuatan peta sebaran biomassa dengan menggunakan citra ALOS PALSAR 12,5 m, yaitu pada pembuatan peta sebaran biomassa dengan menggunakan model 8. Peta sebaran biomassa yang dibuat berdasarkan Model 3 disajikan pada Gambar 15. Berdasarkan Gambar 15a hasil perhitungan akurasi pada peta sebaran menghasilkan overall accuracy 53,98% dan kappa accuracy 26,77%. Kesalahan klasifikasi pada kelas biomassa 3 menjadi kelas biomassa 2 sebanyak 13 titik, kelas biomassa 1 menjadi kelas biomassa 2 sebanyak 7 titik. Kesalahan interpretasi dari klasifikasi yang dilakukan pada peta sebaran biomassa model 5 sebanyak 29 titik dari 63 titik. Pada peta sebaran ini biomassa menyebar di kelas biomassa 1 dan kelas biomassa 2.
Tabel 16 Rekapitulasi hasil perhitungan overall accuracy dan kappa accuracy pada citra ALOS PALSAR resolusi 12,5 m Akurasi (%) Model 3 4
Bentuk Persamaan Y = Exp(6.676 + 0.274*BS_HH125) Y = Exp(8.811 + 0.302*BS_HV125)
Piksel OA 53.97 55.56
KA 26.77 29.75
Petak OA KA 57.14 31.49 60.32 37.15
51
Pada perhitungan akurasi untuk Gambar 15b menghasilkan overall accuracy sebesar 57,14% dan kappa accuracy sebesar 31,49%. Kesalahan klasifikasi pada kelas biomassa 3 menjadi kelas biomassa 2 sebanyak 14 titik. Sedangkan kesalahan interpretasi klasifikasi kelas biomassa 1 menjadi kelas biomassa 2 sebanyak 7 titik. Kesalahan interpretasi total pada peta sebaran biomassa model 3 dengan menggunakan peta areal kerja KPH Kebonharjo sebanyak 27 titik dari 63 titik. Pada peta sebaran Gambar 15 terlihat bahwa kelas biomassa 3 terkonsentrasi pada Bagian Hutan Tuder.
a)
b) Gambar 15 Peta sebaran biomassa model 3, a) peta sebaran biomassa per piksel, b) peta sebaran biomassa berdasarkan peta areal kerja KPH Kebonharjo.
52
Peta sebaran dengan menggunakan Model 4 disajikan pada Gambar 16, Perhitungan akurasi pada peta sebaran biomassa Gambar 16a menghasilkan overall accuracy 55,56% dan kappa accuracy 29,75%. Kesalahan klasifikasi kelas biomassa 3 menjadi kelas biomassa 2 sebanyak 13 titik, kesalahan klasifikasi kelas biomassa 1 menjadi kelas biomassa 2 sebanyak 6 titik. Kesalahan interpretasi pada klasifikasi dalam peta sebaran biomassa model 4 Gambar 16a sebanyak 28 titik dari 63 titik. Pada peta sebaran ini biomassa menyebar pada kelas biomassa 1 dan kelas biomassa 2. Hasil perhitungan akurasi pada peta sebaran biomassa Gambar 16b, yaitu overall accuracy 60,32% dan kappa accuracy 37,15%. Sebaran biomassa berada pada kelas biomassa 1 dan kelas biomassa 2, kelas biomassa 3 ditemukan lebih banyak terkonsentrasi di Bagian Hutan Tuder. Kesalahan klasifikasi Kelas biomassa 3 yang diinterpretasikan menjadi kelas biomassa 2 sebanyak 14 titik. Kesalahan klasifikasi kelas biomassa 1 yang diinterpretasikan menjadi kelas biomassa 2 sebanyak 5 titik. Dari 4 model yang digunakan dalam pembuatan peta sebaran biomassa, model 4 memiliki nilai akurasi paling besar dengan bentuk persamaan Y = Exp(8.811 + 0.302×BS_HV125). Gambar 16b memiliki nilai akurasi paling besar. Pada peta sebaran biomassa Gambar 16b dengan menggunakan peta areal kerja KPH Kebonharjo dapat dilihat bahwa kelas biomassa 1, kelas biomassa 2, dan kelas biomassa 3 memiliki sebaran yang merata, sehingga kesalahan interpretasi pada kelas biomassa 3 dan kelas biomassa 1 menjadi lebih sedikit. Beberapa daerah yang memiliki nilai biomassa yang tinggi (kelas biomassa 3) terklasifikasi menjadi daerah dengan biomassa yang rendah (kelas biomassa 1 atau kelas biomassa 2). Hal ini dikarenakan penggunaan peta petak sebagai area of interest menghasilkan informasi nilai tengah yang tidak hanya mewakili nilai piksel tertentu. Pada pendugaan biomassa per piksel pada setiap titik pengamatan tidak selalu mewakili keadaan dilapangan. Karena nilai backscatter pada setiap piksel dipengaruhi oleh
banyak faktor, antara lain topografi, arah sensor,
konstanta dielektrik, dan kerapatan vegetasi. nilai kappa akurasi yang lebih besar ini terjadi karena pada pengkelasan dengan model 4 selain jumlah kelas yang benar diinterpretasikan di temukan lebih banyak, penyebaran kelas biomassa
53
cenderung lebih menyebar jika dibandingkan peta sebaran biomassa dengan model lain.
a
b Gambar 16 Peta sebaran biomassa model 4, a) peta sebaran biomassa per piksel, b) peta sebaran biomassa berdasarkan peta areal kerja KPH Kebonharjo.
BAB VI KESIMPULAN DAN SARAN 6.1 Kesimpulan 1. Didapatkan 4 model terbaik dari 16 model yang telah dibuat, dua model pendugaan dengan variabel bebas backscatter citra ALOS PALSAR resolusi 50 m dan dua model pendugaan biomassa dengan manggunakan variabel backscatter citra ALOS PALSAR 12,5 m. 2. Model pendugaan biomassa terbaik untuk hubungan biomassa Hendri dengan backscatter citra ALOS PALSAR 50 m untuk polarisasi HH adalah Y = 318,289 + 1/0.030×BS_HH dan untuk polarisasi HV adalah Y = Exp (9,291 + 0.38×BS_HV). Sedangkan model pendugaan terbaik untuk hubungan biomassa Hendri dengan backscatter citra ALOS PALSAR 12,5 m untuk polarisasi HH adalah Y = Exp(6.676 + 0.274×BS_HH125) dan untuk polarisasi HV adalah Y = Exp(8.811 + 0.302×BS_HV125). 3. Hubungan antara backscatter HH dengan biomassa tidak terlalu erat, hal ini dapat dilihat dari nilai koefisien determinasi disesuaikan yang rendah dan error yang lebih tinggi jika dibandingkan dengan hubungan antara biomassa dengan backscatter HV. 4. Uji t-berpasangan menunjukkan bahwa empat model terpilih memiliki nilai pendugaan biomassa yang tidak berbeda nyata dengan biomassa aktual di lapangan sehingga keempat model tersebut dapat digunkan dalam pendugaan biomassa di daerah penelitian untuk jenis tegakan yang sama dan kondisi lapang yang sama. 5. Pembuatan peta sebaran menghasilkan 8 peta dengan akurasi yang beragam, overall accuracy pada setiap peta berkisar pada kurang lebih 50% sedangkan nilai kappa accuracy pada setiap peta berkisar pada 20% - 40%. 6. Peta sebaran biomassa dengan nilai akurasi tertinggi adalah peta sebaran dengan model Y = Exp(8.811 + 0.302×BS_HV125) (model 4).
56
6. 2 Saran 1. Peta sebaran dengan menggunakan peta kerja KPH akan jauh lebih baik untuk menduga biomassa karena mempertimbangkan karakteristik tegakan pada setiap anak petak, tetapi perlu dilakukan analisis akurasi dengan pembuatan kelas biomassa pada selang yang beragam. 2. Perlu dilakukan uji statistik pada 4 model pendugaan pada daerah lain dengan kondisi lapangan dan kondisi tegakan yang sama, atau pada daerah lain dengan keadaan yang berbeda tetapi memiliki jenis tegakan yang sama. 3. Untuk pendugaan biomassa dengan menggunakan koefisien BEF sebaiknya dilakukan dengan menghitung koefisien BEF di daerah penelitian terlebih dahulu, sehingga hasil pendugaan biomassa yang di harapkan dapat tercapai. 4. Digunakan variabel lain selain backscatter untuk menduga biomassa misalnya dengan dilakukan skoring pada komponen yang mempengaruhi nilai backscatter (contoh : topografi, kelembaban, bonita).
DAFTAR PUSTAKA Austin JM, Mackey BG, Van Niel KP. 2003. Estimating forest biomass using satellite radar : an exploratory study in a temperate Australian Eucalyptus forest. [Jurnal] Forest Ecology and Manajement Vol. 176 Hal. 575 – 583. Baker S L. 2006. Non-Linier Regression. http://hspm.sph.sc.edu/courses/J716 /pdf/716-5%20Non-linear%20regression.pdf. [13 Januari 2011]. Bergen MK, Dobson MC. 1999. Integration of remotely sensed Radar imagery in modelling and mapping of Forest Biomassa and Net Primary Production. [Jurnal] Ecological Modelling ELSEVIER. No. 122 Hal. 257-274. Brown S. 1997. Estimating biomass and biomass change of tropical forests. FAO Forest Resources Assessment Publication No.134 Hal 55. Roma. Brown S, Gillespie AJR, Lugo AE.1989. Biomass Estimation Methods for Tropical Forest with Applications to Forest Inventory Data. [Jurnal] Forest Science, Vol 3 No.4 Hal. 881-902. Chave J, Riera B, Dubois MA. Estimation of biomass in neotropical forest of Frence Guiana : spatial and temporal variability. [Jurnal] Journal of Tropical Ecology 2001 No. 17 Hal 79 - 96. Dahlan. 2005. Pendugaan Kandungan Karbon Tegakan Acacia Mangium Willd Menggunakan Citra Landsat ETM+ ean SPOT-5: Studi Kasus di BKPH Parung Panjang KPH Bogor [tesis]. Bogor: Istitut Pertanian Bogor. Dobson MC, Ulaby F, Le Toan T, Beaudoin A, Kasischkle ES, Christensen N. 1992. Dependence of radar backscatter on coniferous forest biomass. [Jurnal] IEEE Transaction on Geoscience and Remote Sensing Vol. 10 No. 2 Hal. 412-415. Dong J, Kaufmann RK, Myneni RB, Tucker CJ, Kauppi P, Liski J. 2003. Remote sensing estimates of boreal and temperate forest woody biomass: Carbon pools, sources, and sinks. [Jurnal] Remote Sensing of Environment Vol. 83 Hal 393-410. [ESA] European Space Agency. 2010. RADAR (Radio Detection and Raging). http://earth.esa.int/applications/data_util/SARDOCS/spaceborne/Radar_ Courses/Radar_Course_III/parameters_affecting.htm . [31 Januari 2011] Fransson JES, Israelsson H. 1999. Estimation of stem volume in boreal forests using ERS-1 C- and JERS-1 L-band SAR data. [Jurnal] International Journal of Remote Sensing Vol. 20 No. 1 Hal 123-137. Fussel J, Rundquist D, Harrington JA Jr. 1986. On defining remote sensing, photogrametry. [Jurnal] Eng & Remote Sensing, Vol 52 Hal. 1507-1511.
58
Hendri. 2001. Analysis of emission and absorption of green hous gasses and evaluation of carbon dioxide mitigation technologies in Perum Perhutani [tesis]. Bogor: Institut Pertanian Bogor. Husch B, Beers TW, Kershaw JA. 2003. Forest measurement : 4th ed. New Jersey, United States : John Wiley 7 Sons. [IPCC] Intergovernmental Panel on Climate Change. 1995. Greenhouse Gas inventory reference manual IPCC WGI technical support unit. Hardley Center, Meteorology Office. London Road, Braknell, RG 122 NY, United Kingdom. [JAXA] Japan Aerospace Exploration Agency. 2010. PALSAR (Phased Array type L-band Synthetic Aperture Radar). http://www.eorc.jaxa.jp/ALOS/en/abo ut/palsar.htm [31 Maret 2010]. Jaya INS. 2002. Aplikasi SIG untuk Kehutanan. Bogor: Fakultas Kehutanan Institut Pertanian Bogor. Johnson GO. 1992. GIS Applications in Emergency Management. [Jurnal] URISA Journal. Vol. 4 No. 1 Hal 66-72 Khaerudin. 1994. Pembibitan Tanaman HTI. Jakarta : Penebar Suadaya. KPH Kebonharjo. 2000. Tabel Volume Lokal Tebangan A.2 Jati. Jawa Tengah : Biro Perencanaan Perum Perhutani Unit I. KPH Kebonharjo. 2007. RPKH Kebonharjo 2007-2016. Jawa Tengah : Biro Perencanaan Perum Perhutani Unit I. Kraenzel M, Castillo A, Moore T, Potvin C. 2003. Carbon Storage of Harvest-Age Teak (Tectona grandis) Plantation, Panama. [Jurnal] Forest Ecology and Management Vol. 173 Hal 213 – 225. Kumar S. 2009. Retrieval of Forest Parameters from Envisat ASAR Data for Biomass Inventory in Dudhwa National Park, UP, India. [disertasi] Indian Institute of Remote Sensing. Lillesand TM, Kiefer RW. 1979. Image Interpretation and Remote Sensing. New York : John Wiley&Son Inc. Lu D. 2005. The potential and challenge of remote sensing-based biomass estimation. [Jurnal] International Journal of Remote Sensing, Vol. 27 No.7 Hal. 1297-1328. Lucas RM, Cronin N, Lee A, Moghaddam M, Witte C, Tickle P. 2006. Empirical relationship between AIRSAR backscatter and LiDAR-derived forest biomass, Queensland, Australia. [Jurnal]. Remote sensing Environtment Vol. 100 Hal. 407 – 425. Lucas RM, Moghaddam M, Cronin N. 2004. Microwave scattering from mixed species woodlands, central Queensland, Australia. [Jurnal]. IEEE
59
Transactions on Geoscience and Remote Sensing Vol. 42 Hal : 2141– 2159. Lugo AE, Snedaker SC. 1974. The Ecology Of Mangroves. [Jurnal] Ann. Rev. Ecol. Syst 5 : 39-65. Mattjik A, Sumertajaya M. 2000. Perancangan Percobaan: Aplikasi SAS dan Minitab. Jilid I. Bogor: IPB Press. Mitchard ETA, Saatchi SS, Woodhouse IH, Nangendo G, Ribeiro NS, Williams M, Ryan CM, Lewis SL, Feldpausch TR, Meir P. 2009. Using Sattelite Radar Backscatter to Predict Above-ground Woody Biomass : A Consistent Relationship Across Four Different African Landscapes. [Jurnal] Geophysical Research Letters. Vol. 36 [NASA] National Aeronautics and Space Administration. Radar Data Center. Freeman Tony. 1996. What is Radar image. Jet Propulsion Laboratory. http://southport.jpl.nasa.gov/. [1 Januari 2011]. Onrizal. 2004. Model penduga biomassa dan karbon tegakan hutan kerangan di Taman Nasional Danau Sentarum, Kalimantan Barat. [Tesis]. Bogor: Institut Pertanian Bogor. Paryono P. 1994. Sistem Informasi Geografis. Yogyakarta : Andi Offset. Prahasta E. 2009. Sistem Informasi Geografis: Konsep-Konsep Dasar (Perspektif Geodesi dan Geomatika ). Bandung : Informatika Bandung. Purwadhi FSH. 2001. Interpretasi Citra Digital. Jakarta : Grasindo. Rosenqvist A, Masanobu S, Manabu W. 2004. ALOS PALSAR: Technical outline and mission concepts Samsuri. 2004. Aplikasi Penginderaan Jauh dalam Pengelolaan Sumber Daya Hutan. [tesis] Jurusan Kehutanan, Fakultas Pertanian, Universitas Sumatera Utara. Sarker MLR, Nichol JE. 2010. Biomass Estimation using Texture Measurements of Dual Polarization SAR Data. Proceeding of the remote sensing and photogrammetry society conference. London. Shimada M, Isoguchi O, Tadono T, Isono K. 2009. PALSAR Calibration Factor Updated. https://auig.eoc.jaxa.jp/auigs/en/doc/an/20090109en_3.html [5 November 2010] Smyth GK. 2002. Nonlinear regression. Encyclopedia of Environmetrics. John Wiley&Sons, Ltd, Chichester. Telenius B, Verwijst T. 1994. The influence of allometric variation, vertical biomass distributions and sampling procedure on biomass estimates in
60
commercial short-rotation forests. [Jurnal] Bioresource Technology No. 51 (1995) Hal. 247-253. Tiryana T, Tatsuhara S, Shiraishi N. 2011. Empirical model for estimating the stand biomass of Teak Plantation in Java, Indonesia. [Jurnal] Journal of Forest Planning No. 16 Hal. 177 – 188. Wijaya A. 2009. Evalusai data mosaik ALOS PALSAR untuk estimasi volume dan biomassa batang (studi kasus dari hutan hujan tropis Indonesia Bagian Tengah). [Jurnal] Jurnal Geografi No. 2 [1 Januari 2009].
LAMPIRAN
63
Lampiran 1 Hasil Perhitungan Biomassa pada Plot Pengamatan dengan menggunakan Alometrik Tegakan Jati yang digunakan sebagai model maupun validasi. No. Plot_ID Kelas Umur Biomassa (Ton/ha) 1 50032A VIII 220.3224 2 50039A II 135.2635 3 50011 VI 164.9799 4 50012A IV 174.8776 5 50014 VIII 177.6879 6 50018A III 125.0429 7 50019A III 120.1237 8 50020 II 78.45688 9 50021 V 149.5169 10 50024A I 39.82471 11 50028A II 94.92379 12 50033 VIII 175.2259 13 50034 III 100.8341 14 50038 II 116.0863 15 50040 VI 169.5666 16 50041 VII 156.0531 17 50042 VI 86.7726 18 50043 II 111.1655 19 50048 VII 184.1373 20 50055 IX 199.0695 21 50069 III 117.2698 22 50071A IV 106.8484 23 50073 I 12.28208 24 50075 V 178.9403 25 50077 I 40.26815 26 50078 I 31.20841 27 50079 I 82.58772 28 50080 I 56.11469 29 50082 I 32.86854 30 50084 I 7.51783 31 50089 II 159.9619 32 50092 VI 127.5499 33 50094 I 69.04149 34 50096 I 13.48198 35 50052A I 7.799567 36 50001A XI 328.6952 37 50005A* I 31.9759 38 50010* V 90.88745 39 50013A* IV 151.28 40 50015* II 116.1468
64 Lampiran 1 (Lanjutan) No. Plot_ID Kelas Umur 41 50016A* VII 42 50025* III 43 50026A* V 44 50027A* II 45 50029A* III 46 50035* III 47 50051A* II 48 50053* VI 49 50070* III 50 50076* I 51 50081* I 52 50083* II 53 50085* I 54 50086* IV 55 50087* III 56 50088* II 57 50091* II 58 50093* IV 59 50095* I 60 50097* II 61 50056A* VIII Keterangan : *) data yang digunakan untuk validasi
Biomassa (Ton/ha) 235.0398 129.85 99.6854 112.9509 197.2791 62.41236 173.1695 153.8217 104.3491 67.94988 89.35894 96.75522 40.76888 182.3471 186.559 88.69974 85.15532 80.49109 62.57825 84.21514 271.6375
65
Lampiran 2 Hasil Perhitungan Biomassa pada Plot Pengamatan dengan menggunakan Biomass Expansion Factor (BEF) yang digunakan sebagai model maupun validasi. No. Plot_ID KU BEF Biomassa BEF (Ton/ha) 1
50039A
II
1.53
133.15
2
50011
VI
1.53
198.48
3
50014
VIII
1.53
222.61
4
50018A
III
1.53
117.4
5
50019A
III
1.53
114.21
6
50020
II
1.53
68.66
7
50021
V
1.53
168.63
8
50024A
I
1.53
31.84
9
50028A
II
1.53
90.6
10
50033
VIII
1.53
230.83
11
50034
III
1.53
96.64
12
50038
II
1.53
114.38
13
50040
VI
1.53
223.85
14
50041
VII
1.53
205.05
15
50042
VI
1.53
105.37
16
50043
II
1.53
106.35
17
50048
VII
1.53
232.99
18
50069
III
1.53
116.56
19
50071A
IV
1.53
108.83
20
50073
I
1.53
8.45
21
50077
I
1.53
30.15
22
50078
I
1.53
24.22
23
50079
I
1.53
76.85
24
50080
I
1.53
45.51
25
50082
I
1.53
24.76
26
50084
I
1.53
4.95
27
50089
II
1.53
156.78
28
50092
VI
1.53
151.38
29
50094
I
1.53
63.28
30
50096
I
1.53
9.96
31
50052A
I
1.53
5.18
32
50074
V
1.53
124.84
33
50070*
III
1.53
97.92
34
50053*
VI
1.53
203.1675
35
50056A*
VIII
1.53
403.6516
36
50051A*
II
1.53
397.6085
37
50035*
III
1.53
58.95235
38
50076*
I
1.53
64.23962
66 Lampiran 2 (Lanjutan) No. Plot_ID KU
BEF
39
50029A*
III
1.53
40
50027A*
II
1.53
116.259
41
50010*
V
1.53
99.62628
42
50013A*
IV
1.53
169.4382
43
50025*
III
1.53
129.1012
44
50026A*
V
1.53
109.4966
45
50001A*
XI
1.53
461.1759
46
50081*
I
1.53
77.47637
47
50083*
II
1.53
88.05862
48
50015*
II
1.53
102.9446
49
50016A*
VII
1.53
299.4584
50
50005A*
I
1.53
24.24776
51
50085*
I
1.53
32.02931
52
50086*
IV
1.53
215.6581
53
50087*
III
1.53
194.3402
54
50088*
II
1.53
85.73939
55
50091*
II
1.53
76.07846
56
50093*
IV
1.53
80.54343
57
50095*
I
1.53
51.97796
58
50097*
II
1.53
79.40788
Keterangan : *) data yang digunakan untuk validasi
Biomassa BEF (Ton/ha) 210.8615
67
Lampiran 3 Hasil Perhitungan Nilai Backscatter pada Plot Pengamatan dengan menggunakan citra ALOS PALSAR resolusi 50 m. Digital Number Backscatter (dB) No. Plot_ID HH HV HH HV 1
50032A
8443
4201.306
-4.47006
-10.5323
2
50039A
8596.472
3915.056
-4.31359
-11.1452
3
50001A
7295.556
3594.917
-5.73883
-11.8862
4
50005A
5381.056
2384.833
-8.38265
-15.4508
5
50010
7614
3573.361
-5.36774
-11.9385
6
50011
7542.694
4008.833
-5.44947
-10.9396
7
50012A
7127.833
3489.306
-5.94085
-12.1452
8
50013A
7754.262
3847.548
-5.20919
-11.2963
9
50014
7137.333
3921.389
-5.92928
-11.1312
10
50015
5649.861
3307.889
-7.95924
-12.609
11
50016A
7094.556
3576.611
-5.9815
-11.9306
12
50018A
6547.917
3742.639
-6.67794
-11.5364
13
50019A
6252.583
3567.833
-7.07881
-11.9519
14
50020
5755.778
3114.056
-7.79792
-13.1335
15
50021
8164.861
3946.75
-4.76102
-11.0752
16
50024A
5996.25
2929.528
-7.44241
-13.664
17
50025
8446.944
3873.222
-4.46601
-11.2386
18
50026A
7896.111
3827.278
-5.05173
-11.3422
19
50027A
6344.139
3080.056
-6.95255
-13.2288
20
50028A
6015.778
3261.583
-7.41416
-12.7314
21
50029A
5209.472
2721.694
-8.66413
-14.3032
22
50033
8084.111
4147.25
-4.84735
-10.6448
23
50034
7963.75
3826.25
-4.97765
-11.3445
24
50035
8201.028
3997.222
-4.72263
-10.9648
25
50038
8350.528
3944.667
-4.56572
-11.0798
26
50040
7961.444
3860.25
-4.98016
-11.2677
27
50041
7647
3798.861
-5.33018
-11.4069
28
50042
7065.472
3423.361
-6.01718
-12.3109
29
50043
7598.639
3659.889
-5.38528
-11.7306
30
50048
7593.972
4134.556
-5.39062
-10.6714
31
50051A
7155.306
3781.194
-5.90744
-11.4474
32
50053
7367.889
3698.083
-5.65314
-11.6405
33
50055
6889.528
3626.667
-6.23621
-11.8098
34
50069
8244.944
3873.306
-4.67625
-11.2384
35
50070
7952.611
3952.5
-4.98981
-11.0626
36
50071A
6784.361
3490.917
-6.36982
-12.1412
37
50073
4642.528
1952.528
-9.66491
-17.1881
38
50075
6750.056
3645.139
-6.41385
-11.7657
39
50076
6238.917
3603
-7.09782
-11.8667
40
50077
6439.25
2910.167
-6.82329
-13.7216
68 Lampiran 3 (Lanjutan) No.
Plot_ID
Digital Number HH HV
Backscatter (dB) HH HV
41
50078
5901.444
2502.806
-7.58083
-15.0315
42
50079
6790.361
3380.528
-6.36214
-12.4203
43
50080
6505.583
2987.222
-6.73428
-13.4946
44
50081
6424
2942
-6.84389
-13.6271
45
50082
5357.476
2297.881
-8.4208
-15.7734
46
50083
7642.389
3537.167
-5.33542
-12.0269
47
50084
5213.139
2461.167
-8.65801
-15.1772
48
50085
6653.083
3226.833
-6.53954
-12.8245
49
50086
8074.861
3542.611
-4.8573
-12.0135
50
50087
7505
3507.917
-5.49299
-12.099
51
50088
8189.048
4091.833
-4.73533
-10.7616
52
50089
7461.417
3714
-5.54357
-11.6032
53
50090
7343.5
3434.917
-5.68194
-12.2817
54
50091
7528.639
3695.028
-5.46567
-11.6476
55
50092
7943.306
3744.333
-4.99997
-11.5325
56
50093
7613.222
3936.806
-5.36863
-11.0971
57
50094
6637
3365.694
-6.56056
-12.4585
58
50095
7126.944
3397.472
-5.94193
-12.3769
59
50096
5295.278
2475.333
-8.52223
-15.1273
60
50097
7179.444
3677.722
-5.87818
-11.6884
61
50052A
4780.639
1998.556
-9.41028
-16.9857
62
50056A
6889.472
3712.972
-6.23628
-11.6056
63
50074
6717.806
3269.111
-6.45545
-12.7114
69
Lampiran 4 Hasil Perhitungan Nilai Backscatter pada Plot Pengamatan dengan menggunakan citra ALOS PALSAR resolusi 12,5 m. No. Plot_ID Digital Number Backscatter (dB) HH HV HH HV 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52
50032A 50039A 50001A 50005A 50010 50011 50012A 50013A 50014 50015 50016A 50018A 50019A 50020 50021 50024A 50025 50026A 50027A 50028A 50029A 50033 50034 50035 50038 50040 50041 50042 50043 50048 50051A 50053 50055 50069 50070 50071A 50073 50075 50076 50077 50078 50079 50080 50081 50082 50083 50084 50085 50086 50087 50088 50089
7051.55 6924.35 6473.70 5617.61 7275.00 7742.28 7046.49 5950.35 7749.85 6460.44 6817.68 6494.71 6013.26 5879.42 7854.01 6539.49 7831.90 7419.04 4448.62 5611.15 5076.99 7121.85 6409.07 6422.57 6743.55 6223.65 6047.45 5801.71 5902.68 6314.63 5763.74 5927.24 6803.75 6476.54 6529.20 5945.11 2927.96 6481.68 6128.93 5197.21 3697.67 5546.17 6654.61 5672.28 4471.53 5935.38 4773.58 6230.93 6477.68 6658.54 6014.94 6579.00
3246.57 3165.07 3056.72 1947.88 3194.17 3421.20 3240.28 2971.19 3718.83 3510.86 3041.53 3539.13 3417.86 3290.91 3416.46 2697.07 3406.58 3387.65 2364.79 3138.97 2738.25 3411.02 3163.32 2999.70 3015.69 3120.03 2988.54 2608.95 2867.93 3231.04 3209.98 2814.07 3334.16 3107.36 3008.99 2973.42 1174.26 3197.94 3085.76 2294.43 1559.40 2663.66 2610.97 2499.34 2086.21 2762.33 1982.90 3057.13 2936.97 3158.36 2850.43 3095.23
-6.03431 -6.19242 -6.77695 -8.00896 -5.76334 -5.22262 -6.04054 -7.50915 -5.21414 -6.79476 -6.32727 -6.74880 -7.41780 -7.61331 -5.09817 -6.68912 -5.12265 -5.59305 -10.03549 -8.01897 -8.88788 -5.94814 -6.86409 -6.84582 -6.42223 -7.11910 -7.36856 -7.72888 -7.57901 -6.99304 -7.78591 -7.54295 -6.34503 -6.77313 -6.70280 -7.51680 -13.66871 -6.76625 -7.25231 -8.68460 -11.64144 -8.12013 -6.53754 -7.92484 -9.99087 -7.53103 -9.42311 -7.10894 -6.77161 -6.53242 -7.41537 -6.63680
-12.77149 -12.99234 -13.29489 -17.20875 -12.91284 -12.31643 -12.78834 -13.54138 -11.59187 -12.09173 -13.33817 -12.02208 -12.32492 -12.65368 -12.32848 -14.38216 -12.35364 -12.40203 -15.52415 -13.06427 -14.25053 -12.34230 -12.99714 -13.45845 -13.41226 -13.11681 -13.49082 -14.67070 -13.84862 -12.81316 -12.86996 -14.01329 -12.54027 -13.15217 -13.43157 -13.53488 -21.60475 -12.90259 -13.21277 -15.78649 -19.14087 -14.49044 -14.66395 -15.04350 -16.61284 -14.17448 -17.05397 -13.29373 -13.64201 -13.01078 -13.90180 -13.18616
70 Lampiran 4 (Lanjutan) No. Plot_ID 53 54 55 56 57 58 59 60 61 62 63
50090 50091 50092 50093 50094 50095 50096 50097 50052A 50056A 50074
Digital Number HH HV 5107.60 6233.28 6838.34 6334.77 5865.24 5407.60 3025.36 6301.30 3139.05 6266.76 6251.43
2345.64 3012.85 3025.49 3064.41 2751.42 2512.88 1242.38 2983.44 1349.52 3407.47 3033.03
Backscatter (dB) HH HV -8.83566 -7.10566 -6.30098 -6.96538 -7.63428 -8.33990 -13.38446 -7.01140 -13.06404 -7.05914 -7.08041
-15.59477 -13.42044 -13.38408 -13.27307 -14.20886 -14.99655 -21.11492 -13.50565 -20.39642 -12.35135 -13.36245
71
Lampiran 5 Hasil pendugaan biomassa dengan menggunakan model hubungan yang disusun berdasarkan alometrik Hendri dan backscatter citra ALOS PALSAR 50 m. Backsatter (dB) Biomassa (Ton/ha) Plot_ID HH HV Aktual Model HH Model HV 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26
50001A 50005A 50010 50013A 50015 50016A 50025 50026A 50027A 50029A 50035 50051A 50053 50070 50076 50081 50083 50085 50086 50087 50088 50091 50093 50095 50097 50056A
-5.7388 -8.3827 -5.3677 -5.2092 -7.9592 -5.9815 -4.4660 -5.0517 -6.9525 -8.6641 -4.7226 -5.9074 -5.6531 -4.9898 -7.0978 -6.8439 -5.3354 -6.5395 -4.8573 -5.4930 -4.7353 -5.4657 -5.3686 -5.9419 -5.8782 -6.2363
-11.8862 -15.4508 -11.9385 -11.2963 -12.6090 -11.9306 -11.2386 -11.3422 -13.2288 -14.3032 -10.9648 -11.4474 -11.6405 -11.0626 -11.8667 -13.6271 -12.0269 -12.8245 -12.0135 -12.0990 -10.7616 -11.6476 -11.0971 -12.3769 -11.6884 -11.6056
328.70 31.98 90.89 151.28 116.15 235.04 129.85 99.69 112.95 197.28 62.41 173.17 153.82 104.35 67.95 89.36 96.76 40.77 182.35 186.56 88.70 85.16 80.49 62.58 84.22 271.64
126.99 38.87 139.36 144.65 52.98 118.91 169.42 149.90 86.54 29.48 160.87 121.37 129.85 151.96 81.70 90.16 140.44 100.30 156.38 135.19 160.44 136.10 139.33 120.22 122.35 110.41
118.42 30.56 116.09 148.18 89.98 116.44 151.46 145.61 71.10 47.27 168.07 139.91 130.01 161.94 119.30 61.11 112.26 82.91 112.83 109.22 181.56 129.66 159.83 98.28 127.66 131.75
72
Lampiran 6 Hasil pendugaan biomassa dengan menggunakan model hubungan yang disusun berdasarkan biomassa BEF dan backscatter citra ALOS PALSAR 50 m. No.
Plot_ID
Backsatter (dB) HH HV
Biomassa (Ton/ha) Aktual
Model HH
Model HV
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26
50001A 50005A 50010 50013A 50015 50016A 50025 50026A 50027A 50029A 50035 50051A 50053 50070 50076 50081 50083 50085 50086 50087 50088 50091 50093 50095 50097 50056A
-5.7388 -8.3827 -5.3677 -5.2092 -7.9592 -5.9815 -4.4660 -5.0517 -6.9525 -8.6641 -4.7226 -5.9074 -5.6531 -4.9898 -7.0978 -6.8439 -5.3354 -6.5395 -4.8573 -5.4930 -4.7353 -5.4657 -5.3686 -5.9419 -5.8782 -6.2363
461.18 24.25 99.63 169.44 102.94 299.46 129.10 109.50 116.26 210.86 58.95 397.61 203.17 97.92 64.24 77.48 88.06 32.03 215.66 194.34 85.74 76.08 80.54 51.98 79.41 403.65
129.35 27.67 143.62 149.72 43.95 120.02 178.31 155.78 82.67 16.84 168.44 122.87 132.65 158.16 77.08 86.85 144.87 98.56 163.26 138.81 167.95 139.86 143.59 121.54 123.99 110.22
116.21 18.08 113.08 158.11 79.69 113.55 162.96 154.37 57.66 32.91 187.98 146.12 132.12 178.63 117.40 46.83 107.98 71.21 108.74 103.99 209.02 131.62 175.44 89.95 128.85 134.54
-11.8862 -15.4508 -11.9385 -11.2963 -12.6090 -11.9306 -11.2386 -11.3422 -13.2288 -14.3032 -10.9648 -11.4474 -11.6405 -11.0626 -11.8667 -13.6271 -12.0269 -12.8245 -12.0135 -12.0990 -10.7616 -11.6476 -11.0971 -12.3769 -11.6884 -11.6056
73
Lampiran 7 Hasil pendugaan biomassa dengan menggunakan model hubungan yang disusun berdasarkan Alometrik Hendri dan backscatter citra ALOS PALSAR 12,5 m. Backsatter (dB) Biomassa (Ton/ha) No. Plot_ID HH HV Aktual Model HH Model HV 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26
50001A 50005A 50010 50013A 50015 50016A 50025 50026A 50027A 50029A 50035 50051A 50053 50070 50076 50081 50083 50085 50086 50087 50088 50091 50093 50095 50097 50056A
-7.23175 -8.43285 -6.06309 -7.77862 -7.12956 -6.70649 -5.41553 -5.85515 -10.3725 -9.22457 -7.17405 -8.11739 -7.84003 -6.97592 -7.71618 -8.34602 -7.84856 -7.48381 -7.1108 -6.87783 -7.71302 -7.39711 -7.2736 -8.64526 -7.31446 -7.34754
-13.6251 -17.6619 -13.1384 -13.7962 -12.3777 -13.7196 -12.5498 -12.6032 -15.7688 -14.5111 -13.6954 -13.1118 -14.2917 -13.6338 -13.582 -15.3627 -14.4591 -13.6474 -13.8922 -13.2071 -14.1216 -13.6356 -13.4865 -15.2238 -13.7024 -12.5852
328.7 31.98 90.89 151.28 116.15 235.04 129.85 99.69 112.95 197.28 62.41 173.17 153.82 104.35 67.95 89.36 96.76 40.77 182.35 186.56 88.7 85.16 80.49 62.58 84.22 271.64
109.34 78.68 150.61 94.13 112.45 126.27 179.85 159.44 46.24 63.34 111.09 85.78 92.56 117.28 95.75 80.58 92.34 102.05 113.03 120.48 95.84 104.50 108.10 74.23 106.89 105.93
109.53 32.37 126.88 104.02 159.64 106.45 151.56 149.13 57.33 83.82 107.23 127.90 89.56 109.24 110.97 64.81 85.14 108.80 101.04 124.27 94.28 109.19 114.21 67.59 107.01 149.95
74
Lampiran 8 Hasil pendugaan biomassa dengan menggunakan model hubungan yang disusun berdasarkan biomassa BEF dan backscatter citra ALOS PALSAR 12,5 m. Backsatter (dB) Biomassa (Ton/ha) No. Plot_ID HH HV Aktual Model HH Model HV 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26
50001A 50005A 50010 50013A 50015 50016A 50025 50026A 50027A 50029A 50035 50051A 50053 50070 50076 50081 50083 50085 50086 50087 50088 50091 50093 50095 50097 50056A
-7.23175 -8.43285 -6.06309 -7.77862 -7.12956 -6.70649 -5.41553 -5.85515 -10.3725 -9.22457 -7.17405 -8.11739 -7.84003 -6.97592 -7.71618 -8.34602 -7.84856 -7.48381 -7.1108 -6.87783 -7.71302 -7.39711 -7.2736 -8.64526 -7.31446 -7.34754
-13.6251 -17.6619 -13.1384 -13.7962 -12.3777 -13.7196 -12.5498 -12.6032 -15.7688 -14.5111 -13.6954 -13.1118 -14.2917 -13.6338 -13.582 -15.3627 -14.4591 -13.6474 -13.8922 -13.2071 -14.1216 -13.6356 -13.4865 -15.2238 -13.7024 -12.5852
461.18 24.25 99.63 169.44 102.94 299.46 129.1 109.5 116.26 210.86 58.95 397.61 203.17 97.92 64.24 77.48 88.06 32.03 215.66 194.34 85.74 76.08 80.54 51.98 79.41 403.65
109.66 75.66 157.36 92.61 113.18 128.98 192.21 167.80 41.55 59.24 111.63 83.41 90.87 118.68 94.42 77.72 90.63 101.44 113.84 122.33 94.51 104.20 108.25 70.85 106.89 105.81
109.56 27.77 129.28 103.37 167.44 106.09 157.92 155.07 52.86 81.06 106.97 130.45 87.34 109.23 111.18 60.68 82.51 108.73 100.05 126.29 92.54 109.17 114.85 63.62 106.72 156.03
75
Lampiran 9 Bentuk
Kelas 1 Kelas 1 10 Model Kelas 2 7 Kelas 3 1 Total 18 PA % 55.56 Overall Accuracy Kappa Accuracy
Bentuk
Kelas 1 Kelas 1 11 Model Kelas 2 7 Kelas 3 0 Total 18 PA % 61.11 Overall Accuracy (%) Kappa Accuracy (%) Bentuk Kelas 1 Model Kelas 1 11 Kelas 2 6 Kelas 3 1 Total 18 PA % 61.11 Overall Accuracy (%) Kappa Accuracy (%)
Bentuk
Kelas 1 Kelas 1 12 Model Kelas 2 5 Kelas 3 1 Total 18 PA % 66.67 Overall Accuracy (%) Kappa Accuracy (%)
Titik Kelas 2 2 21 3 26 80.77
Titik Kelas 2 3 22 1 26 84.62
Titik Kelas 2 2 19 5 26 73.08
UA % 76.92 46.67 20.00
Total Kelas 3 1 13 17 45 1 5 19 63 5.26 50.79 20.93
Total
Kelas 3 1 16 2 19 10.53 55.56 28.64
15 45 3 63
Total Kelas 3 1 14 4 19 21.05
UA % 73.33 48.89 66.67
UA % 78.57 48.72 40.00
14 39 10 63 53.97 27.30
Titik Kelas 2 3 19 4 26 73.08
Kelas 3 1 14 4 19 21.05 55.56 30.06
Total 16 38 9 63
UA % 75.00 50.00 44.44