R语言求助：创建以经纬度为中心的网格化地图（分辨率0.5）

1. 准备工作：安装并加载所需包 #

先确保安装了必要的R包，执行以下命令：

install.packages(c("raster", "readr"))
library(raster)
library(readr)

2. 读取数据 #

如果是本地CSV文件，用read_csv读取；如果直接用你提供的示例数据，可直接构造数据框：

# 读取本地CSV（替换为你的文件路径）
# df <- read_csv("your_data.csv")

# 示例数据构造
df <- data.frame(
  Lon = c(78,78,79,79,83,92,93,94,94),
  Lat = c(14.5,15.5,14.5,15.5,18.5,27.5,27.5,27.5,28.5),
  P = c(14.64,15.06,20.25,23.08,27.04,0.12,15.97,15.14,1.29)
)

3. 创建匹配分辨率的栅格模板 #

因为你需要0.5分辨率且以给定经纬度为中心的栅格，每个栅格的经纬度范围是中心±0.25。先计算整体数据的边界范围，再创建空栅格：

# 计算栅格边界（中心±0.25）
x_min <- min(df$Lon) - 0.25
x_max <- max(df$Lon) + 0.25
y_min <- min(df$Lat) - 0.25
y_max <- max(df$Lat) + 0.25

# 创建空栅格，设置分辨率为0.5
r_template <- raster(
  xmn = x_min, xmx = x_max,
  ymn = y_min, ymx = y_max,
  res = 0.5
)

4. 将P值填充到栅格中 #

用rasterize函数把数据点的P值匹配到对应栅格单元格（你的数据中每个点对应唯一栅格中心，直接赋值即可）：

# 填充P值到栅格
r_p <- rasterize(df[, c("Lon", "Lat")], r_template, field = df$P)

5. 绘制网格化地图 #

用plot函数绘制栅格，可添加颜色和标注优化效果：

# 绘制栅格地图
plot(r_p, 
     main = "P值网格化地图",
     col = terrain.colors(10), # 使用地形色阶
     xlab = "经度", ylab = "纬度")

# 可选：叠加原始数据点验证
points(df$Lon, df$Lat, pch = 16, cex = 0.8, col = "black")

常见报错原因排查 #

• 栅格范围错误：如果直接用原始经纬度的最小/最大值作为栅格边界，会导致部分栅格中心不匹配，必须加上±0.25的偏移量。
• 列顺序错误：rasterize要求坐标列必须是经度（x）在前，纬度（y）在后，如果顺序反了会导致匹配失败。
• 包版本问题：确保raster包是最新版本，旧版本可能存在函数参数兼容问题。

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