35 lines
1.0 KiB
Python
35 lines
1.0 KiB
Python
import nixio as nix
|
|
import os
|
|
from IPython import embed
|
|
|
|
def read_baseline_eod(dataset):
|
|
base = dataset.split(os.path.sep)[-1] + ".nix"
|
|
nix_file = nix.File.open(os.path.join(dataset, base), nix.FileMode.ReadOnly)
|
|
b = nix_file.blocks[0]
|
|
t = b.tags["BaselineActivity_1"]
|
|
eod_da = b.data_arrays["LocalEOD-1"]
|
|
eod = t.retrieve_data("LocalEOD-1")
|
|
time = eod_da.dimensions[0].axis(len(eod))
|
|
nix_file.close()
|
|
return time, eod
|
|
|
|
|
|
def read_baseline_spikes(dataset):
|
|
base = dataset.split(os.path.sep)[-1] + ".nix"
|
|
nix_file = nix.File.open(os.path.join(dataset, base), nix.FileMode.ReadOnly)
|
|
b = nix_file.blocks[0]
|
|
t = b.tags["BaselineActivity_1"]
|
|
embed()
|
|
spike_times = t.retrieve_data("Spikes-1")
|
|
|
|
nix_file.close()
|
|
return spike_times
|
|
|
|
|
|
if __name__ == "__main__":
|
|
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))
|
|
|