如何使用ggplotly()将ggplot地图转为plotly并添加时间交互滑块？

你现有代码存在两个直接导致崩溃/功能失效的问题：

• 变量名前后不一致：inner_join后生成的对象是data_cleaned，ggplot调用的是不存在的data1
• 映射写法错误：aes(fill = "wage")是将填充色映射到固定字符串"wage"，没有关联实际工资字段，也没有绑定时间维度字段，转成plotly后没有帧属性，后续手动加滑块时会重复加载sf几何数据，内存溢出直接导致R崩溃。

可行实现方案 #

不要在单帧ggplotly对象上硬加滑块，优先用原生帧映射实现，能避免几何数据重复存储的问题：

1. 先完成数据预处理
   确保你的最低工资表包含时间维度字段（示例用year字段，可替换为你实际的时间列名），同一州在不同时间点对应不同的最低工资数值：

   library(sf)
   library(readr)
   library(plotly)
   library(dplyr)
   library(ggplot2)
   
   map <- read_sf("map.shp")
   minimumwage <- read_csv("minimumwage.csv")
   data_cleaned <- inner_join(map, minimumwage, by = "States")
   # 将时间字段转为有序因子，保证滑块顺序正确
   data_cleaned$year <- as.ordered(data_cleaned$year)
   

2. 推荐方案：直接用plotly原生接口绘制带动画的地图，稳定性最高，内存占用最低
   这种方式所有时间帧共用一套sf几何数据，不会重复存储几何信息，不会出现崩溃问题：

   # 固定全时段色阶范围，避免切换年份时颜色标尺跳变
   wage_limits <- range(data_cleaned$wage, na.rm = TRUE)
   
   p <- plot_ly() %>%
     add_sf(
       data = data_cleaned,
       # 核心：将frame映射到时间字段，plotly会自动生成时间滑块
       frame = ~year,
       fill = ~wage,
       fillcolor = ~colorRampPalette(c("#f7fbff", "#08306b"))(100)[
         scales::rescale(wage, to = c(1, 100), from = wage_limits)
       ],
       stroke = I("#333"),
       linewidth = I(0.5),
       hoverinfo = "text",
       text = ~paste0("地区：", States, "<br>最低工资：", wage, "<br>年份：", year)
     ) %>%
     # 动画参数配置，关闭不必要的重绘大幅降低内存占用
     animation_opts(
       frame = 500,
       transition = 200,
       redraw = FALSE
     ) %>%
     # 滑块样式配置
     animation_slider(
       currentvalue = list(prefix = "当前年份：", font = list(size = 14))
     ) %>%
     # 播放/暂停按钮配置
     animation_button(
       x = 1, xanchor = "right", y = 0.2, yanchor = "bottom",
       label = "播放/暂停"
     ) %>%
     layout(showlegend = FALSE, title = "各州历年最低工资变化")
   
   p
   

3. 如果你习惯用ggplot语法，必须在ggplot构建阶段就指定帧映射，不要转成ggplotly后再追加

   g <- ggplot(data = data_cleaned, aes(frame = year)) +
     # 注意wage不要加引号，要映射到数据中的实际列
     geom_sf(aes(fill = wage)) +
     # 固定色阶范围
     scale_fill_continuous(limits = wage_limits)
   
   ggplotly(g) %>%
     animation_opts(redraw = FALSE) %>%
     animation_slider(currentvalue = list(prefix = "年份："))
   

避坑说明 #

• 禁止在生成ggplotly对象后手动逐帧绑定sf几何数据，否则几何数据会按时间节点数N复制N份，数据量膨胀几十倍直接触发内存崩溃
• 所有动画帧共用同一套几何边界，仅切换填充属性，能把内存占用降到原来的1/N
• 时间字段不要用无序字符串，转为有序因子或数值类型，避免滑块顺序错乱
• 如果边界数据顶点过多，可先用rmapshaper::ms_simplify()简化几何复杂度，进一步降低内存占用

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