def importAccousticData(index, resampling):
    tdms_file = TdmsFile(dataPath+'_%i'%index+'.tdms')
    accoustic_channel = tdms_file.object('Untitled', 'Accelerometer')
    time=accoustic_channel.time_track()[::resampling]
    signal=accoustic_channel.data[::resampling]
    return time,signal

def importShearData(index, resampling):
    tdms_file = TdmsFile(dataPath+'_%i'%index+'.tdms')
    try:
        shear_load_channel = tdms_file.object('Untitled', 'Shear load')
    except:
        shear_load_channel = tdms_file.object('Untitled', 'shear load')
    time=shear_load_channel.time_track()[::resampling]
    signal=shear_load_channel.data[::resampling]
    return time,signal

def aggregateStrainSignal(start,end, resampling):
    time,signal=importAccousticData(start, resampling)
    for i in range (1+start,end):
        time_temp,signal_temp=importAccousticData(i, resampling)
        time_temp=time_temp+time[-1]+time[1]-time[0]
        time=np.append(time,time_temp,axis=0)
        signal=np.append(signal,signal_temp,axis=0)
    #Convert acceleration data to strain data: strain=velocity/speed of sound
    #with velocity=acceleration/angular frequency and angular frequency=2*pi*f
    #and f=40.3kHz
    signal=(31.6/0.0449)*signal/(700*2*np.pi*40.3E3)
    time=time+3.6*start
    return time,signal


def aggregateShearSignal(start,end, resampling):
    time,signal=importShearData(start, resampling)
    for i in range (1+start,end):
        time_temp,signal_temp=importShearData(i, resampling)
        time_temp=time_temp+time[-1]+time[1]-time[0]
        time=np.append(time,time_temp,axis=0)
        signal=np.append(signal,signal_temp,axis=0)
    time=time+3.6*start
    return time,signal