我正在尝试绘制一些空间数据,我得到了一个奇怪的结果,我之前没有见过问及回答过。
我像这样拉入形状数据(使用中间步骤将区域ID与ObjectIDs相加):
shp<-readOGR('./2016 precincts', '2016_Precincts')
shp<-spTransform(shp, CRS("+proj=longlat +datum=WGS84"))
shp<-tidy(shp)
当我执行以下操作时:
ggmap(map)+geom_polygon(aes(x=long, y=lat,group=group), size=.1, color=borderColor, data=testShp, alpha=0.5)
我得到了一张漂亮的空白地图。但是,一旦我合并了任何数据并尝试填写区域,就像这样:
shp<-merge(shp, precinctTotals, by.x='vtds', by.y='precinct',all.x=T)
边框线正确停止绘制。我无法弄清楚为什么。似乎唯一改变的是数据框中有额外的数据。有没有人经历过这个?