在R中用ggplot/ggmap绘制Ames数据集社区凸包遇坐标匹配问题求助

错误核心原因 #

这个报错本质是ggmap图层默认传递lon/lat美学映射，而geom_sf不需要这些字段；同时如果混用不匹配的坐标系投影（比如ggmap默认的快速投影和EPSG:3857），会进一步导致图层无法对齐。

修正后的完整代码 #

# 加载所需工具包
library(tidyverse)
library(sf)
library(ggmap)
library(modeldata) # 提供Ames数据集

# 1. 构建并转换Ames的SF对象到EPSG:3857
data(ames)
ames_sf <- ames %>%
  # 原始经纬度为WGS84（EPSG:4326），先转为SF对象
  st_as_sf(coords = c("Longitude", "Latitude"), crs = 4326) %>%
  # 转换到谷歌地图使用的Web墨卡托投影（EPSG:3857）
  st_transform(crs = 3857)

# 2. 生成社区凸包多边形
neighborhood_hulls <- ames_sf %>%
  group_by(Neighborhood) %>%
  # 合并同社区的点位，再生成凸包
  summarize(geometry = st_union(geometry)) %>%
  st_convex_hull()

# 3. 获取匹配坐标系的谷歌地图背景
# 从SF对象的边界框计算地图中心
bbox <- st_bbox(ames_sf)
map_center <- c(mean(c(bbox[["xmin"]], bbox[["xmax"]])), 
                mean(c(bbox[["ymin"]], bbox[["ymax"]])))

# 获取谷歌地图，指定CRS为3857
ames_map <- get_googlemap(
  center = map_center,
  zoom = 13, # 根据实际范围调整缩放级别
  maptype = "roadmap",
  crs = 3857
)

# 4. 叠加图层并可视化
ggmap(ames_map) +
  # 关键：添加inherit.aes=FALSE，避免继承ggmap的lon/lat美学映射
  geom_sf(data = neighborhood_hulls,
          aes(fill = Neighborhood),
          alpha = 0.3, # 透明化便于观察社区重叠
          color = "black",
          inherit.aes = FALSE) +
  # 强制使用EPSG:3857坐标系，确保图层对齐
  coord_sf(crs = 3857) +
  theme_minimal() +
  labs(title = "Ames社区凸包与谷歌地图叠加")

关键修正细节 #

• inherit.aes = FALSE：添加到geom_sf中，彻底规避ggmap的lon/lat字段要求对SF图层渲染的干扰。
• 统一坐标系：确保SF对象和谷歌地图都使用EPSG:3857，同时用coord_sf()替代ggmap默认的coord_quickmap()——后者仅适配WGS84经纬度投影，无法匹配Web墨卡托。
• 透明化凸包：设置alpha = 0.3，完美实现《Tidy-modelling with R》中易观察社区重叠的可视化效果。

内容的提问来源于stack exchange，提问作者TheGoat