biomechzoo 0.4.9__tar.gz → 0.4.11__tar.gz

{biomechzoo-0.4.9/src/biomechzoo.egg-info → biomechzoo-0.4.11}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: biomechzoo
-Version: 0.4.9
+Version: 0.4.11
 Summary: Python implementation of the biomechZoo toolbox
 License-Expression: MIT
 Project-URL: Homepage, https://github.com/mcgillmotionlab/biomechzoo
@@ -37,6 +37,7 @@ See also http://www.github.com/mcgillmotionlab/biomechzoo or http://www.biomechz
 ### Installing a dev environment
 conda create -n biomechzoo-dev python=3.11
 conda activate biomechzoo-dev
+cd biomechzoo root folder
 pip install -e ".[dev]"
 
 ### import issues

{biomechzoo-0.4.9 → biomechzoo-0.4.11}/README.md

@@ -20,6 +20,7 @@ See also http://www.github.com/mcgillmotionlab/biomechzoo or http://www.biomechz
 ### Installing a dev environment
 conda create -n biomechzoo-dev python=3.11
 conda activate biomechzoo-dev
+cd biomechzoo root folder
 pip install -e ".[dev]"
 
 ### import issues

{biomechzoo-0.4.9 → biomechzoo-0.4.11}/pyproject.toml

@@ -1,6 +1,6 @@
 [project]
 name = "biomechzoo"
-version = "0.4.9"
+version = "0.4.11"
 description = "Python implementation of the biomechZoo toolbox"
 readme = "README.md"
 requires-python = ">=3.11,<3.12"   # max version for opencap-process tools is 3.11

{biomechzoo-0.4.9 → biomechzoo-0.4.11}/src/biomechzoo/biomech_ops/filter_data.py

@@ -13,10 +13,10 @@ def filter_data(data, ch, filt=None):
         The name(s) of the channel(s) to filter.
     filt : dict, optional
         Dictionary specifying filter parameters. Keys may include:
-        - 'type': 'butter' (default)
+        - 'ftype': 'butter' (default)
         - 'order': filter order (default: 4)
         - 'cutoff': cutoff frequency or tuple (Hz)
-        - 'btype': 'low', 'high', 'bandpass', 'bandstop' (default: 'low')
+        - 'btype': 'low', 'high', 'bandpass', 'bandstop' (default: 'lowpass')
 
     Returns
     -------
