vz2.add_layer(maps, [cov2, 'h', '-v', mz2], ac_display_options)
vz2.add_hmark(mx2)
vz2.add_vmark(my2)
cov_display_options = {}
cov_display_options['interpolation'] = 'nearest'
cov_display_options['cmap'] = pl.cm.RdBu
cov_display_options['vmax'] = 1.
cov_display_options['vmin'] = -1.
vcov = MatshowView(ax_cov, -cov_mat, ['h', 'v'], cov_display_options)
vcov.add_hmark(cov2)
vcov.add_vmark(cov1)
vseries = PlotView(ax_series,
time_series, [cov1, 'h'], {'color': 'b'},
autoscale=True)
ax_series.set_ylim(time_series.min(), time_series.max())
vseries.add_layer(time_series, [cov2, 'h'], {'color': 'r'})
vx1.draw()
vy1.draw()
vz1.draw()
vx2.draw()
vy2.draw()
vz2.draw()
vcov.draw()
vseries.draw()
pl.show()
# Views
vx = ImshowView(ax_x, data, [mx, 'h', '-v', mn], display_options)
vx.add_hmark(my)
vx.add_vmark(mz)
vy = ImshowView(ax_y, data, ['h', my, '-v', mn], display_options)
vy.add_hmark(mx)
vy.add_vmark(mz)
vz = ImshowView(ax_z, data, ['h', '-v', mz, mn], display_options)
vz.add_hmark(mx)
vz.add_vmark(my)
vn = PlotView(ax_n, data, [mx, my, mz, 'h'], {'color': 'gray'})
pl.ylim(0, np.max(data))
vn.add_hmark(mn)
ac_display_options = {}
ac_display_options['interpolation'] = 'nearest'
ac_display_options['cmap'] = pl.cm.autumn
vx.add_layer(data_act, [mx, 'h', '-v', mn], display_options=ac_display_options)
vy.add_layer(data_act, ['h', my, '-v', mn], display_options=ac_display_options)
vz.add_layer(data_act, ['h', '-v', mz, mn], display_options=ac_display_options)
vn.add_layer(data_act, [mx, my, mz, 'h'], display_options={'color': 'r'})
vx.draw()
vy.draw()
vz.draw()
