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

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

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

ArcGIS Tile Package Conversion to Mbtiles using tpkutils (Basemap series part 1)

-
Image
Imagine if all the published basemaps over the internet is not suit to your needs, what you can do?. The answer is simply, just make your own basemaps. And how to make them ? by convert the created map into raster tile package of course. From the created tile package, we can feed it off to map service provider like ArcGIS.com, Mapbox, etc, or just put it on your own server and webgis application/geoportal. Basemap can be authored using ArcGIS or QGIS (just to name the BIG TWO software these days). I am going to talk about ArcGIS first for this post. in ESRI ArcGIS environment, the tile package is called ArcGIS Tile Layer Package, and has *.tpk file extension. How to make them?, simply 1. Make your own map in ArcMap (symbolization, labelling, etc and save it on MXD. 2. Boot ArcToolbox, go to Data Management Tools > Packages > Create Map Tile Package 3. Set the tiling scheme (just use arcgis online/google tiling scheme), tiling format, level or details, tag, summary, ex
Read more