删除列表中包含的 numpy 行？

Question

我有一个 numpy 数组和一个列表。我想删除列表中包含的行。

a = np.zeros((3, 2))
a[0, :] = [1, 2]
l = [(1, 2), (3, 4)]

目前我尝试通过创建一组a 的行来做到这一点，然后排除从l 创建的set，类似于：

sa = set(map(tuple, a))
sl = set(l)
np.array(list(sa - sl))

或者更简单

sl = set(l)
np.array([row for row in list(map(tuple, a)) if row not in sl]

当每一行都很短时，这些效果很好。

有没有更快的方法？我需要优化速度。

Answer 1

方法#1：这是views（将每一行视为一个元素，每个元素都具有扩展的dtype）-

# https://stackoverflow.com/a/45313353/ @Divakar
def view1D(a, b): # a, b are arrays
    a = np.ascontiguousarray(a)
    b = np.ascontiguousarray(b)
    void_dt = np.dtype((np.void, a.dtype.itemsize * a.shape[1]))
    return a.view(void_dt).ravel(),  b.view(void_dt).ravel()

a1D,l1D = view1D(a,l)
out = a[np.in1d(a1D,l1D,invert=True)]

如果您只需要像set 那样在输出中拥有唯一行，请在获得的输出上使用np.unique -

np.unique(out,axis=0)

样本运行输出 -

In [72]: a
Out[72]: 
array([[1, 2],
       [0, 0],
       [0, 0]])

In [73]: l
Out[73]: [(1, 2), (3, 4)]

In [74]: out
Out[74]: 
array([[0, 0],
       [0, 0]])
In [75]: np.unique(out,axis=0)
Out[75]: array([[0, 0]])

方法 #2： 具有相同的降维理念，这里是针对 int dtype 数据的矩阵乘法 -

l = np.asarray(l)
shp = np.maximum(a.max(0)+1,l.max(0)+1)
s = np.r_[shp[::-1].cumprod()[::-1][1:],1]
l1D = l.dot(s)
a1D = a.dot(s)
l1Ds = np.sort(l1D)
out = a[l1D[np.searchsorted(l1Ds,a1D)] != a1D]