$type=ticker$count=10$cols=4$cate=0

Gap fill citra menggunakan Ermapper

Gap Fill dilakukan pada citra satelit Landsat 7 dikarenakan sejak 2003 terjadi kerusakan pada Scan Line Corrector (SLC), satelit L7 berad...


Gap Fill dilakukan pada citra satelit Landsat 7 dikarenakan sejak 2003 terjadi kerusakan pada Scan Line Corrector (SLC), satelit L7 berada pada mode SLC OFF. Meskipun pada saat sekarang sudah ada penggantinya yaitu Landsat 8, untuk kepentingan analisa temporal pada tahun-tahun sebelum ada Landsat 8 tetap saja diperlukan citra Landsat 7. Beberapa software sudah menyediakan tool untuk Gap Fill citra Landsat 7. Tutorial ini hanya salah satu cara saja dengan menggunakan Ermapper. Pun demikian, dengan Ermapper pun banyak cara yang bisa dilakukan untuk melakukan Gap Fill.

Tutorial ini menggunakan dua citra dengan ID LE71220642009202SGS00 dan LE71220642009218SGS00 liputan tahun 2009. Kedua citra tersebut sudah diinspeksi terlebih dahulu dan telah yakin bahwa gap pada citra yang satu dapat diisi dengan data pada citra yang lainnya. Pada banyak kasus mungkin kita perlu 3 atau lebih liputan. Tutorial ini memanfaatkan citra GAP_MASK yang disediakan pada saat kita mendownload data citra Landsat 7. Gap Fill hanya dilakukan pada band 5, 4, dan 3 saja.

1. Buat Dataset

Untuk menjalankan perhitungan dalam algorithma, kita harus membuat dataset terlebih dahulu, yaitu menggabungkan semua data yang ingin dianalisa dalam satu dataset. Buat ermapper, stack layer seperti pada gambar.

Gambar 1.


Nama B31 artinya band 3 pada liputan 1, B31m artinya band 3 liputan 1 gap_mask. Jika ada yg perlu ditanyakan bagaimana membuat dataset seperti pada gambar, silakan komentar.

2. Simpan Algorithm

Algorithma disimpan ke dalam satu file dataset (ERS)

Gambar 2.

3. Buka data hasil export pada Langkah 2 di atas.


Perhatikan nilai pixel citra seperti pada gambar berikut.

Gambar 3.

Perhatikan gambar di atas. Saat di-Klik pada strip, tampil nilai2 pixel pada setiap layer (12 layer). Jumlah layer adalah jumlah liputan (2) dikali dengan jumlah band (3 bands yaitu B3, B4, B5) dan jumlah data/mask (2), sehingga total ada 2 x 3 x 2 = 12 layer.

Kita dapat memanfaatkan fungsi IF THEN ELSE dengan data pada layer dataset. IF THEN ELSE dapat dikembangkan dari logika dengan contoh sebagai berikut
> Untuk layer RED, jika masking band 5 pada liputan pertama adalah SATU, maka diambilah data band 5 pada liputan pertama, jika tidak SATU, maka diambil band 5 pada liputan kedua.

4. Membuat komposit B543

Pada jendela algorithma, Klik pada EMC2 untuk setiap layer. Isikan formula di bawah ini untuk setiap layer
RED dengan formula IF INPUT1 = 1 THEN INPUT2 ELSE INPUT3
GREEN dengan formula IF INPUT1 = 1 THEN INPUT2 ELSE INPUT3
BLUE dengan formula IF INPUT1 = 1 THEN INPUT2 ELSE INPUT3

Gambar 4.



Pada setiap layer, tentukan INPUT1, INPUT2 dan INPUT3 seperti pada gambar di bawah ini

Gambar 5. Formula untuk RED (B5)



Gambar 6. Formula untuk GREEN (B4)


Gambar 7. Formula untuk BLUE (B3)]

 5. Hasil

Simpan algorithma ke file baru. Dapat dipilih OUTPUT TYPE RGB

