Pada sesi kali ini, penulis akan membahas tentang tata cara pengolahan data vektor yaitu pemotongan data vektor pada software ArcGIS. Pemotongan data vektor pada software ArcGIS dapat menggunakan 2 cara yaitu: Select by Attributes dan Clip.

Select by Attributes dapat dilakukan dengan cara memasukkan algoritma pada attribute table sehingga memotong bagian yang diinginkan atau yang tidak diinginkan. Metode Clip dapat digunakan apabila anda ingin memotong data vektor berdasarkan suatu data vektor lainnya.

Berikut adalah tata cara melakukan Select by Attributes pada Software ArcGIS:

1. Siapkan data vektor yang anda ingin potong (dalam kasus ini penulis menggunakan data vektor fishnet yang penulis gunakan dalam post sebelumnya)
2. Pada toolbar Editor > Start Editing
3. Pada Table of Content klik kanan pada data vektor yang ingin dipotong > Open Attribute Table
4. Pada Attribute Table > Table Option > Select by Attributes
5. Masukkan algoritma yang anda inginkan (dalam kasus ini penulis menggunakan  “DEM”=0 or “DEM”=-9999 algoritma ini digunakan untuk menghilangkan data elevasi yang bernilai 0 dan -9999)
6. Klik Delete Selected
7. Berikut adalah hasil pemotongan data vektor menggunakan metode Select by Attributes:

 

Berikut adalah tatacara melakukan Clip pada software ArcGIS:

1. Siapkan data vektor dan data poligon yang akan anda gunakan
2. Pada Toolbar > Geoprocessing > Clip
3. Masukkan Input Features dengan data vektor anda > Clip features dengan shp garis pantai > Simpan pada lokasi yang anda inginkan
4. Berikut adalah hasil pemotongan data vektor menggunakan metode Clip:

 

Berikut adalah cara memotong data raster menggunakan data vektor dengan metode Extract by Mask pada software ArcGIS:

1. Siapkan data vektor yang telah terpotong
2. Pada ArcToolbox > 3D Analyst Tools > Raster Interpolation > Kriging
3. Masukkan Input point features dengan data vektor anda dan masukan z falue dengan data elevasi pada data vektor anda lalu OK
4. Berikut adalah hasil kriging yang telah dilakukan
5. Pada ArcToolbox > Spatial Analyst Tools > Extraction > Extract by Mask
6. Masukkan Input raster dengan hasil kriging yang telah anda lakukan dan input raster or feature mask data dengan data yang akan anda gunakan untuk memotong raster tersebut, lalu OK
7. Berikut adalah hasil pemotongan menggunakan Extract by Mask:

Pada sesi kali ini, penulis akan menunjukkan perbandingan hasil koreksi atmosfer yang telah dilakukan pada sesi-sesi sebelumya. Adapun metode yang telah dilakukan adalah 6SV, Apparent Reflectance, Dark Object Subtraction dan FLAASH Atmospheric Correction.

Berikut adalah perbandingan tampilan citra hasil koreksi atmosfer:

1. Landsat Level-2
2. 6SV
3. Apparent Reflectance
4. Dark Object Subtraction
5. FLAASH Atmospheric Correction

Berikut adalah perbandingan statistik hasil koreksi atmosfer:

1. Landsat Level-2
   

   
2. 6SV
   

   
3. Apparent Reflectance
   

   
4. Dark Object Subtraction
   

   
5. FLAASH Atmospheric Correction
   

   

Effectiveness of DOS (Dark-Object Subtraction) Method and Water Index Techniques to Map Wetlands in Rapidly Urbanising Megacity with Landsat 8 Data

Sean Gilmore, Ashty Saleem, Ashraf Dewan

PENDAHULUAN

Lahan basah pada permukaan bumi kurang lebih adalah 6-9% (Zedler dan Kercher,2005). Peran lahan antara lain adalah menjaga kualitas lingkungan dan termasuk penyimpanan karbon terestrial di dunia. Selain itu, lahan basah juga berperan sebagai oasis di area perkotaan yang sangan penting guna mengurangi suhu udara pemukaan. Perubahan distribusi lahan basah baik karena faktor alam maupun buatan dapat memberikan efek yang signifikan pada ekosistem. Pemetaan yang akurat dan statistik pada area yang tepat sangat penting guna mencegah dalam pengelolan lahan basah.

Data satelit penginderaan jauh banyak digunakan untuk menggambarkan lahan basah di seluruh dunia dengan berbagai teknik, termasuk klasifikasi per-piksel (supervised classification), metode semi otomatis (image segmentation) dan spectral water index. Metode koreksi atmosfer memerlukan informasi in-situ, relative scattering correction berguna untuk mengatasi masalah hamburan (scattering) pada citra. Metode Dark Object Substraction (DOS) adalah teknik berbasis image untuk menghilangkan komponen kabut yang disebabkan oleh hamburan aditif (additive scattering) dari data penginderaan jauh (Chavez Jr, 1988).

