pch

1.4. Plot Titik pada Peta

Syarat: lakukan Project_01:Peta2 (klik disini).

Tujuan: Memplot titik pada peta secara random dengan koordinat terbatas.

Gunakan script peta2.R lalu simpan sebagai peta4.R pada direktori kerja.

Untuk menempatkan titik digunakan fungsi points() yang dapat diisi dengan parameter koordinat x, koordinat y, bentuk karakter titik, dan warnanya. Sebagai latihan, tuliskan perintah berikut sebelum baris plot fitur polygon:


# Plot titik secara random runif
points(x = runif(500, 126, 132) * 10000, y = runif(500, 103, 108) * 10000, pch=16, col='red')

Pada perintah diatas nilai koordinat x dan y dibuat random melalui fungsi runif dengan batasan x antara 126000 dan 132000, dan koordinat y antara 1030000 dan 1080000 (perhatikan nilai koordinat peta untuk x dan y). Nilai 500 (atau n) disini adalah jumlah observasi random. Nilai pch (plotting character) adalah antara 1 hingga 18.

Pada script lengkap, posisi layer points ini diletakkan dibawah layer polygon (artinya: ditulis duluan), kemudian nilai transparan polygon diubah menjadi 0.8 untuk memperlihatkan perbedaan antara warna titik yang tertutup polygon dan yang tidak tertutup.

Gambar yang dihasilkan adalah selalu berbeda karena nilai x dan y yang random, misalnya sebagai berikut:

Hasil penempatan titik secara random.

Hasil plot titik secara random.

Script lengkap dari pekerjaan ini adalah sebagai berikut:


# PROJECT-1: Mengambar Peta Pertama
# Tujuan: Plot titik
# File: peta4.R

# === MEMULAI R
rm(list=ls()) # Menghapus semua data
graphics.off() # Menutup semua grafik

# Menentukan direktori kerja
# Jangan gunakan spasi! gunakan "_" atau "." jika diperlukan
setwd("D:/_randmm/project01")
getwd() # Periksa isi direktori kerja
# ===

# Unggah paket GISTools
library(GISTools)

# Unggah data spasial Georgia
data("georgia")

# Pilih elemen pertama, simpan sebagai appling
appling <- georgia.polys[[1]]

# Tentukan batas luar plot
plot(appling, asp=1, type='n', xlab="Easting", ylab="Northing")

# Plot titik secara random runif
points(x = runif(500, 126, 132) * 10000, y = runif(500, 103, 108) * 10000, pch=16, col='red')

# Plot fitur yang telah dipilih dengan warna dan tingkat transparan
polygon(appling, col=rgb(0, 0.7, 0.5, 0.8))

Advertisements