Program Kegiatan RISET PKPP 2012 Kementerian Riset dan Teknologi (RISTEK)
MODUL BIMBINGAN TEKNIS PENGOLAHAN DATA CITRA SATELIT SOFTWARE OPEN SOURCE ILWIS 3.7 DAN QUANTUM GIS 1.5 Program Kegiatan RISET PKPP 2012 Kementerian Riset dan Teknologi (RISTEK) Koridor Non-Ekonomi / Mendukung Tupoksi Lembaga
PENGUATAN KAPASITAS DAERAH DALAM PEMANFAATAN DATA SYNTHETIC APERTURE RADAR (SAR) UNTUK PENGURANGAN RESIKO DAN MITIGASI BENCANA
Disusun Oleh: Fajar Yulianto, S.Si. Parwati, M.Sc. Dra. Any Zubaidah, M.Si. Ir. Dede Dirgahayu, M.Si. Kusumaning Ayu D.S, S.T. Junita Monika Pasaribu, S.Si. M. Nur Satrio, Amd. Tri Arso Sunugroho, Amd. Widya Putri K, Amd.
Pusat Pemanfaatan Penginderaan Jauh Kedeputian Penginderaan Jauh LEMBAGA PENERBANGAN DAN ANTARIKSA NASIONAL (LAPAN) Jalan LAPAN no. 70, Pekayon, Pasar Rebo, Jakarta 13710, Indonesia Tel./Fax : (021) 8722733
BIMTEK – Kab. Karawang, Prov. Jawa Barat, 2 – 6 Juli 2012
Page 1
Program Kegiatan RISET PKPP 2012 Kementerian Riset dan Teknologi (RISTEK)
BIMTEK – Kab. Karawang, Prov. Jawa Barat, 2 – 6 Juli 2012
Page 2
Program Kegiatan RISET PKPP 2012 Kementerian Riset dan Teknologi (RISTEK)
BIMTEK – Kab. Karawang, Prov. Jawa Barat, 2 – 6 Juli 2012
Page 3
Program Kegiatan RISET PKPP 2012 Kementerian Riset dan Teknologi (RISTEK)
BIMTEK – Kab. Karawang, Prov. Jawa Barat, 2 – 6 Juli 2012
Page 4
Program Kegiatan RISET PKPP 2012 Kementerian Riset dan Teknologi (RISTEK)
BIMTEK – Kab. Karawang, Prov. Jawa Barat, 2 – 6 Juli 2012
Page 5
Program Kegiatan RISET PKPP 2012 Kementerian Riset dan Teknologi (RISTEK)
BIMTEK – Kab. Karawang, Prov. Jawa Barat, 2 – 6 Juli 2012
Page 6
Program Kegiatan RISET PKPP 2012 Kementerian Riset dan Teknologi (RISTEK)
BIMTEK – Kab. Karawang, Prov. Jawa Barat, 2 – 6 Juli 2012
Page 7
Program Kegiatan RISET PKPP 2012 Kementerian Riset dan Teknologi (RISTEK)
BIMTEK – Kab. Karawang, Prov. Jawa Barat, 2 – 6 Juli 2012
Page 8
Program Kegiatan RISET PKPP 2012 Kementerian Riset dan Teknologi (RISTEK)
BIMTEK – Kab. Karawang, Prov. Jawa Barat, 2 – 6 Juli 2012
Page 9
Program Kegiatan RISET PKPP 2012 Kementerian Riset dan Teknologi (RISTEK)
BIMTEK – Kab. Karawang, Prov. Jawa Barat, 2 – 6 Juli 2012
Page 10
Program Kegiatan RISET PKPP 2012 Kementerian Riset dan Teknologi (RISTEK)
BIMTEK – Kab. Karawang, Prov. Jawa Barat, 2 – 6 Juli 2012
Page 11
Program Kegiatan RISET PKPP 2012 Kementerian Riset dan Teknologi (RISTEK)
BIMTEK – Kab. Karawang, Prov. Jawa Barat, 2 – 6 Juli 2012
Page 12
Program Kegiatan RISET PKPP 2012 Kementerian Riset dan Teknologi (RISTEK)
Import Data Citra Pada bagian ini dilakukan import data citra satelit dari format GeoTIFF ke format ILWIS, untuk selanjutnya digunakan dalam analisis dan proses pengolahan data lainnya. Berikut langkah dan proses yang dilakukan. File – Import – GDAL - Raster untuk import file format data ke format data ILWIS Pilih type data Tagged Image File Format (*.TIFF) karena data yang akan diimport memiliki format TIFF
Pilih nama file citra yang akan diimport dalam hal ini adalah data ALOS PALSAR dengan nama file IMG_HH_ALPSRP219617050_H15GUA_48S.tif untuk mulai mengimport data citra Pastikan lokasi output file berada pada direktori kerja yang telah dibuat. Bila belum, dapat dilakukan dengan mengkopi path-nya pada windows explorer dan mem-paste nya pada isian lokasi output, kemudian ditambah strip miring (“/”) dan nama filenya. Nama file sebagai output dalam format ILWIS dalam hal ini adalah IMG_HH_ALPSRP219617050_H15GUA_48S.mpr Bila sudah, dapat klik OK untuk mulai proses data Dengan menggunakan langkah dan proses yang sama dapat dicoba untuk import dan open data ALOS PALSAR scene ID berikut: • • •
IMG_HH_ALPSRP219617060_H15GUA_48S.tif IMG_HH_ALPSRP222097050_H15GUA_48S.tif IMG_HH_ALPSRP222097060_H15GUA_48S.tif
BIMTEK – Kab. Karawang, Prov. Jawa Barat, 2 – 6 Juli 2012
Page 13
Program Kegiatan RISET PKPP 2012 Kementerian Riset dan Teknologi (RISTEK)
Gambar berikut adalah hasil import IMG_HH_ALPSRP219617050_H15GUA_48S.mpr
dan
open
data
Sistem Georeference Pembuatan system georeference diperlukan untuk mempermudah tahap pengolahan data selanjutnya. Hal ini dilakukan dengan tujuan untuk menyamakan sistem koordinat data yang digunakan pada citra satelit. Berikut langkah dan proses yang dilakukan. Pastikan posisi navigator masih pada direktori kerja (nama_user) Klik File – Create – Coordinat System Pilih CoordSystem Projection dan isikan nama file: WGS1984_UTM48S kemudian klik OK Klik tombol Projection untuk menentukan system proyeksi pilih UTM Non-Aktifkan tanda cek pada pilihan Northern Hemisphere dan isikan 48 sebagai kode zona (wilayah Karawang dan sekitarnya termasuk dalam Zone 48S pada sistem koordinat UTM) Klik tombol Ellipsoid dan masukkan WGS84 sebagai definisi ellipsoid Klik tombol Datum dan masukkan WGS1984 sebagai definisi Datum Klik OK dan pastikan isian parameternya sama dengan contoh gambar berikut
BIMTEK – Kab. Karawang, Prov. Jawa Barat, 2 – 6 Juli 2012
Page 14
Program Kegiatan RISET PKPP 2012 Kementerian Riset dan Teknologi (RISTEK)
Klik GeoReference IMG_HH_ALPSRP219617050_H15GUA_48S Pilih Coordinate System yang telah dibuat dengan nama WGS_1984_UTM_49S
BIMTEK – Kab. Karawang, Prov. Jawa Barat, 2 – 6 Juli 2012
Page 15
Program Kegiatan RISET PKPP 2012 Kementerian Riset dan Teknologi (RISTEK)
Geo-Reference Geo-referensi (geo-reference) dalam ILWIS adalah suatu service object yang berfungsi untuk mendefinisikan hubungan antara baris dan kolom suatu data raster dengan posisinya pada sistem koordinat kartesian (X,Y). Data citra raster yang sudah di import pada tahap sebelumnya memiliki informasi koordinat sebatas cakupan scene citra. Geo-referensi dapat difungsikan untuk penggabungan 2 (dua) scene citra atau lebih dengan tujuan memberikan cakupan daerah/wilayah yang lebih luas (daerah yang memiliki cakupan scene citra lebih dari 1 (satu) scene). pengguna dapat menambahkan informasi tersebut dengan proses dan langkah sebagai berikut:
Dari menu bar, pilih File – Create – Georeference Masukkan nama GeoReference: GeoRef_Corners_UTM_49S Pilih GeoRef Corners untuk tipe georeferensinya Coordinate System : pilih WGS_1984_UTM_49S Pixel size: 6.25 (Tergantung resolusi spasial dari piksel citra) Min X,Y : 719165; 9377403 (Batas scene untuk wilayah Karawang) Max X,Y : 834074; 9239190
BIMTEK – Kab. Karawang, Prov. Jawa Barat, 2 – 6 Juli 2012
Page 16
Program Kegiatan RISET PKPP 2012 Kementerian Riset dan Teknologi (RISTEK)
Min X,Y : 719165; 9377403
Max X,Y : 834074; 9239190 IMG_HH_ALPSRP222097050_H15GUA_48S
IMG_HH_ALPSRP219617060_H15GUA_48S
Mozaik Scene Data Citra (Penggabungan Scene Data Citra) Proses penggabungan 2 (dua) scene citra di atas dapat dilakukan dengan pembatas file Geo-referensi (geo-reference) yang telah dibuat. Penggabungan tersebut dapat dilakukan dengan langkah dan tahapan sebagai berikut:
Dari menu file, klik kanan file data: IMG_HH_ALPSRP222097050_H15GUA_48S Pilih Spatial Reference Operation – Resample Isikan file IMG_HH_ALPSRP222097050_H15GUA_48S pada menu Raster Map Plotting point Resampling Method : Nearest Neighbour Output Raster Map: IMG_HH_ALPSRP222097050_H15GUA_48S_R Sistem GeoReference: GeoRef_Corners_UTM_49S Klik Show
BIMTEK – Kab. Karawang, Prov. Jawa Barat, 2 – 6 Juli 2012
Page 17
Program Kegiatan RISET PKPP 2012 Kementerian Riset dan Teknologi (RISTEK)
Dengan menggunakan langkah dan proses yang sama dapat dicoba untuk data ALOS PALSAR scene ID berikut: • • •
IMG_HH_ALPSRP219617060_H15GUA_48S.mpr IMG_HH_ALPSRP222097050_H15GUA_48S.mpr IMG_HH_ALPSRP222097060_H15GUA_48S.mpr
Proses penggabungan 2 (dua) scene citra dapat dilakukan dengan langkah: Operation – Raster Operation – Glue Maps
Pada menu Glue Raster Maps pilih plotting point 2 (untuk penggabungan 2 scene) dengan 1st Map : IMG_HH_ALPSRP222097050_H15GUA_48S_R dan 2nd Map : IMG_HH_ALPSRP219617060_H15GUA_48S_R Sistem GeoReference: GeoRef_Corners_UTM_49S Simpan Output Raster Map dengan nama: Pre_Flood_T1_R
BIMTEK – Kab. Karawang, Prov. Jawa Barat, 2 – 6 Juli 2012
Page 18
Program Kegiatan RISET PKPP 2012 Kementerian Riset dan Teknologi (RISTEK)
Dengan menggunakan langkah dan proses yang sama dapat digabungkan data berikut dengan scene ID berikut: • •
IMG_HH_ALPSRP222097050_H15GUA_48S.mpr IMG_HH_ALPSRP222097060_H15GUA_48S.mpr
Simpan Output Raster Map dengan nama: Post Flood_T2_R
Display Komposit Change Detection Banjir Deteksi banjir dilakukan dengan menggunakan metode change detection. Deteksi perubahan lahan akibat banjir tersebut ditunjukkan oleh data pada kondisi sebelum (Pre-Flood) dan sesudah banjir (PostFlood). Tahapan display komposit data sebelum dan sesudah banjir dilakukan sebagai berikut:
Pada menu toolbar pilih Operation – Image Processing – Color Composite. Muncul menu Color Composite pilih plotting point : standart, isikan file data pada Red Band : Pre_Flood_T1_R, Green Band : Post_Flood_T2_R, Blue Band : Post_Flood_T2_R. Pada Output Raster Map dapat diisikan dengan nama file Citra Banjir Karawang. Hasil pengolahan deteksi banjir dapat dilihat pada display gambar berikut. Deteksi pada daerah berwarna merah merupakan genangan banjir. Warna hitam menunjukan daerah genangan seperti sawah dll, yang merupakan BIMTEK – Kab. Karawang, Prov. Jawa Barat, 2 – 6 Juli 2012
Page 19
Program Kegiatan RISET PKPP 2012 Kementerian Riset dan Teknologi (RISTEK)
daerah yang sudah terdapat air sebelum banjir. Warna putih hingga abu – abu menunjukkan daerah kering yang tidak terkena banjir dan tidak mengalami perubahan.
Eksport Data Citra Pada bagian ini dilakukan eksport data citra satelit dari format ILWIS ke format *.LAN, untuk selanjutnya digunakan dalam analisis dan proses pengolahan data lainnya. Berikut langkah dan proses yang dilakukan. Klik kanan file data pada menu file yang akan di eksport, pilih Export – Method ILWIS – Format : Erdas.LAN. Output Filename: Citra Banjir Karawang.
BIMTEK – Kab. Karawang, Prov. Jawa Barat, 2 – 6 Juli 2012
Page 20
Program Kegiatan RISET PKPP 2012 Kementerian Riset dan Teknologi (RISTEK)
Quantum GIS Pada bagian ini software quantum GIS digunakan untuk melakukan overlay/tumpang susun dengan data vector (batas administrasi, jalan, sungai, dll), yang selanjutnya dapat digunakan untuk keperluan layout dan cetak peta. Pada menu toolbar klik Layer – Add Raster Layer. Atau dapat juga dilakukan dengan klik icon add Raster Layer. Open data file : Citra Banjir Karawanglan.lan Display map Citra Banjir Karawanglan.lan dapat terlihat dalam view gambar berikut.
Data vector seperti batas administrasi, jalan, sungai dapat dimasukan untuk menambahkan unsure informasi dalam kebutuhan analisis. Pada menu toolbar klik Layer – Add Vektor Layer. Atau dapat juga dilakukan dengan klik icon add Vektor Layer. Pada menu Add Vektor Layer yang tersedia, pilih source tipe plotting point File. Selanjutnya, klik menu Browser untuk memasukkan data vector: batas administrasi, jalan, sungai, dan landuse.
BIMTEK – Kab. Karawang, Prov. Jawa Barat, 2 – 6 Juli 2012
Page 21
Program Kegiatan RISET PKPP 2012 Kementerian Riset dan Teknologi (RISTEK)
Layout peta dapat dipersiapkan setelah komponen kelengkapan data yang dibutuhkan tersedia. Klik menu File – Composer Manager – Add – Show.
BIMTEK – Kab. Karawang, Prov. Jawa Barat, 2 – 6 Juli 2012
Page 22
Program Kegiatan RISET PKPP 2012 Kementerian Riset dan Teknologi (RISTEK)
BIMTEK – Kab. Karawang, Prov. Jawa Barat, 2 – 6 Juli 2012
Page 23
Program Kegiatan RISET PKPP 2012 Kementerian Riset dan Teknologi (RISTEK)
BIMTEK – Kab. Karawang, Prov. Jawa Barat, 2 – 6 Juli 2012
Page 24
Program Kegiatan RISET PKPP 2012 Kementerian Riset dan Teknologi (RISTEK)
BIMTEK – Kab. Karawang, Prov. Jawa Barat, 2 – 6 Juli 2012
Page 25
Program Kegiatan RISET PKPP 2012 Kementerian Riset dan Teknologi (RISTEK)
BIMTEK – Kab. Karawang, Prov. Jawa Barat, 2 – 6 Juli 2012
Page 26
Program Kegiatan RISET PKPP 2012 Kementerian Riset dan Teknologi (RISTEK)
BIMTEK – Kab. Karawang, Prov. Jawa Barat, 2 – 6 Juli 2012
Page 27
Program Kegiatan RISET PKPP 2012 Kementerian Riset dan Teknologi (RISTEK)
BIMTEK – Kab. Karawang, Prov. Jawa Barat, 2 – 6 Juli 2012
Page 28
Program Kegiatan RISET PKPP 2012 Kementerian Riset dan Teknologi (RISTEK)
BIMTEK – Kab. Karawang, Prov. Jawa Barat, 2 – 6 Juli 2012
Page 29
Program Kegiatan RISET PKPP 2012 Kementerian Riset dan Teknologi (RISTEK)
BIMTEK – Kab. Karawang, Prov. Jawa Barat, 2 – 6 Juli 2012
Page 30
Program Kegiatan RISET PKPP 2012 Kementerian Riset dan Teknologi (RISTEK)
MANUAL PENGOLAHAN DATA QMORPH Per-Jam VERSI GRADS 1. Download data di ftp://ftp.cpc.ncep.noaa.gov/precip/qmorph/30min_8km/
Datanya dalam format winzip/ WinRAR Extract save di c:/QMORPH/
(QMORPH_8km_20120324.Z)
2. Buka explorer dan Open file qmorph_00z.ctl dari jam 00z hingga 23z
Ubah sesuai tanggal
Ubah sesuai tanggal
BIMTEK – Kab. Karawang, Prov. Jawa Barat, 2 – 6 Juli 2012
Page 31
Program Kegiatan RISET PKPP 2012 Kementerian Riset dan Teknologi (RISTEK)
3. Buka Window Explorer, Open file plot_qmorph_00z.gs Edit didalamnya dengan tanggal yang baru. .
Ubah sesuai tanggal
4. Buka software grads - GradsNC/GradsC
Tekan enter
BIMTEK – Kab. Karawang, Prov. Jawa Barat, 2 – 6 Juli 2012
Page 32
Program Kegiatan RISET PKPP 2012 Kementerian Riset dan Teknologi (RISTEK)
5. Kemudian di run dengan mengetik perintah: run c:/qmorph/plot_qmorph_00.gs
Tekan enter
6. Lalu enter, maka akan muncul jendela/window baru sbb:
BIMTEK – Kab. Karawang, Prov. Jawa Barat, 2 – 6 Juli 2012
Page 33
Program Kegiatan RISET PKPP 2012 Kementerian Riset dan Teknologi (RISTEK)
7. Save hasilnya dengan command: ! gxps –c
Ketik dan Enter
8. Enter input file: c:/qmorph/tem
Ketik dan Enter
BIMTEK – Kab. Karawang, Prov. Jawa Barat, 2 – 6 Juli 2012
Page 34
Program Kegiatan RISET PKPP 2012 Kementerian Riset dan Teknologi (RISTEK)
9. Enter ouput file: c:/qmorph/QMORPH_8km_2012032400.ps
Ketik dan Enter
Tunggu sebentar sampai proses selesai lalu ulangi kedata berikutnya.
10. Open hasilnya di adobe photoshop, ganti image sizenya : Width : 100 percent ganti 75 percent Height: 100 percent ganti 75 percent Pilih Resolution: 72 Mode: RGB
BIMTEK – Kab. Karawang, Prov. Jawa Barat, 2 – 6 Juli 2012
Page 35
Program Kegiatan RISET PKPP 2012 Kementerian Riset dan Teknologi (RISTEK)
11. Rotate imagenya: Klik Image Rotate canvas 90° CW
12. Save hasilnya ke format JPG, pilih low quality. 13. Kerjakan untuk data per-jam yang lainnya hingga selesai 24 data jam-an.
BIMTEK – Kab. Karawang, Prov. Jawa Barat, 2 – 6 Juli 2012
Page 36
Program Kegiatan RISET PKPP 2012 Kementerian Riset dan Teknologi (RISTEK)
MANUAL PENGOLAHAN DATA TRMM VERSI GRADS 1. Download data di server ftp://trmmopen.gsfc.nasa.gov/pub/merged/mergeIRMicro/
Datanya dalam format winzip/ WinRAR Extract save di c:/TRMM/
(3B42RT.20120324.6.bin.gz)
2. Buka explorer dan Open file 3B42rt_00z.ctl dari jam 00z hingga 21z Ganti sesuai tanggal
Ganti sesuai tanggal
BIMTEK – Kab. Karawang, Prov. Jawa Barat, 2 – 6 Juli 2012
Page 37
Program Kegiatan RISET PKPP 2012 Kementerian Riset dan Teknologi (RISTEK)
3. Buka Window Explorer, Open file plot_trmm_00.gs Edit didalamnya dengan tanggal yang baru.
Ganti sesuai tanggal
4. Buka software grads - GradsNC/GradsC
Tekan enter
BIMTEK – Kab. Karawang, Prov. Jawa Barat, 2 – 6 Juli 2012
Page 38
Program Kegiatan RISET PKPP 2012 Kementerian Riset dan Teknologi (RISTEK)
5. Kemudian di run dengan mengetik perintah: run c:/trmm/plot_trmm_00.gs
Tulis utk running data lalu enter
6. Lalu enter, maka akan muncul jendela/window baru sbb:
BIMTEK – Kab. Karawang, Prov. Jawa Barat, 2 – 6 Juli 2012
Page 39
Program Kegiatan RISET PKPP 2012 Kementerian Riset dan Teknologi (RISTEK)
7. Kemudian Ketikan perintah: !gxps(spasi)-c Lalu enter kembali
Ketik dan enter
8. Enter input file: c:/trmm/tem
Ketik dan enter
BIMTEK – Kab. Karawang, Prov. Jawa Barat, 2 – 6 Juli 2012
Page 40
Program Kegiatan RISET PKPP 2012 Kementerian Riset dan Teknologi (RISTEK)
9. Enter ouput file: c:/trmm/3B42RT.2012032400.ps
Ketik dan enter
Tunggu sebentar sampai proses selesai lalu ulangi kedata berikutnya.
10. Open hasilnya di adobe photoshop, ganti image sizenya : Width : 100 percent ganti 75 percent Height: 100 percent ganti 75 percent Pilih Resolution: 72 Mode: RGB
Ketik sesuaikan dgn tanggal data dan enter
BIMTEK – Kab. Karawang, Prov. Jawa Barat, 2 – 6 Juli 2012
Page 41
Program Kegiatan RISET PKPP 2012 Kementerian Riset dan Teknologi (RISTEK)
11. Rotate imagenya: Klik Image Rotate canvas 90° CW
12. Save hasilnya ke format JPG, pilih low quality 13. Kerjakan untuk data per-jam yang lainnya hingga selesai 8 data jam-an.
BIMTEK – Kab. Karawang, Prov. Jawa Barat, 2 – 6 Juli 2012
Page 42