>>> a = np.arange(16).reshape(4, 4)
>>> a
array([[ 0, 1, 2, 3],
[ 4, 5, 6, 7],
[ 8, 9, 10, 11],
[12, 13, 14, 15]])
>>> iu1 = np.triu_indices(4)
>>> iu1
(array([0, 0, 0, 0, 1, 1, 1, 2, 2, 3]), array([0, 1, 2, 3, 1, 2, 3, 2, 3, 3])) #表示 哪行那列的 矩阵数据
>>> a[iu1]
array([ 0, 1, 2, 3, 5, 6, 7, 10, 11, 15]) #按照索引取值
>>> iu1 = np.triu_indices(4,1)
>>> iu1
(array([0, 0, 0, 1, 1, 2]), array([1, 2, 3, 2, 3, 3]))
>>> a[iu1]
array([ 1, 2, 3, 6, 7, 11])
eucl_dist = euclidean_distances(norm, norm)[np.triu_indices(T, k=1)] # 100 *100的距离矩阵
数量应该 是 (100*100 - 主对角元素100)/2 = 4950 的一维向量