Gambar 8.

6. Gap fill untuk 3 liputan citra

Tahapan 1 - 5 ditujukan untuk Gap fill citra dengan 2 liputan saja, Jika berdasarkan inspeksi pendahuluan ditentukan perlu 3 buah liputan, maka lakukan langkah 1 - 5 dengan tambahan citra liputan baru.

Khusus untuk Langkah 4, algorithma akan menjadi seperti gambar di bawah ini.

Gambar 9. Contoh algorithma menggunakan 3 buah liputan untuk Band 5

Catatan:1. Gap fill untuk analisa spektral seharusnya dilakukan pada citra yg sudah dikoreksi. Namun untuk analisis non-spektral, misalnya analisis visual, Gap fill dapat dilakukan terhadap citra belum terkoreksi.
2. Tutorial ini memanfatkan citra GAP_MASK yang disediakaan. Pilihan algorithma lain tentu ada misalnya menggunakan definisi nilai GAP itu sendiri. Namun hal tersebut tidak dibahas pada tutorial ini.

/Have fun.

COMMENTS

BLOGGER: 4
  1. mau nanya mas, dalam pembuatan data set untuk komposit 543 seperti diatas, darimana mendapatkan band 32,42 dan 52?
    makasih mas

    BalasHapus
  2. Band 52, 42 dan 32 adalah merujuk kepada band 543 citra pengisi, sedangkan 51, 41 dan 31 adalah merujuk kepada band 543 citra master (yang diisi).
    Menentukan mana citra master (1) dan mana citra pengisi (1) ditentukan secara subyektif.

    BalasHapus
  3. cara menggabung dataset seperti di langkah 1 itu bgaimana mas? terimma kasih sebelumnya

    BalasHapus
  4. mohon maaf, ada modul untuk pengolahan garis pantai menggunakan er mapper atau tidak ya. kalau ada mohon berbagi mas.

    BalasHapus

Nama

Data,56,Geospasial,147,GIS,98,Pustaka,20,RS,96,Software,31,Tutorial,69,
ltr
item
RSGIS INDONESIA: Gap fill citra menggunakan Ermapper
Gap fill citra menggunakan Ermapper
https://2.bp.blogspot.com/-Xxeg-l9Za78/XJG0dqYD4jI/AAAAAAABwWM/qTyk9VttZfwtbBFQkdPt4N-bOL543Q00wCEwYBhgL/s640/gap05.png
https://2.bp.blogspot.com/-Xxeg-l9Za78/XJG0dqYD4jI/AAAAAAABwWM/qTyk9VttZfwtbBFQkdPt4N-bOL543Q00wCEwYBhgL/s72-c/gap05.png
RSGIS INDONESIA
https://www.rsgis.info/2015/03/gap-fill-citra-menggunakan-ermapper_6.html
https://www.rsgis.info/
https://www.rsgis.info/
https://www.rsgis.info/2015/03/gap-fill-citra-menggunakan-ermapper_6.html
true
3008467515520484700
UTF-8
Loaded All Posts Not found any posts VIEW ALL Readmore Reply Cancel reply Delete By Home PAGES POSTS View All RECOMMENDED FOR YOU LABEL ARCHIVE SEARCH ALL POSTS Not found any post match with your request Back Home Sunday Monday Tuesday Wednesday Thursday Friday Saturday Sun Mon Tue Wed Thu Fri Sat January February March April May June July August September October November December Jan Feb Mar Apr May Jun Jul Aug Sep Oct Nov Dec just now 1 minute ago $$1$$ minutes ago 1 hour ago $$1$$ hours ago Yesterday $$1$$ days ago $$1$$ weeks ago more than 5 weeks ago Followers Follow THIS PREMIUM CONTENT IS LOCKED STEP 1: Share. STEP 2: Click the link you shared to unlock Copy All Code Select All Code All codes were copied to your clipboard Can not copy the codes / texts, please press [CTRL]+[C] (or CMD+C with Mac) to copy