来自shapefile的奇异多边形

Question

我正在尝试将shapefile应用于ggmaps地图,但它给了我非常奇怪的结果.有问题的shape文件是"统计局域网"(组类似邮编)shape文件从可用统计澳大利亚统计局这里.

通常情况下,我可能认为这是一个切断边缘点的问题,但即使在缩放级别1我也会遇到它(事实上它看起来更糟):

以下是我用于生成上述图表的一些代码:

library(tidyverse)
library(ggmap)
library(rgdal)

slas <- readOGR(dsn="SLA",layer="SLA11aAust")

aus4 <- get_map("Australia",zoom=4)
ggmap(aus4)

ggmap(aus4)+
  geom_polygon(data=slas, aes(x=long,y=lat))

aus1 <- get_map("Australia",zoom=1)
ggmap(aus1)

ggmap(aus1)+
  geom_polygon(data=slas, aes(x=long,y=lat))

我做错了什么,或者shapefile是否以某种方式错误配置？

Answer 1

我认为你只需要(可选)fortify变量slas,不要忘记group并用颜色使边界可见:

slas <- fortify(slas, region = "SLA_CODE11")
ggmap(aus4) + 
  geom_polygon(data = slas2, color = "white", aes(x = long, y = lat, group = group))