使用geom_sf与gganimate制作地图动画时遇CRS不匹配错误

核心问题分析 #

静态分面图正常但动画报错，本质是合并后的sf数据在帧迭代时出现CRS不一致，NA警告则是年度数据存在区域缺失导致的。以下是针对性解决步骤：

1. 强制统一CRS（必做） #

即使静态图显示正常，gganimate在逐帧渲染时会严格检查CRS一致性，需确保所有年度的空间数据CRS完全匹配：

# 先确认原始UK地图的CRS（通常为EPSG:27700，英国国家坐标系）
st_crs(uk_map)

# 确保合并后的动画数据是标准sf对象，并重设CRS
uk_animated_data <- st_as_sf(uk_animated_data)
uk_animated_data <- st_set_crs(uk_animated_data, st_crs(uk_map))

# 验证各年度子集的CRS是否一致（输出应为单一EPSG码）
unique(sapply(split(uk_animated_data, uk_animated_data$year), function(x) st_crs(x)$epsg))

如果验证输出有多个值，说明部分年度的空间数据CRS丢失，需重新合并或修复。

2. 处理NA值消除警告 #

transition_states的NA警告通常是部分年度缺失区域数据导致的，可通过补全或过滤解决：

补全缺失数据（推荐）

用tidyr::complete补全所有年度的所有区域，缺失计数用0填充：

library(tidyr)
# 先转为tibble补全，再转回sf对象
uk_animated_data <- uk_animated_data %>%
  as_tibble() %>%
  complete(year, objectid, fill = list(count = 0)) %>%
  st_as_sf() %>%
  st_set_crs(st_crs(uk_map)) # 转回后需重新设置CRS

过滤缺失数据

如果不需要保留缺失区域，直接过滤：

uk_animated_data <- uk_animated_data %>%
  filter(!is.na(count), !is.na(objectid))

3. 修正gganimate代码写法 #

避免在动画中重复设置CRS或出现类型不兼容问题：

用transition_time的正确示例

library(gganimate)
p <- ggplot(uk_animated_data) +
  geom_sf(aes(fill = count)) +
  coord_sf(crs = st_crs(uk_map)) + # 全局统一设置CRS
  scale_fill_viridis_c(option = "plasma") +
  labs(title = "年度计数: {frame_time}") +
  transition_time(year) +
  ease_aes("linear")

# 渲染动画
animate(p, renderer = gifski_renderer(), width = 800, height = 600)

用transition_states的正确示例（需将year转为因子）

# 先将year转为因子，避免NA警告
uk_animated_data$year <- factor(uk_animated_data$year)

p <- ggplot(uk_animated_data) +
  geom_sf(aes(fill = count)) +
  coord_sf(crs = st_crs(uk_map)) +
  labs(title = "年度计数: {closest_state}") +
  transition_states(year, transition_length = 1, state_length = 1)

animate(p)

4. 排查合并过程的潜在问题 #

确保objectid在两个数据框中类型一致，避免合并时出现隐式不匹配：

# 统一objectid为整数类型
uk_map$objectid <- as.integer(uk_map$objectid)
yearCountsAllRegions$objectid <- as.integer(yearCountsAllRegions$objectid)

# 重新执行合并
uk_animated_data <- left_join(uk_map, yearCountsAllRegions, by = "objectid")

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