R 形状文件中的奇异多边形

R 形状文件中的奇异多边形,r,R,我试图将一个shapefile应用到ggmaps地图,但它给了我非常奇怪的结果。讨论中的shapefile是澳大利亚统计局提供的“统计本地区”(类似于邮政编码的组)shapefile 通常我可能会认为这是一个切割边缘点的问题,但即使在缩放级别1,我也会遇到这个问题(事实上,它看起来更糟): 下面是我用来制作上述图表的一些代码: library(tidyverse) library(ggmap) library(rgdal) slas <- readOGR(dsn="SLA",laye

我试图将一个shapefile应用到ggmaps地图,但它给了我非常奇怪的结果。讨论中的shapefile是澳大利亚统计局提供的“统计本地区”(类似于邮政编码的组)shapefile

通常我可能会认为这是一个切割边缘点的问题,但即使在缩放级别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))
库(tidyverse)
图书馆(ggmap)
图书馆(rgdal)
sla我认为您只需要(可选地)
加强变量
sla
,不要忘了
分组
,并使用颜色显示边界:

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

sla你说得对,看起来好多了。
fortify
做什么?文档现在只是说“改用
broom
包”,这。。。没有帮助。
fortify
sla
空间多边形数据帧
转换为常规的
数据帧
。严格地说,这是没有必要的;缺少的关键是
group=group
参数。