ggplot2按SpatialPolygonsDataFrame的task_color列为网格地图单元格着色

核心错误原因 #

你遇到的报错和无法正确填充的问题来自三个错误用法：

1. SpatialPolygonsDataFrame属于sp包的空间对象，ggplot的几何图层无法直接读取对象中@data槽里的属性列，需要先转换为带属性关联的表格数据
2. 你把fill = task_color写在了aes()映射外，ggplot会把它识别为全局固定填充值，不会按每行数据匹配对应的属性值
3. 如果你要直接使用task_color列中存储的颜色字符串作为填充色，需要额外加scale_fill_identity()声明，否则ggplot会把颜色字符串识别为分类标签，映射默认调色板的颜色

实现代码 #

方法1：基于传统sp对象转换（兼容旧代码） #

先把空间对象的坐标和属性合并为绘图用的完整表格：

# 转换空间坐标为长表
grid_df <- fortify(grid)
# 提取属性表并添加匹配用的id列
grid_attr <- grid@data
grid_attr$id <- rownames(grid_attr)
# 合并坐标与属性
grid_full <- merge(grid_df, grid_attr, by = "id")

绘图代码：

ggplot() +
  # 多边形填充，把fill放在aes内映射task_color列
  geom_polygon(data = grid_full, aes(x = long, y = lat, group = group, fill = task_color), alpha = 0.3) +
  # 多边形描边
  geom_path(data = grid_full, aes(x = long, y = lat, group = group), size = 2.5, alpha = 0.3, col = 'pink') +
  # 声明直接使用列中的颜色值作为填充
  scale_fill_identity() +
  # 空间图固定坐标比例
  coord_equal()

方法2：基于sf包简化实现（推荐） #

现在ggplot2原生支持sf格式的简单要素对象，不需要手动合并坐标和属性：

library(sf)
# 直接将sp对象转换为sf对象，自动保留所有属性
grid_sf <- st_as_sf(grid)

# 绘图代码，一行即可同时实现填充和描边
ggplot() +
  geom_sf(data = grid_sf, aes(fill = task_color), alpha = 0.3, col = 'pink', linewidth = 2.5) +
  scale_fill_identity()

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