jgrewe/gp_neurobio
1
0
Fork 0
Code Issues Pull Requests Releases Wiki Activity

nee

Browse Source
This commit is contained in:
efish 2018-11-16 17:54:53 +01:00
parent 1c875aa039
commit a77c5a0e90
2 changed files with 46 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

21
code/base_eod.py Normal file
View File

@ -0,0 +1,21 @@
from read_baseline_data import *
#import nix_helpers as nh
import matplotlib.pyplot as plt
import numpy as np
from IPython import embed #Funktionen importieren
data_dir = "../data"
dataset = "2018-11-09-aa-invivo-1"
#data = ("2018-11-09-aa-invivo-1", "2018-11-09-ab-invivo-1", "2018-11-09-ac-invivo-1", "2018-11-09-ad-invivo-1", "2018-11-13-aa-invivo-1", "2018-11-13-ab-invivo-1", "2018-11-13-ad-invivo-1", "2018-11-09-af-invivo-1", "2018-11-09-ag-invivo-1", "2018-11-13-ah-invivo-1", "2018-11-13-ai-invivo-1", "2018-11-13-aj-invivo-1", "2018-11-13-ak-invivo-1", "2018-11-13-al-invivo-1", "2018-11-14-aa-invivo-1", "2018-11-14-ab-invivo-1", "2018-11-14-ac-invivo-1", "2018-11-14-ad-invivo-1", "2018-11-14-ae-invivo-1", "2018-11-14-af-invivo-1", "2018-11-14-ag-invivo-1", "2018-11-14-aa-invivo-1", "2018-11-14-aj-invivo-1", "2018-11-14-ak-invivo-1", "2018-11-14-al-invivo-1", "2018-11-14-am-invivo-1", "2018-11-14-an-invivo-1")
time,eod = read_baseline_eod(os.path.join(data_dir, dataset))
zeit = np.asarray(time)
plt.plot(zeit[0:1000], eod[0:1000])
plt.title('A.lepto EOD')#Plottitelk
plt.xlabel('Time [ms]', fontsize = 12)#Achsentitel
plt.ylabel('Amplitude[mv]', fontsize = 12)#Achsentitel
plt.xticks(fontsize = 12)
plt.yticks(fontsize = 12)
plt.show()

25
code/vector_phase.py Normal file
View File

@ -0,0 +1,25 @@
from read_baseline_data import *
from utility import *
#import nix_helpers as nh
import matplotlib.pyplot as plt
import numpy as np
from IPython import embed #Funktionen importieren
#Zeitpunkte einer EOD über Zero-crossings finden, die in einer Steigung liegen
data_dir = "../data"
dataset = "2018-11-09-ad-invivo-1"
time,eod = read_baseline_eod(os.path.join(data_dir, dataset))
spike_times = read_baseline_spikes(os.path.join(data_dir, dataset))
print(len(spike_times))
eod_times = zero_crossing(eod,time)
eod_durations = np.diff(eod_times)
print(len(spike_times))
print(len(eod_durations))
#for st in spike_times:
#et = eod_times[eod_times < st]
#dt = st - et
#vs = vector_strength(spike_times, eod_durations)