r - 来自 shapefile 的奇异多边形

Question

我正在尝试将 shapefile 应用于 ggmaps map ，但它给了我非常奇怪的结果。有问题的 shapefile 是澳大利亚统计局提供的“统计局部区域”(类似于邮政编码的组)shapefile here .

通常我可能会认为这是边缘点被截断的问题，但即使在缩放级别 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

我认为你只需要(可选)强化变量slas，不要忘记group并使边界可见一种颜色:

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