From 28a0b09e3aa01f69cea46f13b38c9dbbb222610a Mon Sep 17 00:00:00 2001
From: AnnaCastelnovo <castelnovo.anna@gmail.com>
Date: Sun, 7 Jan 2018 17:50:04 +0100
Subject: [PATCH] Update swa_FindSWTravelling.m

---
 SW/swa_FindSWTravelling.m | 23 +++++++++++++++++++++++
 1 file changed, 23 insertions(+)

diff --git a/SW/swa_FindSWTravelling.m b/SW/swa_FindSWTravelling.m
index fba0d0e..1ea14b1 100644
--- a/SW/swa_FindSWTravelling.m
+++ b/SW/swa_FindSWTravelling.m
@@ -166,6 +166,29 @@
         SW(nSW).Travelling_Streams{end+1} = Streams{maxAngleId};
     end
 
+    %% ADDED BY ANNA
+    % find stream delays   % TODO (optimize the code)
+    y = []; x1 = []; x2 = [];
+    for x = 1 : Info.Parameters.Travelling_GS
+        for k = 1 : Info.Parameters.Travelling_GS
+            y(end + 1) = SW(nSW).Travelling_DelayMap(x, k);
+            x1(end + 1) = x;
+            x2(end + 1) = k;
+        end
+    end
+    tw = SW(nSW).Travelling_Streams{1}';
+    x1(isnan(y)) = [];
+    x2(isnan(y)) = [];
+    y (isnan(y)) = [];
+    Function = scatteredInterpolant(x1',x2',y');
+    SW(nSW).Stream_Travelling_Delay = Function(tw);
+    
+    % 3D projection & compute local speed
+    % TODO: scaleFactor: find relationship between EEG net model and real head size
+    scaleFactor = 2;
+    [SW(nSW).Stream_Coord, SW(nSW).Stream_Speed] = swa_Project3DSpeed(Info, SW(nSW), scaleFactor, 0);
+    
+
     % Update waitbar
     if flag_wait
         waitbar(nSW/length(SW),h,sprintf('Slow Wave %d of %d',nSW, length(SW)))