def rastercountplot(spiketrain, nbins=50, **kwargs):
fig = plt.figure(figsize=(14, 6))
gs = gridspec.GridSpec(2, 1, hspace=0.01, height_ratios=[0.2,0.8])
ax1 = plt.subplot(gs[0])
ax2 = plt.subplot(gs[1])
color = kwargs.get('color', None)
if color is None:
color = '0.4'
ds = (spiketrain.support.stop - spiketrain.support.start)/nbins
flattened = spiketrain.bin(ds=ds).flatten()
steps = np.squeeze(flattened.data)
stepsx = np.linspace(spiketrain.support.start, spiketrain.support.stop, num=flattened.n_bins)
# ax1.plot(stepsx, steps, drawstyle='steps-mid', color='none');
ax1.set_ylim([-0.5, np.max(steps)+1])
rasterplot(spiketrain, ax=ax2, **kwargs)
utils.clear_left_right(ax1)
utils.clear_top_bottom(ax1)
utils.clear_top(ax2)
ax1.fill_between(stepsx, steps, step='mid', color=color)
utils.sync_xlims(ax1, ax2)
return ax1, ax2
评论列表
文章目录