Рисование границ провинции в R

Question

Я хочу нарисовать карту choropleth для моего отчета. Я загрузил данные из GADM. Мне удалось создать карту choropleth, но когда я добавил границы, произошло что-то неправильное. Появляется много синей линии.

Мой код:

#Read Database 
HCM <- read.csv("HCM.csv") 

#Load Packages 
library(plyr) 
library(sp) 
library(ggplot2) 
library(rgeos) 
library(maptools) 
library(scales) 

#Read Map (Map downloaded from GADM) 
vie_map1 <- readRDS('VNM_adm1.rds') 
vie_map <- fortify(vie_map1, region = "ID_1") 

#Draw the map 
ggplot() + 
    geom_map(data = HCM, aes(map_id = ID, fill = WeightRange), 
      map = vie_map) + expand_limits(x = vie_map$long, y = vie_map$lat) 

Но когда я хочу добавить границы, что-то неладное произошло (прямая линия в синий цвет). Код для добавления границ:

# Add the boundaries 
ggplot() + 
    geom_map(data = HCM, aes(map_id = ID, fill = WeightRange), 
      map = vie_map) + expand_limits(x = vie_map$long, y = vie_map$lat)+ 
geom_path() + 
    geom_path(data = vie_map, aes(x = long, y = lat), color = "blue") + 
    theme_bw() 

Скажите, пожалуйста, как это исправить (удалите прямую линию).

Answer 1

Также установлено group в пределах эстетики отображений geom_path не только x и y:

library(ggplot2) 
library(raster) 
library(maptools) 
vie_map1 <- getData("GADM", country="VNM", level=1) 
# trying URL 'http://biogeo.ucdavis.edu/data/gadm2.8/rds/VNM_adm1.rds' 
# Content type '' length 2765340 bytes (2.6 MB) 
# downloaded 2.6 MB 
vie_map <- fortify(vie_map1, region = "ID_1") 
ggplot() + 
    geom_polygon(data = vie_map, aes(x = long, y = lat, group = group), fill = "grey90") + 
    geom_path(data = vie_map, aes(x = long, y = lat, group = group), color = "blue") + 
    theme_bw()