예제 #1
0
def ahurst(settings, signal, data):
    object = wla.toolbox()
    transtype = settings.transtype
    l = settings.l
    type = settings.type
    refl = settings.refl
    j_0 = settings.j_0
    j_min = settings.j_min
    j_max = settings.j_max#Must set j_0 wvar needs this
    bias = settings.bias
    p = settings.p
    n = len(signal.e1) #Assume same length of all epochs!!!!
    print l, type, refl, j_0, j_min, j_max, bias, p, n
    
    for i in range(0,j_0):
        data.y1.append(0)
        data.y2.append(0)
        data.y3.append(0)
    for j in range (0, j_0):
        data.y_hat1.append(0)
        data.y_hat2.append(0)
        data.y_hat3.append(0)
    for k in range(0,2*j_0):
        data.ci1.append(0)
        data.ci2.append(0)
        data.ci3.append(0)
    #e4 = list_to_array(data.e4)
    e1 = list_to_array(signal.e1)
    e2 = list_to_array(signal.e2)
    e3 = list_to_array(signal.e3)
    y1 = list_to_array(data.y1)
    y2 = list_to_array(data.y2)
    y3 = list_to_array(data.y3)
    #y4 = list_to_array(data.y4)
    y_hat1 = list_to_array(data.y_hat1)
    y_hat2 = list_to_array(data.y_hat2)
    y_hat3 = list_to_array(data.y_hat3)
    #wvar4 = list_to_array(data.y_hat4)
    ci1 = list_to_array(data.ci1)
    ci2 = list_to_array(data.ci2)
    ci3 = list_to_array(data.ci3)
    #ci4 = list_to_array(data.ci4)
    out1 = list_to_array(data.out1)
    out2 = list_to_array(data.out2)
    out3 = list_to_array(data.out3)
    #out4 = list_to_array(data.out4)

    #if settings.filetype is 1:
    #object.ahurst(transtype, l, type, bias, refl, n, j_0, j_min, j_max, p, e1, y1, y_hat1, ci1, out1)#What about j_0???
    #object.ahurst(transtype, l, type, bias, refl, n, j_0, j_min, j_max, p, e2, y2, y_hat2, ci2, out2)
    #object.ahurst(transtype, l, type, bias, refl, n, j_0, j_min, j_max, p, e3, y3, y_hat3, ci3, out3)

    #if settings.filetype is 0:
    #    object.ahurst(transtype, l, type, bias, refl, n, j_min, j_max, p, e4, y4, wvar4, out4)#Error with wvar

    
    for i in range(0,3):
        data.out1[i] = carray_get(out1,i)
        data.out2[i] = carray_get(out2,i)
        data.out3[i] = carray_get(out3,i)
예제 #2
0
def wa(settings, signal, data):
    object = wla.toolbox()
    transtype = settings.transtype
    l = settings.l
    type = settings.type
    refl = settings.refl
    j_0 = settings.j_0
    n = len(signal.etime) #Assume constant length!!
    print transtype, l, type, refl, j_0, n

    for i in range(0, 2*(j_0+1)):
        data.boundaries.append(0)
    e1 = list_to_array(signal.signal)
    e2 = list_to_array(signal.e2)
    e3 = list_to_array(signal.e3)
    boundaries = list_to_array(data.boundaries)
    wv1 = new_cmatrix(j_0+1,n)
    wv2 = new_cmatrix(j_0+1,n)
    wv3 = new_cmatrix(j_0+1,n)

    object.wa(transtype, l, type, refl, j_0, n, boundaries, e1, wv1)
    object.wa(transtype, l, type, refl, j_0, n, boundaries, e2, wv2)
    object.wa(transtype, l, type, refl, j_0, n, boundaries, e3, wv3)
    #Check for filetype and do extra round...
    if transtype is 0:
        for i in range(j_0+1,0,-1):
            level = i-1
            dummy = 0
            tmp_list1 = []
            for j in range(0,n):
                tmp_list1.append(0)
            for k in range(0,n,2**i):
                tmp_list1[k] = (cmatrix_get(wv1,level,dummy))
                dummy +=1
            data.wv1.append(tmp_list1)
        #boundaries too...
    elif transtype is 1:
        for i in range(j_0+1,0,-1): #Now it plots only the first half of V_j0 
            j = i-1
            tmp_list1 = []
            tmp_list2 = []
            tmp_list3 = []
            for k in range(0,n):
                tmp_list1.append(cmatrix_get(wv1,j,k))
                tmp_list2.append(cmatrix_get(wv2,j,k))
                tmp_list3.append(cmatrix_get(wv3,j,k))
            data.wv1.append(tmp_list1)
            data.wv2.append(tmp_list2)
            data.wv3.append(tmp_list3)

    data.boundaries[0] = carray_get(boundaries, 2*j_0)
    data.boundaries[j_0+1] = carray_get(boundaries, 2*j_0+1)
    j = 1
    for i in range(j_0,0,-1):
        data.boundaries[j] = carray_get(boundaries, i-1)
        j += 1
    j = j_0+2
    for i in range(2*j_0,j_0,-1):
        data.boundaries[j] = carray_get(boundaries, i-1)
        j += 1