@@ -25,32 +25,34 @@ def filter_data(data, ch, filt=None):
     """
 
     if filt is None:
-        filt = {}
+        filt = {'ftype': 'butter',
+                'order': 4,
+                'cutoff': 10,
+                'btype': 'lowpass',
+                'filtfilt': True}
 
     if isinstance(ch, str):
         ch = [ch]
 
-    analog_channels = data['zoosystem']['Analog']['Channels']
-    if analog_channels:
-        analog_freq = data['zoosystem']['Analog']['Freq']
-    video_channels = data['zoosystem']['Video']['Channels']
-    if video_channels:
-        video_freq = data['zoosystem']['Video']['Freq']
-
+    # loop through all channels and filter
     for c in ch:
         if c not in data:
             raise KeyError('Channel {} not found in data'.format(c))
 
         if 'fs' not in filt:
+
+            video_channels = data['zoosystem']['Video']['Channels']
+            analog_channels = data['zoosystem']['Analog']['Channels']
+
             if c in analog_channels:
-                filt['fs'] = analog_freq
-            elif c in video_freq:
-                filt['fs'] = video_freq
+                filt['fs'] = data['zoosystem']['Analog']['Freq']
+            elif c in video_channels:
+                filt['fs'] =  data['zoosystem']['Video']['Freq']
             else:
-                raise ValueError('frequency not provided and cannot be inferred from zoosystem for channel'.format(c))
+                raise ValueError('Channel not analog or video')
 
         signal_raw = data[c]['line']
-        signal_filtered = filter_line(signal_raw, filt)
+        signal_filtered = filter_line(signal_raw=signal_raw, filt=filt)
         data[c]['line'] = signal_filtered
 
     return data

biomechzoo-0.4.11/src/biomechzoo/biomech_ops/filter_line.py

@@ -0,0 +1,85 @@
+import numpy as np
+import scipy.signal as sgl
+
+
+def filter_line(signal_raw, filt=None, fs=None):
+    """Filter an array using a Butterworth filter."""
+    #todo: verify that filter is working correctly
+    #todo add more filters
+    #todo: consider using kineticstoolkit
+
+    if filt is None:
+        filt = {'ftype': 'butter',
+                'order': 4,
+                'cutoff': 10,
+                'btype': 'lowpass',
+                'filtfilt': True}
+        if fs is None:
+            raise ValueError('fs is required if no filt is specified')
+
+    else:
+        if 'fs' not in filt:
+            raise ValueError('fs is a required key of filt')
+
+    # Normalize filter type strings
+    if filt['ftype'] == 'butterworth':
+        filt['ftype'] = 'butter'
+    if filt['btype'] is 'low':
+        filt['btype'] = 'lowpass'
+    if filt['btype'] is 'high':
+        filt['btype'] = 'highpass'
+
+    # Extract parameters
+    ftype = filt['ftype']
+    order = filt['order']
+    cutoff = filt['cutoff']
+    btype = filt['btype']
+    filtfilt = filt['filtfilt']
+    fs = filt['fs']
+
+    # prepare normalized cutoff(s)
+    nyq = 0.5 * fs
+    norm_cutoff = np.atleast_1d(np.array(cutoff) / nyq)
+
+    if ftype is 'butter':
+        signal_filtered = kt_butter(ts=signal_raw, fc=norm_cutoff, fs=fs, order=order, btype=btype, filtfilt=filtfilt)
+    else:
+        raise NotImplementedError(f"Filter type '{ftype}' not implemented.")
+
+    return signal_filtered
+
+
+def kt_butter(ts, fc, fs, order=2, btype='lowpass', filtfilt=True):
+    """
+    Apply a Butterworth filter to data.
+
+    Parameters
+    ----------
+    ts, ndarray, 1d.
+    fc, Cut-off frequency in Hz. This is a float for single-frequency filters
+        (lowpass, highpass), or a tuple of two floats (e.g., (10., 13.)
+        for two-frequency filters (bandpass, bandstop)).
+    order, Optional. Order of the filter. Default is 2.
+    btype, Optional. Can be either "lowpass", "highpass", "bandpass" or
+        "bandstop". Default is "lowpass".
+    filtfilt, Optional. If True, the filter is applied two times in reverse direction
+        to eliminate time lag. If False, the filter is applied only in forward
+        direction. Default is True.
+
+    Returns
+    -------
+    ts_f,  A copy of the input data which each data being filtered.
+    
+    Notes: 
+    - This code was adapted from kineticstoolkit Thanks @felxi
+    """
+
+    sos = sgl.butter(order, fc, btype, analog=False, output="sos", fs=fs)
+
+    # Filter
+    if filtfilt:
+        ts_f = sgl.sosfiltfilt(sos, ts, axis=0)
+    else:
+        ts_f = sgl.sosfilt(sos,ts, axis=0)
+
+    return ts_f

{biomechzoo-0.4.9 → biomechzoo-0.4.11}/src/biomechzoo/biomechzoo.py

@@ -8,7 +8,7 @@ from biomechzoo.utils.batchdisp import batchdisp
 from biomechzoo.utils.get_split_events import get_split_events
 from biomechzoo.utils.split_trial import split_trial
 from biomechzoo.conversion.c3d2zoo_data import c3d2zoo_data
-from biomechzoo.conversion.table2zoo import table2zoo_data
+from biomechzoo.conversion.table2zoo_data import table2zoo_data
 from biomechzoo.conversion.mvnx2zoo_data import mvnx2zoo_data
 from biomechzoo.processing.removechannel_data import removechannel_data
 from biomechzoo.processing.renamechannel_data import renamechannel_data
@@ -19,6 +19,7 @@ from biomechzoo.processing.renameevent_data import renameevent_data
 from biomechzoo.biomech_ops.normalize_data import normalize_data
 from biomechzoo.biomech_ops.phase_angle_data import phase_angle_data
 from biomechzoo.biomech_ops.continuous_relative_phase_data import continuous_relative_phase_data
+from biomechzoo.biomech_ops.filter_data import filter_data
 
 class BiomechZoo:
     def __init__(self, in_folder, inplace=False, subfolders=None, name_contains=None, verbose=0):
@@ -306,7 +307,7 @@ class BiomechZoo:
         # Update self.folder after  processing
         self._update_folder(out_folder, inplace, in_folder)
 
-    def addevent(self, ch, evt_type, evt_name, out_folder=None, inplace=None):
+    def addevent(self, ch, event_type, event_name, out_folder=None, inplace=None):
         """ adds events of type evt_type with name evt_name to channel ch """
         start_time = time.time()
         verbose = self.verbose
@@ -316,9 +317,9 @@ class BiomechZoo:
         fl = engine(in_folder, extension='.zoo', name_contains=self.name_contains, subfolders=self.subfolders)
         for f in fl:
             if verbose:
-                batchdisp('adding event {} to channel {} for {}'.format(evt_type, ch, f), level=2, verbose=verbose)
+                batchdisp('adding event {} to channel {} for {}'.format(event_type, ch, f), level=2, verbose=verbose)
             data = zload(f)
-            data = addevent_data(data, ch, evt_type, evt_name)
+            data = addevent_data(data, ch, event_type, event_name)
             zsave(f, data, inplace=inplace, out_folder=out_folder, root_folder=in_folder)
         method_name = inspect.currentframe().f_code.co_name
         batchdisp('{} process complete for {} file(s) in {:.2f} secs'.format(method_name, len(fl), time.time() - start_time), level=1, verbose=verbose)
@@ -346,28 +347,22 @@ class BiomechZoo:
         self._update_folder(out_folder, inplace, in_folder)
 
     def filter(self, ch, filt=None, out_folder=None, inplace=None):
-        raise NotImplementedError
-        # verbose = self.verbose
-        # in_folder = self.in_folder
-        # if inplace is None:
-        #     inplace = self.inplace
-        #
-        # # set filter type
-        # if filt is None:
-        #     filt = {'type': 'butterworth',
-        #             'order': 3,
-        #             'pass': 'lowpass'}
-        #
-        # fl = engine(in_folder)
-        # for f in fl:
-        #     batchdisp('filtering data in channels {} for {}'.format(ch, f), level=2, verbose=verbose)
-        #     data = zload(f)
-        #     data = filter_data(data, ch, filt)
-        #     zsave(f, data, inplace=inplace, root_folder=in_folder, out_folder=out_folder, verbose=verbose)
-        # method_name = inspect.currentframe().f_code.co_name
-        # batchdisp('{} computation complete for {} file(s)'.format(method_name, len(fl)), level=1, verbose=verbose)
-        #
-        # # Update self.folder after  processing
-        # self._update_folder(out_folder, inplace, in_folder)
-
+        """ filter data"""
+        start_time = time.time()
+        verbose = self.verbose
+        in_folder = self.in_folder
+        if inplace is None:
+            inplace = self.inplace
+        fl = engine(in_folder, name_contains=self.name_contains, subfolders=self.subfolders)
+        for f in fl:
+            if verbose:
+                batchdisp('filtering data for channel {} in {}'.format(ch, f), level=2, verbose=verbose)
+            data = zload(f)
+            data = filter_data(data, ch, filt)
+            zsave(f, data, inplace=inplace, out_folder=out_folder, root_folder=in_folder)
+        method_name = inspect.currentframe().f_code.co_name
+        batchdisp('{} process complete for {} file(s) in {:.2f} secs'.format(method_name, len(fl), time.time() - start_time),
+            level=1, verbose=verbose)
+        # Update self.folder after  processing
+        self._update_folder(out_folder, inplace, in_folder)
 

{biomechzoo-0.4.9 → biomechzoo-0.4.11}/src/biomechzoo/conversion/c3d2zoo_data.py

@@ -8,10 +8,14 @@ def c3d2zoo_data(c3d_obj):
     - data (dict): Zoo dictionary with 'line' and 'event' fields per channel.
     """
     data = {}
