DEM Interferometri Menggunakan ALOS PALSAR Data

-
Setelah minggu lalu mencoba interferometri dengan Sentinel 1 dan ESA SNAP, minggu ini saya mencoba kembali percobaan interferometri untuk memperoleh DEM, tetapi menggunakan data ALOS PALSAR dan ENVI SARSCAPE. Percobaan saya lakukan dua kali, pertama menggunakan pasangan data ALOS PALSAR Level 1.1 FBS (Fine Beam Single) SLC, dan percobaan kedua menggunakan ALOS PALSAR Level 1.0 FBD (Fine Beam Double) RAW. Untuk pasangan data pertama selisih temporalnya cukup panjang, hampir 10 bulan, sementara untuk percobaan kedua selisih temporalnya lebih singkat, hanya 46 hari.

Tahapan pemrosesan yang dilakukan di Percobaan pertama meliputi:
1. Data Import (Sarscape>Import Data>SAR Spaceborne>ALOS PALSAR)
2. Simulated Orbit Correction menggunakan DEM(Sarscape>General Tools>Orbit Correction>automatic Orbit Correction). DEM yang dipakai adalah SRTM90
3. Interferogram Generation (Sarscape>interferometry>Phase Processing>Interferometry Workflow>Interferogram Generation)
4. Refinement dan Reflattening menggunakan 7 GCP (Sarscpae>Interferometry>Phase Processing>Interferometry Workflow>Refinement and Reflattening)
5. Phase to Height Conversion and Geocoding (Sarscape>interferometry>phase Processing>inferferometry Workflow>Phase to Height Conversion and Geocoding).

Orbit correction dilakukan karena data PALSAR Level 1.1 sudah tidak memiliki informasi orbit precise, jadi harus disimulasikan parameternya.

Hasil yang diperoleh DEM nya cukup bagus, hanya terdapat artifacts




Saya tidak tahu dari mana artifact ini muncul. Dugaan saya adalah karena orbit parameternya adalah hasil simulasi. Bukan data orbit precise asli dari satelit.

Kemudian percobaan kedua saya lakukan menggunakan data RAW Level 1.0 dengan temporal resolution yang lebih singkat. Saya mendapat informasi bahwa data orbit precise PALSAR tersimpan di metadata dari data Level 1.0, jadi menarik untuk dicoba.

Proses yang dilakukan untuk data level 1.0 sebagai berikut:

1. Focusing (Sarscape > Focusing > PALSAR > JAXA Palsar Single Pol), dengan polarisasi yang dipilih adalah HH, karena backscatternya lebih konsisten (koherensi lebih bagus).
2. Interferogram Generation (Sarscape>interferometry>Phase Processing>Interferometry Workflow>Interferogram Generation)
3. Refinement dan Reflattening menggunakan 7 GCP (Sarscpae>Interferometry>Phase Processing>Interferometry Workflow>Refinement and Reflattening)
4. Phase to Height Conversion and Geocoding (Sarscape>interferometry>phase Processing>inferferometry Workflow>Phase to Height Conversion and Geocoding). DEM yang dipakai adalah SRTM90.

Hasil yang diperoleh lebih baik (dari sisi artifacts yang muncul) dari percobaan pertama, hanya artifactnya masih ada, kali ini di sisi atas.


Minor artifact  pada percobaan kedua (mungkin) disebabkan oleh koherensi yang rendah antara dua pasangan citra.
Secara kualitas hasil, DEM pada percobaan pertama lebih baik dari percobaan kedua, walaupun perbedaan temporalnya lebih panjang. Hal ini dikarenakan perpendicular baselinenya lebih jauh, sehingga lebih sensitif dengan perbedaan elevasi.

Jadi. Kesimpulan sementara, Orbit file berpengaruh dan harus dipertimbangkan dalam interferometri. Untuk pengolahan PALSAR di Sarscape, lebih baik dimulai dari data RAW daripada data SLC. Selain itu, balancing antara temporal differences dan perpendicular baseline juga harus dipertimbangkan

Data ALOS Palsar dapat diunduh cuma cuma di ASF VERTEX.

