Clip Multiple Rasters Using Polygons (Global Mapper)

-
Ada kasus, misalnya, anda diminta memotong (clip) data raster/citra sebanyak "n" file menggunakan grid polygon. Misalnya seperti gambar dibawah, saya mau membuat mosaik peta RBI hasil scanning, tapi saya tidak perlu informasi tepi dan bawah peta yang mengganggu. Jika bagian tidak penting ini tidak diclip, tentu hasil mosaik akan "awut - awut"-an.



Melakukan clip raster untuk kasus diatas sangat mudah jika data yang akan dimosaik cuma satu dua buah, tapi bagaimana jika ratusan??  anda tentu memerlukan sebuah tool yang bisa melakukannya secara otomatis.

Secara default fungsi clip multiple raster ini tidak tersedia di ArcGIS. Fungsionalitas batch processing juga tidak dapat digunakan untuk kasus seperti ini. Ada fungsi lain yang bisa, yaitu "iterate/loop raster" melalui model builder. Disini software akan menelusuri seluruh data raster yang ada dalam direktori, me-load ke software, diclip menggunakan polygon (identifikasi polygon menggunakan field tertentu, yang apesnya harus diiterasi juga). Dikarenakan dalam sebuah model geoprocessing iterasi tidak dapat dilakukan dua kali, maka untuk pekerjaan ini mau nggak mau harus dibuat dua model, seperti gambar di bawah.


Ribet, adakah cara lain???

setelah menelusuri dengan dipandu mbah google, ternyata prosedur serupa dapat dilakukan dengan relatif lebih mudah menggunakan Gombal (global) mapper melalui script tertentu yang berfungsi sebagai batch processing dan Looping. 

Here the Script. 

GLOBAL_MAPPER_SCRIPT VERSION=1.00
UNLOAD_ALL
DIR_LOOP_START DIRECTORY="E:\Workspace2014\Supervisi_PengelolaanDataSpasial\DATA RBI\PETA RBI SKALA 1 ; 25.000\Test" FILENAME_MASKS="*.tif" RECURSE_DIR=NO
IMPORT FILENAME="%FNAME_W_DIR%" CONTRAST_SHARED=YES CONTRAST_MODE=MIN_MAX
EXPORT_RASTER FILENAME="C:\GlobalMapper12\Result\%FNAME_WO_EXT%_ Cropped.tif " BG_TRANSPARENT=YES \
OVERWRITE_EXISTING=YES TYPE=GEOTIFF GEN_WORLD_FILE=YES \
POLYGON_CROP_FILE="E:\Workspace2014\Supervisi_PengelolaanDataSpasial\DATA RBI\PETA RBI SKALA 1 ; 25.000\test_shp.shp" \
POLYGON_CROP_USE_EACH=YES \
POLYGON_CROP_BBOX_ONLY=YES \
POLYGON_CROP_NAME_ATTR="ID_unique"
UNLOAD_ALL
DIR_LOOP_END

Ganti parameter DIRECTORY dengan folder anda menyimpan data raster yang akan dipotong.
Ganti parameter POLYGON_CROP_FILE dengan lokasi polygon pemotong
isikan POLYGON_CROP_NAME_ATTR dengan nama field yang akan menjadi dasar pemotongan data raster


Simpan script diatas di Notepad dan save dengan ekstensi *GMS.

Jalankan Script dari Menu File>Run Script

And, Here are the results


memang, masih ada bagian yang cukup menyebalkan, yaitu menghapus potongan - potongan kecil dari data raster yang tidak dipakai, tapi saya pikir ini lebih baik daripada harus clipping manual. 

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