-
+    data['zoosystem'] = set_zoosystem()
+    video_freq = None
+    analog_freq = None
+    # extract "video" data
     if 'points' in c3d_obj['data']:
         points = c3d_obj['data']['points']  # shape: (4, n_markers, n_frames)
-        labels = c3d_obj['parameters']['POINT']['LABELS']['value']
+        labels = list(c3d_obj['parameters']['POINT']['LABELS']['value'])
+        video_freq = int(c3d_obj['parameters']['POINT']['RATE']['value'][0])
         for i, label in enumerate(labels):
             line_data = points[:3, i, :].T  # shape: (frames, 3)
             data[label] = {
@@ -19,9 +23,32 @@ def c3d2zoo_data(c3d_obj):
                 'event': {}  # empty for now
             }
 
-        params = c3d_obj['parameters']
-        video_freq = c3d_obj['parameters']['POINT']['RATE']['value'][0]
-        if 'EVENT' in params and 'TIMES' in params['EVENT']:
+        data['zoosystem']['Video']['Freq'] = video_freq
+        data['zoosystem']['Video']['Channels'] = labels
+
+    if 'analogs' in c3d_obj['data']:
+        analog_data = c3d_obj['data']['analogs']  # shape: (subframes, n_analog_channels, n_frames)
+        analog_labels = list(c3d_obj['parameters']['ANALOG']['LABELS']['value'])
+        analog_freq = int(c3d_obj['parameters']['ANALOG']['RATE']['value'][0])
+        # Flatten to 2D: (n_samples, n_channels)
+        # ezc3d stores analogs as subframes per frame, so we flatten across all
+        n_subframes, n_channels, n_frames = analog_data.shape
+        analog_data = analog_data.reshape(n_subframes * n_frames, n_channels)
+
+        for i, label in enumerate(analog_labels):
+            line_data = analog_data[:, i].reshape(-1, 1)  # shape: (samples, 1)
+            data[label] = {
+                'line': line_data,
+                'event': {},
+            }
+
+        data['zoosystem']['Analog']['Freq'] = analog_freq
+        data['zoosystem']['Analog']['Channels'] = analog_labels
+
+    # extract event information
+    params = c3d_obj['parameters']
+    if 'EVENT' in params and 'TIMES' in params['EVENT']:
+        if 'points' in c3d_obj['data']:
             times_array = params['EVENT']['TIMES']['value']
             frames = times_array[1]  # should be time depending on C3D file
 
