[tools] numpy: adds intersect2d

numpy.py

@@ -75,3 +75,15 @@ def template_to(inp, template):
         xsel(np.argsort(inp, axis=-1), 1, np.argsort(np.argsort(template, axis=-1), axis=-1), 1),
         axis=-1,
     )
+
+
+def intersect2d(A, B, *args, **kwargs):
+    nrows, ncols = A.shape
+    dtype = {"names": ["f{}".format(i) for i in range(ncols)], "formats": ncols * [A.dtype]}
+    C = np.intersect1d(A.view(dtype), B.view(dtype), *args, **kwargs)
+    if len(C) > 0:
+        c = C[0]
+    else:
+        c = C
+    c = c.view(A.dtype).reshape(-1, ncols)
+    return C