Нет проблем:

>>> t = np.array([[1,1,1,1,1],[2,2,2,2,2],[3,3,3,3,3],[4,4,4,4,4],[5,5,5,5,5]])
>>> x = np.arange(5).reshape((-1,1)); y = np.arange(5)
>>> print (t[[x]],t[[y]])

Большая проблема:

>>> s = scipy.sparse.csr_matrix(t)
>>> print (s[[x]].toarray(),s[[y]].toarray())
Traceback (most recent call last):
  File "<pyshell#22>", line 1, in <module>
:               :
:               :
ValueError: data, indices, and indptr should be rank 1

s.toarray()[[x]] прекрасно работает, но побеждает всю мою цель, используя разреженные матрицы, так как мои массивы слишком большие. Я проверил Атрибуты и Методы, связанные с некоторыми из разреженных матриц для всего, что ссылается на Расширенное индексирование, но не играло в кости. Есть идеи?