3D GIS Visualization and Application using R Program

Membuat tampilan 3 dimensi data GIS seperti ketinggian, slope, keterjalan tapak dan effek hillshade menggunakan Program R

Morfologi rupa bumi dapat digambarkan dalam bentuk digital elevation model (DEM). Data ini diperoleh dari SRTM atau Shuttle Radar Topography Mission yang merupakan salah satu produk hasil kerjasama antara NASA, the National Geospatial-Intelligence Agency, and the German and Italian Space Agencies. Dengan adanya data ketinggian ini, kita dapat menghitung atau konversi data elevasi menjadi data kontur, kemiringan lereng (slope), arah aliran sungai (flow direction), beberapa indeks keterjalan tapak seperti roughness, Terrain Ruggedness Index (TRI), hingga membuat efek tampilan hillshade untuk menggambarkan keadaan rupa bumi pada peta.

Informasi tersebut diatas sangat penting untuk berbagai aplikasi khususnya di bidang ekologi, misalnya untuk melihat hubungan antara spesies dengan lingkungannya. Ahli ekologi berpendapat bahwa sebaran organisme sangat dipengaruhi oleh suhu dan curah hujan. Sementara suhu dan curah hujan sangat dipengaruhi oleh bentuk rupa bumi dan posisinya terhadap garis khatilistiwa. Artinya sebaran spesies juga dipengaruhi oleh ketinggian, keterjalan tapak, kemiringan lereng, jarak dari sungai, dan sebagainya. Parameter-parameter ini biasanya digunakan dalam pemodelan distribusi spesies.

Pada tutorial sebelumnya, kita telah mebahas bagaimana cara mendapatkan data DEM ini baik secara langsung menggunakan program R (lihat tutorial) atau dari website USGS (lihat tutorial). Selain itu, kita juga sudah membahas bagaimana membuat peta kontur dari data DEM (lihat tutorial). Kali ini kita akan membahas bagaimana membuat peta kemiringan lahan (slope), keterjalan tapak (roughness/terrain ruggedness) dan membuat efek hillshade. Dengan menggunakan Program R, kita hanya perlu membuat satu baris perintah (script) sederhana untuk satu fungsi, jadi cukup ringkas dan cepat. Anda dapat mengunduh data yang akan digunakan pada tutorial ini (download data), atau jika terhubung dengan internet, Anda bisa langsung mengunduhnya melalui program R (lihat tutorial). Peta 3D yang akan dihasilkan dapat dilihat pada gambar berikut.

Gambar 1. Tampilan data ketinggian, slope, keterjalan, dan arah muka bumi menggunkan aplikasi R

Gambar 2. Tampilan hillshade pada data ketinggian di sekitar TNGP

Download Global Administrative, Topography and Bioclimatic Database Through R Program

Mengunduh data administrasi, topografi dan bio-iklim langsung melalui Program R

Salah satu keunggulan menggunakan program R adalah kita bisa langsung menunduh data dengan menggunakan satu atau dua baris perintah di program R tanpa harus membuka portal website tertentu untuk mencari dan mengunduhnya. Keunggulan ini dimiliki oleh pustaka “raster” yang sudah sering dipakai pada tutorial-tutorial sebelumnya. Data-data yang bisa diunduh langsung seperti batas administrasi suatu daerah/negara (Global Administrative Database-GADM)pada beberapa level, data ketinggian (Digital Elevation Model) resolusi 90m, data iklim global seperti suhu tahunan rata-rata/minimum/maksimum, curah hujan, dan data bioclime yang biasanya dipakai untuk analisis pemodelan distribusi spesies (Species Distribution Model).

Pada tutorial kali ini kita akan mencoba mengunduh langsung data-data tersebut yang kemudian bisa dipakai langsung untuk kebutuhan analisis data. Secara bertahap, kita akan memulainya dengan menyiapkan folder dimana output akan disimpan secara otomatis. Setelah itu memanggil data wilayah yang akan dijadikan acuan (area of interest) untuk mengunduh data. Lalu memanggil pustaka raster yang menyediakan perintah yang dibutuhkan untuk mengunduh data-data tersebut.

Data yang akan digunakan dalam tutorial ini dapat diunduh disini. Salah satu data yang akan dihasilkan dapat dilihat pada Gambar 1.

Gambar 1. Peta Bio-Iklim 1 yang merupakan data suhu tahunan rata-rata

## Mengaktifkan folder kerja, dimana semua hasil akan disimpan
setwd("D:/LUBIS_PRIVATE_DATA/R_GIS_PROJECT/04_Hasil")

Create distance layers using R program

Membuat peta jarak dari lokasi menggunakan software R

