使用ggmap绘制形状文件：当形状文件大于ggmap时进行裁剪

当我尝试将 ggmap 与形状文件组合时，我遇到了剪切问题。 Kahle 和 Wickham (2013: 158) 中的示例工作正常，因为来自 ggmap 的光栅图像覆盖了整个形状文件。下面是当我尝试在覆盖较小区域的 ggmap 图上绘制美国各州的形状文件时发生的情况的示例。 ggmap 显示纽约市，我想将其与美国各州的边界重叠（仅作为示例）。由此产生的地图没有任何意义。问题是形状文件被剪裁并且 ggplot 连接未剪裁的点。下面是代码。形状文件来自here http://www.arcgis.com/home/item.html?id=f7f805eb65eb4ab787a0a3e1116ca7e5。我只是在这里展示最后一个情节。

我怎么解决这个问题？

path <- "PATH TO SHAPEFILE"
library("ggmap")
library("rgdal")

# shapefile
states <- readOGR(dsn = path, layer = "states")
states_df <- fortify(states)
# plot shapefile
plot(states, lwd = 0.1)
ggplot(states_df, aes(long, lat, group = group)) +
    geom_polygon(colour = "black", fill = NA, size = 0.1)


# combine ggmap with shapefile
map <- get_map("new york city", zoom = 10, source = "stamen")
ggmap(map, extent = "device")

ggmap(map, extent = "device") +
    geom_polygon(aes(long, lat, group=group), data = states_df, colour = "red", fill = NA, size = 1)

卡勒、大卫和哈德利·威克姆。 2013.“Ggmap：使用 ggplot2 进行空间可视化。” R 期刊 5(1)：144–61。

这是我的尝试。我经常使用 GADM shapefile，您可以使用raster包裹。我对 NY、NJ 和 CT 的形状文件进行了子集化。您最终可能不必这样做，但减少数据量可能会更好。当我绘制地图时，ggplot 自动删除了 ggmap 图像 bbox 之外的数据点。因此，我不需要做任何额外的工作。我不确定你使用的是哪个形状文件。但是，GADM 的数据似乎与 ggmap 图像配合得很好。希望这对您有帮助。

library(raster)
library(rgdal)
library(rgeos)
library(ggplot2)

### Get data (shapefile)
us <- getData("GADM", country = "US", level = 1)

### Select NY and NJ
states <- subset(us, NAME_1 %in% c("New York", "New Jersey", "Connecticut"))

### SPDF to DF
map <- fortify(states)

## Get a map
mymap <- get_map("new york city", zoom = 10, source = "stamen")


ggmap(mymap) +
geom_map(data = map, map = map, aes(x = long, y = lat, map_id = id, group = group))

如果您只想要线条，那么以下就是您所追求的。

ggmap(mymap) +
geom_path(data = map, aes(x = long, y = lat, group = group))