Pada penelitian ini menggunakan daerah Dhaka, Bangladesh. Penelitian ini juga berguna untuk menganalisis kesesuaian indeks berbasis air dalam menilai spasial lahan basah di wilayah perkotaan yang mengalami urbanisasi yang cepat. Penelitian ini diharapkan dapat berkontribusi secara signifikan terhadap basis pengetahuan yang ada di lokasi spasial lahan basah di area Rencana Pengembangan Metropolitan Dhaka (DMDP) yang merupakan unit perencanaan yang baru-baru ini dikembangkan oleh para pembuat kebijakan, menegakkan organisasi lokal untuk melestarikan ekosistem lahan basah yang tersisa.

 

PENGOLAHAN DATA

Citra yang digunakan berasal dari USGS Earth Explorer web service, dengan dua citra yang memiliki path-row 137-43 dan 137-44 sebagai area penelitian. Kedua citra tersebut dimosaik dan gambar yang dihasilkan dipotong sesuai area penelitian. Song dkk (2001) menjelaskan, pada citra Landsat, algoritma DOS mengasumsikan keberadaan ‘objek gelap’, yang memiliki piksel nol ke angka reflektansi sangat kecil sehingga nilai DN (digital number) minimum dalam histogram dianggap sebagai efek dari hamburan atmosfer dan dikurangi dari semua piksel pada citra, sehingga menciptakan “objek gelap” dengan nilai DN nol.

METODOLOGI

Dalam penelitian ini digunakan dua buah citra yaitu citra mentah dan citra yang telah tekoreksi DOS.

Normalised Difference Water Index (NDWI) dihitung menggunakan rumus yang diusulkan oleh McFeeters (1996):

NDWI = (Green Band – NIR Band) / (Green Band + NIR Band)

Dimana :
GREEN = Nilai reflektan kanal hijau 0,52 – 0,60 µm (Band 3)
NIR = Nilai reflektan kanal infra merah dekat 0,76 – 0,90 µm (Band 5)

Modified Normalized Difference Water Index (MNDWI) dihitung menggunakan rumus yang diusulkan oleh Xu (2006):

MNDWI = (Green Band – MIR Band) / (Green Band + MIR Band)

Dimana :
GREEN = Nilai reflektan kanal hijau 0,52 – 0,60 µm (Band 3)
MIR = Nilai reflektan kanal Inframerah gelombang menengah 1,55 – 1,75 µm (Band 6)

The Collective Indices classification dilakukan dengan menggunakan versi modifikasi dari formula yang dikembangkan Lu et al (2011). Pertama-tama, MNDWI digunakan karena lebih baik dalam mengekstraksi lahan basah dalam lingkungan yang didominasi perkotaan. Kedua, rumus tersebut dibalik dengan mengurangi NDVI dari MNDWI.

Collective Indices = MNDWI – NDVI

Dimana, NDVI dihitung menggunakan rumus yang diusulkan oleh Tucker (1979):

NDVI = (NIR Band – Red Band) / (NIR Band + Red Band)

Dimana :
NIR = Nilai reflektan kanal Inframerah gelombang dekat 0,76 – 0,90 µm (Band 5)
RED = Nilai reflektan kanal merah 0,63 – 0,69 µm (Band 4)

Automated Water Extraction Index (AWEI) dihitung menggunakan versi modifikasi dari formula non-shadow yang dikembangkan oleh Feyisa dkk (2014). Rumus yang telah dimodifikasi adalah:

AWEI = 4 x (Band 3 – Band 6) – (0.25 x Band 5 + 2.75 x Band 7)

Dimana :
BAND 3 = Nilai reflektan kanal hijau 0,53 – 0,59 µm (Green Band)
BAND 5 = Nilai reflektan kanal Inframerah gelombang dekat 0,85 – 0,88 μm (NIR Band)
BAND 6 = Nilai reflektan kanal  Inframerah gelombang pendek 1.57 – 1.65 µm (SWIR1 Band)
BAND 7 = Nilai reflektan kanal Inframerah gelombang pendek 2.11 – 2.29 µm (SWIR2 Band)

HASIL DAN ANALISIS

MNDWI adalah indeks terbaik dengan akurasi keseluruhan baku / DOS sebesar 97.28 % / 98% dan nilai kappa lebih besar dari 0.9 yang efektif membedakan lahan basah di pusat perkotaan dan juga di daerah rawa.