예제 #3
0
def shurst(data):
    l = 8
    type = 2
    bias = 0
    refl = 0
    j_0 = 11
    j_min = 5
    j_max = 11
    n = len(data.e2)
    n_b = 2
    signal = list_to_array(data.e2)
    object = wla.toolbox()
    object.shurst(l, type, refl, j_0, j_min, j_max, n, n_b, signal)
예제 #4
0
def mra(settings, signal, data):
    object = wla.toolbox()
    transtype = settings.transtype
    l = settings.l
    type = settings.type
    refl = settings.refl
    j_0 = settings.j_0
    n = len(signal.e1)
    for i in range(0, 2*(j_0+1)):
        data.boundaries.append(0)
    e1 = list_to_array(signal.e1)
    e2 = list_to_array(signal.e2)
    e3 = list_to_array(signal.e3)
    boundaries = list_to_array(data.boundaries)
    ds1 = new_cmatrix(j_0+1,n)
    ds2 = new_cmatrix(j_0+1,n)
    ds3 = new_cmatrix(j_0+1,n)

    object.mra(transtype, l, type, refl, j_0, n, boundaries, e1, ds1)
    object.mra(transtype, l, type, refl, j_0, n, boundaries, e2, ds2)
    object.mra(transtype, l, type, refl, j_0, n, boundaries, e3, ds3)

    for i in range(j_0+1,0,-1):
        j = i-1
        tmp_list1 = []
        tmp_list2 = []
        tmp_list3 = []
        for k in range(0,n):
            tmp_list1.append(cmatrix_get(ds1,j,k))
            tmp_list2.append(cmatrix_get(ds2,j,k))
            tmp_list3.append(cmatrix_get(ds3,j,k))
        data.ds1.append(tmp_list1)
        data.ds2.append(tmp_list2)
        data.ds3.append(tmp_list3)

    data.boundaries[0] = carray_get(boundaries, 2*j_0)
    data.boundaries[j_0+1] = carray_get(boundaries, 2*j_0+1)
    j = 1
    for i in range(j_0,0,-1):
        data.boundaries[j] = carray_get(boundaries, i-1)
        j += 1
    j = j_0+2
    for i in range(2*j_0,j_0,-1):
        data.boundaries[j] = carray_get(boundaries, i-1)
        j += 1
예제 #5
0
def ihurst(settings, signal, data):
    object = wla.toolbox()
    l = settings.l
    type = settings.type
    refl = settings.refl 
    j_0 = settings.j_0
    j_min = settings.j_min
    j_max = settings.j_max
    n = len(signal.e1)
    print l, type, refl, j_0, j_min, j_max, n
    for i in range(0, n):
        data.hurst1.append(0)
        data.hurst2.append(0)
        data.hurst3.append(0)

    e1 = list_to_array(signal.e1)
    e2 = list_to_array(signal.e2)
    e3 = list_to_array(signal.e3)
    hurst1 = list_to_array(data.hurst1)
    hurst2 = list_to_array(data.hurst2)
    hurst3 = list_to_array(data.hurst3)
    out1 = list_to_array(data.out1)
    out2 = list_to_array(data.out2)
    out3 = list_to_array(data.out3)
    
    object.ihurst(l, type, refl, j_0, j_min, j_max, n, e1, hurst1, out1)
    object.ihurst(l, type, refl, j_0, j_min, j_max, n, e2, hurst2, out2)
    object.ihurst(l, type, refl, j_0, j_min, j_max, n, e3, hurst3, out3)

    for i in range(0,3):
        data.out1[i] = carray_get(out1,i)
        data.out2[i] = carray_get(out2,i)
        data.out3[i] = carray_get(out3,i)
    for i in range(0,n):
        data.hurst1[i] = carray_get(hurst1,i)
        data.hurst2[i] = carray_get(hurst2,i)
        data.hurst3[i] = carray_get(hurst3,i)
예제 #6
0
def analysis():
    object = wla.toolbox()
    signal_list = []#Set signal as "global" here
    y_list = [] #Get these as input
    yhat_list = []#Convert to arrays here
    ci_list = []
    output_list = [0,0,0]
    hurst_list = []
    boundaries_list = []
