来自 netCDF 文件的 numpy 数组中的无效值（打印为“--”）

Question

我通过netCDF4 将netCDF 文件读入python。 我得到具有无效值的数组，这些数组打印为--，例如：

[[[-- -- -- ..., 0.0 0.0 0.0]
  [-- -- -- ..., 0.0 0.0 0.0]
  [-- -- -- ..., 0.0 0.0 0.0]]

 [[-- -- -- ..., 0.0 0.0 0.0]
  [-- -- -- ..., 0.0 0.0 0.0]
  [-- -- -- ..., 0.0 0.0 0.0]]]

这些值从何而来，我如何识别（和替换）它们？诸如numpy.isnan 之类的东西也会返回给我--。我可以稍后将这些数组与-- 一起使用而不会崩溃（例如，我可以绘制它们），但它们似乎对某些操作有问题。

例如，我曾经做过类似的事情

numpy.mean(myarray, axis=(1,2)) # the tuple for `axis` is not very regular

在myarray 上源自我的netCDF 文件之前没有任何问题，但我得到TypeError: tuple indices must be integers, not tuple 并且这些奇怪的-- 填充数组崩溃...

我正在使用python 2.7.9（如果与打印到屏幕格式相关，PyCharm 用于代码编辑）。

Answer 1

这里有一个masked array，其中-- 代表被屏蔽的元素。

只需使用myarray.data 和myarray.mask 分别获取数据和无效元素。

很可能，一个掩码的 numpy 数组被故意保存在 netcdf4 文件中，这可能是有原因的。所以我会简单地继续使用这个掩码数组。正如您所注意到的，大多数 Numpy/Scipy/Matplotlib 函数（例如 numpy.mean）都有特定的逻辑来处理它们。