'DsInh_b': 0., 'DsInh_tauw': 1e9 } if __name__ == '__main__': import sys sys.path.append('../..') import main as ntwk if sys.argv[-1] == 'plot': ## load file data = ntwk.load_dict_from_hdf5('CellRep2019_data.h5') print(data) # ## plot fig, _ = ntwk.raster_and_Vm_plot(data, smooth_population_activity=10.) ntwk.show() else: import numpy as np from analyz.processing.signanalysis import smooth NTWK = ntwk.build_populations(Model, ['RecExc', 'RecInh', 'DsInh'], AFFERENT_POPULATIONS=['AffExc'], with_raster=True, with_Vm=4, verbose=True) ntwk.build_up_recurrent_connections(NTWK, SEED=5, verbose=True) #######################################
# ###################### # ## ----- Plot ----- ## # ###################### # vision_model = ntwk.earlyVis_model(from_file='drifting-grating-data.npz') # ps = ntwk.vision_plot(model=vision_model) # fig0 = ps.protocol_plot() ## load file data = ntwk.load_dict_from_hdf5('visual_input_data.h5') print(data['iRASTER_AffExc']) print(data['tRASTER_AffExc']) # ## plot fig, _ = ntwk.raster_and_Vm_plot(data, ['Exc', 'Inh', 'DsInh', 'AffExc'], smooth_population_activity=10.) ntwk.show() else: if False: params = ntwk.vision_params.copy() params['rf_theta'] = [3 * np.pi / 12., 5. * np.pi / 12.] params['rf_psi'] = [0., 0.] params['rf_freq'] = [0.06, 0.08] params['Ncells'] = 50 vision_model = ntwk.earlyVis_model(params=params) vision_model.full_process('drifting-grating', '') vision_model.save_data('drifting-grating-data.npz')