如何在R中按经纬度重塑向量数据为网格格式并写入NetCDF

现有R语言生成的数据框df，需要将其中的m_1、m_2、m_3三个变量转换为**450（纬度，对应Y）×350（经度，对应X）**的网格矩阵格式，同时得到长度为350的唯一经度向量X、长度为450的唯一纬度向量Y，最终格式需匹配array(1:100, c(450,350))，用于写入NetCDF文件。

可复现代码（修正后） #

注意：原始代码中length(5000)为错误写法（返回值为1），已修正为匹配网格维度的总数据点数：

# 定义网格维度：纬度450，经度350
n_lat <- 450
n_lon <- 350

# 生成经纬度网格并展开为向量
X <- rep(seq(from = -10, to = -5, length.out = n_lon), each = n_lat)
Y <- rep(seq(from = 51, to = 55, length.out = n_lat), times = n_lon)

# 生成变量数据
m_1 <- seq(from = 1, to = 100, length.out = n_lat * n_lon)
m_2 <- seq(from = 1, to = 100, length.out = n_lat * n_lon)
m_3 <- seq(from = 1, to = 100, length.out = n_lat * n_lon)

# 构建数据框
df <- data.frame(m_1, m_2, m_3, X, Y)

转换方法 #

方法1：直接矩阵重塑（高效简洁）

前提是数据框中的X和Y已按网格顺序排列，若未排序需先执行df <- df[order(df$Y, df$X), ]：

# 获取唯一的经纬度向量
unique_X <- unique(df$X)
unique_Y <- unique(df$Y)

# 将变量转换为450×350矩阵
m1_matrix <- matrix(df$m_1, nrow = length(unique_Y), ncol = length(unique_X))
m2_matrix <- matrix(df$m_2, nrow = length(unique_Y), ncol = length(unique_X))
m3_matrix <- matrix(df$m_3, nrow = length(unique_Y), ncol = length(unique_X))

# 验证维度
dim(m1_matrix)  # 输出：450 350

方法2：使用reshape2包重塑数据

适合数据顺序混乱的场景，自动按经纬度分组：

library(reshape2)

# 转换m_1为网格矩阵（去掉第一列的Y值，保留经度列作为矩阵的列）
m1_matrix <- dcast(df, Y ~ X, value.var = "m_1")[-1] %>% as.matrix()
m2_matrix <- dcast(df, Y ~ X, value.var = "m_2")[-1] %>% as.matrix()
m3_matrix <- dcast(df, Y ~ X, value.var = "m_3")[-1] %>% as.matrix()

关键说明 #

• 必须确保原始数据的总点数等于450×350=157500，否则无法完整填充网格。
• 转换后的矩阵行对应纬度（Y的唯一值），列对应经度（X的唯一值），符合NetCDF的常规数据格式要求。

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