@@ -61,8 +88,8 @@ def c3d2zoo_data(c3d_obj):
                     else:
                         data[labels[0]]['event'][key_name] = [frame-1, 0, 0]    # remove 1 to follow python
 
-    # todo add relevant meta data to zoosystem
-    data['zoosystem'] = set_zoosystem()
-    data['zoosystem']['Analog']['Freq'] = int(params['ANALOG']['RATE']['value'][0])
+    # add more zoosystem
+    if analog_freq is not None and video_freq is not None:
+        data['zoosystem']['AVR'] = analog_freq/video_freq
 
     return data

{biomechzoo-0.4.9 → biomechzoo-0.4.11}/src/biomechzoo/conversion/mvnx2zoo_data.py

@@ -10,6 +10,7 @@ def mvnx2zoo_data(fl):
 
     # create zoo data dict
     data = {}
+    data['zoosystem'] = set_zoosystem()
 
     # extract joint angle data (All JOINTS may not exist in a given dataset)
     for key, val in JOINTS.items():
@@ -58,7 +59,6 @@ def is_valid_for_zoo(val):
 
 def _get_meta_info(fl, mvnx_file, data):
     # todo: add more, see mvnx_file object
-    data['zoosystem'] = set_zoosystem(fl)
     data['zoosystem']['Video']['Freq'] = int(mvnx_file.frame_rate)
     data['zoosystem']['mvnx_version'] = mvnx_file.version
     data['zoosystem']['mvnx_configuration'] = mvnx_file.configuration

