def opto_get_lfp_filtered(opto_dataset,channel,fa,fb,order=4):
'''
Retrieves channel or channels from opto LFP dataset
Channels are 1-indexed
Parameters
----------
opto_dataset : string
path or string identifier for a dataset
channel:
1-indexed channel ID or None to return a NTimes x NChannel array
of all LFP data
fa:
low frequency of band-pass, or 'None' to use a low-pass filter.
if fb is 'None' then this is the cutoff for a high-pass filter.
fb:
high-frequency of band-pass, or 'None to use a high-pass filter.
if fa is 'None' then this is the cutoff for a low-pass filter
'''
Fs = opto_get_Fs(opto_dataset)
if channel is None:
data = metaloadmat(opto_dataset)['LFP'].T
#return array([bandfilter(x,fa,fb,Fs,order) for x in data])
problems = [(i,x,fa,fb,Fs,order) for i,x in enumerate(data)]
data = np.array(parmap(__opto_get_lfp_filtered_helper__,problems))
return squeeze(data)
else:
assert channel>=1
data = metaloadmat(opto_dataset)['LFP'][:,channel-1]
return bandfilter(data,fa,fb,Fs,order)
def opto_get_all_lfp_analytic_quick(opto_dataset,fa,fb):
Fs = 1000.0
order = 4
data = metaloadmat(opto_dataset_compact)['lfp']
data = data.transpose((0,2,1))
data = bandfilter(hilbert(data),fa,fb,Fs,order)
return data
def opto_get_lfp(opto_dataset,channel):
'''
Retrieves channel or channels from opto LFP dataset
Channels are 1-indexed
Parameters
----------
opto_dataset : string
path or string identifier for a dataset
channel:
1-indexed channel ID or None to return a NTimes x NChannel array
of all LFP data
'''
if channel is None:
return metaloadmat(opto_dataset)['LFP'].T
else:
assert channel>=1
return metaloadmat(opto_dataset)['LFP'][:,channel-1]
def opto_get_events_passive(opto_dataset):
start,stop = metaloadmat(opto_dataset)['events']
return start, stop
def opto_get_laser(opto_dataset):
return metaloadmat(opto_dataset)['laser'][0]
def opto_get_Fs(opto_dataset):
return metaloadmat(opto_dataset)['Fs'][0,0]
def opto_get_map(opto_dataset):
return np.array(metaloadmat(opto_dataset_compact)['arrayChannelMap'])
def opto_get_all_lfp_quick(opto_dataset):
data = metaloadmat(opto_dataset_compact)['lfp']
return data
