source: altifloat/matlab_toolbox/twinexp4f.m @ 160

% 1 floater 2 pts
figure
ls
load obs_float_test.dat
x1=obs_float_test (1:2:40,4)
y1=obs_float_test (2:2:40,4)
x1=x1+1;
y1=y1+1;
plot(x1(1),y1(1),'+')
hold on
plot(x1(20),y1(20),'>')
hold on
plot(x1(1:20),y1(1:20),'-r')
hold on
x2=obs_float_test (41:2:80,4)
y2=obs_float_test (42:2:80,4)
x2=x2+1;
y2=y2+1;
plot(x2(1),y2(1),'+')
hold on
plot(x2(20),y2(20),'>')
hold on
plot(x2(1:20),y2(1:20),'-r')
hold on
x3=obs_float_test (81:2:120,4)
y3=obs_float_test (82:2:120,4)
x3=x3+1;
y3=y3+1;
plot(x3(1),y3(1),'+')
hold on
plot(x3(20),y3(20),'>')
hold on
plot(x3(1:20),y3(1:20),'-r')
hold on
x4=obs_float_test (121:2:160,4)
y4=obs_float_test (122:2:160,4)
x4=x4+1;
y4=y4+1;
plot(x4(1),y4(1),'+')
hold on
plot(x4(20),y4(20),'>')
hold on
plot(x4(1:20),y4(1:20),'-r')
hold on



load rfloat_total.dat;
whos;
rfloat_total(:,:)
x5=rfloat_total(1:20,4)
y5=rfloat_total(1:20,3)
x5=x5+1;
y5=y5+1;
plot(x5(2),y5(2),'+')
hold on
plot(x5(19),y5(19),'>')
hold on
plot(x5(2:19),y5(2:19),'-g')
hold on
x6=rfloat_total(21:40,4)
y6=rfloat_total(21:40,3)
x6=x6+1;
y6=y6+1;
plot(x6(2),y6(2),'+')
hold on
plot(x6(19),y6(19),'>')
hold on
plot(x6(2:19),y6(2:19),'-g')
hold on
ls;
x7=rfloat_total(41:60,4)
y7=rfloat_total(41:60,3)
x7=x7+1;
y7=y7+1;
plot(x7(2),y7(2),'+')
hold on
plot(x7(19),y7(19),'>')
hold on
plot(x7(2:19),y7(2:19),'-g')
hold on
ls;
x8=rfloat_total(61:80,4)
y8=rfloat_total(61:80,3)
x8=x8+1;
y8=y8+1;
plot(x8(2),y8(2),'+')
hold on
plot(x8(19),y8(19),'>')
hold on
plot(x8(2:19),y8(2:19),'-g')
hold on
ls;


load uv_total.dat
X=uv_total(:,end-3);
Y=uv_total(:,end-2);
X=X+1;
Y=Y+1;
U1=uv_total(:,end-1);
V1=uv_total(:,end);
quiver(X,Y,U1,V1,'r');


load vel_bck_all.dat;
whos;
X=vel_bck_all(:,end-3);
Y=vel_bck_all(:,end-2);
U2=vel_bck_all(:,end-1);
V2=vel_bck_all(:,end);
quiver(X,Y,U2,V2,'b');

load vel_next_all.dat;
whos;
X=vel_next_all(:,end-3);
Y=vel_next_all(:,end-2);
U3=vel_next_all(:,end-1);
V3=vel_next_all(:,end);
quiver(X,Y,U3,V3,'g');
figure

it=find(vel_next_all(:,1)==1);
T1=reshape(vel_next_all(it,4),58,87);
it=find(vel_next_all(:,1)==2);
T2=reshape(vel_next_all(it,4),58,87);
it=find(vel_next_all(:,1)==3);
T3=reshape(vel_next_all(it,4),58,87);
it=find(vel_next_all(:,1)==4);
T4=reshape(vel_next_all(it,4),58,87);
it=find(vel_next_all(:,1)==5);
T5=reshape(vel_next_all(it,4),58,87);
it=find(vel_next_all(:,1)==6);
T6=reshape(vel_next_all(it,4),58,87);
it=find(vel_next_all(:,1)==7);
T7=reshape(vel_next_all(it,4),58,87);
it=find(vel_next_all(:,1)==8);
T8=reshape(vel_next_all(it,4),58,87);
it=find(vel_next_all(:,1)==9);
T9=reshape(vel_next_all(it,4),58,87);
it=find(vel_next_all(:,1)==10);
T10=reshape(vel_next_all(it,4),58,87);
it=find(vel_next_all(:,1)==11);
T11=reshape(vel_next_all(it,4),58,87);
it=find(vel_next_all(:,1)==12);
T12=reshape(vel_next_all(it,4),58,87);
it=find(vel_next_all(:,1)==13);
T13=reshape(vel_next_all(it,4),58,87);
it=find(vel_next_all(:,1)==14);
T14=reshape(vel_next_all(it,4),58,87);
it=find(vel_next_all(:,1)==15);
T15=reshape(vel_next_all(it,4),58,87);
it=find(vel_next_all(:,1)==16);
T16=reshape(vel_next_all(it,4),58,87);
it=find(vel_next_all(:,1)==17);
T17=reshape(vel_next_all(it,4),58,87);
it=find(vel_next_all(:,1)==18);
T18=reshape(vel_next_all(it,4),58,87);
it=find(vel_next_all(:,1)==19);
T19=reshape(vel_next_all(it,4),58,87);

