dataTimesCB.m

action = get(gcbo,'Tag');
hfig = get(gcbo,'Parent');
UD = get(hfig,'UserData');
DataFiles = UD.DataFiles;
NIMSinfo = UD.NIMSinfo;
Sites = UD.Sites;
maxSites = UD.MaxSites;
HZ = UD.HZ;
getOS;

dataHome = [ getenv('MTdata') filesep DataFiles{1}.surveyName ];
dataDir = [ dataHome filesep 'data' filesep];
binDir = getenv('MTbin');
header = '';

switch action
   case 'Segment'
      v = get(gcbo,'Value');
      UD = get(gcbo,'UserData');
      if ~isempty(UD.hSeg)
          if v
              set(gcbo,'ForeGroundColor',[1,0,0]);
              if isfield(UD.hSeg,'Color') || isprop(UD.hSeg,'Color')
                set(UD.hSeg,'Color',[1,0,0]);
              end
          else
              set(gcbo,'ForeGroundColor',UD.defClr);
              if isfield(UD.hSeg,'Color') || isprop(UD.hSeg,'Color')
                set(UD.hSeg,'Color',UD.defClr);
              end
          end
      end
      H = findobj('Parent',hfig,'Tag','Segment');
      SelectedDataFiles = {};
      FileNames = {};
      nFiles = 0;
      allFiles = [];
      for l = 1:length(H)
         if get(H(l),'value')
            UD = get(H(l),'UserData');
            nFiles = nFiles+1;
            SelectedDataFiles{nFiles} = DataFiles{UD.iFile};
            if length(UD.Name) == 7
              sName = [UD.Name(3) UD.Name(5:7)];
            else
              sName = UD.Name;
            end
            FileNames{nFiles} = UD.Name;
            allFiles = [allFiles sName ';'];
         end
      end
      %  set text window
      if(length(allFiles) > 0)
         allFiles = allFiles(1:end-1);
      end
      set(findobj('Parent',hfig,'Tag','FileNames'),...
		'String',allFiles,'UserData',SelectedDataFiles);
      % Enable buttons checks number of files in selected group,
      %   enables shift, proecessing buttons accordingly
      SelectedSites = uniqueSites(SelectedDataFiles);
      EnableButtons

   case 'decimated'
      %  no need to do anything for this case ... just
      %    check value before loading files
   case 'Load'
      set(hfig,'Pointer','Watch')
      %  make list of selected file names
      H = findobj('Parent',hfig,'Tag','Segment');

      LDF = {};
      nFiles = 0;
      scaleFac = [];
      rescaleData = 1;
      for l = 1:length(H)
         if get(H(l),'value')
            UD = get(H(l),'UserData');
            nFiles = nFiles+1;
            LDF{nFiles} = DataFiles{UD.iFile};
            LDF{nFiles}.FileName = UD.Name;
            if length(LDF{nFiles}.scaleFactors)>0
               scaleFac = [scaleFac LDF{nFiles}.scaleFactors];
            else
               rescaleData = 0;
            end
         end
      end

      plotStart = get(findobj('Tag','Plot','Parent',hfig),'UserData');
      DEC = get(findobj('Tag','decimated','Parent',hfig),'value') ;
      if(DEC)
         for k = 1:length(LDF)
            FileNames{k} =  [LDF{k}.FileName '.dbn'];
            LDF{k}.FileName = FileNames{k};
         end
      else
         for k = 1:length(LDF)
            FileNames{k} = [LDF{k}.FileName '.bin'];
            LDF{k}.FileName = FileNames{k};
         end
      end
      %  add list of selected file names to list of LOADED
      LOADED = get(gcbo,'UserData');
      nLoaded = length(LOADED);
      if nLoaded > 0
         LOADEDnames = get(findobj('Tag','Plot List','Parent',hfig),'String');
      else
         LOADEDnames = {};
      end
      iA = nLoaded+1 ;
      LOADEDnames{iA} = get(findobj('Parent',hfig,'Tag','FileNames'),...
		'String');
      set(findobj('Tag','Plot List','Parent',hfig),...
		'String',LOADEDnames,'Enable','on','Value',iA);
      set(findobj('Tag','Plot','Parent',hfig),'Enable','on')

      %   load NIMS data  ...  no decimation, scale data to nT, mV
      %    (NOT mV/km)
      decFac = 1;
      SCALE = 1;
      currentDir = pwd;
      eval(['cd ''' dataDir '''']);
      [y,SysTF,dt,t0,dataZero,ch,sta,error] = readTSnims(FileNames,decFac,SCALE);
      eval(['cd ''' currentDir '''']);

      %  rescale data channels to physical units
      if rescaleData
         [nch,npts] = size(y);
         for k = 1:nch
            y(k,:) = scaleFac(k)*y(k,:);
         end
      end

      %   set up TSd
      TSd1 = setTSd(y,t0,dt,ch,sta,SysTF,dataZero);
      TSd1.sta = char(TSd1.sta);
      data0 = (TSd1.range(:,1)+TSd1.range(:,2))/2;
      rng = (TSd1.range(:,2)-TSd1.range(:,1))/2;
      rng = plotStart.rangeScale*rng;
      TSd1.range(:,1:2) = [ data0-rng data0+rng];
      TSd1.DataDir = dataDir;
      TSd1.arrayID = LDF{1}.Name;

      %   copy into cell arrays
      if iA == 1
         TSdata = {};
         TSd = {};
       end
      TSdata{iA} = y;
      clear y
      TSd{iA} = TSd1;
      clear TSd1
      LDF{1}.dt = dt;
      LOADED{iA}=LDF;
      set(gcbo,'UserData',LOADED);
      set(hfig,'Pointer','Arrow')
      % Enable buttons checks number of files in selected group,
      %   enables shift, proecessing buttons accordingly
       EnableButtons
      %  now set based on what is in selction window

      %   need sampling rate for shifting ... we are assuming that 
      %    the same sampling rate is used for all files loaded

   case 'Plot List'
      iA = get(gcbo,'Value');
      % Enable buttons checks number of files in selected group,
      %   enables shift, proecessing buttons accordingly
      EnableButtons

   case 'Plot'
      iA = get(findobj('Tag','Plot List','Parent',hfig),'Value');

      plotStart = get(gcbo,'UserData');
      %  reset nPw (# of points to plot initially on one page) now
      %   that dt has been read from data file  (NOTE: dt returned
      %   by readTSnims accounts for subsampling by a factor of decFac)
      dt = TSd{iA}.dt;
      nPw = fix(plotStart.tFac*plotStart.nTunits/dt);
      plotStart.nPw = nPw;

      rect = plotStart.rect;
      [hTSw,TSw1] = plotTSwin(TSd{iA},plotStart);
      TSfigInd(hTSw.Number) = iA;
      % set(hTSw,'Name',DataArrays{iA});
      i0 = TSw1.i0; i1 = TSw1.i1; di = TSw1.di;
      yFac = TSw1.yFac;
      ddt = TSd{iA}.dt/TSw1.tFac;
      t0 = (TSd{iA}.t0)*(1-TSw1.t0Eq0)/TSw1.tFac;
      I = i0:di:min(i1,size(TSdata{iA},2));
      [TSw1] = plotPage(ddt*I+t0,...
                yFac*TSdata{iA}(:,I),...
                TSw1,TSd{iA},1);
      TSw{iA} = TSw1;
      TSc = {};
      TSp= {};

       % LOADED = get(findobj('Tag','Load','Parent',hfig),'UserData');
       % FileNames = LOADED{iPlot};
   case 'Map'
      nSites = length(SelectedSites);
      if nSites==1 
          % do nothing 
      else
      Lats = zeros(nSites,1);Lons = Lats;
      for k = 1:nSites
        Lats(k) = SelectedSites{k}.lat;
        Lons(k) = SelectedSites{k}.lon;
      end
      minLat = min(Lats);
      minLon = min(Lons);
      maxLat = max(Lats);
      maxLon = max(Lons);
      midLat = (maxLat+minLat)/2;
      midLon = (maxLon+minLon)/2;
      X = 111*cos(pi*midLat/180)*(Lons-midLon);
      Y = 111*(Lats-midLat);
      xR = max(X)-min(X);
      yR = max(Y)-min(Y);
      asp = yR/xR;
      if(asp < 1)
         Rect = [100,100,400,400*asp];
         orient = 'LandScape';
      else
         Rect = [100,100,400/asp,400];
         orient = 'Portrait';
      end
      figure('Position',Rect,'PaperOrientation',orient,...
	'PaperPosition',Rect/100,'Name','Site Map');
      plot(X,Y,'r+','MarkerSize',14,'LineWidth',2)
      set(gca,'FontWeight','demi','FontSize',12)
      xlabel('km E')
      ylabel('km N')
      for k = 1:nSites
         ShortName = SelectedSites{k}.Name([3 end-1:end]);
         text(X(k)+5,Y(k)-5,ShortName,'FontWeight','demi');
      end
      end
      
   case 'Shift'
      iA = get(findobj('Tag','Plot List','Parent',hfig),'Value');
      LDF = LOADED{iA};
      set(hfig,'Pointer','watch')
      shiftPlotStart(TSdata{iA},LDF)
      set(hfig,'Pointer','arrow')
      
   case 'xmledit'
      SelectedDataFiles = get(findobj('Parent',hfig,'Tag','FileNames'),...
		'UserData');
      xmleditfig = checkHeaderGUI(SelectedDataFiles{1}.metaData,'xmlEdit');
      %disp(SelectedDataFiles{1}.metaData);
   case 'dnff'
      SelectedDataFiles = get(findobj('Parent',hfig,'Tag','FileNames'),...
		'UserData');
      currentDir = pwd;
      eval(['cd ''' dataHome ''''])
      binFile = [SelectedDataFiles{1}.Name '.bin'];
      dnffRunFile = ['tmp' filesep SelectedDataFiles{1}.Name 'FT.cfg'];
      logFile =  ['tmp' filesep SelectedDataFiles{1}.Name 'FT.log'];
      fid = fopen(dnffRunFile,'w');            
      fprintf(fid,'%s\n',binFile);
      fprintf(fid,'%s\n','n');
      fclose(fid);
      if exist(['FC' filesep SelectedDataFiles{1}.Name '.f5'],'file')
          system(['rm FC' filesep SelectedDataFiles{1}.Name '.f5']);
      end
      system([binDir 'dnff -B < ' dnffRunFile '>! ' logFile ' &'])
      eval(['cd ''' currentDir '''']);
   case 'RR'
      SelectedDataFiles = get(findobj('Parent',hfig,'Tag','FileNames'),...
		'UserData');
      optionsFile = get(gcbo,'UserData');
      ComputeTF(SelectedDataFiles,'RR',optionsFile);
   case 'SS'
      SelectedDataFiles = get(findobj('Parent',hfig,'Tag','FileNames'),...
		'UserData');
      optionsFile = get(gcbo,'UserData');
      ComputeTF(SelectedDataFiles,'SS',optionsFile);
   case 'MMT'
      SelectedDataFiles = get(findobj('Parent',hfig,'Tag','FileNames'),...
		'UserData');
      optionsFile = get(gcbo,'UserData');
      ComputeTF(SelectedDataFiles,'MMT',optionsFile);
    
   case 'MaxTimeString'
      hMax = gcbo;  
      hMin = findobj('Parent',hfig,'Tag','MinTimeString');
      tL2 = round(str2num(get(hMax,'String')));
      tL1 = round(str2num(get(hMin,'String')));
      Tdata = get(hMax,'UserData');
      %  check that tL2 is in range, adjust if not
      if tL2 > Tdata(2)
         tL2 = Tdata(2);
         set(hMax,'String',num2str(tL2));
      elseif tL2 < Tdata(1) + Tdata(3)
         tL2 = round(Tdata(1) + Tdata(3));
         set(hMax,'String',num2str(tL2));
      end
      %  set slider to adjusted tL2 value
      set(findobj('Parent',hfig,'Tag','MaxTimeSlider'),'Value',tL2);

      %  Then check to see if lower limit should be adjusted
      if tL1 > tL2 - Tdata(3)
         tL1 = round(tL2-Tdata(3));
         %  set text window, slider to adjusted tL1 value
         set(hMin,'String',num2str(tL1));
         set(findobj('Parent',hfig,'Tag','MinTimeSlider'),'Value',tL1);
      end
   case 'MinTimeString'
      hMin = gcbo;  
      hMax = findobj('Parent',hfig,'Tag','MaxTimeString');
      tL2 = str2num(get(hMax,'String'));
      tL1 = round(str2num(get(hMin,'String')));
      Tdata = get(hMin,'UserData');
      %  check that tL1 is in range, adjust if not
      if tL1 < Tdata(1)
         tL1 = Tdata(1);
         set(hMin,'String',num2str(tL1));
      elseif tL1 > Tdata(2) - Tdata(3)
         tL1 = Tdata(2) - round(Tdata(3));
         set(hMin,'String',num2str(tL1));
      end
      %  set slider to adjusted tL2 value
      set(findobj('Parent',hfig,'Tag','MinTimeSlider'),'Value',tL1);

      %  Then check to see if upper limit should be adjusted
      if tL2 < tL1 + Tdata(3)
         tL2 = tL1+round(Tdata(3));
         %  set text window, slider to adjusted tL2 value
         set(hMax,'String',num2str(tL2));
         set(findobj('Parent',hfig,'Tag','MaxTimeSlider'),'Value',tL2);
      end
   case 'MaxTimeSlider'
      hMax = gcbo;
      hMin = findobj('Parent',hfig,'Tag','MinTimeSlider');
      tL2 = round(get(hMax,'Value')); 
      tL1 = round(get(hMin,'Value')); 
      Tdata = get(hMax,'UserData');
      %  set edit window to new tL2 value
      set(findobj('Parent',hfig,'Tag','MaxTimeString'),...
			'String',num2str(tL2));
      %  Then check to see if lower limit should be adjusted
      if tL1 > tL2 - Tdata(3)
         tL1 = tL2-round(Tdata(3));
         %  set text window, slider to adjusted tL1 value
         set(hMin,'Value',tL1);
         set(findobj('Parent',hfig,'Tag','MinTimeString'),...
			'String',num2str(tL1));
      end
   case 'MinTimeSlider'
      hMin = gcbo;
      hMax = findobj('Parent',hfig,'Tag','MaxTimeSlider');
      tL2 = round(get(hMax,'Value')); 
      tL1 = round(get(hMin,'Value')); 
      Tdata = get(hMin,'UserData');
      %  set edit window to new tL1 value
      set(findobj('Parent',hfig,'Tag','MinTimeString'),...
			'String',num2str(tL1));
      %  Then check to see if upper limit should be adjusted
      if tL2 < tL1 + Tdata(3)
         tL2 = tL1+round(Tdata(3));
         %  set text window, slider to adjusted tL2 value
         set(hMax,'Value',tL2);
         set(findobj('Parent',hfig,'Tag','MaxTimeString'),...
			'String',num2str(tL2));
      end
   case 'Replot Data Times'
      tL1 = str2num(get(findobj('Parent',hfig,...
		'Tag','MinTimeString'),'String'));
      tL2 = str2num(get(findobj('Parent',hfig,...
		'Tag','MaxTimeString'),'String'));
      plotStart = get(findobj('Tag','Plot','Parent',hfig),'UserData');
      delete(hfig)
      clear TSd
      clear TSw
      clear TSdata
      dataTimes(DataFiles,NIMSinfo,Sites,plotStart,[tL1,tL2],HZ)
    case 'Sort Alph'
      tL1 = get(findobj('Parent',hfig,...
		'Tag','MinTimeSlider'),'Value');
      tL2 = get(findobj('Parent',hfig,...
		'Tag','MaxTimeSlider'),'Value');
      plotStart = get(findobj('Tag','Plot','Parent',hfig),'UserData');
      delete(hfig)
      clear TSd
      clear TSw
      clear TSdata
      dataQC(DataFiles,NIMSinfo,Sites,plotStart,[round(tL1),round(tL2)],HZ,1)
   case 'All Data Times'
      plotStart = get(findobj('Tag','Plot','Parent',hfig),'UserData');
      clear TSd
      clear TSw
      clear TSdata
      delete(hfig)
      dataTimes(DataFiles,NIMSinfo,Sites,plotStart,maxSites,HZ)   
    case 'Replot QC'
      tL1 = str2num(get(findobj('Parent',hfig,...
		'Tag','MinTimeString'),'String'))
      tL2 = str2num(get(findobj('Parent',hfig,...
		'Tag','MaxTimeString'),'String'))
      plotStart = get(findobj('Tag','Plot','Parent',hfig),'UserData');
      delete(hfig)
      clear TSd
      clear TSw
      clear TSdata
      dataQC(DataFiles,NIMSinfo,Sites,plotStart,[tL1,tL2],HZ)
   case 'QC'
      plotStart = get(findobj('Tag','Plot','Parent',hfig),'UserData');
      clear TSd
      clear TSw
      clear TSdata
      delete(hfig)
      dataQC(DataFiles,NIMSinfo,Sites,plotStart,maxSites,HZ)
   case 'Site Name'
      siteName = get(gcbo,'String');
      MTdir = [dataHome filesep 'MT'];
      cd(MTdir);
      zid = [ siteName '*.z*; ' lower(siteName) '*.z*; *.zmm' ];
      apresplt
      
   % Callbacks for the Summary function
   case 'Best TF'
      l = get(gcbo,'UserData');
      options = get(gcbo,'String');
      Sites{l}.bestTF = char(options(get(gcbo,'Value')));
      UD.Sites = Sites;
      set(hfig,'UserData',UD);
   case 'Good From'
      l = get(gcbo,'UserData');
      Sites{l}.period(1) = str2double(get(gcbo,'String'));
      UD.Sites = Sites;
      set(hfig,'UserData',UD);
   case 'Good To'
      l = get(gcbo,'UserData');
      Sites{l}.period(2) = str2double(get(gcbo,'String'));
      UD.Sites = Sites;
      set(hfig,'UserData',UD);
   case 'Data Quality'
      l = get(gcbo,'UserData');
      options = get(gcbo,'String');
      Sites{l}.quality = str2double(options(get(gcbo,'Value')));
      UD.Sites = Sites;
      set(hfig,'UserData',UD);
    case 'Load From File'
        [qcfile,qcdir] = uigetfile('*.qc');
        fid = fopen([qcdir filesep qcfile],'r');
        header = fgetl(fid);
        %qcinfo = textscan(fid,'%s %f %f %s %f %s %f %f %f %f %q\n',[11 inf]);
        qcinfo = textscan(fid,'%s %f %f %s %f %s %f %f %f %f %q\n',inf);
        fclose(fid);
        for l = 1:length(Sites)
            k = strmatch(Sites{l}.Name,qcinfo{1});
            if ~isempty(k)
                Sites{l}.lat = qcinfo{2}(k);
                Sites{l}.lon = qcinfo{3}(k);
                %Sites{l}.suffix = char(qcinfo{4}(k)); % skip suffix array
                %Sites{l}.days = qcinfo{5}(k);
                Sites{l}.bestTF = char(qcinfo{6}(k));
                Sites{l}.period(1) = qcinfo{7}(k);
                Sites{l}.period(2) = qcinfo{8}(k);
                Sites{l}.quality = qcinfo{9}(k);
                Sites{l}.flag = qcinfo{10}(k);
                Sites{l}.comments = char(qcinfo{11}(k));
            end
        end
        plotStart = get(findobj('Tag','Plot','Parent',hfig),'UserData');
        clear TSd
        clear TSw
        clear TSdata
        delete(hfig)
        dataQC(DataFiles,NIMSinfo,Sites,plotStart,maxSites)
    case 'Save To File'
        [qcfile,qcdir] = uiputfile('*.qc');
        fid = fopen([qcdir filesep qcfile],'wt');
        fprintf(fid,['Site ID\tLatitude\tLongitude\tRuns\tTime (days)\t' ...
            'Best TF\tGood From Period (secs)\tGood To Period (secs)\t' ...
            'Quality Rating\tWarning Flag\tComments\n'],header);
        for l = 1:length(Sites)
            id = Sites{l}.Name;
            lat = Sites{l}.lat;
            lon = Sites{l}.lon;
            suffix = Sites{l}.suffix;
            days = sum(Sites{l}.days);
            quality = Sites{l}.quality;
            bestTF = char(Sites{l}.bestTF);
            goodfrom = round(Sites{l}.period(1));
            goodto = round(Sites{l}.period(2));
            flag = Sites{l}.flag;
            comments = char(Sites{l}.comments);
            fprintf(fid,'%s\t%f\t%f\t%s\t%.1f\t%s\t%d\t%d\t%d',...
                id,lat,lon,suffix,days,bestTF,goodfrom,goodto,quality);
            if ~isnan(flag)
                fprintf(fid,'\t%d\t"%s"',flag,comments);
            end
            fprintf(fid,'\n');
        end
        fclose(fid);
    case 'Copy TFs'
        outdir = uigetdir(dataHome,'Copy TFs to this directory:');
        MTdir = [dataHome filesep 'MT'];
        for l = 1:length(Sites)
            zfile = [ MTdir filesep char(Sites{l}.bestTF) ];
            system(['cp ' zfile ' ' outdir]);
        end
end