Saya telah membuat contoh video tutorial pengolahan dan ekstraksi DEM menggunakan pendekatan interferometri, dengan contoh data Sentinel-1 SAR, dan menggunakan software ESA SNAP di bawah ini.


Comments

Post a Comment

Popular posts from this blog

Daftar Alamat dan Link Geoportal dan GIS Portal Kementerian, Lembaga, Badan dan Pemerintah Daerah di Indonesia

-
 Seiring dengan semakin populernya penggunaan data spasial di berbagai institusi, baik pemerintah maupun swasta di Indonesia, banyak lembaga pemerintah telah mengembangkan Geoportal sebagai sarana berbagi pakai data geospasial resmi pemerintah. Berbagai macam peta dasar dan peta tematik dengan berbagai format telah dipublikasikan secara resmi oleh geoportal yang dikelola masing - masing instansi.  Berikut ini adalah daftar alamat web geoportal milik instansi pemerintah indonesia, baik kementerian, badan, lembaga, maupun pemerintah daerah yang telah dibangun dan diresmikan. Daftar ini akan diperbarui sewaktu - waktu apabila tercatat ada lembaga pemerintah yang telah membangun dan menerbitkan geoportal sebagai salah satu sarana berbagi pakai data dan informasi geospasial.  Kementerian / Lembaga / Pemerintah Daerah Pengelola Alamat URL Portal INA SDI ">https://tanahair.indonesia.go.id/portal-web/# Kementerian Kelautan dan Perikanan ">http://satupeta.kkp.go.id/ ...
Read more

HAE (Height Above Ellipsoid) and MSL (Mean Sea Level) Conversion Using Pathfinder Office

-
Image
Di tepi pantai kok GPS-nya nilai ketinggiannya bukan nol ? mending bro yang saya alami malah nilai ketinggiannya negatif!, hmmmmm pasti GPSnya rusak, mari kita "lembiru" (lempar ganti baru)...... :) GNSS/GPS receiver secara default menggunakan acuan ketinggian (datum) berdasarkan ellipsoid WGS-84 (datum koordinat global), sehingga nilai ketinggian/elevasi yang diperoleh merupakan ketinggian diatas ellipsoid (HAE : Height Above Ellipsoid). Jika GPS yang kita gunakan menggunakan setting ini, maka tidak heran apabila ketika kita survei di pinggir pantai, ketinggiannya tidak nol atau malah mungkin negatif. Hal ini dikarenakan lokasi yang kita survei menurut WGS-84 ketinggiannya bukan nol. Ketinggian yang mendasarkan pada nilai nol diatas garis pantai rata rata disebut dengan MSL (Mean Sea Level). Kenapa bisa ada dua acuan ketinggian yang berbeda?. Hal ini dikarenakan bentuk bumi tidak sepenuhnya bulat, karena itu ketinggian nol di pantai ancol jakarta sejatinya tidak sama de...
Read more

TUTORIAL ORTHOREKTIFIKASI CITRA SATELIT RESOLUSI SEDANG (CITRA ASTER)

-
Image
Orthorektifikasi adalah proses koreksi geometrik citra satelit atau foto udara untuk memperbaiki kesalahan geometrik citra yang bersumber dari pengaruh topografi, geometri sensor dan kesalahan lainnya. Hasil dari orthorektifikasi adalah citra tegak (planar) yang mempunyai skala seragam di seluruh bagian citra. Orthorektifikasi sangat penting untuk dilakukan apabila citra akan digunakan untuk memetakan dan mengekstrak informasi dimensi, seperti lokasi, jarak, panjang, luasan, dan volume. Citra tegak merupakan citra) yang telah dikoreksi segala kesalahan geometriknya, sebagai akibat dari mekanisme perekaman citra. Kesalahan geometrik citra dapat berasal dari sumber internal satelit dan sensor (sensor miring/off nadir) ataupun sumber eksternal, yang dalam hal ini adalah topografi permukaan bumi. Perekaman off nadir dan perbedaan ketinggian berbagai obyek di permukaan bumi menyebabkan adanya kesalahan citra yang disebut relief displacement. Relief displacement sendiri dapat didefinisi...
Read more