结构化数组的掩码是否应该自行结构化？

我正在调查numpy 问题 2972 https://github.com/numpy/numpy/issues/2972以及几个相关的问题。事实证明，所有这些问题都与数组本身是结构化的情况有关，但它的掩码不是：

In [38]: R = numpy.zeros(10, dtype=[("A", "<f2"), ("B", "<f4")])

In [39]: Rm = numpy.ma.masked_where(R["A"]<5, R)

In [41]: Rm.dtype
Out[41]: dtype([('A', '<f2'), ('B', '<f4')])

In [42]: Rm.mask.dtype
Out[42]: dtype('bool')

# Now, both `__getitem__` and `__repr__` will result in errors — see issue #2972

如果我以不同的方式创建掩码数组，则掩码数据类型的结构类似于数组本身的数据类型：

In [44]: Q.dtype
Out[44]: dtype([('A', '<f4'), ('B', '<f4')])

In [45]: Q.mask.dtype
Out[45]: dtype([('A', '?'), ('B', '?')])

前一种情况暴露出几个问题。例如，Rm.__repr__() and Rm["A"]两者都导致IndexError，虽然这是一个ValueError在过去。

根据设计，该模式应该是可能的，其中A.dtype是结构化的，但是A.mask.dtype是不是结构化？

换句话说：是错误__repr__ and __getitem__中的方法numpy.ma.core.MaskedArray，或者真正的错误是否发生在之前——通过允许这样一个掩码结构化数组存在？

第一种情况中的错误表明这些方法期望掩码具有与基本数组相同的字段数量（和名称）

__getitem__:  dout._mask = _mask[indx]
_recursive_printoption: (curdata, curmask) = (result[name], mask[name])

如果掩码数组是使用“main”构造函数创建的，则掩码具有相同的结构

Rn = np.ma.masked_array(R, mask=R['A']>5)
Rn.mask.dtype: dtype([('A', '?'), ('B', '?')])

换句话说，每个元素的每个字段都有一个掩码值。

The masked_array医生显然打算将“相同形状”包括在内dtype结构。Mask: Must be convertible to an array of booleans with the same shape as 'data'.

如果我尝试以相同的方式设置蒙版masked_where does

Rn._mask=R['A']>5

我得到同样的打印错误。结构化掩码被新的布尔值覆盖，从而更改其数据类型。相反，如果我使用

Rn.mask=R['A']<5

Rn打印效果很好。.mask是一个属性，其set方法显然可以正确处理结构化掩模。

如果不深入研究代码历史（在 github 上），我的猜测是masked_where是一个便利函数，当结构数据类型添加到结构的其他部分时，该函数不会更新ma代码。相比ma.masked_array这是一个简单的函数，根本不考虑数据类型。其他便利功能，例如ma.masked_greater use masked_where。改变result._mask = cond to result.mask = cond可能是纠正此问题所需的全部内容。

您对非结构化掩模的后果进行了多彻底的测试？

Rm.flatten()

返回一个带有结构化掩码的数组，即使它以非结构化掩码开头。那是因为它使用Rm.__setmask__，对字段敏感。这就是set函数为mask财产。

Rm.tolist()  # same error as str()

masked_where以。。开始：

cond = make_mask(condition)

make_mask返回简单的“bool”数据类型。也可以使用 dtype 调用它，生成结构化掩码：np.ma.make_mask(R['A']<5,dtype=R.dtype)。但这种结构化掩模在使用时会变平masked_where. masked_where不仅允许非结构化掩码，而且还强制其非结构化。

您的非结构化掩码已经部分实现，recordmask财产：

recordmask = property(fget=_get_recordmask)

我说部分是因为它有一个get方法，但是set方法尚未实现。看def _set_recordmask(self):

我越看这个就越确信masked_where是错的。它可以更改为设置结构化掩码，但它与masked_array。如果在构造数组时引发错误（有dtype.names）。那样masked_where对于非结构化数值数组仍然有用，同时防止误用到结构化数组。

我还应该看看测试代码。