c3d-parser 1.5.4__tar.gz → 1.6.1__tar.gz

{c3d_parser-1.5.4/src/c3d_parser.egg-info → c3d_parser-1.6.1}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: c3d-parser
-Version: 1.5.4
+Version: 1.6.1
 Summary: C3D parser for gait data harmonisation.
 Author-email: Timothy Salemink <tim.nicolas@outlook.com>, Sally Jack <sallyjaack@gmail.com>
 License-Expression: Apache-2.0
@@ -12,7 +12,7 @@ Requires-Dist: numpy==1.26.4
 Requires-Dist: pandas
 Requires-Dist: c3d==0.5.2
 Requires-Dist: trc-data-reader>=0.4.1
-Requires-Dist: opensim-model-creator==0.8.3
+Requires-Dist: opensim-model-creator==0.10.1
 Requires-Dist: PySide6
 Requires-Dist: matplotlib
 Requires-Dist: mplcursors

{c3d_parser-1.5.4 → c3d_parser-1.6.1}/requirements.txt

@@ -3,7 +3,7 @@ numpy==1.26.4
 pandas
 c3d==0.5.2
 trc-data-reader>=0.4.1
-opensim-model-creator==0.8.3
+opensim-model-creator==0.10.1
 PySide6
 matplotlib
 mplcursors

c3d_parser-1.6.1/src/c3d_parser/__init__.py

@@ -0,0 +1,2 @@
+
+__version__ = "1.6.1"

{c3d_parser-1.5.4 → c3d_parser-1.6.1}/src/c3d_parser/core/c3d_parser.py

@@ -16,7 +16,6 @@ from decimal import Decimal, ROUND_CEILING, ROUND_FLOOR
 
 from trc import TRCData
 from opensim_model_creator.Create_Model import create_model
-from ll_visualiser.utils import get_fit_metrics, load_landmarks, define_measurements
 
 from c3d_parser.core.c3d_patch import c3d
 from c3d_parser.core.utils import clear_directory
@@ -38,11 +37,15 @@ torso_markers = ["C7", "T2", "T10", "MAN"]
 required_markers = [{"LASI", "RASI"}, {"LKNE", "RKNE"}, {"LANK", "RANK"}, {"LMED", "RMED"}, {"LHEE", "RHEE"},
                     {"LTOE", "RTOE"}, ({"LPSI", "RPSI"}, {"SACR"}), ({"LKNEM", "RKNEM"}, {"LKAX", "RKAX"})]
 
+required_measurements = ['Sex', 'Age', 'Height', 'Weight', 'ASIS Width', 'Left Knee Width', 'Right Knee Width',
+                         'Left Ankle Width', 'Right Ankle Width']
+
 
 def parse_session(static_trial, dynamic_trials, input_directory, output_directory, lab, marker_diameter, static_data,
                   left_foot_flat, right_foot_flat, toe_marker_proximal, optimise_knee_axis, filter_trc, filter_grf,
                   ik_task_set, running_gait, progress_tracker):
 
+    validate_static_data(static_data)
     clear_directory(output_directory)
 
     logger.info(f"Processing session {os.path.normpath(input_directory)}.")
@@ -89,7 +92,8 @@ def parse_session(static_trial, dynamic_trials, input_directory, output_director
                                    output_directory, left_foot_flat, right_foot_flat, toe_marker_proximal,
                                    optimise_knee_axis, progress_tracker)
 
-    write_c3d_parser_history(input_directory, output_directory, static_trial, deidentified_file_names, static_data)
+    write_c3d_parser_history(input_directory, output_directory, lab, static_trial, deidentified_file_names, static_data,
+                             left_foot_flat, right_foot_flat, toe_marker_proximal)
 
     if not trc_file_paths:
         raise CancelException("No dynamic trials found.")
@@ -127,16 +131,26 @@ def approximate_anthropometrics(c3d_file, lab, marker_diameter):
 
 def calculate_anthropometrics(frame_data, marker_diameter):
 
