Portanto, em matrizes numpy existe a função interna para obter os índices diagonais, mas não consigo descobrir como obter a diagonal começando do canto superior direito em vez do superior esquerdo.

Este é o código normal para começar a partir do canto superior esquerdo:

>>> import numpy as np
>>> array = np.arange(25).reshape(5,5)
>>> diagonal = np.diag_indices(5)
>>> array
array([[ 0,  1,  2,  3,  4],
   [ 5,  6,  7,  8,  9],
   [10, 11, 12, 13, 14],
   [15, 16, 17, 18, 19],
   [20, 21, 22, 23, 24]])
>>> array[diagonal]
array([ 0,  6, 12, 18, 24])

então o que eu uso se eu quiser retornar:

array([ 4,  8, 12, 16, 20])