From b34467971868f03c2e3b0d3130bd676c0b16ac11 Mon Sep 17 00:00:00 2001
From: LexvanGastel <l.vangastel@studen.utwente.nl>
Date: Wed, 2 Apr 2025 11:05:44 +0200
Subject: [PATCH] vraag 18
---
Q2main.py | 26 ++++++++++++++++++++++++++
main.py | 21 ++++++++++++++++++++-
2 files changed, 46 insertions(+), 1 deletion(-)
diff --git a/Q2main.py b/Q2main.py
index e69de29..af4ad41 100644
--- a/Q2main.py
+++ b/Q2main.py
@@ -0,0 +1,26 @@
+
+class ImageProcessor:
+ def __init__(self, image_path):
+ self.image
+ self.image_rgb
+ self.image_hsv
+
+ def apply_median_filter(self, kernel_size):
+ return filterd_image
+
+ def get_color_channels(self):
+ return red_channel, green_channel, blue_channel
+
+ def convert_to_hsv(self):
+ return hsv_image
+
+class Plotter:
+ def __init__(self):
+ pass
+
+ def plot_image_channels(self, image):
+ return None
+
+ def plot_histograms(self, image):
+ return None
+
\ No newline at end of file
diff --git a/main.py b/main.py
index 035bb00..76c6194 100644
--- a/main.py
+++ b/main.py
@@ -1,6 +1,6 @@
import scipy.io
import numpy as np
-from scipy.signal import butter, filtfilt
+from scipy.signal import butter, filtfilt, find_peaks
import matplotlib.pyplot as plt
def load_emg_data(file_path):
@@ -130,6 +130,25 @@ def plot_emg_data(file_path):
freq_rf, fft_rf = compute_fft(envelope_rf, sampling_rate)
freq_ta, fft_ta = compute_fft(envelope_ta, sampling_rate)
+ peaks_rf, _ = find_peaks(envelope_rf)
+ peaks_ta, _ = find_peaks(envelope_ta)
+
+ num_peaks_rf = len(peaks_rf)
+ num_peaks_ta = len(peaks_ta)
+
+ max_value_rf = np.max(envelope_rf)
+ max_value_ta = np.max(envelope_ta)
+
+ area_rf = np.trapz(envelope_rf, time_vector)
+ area_ta = np.trapz(envelope_ta, time_vector)
+
+ print(f"1.Numbers of EMG peaks Rectus Femoris: {num_peaks_rf}")
+ print(f"1.Numbers of EMG peaks Tibialis Anterior: {num_peaks_ta}")
+ print(f"2. Maximale peak Rectus Femoris: {max_value_rf}")
+ print(f"2. Maximale peak Tibialis Anterior: {max_value_ta}")
+ print(f"3. Area under Curve Rectus Femoris {area_rf}")
+ print(f"3. Area under Curve Tibialis Anterior {area_ta}")
+
plt.figure(figsize=(13,8))
plt.subplot(2, 2, 1)
--
GitLab