Peta koroplet
Sekarang setelah Anda memahami cara menggambar poligon, mari tampilkan poligon Anda di peta. Ingat, Anda mengganti pemanggilan ggplot() dengan ggmap() dan pemanggilan ggplot() asli dipindahkan ke argumen base_layer(), lalu Anda menambahkan layer poligon seperti biasa:
ggmap(corvallis_map_bw,
base_layer = ggplot(ward_sales,
aes(lon, lat))) +
geom_polygon(aes(group = group, fill = ward))
Cobalah di konsol sekarang!
Waduh, ada yang tidak beres. Ward 1, 3, dan 8 terlihat bergerigi dan salah. Apa yang terjadi? Sebagian batas ward berada di luar batas peta. Karena pengaturan bawaan di ggmap(), data apa pun yang berada di luar peta akan dibuang sebelum pemetaan, sehingga beberapa batas poligon terhapus dan ketika titik yang tersisa dihubungkan, bentuknya jadi keliru.
Jangan khawatir, ada solusinya: ggmap() menyediakan beberapa argumen untuk mengendalikan perilaku ini. Argumen extent = "normal" bersama dengan maprange = FALSE memaksa plot menggunakan rentang data alih-alih rentang peta untuk menentukan batas pemetaan.
Latihan ini adalah bagian dari kursus
Memvisualisasikan Data Geospasial di R
Latihan interaktif praktis
Cobalah latihan ini dengan menyelesaikan kode contoh berikut.
# Fix the polygon cropping
ggmap(corvallis_map_bw,
base_layer = ggplot(ward_sales, aes(lon, lat))) +
geom_polygon(aes(group = group, fill = ward))