Question

我想仅绘制this data set中ethnic变量的外边界。您在此链接下找到的数据集plot_data是来自苏丹PRIO GRID的shapefile的强化数据集。

我目前的代码如下：

plot_data <- load("plot_data.rdata")

ggplot(plot_data, 
       aes(x= long, 
           y = lat, 
           group = id)) + 
  geom_polygon() + 
  geom_polygon(data = plot_data %>% 
                 filter(!is.na(ethnic)) %>% # subset data with ethnicity only
                 as.data.frame(), 
               aes(color = ethnic)) + 
  coord_equal()

这给了我以下输出：

但是，我想删除指定种族群体的居住区域的所有内部线路，并仅绘制该区域的外边界。

我没有得到this similar problem的解决方案。

我可能需要在之前组合单元格我强化shapefile;但也许在强化之后还有一种方法。我试图删除重复的坐标，但这不起作用。任何建议将不胜感激。感谢。

Answer 1

基本思路很简单：你需要将填充和颜色映射到整个绘图的美学，然后将颜色映射到较小多边形的美学，这样你就可以用相同的颜色为较小的多边形着色。更大的多边形。这就是摆脱较小多边形内线的原因。然后我在较小多边形的美学之外添加了尺寸，以便将红线的宽度增加到可感知的水平。

plot_data <- load("plot_data.rdata")   

library(ggplot2) 

df<-plot_data %>% 
  filter(!is.na(ethnic)) %>% # subset data with ethnicity only
  as.data.frame()    

ggplot(plot_data, 
       aes(x= long, y = lat, group = id,fill="", color="")) + 
  geom_polygon() + 
  geom_polygon(data = df, 
               aes(color = ethnic), size=1) + 
  geom_polygon(data=df, aes(x= long, y = lat, group = id))+
  scale_fill_manual(values="black", guide=F)+
  scale_color_manual(name="ethnic",
               labels=c("","Shaygiyya, Ja'aliyyin and Danagla (Arab)"),
               values=c("black","red"), 
               breaks = c("NA","Shaygiyya, Ja'aliyyin and Danagla (Arab)")) +
  coord_equal()

我的会话信息：

sessionInfo()
R version 3.3.1 (2016-06-21)
Platform: x86_64-apple-darwin13.4.0 (64-bit)
Running under: OS X 10.11.5 (El Capitan)

locale:
[1] en_US.UTF-8/en_US.UTF-8/en_US.UTF-8/C/en_US.UTF-8/en_US.UTF-8

attached base packages:
[1] stats     graphics  grDevices utils     datasets  methods   base     

loaded via a namespace (and not attached):
 [1] Rcpp_0.12.7        XML_3.98-1.4       proj4_1.0-8        bitops_1.0-6      
 [5] MASS_7.3-45        grid_3.3.1         plyr_1.8.3         gtable_0.2.0      
 [9] scales_0.4.0       KernSmooth_2.23-15 ggplot2_2.1.0      ash_1.0-15        
[13] RColorBrewer_1.1-2 RJSONIO_1.3-0      tools_3.3.1        RSelenium_1.4.2   
[17] munsell_0.4.3      RCurl_1.95-4.8     maps_3.1.0         colorspace_1.2-6  
[21] caTools_1.17.1

更新： OP使用我发布的代码制作相同的情节令人不安。我有一个朋友使用我在这里提供的相同代码在Windows系统（我在El Capitan上）重现该情节，他得到了相同的结果（下图）。

仅在ggplot2 map / geom_polygon

1 个答案: