'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')