R语言二维数组条件切片：提取指定经纬度区域位置

核心需求 #

从二维经纬度矩阵中筛选满足60 < lat < 70且-60 < long < -50的位置，再用这些位置提取三维温度数组（第三维度为深度）的对应数据。

你的代码存在的问题 #

1. filter函数误用：filter是tidyverse中处理数据框的工具，不能直接用于数组/矩阵，会导致逻辑错误。
2. between参数顺序错误：between的正确语法是between(目标变量, 下限, 上限)，你写的between(60,70)把数值顺序搞反了，完全不符合筛选逻辑。
3. 不必要的数组合并：其实不需要把经纬度合并成三维数组，直接用二维矩阵的逻辑筛选更高效。

正确解法 #

1. 处理真实NetCDF数据场景

假设你的温度数组名为temp_data（维度为1442x398x深度），步骤如下：

# 读取经纬度矩阵
nav_lat <- ncvar_get(nc,"nav_lat")
nav_long <- ncvar_get(nc,"nav_lon")

# 创建筛选掩码：标记所有满足条件的经纬度位置
region_mask <- nav_lat > 60 & nav_lat < 70 & nav_long > -60 & nav_long < -50

# 提取符合条件的经纬度坐标
selected_lat <- nav_lat[region_mask]
selected_long <- nav_long[region_mask]

# 提取对应区域的温度数据（保留所有深度维度）
selected_temp <- temp_data[region_mask, ]

2. 修正你的小型可复现示例

先修正数组构造逻辑（确保每个(i,j)位置对应一对经纬度），再执行筛选：

library(tidyverse)

# 构造经纬度矩阵（每个位置对应一对坐标）
lat_mat <- matrix(c(5,9,3, 4,11,2, 6,8,1), nrow=3, ncol=3)
long_mat <- matrix(c(10,11,12, 9,13,10, 11,12,9), nrow=3, ncol=3)
# 合并为三维数组：[行,列,1]是纬度，[行,列,2]是经度
result <- array(c(lat_mat, long_mat), dim=c(3,3,2))

# 创建筛选掩码：纬度在5-10之间，经度在10-12之间
mask <- between(result[,,1], 5, 10) & between(result[,,2], 10, 12)

# 提取符合条件的坐标对
selected_coords <- cbind(
  latitude = result[,,1][mask],
  longitude = result[,,2][mask]
)

# 查看结果
selected_coords

运行后会得到符合条件的坐标对：

latitude longitude
[1,]        5        10
[2,]        9        11
[3,]        6        11
[4,]        8        12

关键逻辑说明 #

• 掩码（mask）：是一个和经纬度矩阵维度相同的逻辑矩阵，TRUE代表该位置满足筛选条件，FALSE则不满足。
• 数组索引：用逻辑矩阵作为索引，可以直接提取所有满足条件的元素，温度数组的前两个维度和经纬度矩阵一一对应，因此用同一个掩码就能提取对应位置的所有深度数据。

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