Верхняя диагональ двумерного массива numpy

Это кажется простым (и тривиально написать трехстрочный цикл для), но как я могу использовать numpy-срез, составьте список местоположений индекса верхней диагонали массив numpy? Т.е.

Учитывая массив 4x4, я бы хотел, чтобы индексные позиции были у X:

[ X X X X ]
[ 0 X X X ]
[ 0 0 X X ]
[ 0 0 0 X ]

Давать:

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