From 3f9f27f31d8328a28133c38603b460deabb55f81 Mon Sep 17 00:00:00 2001
From: Théo de la Hogue
Date: Thu, 24 Nov 2022 08:42:03 +0100
Subject: Testing dispersion based fixation and saccade identification.

---
 .../DispersionBasedGazeMovementIdentifier.py       | 117 +++++++++++++++++++++
 1 file changed, 117 insertions(+)
 create mode 100644 src/argaze.test/GazeAnalysis/DispersionBasedGazeMovementIdentifier.py

(limited to 'src/argaze.test/GazeAnalysis/DispersionBasedGazeMovementIdentifier.py')

diff --git a/src/argaze.test/GazeAnalysis/DispersionBasedGazeMovementIdentifier.py b/src/argaze.test/GazeAnalysis/DispersionBasedGazeMovementIdentifier.py
new file mode 100644
index 0000000..6d84333
--- /dev/null
+++ b/src/argaze.test/GazeAnalysis/DispersionBasedGazeMovementIdentifier.py
@@ -0,0 +1,117 @@
+#!/usr/bin/env python
+
+import unittest
+import random
+import time
+
+from argaze import GazeFeatures
+from argaze.GazeAnalysis import DispersionBasedGazeMovementIdentifier
+
+import numpy
+
+def build_gaze_fixation(size: int, center: tuple, dispersion: float, start_time: float, min_time: float, max_time: float):
+    """ Generate N TimeStampedGazePsoitions dispersed around a center point for testing purpose.
+        Timestamps are current time after random sleep (second).
+        GazePositions are random values.
+    """
+    ts_gaze_positions = GazeFeatures.TimeStampedGazePositions()
+
+    for i in range(0, size):
+
+        # Edit gaze position
+        random_x = center[0] + dispersion * (random.random() - 0.5)
+        random_y = center[1] + dispersion * (random.random() - 0.5)
+        random_gaze_position = GazeFeatures.GazePosition((random_x, random_y))
+
+        # Store gaze position
+        ts = time.time() - start_time
+        ts_gaze_positions[ts] = random_gaze_position
+
+        #print(ts, random_x, random_y)
+
+        # Sleep a random time
+        sleep_time = random.random() * (max_time - min_time) + min_time
+        time.sleep(sleep_time)
+
+    return ts_gaze_positions
+
+class TestDispersionBasedGazeMovementIdentifierClass(unittest.TestCase):
+    """Test DispersionBasedGazeMovementIdentifier class."""
+
+    def test_fixation_identification(self):
+        """Test DispersionBasedGazeMovementIdentifier fixation identification."""
+
+        size = 10
+        center = (0, 0)
+        dispersion = 10
+        start_time = time.time()
+        min_time = 0.01
+        max_time = 0.1
+        
+        ts_gaze_positions = build_gaze_fixation(size, center, dispersion, start_time, min_time, max_time)
+        gaze_movement_identifier = DispersionBasedGazeMovementIdentifier.GazeMovementIdentifier(dispersion_threshold=dispersion, duration_threshold=min_time*2)
+        ts_fixations, ts_saccades, ts_status = gaze_movement_identifier.identify(ts_gaze_positions)
+
+        # Check result size
+        self.assertEqual(len(ts_fixations), 1)
+        self.assertEqual(len(ts_saccades), 0)
+        self.assertEqual(len(ts_status), size)
+
+        # Check fixation
+        ts, fixation = ts_fixations.pop_first()
+
+        self.assertEqual(len(fixation.positions.keys()), size)
+        self.assertLessEqual(fixation.dispersion, dispersion)
+        self.assertGreaterEqual(fixation.duration, size * min_time)
+        self.assertLessEqual(fixation.duration, size * max_time)
+
+    def test_fixation_and_saccade_identification(self):
+        """Test DispersionBasedGazeMovementIdentifier fixation and saccade identification."""
+
+        size = 10
+        center_A = (0, 0)
+        center_B = (50, 50)
+        dispersion = 10
+        start_time = time.time()
+        min_time = 0.01
+        max_time = 0.1
+        
+        ts_gaze_positions_A = build_gaze_fixation(size, center_A, dispersion, start_time, min_time, max_time)
+        ts_gaze_positions_B = build_gaze_fixation(size, center_B, dispersion, start_time, min_time, max_time)
+
+        ts_gaze_positions = ts_gaze_positions_A.append(ts_gaze_positions_B)
+
+        gaze_movement_identifier = DispersionBasedGazeMovementIdentifier.GazeMovementIdentifier(dispersion_threshold=dispersion, duration_threshold=min_time*2)
+        ts_fixations, ts_saccades, ts_status = gaze_movement_identifier.identify(ts_gaze_positions)
+
+        # Check result size
+        self.assertEqual(len(ts_fixations), 2)
+        self.assertEqual(len(ts_saccades), 1)
+        self.assertEqual(len(ts_status), size*2)
+
+        # Check first fixation
+        ts, fixation = ts_fixations.pop_first()
+
+        self.assertEqual(len(fixation.positions.keys()), size)
+        self.assertLessEqual(fixation.dispersion, dispersion)
+        self.assertGreaterEqual(fixation.duration, size * min_time)
+        self.assertLessEqual(fixation.duration, size * max_time)
+
+        # Check first saccade
+        ts, saccade = ts_saccades.pop_first()
+
+        self.assertEqual(len(saccade.positions.keys()), 2)
+        self.assertGreaterEqual(saccade.duration, min_time)
+        self.assertLessEqual(saccade.duration, max_time)
+
+        # Check second fixation
+        ts, fixation = ts_fixations.pop_first()
+
+        self.assertEqual(len(fixation.positions.keys()), size)
+        self.assertLessEqual(fixation.dispersion, dispersion)
+        self.assertGreaterEqual(fixation.duration, size * min_time)
+        self.assertLessEqual(fixation.duration, size * max_time)
+
+if __name__ == '__main__':
+
+    unittest.main()
\ No newline at end of file
-- 
cgit v1.1