##
#
#    for value in open('/Users/vegard/Master/Hurst05.dat'): #('/Users/vegard/Master/ocean2.dat'):
#        signal_list.append(float(value.split()[1])) #For reading Ingve's files...
#        hurst_list.append(0)
#       	#signal_list.append(float(value))
#    for i in range(j_min-1,j_max):
#        y_list.append(0)
#    for j in range(0, 3*j_max):
#        wvar_list.append(0)
#    for k in range(0, 2*(j_max+1)):
#        boundaries_list.append(0)
#    n = len(signal_list)
#    signal = list_to_array(signal_list)
#    y = list_to_array(y_list)
#    wvar = list_to_array(wvar_list)
#    output = list_to_array(output_list)
#    hurst = list_to_array(hurst_list)
#    boundaries = list_to_array(boundaries_list)
#    wv = new_cmatrix(j_max+1,n)
#    ds = new_cmatrix(j_max+1,n)
#    object = wla.toolbox()
#    #
#
#    if(wa.get_active() == 1):
#        transtype = wa.get_transtype()
#        l = wa.get_l()
#        type = wa.get_type()
#        refl = wa.get_refl()
#        j_max = wa.get_jmax()
#
#        #Get epochs
#        object.wa(transtype, l, type, refl, j_max, n, boundaries, signal, w)
#    
#    if(mra.get_active() == 1):
#        transtype = mra.get_transtype()
#        l = mra.get_l()
#        type = mra.get_type()
#        refl = mra.get_refl()
#        j_max = mra.get_jmax()
#        
#        #Get epochs
#        object.mra(transtype, l, type, refl, j_max, n, boundaries, signal, ds)
#    
#    if(ahurst.get_active() == 1):
#        transtype = ahurst.get_transtype()
#        l = ahurst.get_l()
#        type = ahurst.get_type()
#        refl = ahurst.get_refl()
#        j_min = ahurst.get_jmin()
#        j_max = ahurst.get_jmax()
#        bias = ahurst.get_bias()
#        p = ahurst.get_p()
#    
#        #Get epochs
#        object.ahurst(transtype, l, type, bias, refl, n, j_min, j_max, p, signal, y, wvar, output) #
#
#    if(ihurst.get_active() == 1):
#        l = ihurst.get_l()
#        type = ihurst.get_type()
#        refl = ihurst.get_refl() #Not implemented yet in C++
#        j_min = ihurst.get_jmin()
#        j_max = ihurst.get_jmax()
#    
#        #Get epochs
#        object.ihurst(l, type, j_min, j_max, len(signal_list), signal, hurst, output)
#        
        
    transtype = 1
    l = 8
    type = 2
    bias = 0
    refl = 0
    j_0 = 6
    j_min = 2
    j_max = 6
    p = 95
    
    for value in open('/Users/vegard/Master/ecg.dat'): #('/Users/vegard/Master/ocean2.dat'):
        #signal_list.append(float(value.split()[1])) #For reading Ingve's files...
        hurst_list.append(0)
      	signal_list.append(float(value))
    for i in range(0,j_0):
        y_list.append(0)
    for j in range(0, j_0):
        yhat_list.append(0)
    for k in range(0,2*j_0):
        ci_list.append(0)
    for k in range(0, 2*(j_max+1)):
        boundaries_list.append(0)
    n = len(signal_list)
    signal = list_to_array(signal_list)
    y = list_to_array(y_list)
    yhat = list_to_array(yhat_list)
    ci = list_to_array(ci_list)
    output = list_to_array(output_list)
    hurst = list_to_array(hurst_list)
    boundaries = list_to_array(boundaries_list)
    wv = new_cmatrix(j_max+1,n)
    #ds = new_cmatrix(j_max+1,n)
    #object = wla.toolbox()
   #For iHurst use j_min:1, j_max:10
    object.ahurst(transtype, l, type, bias, refl, n, j_0, j_min, j_max, p, signal, y, yhat, ci, output) 
예제 #7
0
    
    signal_list = []
    #for value in open('/Users/vegard/Master/ecg.dat'): #('/Users/vegard/Master/ocean2.dat'):
    #    signal_list.append(float(value.split()[1]))
        #signal_list.append(float(value))
    l = 8
    type = 2
    bias = 0
    refl = 0
    j_0 = 6
    j_min = 1
    j_max = 6
    n = len(signal_list)
    n_b = 6
    signal = list_to_array(signal_list)
    object = wla.toolbox()
#    object.shurst(l, type, refl, j_0, j_min, j_max, n, n_b, signal)

    analysis()
    #ahurstdata = Data()
    #ihurstdata = Data()
    #genset = Settings()
    #ahurstset = Settings()
    #ihurstset = Settings()
    #test = [1,2,3]
    #tittel = "Testfil"
    #ofilename = "/Users/vegard/Master/saves/testfile"
    #ahurstdata.out1 = test
    #ahurstdata.out2 = test
    #ahurstdata.out3 = test
    #save(ofilename, tittel, genset, ahurstset, ihurstset, ahurstdata, ihurstdata)