{biomechzoo-0.4.9 → biomechzoo-0.4.11}/src/biomechzoo/conversion/table2zoo_data.py

@@ -32,13 +32,15 @@ def table2zoo_data(csv_path, type='csv', skip_rows=0, freq=None):
         raise ValueError('type must be csv or parquet')
 
     # Use all columns
-    data = df.iloc[:, 0:]
+    df_data = df.iloc[:, 0:]
 
     # assemble zoo data
-    zoo_data = {}
-    for ch in data.columns:
-        zoo_data[ch] = {
-            'line': data[ch].values,
+    data = {}
+    data['zoosystem'] = set_zoosystem()
+
+    for ch in df_data.columns:
+        data[ch] = {
+            'line': df_data[ch].values,
             'event': []
         }
 
@@ -59,14 +61,12 @@ def table2zoo_data(csv_path, type='csv', skip_rows=0, freq=None):
             freq = compute_sampling_rate_from_time(time_data)
 
     # add metadata
-    # todo update zoosystem to match biomechzoo requirements
-    zoo_data['zoosystem'] = set_zoosystem(csv_path)
-    zoo_data['zoosystem']['Video']['Freq'] = freq
-    zoo_data['zoosystem']['Analog']['Freq'] = 'None'
+    data['zoosystem']['Video']['Freq'] = freq
+    data['zoosystem']['Analog']['Freq'] = 'None'
     if 'version' in metadata:
-        zoo_data['zoosystem']['collection_system_version'] = metadata['version']
+        data['zoosystem']['collection_system_version'] = metadata['version']
 
-    return zoo_data
+    return data
 
 
 def _parse_metadata(header_lines):

{biomechzoo-0.4.9 → biomechzoo-0.4.11}/src/biomechzoo/processing/addevent_data.py

@@ -37,6 +37,11 @@ def addevent_data(data, ch, ename, etype):
         elif etype == 'rom':
             eyd = float(np.max(yd) - np.min(yd))
             exd = 0  # dummy index (like MATLAB version)
+        elif etype == 'max_stance':
+            # special event for gait and running
+            exd = max_stance(yd)
+            eyd = float(yd[exd])
+            eyd = float(yd[exd])
         else:
             raise ValueError(f'Unknown event type: {etype}')
 
@@ -44,3 +49,8 @@ def addevent_data(data, ch, ename, etype):
         data[channel]['event'][ename] = [exd, eyd, 0]
 
     return data
+
+def max_stance(yd):
+    """ extracts max from first 40% of the gait cycle"""
+    raise NotImplementedError
+    return exd

{biomechzoo-0.4.9 → biomechzoo-0.4.11}/src/biomechzoo/processing/explodechannel_data.py

@@ -1,4 +1,5 @@
 import copy
+import numpy as np
 
 def explodechannel_data(data, channels=None):
     """ Explodes 3D channels (n x 3 arrays) into separate X, Y, Z channels.
@@ -6,13 +7,23 @@ def explodechannel_data(data, channels=None):
     Arguments:
         data (dict): Zoo data loaded from a file
         channels (list of str or None): Channels to explode.
-        If None, explode all channels with 'line' shaped (n x 3).
+            If None, explode all channels with 'line' shaped (n x 3).
 
     Returns:
         data_new (dict): Modified zoo dictionary with exploded channels.
     """
-
     data_new = copy.deepcopy(data)
+
+    # Ensure zoosystem channel lists are Python lists
+    for sys in ['Video', 'Analog']:
+        if sys in data_new.get('zoosystem', {}):
+            ch_list = data_new['zoosystem'][sys].get('Channels', [])
+            if isinstance(ch_list, np.ndarray):
+                ch_list = ch_list.tolist()
+            # strip whitespace
+            ch_list = [str(ch).strip() for ch in ch_list]
+            data_new['zoosystem'][sys]['Channels'] = ch_list
+
     # Find default channels if none provided
     if channels is None:
         channels = []
@@ -39,9 +50,20 @@ def explodechannel_data(data, channels=None):
             key = ch + axis
             data_new[key] = {
                 'line': line,
-                'event': data[ch]['event']}
+                'event': data_new[ch]['event']}
 
         # Remove original channel
         del data_new[ch]
 
