%For bugs, comments, etc. please contact Kapil Krishan :
%kkrishan@uci.edu
%www.physics.uci.edu/~foams
%version: I

clear all;close all;

filename='image0.bmp';pathname='../images/';
[P_0,filename,pathname,Params]=new_centers_list(filename,pathname,struct([]),'ON');

filename='image1.bmp';pathname='../images/';
[P_1,filename,pathname,Params]=new_centers_list(filename,pathname,Params,'ON');

reg_of_motion=10*[1 1];
mapping=two_frame_mapping(P_0,P_1,reg_of_motion);

velocities_mapping=P_1(mapping(:,2),[2 3])-P_0(mapping(:,1),[2 3]);
quiver(P_1(mapping(:,2),2),P_1(mapping(:,2),3),velocities_mapping(:,1),velocities_mapping(:,2));
axis equal;axis image;