Differences

This shows you the differences between two versions of the page.

--- visual3d:documentation:emg:processing:onset [2024/06/14 17:11] – created sgranger
+++ visual3d:documentation:emg:processing:onset [2024/07/17 15:45] (current) – created sgranger
@@ Line 1: / Line 1: @@
-|**Language:**|** English**  • [[index.php?title=EMG_Onset/fr&action=edit&redlink=1|français]] • [[index.php?title=EMG_Onset/it&action=edit&redlink=1|italiano]] • [[EMG_Onset/pt|português]] • [[index.php?title=EMG_Onset/es&action=edit&redlink=1|español]] ****|
+====== Onset ======
-|===== Contents =====\\ \\ \\ \\ * [[#Count_the_Number_of_Frames_the_signals_is_ON|1 Count the Number of Frames the signals is ON]]|
+||
-Computing the onset of an EMG signal can be quite tricky. In this example, we describe a pipeline for a straightforward estimation of EMG onset based on the EMG envelope and [[Event_Threshold|Event_Threshold]]. This is not a definitive definition, but simply a starting point for the user.
+Computing the onset of an EMG signal can be quite tricky. In this example, we describe a pipeline for a straightforward estimation of EMG onset based on the the [[Visual3D:Documentation:Pipeline:Signal_Commands:Teager_Kaiser_Energy|Teager-Kaiser Energy]] Operator and [[Visual3D:Documentation:Pipeline:Event_Commands:Event_Threshold|Event_Threshold]]. This is not a definitive definition, but simply a starting point for the user.
+{{:Teager_Kaiser.jpg}}
+\\
-[[File:Emg2.jpg|{{/images/b/b8/Emg2.jpg?454x327}}]]
-Three event labels are identified on this figure. The yellow and red events specify the "rest period". This is arbitrary and is specified by the user explicitly. The blue event label is the first threshold crossing above the "median + 3*standard deviation".
 <code>
@@ Line 14: / Line 15: @@
 Highpass_Filter
 /SIGNAL_TYPES=ANALOG
-/SIGNAL_NAMES=EMG2
+/SIGNAL_NAMES=EMG_C
 ! /SIGNAL_FOLDER=ORIGINAL
 ! /RESULT_SUFFIX=
-/RESULT_FOLDER=RMS
+/RESULT_FOLDER=PROCESSED
 ! /FILTER_CLASS=BUTTERWORTH
 /FREQUENCY_CUTOFF=20
@@ Line 25: / Line 26: @@
 ;
-Lowpass_Filter
+!Apply the Teager-Kaiser Operator
-/SIGNAL_TYPES=ANALOG
-! /SIGNAL_NAMES=
-/SIGNAL_FOLDER=RMS
-! /RESULT_SUFFIX=
-/RESULT_FOLDER=RMS
-! /FILTER_CLASS=BUTTERWORTH
-/FREQUENCY_CUTOFF=500
-/NUM_REFLECTED=0
-/TOTAL_BUFFER_SIZE=100
-! /NUM_BIDIRECTIONAL_PASSES=1
-;
-!Compute the RMS of the emg signal based on a 100 ms window
+Teager_Kaiser_Energy
-Moving_RMS
 /SIGNAL_TYPES=ANALOG
-/SIGNAL_NAMES=EMG2
+/SIGNAL_NAMES=EMG_C
-/SIGNAL_FOLDER=RMS
+/SIGNAL_FOLDER=PROCESSED
-! /RESULT_SUFFIX=
+! /SIGNAL_COMPONENTS=
-/RESULT_FOLDER=RMS
+! /RESULT_FOLDER=PROCESSED
-/NUM_WINDOW_FRAMES=1+0.1*PARAMETER::ANALOG::RATE
+! /EVENT_SEQUENCE=
+! /EXCLUDE_EVENTS=
 ;
@@ Line 63: / Line 52: @@
 ;
-!Compute the median and standard deviation of the EMG envelope during rest
+!Compute the median and standard deviation of the TKO envelope during rest
 Metric_Median
@@ Line 70: / Line 59: @@
 /RESULT_METRIC_FOLDER=EMG
 /SIGNAL_TYPES=ANALOG
-/SIGNAL_NAMES=EMG2
+/SIGNAL_NAMES=EMG_C
-/SIGNAL_FOLDER=RMS
+/SIGNAL_FOLDER=PROCESSED
 ! /SIGNAL_COMPONENTS=ALL_COMPONENTS
 /EVENT_SEQUENCE=REST1+REST2
@@ Line 84: / Line 73: @@
 /RESULT_METRIC_FOLDER=EMG
 /SIGNAL_TYPES=ANALOG
-/SIGNAL_NAMES=EMG2
+/SIGNAL_NAMES=EMG_C
-/SIGNAL_FOLDER=RMS
+/SIGNAL_FOLDER=PROCESSED
 ! /SIGNAL_COMPONENTS=ALL_COMPONENTS
 /EVENT_SEQUENCE=REST1+REST2
@@ Line 93: / Line 82: @@
 ;
-!Identify the event EMG2_ON based on a threshold crossing
+!Identify the event EMG_C_ON based on a threshold crossing
 !of the median + 3 * standard deviation
 Event_Threshold
 /SIGNAL_TYPES=ANALOG
-/SIGNAL_NAMES=EMG2
+/SIGNAL_NAMES=EMG_C
-/SIGNAL_FOLDER=RMS
+/SIGNAL_FOLDER=PROCESED
-/EVENT_NAME=EMG2_ON
+/EVENT_NAME=EMG_C_ON
 /SELECT_X=TRUE
 ! /SELECT_Y=FALSE
@@ Line 118: / Line 107: @@
 </code>
-\\
-===== Count the Number of Frames the signals is ON =====
-<code>
-! Create a signal that is zero when the EMG is off, and one when it is ON
-! Start with creating a signal that is zero at all frames
-Evaluate_Expression
-/EXPRESSION=FRAME_NUMBERS::ORIGINAL::FRAMES * 0
-! /SIGNAL_TYPES=
-! /SIGNAL_FOLDER=ORIGINAL
-! /SIGNAL_NAMES=
-/RESULT_TYPES=DERIVED
-/RESULT_FOLDERS=EMG_COUNT
-/RESULT_NAME=EMG2
-! /APPLY_AS_SUFFIX_TO_SIGNAL_NAME=FALSE
-;
-! Set all frames between the ON/OFF events to one
-Set_Data_To_New_Values
-/SIGNAL_TYPES=DERIVED
-/SIGNAL_FOLDER=EMG_COUNT
-/SIGNAL_NAMES=EMG2
-! /SIGNAL_COMPONENTS=ALL_COMPONENTS
-/RESULT_TYPES=DERIVED
-/RESULT_FOLDERS=EMG_COUNT
-/RESULT_SUFFIX=
-/EVENT_SEQUENCE=EMG2_ON+EMG2_OFF
-! /EXCLUDE_EVENTS=
-! /EVENT_INSTANCE=0
-! /START_FRAME=
-! /END_FRAME=
-! /USE_POINT_RATE=TRUE
-/REPLACEMENT_VALUES=1
-! /THRESHOLD_HIGH=NO_DATA
-! /THRESHOLD_LOW=NO_DATA
-! /SET_DATA_AT_SINGLE_EVENT=FALSE
-;
-! Add the signal during the gait cycle
-!   The result will be the number of frames the signal is ON during the gait cycle
-Metric_Sum
-/RESULT_METRIC_FOLDER=EMG_COUNT
-/RESULT_METRIC_NAME=
-/APPLY_AS_SUFFIX_TO_SIGNAL_NAME=TRUE
-/SIGNAL_TYPES=DERIVED
-/SIGNAL_FOLDER=EMG_COUNT
-/SIGNAL_NAMES=EMG2
-/COMPONENT_SEQUENCE=ALL
-/EVENT_SEQUENCE=RHS+RHS
-/EXCLUDE_EVENTS=
-/SEQUENCE_PERCENT_START=
-/SEQUENCE_PERCENT_END=
-! /GENERATE_MEAN_AND_STDDEV=TRUE
-! /APPEND_TO_EXISTING_VALUES=FALSE
-;
-</code>
-Retrieved from ""