ggplot2与ggmap奇形多边形R

Question

我正在尝试使用 ggplot2 绘制多边形 shapefile，但得到了一些奇怪的结果。我使用以下代码读取具有单个多边形的 shapefile：

zctaSp<-readShapePoly("zctaSp.shp")

然后我使用几种不同的方法进行绘图，其中两种可行，但我需要的那一种不行。

情节 1：形状看起来正确：

plot(zctaSp)

情节 2：形状看起来也正确，这似乎与情节 2 几乎相同：

ggplot(data=zctaSp, aes(x=long, y=lat, group=group)) + geom_polygon()

情节 3：但形状被这样弄乱了：

ggplot(data=zctaSp, aes(x=long, y=lat, group=group)) + geom_polygon()    
atl <- qmap('atlanta', zoom=11, color="bw")
atl + geom_polygon(data=zctaSp, aes(x=long, y=lat, group=group), alpha=1)

我将 shapefile 放在http://bit.ly/1nnlAg3。

请注意，在按照 Hadley Wickham 在 this link 的说明运行 fortify 命令后，我也尝试过绘图，但这并没有改善问题。

Answer 1

尝试：

library(ggmap)
library(rgdal)
# Data using NAD83 - epsg: 4269
zct <- readOGR(dsn = 'D:/Programacao/R/Stackoverflow/22387624',
               layer = 'zctaSp')
zctdf <- fortify(zct)

# Project to wgs84
wgs84proj <- CRS('+init=epsg:4326')
zct_g <- spTransform(zct, wgs84proj)
zctgdf <- fortify(zct_g)
map_loc <- get_map(location = c(lon = mean(zctgdf$lon), mean(zctgdf$lat)),
                   source = 'google', zoom = 11)
map <- ggmap(map_loc, extent = 'device')
map + 
  geom_polygon(data=zctgdf, aes(x=long, y=lat, group=group), alpha=.8)

来自 Qgis 的地图

关于本练习中使用的预测，NAD 83 / WGS84 可能是一个陷阱。 NAD83 和 WGS84 之间存在细微差别。 NAD83 依赖于 GRS80 基准，其实现与 WGS84 非常相似，但又不一样。

对于这些 ggmap 应该始终使用未投影的 WGS84 (epsg4326)。