T=(T1+T2+T3+T4+T5+T6+T7+T8+T9+T10+T11+T12+T13+T14+T15+T16+T17+T18+T19)/19;

it=find(vel_bck_all(:,1)==1);
A1=reshape(vel_bck_all(it,4),58,87);
it=find(vel_bck_all(:,1)==2);
A2=reshape(vel_bck_all(it,4),58,87);
it=find(vel_bck_all(:,1)==3);
A3=reshape(vel_bck_all(it,4),58,87);
it=find(vel_bck_all(:,1)==4);
A4=reshape(vel_bck_all(it,4),58,87);
it=find(vel_bck_all(:,1)==5);
A5=reshape(vel_bck_all(it,4),58,87);
it=find(vel_bck_all(:,1)==6);
A6=reshape(vel_bck_all(it,4),58,87);
it=find(vel_bck_all(:,1)==7);
A7=reshape(vel_bck_all(it,4),58,87);
it=find(vel_bck_all(:,1)==8);
A8=reshape(vel_bck_all(it,4),58,87);
it=find(vel_bck_all(:,1)==9);
A9=reshape(vel_bck_all(it,4),58,87);
it=find(vel_bck_all(:,1)==10);
A10=reshape(vel_bck_all(it,4),58,87);
it=find(vel_bck_all(:,1)==11);
A11=reshape(vel_bck_all(it,4),58,87);
it=find(vel_bck_all(:,1)==12);
A12=reshape(vel_bck_all(it,4),58,87);
it=find(vel_bck_all(:,1)==13);
A13=reshape(vel_bck_all(it,4),58,87);
it=find(vel_bck_all(:,1)==14);
A14=reshape(vel_bck_all(it,4),58,87);
it=find(vel_bck_all(:,1)==15);
A15=reshape(vel_bck_all(it,4),58,87);
it=find(vel_bck_all(:,1)==16);
A16=reshape(vel_bck_all(it,4),58,87);
it=find(vel_bck_all(:,1)==17);
A17=reshape(vel_bck_all(it,4),58,87);
it=find(vel_bck_all(:,1)==18);
A18=reshape(vel_bck_all(it,4),58,87);
it=find(vel_bck_all(:,1)==19);
A19=reshape(vel_bck_all(it,4),58,87);

A=(A1+A2+A3+A4+A5+A6+A7+A8+A9+A10+A11+A12+A13+A14+A15+A16+A17+A18+A19)/19;
F=A.^2-T.^2
imagesc(F)

figure
it=find(uv_total(:,1)==1);
I1=reshape(uv_total(it,4),58,87);
it=find(uv_total(:,1)==2);
I2=reshape(uv_total(it,4),58,87);
it=find(uv_total(:,1)==3);
I3=reshape(uv_total(it,4),58,87);
it=find(uv_total(:,1)==4);
I4=reshape(uv_total(it,4),58,87);
it=find(uv_total(:,1)==5);
I5=reshape(uv_total(it,4),58,87);
it=find(uv_total(:,1)==6);
I6=reshape(uv_total(it,4),58,87);
it=find(uv_total(:,1)==7);
I7=reshape(uv_total(it,4),58,87);
it=find(uv_total(:,1)==8);
I8=reshape(uv_total(it,4),58,87);
it=find(uv_total(:,1)==9);
I9=reshape(uv_total(it,4),58,87);
it=find(uv_total(:,1)==10);
I10=reshape(uv_total(it,4),58,87);
it=find(uv_total(:,1)==11);
I11=reshape(uv_total(it,4),58,87);
it=find(uv_total(:,1)==12);
I12=reshape(uv_total(it,4),58,87);
it=find(uv_total(:,1)==13);
I13=reshape(uv_total(it,4),58,87);
it=find(uv_total(:,1)==14);
I14=reshape(uv_total(it,4),58,87);
it=find(uv_total(:,1)==15);
I15=reshape(uv_total(it,4),58,87);
it=find(uv_total(:,1)==16);
I16=reshape(uv_total(it,4),58,87);
it=find(uv_total(:,1)==17);
I17=reshape(uv_total(it,4),58,87);
it=find(uv_total(:,1)==18)
I18=reshape(uv_total(it,4),58,87);
it=find(uv_total(:,1)==19)
I19=reshape(uv_total(it,4),58,87);

I=(I1+I2+I3+I4+I5+I6+I7+I8+I9+I10+I11+I12+I13+I14+I15+I16+I17+I18+I19)/19;
F1=A.^2-I.^2;
imagesc(F1)

F2=I.^2-T.^2;
figure
imagesc(F2)