Metode collective indices terbaik kedua dengan akurasi keseluruhan baku / DOS 96.57 / 97.86% dan nilai kappa lebih besar dari 0.89 serta mampu mengekstrak mayoritas daerah lahan basah rawa, namun salah mengklasifikasi beberapa piksel dalam area perkotaan.

AWEI adalah indeks berkinerja terbaik ketiga dengan akurasi keseluruhan mentah / DOS dari 94.85 / 96.86% dan nilai kappa lebih besar dari 0.84 tetapi salah mengklasifikasikan beberapa daerah perkotaan sebagai lahan basah.

NDWI adalah indeks berkinerja terburuk dengan akurasi baku / DOS 92.7 % / 92.86% dan nilai kappa lebih besar dari 0.7 cukup buruk dalam memisahkan daerah rawa daerah perkotaan.

Lokasi spasial lahan basahs menggunakan teknik rasio band yang berbeda:

(a) NDWI diproduksi dari raw citra

(b) NDI yang diklasifikasikan menggunakan ambang optimal

(c) Klasifikasi NDWI dari citra DOS

(d) Klasifikasi DOS NDWI menggunakan optimal threshold

(e) MNDWI yang dihasilkan dari citra raw

(f) Klasifikasi MNDWI menggunakan optimal threshold

(g) MNDWI dihasilkan dari citra DOS

(h) Klasifikasi DOS MNDWI menggunakan optimal threshold

(i) Metode Collective indices yang dihasilkan dari citra raw

(j) Metode klasifikasi Collective indices menggunakan optimal threshold

(k) Metode Collective indices yang dihasilkan dari citra DOS

(l) Metode Klasifikasi Collective indices DOS yang menggunakan optimal threshold

(m) AWEI dihasilkan dari raw image

(n) Classified AWEI menggunakan optimal threshold

(o) AWEI yang dihasilkan dari citra DOS

(p) Collective Indices DOS yang diklasifikasikan menggunakan optimal threshold

KESIMPULAN

Dari penelitian ini dapat disimpulkan bahwa efek hamburan atmosfer harus dihilangkan sebelum analisis dilakukan pada data Landsat 8. MNDWI ditemukan menjadi indeks yang cocok dan secara akurat menentukan lokasi spasial lahan basah di dalam area DMDP.

Pada sesi kali ini, penulis akan menjelaskan cara melakukan koreksi atmosfer metode 6SV (Seconds Simulation of a Satellite Signal in the Solar Spectrum – Vector) pada citra landsat 8. Metode 6SV ini dapat dilakukan melalui website http://6s.ltdri.org/. Berikut adalah cara melakukan koreksi atmosfer menggunakan metode 6SV pada citra landsat 8:

