以矩阵为元素的数组的块状相交

Question

我有两个数组，一个是形状(200000, 28, 28)，另一个是shape (10000, 28, 28)，所以实际上是两个以矩阵为元素的数组。现在，我想计算并获得(N, 28, 28)在两个数组中都重叠的所有元素（以形式）。使用普通的for循环，这是减慢速度的方法，因此我尝试使用numpys
intersect1d方法进行尝试，但我不知道如何将其应用于此类数组。

Answer 1

使用关于唯一行的此问题中的方法)

def intersect_along_first_axis(a, b):
    # check that casting to void will create equal size elements
    assert a.shape[1:] == b.shape[1:]
    assert a.dtype == b.dtype

    # compute dtypes
    void_dt = np.dtype((np.void, a.dtype.itemsize * np.prod(a.shape[1:])))
    orig_dt = np.dtype((a.dtype, a.shape[1:]))

    # convert to 1d void arrays
    a = np.ascontiguousarray(a)
    b = np.ascontiguousarray(b)
    a_void = a.reshape(a.shape[0], -1).view(void_dt)
    b_void = b.reshape(b.shape[0], -1).view(void_dt)

    # intersect, then convert back
    return np.intersect1d(b_void, a_void).view(orig_dt)

请注意，使用void浮点数是不安全的，因为它将导致-0不等于0