+        # --- Update zoosystem lists ---
+        for sys in ['Video', 'Analog']:
+            if sys in data_new['zoosystem']:
+                ch_list = data_new['zoosystem'][sys]['Channels']
+                if ch in ch_list:
+                    # Remove original channel
+                    ch_list = [c for c in ch_list if c != ch]
+                    # Add exploded channels
+                    ch_list.extend([ch + '_x', ch + '_y', ch + '_z'])
+                    data_new['zoosystem'][sys]['Channels'] = ch_list
+
     return data_new

{biomechzoo-0.4.9 → biomechzoo-0.4.11}/src/biomechzoo/processing/removechannel_data.py

@@ -13,7 +13,6 @@ def removechannel_data(data, channels, mode='remove'):
     if mode not in ['remove', 'keep']:
         raise ValueError("mode must be 'remove' or 'keep'.")
 
-    zoosystem = data.get('zoosystem', {})
     all_channels = [ch for ch in data if ch != 'zoosystem']
 
     # Check for missing channels
@@ -28,9 +27,20 @@ def removechannel_data(data, channels, mode='remove'):
     else:
         raise ValueError("Mode must be 'remove' or 'keep'.")
 
-    # Build new zoo dictionary
-    data_new = {'zoosystem': zoosystem}
-    for ch in keep_channels:
-        data_new[ch] = data[ch]
+    # --- Compute channels to remove ---
+    remove_channels = [ch for ch in all_channels if ch not in keep_channels]
 
-    return data_new
+    if remove_channels:
+        print('Removing channels: {}'.format(remove_channels))
+    else:
+        print('No channels to remove')
+
+    # Remove from main data dict ---
+    for ch in remove_channels:
+        data.pop(ch, None)
+        if ch in data['zoosystem']['Video']['Channels']:
+            data['zoosystem']['Video']['Channels'] = [c for c in data['zoosystem']['Video']['Channels'] if c != ch]
+        if ch in data['zoosystem']['Analog']['Channels']:
+            data['zoosystem']['Analog']['Channels'] = [c for c in data['zoosystem']['Analog']['Channels'] if c != ch]
+
+    return data

{biomechzoo-0.4.9 → biomechzoo-0.4.11}/src/biomechzoo/utils/zload.py

@@ -1,5 +1,6 @@
 from scipy.io import loadmat
 import os
+import numpy as np
 
 
 def zload(filepath):
@@ -29,6 +30,16 @@ def zload(filepath):
     if 'data' in data:
         data = data['data']
 
+    # Convert Video and Analog channel arrays to Python lists
+    for sys in ['Video', 'Analog']:
+        if 'zoosystem' in data and sys in data['zoosystem']:
+            channels = data['zoosystem'][sys].get('Channels', [])
+            # Convert to list and strip spaces
+            if isinstance(channels, np.ndarray):
+                channels = channels.tolist()
+            channels = [str(ch).strip() for ch in channels]
+            data['zoosystem'][sys]['Channels'] = channels
+
     return data
 
 

