R语言ggplot2地图自定义颜色及NA值问题解决问询

一、手动绑定fit.mean刻度与指定颜色的实现方法 #

完全可以实现，以ggplot2为例（假设你使用该绘图包），操作步骤如下：

1. 按0、0.05…1的顺序定义对应21个刻度的颜色向量，示例如下：

# 替换为你实际需要的颜色值
custom_colors <- c(
  "black", "#2D2D2D", "#4A4A4A", "#686868", "#868686",
  "#A3A3A3", "#C1C1C1", "#DCDCDC", "#EAEAEA", "#F2F2F2",
  "#F5F5F5", "#F8F8F8", "#FAFAFA", "#FCFCFC", "#FEFEFE",
  "#ECECEC", "#E0E0E0", "#D3D3D3", "#C6C6C6", "#B9B9B9",
  "#ACACAC"
)
# 对应的刻度值序列
custom_breaks <- seq(0, 1, by = 0.05)

2. 修改绘图代码，将原有连续色阶标度替换为scale_fill_gradientn，绑定自定义颜色与刻度：

library(ggplot2)
library(sf)

ggplot() +
  geom_sf(data = SHP_0, aes(fill = fit.mean), color = "white") +
  # 替换原scale_fill_continuous等标度
  scale_fill_gradientn(
    colors = custom_colors,
    breaks = custom_breaks,
    limits = c(0, 1),
    guide = guide_colorbar(ticks = TRUE, ticks.colour = "black")
  ) +
  theme_minimal()

若需要分段着色而非连续过渡，可改用scale_fill_stepsn，参数用法与上述一致。

二、NA值排查与着色解决方法 #

排查原因 #

1. 数据缺失验证：先确认意大利对应的fit.mean是否真为NA，执行代码查看：

# 假设国家名称列名为"country_name"，根据实际列名调整
italy_data <- subset(d.fit.spatial.1805, country_name == "Italy")
print(italy_data$fit.mean)

若输出为NA，说明1805年意大利的拟合数据确实缺失；若输出不为NA，大概率是空间数据与拟合数据的连接键不匹配（如国家名称拼写不一致），需检查数据合并逻辑。

2. 负值范围确认：德国的负值正常着色，说明色标范围已覆盖负值区间，意大利的NA与负值无关，属于数据本身缺失或匹配错误。

NA值着色解决方法 #

在填充色标度参数中添加na.value，指定NA区域的颜色，示例：

ggplot() +
  geom_sf(data = d.fit.spatial.1805, aes(fill = fit.mean), color = "white") +
  scale_fill_continuous(
    na.value = "lightgray",  # 替换为你需要的NA颜色
    limits = c(min(d.fit.spatial.1805$fit.mean, na.rm = TRUE), max(d.fit.spatial.1805$fit.mean, na.rm = TRUE))
  ) +
  theme_minimal()

无论使用scale_fill_gradientn还是其他色标函数，均可通过na.value参数设置NA区域的颜色。

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