How to destrip (gap filling) Landsat 7 SLC-off

Cara mengisi informasi yang hilang pada Citra Landsat 7 SLC-off karena adanya Stripping

Berbagai metode atau software mungkin memiliki kemampuan untuk mengisi data atau informasi yang hilang pada Citra Landsat 7 yang rusak karena adanya stripping. Namun beberapa software atau metode tersebut belum tentu valid untuk digunakan untuk scientific research, karena kebanyakan metode tersebut digunakan hanya untuk memperbaiki kualitas tampilan dan bukan memperbaiki kualitas informasi.

Namun, ada beberapa metode yang bisa digunakan untuk tujuan memperbaiki data yang hilang dengan menggunakan informasi pada Landsat images lainnya yang berbeda waktu. Artinya, data yang hilang tersebut dapat diisi dengan nilai yang ada pada citra lainnya yang memiliki waktu perekaman yang berdekatan, dan dilakukan pada setiap band yang ada atau band yang diperlukan. Hal ini mungkin dapat dilakukan pada berbagai images processing software, namun disini saya menggunakan Erdas 9.1.

Sebelum memulai, siapkan data Citra Landsat Utama yang akan ditambal, dan juga data satu atau lebih citra sebagai citra penambal yang memiliki waktu perekaman yang tidak jauh berbeda dengan waktu perekaman citra utama. Hasil penambalan tersebut dapat dilihat pada Gambar 1. 

Gambar 1. Tampilan citra sebelum (Kiri) dan sesudah ditambal (Kanan)

Adapun langkahnya adalah sebagai berikut:

1. Siapkan citra utama dan citra penambal (satu atau lebih). Disini saya menggunakan 1 citra Landsat utama (le71220652010205edc00_stack) dengan 1 buah citra penambal (le71220652010173sgs00_stack). Kedua citra tersebut sudah di-layerstack dengan jumlah dan urutan band yang sama pada kedua citra tersebut.
2. Buka ERDAS Imagine 9.1.
3. Buka citra *.tif file untuk citra utama dan citra penambal.

   

   Gambar 2. Citra Utama (Kiri) dan citra penambal (Kanan)

4. 
   

   
   

   Buka Model maker, lalu buat ekspresi yang menyatakan bahwa " Jika Image Utama > 0, maka gunakan data Image Utama, jika tidak, maka gunakan data Image 2". 
   
5. Ekspresi tersebut dapat ditulisakan dengan menggunakan model maker dengan sintaks berikut (contoh untuk Band 1) dalam Erdas Model Maker: EITHER $ n1_l71015033_03320070515_b10 IF ($ n1_l71015033_03320070515_b10> 0) Or $ n2_l71015033_03320070531_b50 OTHERWISE
   
6. Lakukan hal tersebut berulang-ulang pada band yang berbeda sehingga semua band tertambal dengan baik dan siap untuk di layerstack ulang.
7. Tambahan::Jika hasil tambalan masih kurang, silahkan tambahkan Citra Penambal kedua, setelah itu lakukan hal yang sama seperti pada point ke-4-5 diatas. Untuk diperhatikan, jika menggunakan 2 citra penambal, maka hasil penambalan pertama harus di setting "Temporary Raster Only". Untuk kemudian ditambal lagi dengan Citra ketiga. Begitu seterusnya. 

   

   Gambar 3. Tampilan model maker dengan dua citra penambal

8. Lalu band-band tersebut siap untuk di-layerstack.

Sumber referensi >> http://landsat.usgs.gov/sci_an.php