Beberapa Package R yang Digunakan untuk Membuat Peta
library(maps)
library(ggplot2)
library(mapdata)
library(ggthemes)
library(ggrepel)
library(raster)
library(plyr)
library(dplyr)
Membuat Peta Dunia dengan Fungsi map_data()
peta_dunia <- map_data("world")
peta_dunia
gambar1 <- ggplot() +
geom_polygon(data = peta_dunia, aes(x=long, y = lat, group = group),
fill = "gray85", color = "gray80")
print(gambar1)
Membuat Peta Indonesia
peta_indonesia <- gambar1 + xlim(94.5,141.5) + ylim(-11.5,6.5)
print(peta_indonesia)
Menggunakan fungsi getData(), Mendapatkan Nama-Nama Provinsi
peta_indonesia_level_1 <- getData('GADM', country = 'IDN', level = 1)
## Warning in getData("GADM", country = "IDN", level = 1): getData will be removed in a future version of raster
## . Please use the geodata package instead
nama_provinsi = peta_indonesia_level_1$NAME_1
nama_provinsi
Melakukan Pengelompokkan Wilayah Provinsi pada Peta Indonesia
label_angka <- seq(1, 33)
NAME_1 = peta_indonesia_level_1$NAME_1
gabung_label_nama <- data.frame(NAME_1, label_angka)
indo1 = peta_indonesia_level_1
indo1@data$id <- rownames(indo1@data)
indo1@data <- plyr::join(indo1@data, gabung_label_nama, by = "NAME_1")
new_indo1 <- fortify(indo1)
## Regions defined for each Polygons
new_indo1 <- plyr::join(new_indo1, indo1@data, by="id")
new_indo1
warna <- rainbow(33)
Indo <- peta_indonesia +
scale_color_manual(values= c(warna) )+
geom_polygon(data = new_indo1, aes(x=long, y = lat, group = group, fill = NAME_1), color = "grey") +
scale_fill_manual(values=c(warna)) +
theme(legend.position = "none")
Indo
Ekstrak Label
label_kota = new_indo1 %>% distinct(NAME_1, .keep_all = TRUE)
label_kota
ambil_label = read.table("label.txt", header = T, sep = "\t")
ambil_label
## woe_label postal latitude longitude
## 1 Aceh AC 4.4153300 96.9956
## 2 North Sumatra SU 2.3630400 99.2161
## 3 West Sumatra SB -0.6426110 100.6110
## 4 Riau RI 0.3968920 101.7450
## 5 Jambi JA -1.6549700 102.8230
## 6 South Sumatra SL -3.3824000 104.0730
## 7 Bengkulu BE -3.4860600 102.3680
## 8 Lampung LA -4.9190000 105.0590
## 9 Bangka-Belitung BB -2.9581700 106.8190
## 10 Kepulauan Riau KR -0.1426390 104.6010
## 11 DKI Jakarta JK -6.2246200 106.8370
## 12 West Java JR -6.9076300 107.6380
## 13 Central Java JT -7.2901000 109.8960
## 14 Special Region of Yogyakarta YO -7.8353300 110.4430
## 15 East Java JI -7.8812900 112.6160
## 16 Banten BT -6.2579400 106.1670
## 17 Bali BA -8.3412000 115.1790
## 18 West Nusa Tenggara NB -8.5710600 116.2940
## 19 East Nusa Tenggara NT -8.5844300 120.6890
## 20 West Kalimantan KB -0.0636474 111.3040
## 21 Central Kalimantan KT -1.8421700 113.2860
## 22 South Kalimantan KS -3.0071300 115.4510
## 23 East Kalimantan KI 1.2891500 116.3540
## 24 North Kalimantan KU 3.0730930 116.0414
## 25 North Sulawesi SW 0.8530390 124.4460
## 26 Central Sulawesi ST -1.3009700 120.1530
## 27 South Sulawesi SE -3.7444000 119.9900
## 28 South East Sulawesi SG -4.0746000 122.1190
## 29 Gorontalo GO 0.7605010 122.3310
## 30 West Sulawesi SR -2.6961500 119.3160
## 31 Maluku MA -8.2481400 128.1610
## 32 North Maluku LA -2.2894500 125.9760
## 33 West Papua IB -1.3252500 132.8250
## 34 Papua PA -4.1082500 138.6890
Indo <- peta_indonesia +
scale_color_manual(values= c(warna) )+
geom_polygon(data = new_indo1, aes(x=long, y = lat, group = group, fill = NAME_1), color = "grey") +
scale_fill_manual(values=c(warna)) +
theme(legend.position = "none") +
geom_point(data = ambil_label,aes(x = longitude, y = latitude ),
color = "black", size=2.5, alpha = 0.5, show.legend = F) +
geom_text_repel(data = ambil_label, aes(longitude, latitude, label = woe_label),
size = 3)
Indo
