pymvpa——如何计算一个普通列表的相似性矩阵

原创，转载请指明出处！

## 此函数用来将list 转成 可以画图的 ndarray
#def listtomatrix(list):
#    lenlist=len(list)
#    m=np.array(list)  #将list 转成 ndarray
#    RDM = scipy.spatial.distance.pdist(list, 'euclidean')  
#    这里计算出来的是上三角矩阵的元素，但是全在一行，下面的步骤将数据正真转成上三角
#    rdm_uptriangle=np.zeros(shape=(lenlist,lenlist))        
##   第一行 的最后一列结尾元素是lenth-1 下标      对角线上全是0
#    h=0  #h 用来加在k上面 因为k一直都是从0开始 ，
#    j=0 # j 用来从RDM 中取元素 每取一个元素加1
#    for i in range(lenlist): # 行遍历 这个上三角矩阵 只遍历到最后一个元素即可
#        h+=1
#        for k in range(lenlist-i-1):  
#            k+=h # 每次都跳过对角和下三角区域元素      k+=h 
#            
##            print(k)
#            rdm_uptriangle[i][k]=RDM[j]
#            j+=1
#        
#    rdm=rdm_uptriangle+rdm_uptriangle.T
#    type(rdm)
#    return rdm

转成ndarray 之后

import pylab as pl
b=[[1],[2],[3],[4]]
a=listtomatrix(b)
pl.imshow(a)

结果如图，将一个list  转成 图像就出来了

常用函数查找，看源代码  Ctrl +Alt 大法好