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