Analisis jarak (distance) sangat banyak manfaatnya dan telah banyak digunakan dalam studi atau penelitian dibidang ekologi maupun bidang lainnya. Misalnya pada penelitian mengenai tingkat deforestasi di suatu tempat, biasanya jarak dari jalan ataupun jarak dari pemukiman berkolerasi positif terhadap tingkat deforestasi. Artinya, hutan yang berada dekat dengan jalan atau perkampungan biasanya memiliki tingkat deforestasi yang lebih tinggi dibanding dengan hutan yang berada jauh dari jalan atau pemukiman.

Jadi pada tutorial kali ini, kita akan belajar melakukan analisis jarak dari suatu lokasi baik yang berbentuk titik, garis ataupun area. Pada data yang memiliki format raster, analisis jarak dapat dilakukan dengan menggunakan fungsi distance yang ada pada pustaka raster. Pada pustaka yang sama, kita juga bisa melakukan analisis jarak dari data vektor yang berbentuk titik tanpa harus merubah format ke raster terlebih dahulu dengan menggunakan perintah distanceFromPoints.

Analisis jarak berikutnya adalah membuat matriks jarak antar sepasang titik dengan menggunakan perintah pointDistance, hasil luarannya berupa matriks data, bukan peta. Hal ini biasanya digunakan untuk menghitung jarak rata-rata antar lokasi, misalnya jarak rata-rata antar lokasi survey atau transek. Data vektor atau raster yang digunakan untuk analisis jarak tidak harus memiliki sistem proyeksi UTM, bisa juga dilakukan pada data yang memiliki sistem kordinat geografis.

Kita akan memulai tutorial ini dengan menggunakan fungsi distance dari data vektor jalan yang sudah diubah menjadi raster pada tutorial sebelumnya. Setelah itu, kita akan membuat peta jarak antar desa tanpa harus mengkonversi data titik desa dalam format vektor ke format raster. Yang terakhir, kita akan membuat matriks jarak antar titik lokasi penelitian di TNGP. Data yang akan digunakan dalam tutorial ini dapat diunduh disini.. Hasil dari analisis jarak dapat dilihat pada gambar dibawah.

Gambar 1. Hasil analisis jarak dari batas TNGP dalam bentuk raster

## Mengaktifkan folder kerja, dimana semua hasil luaran akan disimpan
setwd("D:/LUBIS_PRIVATE_DATA/R_GIS_PROJECT/04_Hasil")

Raster to/from vector conversion in R

Mengubah format data raster dari dan ke vektor

Data spasial dalam bentuk vektor dapat diubah menjadi format raster, begitu juga sebaliknya. Seringkali kita dihadapkan pada situasi dimana bekerja dengan format raster lebih mudah dan lebih ringan dibanding bekerja dengan format vektor. Namun penting untuk diperhatikan bahwa perubahan atau transformasi data ini dapat menyebabkan berkurangnya akurasi, dan ini sangat bergantung dengan ukuran piksel yang diinginkan.

Tutorial kali ini kita akan mengubah format data vektor menjadi data raster dan sebaliknya menggunakan program R. Pada dasarnya, proses rasterisasi ini membutuhkan dua input data, yakni data vektor yang akan diubah menjadi raster dan template raster tujuan yang telah memiliki resolusi atau dimensi geografis yang dinginkan. Template raster ini berfungsi sebagai wadah tempat dimana nilai atau informasi yang ada di vektor akan disimpan. Proses rasterasi membutuhkan memori yang cukup banyak, tergantung jenis data vektor yang akan dikonversi. Sebagai contoh, mengubah data vektor seperti jalan yang berkelok-kelok atau sungai yang bercabang-cabang akan membutuhkan memori yang lebih besar dibandingkan dengan mengubah data lokasi atau titik desa dalam bentuk titik.

Ringkasnya, kita akan melakukan beberapa tahapan dimulai dari menentukan folder aktif, memanggil data vektor, memanggil data raster yang akan dijadikan sebagai template/wadah dari raster tujuan, kemudian mengubah kembali dari raster ke vektor, lalu membuat layout dan menyimpan data yang telah dikonversi. Tutorial sebelumnya telah membahas beberapa rangkaian proses ini. Salah satu hasi layout dari tutorial ini terlihat seperti gambar di bawah 1 di bawah.

Data yang digunakan dalam tutorial ini dapat diunduh disini.

Gambar 1. Hasil konversi data (DEM) dan kawasan TNGP dalam format raster menjadi vektor (kontur dan polygon)

# Mengaktifkan folder kerja, dimana semua hasil luaran akan disimpan
setwd("D:/LUBIS_PRIVATE_DATA/R_GIS_PROJECT/04_Hasil")