使用BFAST及bfastSpatial绘制年份颜色图例断点值问题咨询

• 第一步：预处理断点 raster 结果
  先提取断点结果中的数值部分，所有值为浮点型（例如2015.12对应2015年2月），首先过滤掉无断点的无效值（通常为NA或超出研究时段的数值），将有效值统一限制在你的研究时间区间内。
• 第二步：自定义年份色阶映射
  不需要直接调用自带的changeMonth函数，手动生成和年份匹配的连续颜色映射即可，R环境下示例代码如下：

  library(RColorBrewer)
  # 替换为你自己的研究起止年份
  start_year <- 2000
  end_year <- 2020
  # 生成对应年份数量的渐变色，可替换为你需要的色系
  year_color_pal <- colorRampPalette(brewer.pal(11, "Spectral"))(end_year - start_year + 1)
  

• 第三步：出图匹配目标效果
  调用raster::plot绘制断点图层时，设置breaks参数为每年的起始值，col参数为上一步生成的色阶，自定义图例时将刻度标签设置为对应年份即可，最终效果和参考图一致：
  

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