Neuroinflab / kCSD-python

Kernel Current Source Density
https://kcsd-python.readthedocs.io/en/latest/
Other
15 stars 23 forks source link

subplots in figures/npx/traub_data_kcsd_column_figure.py are for different time points #136

Closed mczerwinski closed 3 years ago

mczerwinski commented 3 years ago

In make_column_plot function: 1st, 2nd and 4th subplots show all of the currents 3rd - potentials, show potentials in moment defined by the variable time_pt_interest - that is 3000, in the middle of 6000 long signal.

ax1 = plt.subplot(241, aspect='equal')
plot_all_currents(ax1, xmin, xmax, ymin, ymax, all_x, all_y, all_z, all_val, letter='A')

ax2 = plt.subplot(242, aspect='equal') plot_csd_slice(ax2, xmin, xmax, ymin, ymax, all_x, all_y, all_z, all_val, letter='B')

ax3 = plt.subplot(243, aspect='equal') plot_dense_potentials(ax3, h, pop_names, time_pts, time_pt_interest, letter='C', filt=False)

ax4 = plt.subplot(244, aspect='equal') xx, yy = np.mgrid[xmin:xmax:np.complex(0, true_csd.shape[0]), ymin:ymax:np.complex(0, true_csd.shape[1])] plot_csd_smooth(ax4, xmin, xmax, ymin, ymax, true_csd[:, :], xx, yy, letter='D')

Later, in the same function:

plot_csd_smooth(ax_list1[i], xmin, xmax, ymin, ymax, kcsd[:, :, 750], x, y, letter=letter_list1[i], ele_pos=ele_pos)

The kcsd is plotted in time point of 750.

mczerwinski commented 3 years ago

It's because potentials are downsampled 4 times. Having values written plain was hard to keep track of it all.