如何将叠加SpatialPointsDataFrame到栅格后的地图保存为变量

这是个很常见的基础绘图小坑——R的基础plot()函数本身不会返回可复用的绘图对象，它只是把图形直接绘制到当前的绘图设备里，所以你把它赋值给变量自然会得到NULL。下面给你几个实用的解决方案：

方案1：把绘图逻辑封装成自定义函数（基础绘图友好） #

把你需要的绘图步骤打包成一个函数，之后每次需要生成地图时直接调用函数就行，不用重复写两行plot代码：

# 定义可复用的绘图函数
plot_bird_map <- function(map_raster, bird_points) {
  plot(map_raster)
  plot(bird_points, add = TRUE, col = "blue")
}

# 第一次绘制地图
plot_bird_map(map, yellowbirds)

# 后续需要重绘时，直接调用函数
plot_bird_map(map, yellowbirds)

方案2：改用ggplot2（推荐，支持保存绘图对象） #

如果愿意切换到ggplot2生态，它的绘图逻辑是基于可保存的对象设计的，你可以把整个绘图定义保存为变量，后续随时打印重绘，还能灵活修改样式：

library(ggplot2)
library(rasterVis) # 用于将raster转换为ggplot兼容格式

# 转换数据格式
map_df <- as.data.frame(map, xy = TRUE)  # raster转data.frame
yellowbirds_df <- as.data.frame(yellowbirds)  # SpatialPointsDataFrame转data.frame

# 创建并保存ggplot绘图对象
overlay_map <- ggplot() +
  geom_raster(data = map_df, aes(x = x, y = y, fill = layer)) +  # 绘制底图
  geom_point(data = yellowbirds_df, aes(x = coords.x1, y = coords.x2), color = "blue") +  # 叠加鸟点
  theme_minimal()

# 第一次绘图
print(overlay_map)

# 后续重绘只需打印变量
print(overlay_map)

方案3：用recordPlot()记录基础绘图状态（应急用） #

如果非要用基础绘图且想保存当前绘图状态，可以用recordPlot()记录，之后用replayPlot()重绘，但这个方法局限性较大（比如设备尺寸变化可能导致重绘异常）：

# 先完成绘图
plot(map)
plot(yellowbirds, add = TRUE, col = "blue")

# 记录当前绘图状态到变量
overlay_map <- recordPlot()

# 重绘时调用replayPlot
replayPlot(overlay_map)

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