#%% example play movie
m.playMovie(frate=.03,gain=20.0,magnification=4)
#%% example take first channel of two
channelId=1;
totalChannels=2;
m.makeSubMov(range(channelId-1,m.mov.shape[0],totalChannels))
#%% prtition into two movies up and down
m.crop(crop_top=150,crop_bottom=150,crop_left=150,crop_right=150,crop_begin=0,crop_end=0)
#%% concatenate movies (it will add to the original movie)
# you have to create another movie new_mov=XMovie(...)
m.append(new_mov)
#%% motion correct run 3 times
# WHEN YOU RUN motion_correct YOUR ARE MODIFYING THE OBJECT!!!!
templates=[];
shifts=[];
max_shift=5;
num_iter=10; # numer of times motion correction is executed
template=None # here you can use your own template (best representation of the FOV)
for j in range(0,num_iter):
template_used,shift=m.motion_correct(max_shift=max_shift,template=None,show_movie=False);
templates.append(template_used)
shift=np.asarray(shift)
shifts.append(shift)
