% SONAR_MATRIX - Script for processing sonar data
deployment = 2;

% for Eurostrat first deployment 2002-2003
if(deployment == 1),
  FirstSonarDay = julian(2002,11,7,0);
  LastSonarDay =  julian(2002,12,8,0);
  dnFirst = datenum(2002,11,7);
  dnLast  = datenum(2002,12,8);
  rot = 180-17+43; % correct orientation for ES dep 1
  %RootDataDir = 'E:\EuroStrat\701im\test_images_for_movies_single\'; % don't delete these
  %RootDataDir = 'E:\EuroStrat\701im\test_images_for_movies\'; % don't delete these
  %RootDataDir = 'E:\EuroStrat\701im\images_for_movies\'; % don't delete these
  RootDataDir = 'C:\crs\data\EuroStrat_0203\701-Particle\images_for_movies\';
  %RootDataDir = 'C:\projects\euro\701im\test_images_for_movies_single\';
  %RootDataDir = 'G:\Feb03EuroStrat\701-Particle\701im\';
  %RootDataDir = 'C:\crs\data\EuroStrat_0203\701-Particle\701im\';
  FanSettings = {'348','0','1','500','0.55','2002'};
elseif(deployment == 2),
  FirstSonarDay = julian(2003,2,19,0);
  LastSonarDay = julian(2003,4,21,0);
  dnFirst = datenum(2003,2,19)
  dnLast  = datenum(2003,4,21)
  rot = 180-17+(276-360) % this produced correct orientation for ES dep 2
  RootDataDir = 'C:\crs\data\EuroStrat_0503\711-Particle\711_IM_Movies';
  %RootDataDir = 'C:\crs\data\EuroStrat_0503\711-Particle\Short_movie_test';
  FanSettings = {'348','0','1','500','0.55','2003'};
else
  error('Bad deployment number')
end

settings(1).USGS_mooring_number = 7012;
settings(1).filename = 'none';
settings(1).jd = 'NaN';
settings(1).rot = rot;
settings(1).Sonar_Azimuth = 0.;
settings(1).SectorSweep = 348;
settings(1).AngleSweepAround = 0;
settings(1).StepSize = 1;
settings(1).DataPoints = 500;           % number of 'ranges' recorded in each ping
settings(1).Height = .55; 
settings(1).datayear = 2002;
settings(1).DegPerStep = 'NaN';
settings(1).SonarRange = 5;
if(deployment==2),
  settings(1).USGS_mooring_number = 7112;
  settings(1).datayear = 2003;
end

%avim = zeros(2320,449);
Fanidx = 1;
SubDataDir = dir(RootDataDir);
% go through each and parse
for d=3:length(SubDataDir), % the first two are always . and ..
  disp(SubDataDir(d).name)
  % act only on directories
  if exist(fullfile(RootDataDir,SubDataDir(d).name)) == 7,
    DataFiles = dir(fullfile(RootDataDir,SubDataDir(d).name))
    length(DataFiles)
    for f = 3:length(DataFiles),
      % disp(DataFiles(f).name)
      % name convention is S1111020.F50 = S1MMDDHH.?MM 
      % where ? = F for fan or P for pencil
      a = sscanf(DataFiles(f).name,'S1%c%c%c%c%c%c.%c%c%c')
      mon = str2num(a(1:2));
      day = str2num(a(3:4));
      hr = str2num(a(5:6));
      min = str2num(a(8:9));
      %if(deployment==1),this_year = 2002;,end
      %if(deployment==2),this_year = 2003;,end   
      %datenum_sonar = datenum(this_year,mon,day,hr,min,0);
      WhichSonar = a(7);
      if upper(WhichSonar) == 'F',
	settings(Fanidx) = settings(1);
	%FanDatenum=datenum([settings.datayear mon day hr min 0]);
	fn = fullfile(RootDataDir,SubDataDir(d).name,DataFiles(f).name)
        settings(Fanidx).filename = fn;
        settings(Fanidx).jd = julian([settings(Fanidx).datayear mon day hr min 0]); 
	clf
	[imagedata, Xplot,Yplot,Zi, settings(Fanidx)]=...
	   showfan(fn,'none',settings(Fanidx));

	drawnow;
	shg
	pause
	pfn = sprintf('./dep2_frames/frame%03d.png',Fanidx);
	eval(['print -dpng ',pfn]) 
	Fanidx = Fanidx+1;
      end
    end
  end
end