spikesorting.py 文件源码

def __find_temperature(self,tree,mintemp,maxtemp,tempstep,min_clus):
        num_temp = int(floor(float(maxtemp-mintemp)/tempstep))
        aux = np.diff(tree[:,4])
        aux1 = np.diff(tree[:,5])
        aux2 = np.diff(tree[:,6])
        aux3 = np.diff(tree[:,7])
        temp=0;
        for t in range(0,num_temp-1):
            if(aux[t] > min_clus or aux1[t] > min_clus or aux2[t] > min_clus or aux3[t] >min_clus):
                temp=t+1

        if (temp==0 and tree[temp][5]<min_clus):
            temp=1

        return temp