提取 3D 阵列数据 - 编程之家

如何解决提取 3D 阵列数据

我在 R 方面有一些经验，但这是一个我从未深入研究过的新领域。

简短说明：我有一堆具有 2D 矩阵格式（纬度/经度）的连续月度文件，我需要将它们组合成一个 3D 数组（日期构成第 3 维，我想我可以做这部分）。然后，我需要提取一组特定的坐标，以从我拥有文件的所有月份的纬度/经度中获取连续的 NO2 值列表（基本上将数组侧向并提取一列日期和浓度）。有没有一种简单的（或者至少不是疯狂的技术）方法来做这样的事情？

详细解释：

我有来自 OMI 站点的每月 NO2 浓度矩阵。格式是经度和纬度的 2D 矩阵（分别为 -180 到 180 和 -90 到 90，增量为 0.125 度）。所以 2880 列和 1440 行。具体格式为：

ncols 2880
nrows 1440
xllcorner -180.0
yllcorner -90.0
细胞大小 0.125
nodata_value -999
2.0版本
-999 -999 -999 -999 -999 -999 -999 -999 -999 -999 -999 -999 <..4> -999 -999

在第一个 -999 之前有一个前导“空格”，但我不知道如何格式化它。标题（以 2.0 版结尾）在所有文件上（自始至终完全相同）。我可以在 ArcMap 中单独打开这些文件，并轻松找到特定月份的 NO2 值。但是我有这些文件 15 年的价值，并且将来可能能够获得每日价值（GB 价值的数据......）。我宁愿不必加载每个单独的文件并选择浓度。我很想获取所有文件，将它们连接成一个 3D 数组（纬度 x 长 x 日期），然后选择一个特定坐标从该坐标中吐出日期和 NO2 浓度的列表。听起来这不是一项艰巨的任务，但我已经很长时间没有编程了，并且在我的搜索中找不到解决方案。我也不确定提供更多上下文的其他信息。如果我在正确的论坛中，请告诉我。

谢谢！