added epoch pb & partial sacdet

test_experiment.m

@@ -78,4 +78,130 @@ for i = 1:3:10
     
     i
 end
+%%
+
+fs          = 120;
+duration    = 3;
+samples     = fs*duration - 1;
+%%
+
+for iB = 1:length(Data.timestamps)
+   
+   % Empty traces
+   E.AzChairEpoched   	= [];
+   E.ElChairEpoched    = [];
+   E.AzGazeEpoched     = [];
+   E.ElGazeEpoched     = [];
+   E.AzEyeEpoched      = [];
+   E.ElEyeEpoched      = [];
+   E.AzHeadEpoched     = [];
+   E.ElHeadEpoched     = [];
+   
+   % Interpolate Gaze
+   lsl_opti  	= Data.timestamps(iB).optitrack;
+   Data.timestamps(iB).epoch_interp       = 0:1/120:lsl_opti(end)-lsl_opti(1);
+   Data.position(iB).gaze_interp(:,1)     = interp1(lsl_opti-lsl_opti(iB), Data.position(iB).gaze(:,1), Data.timestamps(iB).epoch_interp,'pchip')';
+   Data.position(iB).gaze_interp(:,2)  	= interp1(lsl_opti-lsl_opti(iB), Data.position(iB).gaze(:,2), Data.timestamps(iB).epoch_interp,'pchip')';
+   
+   % Interpolate Chair
+   CUT_OFF = 203;
+   Data.position(iB).chair_interp(:,1)    = interp1(Data.timestamps(iB).chair, Data.position(iB).chair, Data.timestamps(iB).epoch_interp,'pchip')';
+   Data.position(iB).chair_interp(Data.timestamps(iB).epoch_interp>CUT_OFF,1)  = zeros(1,sum(Data.timestamps(iB).epoch_interp>CUT_OFF));              % Correct for extrapolation == 0;
+   Data.position(iB).chair_interp(:,2)    = zeros(size(Data.position(iB).chair_interp(:,1)));
+   
+   % Interpolate Eye
+   Data.position(iB).eye_interp(:,1)      = interp1(lsl_opti-lsl_opti(iB), Data.position(iB).pupillabs(:,1), Data.timestamps(iB).epoch_interp,'pchip')';
+   Data.position(iB).eye_interp(:,2)      = interp1(lsl_opti-lsl_opti(iB), Data.position(iB).pupillabs(:,2), Data.timestamps(iB).epoch_interp,'pchip')';
+   
+   % Interpolate Head
+   Data.position(iB).head_interp(:,1)      = interp1(lsl_opti-lsl_opti(iB), Data.position(iB).optitrack(:,1), Data.timestamps(iB).epoch_interp,'pchip')';
+   Data.position(iB).head_interp(:,2)      = interp1(lsl_opti-lsl_opti(iB), Data.position(iB).optitrack(:,2), Data.timestamps(iB).epoch_interp,'pchip')';
+   
+   
+   % Select stimuli indices
+   nstim       = length(Data.stimuli(iB).azimuth);
+   ntriggers   = length(Data.timestamps(iB).stimuli);
+   ind         = 1;
+   ext         = 2;
+
+   for iS = 1:nstim
+      % epoch for stimuli
+      start             = Data.timestamps(iB).stimuli(ind) - lsl_opti(1);
+      [~,idx]           = min(abs(Data.timestamps(iB).epoch_interp-start));
+      
+      % Gaze
+      E.AzGazeEpoched  	= [E.AzGazeEpoched, Data.position(iB).gaze_interp(idx:idx+samples,1)'];
+      E.ElGazeEpoched  	= [E.ElGazeEpoched, Data.position(iB).gaze_interp(idx:idx+samples,2)'];
+      
+      % Chair
+      E.AzChairEpoched  = [E.AzChairEpoched, Data.position(iB).chair_interp(idx:idx+samples,1)'];
+      E.ElChairEpoched 	= [E.ElChairEpoched, Data.position(iB).chair_interp(idx:idx+samples,2)'];
+      
+      % Eye
+      E.AzEyeEpoched   	= [E.AzEyeEpoched, Data.position(iB).eye_interp(idx:idx+samples,1)'];
+      E.ElEyeEpoched   	= [E.ElEyeEpoched, Data.position(iB).eye_interp(idx:idx+samples,2)'];
+
+      % Head
+      E.AzHeadEpoched   = [E.AzHeadEpoched, Data.position(iB).head_interp(idx:idx+samples,1)'];
+      E.ElHeadEpoched   = [E.ElHeadEpoched, Data.position(iB).head_interp(idx:idx+samples,2)'];
+      
+      ind = ind + ext;
+   end
+   Data.epoch(iB) = E;
+end
+
+
+fn = strrep('preprocessed_data_JJH-0001-21-04-15-0001.mat','preprocessed','epoched');
+save(fn,'Data');
+%%
+if ~exist('sacdet','dir')
+       mkdir('sacdet')
+end
 
+cfn   = pb_clean('cd','C:\Users\marc_\OneDrive\Documents\Studie\STAGE\MATLAB Code\Data\Recordings\test experiment');
+path  = pb_getdir('cdir',cd);
+cd(path)
+
+%  Load data
+l     = dir('epoched*.mat');
+
+if isempty(l); return; end
+fn    = l(1).name;
+load(fn);
+
+%  Switch to sacdet folder
+path  = [path filesep 'sacdet'];
+cd(path)
+%%	Run sacdet over epoched data
+%  Read converted data file, preprocess all blocks, and store data
+
+l     = dir('*.sac');
+if ~isempty(l); return; end
+
+
+%% Make hv & csv and run sacdet
+
+
+Dlen        = length(Data.epoch);
+fs          = 120;
+duration    = 3;
+samples     = duration * fs;
+
+for iB = 1:Dlen
+
+   % Saving calibrated data
+   fname                   = fcheckext(['sacdet_' fn(14:end-5) '_block_' num2str(iB) '_azel'] ,'.hv');
+   fid                     = fopen([path filesep fname],'w','l');
+   AZEL                    = [Data.epoch(iB).AzGazeEpoched; Data.epoch(iB).ElGazeEpoched];
+   
+   fwrite(fid,AZEL,'float');
+   fclose(fid);
+   
+   fn_csv = [path filesep fname];
+   VC2csv(fn_csv,fs,samples,1:length(Data.epoch(iB).AzGazeEpoched)/samples);
+   
+   if ~isfile(fn_csv(1:end-3)) == 1
+       pa_sacdet;
+       pause;
+   end
+end
\ No newline at end of file