{biomechzoo-0.4.9 → biomechzoo-0.4.11/src/biomechzoo.egg-info}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: biomechzoo
-Version: 0.4.9
+Version: 0.4.11
 Summary: Python implementation of the biomechZoo toolbox
 License-Expression: MIT
 Project-URL: Homepage, https://github.com/mcgillmotionlab/biomechzoo
@@ -37,6 +37,7 @@ See also http://www.github.com/mcgillmotionlab/biomechzoo or http://www.biomechz
 ### Installing a dev environment
 conda create -n biomechzoo-dev python=3.11
 conda activate biomechzoo-dev
+cd biomechzoo root folder
 pip install -e ".[dev]"
 
 ### import issues

biomechzoo-0.4.9/src/biomechzoo/biomech_ops/filter_line.py

@@ -1,88 +0,0 @@
-import numpy as np
-from scipy.signal import butter, filtfilt
-
-
-def filter_line(signal_raw, filt):
-    """ filter an array
-
-    Arguments
-    ----------
-    signal_raw : n, or n x 3 array signal to be filtered
-    filt : dict, optional
-        Dictionary specifying filter parameters. Keys may include:
-        - 'type': 'butter' (default)
-        - 'order': filter order (default: 4)
-        - 'cutoff': cutoff frequency or tuple (Hz)
-        - 'btype': 'low', 'high', 'bandpass', 'bandstop' (default: 'low')
-        - 'fs' frequency
-
-    Returns
-    -------
-    signal_filtered: filtered version of signal_raw"""
-    # todo allow for missing frequency to be obtained from zoosystem metadata
-    if filt is None:
-        filt = {}
-    if filt['type'] is 'butterworth':
-        filt['type'] = 'butter'
-    # Set default filter parameters
-    ftype = filt.get('type', 'butter')
-    order = filt.get('order', 4)
-    cutoff = filt.get('cutoff', None)
-    btype = filt.get('btype', 'low')
-    fs = filt.get('fs', None)
-
-    if ftype != 'butter':
-        raise NotImplementedError(f"Filter type '{ftype}' not implemented.")
-
-    if fs is None:
-        raise ValueError("Sampling frequency 'fs' must be specified in filt.")
-
-    if cutoff is None:
-        raise ValueError("Cutoff frequency 'cutoff' must be specified in filt.")
-
-    nyq = 0.5 * fs
-    norm_cutoff = np.array(cutoff) / nyq
-
-    b, a = butter(order, norm_cutoff, btype=btype, analog=False)
-
-    if signal_raw.ndim == 1:
-        signal_filtered = filtfilt(b, a, signal_raw)
-    else:
-        # Apply filter to each column if multivariate
-        signal_filtered = np.array([filtfilt(b, a, signal_raw[:, i]) for i in range(signal_raw.shape[1])]).T
-
-    return signal_filtered
-
-
-if __name__ == '__main__':
-    """ -------TESTING--------"""
-    import os
-    import matplotlib.pyplot as plt
-    from src.biomechzoo.utils.zload import zload
-    current_dir = os.path.dirname(os.path.abspath(__file__))
-    project_root = os.path.dirname(current_dir)
-    fl = os.path.join(project_root, 'data', 'other', 'HC030A05.zoo')
-    data = zload(fl)
-    data = data['data']
-    signal_raw = data['ForceFz1']['line']
-    filt = {'type': 'butterworth',
-            'order': 3,
-            'cutoff': 20,
-            'btype': 'low',
-            'fs': data['zoosystem']['Analog']['Freq']
-            }
-    signal_filtered = filter_line(signal_raw, filt)
-
-    # now plot
-    plt.figure(figsize=(10, 4))
-    plt.plot(signal_raw, label='Raw', alpha=0.6)
-    plt.plot(signal_filtered, label='Filtered', linewidth=2)
-    plt.xlabel('Frame')
-    plt.ylabel('Amplitude')
-    plt.title('Testing filter_line')
-    plt.legend()
-    plt.grid(True)
-    plt.tight_layout()
-    plt.show()
-
-