 |
| Visualisasi peta dengan R |
Setelah kita mempraktikkan bagaimana cara memvisualisasikan grafik garis dengan R, kali ini kita akan beranjak pada visualisasi berikutnya, yaitu visualisasi peta (map) dengan menggunakan R. Saya rasa, visualisasi jenis ini saat ini banyak disukai khususnya para mereka yang memiliki penelitian di bidang geospasial atau sejenisnya. Saya sendiri belajar otodidak dari beragam sumber dan baru bisa dasarnya setelah setahun lamanya belajar step by step. Sebab, sangat sedikit referensi khususnya di mbah google yang menjelaskan secara terperinci dan step by step mengenai cara visualisasi peta dengan R, kebanyakan di youtube pun juga membahas dengan ArcGIS, QGIS, atau sejenisnya yang memang spesifik untuk pemetaan dan keperluan digitasi peta lapangan, khususnya penelitian perencaan tata ruang atau lingkungan.
Lantas, sebenarnya gimana sih caranya memvisualisasikan peta (map) dengan R? Kali ini coba kita fokus bareng-bareng untuk mempraktikkannya dengan memanfaatkan package rgdal. Adapun syarat utama kita bisa memvisualisasikan dengan package ini adalah tersedianya file berekstensi *dal di dalam file peta kita. Lalu, perlu juga data yang telah diisi dengan data-data terlebih dahulu sebagai bekal visualisasi gradasi atau sejenisnya.
Khusus untuk visualisasi peta gradasi, kita dapat menggunakan fungsi spplot yang tertanam di dalam package rgdal ini. Oke, sebelum mengeksekusinya, kita baiknya mengunduh file peta dulu ya teman-teman, peta yang kita gunakan kali ini adalah peta provinsi Jawa Timur yang terdapat di sini. Visualisasinya kita praktikan dengan beberapa code yang telah saya siapkan berikut hasilnya:
Code:
#install dan mengaktifkan package rgdal
install.packages("rgdal")
library(rgdal)
#Memanggil Peta Jatim, pahami betul cara import data
ptjatim=readOGR(dsn= "E:/peta",layer='jatim')
#Melihat variabel yang ada di dalam ptjatim
names(ptjatim)
Hasil:
[1] "KODE_KAB" "NAMA_KAB" "KODE_PROP" "NAMA_PROP" "KODE" "Y" "X1" "X2" "X3"
[10] "X4" "X5" "X6" "X7" "X8" "X9" "DATA" "Y1"
Code:
#Melihat Variabel Y yang akan kita visualisasikan
ptjatim$Y
Hasil:
[1] 1.1 3.0 1.9 3.2 3.1 4.4 8.0 6.8 4.9 4.3 3.0 1.7 3.3 3.0 5.3 6.0 8.5 1.9 2.7 3.0 5.0
[22] 6.5 1.6 1.0 1.2 2.7 2.1 2.7 6.7 3.6 3.7 9.3 3.1 8.3 16.2 4.6 4.6 8.5
Code:
#Visualisasi Dasar peta jatim
plot(ptjatim)
Hasil:
 |
| Visualisasi 1 |
Code:
#Visualisasi peta dengan warna
plot(ptjatim, col = "green")
Hasil:
 |
| Visualisasi 2 |
Code:
#Visualisasi gradasi variabel Y menurut wilayah
#membuat fungsi Warna dan Gradasi Visualisasi Peta
colfunc=colorRampPalette(c("#006633", "#66ff66", "#ccffcc"))
color=colfunc(16)
#Visualisasi Peta untuk variabel Y
spplot(ptjatim,"Y", col.regions=color, colorkey = T,
main="Sebaran Variabel Y menurut Kabupaten/Kota Jawa Timur")
Hasil:
 |
| Visualisasi 3 |