4

我有两个代表二维坐标的 numpy 数组。每行代表(x, y)对:

a = np.array([[1, 1], [2, 1], [3, 1], [3, 2], [3, 3], [5, 5]])
b = np.array([[1, 1], [5, 5], [3, 2]])

我想删除有效a的元素b。所以结果是:

array([[2, 1], [3, 1], [3, 3]])

我可以通过循环和比较来做到这一点,我希望我能做得更容易。

4

1 回答 1

1

Pythonsets在给出差异方面做得很好。但是,它并不能维持秩序

np.array(list(set(tuple(x) for x in a.tolist()).difference(set(tuple(x) for x in b.tolist()))))

或者要使用布尔索引,使用广播创建一个outer equals, 和sumwith anyandall

A = np.all((a[None,:,:]==b[:,None,:]),axis=-1)
A = np.any(A,axis=0)
a[~A,:]

或使ab复杂:

ac = np.dot(a,[1,1j])
bc = np.dot(b,[1,1j])
A = np.any(ac==bc[:,None],axis=0)
a[~A,:]

或使用setxor1d

xx = np.setxor1d(ac,bc)
# array([ 2.+1.j,  3.+1.j,  3.+3.j])
np.array([xx.real,xx.imag],dtype=int).T

==================

In [222]: ac = np.dot(a,[1,1j])
     ...: bc = np.dot(b,[1,1j])
In [223]: ac
Out[223]: array([ 1.+1.j,  2.+1.j,  3.+1.j,  3.+2.j,  3.+3.j,  5.+5.j])
In [225]: bc
Out[225]: array([ 1.+1.j,  5.+5.j,  3.+2.j])
In [226]: ac == bc[:,None]
Out[226]: 
array([[ True, False, False, False, False, False],
       [False, False, False, False, False,  True],
       [False, False, False,  True, False, False]], dtype=bool)
于 2013-10-01T01:28:28.803 回答