-    def distance(frame, landmark_1, landmark_2):
-        return round(np.linalg.norm(np.array(frame[landmark_1]) - np.array(frame[landmark_2])), 1)
+    def distance(frame, landmark_1, landmark_2, marker_diameter=0):
+        missing = [landmark for landmark in (landmark_1, landmark_2) if landmark not in frame]
+        if missing:
+            logger.warn(f"Unable to approximate distance between {landmark_1} and {landmark_2}. "
+                        f"Missing marker(s): {', '.join(missing)}.")
+            return 0
+
+        point_1 = np.array(frame[landmark_1])
+        point_2 = np.array(frame[landmark_2])
+        length = round(np.linalg.norm(point_1 - point_2), 1)
+
+        return length - marker_diameter
 
     mid_frame = frame_data.iloc[len(frame_data) // 2].copy()
 
     anthropometrics = {
-        'left_knee_width': distance(mid_frame, 'LKNE', 'LKNEM') - marker_diameter,
-        'right_knee_width': distance(mid_frame, 'RKNE', 'RKNEM') - marker_diameter,
-        'left_ankle_width': distance(mid_frame, 'LANK', 'LMED')  - marker_diameter,
-        'right_ankle_width': distance(mid_frame, 'RANK', 'RMED')  - marker_diameter,
+        'left_knee_width': distance(mid_frame, 'LKNE', 'LKNEM', marker_diameter),
+        'right_knee_width': distance(mid_frame, 'RKNE', 'RKNEM', marker_diameter),
+        'left_ankle_width': distance(mid_frame, 'LANK', 'LMED', marker_diameter),
+        'right_ankle_width': distance(mid_frame, 'RANK', 'RMED', marker_diameter),
         'left_leg_length': distance(mid_frame, 'LASI', 'LMED'),
         'right_leg_length': distance(mid_frame, 'RASI', 'RMED'),
         'inter_asis_distance': distance(mid_frame, 'LASI', 'RASI'),
@@ -238,15 +252,18 @@ def parse_dynamic_trial(c3d_file, lab, output_directory, trial_index, marker_dat
     return analog_data, events, s_t_data, trc_file_path, grf_file_path
 
 
-def write_c3d_parser_history(input_directory, output_directory, static_trial, deidentified_file_names, static_data):
+def write_c3d_parser_history(input_directory, output_directory, lab, static_trial, deidentified_file_names, static_data,
+                             left_foot_flat, right_foot_flat, toe_marker_proximal):
     log_path = os.path.join(input_directory, "c3d_parser_history.log")
     timestamp = datetime.now().strftime("%Y-%m-%d %H:%M")
 
     try:
         with open(log_path, "w", encoding="utf-8") as f:
             f.write("C3D-Parser History\n")
-            f.write(f"Last Run: {timestamp}\n\n")
+            f.write(f"Last Run: {timestamp}\n")
+            f.write(f"C3D-Parser Version: {VERSION}\n\n")
 
+            f.write(f"Lab: {lab}\n")
             f.write(f"Static trial: {static_trial}\n\n")
 
             if deidentified_file_names:
@@ -273,6 +290,12 @@ def write_c3d_parser_history(input_directory, output_directory, static_trial, de
                     unit = f" {units[key]}" if key in units else ""
                     f.write(f"{label}: {value}{unit}\n")
 
+            f.write("\n")
+            f.write("Foot Settings:\n")
+            f.write(f"Toe Marker: {'Proximal' if toe_marker_proximal else '2nd/3rd Met'}\n")
+            f.write(f"Left Foot Flat: {left_foot_flat}\n")
+            f.write(f"Right Foot Flat: {right_foot_flat}\n")
+
             metrics, measurements = define_metrics_and_measurements(output_directory)
             if metrics:
                 f.write("\n")
@@ -290,6 +313,8 @@ def write_c3d_parser_history(input_directory, output_directory, static_trial, de
 
 
 def define_metrics_and_measurements(output_directory):
+    from ll_visualiser.utils import get_fit_metrics, load_landmarks, define_measurements
+
     model_directory = os.path.join(output_directory, "Models", "Meshes")
     left_predicted_landmark_file = os.path.join(model_directory, 'predicted_lms_left.txt')
     right_predicted_landmark_file = os.path.join(model_directory, 'predicted_lms_right.txt')
@@ -1700,3 +1725,23 @@ def get_subject_info(static_data):
         'right_malleolar_width'
     ])
     return info
+
+
+def height_to_mm(height):
+    if 0 < height < 300:
+        logger.info(f"Subject height ({height}) assumed to be in cm. Converting to mm.")
+        return height * 10
+
+    return height
+
+
+def validate_static_data(static_data):
+    missing = [key for key in required_measurements if not static_data.get(key)]
+    if missing:
+        raise ParserError("Subject measurements missing:\n- " + "\n- ".join(missing))
+
+    static_data['Height'] = height_to_mm(static_data['Height'])
+
+    # Cap age at 18.
+    # Current shape model is paediatric only, ages over 18 can give strange results.
+    static_data['Age'] = min(static_data['Age'], 18)

{c3d_parser-1.5.4 → c3d_parser-1.6.1}/src/c3d_parser/view/main_window.py

@@ -17,7 +17,7 @@ from ll_visualiser.visualiser import visualise_model
 
 from c3d_parser.core.c3d_parser import (parse_session, extract_static_data, extract_marker_names, is_dynamic,
     CancelException, write_normalised_grfs, write_normalised_kinematics, write_normalised_kinetics,
-    write_spatiotemporal_data, approximate_anthropometrics)
+    write_spatiotemporal_data, approximate_anthropometrics, height_to_mm)
 from c3d_parser.settings.general import (APPLICATION_NAME, VERSION, DEFAULT_STYLE_SHEET, INVALID_STYLE_SHEET,
                                          get_marker_maps_dir)
 from c3d_parser.view.ui.ui_main_window import Ui_MainWindow
@@ -124,7 +124,6 @@ class MainWindow(QMainWindow):
         self._kinetic_curves.signals.cycles_included.connect(self._update_cycle_count)
 
     def _setup_progress_bar(self):
-        self._progress_text = ""
         self._progress_value = 0
         self._ui.progressBar.setVisible(False)
 
@@ -160,7 +159,6 @@ class MainWindow(QMainWindow):
         self._plotter.clear_actors()
 
     def _clear_progress_bar(self):
-        self._progress_text = ""
         self._ui.labelProgress.setText("")
         self._progress_value = 0
         self._ui.progressBar.setValue(0)
@@ -460,33 +458,32 @@ class MainWindow(QMainWindow):
                 self._set_subject_info(static_trial)
                 self._static_trial = static_trial
 
-            # Fill in any anthropometric gaps.
-            if self._approximate_anthropometrics:
-                input_directory = self._ui.lineEditInputDirectory.text()
-                lab = self._ui.comboBoxLab.currentText()
-                marker_diameter = self._ui.doubleSpinBoxMarkerDiameter.value()
-                c3d_file = os.path.normpath(os.path.join(input_directory, static_trial))
-                anthropometrics = approximate_anthropometrics(c3d_file, lab, marker_diameter)
-
-                if not self._ui.doubleSpinBoxASISWidth.value():
-                    self._ui.doubleSpinBoxASISWidth.setValue(anthropometrics['inter_asis_distance'])
-                if not self._ui.doubleSpinBoxLeftLegLength.value():
-                    self._ui.doubleSpinBoxLeftLegLength.setValue(anthropometrics['left_leg_length'])
-                if not self._ui.doubleSpinBoxLeftKneeWidth.value():
-                    self._ui.doubleSpinBoxLeftKneeWidth.setValue(anthropometrics['left_knee_width'])
-                if not self._ui.doubleSpinBoxLeftAnkleWidth.value():
-                    self._ui.doubleSpinBoxLeftAnkleWidth.setValue(anthropometrics['left_ankle_width'])
-                if not self._ui.doubleSpinBoxRightLegLength.value():
-                    self._ui.doubleSpinBoxRightLegLength.setValue(anthropometrics['right_leg_length'])
-                if not self._ui.doubleSpinBoxRightKneeWidth.value():
-                    self._ui.doubleSpinBoxRightKneeWidth.setValue(anthropometrics['right_knee_width'])
-                if not self._ui.doubleSpinBoxRightAnkleWidth.value():
-                    self._ui.doubleSpinBoxRightAnkleWidth.setValue(anthropometrics['right_ankle_width'])
+                # Fill in any anthropometric gaps.
+                if self._approximate_anthropometrics:
+                    input_directory = self._ui.lineEditInputDirectory.text()
+                    lab = self._ui.comboBoxLab.currentText()
+                    marker_diameter = self._ui.doubleSpinBoxMarkerDiameter.value()
+                    c3d_file = os.path.normpath(os.path.join(input_directory, static_trial))
+                    anthropometrics = approximate_anthropometrics(c3d_file, lab, marker_diameter)
+
+                    if not self._ui.doubleSpinBoxASISWidth.value():
+                        self._ui.doubleSpinBoxASISWidth.setValue(anthropometrics['inter_asis_distance'])
+                    if not self._ui.doubleSpinBoxLeftLegLength.value():
+                        self._ui.doubleSpinBoxLeftLegLength.setValue(anthropometrics['left_leg_length'])
+                    if not self._ui.doubleSpinBoxLeftKneeWidth.value():
+                        self._ui.doubleSpinBoxLeftKneeWidth.setValue(anthropometrics['left_knee_width'])
+                    if not self._ui.doubleSpinBoxLeftAnkleWidth.value():
+                        self._ui.doubleSpinBoxLeftAnkleWidth.setValue(anthropometrics['left_ankle_width'])
+                    if not self._ui.doubleSpinBoxRightLegLength.value():
+                        self._ui.doubleSpinBoxRightLegLength.setValue(anthropometrics['right_leg_length'])
+                    if not self._ui.doubleSpinBoxRightKneeWidth.value():
+                        self._ui.doubleSpinBoxRightKneeWidth.setValue(anthropometrics['right_knee_width'])
+                    if not self._ui.doubleSpinBoxRightAnkleWidth.value():
+                        self._ui.doubleSpinBoxRightAnkleWidth.setValue(anthropometrics['right_ankle_width'])
 
         # Check height units.
         height = self._ui.doubleSpinBoxHeight.value()
-        if height < 250:
-            self._ui.doubleSpinBoxHeight.setValue(height * 10)
+        self._ui.doubleSpinBoxHeight.setValue(height_to_mm(height))
 
 
     def _clear_subject_info(self):
@@ -586,10 +583,6 @@ class MainWindow(QMainWindow):
             'Right Leg Length': self._ui.doubleSpinBoxRightLegLength.value(),
         }
 
-        # Cap age at 18.
-        # Current shape model is paediatric only, ages over 18 can give strange results.
-        static_data['Age'] = min(static_data['Age'], 18)
-
         missing = [key for key, value in static_data.items() if not value]
         if 'Left Leg Length' in missing or 'Right Leg Length' in missing:
             reply = QMessageBox.warning(self, "Warning", "Leg length measurements missing. You can continue "
@@ -604,10 +597,6 @@ class MainWindow(QMainWindow):
             QMessageBox.warning(self, "Warning", "Subject measurements missing:\n- " + "\n- ".join(missing))
             return
 
-        if static_data['Height'] < 300:
-            QMessageBox.warning(self, "Warning", "Please ensure subject height is in mm.")
-            return
-
         optimise_knee_axis = self._optimise_knee_axis
         if not dynamic_trials:
             reply = QMessageBox.information(self, "No Dynamic Trials",
@@ -685,7 +674,6 @@ class MainWindow(QMainWindow):
 
     def _update_progress(self, message, color):
         self._ui.labelProgress.setText(message)
-        self._progress_text = message
         self._ui.labelProgress.setStyleSheet(f"color: {color};")
 
         self._progress_value = min(self._progress_value + 20, 100)
@@ -1010,6 +998,7 @@ class MainWindow(QMainWindow):
         dlg.exec()
         self._reset_lab_combo_box()
 
+    @handle_runtime_error
     def _handle_error(self, e):
         raise e
 

{c3d_parser-1.5.4 → c3d_parser-1.6.1/src/c3d_parser.egg-info}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: c3d-parser
-Version: 1.5.4
+Version: 1.6.1
 Summary: C3D parser for gait data harmonisation.
 Author-email: Timothy Salemink <tim.nicolas@outlook.com>, Sally Jack <sallyjaack@gmail.com>
 License-Expression: Apache-2.0
@@ -12,7 +12,7 @@ Requires-Dist: numpy==1.26.4
 Requires-Dist: pandas
 Requires-Dist: c3d==0.5.2
 Requires-Dist: trc-data-reader>=0.4.1
-Requires-Dist: opensim-model-creator==0.8.3
+Requires-Dist: opensim-model-creator==0.10.1
 Requires-Dist: PySide6
 Requires-Dist: matplotlib
 Requires-Dist: mplcursors

{c3d_parser-1.5.4 → c3d_parser-1.6.1}/src/c3d_parser.egg-info/requires.txt

@@ -3,7 +3,7 @@ numpy==1.26.4
 pandas
 c3d==0.5.2
 trc-data-reader>=0.4.1
-opensim-model-creator==0.8.3
+opensim-model-creator==0.10.1
 PySide6
 matplotlib
 mplcursors

c3d_parser-1.5.4/src/c3d_parser/__init__.py

@@ -1,2 +0,0 @@
-
-__version__ = "1.5.4"