如何用numpy数组中的掩码值替换掩码值

Question

使用旁边最接近的未屏蔽值的平均值替换数组中屏蔽值的最有效的numpy方法是什么？

例如：

a = np.array([2,6,4,8])
b = np.ma.masked_where(a>5,a)
print b

masked_array(data = [2 -- 4 --],
         mask = [False  True False  True],
   fill_value = 999999)

我希望b中的掩码值替换为它们旁边的值的平均值。边界可以重复最接近的未屏蔽值。所以在这个例子中，b将是以下内容：

b = [2,3,4,4]

这个问题的主要原因是看看这是否可以在不使用迭代器的情况下有效地完成。

Answer 1

您可以使用np.interp和np.where

import numpy as np

a = np.array([2,6,4,8])
mymask = a>5
b = np.ma.masked_where(mymask,a)

print b
# [2 -- 4 --]

c = np.interp(np.where(mymask)[0],np.where(~mymask)[0],b[np.where(~mymask)[0]])
b[np.where(mymask)[0]] = c

print b
# [2 3 4 4]