使用ggplot2绘制RasterLayer对象报错，求解决方案及示例数据

问题原因 #

geom_raster() 需要输入包含坐标（x、y）和栅格值的数据框，直接传入RasterLayer（S4对象）会触发类型不匹配错误；普通的as.data.frame()转换会丢失坐标信息，导致转换后的数据无法用于绘图。

正确转换与绘图方法 #

单个RasterLayer的处理 #

使用raster::as.data.frame()并指定xy=TRUE参数，保留栅格的坐标信息：

library(raster)
library(ggplot2)

# 加载或定义你的RasterLayer
r <- raster(system.file("external/test.grd", package="raster"))

# 转换为带坐标的数据框
r_df <- as.data.frame(r, xy = TRUE)

# 用geom_raster绘图
ggplot(r_df, aes(x = x, y = y, fill = layer)) +
  geom_raster() +
  scale_fill_viridis_c() # 可选：优化颜色映射

列表中多个RasterLayer的处理 #

如果栅格存储在列表中，可批量转换并绘制（结合purrr和dplyr简化操作）：

library(purrr)
library(dplyr)

# 创建示例栅格列表示例
list_raster <- list(
  raster(system.file("external/test.grd", package="raster")),
  raster(nrow=100, ncol=100, vals=runif(10000))
)

# 批量转换为带坐标的数据框，并添加标识列
list_df <- imap(list_raster, ~mutate(as.data.frame(.x, xy=TRUE), raster_name = paste("栅格", .y)))

# 合并为单个数据框，绘制多子图
combined_df <- bind_rows(list_df)
ggplot(combined_df, aes(x=x, y=y, fill=layer)) +
  geom_raster() +
  facet_wrap(~raster_name) +
  scale_fill_viridis_c()

可复现的栅格示例数据 #

无需外部下载，直接使用R包内置数据或生成随机栅格：

# 方法1：raster包内置测试数据
r_test <- raster(system.file("external/test.grd", package="raster"))

# 方法2：生成随机栅格
r_random <- raster(nrow=50, ncol=50)
values(r_random) <- rnorm(ncell(r_random)) # 正态分布随机值

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