continue # print 'Doing ', exp dat = load_PopData(exp) dat_c = dat['dat_c'].mean(2) dat_c = normalize(dat_c, axis=0) dat_w = dat['dat_w'].mean(2) dat_w = normalize(dat_w, axis=0) active = dat['active'] d = np.where(active[:, 1])[0] lum_mask, con_mask, flow_mask, four_mask, four_mask_shape,\ freq_mask, orient_mask,\ lum_surr, con_surr, flow_surr, four_surr, four_surr_shape,\ freq_surr, orient_surr,\ lum_whole, con_whole, flow_whole, four_whole, four_whole_shape,\ freq_whole, orient_whole = load_parsed_movie_dat(exp, 'POP', None) all_dat = {'Data': {'Centre': dat_c, 'Whole': dat_w}, 'Movie': {}} all_dat['Movie']['Contrast'] = con_mask all_dat['Movie']['Luminence'] = lum_mask all_dat['Movie']['Fourier'] = np.append(four_mask.real, four_mask.imag, axis=2).astype(np.float) all_dat['Movie']['Frequency'] = freq_mask all_dat['Movie']['Orientation'] = orient_mask for d in all_dat['Data']: # print d X = all_dat['Data'][d] for m in all_dat['Movie']: # print m ys = all_dat['Movie'][m]
def get_mov_data(comb, targ_type, src_type, e, cellid, exp_type, four_downsample=None, shift=0, randomise=None): lum_mask, con_mask, flow_mask, four_mask, four_mask_shape,\ freq_mask, orient_mask,\ lum_surr, con_surr, flow_surr, four_surr, four_surr_shape,\ freq_surr, orient_surr,\ lum_whole, con_whole, flow_whole, four_whole, four_whole_shape,\ freq_whole, orient_whole\ = load_parsed_movie_dat(cellid, exp_type, four_downsample) flow_mask = bin_flow(flow_mask) flow_surr = bin_flow(flow_surr) flow_whole = bin_flow(flow_whole) idx_bar = [] idx_flow = [] idx_four = [] idx_freq = [] idx_orient = [] lbl_bar = [] four_shape = [] all_dat = None if randomise is None: if targ_type == 'Center': source = e['psth_c_shift'][shift] elif targ_type == 'Surround': source = e['psth_s_shift'][shift] elif targ_type == 'Whole': source = e['psth_w_shift'][shift] elif randomise == 'random': if targ_type == 'Center': source = e['psth_c_rand'] elif targ_type == 'Surround': source = e['psth_s_rand'] elif targ_type == 'Whole': source = e['psth_w_rand'] elif randomise == 'generated': if targ_type == 'Center': source = e['psth_c_gen'] elif targ_type == 'Surround': source = e['psth_s_gen'] elif targ_type == 'Whole': source = e['psth_w_gen'] source = source.mean(0) edge = e['edge'] if src_type == 'Center': if 'Luminance' in comb: for l in lum_mask: all_dat = append_Nones(all_dat, l[:, np.newaxis], 1) idx_bar += [all_dat.shape[1] - 1] lbl_bar += ['Luminance'] if 'Contrast' in comb: for c in con_mask: all_dat = append_Nones(all_dat, c[:, np.newaxis], 1) idx_bar += [all_dat.shape[1] - 1] lbl_bar += ['Contrast'] if 'Flow' in comb: for f in flow_mask: if all_dat != None: pre_len = all_dat.shape[1] else: pre_len = 0 all_dat = append_Nones(all_dat, f[:, :-1], 1) idx_flow += [range(pre_len, all_dat.shape[1])] for f in flow_mask: all_dat = append_Nones(all_dat, f[:, -1:], 1) idx_bar += [all_dat.shape[1] - 1] lbl_bar += ['Flow Vel'] if 'Fourier' in comb: for f in four_mask: if all_dat != None: pre_len = all_dat.shape[1] else: pre_len = 0 f = np.append(f.real, f.imag, 1) all_dat = append_Nones(all_dat, f, 1) idx_four += [range(pre_len, all_dat.shape[1])] four_shape.append(four_mask_shape) if 'Frequency' in comb: for f in freq_mask: if all_dat != None: pre_len = all_dat.shape[1] else: pre_len = 0 all_dat = append_Nones(all_dat, f, 1) idx_freq += [range(pre_len, all_dat.shape[1])] if 'Orientation' in comb: for f in orient_mask: if all_dat != None: pre_len = all_dat.shape[1] else: pre_len = 0 all_dat = append_Nones(all_dat, f, 1) idx_orient += [range(pre_len, all_dat.shape[1])] elif src_type == 'Surround': if 'Luminance' in comb: for l in lum_surr: all_dat = append_Nones(all_dat, l[:, np.newaxis], 1) idx_bar += [all_dat.shape[1] - 1] lbl_bar += ['Luminance'] if 'Contrast' in comb: for c in con_surr: all_dat = append_Nones(all_dat, c[:, np.newaxis], 1) idx_bar += [all_dat.shape[1] - 1] lbl_bar += ['Contrast'] if 'Flow' in comb: for f in flow_surr: if all_dat != None: pre_len = all_dat.shape[1] else: pre_len = 0 all_dat = append_Nones(all_dat, f[:, :-1], 1) idx_flow += [range(pre_len, all_dat.shape[1])] for f in flow_surr: all_dat = append_Nones(all_dat, f[:, -1:], 1) idx_bar += [all_dat.shape[1] - 1] lbl_bar += ['Flow Vel'] if 'Fourier' in comb: for f in four_surr: if all_dat != None: pre_len = all_dat.shape[1] else: pre_len = 0 all_dat = append_Nones(all_dat, f, 1) idx_four += [range(pre_len, all_dat.shape[1])] four_shape.append(four_surr_shape) if 'Frequency' in comb: for f in freq_surr: if all_dat != None: pre_len = all_dat.shape[1] else: pre_len = 0 all_dat = append_Nones(all_dat, f, 1) idx_freq += [range(pre_len, all_dat.shape[1])] if 'Orientation' in comb: for f in orient_surr: if all_dat != None: pre_len = all_dat.shape[1] else: pre_len = 0 all_dat = append_Nones(all_dat, f, 1) idx_orient += [range(pre_len, all_dat.shape[1])] elif src_type == 'Whole': if 'Luminance' in comb: for l in lum_whole: all_dat = append_Nones(all_dat, l[:, np.newaxis], 1) idx_bar += [all_dat.shape[1] - 1] lbl_bar += ['Luminance'] if 'Contrast' in comb: for c in con_whole: all_dat = append_Nones(all_dat, c[:, np.newaxis], 1) idx_bar += [all_dat.shape[1] - 1] lbl_bar += ['Contrast'] if 'Flow' in comb: for f in flow_whole: if all_dat != None: pre_len = all_dat.shape[1] else: pre_len = 0 all_dat = append_Nones(all_dat, f[:, :-1], 1) idx_flow += [range(pre_len, all_dat.shape[1])] for f in flow_whole: all_dat = append_Nones(all_dat, f[:, -1:], 1) idx_bar += [all_dat.shape[1] - 1] lbl_bar += ['Flow Vel'] if 'Fourier' in comb: for f in four_whole: if all_dat != None: pre_len = all_dat.shape[1] else: pre_len = 0 all_dat = append_Nones(all_dat, f, 1) idx_four += [range(pre_len, all_dat.shape[1])] four_shape.append(four_whole_shape) if 'Frequency' in comb: for f in freq_whole: if all_dat != None: pre_len = all_dat.shape[1] else: pre_len = 0 all_dat = append_Nones(all_dat, f, 1) idx_freq += [range(pre_len, all_dat.shape[1])] if 'Orientation' in comb: for f in orient_whole: if all_dat != None: pre_len = all_dat.shape[1] else: pre_len = 0 all_dat = append_Nones(all_dat, f, 1) idx_orient += [range(pre_len, all_dat.shape[1])] #all_dat = np.tile(all_dat, [source.shape[0], 1, 1]) plot_params = {} plot_params['idx_bar'] = idx_bar plot_params['idx_flow'] = idx_flow plot_params['idx_four'] = idx_four plot_params['lbl_bar'] = lbl_bar plot_params['four_shape'] = four_shape return all_dat, source, plot_params
continue # print 'Doing ', exp dat = load_PopData(exp) dat_c = dat['dat_c'].mean(2) dat_c = normalize(dat_c, axis=0) dat_w = dat['dat_w'].mean(2) dat_w = normalize(dat_w, axis=0) active = dat['active'] d = np.where(active[:, 1])[0] lum_mask, con_mask, flow_mask, four_mask, four_mask_shape,\ freq_mask, orient_mask,\ lum_surr, con_surr, flow_surr, four_surr, four_surr_shape,\ freq_surr, orient_surr,\ lum_whole, con_whole, flow_whole, four_whole, four_whole_shape,\ freq_whole, orient_whole = load_parsed_movie_dat(exp, 'POP', None) all_dat = {'Data': {'Centre': dat_c, 'Whole': dat_w}, 'Movie': {}} all_dat['Movie']['Contrast'] = con_mask all_dat['Movie']['Luminence'] = lum_mask all_dat['Movie']['Fourier'] = np.append(four_mask.real, four_mask.imag, axis=2).astype(np.float) all_dat['Movie']['Frequency'] = freq_mask all_dat['Movie']['Orientation'] = orient_mask for d in all_dat['Data']: # print d X = all_dat['Data'][d] for m in all_dat['Movie']: # print m ys = all_dat['Movie'][m] if len(ys.shape) < 3:
def get_mov_data(targ_type, e, expdate, exp_type, four_downsample=None): lum_mask, con_mask, flow_mask,\ four_mask, four_mask_shape,\ lum_surr, con_surr, flow_surr,\ four_surr, four_surr_shape,\ lum_whole, con_whole, flow_whole,\ four_whole, four_whole_shape \ = load_parsed_movie_dat(expdate, exp_type, four_downsample) flow_mask = bin_flow(flow_mask) flow_surr = bin_flow(flow_surr) flow_whole = bin_flow(flow_whole) idx_bar = [] idx_flow = [] idx_four = [] lbl_bar = [] four_shape = [] all_dat = None if targ_type == 'Center': source = e['psth_c'] for f in four_mask: if all_dat != None: pre_len = all_dat.shape[1] else: pre_len = 0 all_dat = append_Nones(all_dat, f, 1) idx_four += [range(pre_len, all_dat.shape[1])] four_shape.append(four_mask_shape) elif targ_type == 'Surround': source = e['psth_s'] for f in four_surr: if all_dat != None: pre_len = all_dat.shape[1] else: pre_len = 0 all_dat = append_Nones(all_dat, f, 1) idx_four += [range(pre_len, all_dat.shape[1])] four_shape.append(four_surr_shape) elif targ_type == 'Whole': source = e['psth_w'] for f in four_mask: if all_dat != None: pre_len = all_dat.shape[1] else: pre_len = 0 all_dat = append_Nones(all_dat, f, 1) idx_four += [range(pre_len, all_dat.shape[1])] four_shape.append(four_mask_shape) # all_dat = np.tile(all_dat, [source.shape[0], 1, 1]) plot_params = {} plot_params['idx_bar'] = idx_bar plot_params['idx_flow'] = idx_flow plot_params['idx_four'] = idx_four plot_params['lbl_bar'] = lbl_bar plot_params['four_shape'] = four_shape return all_dat, source, plot_params
def get_mov_data(comb, targ_type, src_type, e, expdate, exp_type, four_downsample=None): lum_mask, con_mask, flow_mask,\ four_mask, four_mask_shape,\ lum_surr, con_surr, flow_surr,\ four_surr, four_surr_shape,\ lum_whole, con_whole, flow_whole,\ four_whole, four_whole_shape \ = load_parsed_movie_dat(expdate, exp_type, four_downsample) flow_mask = bin_flow(flow_mask) flow_surr = bin_flow(flow_surr) flow_whole = bin_flow(flow_whole) idx_bar = [] idx_flow = [] idx_four = [] lbl_bar = [] four_shape = [] all_dat = None if targ_type == 'Center': source = e['psth_c'] elif targ_type == 'Surround': source = e['psth_s'] elif targ_type == 'Whole': source = e['psth_w'] if src_type == 'Center': if 'Luminance' in comb: for l in lum_mask: all_dat = append_Nones(all_dat, l[:, np.newaxis], 1) idx_bar += [all_dat.shape[1] - 1] lbl_bar += ['Luminance'] if 'Contrast' in comb: for c in con_mask: all_dat = append_Nones(all_dat, c[:, np.newaxis], 1) idx_bar += [all_dat.shape[1] - 1] lbl_bar += ['Contrast'] if 'Flow' in comb: for f in flow_mask: if all_dat != None: pre_len = all_dat.shape[1] else: pre_len = 0 all_dat = append_Nones(all_dat, f[:, :-1], 1) idx_flow += [range(pre_len, all_dat.shape[1])] for f in flow_mask: all_dat = append_Nones(all_dat, f[:, -1:], 1) idx_bar += [all_dat.shape[1] - 1] lbl_bar += ['Flow Vel'] if 'Fourier' in comb: for f in four_mask: if all_dat != None: pre_len = all_dat.shape[1] else: pre_len = 0 all_dat = append_Nones(all_dat, f, 1) idx_four += [range(pre_len, all_dat.shape[1])] four_shape.append(four_mask_shape) elif src_type == 'Surround': if 'Luminance' in comb: for l in lum_surr: all_dat = append_Nones(all_dat, l[:, np.newaxis], 1) idx_bar += [all_dat.shape[1] - 1] lbl_bar += ['Luminance'] if 'Contrast' in comb: for c in con_surr: all_dat = append_Nones(all_dat, c[:, np.newaxis], 1) idx_bar += [all_dat.shape[1] - 1] lbl_bar += ['Contrast'] if 'Flow' in comb: for f in flow_surr: if all_dat != None: pre_len = all_dat.shape[1] else: pre_len = 0 all_dat = append_Nones(all_dat, f[:, :-1], 1) idx_flow += [range(pre_len, all_dat.shape[1])] for f in flow_surr: all_dat = append_Nones(all_dat, f[:, -1:], 1) idx_bar += [all_dat.shape[1] - 1] lbl_bar += ['Flow Vel'] if 'Fourier' in comb: for f in four_surr: if all_dat != None: pre_len = all_dat.shape[1] else: pre_len = 0 all_dat = append_Nones(all_dat, f, 1) idx_four += [range(pre_len, all_dat.shape[1])] four_shape.append(four_surr_shape) elif src_type == 'Whole': if 'Luminance' in comb: for l in lum_whole: all_dat = append_Nones(all_dat, l[:, np.newaxis], 1) idx_bar += [all_dat.shape[1] - 1] lbl_bar += ['Luminance'] if 'Contrast' in comb: for c in con_whole: all_dat = append_Nones(all_dat, c[:, np.newaxis], 1) idx_bar += [all_dat.shape[1] - 1] lbl_bar += ['Contrast'] if 'Flow' in comb: for f in flow_whole: if all_dat != None: pre_len = all_dat.shape[1] else: pre_len = 0 all_dat = append_Nones(all_dat, f[:, :-1], 1) idx_flow += [range(pre_len, all_dat.shape[1])] for f in flow_whole: all_dat = append_Nones(all_dat, f[:, -1:], 1) idx_bar += [all_dat.shape[1] - 1] lbl_bar += ['Flow Vel'] if 'Fourier' in comb: for f in four_whole: if all_dat != None: pre_len = all_dat.shape[1] else: pre_len = 0 all_dat = append_Nones(all_dat, f, 1) idx_four += [range(pre_len, all_dat.shape[1])] four_shape.append(four_whole_shape) all_dat = np.tile(all_dat, [source.shape[0], 1, 1]) plot_params = {} plot_params['idx_bar'] = idx_bar plot_params['idx_flow'] = idx_flow plot_params['idx_four'] = idx_four plot_params['lbl_bar'] = lbl_bar plot_params['four_shape'] = four_shape return all_dat, source, plot_params