基于R语言raster包按纬度区间为栅格地图添加分类标签的方法

嘿，我来给你分享几个可行的方案，帮你给WorldClim的最低气温栅格图按纬度区间添加分类标签：

这个方案直接用你已经在使用的raster包工具，不需要额外加载太多依赖，适合快速实现需求：

1. 先准备好数据和纬度区间规则

library(raster)
# 获取你需要的最低气温数据
dat <- getData('worldclim', var='tmin', res=10)
tmin1 <- dat$tmin1

# 定义纬度区间（你可以根据需求自由调整边界）
lat_intervals <- c(-90, -10, 0, 10, 20, 90)
# 给每个区间设置对应的标签文本
interval_labels <- c("纬度 < -10°", "-10° ~ 0°", "0° ~ 10°", "10° ~ 20°", "纬度 > 20°")

2. 计算标签的放置位置，然后绘图加标签

# 计算每个纬度区间的中心纬度，用来放置标签
centers_lat <- (lat_intervals[-1] + lat_intervals[-length(lat_intervals)])/2
# 选一个合适的经度位置放标签（这里选栅格左侧偏内的位置，避免贴边）
lon_pos <- xmin(tmin1) + 5

# 绘制栅格图
plot(tmin1, main="1月最低气温（按纬度区间分类）")
# 添加分类标签，用白色粗体让标签更显眼
text(x=lon_pos, y=centers_lat, labels=interval_labels, col="white", font=2)

进阶优化：如果觉得标签不够清晰，可以先给每个纬度区间加半透明背景框

plot(tmin1, main="1月最低气温（按纬度区间分类）")
# 给每个区间绘制半透明背景
for(i in 1:length(interval_labels)){
  rect(xleft=xmin(tmin1), ybottom=lat_intervals[i], 
       xright=xmax(tmin1), ytop=lat_intervals[i+1], 
       col=adjustcolor("black", alpha.f=0.2), border=NA)
}
# 再添加标签
text(x=lon_pos, y=centers_lat, labels=interval_labels, col="white", font=2, cex=0.9)

如果你想要更美观、可定制性更强的图表，可以结合ggplot2和sf包来实现，语法更直观，调整样式也更方便：

library(raster)
library(ggplot2)
library(sf)

# 获取数据
dat <- getData('worldclim', var='tmin', res=10)
tmin1 <- dat$tmin1

# 将栅格转换为SF对象，方便ggplot处理
tmin_sf <- st_as_sf(rasterToPoints(tmin1, spatial=TRUE))
# 添加纬度分类列，用cut函数划分区间
tmin_sf$lat_category <- cut(tmin_sf$y, 
                            breaks=c(-90, -10, 0, 10, 20, 90), 
                            labels=c("纬度 < -10°", "-10° ~ 0°", "0° ~ 10°", "10° ~ 20°", "纬度 > 20°"))

# 绘制可视化图表
ggplot() +
  # 绘制气温栅格
  geom_sf(data=tmin_sf, aes(fill=tmin1), color=NA) +
  # 设置颜色标尺，用viridis配色更友好
  scale_fill_viridis_c(name="1月最低气温(℃)") +
  # 添加纬度区间标签，位置可以自由调整
  annotate("text", x=xmin(tmin1)+5, 
           y=c(-50, -5, 5, 15, 50), 
           label=c("纬度 < -10°", "-10° ~ 0°", "0° ~ 10°", "10° ~ 20°", "纬度 > 20°"),
           color="white", fontface="bold", size=3) +
  # 可选：添加区间分隔虚线，让分类更清晰
  geom_hline(yintercept=c(-10, 0, 10, 20), color="white", linetype="dashed") +
  theme_minimal() +
  labs(title="1月最低气温（按纬度区间分类）")

小提示 #

• 你可以完全自定义lat_intervals的边界值，比如增加更多细分区间或者调整现有边界
• 标签的x/y坐标可以根据你的图表布局微调，确保标签不会遮挡关键数据区域

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