def parfor_calc_pixel(arr, idx, r, ecs_list, pt, tm, rs, dojit=False): sr = ec2b.calc_pixel(idx[1], idx[0], r, ecs_list, pt, tm, rs, dojit) return sr.data
def parfor_calc_pixel(arr, idx, r, ecs_list, pt, tm, rs, dojit=False): sr=ec2b.calc_pixel(idx[1], idx[0], r, ecs_list, pt, tm, rs, dojit) return sr.data
interphase_dur=.5 / 1000., tsample=.1 / 1000.) else: ptrain = e2cm.Psycho2Pulsetrain(current_amplitude=0, dur=0.25, pulse_dur=.5 / 1000., interphase_dur=.5 / 1000., tsample=.1 / 1000.) pt.append(ptrain) [ecs_list, cs_list] = r.electrode_ecs(e_all) tm1 = ec2b.TemporalModel() rs = 1 / (fps * pt[0].tsample) sr_tmp = ec2b.calc_pixel(0, 0, r, ecs_list, pt, tm1, rs, dojit=False) brightness_movie = np.zeros( (r.gridx.shape[0], r.gridx.shape[1], sr_tmp.shape[0])) #brightness_movie = np.zeros((4, 4, sr_tmp.shape[0])) def parfor_calc_pixel(arr, idx, r, ecs_list, pt, tm, rs, dojit=False): sr = ec2b.calc_pixel(idx[1], idx[0], r, ecs_list, pt, tm, rs, dojit) return sr.data brightness_movie = utils.parfor(brightness_movie, parfor_calc_pixel, r, ecs_list, pt, tm1,
for ct, rf in enumerate(e_rf): if ct==ee: ptrain=e2cm.Psycho2Pulsetrain(current_amplitude=3, dur=0.25, pulse_dur=.5/1000., interphase_dur=.5/1000., tsample=.1/1000.) else: ptrain=e2cm.Psycho2Pulsetrain(current_amplitude=0, dur=0.25, pulse_dur=.5/1000., interphase_dur=.5/1000., tsample=.1/1000.) pt.append(ptrain) [ecs_list, cs_list] = r.electrode_ecs(e_all) tm1 = ec2b.TemporalModel() rs=1/(fps*pt[0].tsample) sr_tmp=ec2b.calc_pixel(0, 0, r, ecs_list, pt, tm1, rs, dojit=False) brightness_movie = np.zeros((r.gridx.shape[0], r.gridx.shape[1], sr_tmp.shape[0])) #brightness_movie = np.zeros((4, 4, sr_tmp.shape[0])) def parfor_calc_pixel(arr, idx, r, ecs_list, pt, tm, rs, dojit=False): sr=ec2b.calc_pixel(idx[1], idx[0], r, ecs_list, pt, tm, rs, dojit) return sr.data brightness_movie = utils.parfor(brightness_movie, parfor_calc_pixel, r, ecs_list, pt, tm1, rs, dojit=False, n_jobs=12, axis=-1) # brightnessmovie[yy, xx, :] = sr_rs filename='SE_' + retinaname + '_E' + str(ee) np.save(filename, brightness_movie) scipy.misc.imsave(filename +'.jpg', brightness_movie[:, :, 10]) # moviefilename='singleelectrode_' + retinaname + str(ee)