R.W. van Bemmelen adalah geologist paling popular di Indonesia. Karya terutamanya, kitab the Geology of Indonesia adalah bacaan wajib seluruh mahasiswa geologi di Indonesia, sejak almarhum Katili, hingga sekarang. Selain buku ini, van Bemmelen juga aktif melakukan pemetaan. Salah satu karyanya adalah Peta Geologi Lembar Bandung skala 1:100.000 yang diterbitkan pada tahun 1934.
Sejak memahami tentang teknik menggeoreferensi, saya ingin mengetahui apakah saya bisa menggeoreferensi peta van Bemmelen ini, sehingga kita tahu bagaimana sebaran batuan menurut van Bemmelen pada kondisi yang sekarang. Tidak terlalu penting sebenarnya, tapi kalau bisa ya kenapa tidak.
Ok pertama kita akan mengunduh peta geologi Lembar Bandung tahun 1934. Kita bisa mengunduh data ini di website koleksi digital Universitas Leiden https://digitalcollections.universiteitleiden.nl/ gunakan kata kunci “Bandoeng” dan cari di bagian “Maps (KITLV) (under construction)“. Kalau sulit, bisa pakai tautan ini.
Di tautan tersebut, Anda akan dapat 11 peta geologi jadul, keluaran tahun 1930an. Kita akan pilih nomor D E 14,3 yang merupakan kode Peta Geologi Lembar Bandung. Di bagian bawah peta, terdapat ikon untuk menguduh. Bentuk ikonnya seperti kardus pipih dengan mata panah mengarah ke kardus tersebut. Letaknya di samping ikon rantai, di atas tulisan In Collections. Silakan mengunduh Original Master dengan ukuran file 121 MB.
Untuk menggeoreferensi, kita akan menggunakan piranti lunak QGIS 3.10.5. Jika tidak punya software ini, maka bisa mengunduh dan menginstalnya di komputer masing-masing. Dalam tulisan ini, saya berasumsi bahwa pembaca sudah paham prinsip dasar menggeoreferensi.
Hal pertama yang harus dilakukan sebelum menggeoreferensi adalah mengetahui sistem koordinat. Bagaimana caranya? Yaitu dengan melihat pojok-pojok peta. Pada Peta Geologi Lembar Bandung ini kita melihat terdapat sistem koordinat, terutama pada bagian pojok kiri bawah di mana terdapat informasi sebagai berikut:
Di dalam peta ini kita dapat melihat bahwa Van Bemmelen sudah meletakkan sistem koordinat. Informasi di atas berarti bahwa titik pojok kiri bawah pada peta ini berjarak 0 derajat 40 menit ke arah timur dari Meridian Batavia, dan 7 derajat ke arah selatan dari ekuator. Meridian dari Batavia sendiri didefinisikan berjarak 106 derajat 48 menit 27.8 detik ke arah timur dari Greenwich. Informasi ini sangat penting, karena akan kita bandingkan dengan informasi yang kita miliki di QGIS.
Sekarang kita akan membuka QGIS. Ketika kita membuka QGIS, di bagian pojok kanan bawah, terdapat tulisan EPSG:XXXX yang merupakan sistem koordinat projek kita. Kita bisa klik gambar globe di samping tulisan EPSG untuk tahu sistem koordinat apa saja yang tersedia.
Selanjutnya, karena kita akan menggeoreferensi peta lama Indonesia, maka kita gunakan kata kunci “Batavia”. Maka akan ditemukan beberapa sistem koordinat (lihat baris Predefined Coordinate Reference Systems), yaitu sistem koordinat geografis dan terproyeksi (Geographic Coordinate Systems dan Projected Coordinate Systems).
Karena data kita menggunakan sistem koordinat geografis, maka kita cari di bagian Geographic Coordinate Systems. Perhatikan bahwa ada dua sistem, yaitu EPSG:4211 dan EPSG:4813. Apa bedanya? Kita akan tinjau.
Perbedaannya adalah EPSG:4211 berprime-meridien di Greenwich, artinya titik nol garis bujur berada di Greenwich. Sementara EPSG:4813 berprime-meridien di Batavia, artinya titik nol garis bujur pada sistem koordinat ini berada di Batavia. Namun terdapat keterangan bahwa Jakarta(Batavia) berjarak 106.807719 derajat dari Greenwich. Angka ini jika diubah menjadi derajat menit detik akan sama dengan 106 derajat 48 menit 27.8 detik. Artinya sistem koordinat ini akan sesuai dengan peta kita. Berikut adalah informasi geografis dari sistem koordinat EPSG:4813:
GEOGCRS["Batavia (Jakarta)",
DATUM["Batavia (Jakarta)",
ELLIPSOID["Bessel 1841",6377397.155,299.1528128,
LENGTHUNIT["metre",1]]],
PRIMEM["Jakarta",106.807719444444,
ANGLEUNIT["degree",0.0174532925199433]],
CS[ellipsoidal,2],
AXIS["geodetic latitude (Lat)",north,
ORDER[1],
ANGLEUNIT["degree",0.0174532925199433]],
AXIS["geodetic longitude (Lon)",east,
ORDER[2],
ANGLEUNIT["degree",0.0174532925199433]],
USAGE[
SCOPE["unknown"],
AREA["Indonesia - Bali, Java and western Sumatra onshore"],
BBOX[-8.91,95.16,5.97,115.77]],
ID["EPSG",4813]]
Oleh karena itu kita akan pilih sistem koordinat ini di peta kita.
Setelah itu kita akan masuk di bagian georeferensi.
Pada toolbar, klik Raster, kemudian pilih Georeferencer.
Di pojok kiri atas windows Georeferencer, pilih Open Raster. Pilih Raster yang ingin digeoreferensi. Dalam hal ini, raster yang sudah kita download, yaitu raster Peta Geologi Lembar Bandung. Hasilnya kurang lebih akan seperti gambar di bawah ini.
Di bawah tulisan settings, agak ke kanan sedikit, terdapat ikon berwarna kuning yang berbentuk seperti gerigi, yaitu ikon Transformation Settings. Kita akan mengatur Transformation pada pengaturan ini.
Pada pengaturan ini, pilih tipe transformasi menjadi Helmert. Metode Resampel menjadi Nearest neightbour dan Target SRS menjadi Project CRS, yaitu EPSG:4813. Atur output raster, sisanya tinggalkan sebagaimana default.
Pilih empat titik di pojok kiri atas, kanan atas, kanan bawah, dan kiri bawah sebagai titik ikat. Pastikan bahwa Anda menggunakan bagian dalam dari frame sebagai lokasi koordinat. Sebagai contoh saya menggeoreferensi pojok kiri atas peta, maka saya mengisi X/East sebagai 0 dan 40, artinya 0 derajat dan 40 menit. QGIS secara otomatis membacanya demikian. Sementara pada Y/North Anda akan mengisikan sebagai -6 dan 40 karena kita berada di selatan ekuator. Lakukan hal yang sama untuk sudut-sudut peta lainnya.
Setelah empat sudut terpilih sebagai GCP, maka kita bisa mulai menggeoreferensi dengan mengeklik ikon segitiga miring berwarna hijau, yang mirip dengan ikon play. Setelah itu kita bisa mengecek hasil georeferensi kita. Caranya dengan membuka file raster yang kita buat ke dalam QGIS.
Untuk mengecek apakah peta kita tergeoreferensi dengan benar, kita bisa bandingkan dengan peta OSM. Kita atur transparansi dari layer OSM untuk kita bandingkan bagaimana kedua layer ini. Saya pikir layer ini tergeoreferensi dengan baik sekali.
Demikian tutorial singkat menggeoreferensi Peta Geologi Van Bemmelen tahun 1934.