1. Buka website http://6s.ltdri.org/ > Run 6SV
2. Submit > pada make your own atmospheric correction pilih Geometrical Condition
3. Step 1 > Geometrical Condition > User’s > Submit
4. Isi User’s Geometrical Condition sesuai dengan metadata citra anda, untuk bulan dan hari dapat dilihat pada Date Acquired dan Scene Center Time, Solar Zenithal adalah 90 – Sun Elevation, Solar Azimuthal adalah Sun Azimuthal, Sensor Zenithal adalah Roll Angle, dan Sensor Azimuthal dapat dilakukan pada website www.fcc.gov dengan mencari azimuth pada Corner UL Lat Long dan Corner LL Lat Long > submit
5. Berikut adalah data Roll angle, Sun azimuthal dan Sun Elevation
6. Berikut adalah Corner UL Lat Long dan Corner LL Lat Long
7. Berikut adalah website http://www.fcc.gov
8. Step 2 > Atmospherical Model > Atmospherical Profile = Dipilih sesuai daerah yang anda pilih (dalam kasus ini penulis menggunakan Tropical karena penulis menggunakan citra Lombok) > Aerosol Model = Dipilih sesuai daerah yang anda pilih (dalam kasus ini penulis menggunakan Maritime Model karena penulis menggunakan citra Lombok) > Submit
9. Pada maritime model > Visibility (km) > nilai Visibility dapat dilihat pada www.wunderground.com > Submit
10. Konversikan nilai visibility yang anda dapat dalam satuan km
11. Step 3 > Target Sensor Altitude > Target altitude = Altitude (nilai altitude dapat dilihat pada BPS wilayah anda) > Sensor altitude = Satellite Level > Submit
12. Step 4 > Spectral Condition > Spectral Condition = same filter function=1 > Band = 1st band of tm (landsat 5) (0.430-0.560) > Submit
13. Isikan wavelength inf dengan nilai minimum wavelength band yang anda inginkan (dalam kasus ini band 1 landsat 8 adalah 0.435) > isikan wavelength sup dengan nilai maksimum band yang anda inginkan (dalam kasus ini band 1 landsat 8 adalah 0.451) > Submit
14. Step 5 > Ground Reflectance > Ground Reflectance = homogeneous surface > Directional Effects = non directional effects > Submit
15. Specify surface reflectance = dipilih sesuai lokasi anda (dalam kasus ini  penulis menggunakan mean spectral value of green vegetation karena sebagian besar area yang penulis gunakan adalah vegetasi) > Submit > Submit
16. Step 6 > Signal > Atmospheric correction mode = Atmospheric correction with Lambertian assumption > Pilih satuan reflectance = [w/m2/str/mic] > tuliskan 0 > Submit
17. Step 7 > Result > tunggu hinggal tombol Output File muncul > Klik
18. Setelah muncul hasilnya scroll hingga kebawah dan anda akan menemukan atmospheric correction result > coefficients xa xb xc > Simpan nilai tersebut
19. Kembali pada tampilan result > ulangi Step 4 > Lakukan dengan memperhatikan nilai wavelength setiap band dari 1 – 8
20. Berikut adalah hasil Xa Xb Xc band 1 hingga band 8
21. Buka software Envi anda
22. Pada toolbox > band math
23. Hitunglah nilai Radiance masing-masing band dengan menggunakan rumus berikut
24. Nilai ML dan AL yang dapat anda lihat pada metadata citra anda
25. Pada band math > masukkan rumus (ML*float(B1)+(AL) >Add to list > OK
26. Pilih band 2 citra anda > Simpan > OK
27. Lakukan pada band 2, 3 dan 4
28. Pada toolbox > band math
29. Hitunglah nilai Y masing-masing band dengan menggunakan rumus berikut
30. Pada band math > masukkan rumus (Xa*float(B1)-Xb) > Add to list > OK
31. Pilih radiance band 2 citra anda > Simpan > OK
32. Lakukan pada radiance band 2, 3 dan 4
33. Pada toolbox > band math
34. Hitunglah nilai acr (reflektan-permukaan) masing-masing band dengan menggunakan rumus berikut
35. Pada band math > masukkan rumus (float(B1)/(1+(Xc*float(B1))) > Add to list > OK
36. Pilih Y band 2 citra anda > Simpan > OK
37. Lakukan pada Y band 2, 3 dan 4
38. Pada toolbox > Layer Stacking
39. Pada Layer Stacking Parameter > Import file = hasil acr band 4, 3 dan 2 > Simpan > OK
40. Berikut adalah hasil Layer Stacking metode 6SV
41. Pada toolbox > Compute Statistics
42. Pilih hasil Layer Stacking metode 6SV anda pada kolom select input file > OK
43. Berikut adalah hasil Compute Statistic metode 6SV
44. Karena pada nilai Min < 0 dan Max > 1 maka perlu dilakukan perhitungan pada band math
45. Pada toolbox > band math
46. Pada band math > masukkan rumus (b1 le 0)*0+(b1 ge 10000)*1+(b1 gt 0 and b1 lt 10000)*float(b1)/10000 > Add to list > OK
47. Pilih acr band 2 citra anda > Simpan OK
48. Lakukan pada acr band 2, 3 dan 4
49. Pada toolbox > Layer Stacking
50. Pada Layer Stacking Parameter > Import file = hasil acr terkoreksi band 4, 3 dan 2 > Simpan > OK
51. Berikut adalah hasil Layer Stacking metode 6SV terkoreksi
52. Pada toolbox > Compute Statistics
53. Pilih hasil Layer Stacking metode 6SV terkoreksi anda pada kolom select input file > OK
54. Berikut adalah hasil Compute Statistic metode 6SV terkoreksi

Pada sesi kali ini metode Apparent Reflectance menggunakan software ArcGIS. Berikut adalah tutorial metode apparent yang telah penulis lakukan :

1. Pada menu bar klik Add Data > citra Landsat 8 anda dalam format .MTL > Multi Spectral > Add
2. Klik kanan pada layer citra anda > Properties > Function
3. Klik kanan pada Stretch Function lalu remove
4. Pada Composite Band Function > Insert Function > Apparent Reflectance Function
5. Setelah jendela Raster Function Properties muncul > centang Albedo lalu OK > pada Layer Properties > Apply > OK
6. Hasilnya adalah seperti berikut
7. Klik kanan pada layer citra anda > Data > Export Data
8. Pada jendela Export Raster Data > Setting seperti gambar berikut > Save
9. Pada jendela Output Raster klik Yes > Hasilnya adalah sebagai berikut

Pada metode ini membutuhkan citra yang telah terkalibrasi radiometrik yang dapat dilihat pada post yang telah saya upload sebelumnya DOS. Setelah anda memiliki data radiance, selanjutnya saya akan memberikan tutorial koreksi atmosfer metode FLAASH.

Berikut adalah tutorial koreksi atmosfer metode FLAASH:

1. Pada menu toolbox cari FLAASH lalu pilih FLAASH Atmospheric Correction
2. Setelah menu FLAASH terbuka isikan Input Radiance Image dengan data radiance yang telah anda lakukan, lalu isikan lokasi penyimpanan anda hingga Rootname FLAASH file anda, lalu isikan Flight Date dan Flight Time sesuai metadata anda, pilih sensor type anda sesuai citra anda (dalam kasus ini penulis menggunakan Landsat-8 OLI), untuk memperoleh Ground Elevation (km) anda dapat membuka BPS daerah yang anda inginkan, klik Multispectral Settings
3. Pada jendela Multispectral Setting > Kaufman-Tanre Aerosol Retrieval > Assign Default Value Based on Retrieval Conditions pilih Over-Land Retrieval Standard (660:2100 nm) lalu OK > setelah muncul jendela FLAASH klik OK
4. Hasilnya adalah sebagai berikut
5. Untuk pengecekan pada toolbox cari Compute lalu pilih Compute Statistics
6. Pilih citra FLAASH anda lalu OK
7. Setelah muncul jendela Compute Statistics Parameters klik OK
8. Berikutlah hasil Compute Statistics tetapi hasilnya tidak sesuai karena nilai min tidak 0 dan maks tidak 1, maka perlu dilakukan perhitungan algoritma menggunakan band math
9. Pada toolbox cari Band lalu pilih Band Math
10. Masukkan rumus “(b1 le 0)*0+(b1 ge 10000)*1+(b1 gt 0 and b1 lt 10000)*float(b1)/10000” lalu OK
11. Pilih citra FLAASH anda > pilih lokasi penyimpanan anda > OK
12. Pada toolbox cari Compute lalu pilih Compute Statistics
13. Pilih citra hasil band math yang telah dilakukan anda lalu OK
14. Setelah muncul jendela Compute Statistics Parameters klik OK
15. Berikut adalah hasil Compute Statistics setelah dihitung menggunakan Band Math

Pada sesi kali ini penulis akan memanfaatkan data kontur yang telah kita buat pada post pembuatan kontur. Berikut adalah tutorial pembuatan Fishnet menggunakan software ArcGIS:

1. Munculkan citra DEM-SRTM dan Kontur yang telah dibuat
2. Pada Table of Content, klik kanan pada citra DEM-SRTM lalu pilih Properties > Symbology > atur Display Background Value : 0 dan ubah warnanya menjadi Hitam
3. Pada ArcToolbox pilih Data Management Tools > Sampling > Create Fishnet
4. Pada jendela Create Fishnet, pilih lokasi penyimpanan anda >pada kolom Template Extent (optional) pilih same as layer DEM.dat > isikan Number of Rows dan Number of Coloums sesuai yang anda inginkan (dalam kasus ini penulis menggunakan 50 pada keduanya), lalu OK
5. Berikut adalah hasil dari Fishnet
6. Pada ArcToolbox pilih Spatial Analyst Tools > Extraction > Extract Multi Values to Points
7. Pada jendela Extract Multi Values to Points, pilih hasil Fishnet_label anda untuk input point features > pilih DEM.dat untuk input raster lalu OK
8. Berikut adalah hasil Extract Multi Values to Points dapat dilihat pada Open Attribute Table pada Fishnet_label
9. Selanjutnya Start Editing dan hapus semua nilai ketinggian 0 untuk menghilangkan Fishnet diluar pulau lalu save editing
10. Berikut adalah hasil dari Editing pada Fishnet :

Berikut ini adalah tutorial pembuatan TIN menggunakan software ArcGIS :

1. Pada menu ArcToolbox pilih 3D Analyst Tools > Data Management > Tin > Create TIN
2. Pada jendela Create TIN, pilih lokasi penyimpanan anda > Pilih sistem koordinat yang anda gunakan (dalam kasus ini penulis menggunakan GCS_WGS_84 karena penulis ini menyamakan sistem koordinat dengan data DEM) > masukkan data kontur pada Input Feature Class (optional) lalu OK
3. berikut adalah hasil dari TIN menggunakan software ArcGIS :

Berikut ini adalah tutorial pembuatan Kriging menggunakan tools Geostatistical Analyst pada software ArcGIS :

1. Klik kiri pada Customize > Toolbars > Geostatistical Analyst
2. Letakkan tools Geostatistical Analyst sesuai keinginan anda > klik kiri pada Geostatistical Analyst > pilih Geostatistical Wizard
3. Pada Jendela Geostatistical Wizard lihat pada Geostatistical Methods > Kriging / CoKriging > isi data field dengan DEM.dat lalu next
4. Pada step 2 of 5 pilih Ordinary pada Kriging Type dan Prediction pada Output Surface Type lalu next
5. Pada step 3 of 5 lihat pada General > Variable diganti Covariance > lihat pada Model Nugget > Enable diganti True > Calculate Nugget > True > lihat pada Model #1 > Type diganti Spherical > Anisotropy diganti True > Calculate Partial Sill diganti True > lihat pada View Settings > Show search direction diganti True > Show all lines diganti True > next
6. Pada step 4 of 5 langsung saja klik next
7. Pada step 5 of 5 langsung saja klik next
8. Setelah mengisi kelima step tadi maka akan muncul Method Report yang berisikan tentang ringkasan dari pengaturan yang telah kita lakukan, apabila telah sesuai klik OK
9. Berikut adalah hasil dari Kriging menggunakan tools Geostatistical Analyst pada software ArcGIS :

Berikut adalah tutorial pembuatan Hillshade menggunakan software ArcGIS :

1. Pada menu ArcToolbox pilih 3D Analyst Tools > Raster Surface > Hillshade
2. Pada jendela Hillshade > isikan DEM.dat pada input raster > simpan sesuai keinginan anda (apabila terjadi kegagalan maka pindah lokasi penyimpanan anda) > OK
3. Berikut adalah hasil dari Hillshade menggunakan software ArcGIS :

Berikut adalah tutorial pembuatan Slope menggunakan software ArcGIS

1. Pada menu ArcToolbox pilih 3D Analyst Tools > Raster Surface > Slope
2. Pada jendela Slope > isikan DEM.dat pada input raster > simpan sesuai keinginan anda (apabila terjadi kegagalan maka pindah lokasi penyimpanan anda) > isikan PERCENT_RISE pada Output Measurement (optional) > OK
3. Berikut adalah hasil dari Slope menggunakan software ArcGIS :

Berikut adalah contoh pembuatan tampilan 3D menggunakan data Hillshade yang telah kita buat caranya adalah pada Table of Content > klik kanan pada layer Hillshade > Properties > Display > atur transparancy menjadi 60% >Apply > OK

Pastikan layer Hillshade diatas layer DEM, Slope, Kriging dan TIN > nyalakan hanya 2 kombinasi Layer (Hillshade + DEM, Hillshade + Slope, Hillshade + Kriging dan Hillshade + TIN). Sebagai hasilnya akan ditunjukkan pada puncak Gunung Rinjani tepatnya pada Danau Segara Anak seperti gambar-gambar dibawah ini :

Pada sesi kali ini akan dibahas tentang pembuatan kontur otomatis menggunakan DEM-SRTM yang telah di download pada USGS. DEM-SRTM tersebut memiliki resolusi sebesar 30 meter. Selain dapat menggunakan DEM-SRTM milik USGS kita juga dapat menggunakan DEMNAS yang dapat kita download pada DEMNAS-BIG. DEMNAS adalah data DEM yang dimiliki oleh BIG dengan beberapa sumber data mulai dari IFSAR yang memiliki resolusi 5 m, Terrasar X yang memiliki resolusi 5 m dan ALOS PALSAR yang memiliki resolusi 11.25 m. Untuk info lebih lengkap mengenai DEMNAS dapat anda akses pada DEMNAS-BIG. Sebelum kita mencoba membuat kontur dari DEM-SRTM, mari kita pelajari terlebih dahulu apa itu DEM, DTM, DSM.

• DEM (Digital Elevation Model)

DEM adalah data digital yang menggambarkan geometri dari bentuk permukaan bumi atau bagiannya yang terdiri dari himpunan titik-titik koordinat hasil sampling dari permukaan dengan algoritma yang mendefinisikan permukaan tersebut menggunakan himpunan koordinat (Tempfli, 1991).

DEM merupakan suatu sistem, model, metode, dan alat dalam mengumpulkan, prosessing, dan penyajian informasi medan. Susunan nilai-nilai digital yang mewakili distribusi spasial dari karakteristik medan, distribusi spasial di wakili oleh nilai sistem koordinat horisontal X Y dan karakteristik medan diwakili oleh ketinggian medan dalam sistem koordinat Z (Frederic J. Doyle, 1991).

DEM khususnya digunakan untuk menggambarkan relief medan. Gambaran model relief rupabumi tiga dimensi (3 dimensi yang menyerupai keadaan sebenarnya di dunia nyata (real world) divisualisaikan dengan bantuan teknologi komputer grafis dan teknologi virtual reality (Mogal, 1993).

DEM dapat diekstrak datanya untuk menampilkan perspektif 3D, kontur, kelas elevasi, profil, efek bayangan (hillshading), kemiringan lereng (slope), aspek permukaan, galian-timbunan, dll. DEM dapat digunakan untuk analisis medan (topografi), koreksi citra radar karena pengaruh layover pada medan perbukitan, dan visualisasi untuk penggambaran medan.

• DTM (Digital Terrain Model)

DTM adalah singkatan dari Digital Terrain Model atau bentuk digital dari terrain (permukaan tanah, tidak termasuk objek diatasnya) DTM menampilkan data yang lebih lengkap dari DEM. DTM digambarkan sebagai tiga representasi dimensi permukaan medan yang terdiri dari X,Y, Z koordinat disimpan dalam bentuk digital yang tidak hanya mencakup ketinggian dan elevasi unsur – unsur geografis lainnya dan fitur alami seperti sungai, jalur punggungan, dll DTM secara efektif DEM yang telah ditambah dengan unsur-unsur seperti breaklines dan pengamatan selain data asli untuk mengoreksi artefak yang dihasilkan dengan hanya menggunakan data asli. Dengan meningkatnya penggunaan komputer dalam rekayasa dan pengembangan cepat tiga dimensi grafis komputer DTM menjadi alat yang ampuh untuk sejumlah besar aplikasi di bumi dan ilmu teknik.

• DSM (Digital Surface Model)

DSM adalah digital surface model atau dapat diartikan sebagai model permukaan digital. DSM juga merupakan model elevasi yang menampilkan ketinggian permukaan, jika DTM hanya menampilkan ground (permukaan tanah tanpa apapun yang diatasnya) maka DSM menampilkan bentuk permukaan apapun yang ada seperti ketinggian pohon, bangunan dan objek apapun yang ada diatas tanah.

• SRTM (Shuttle Radar Topography Mission)

SRTM adalah sebuah citra yang memiliki data elevasi resolusi tinggi merepresentasikan topografi bumi dengan cakupan global (80% luasan dunia). Data SRTM adalah data elevasi muka bumi yang dihasilkan dari satelit yang diluncurkan NASA. Data ini dapat digunakan untuk melengkapi informasi ketinggian dari produk peta 2D, seperti kontur dan profil. Ketelitian ini dapat mencapai 15 m dan berguna untuk pemetaan skala menengah sampai skala tinggi (Lili somantri). SRTM 30 meter dapat kita download secara gratis dan dapat digunakan untuk berbagai pekerjaan yang berkaitan dengan ketinggian dengan sangat baik, misal membuat slope, kontur, TIN, aspec, flow accumulation, sungai, DAS dan lain-lain.

Setelah mempelajari dasar ilmu DEM, DEMNAS, DSM, DTM dan SRTM sekarang kita akan memulai untuk membuat kontur menggunakan data DEM-SRTM 30 menggunakan software ENVI dan dilanjutkan oleh ArcGIS dalam pembuatan konturnya.

Berikut adalah tahapan untuk melakukan percobaan ini :

1. Buka citra DEM-SRTM yang telah anda download pada software ENVI, (dalam kasus ini penulis menggunakan 2 citra DEM SRTM yang bersebelahan karena lokasi yang menjadi subjek penulis berada di 2 citra DEM-SRTM yang berbeda sehingga perlu disatukan terlebih dahulu kedua citra tersebut) apabila anda tidak mengalami masalah tersebut anda bisa langsung lewati hingga nomor 9
2. Pada menu Toolbox tuliskan Mosaicking, lalu pilih Pixel Based Mosaicking
3. Klik import > import file
4. Pilih 2 citra SRTM yang bersebelahan tersebut dengan cara tahan shift lalu klik kiri pada citra 1 lalu klik kiri lagi pada citra 2, lalu OK
5. Setelah muncul jendela berikut, isikan mosaic Xsize dengan 7202 dan mosaic Ysize 3601, lalu OK
6. Sesuaikan posisi citra dengan cara tahan klik kiri pada salah satu citra lalu ditarik ke kanan sehingga terbentuk sesuai yang diinginkan
7. Jika sudah sesuai, klik File > Apply > Simpan sesuai keinginan > OK
8. Hasilnya adalah sebagai berikut :
9. Masukkan vektor batas wilayah/ batas administrasi/ garis pantai (dalam kasus ini penulis menggunakan garis pantai) sehingga berpotongan dengan citra SRTM tersebut
10. Buat ROI dengan cara klik File > New > Region of Interest
11. Setelah muncul jendela berikut Klik File > Import
12. Pada bagian Select Input File pilih vektor yang anda gunakan lalu OK
13. Setelah terbentuk ROI pada citra SRTM anda, selanjutnya pada menu Toolbox tuliskan Subset Data dan pilih Subset Data from ROIs
14. Pilih citra SRTM anda (dalam kasus ini penulis menggunakan citra yang penulis telah gabungkan pada tahap-tahap awal) lalu OK
15. Pada jendela Spatial Subset via ROI Parameter, pilih ROI yang telah anda buat > ubah Mask pixel output of ROI menjadi Yes > Simpan sesuai keinginan > OK
16. Hasilnya adalah sebagai berikut
17. Buka software ArcGIS dan Inputkan hasil Subset ROI yang telah dilakukan sebelumnya, lalu pada menu ArcToolbox > 3D Analyst Tools > Raster Surface > Contour
18. Pada bagian input, isikan hasil Subset ROI anda, lalu pilih lokasi penyimpanan anda dan isikan Contour Interval sesuai kebutuhan anda (dalam kasus ini penulis melakukan percobaan dengan Contour Interval 15 dan 30) lalu OKCon

    

    

19. Berikut adalah kontur yang dihasilkan :
    

    Hasil Contour Interval 15m

    

    Hasil Contour Interval 30m

     

    

    Perbandingan 15 m dan 30 m

    

    Perbandingan Secara Detil

     

 

Berikut adalah tahapan untuk melakukan koreksi atmosfer menggunakan metode DOS pada citra landsat 8

1. Buka citra .MTL pada citra landsat anda
2. Pada menu toolbox, carilah Radiometric Calibration lalu klik 2 kali
3. Pilih citra multispektral anda lalu OK
4. Pilih Calibration Type anda dengan reflectance, lalu pilih Output interleave anda dengan BIL atau BIP dan simpan output pada lokasi yang anda inginkan, lalu OK
5. Lakukan step 2-5 dengan memilih radiance pada Calibration Typenya
6. Hasilnya akan seperti berikut

   

   

7. Pada menu toolbox, carilah Dark Subtraction lalu klik 2 kali
8. Pada kolom input file, pilih hasil radiance anda lalu OK
9. Pada Subtraction Method pilih band minimum, lalu simpan sesuai keinginan dan klik OK
10. Lakukan step 8-10 dengan memilih reflectance pada input filenya
11. Hasilnya akan seperti berikut

    

    

     

Dalam pembuatan peta skala 1 : 5000 menggunakan citra satelit, membutuhkan citra yang minimal resolusi spasialnya adalah 50cm. Data yang digunakan dalam pembuatan peta tersebut haruslah data yang posisinya berada di atas permukaan bumi. Sedangkan citra satelit sendiri adalah data yang posisinya berada di luar angkasa. Untuk mendapatkan data yang posisinya berada diatas permukaan bumi maka perlu dilakukan koreksi radiometrik.

Koreksi radiometrik adalah untuk memperbaiki nilai piksel supaya sesuai dengan yang seharusnya, biasanya mempertimbangkan faktor gangguan atmosfer sebagai sumber kesalahan utama. Efek atmosfer menyebabkan nilai pantulan obyek dipermukaan bumi yang terekam oleh sensor menjadi bukan merupakan nilai aslinya, tetapi menjadi lebih besar oleh karena adanya hamburan atau lebih kecil karena proses serapan. Kondisi kecerahan data citra selain pengaruh dari kondisi dan efek atmosfir, juga dipengaruhi oleh sudut sinar matahari dan sensitifitas sensor.

Gambar diatas menjelaskan bahwa garis putus-putus biru adalah data yang diperoleh oleh satelit (Top of Atmosphere), sedangkan garis hijau adalah data yang diinginkan (Bottom of Atmosphere). Untuk merubah data ToA menjadi data BoA maka diperlukan koreksi atmosfer sehingga efek atmosfer dapat dihilangkan.

Dalam proses koreksi radiometrik, terdapat 2 tahapan yaitu kalibrasi radiometrik dan koreksi atmosfer. Kalibrasi radiometrik adalah proses mengubah data pada citra yang (pada umumnya) disimpan dalam bentuk Digital Number (DN) menjadi radiance dan/atau reflectance. Sedangkan koreksi atmosfer adalah proses menghilangkan pengaruh atmosfer (molekul dan partikel) yang ikut serta menghamburkan sinyal sebelum direkam oleh sensor penginderaan jauh (pasif).

Dalam melakukan koreksi atmosfer sederhana, kita dapat menggunakan 2 metode yaitu metode DOS (Dark Object Subtraction) dan metode Radiative Transfer 6SV (Second Simulation of a Satellite Signal in the Solar Spectrum – Vector).