Afile='';
Bfile=''; 


%  6 passes starting with 64*64 windows and ending with 16*16.
winsize=[64 64;64 64;32 32;32 32;16 16;16 16];
winsize=[64 64;32 32;32 32;16 16;16 16];
winsize=[64 64;64 64;32 32;32 32];
winsize=[64 64;64 64];

% 5 passes ending up with 32*16 pixels large interrogation regions
%winsize=[64 64;64 32;32 32;32 16;32 16];
%winsize=[32 32;32 32;16 16;16 16;8 8;8 8];
delT=0.078;
overlap=0.75;
method='multin';
wocofile='worldco1.mat';
%wocofile='worldco4.mat';
%wocofile='';
mask='polymask.mat';





[filelist]=textread('filelist.txt','%[^\n]');

for pair=1:2:(size(filelist))

Afile=cell2mat(filelist(pair))
Bfile=cell2mat(filelist(pair+1))
%Cfile=cell2mat(filelist(pair+2))


[Atick level]=strread(Afile,'frame.%u.%[^.].jpg');
[Btick level]=strread(Bfile,'frame.%u.%[^.].jpg');
%[Ctick level]=strread(Cfile,'frame.%u.%[^.].jpg');
delT1=(Btick-Atick)/1000;
%delT2=Ctick-Atick;
%delT2=(Ctick-Btick)/1000;


  [x1,y1,u1,v1,snr1,pkh1]=matpiv(Afile,Bfile,winsize,delT1,...
     overlap,method,wocofile,mask);
%  [x2,y2,u2,v2,snr2,pkh2]=matpiv(Bfile,Cfile,winsize,delT2,...
%     overlap,method,wocofile,mask);
%[x,y,u,v]=matpiv(Afile,Bfile,winsize,delT,overlap,method);
%[su,sv]=snrfilt(x,y,u,v,snr,1.3);
%[pu,pv]=peakfilt(x,y,su,sv,pkh,0.5);
%[gu,gv]=globfilt(x,y,pu,pv,3);
%[mu,mv]=localfilt(x,y,gu,gv,2,'median',3,'polymask.mat');
%[fu,fv]=naninterp(mu,mv,'linear','polymask.mat',x,y);



field=(pair+1)/2

if 1

SNR(:,:,field)=snr1;

[su1,sv1]=snrfilt(x1,y1,u1,v1,snr1,1.2);
SU1(:,:,field)=su1;
SV1(:,:,field)=sv1;
                                                                                
%[su2,sv2]=snrfilt(x2,y2,u2,v2,snr2,1.3);
%SU2(:,:,field)=su2;
%SV2(:,:,field)=sv2;

if 1                                                                                 
[pu1,pv1]=peakfilt(x1,y1,su1,sv1,pkh1,0.3);
PU1(:,:,field)=pu1;
PV1(:,:,field)=pv1;
                                                                                
%[pu2,pv2]=peakfilt(x2,y2,su2,sv2,pkh2,0.4);
%PU2(:,:,field)=pu2;
%PV2(:,:,field)=pv2;
else
pu1=su1;
pv1=sv1;
pu2=su2;
pv2=sv2;
end
                                                                                
[gu1,gv1]=globfilt(x1,y1,pu1,pv1,4);
GU1(:,:,field)=gu1;
GV1(:,:,field)=gv1;
                                                                                
XX1(:,:,field)=x1;
YY1(:,:,field)=y1;
                                                                                
%[gu2,gv2]=globfilt(x2,y2,pu2,pv2,4);
%GU2(:,:,field)=gu2;
%GV2(:,:,field)=gv2;

if 0                                                                                
[mu1,mv1]=localfilt(x1,y1,gu1,gv1,2,'median',3,mask);
MU1(:,:,field)=mu1;
MV1(:,:,field)=mv1;
                                                                                
[mu2,mv2]=localfilt(x2,y2,gu2,gv2,2,'median',3,mask);
MU2(:,:,field)=mu2;
MV2(:,:,field)=mv2;
                                                                                
[fu1,fv1]=naninterp(mu1,mv1,'linear',mask,x1,y1);
FU1(:,:,field)=fu1;
FV1(:,:,field)=fv1;

[fu2,fv2]=naninterp(mu2,mv2,'linear',mask,x2,y2);
FU2(:,:,field)=fu2;
FV2(:,:,field)=fv2;
end;

end

end;

save errsuite.mat GU1 GV1 XX1 YY1;