From 6dd4a4f5de57c1af9f4c680dc70bb23dc6b730df Mon Sep 17 00:00:00 2001
From: Jan Grewe <jan.grewe@g-node.org>
Date: Wed, 14 Sep 2022 09:53:17 +0200
Subject: [PATCH] [io] adaptor for nixtrack data

---
 etrack/io/nixtrack_data.py | 63 ++++++++++++++++++++++++++++++++++++++
 1 file changed, 63 insertions(+)
 create mode 100644 etrack/io/nixtrack_data.py

diff --git a/etrack/io/nixtrack_data.py b/etrack/io/nixtrack_data.py
new file mode 100644
index 0000000..9eeab79
--- /dev/null
+++ b/etrack/io/nixtrack_data.py
@@ -0,0 +1,63 @@
+import os
+import numpy as np
+import pandas as pd 
+import numbers as nb
+import nixtrack as nt
+
+from ..tracking_data import TrackingData
+from IPython import embed
+
+
+class NixtrackData(object):
+    
+    def __init__(self, results_file, crop=(0, 0)) -> None:
+        """
+        If the video data was cropped before tracking and the tracked positions are with respect to the cropped images, we may want to correct for this to convert the data back to absolute positions in the video frame.
+
+        Parameters
+        ----------
+        crop : 2-tuple
+            tuple of (xoffset, yoffset)
+
+        Raises
+        ------
+        ValueError if crop value is not a 2-tuple
+        """
+        if not os.path.exists(results_file):
+            raise ValueError("File %s does not exist!" % results_file)
+        if not isinstance(crop, tuple) or len(crop) < 2:
+            raise ValueError("Cropping info must be a 2-tuple of (x, y)")
+        self._file_name = results_file
+        self._crop = crop
+        self._dataset = nt.Dataset(self._file_name, nt.util.FileMode.ReadOnly)
+        if not self._dataset.is_open:
+            raise ValueError(f"An error occurred opening file {self._file_name}! File is not open!")
+
+    @property
+    def filename(self):
+        return self._file_name
+
+    @property
+    def bodyparts(self):
+        return self._dataset.nodes
+    
+    def _correct_cropping(self, orgx, orgy):
+        x = orgx + self._crop[0]
+        y = orgy + self._crop[1]
+        return x, y
+    
+    @property
+    def tracks(self):
+        return self._dataset.tracks
+
+    def track(self, track=None, bodypart=0):
+        if isinstance(bodypart, nb.Number):
+            bp = self.bodyparts[bodypart]
+        elif isinstance(bodypart, str) and bodypart in self.bodyparts:
+            bp = bodypart
+        else:
+            raise ValueError(f"Body part {bodypart} is not a tracked node!")
+        positions, time, iscore, nscore = self._dataset.positions(node=bp, axis_type=nt.AxisType.Time)
+        embed()
+
+        return TrackingData(positions[:, 0], positions[:, 1], time, nscore, bp, fps=self._dataset.fps)
\ No newline at end of file