From 992b84ea72e1d20b395ab8d3d50abbd494c1a749 Mon Sep 17 00:00:00 2001
From: Théo de la Hogue
Date: Tue, 10 Oct 2023 13:51:11 +0200
Subject: Fixing KCoefficient formula.

---
 src/argaze/GazeAnalysis/KCoefficient.py | 29 ++++++++++++++++++++---------
 1 file changed, 20 insertions(+), 9 deletions(-)

diff --git a/src/argaze/GazeAnalysis/KCoefficient.py b/src/argaze/GazeAnalysis/KCoefficient.py
index 80fe1fd..c50bc3a 100644
--- a/src/argaze/GazeAnalysis/KCoefficient.py
+++ b/src/argaze/GazeAnalysis/KCoefficient.py
@@ -52,19 +52,24 @@ class ScanPathAnalyzer(GazeFeatures.ScanPathAnalyzer):
         duration_std = numpy.std(durations)
         amplitude_std = numpy.std(amplitudes)
 
-        Ks = []
-        for scan_step in scan_path:
+        if duration_std > 0. and amplitude_std > 0.:
+
+            Ks = []
+            for scan_step in scan_path:
+
+                Ks.append((abs(scan_step.duration - duration_mean) / duration_std) - (abs(scan_step.last_saccade.amplitude - amplitude_mean) / amplitude_std))
+
+            self.__K = numpy.array(Ks).mean()
 
-            Ks.append(((scan_step.duration - duration_mean) / duration_std) - ((scan_step.last_saccade.amplitude - amplitude_mean) / amplitude_std))
+        else:
 
-        self.__K = numpy.array(Ks).mean()
+            self.__K = 0.
 
     @property
     def K(self) -> float:
         """K coefficient."""
 
         return self.__K
-    
 
 @dataclass
 class AOIScanPathAnalyzer(GazeFeatures.AOIScanPathAnalyzer):
@@ -104,12 +109,18 @@ class AOIScanPathAnalyzer(GazeFeatures.AOIScanPathAnalyzer):
         duration_std = numpy.std(durations)
         amplitude_std = numpy.std(amplitudes)
 
-        Ks = []
-        for aoi_scan_step in aoi_scan_path:
+        if duration_std > 0. and amplitude_std > 0.:
+
+            Ks = []
+            for aoi_scan_step in aoi_scan_path:
+
+                Ks.append((abs(aoi_scan_step.duration - duration_mean) / duration_std) - (abs(aoi_scan_step.last_saccade.amplitude - amplitude_mean) / amplitude_std))
+
+            self.__K = numpy.array(Ks).mean()
 
-            Ks.append(((aoi_scan_step.duration - duration_mean) / duration_std) - ((aoi_scan_step.last_saccade.amplitude - amplitude_mean) / amplitude_std))
+        else:
 
-        self.__K = numpy.array(Ks).mean()
+            self.__K = 0.
 
     @property
     def K(self) -> float:
-- 
cgit v1.1