Differences

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

--- visual3d:documentation:emg:processing:detecting_muscle_onset [2024/06/14 17:11] – created sgranger
+++ visual3d:documentation:emg:processing:detecting_muscle_onset [2024/07/17 15:45] (current) – created sgranger
@@ Line 1: / Line 1: @@
-|**Language:**|** English**  • [[index.php?title=EMG_Detecting_Muscle_Onset/fr&action=edit&redlink=1|français]] • [[index.php?title=EMG_Detecting_Muscle_Onset/it&action=edit&redlink=1|italiano]] • [[index.php?title=EMG_Detecting_Muscle_Onset/pt&action=edit&redlink=1|português]] • [[index.php?title=EMG_Detecting_Muscle_Onset/es&action=edit&redlink=1|español]] ****|
+====== Detecting Muscle Onset ======
-|===== Contents =====\\ \\ \\ \\ * [[#Introduction|1 Introduction]]\\ * [[#EMG_Onset_using_Event_Threshold|2 EMG Onset using Event_Threshold]]\\   * [[#Pipeline|2.1 Pipeline]]\\   * [[#Pipeline_for_computing_the_number_of_frames_the_signal_is_on|2.2 Pipeline for computing the number of frames the signal is on]]\\ * [[#EMG_Onset_based_on_TKEO|3 EMG Onset based on TKEO]]\\   * [[#Pipeline_2|3.1 Pipeline]]\\ * [[#EMG_Onset_based_on_Automatic_Detection|4 EMG Onset based on Automatic Detection]]|
+||
-===== Introduction =====
+Surface EMG signals typically have a frequency content between 50 and 500 Hz. These signals may have a DC bias. It is common to apply a high pass filter with a cutoff frequency of 50 Hz to EMG signals prior to any other processing to remove the DC bias.
-To determine when a muscle is activated a threshold needs to be chosen from a period of inactivity. When the signal surpasses the threshold, the muscle is turned on. There are three different ways to do so.
+**Note: Many laboratories assume a frequency content between 20 and 500 Hz**
-  - EMG Onset using Event_Threshold
+**Note: Applying a high pass filter will effectively remove any DC bias in the signal.**
-  - EMG Onset based on TKEO
-  - EMG Onset based on automatic detection
-===== EMG Onset using Event_Threshold =====
-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. Two events should be created at the start and beginning of the rest period. 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".
-[[File:Emg2.jpg|{{/images/b/b8/Emg2.jpg?454x327}}]]
-==== Pipeline ====
-The following pipeline first includes a low pass filter, a high pass filter and a moving RMS. The first event named REST1 is created by selecting the first frame indicating the start of the rest. The second event named REST2 is created by selecting the frame number at the end of the rest. A typical threshold is median + 3*standard deviation. The median and standard deviation are then calculated by using the [[Metric_Mean|Metric_Median]] and [[Metric_StdDev|Metric_StdDev]] functions. The [[Event_Threshold|Event_Threshold]] pipeline command is then used to compute the threshold by computing the median + 3*standard deviation.
 <code>
-!Apply a high pass and low pass filter assuming that the surface emg signal contains information
+! By default, Visual3D uses a second order Butterworth Bidirectional Filter (resulting in a
-!from 20 to 500 Hz
+! fourth order filter) for Lowpass/Highpass filtering.
-Highpass_Filter
+! High pass filter
-/SIGNAL_TYPES=ANALOG
+Highpass_Filter
-/SIGNAL_NAMES=EMG2
+/Signal_Types= ANALOG
-! /SIGNAL_FOLDER=ORIGINAL
+/SIgnal_Names= EMG1+EMG2+EMG3
-! /RESULT_SUFFIX=
+/Signal_Folder= ORIGINAL
-/RESULT_FOLDER=RMS
+! /Result_Suffix=
-! /FILTER_CLASS=BUTTERWORTH
+! /Result_Folder=PROCESSED
-/FREQUENCY_CUTOFF=20
+! /Filter_Class=BUTTERWORTH
-/NUM_REFLECTED=0
+/Frequency_Cutoff= 50
-/TOTAL_BUFFER_SIZE=100
+/Num_Reflected= 0
-! /NUM_BIDIRECTIONAL_PASSES=1
+/Total_Buffer_Size=100
-;
+/Num_Bidirectional_Passes=1
-Lowpass_Filter
-/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
-Moving_RMS
-/SIGNAL_TYPES=ANALOG
-/SIGNAL_NAMES=EMG2
-/SIGNAL_FOLDER=RMS
-! /RESULT_SUFFIX=
-/RESULT_FOLDER=RMS
-/NUM_WINDOW_FRAMES=1+0.1*PARAMETER::ANALOG::RATE
-;
-!Create events that define the resting baseline of the muscle
-Event_Explicit
-/EVENT_NAME=REST1
-/FRAME=1
-! /TIME=
-;
-Event_Explicit
-/EVENT_NAME=REST2
-/FRAME=150
-! /TIME=
-;
-!Compute the median and standard deviation of the EMG envelope during rest
-Metric_Median
-/RESULT_METRIC_NAME=_MED
-/APPLY_AS_SUFFIX_TO_SIGNAL_NAME=TRUE
-/RESULT_METRIC_FOLDER=EMG
-/SIGNAL_TYPES=ANALOG
-/SIGNAL_NAMES=EMG2
-/SIGNAL_FOLDER=RMS
-! /SIGNAL_COMPONENTS=ALL_COMPONENTS
-/EVENT_SEQUENCE=REST1+REST2
-/EXCLUDE_EVENTS=
-/GENERATE_MEAN_AND_STDDEV=FALSE
-! /APPEND_TO_EXISTING_VALUES=FALSE
-;
-Metric_StdDev
-/RESULT_METRIC_NAME=_SD
-/APPLY_AS_SUFFIX_TO_SIGNAL_NAME=TRUE
-/RESULT_METRIC_FOLDER=EMG
-/SIGNAL_TYPES=ANALOG
-/SIGNAL_NAMES=EMG2
-/SIGNAL_FOLDER=RMS
-! /SIGNAL_COMPONENTS=ALL_COMPONENTS
-/EVENT_SEQUENCE=REST1+REST2
-/EXCLUDE_EVENTS=
-/GENERATE_MEAN_AND_STDDEV=FALSE
-! /APPEND_TO_EXISTING_VALUES=FALSE
-;
-!Identify the event EMG2_ON based on a threshold crossing
-!of the median + 3 * standard deviation
-Event_Threshold
-/SIGNAL_TYPES=ANALOG
-/SIGNAL_NAMES=EMG2
-/SIGNAL_FOLDER=RMS
-/EVENT_NAME=EMG2_ON
-/SELECT_X=TRUE
-! /SELECT_Y=FALSE
-! /SELECT_Z=FALSE
-! /SELECT_RESIDUAL=FALSE
-/THRESHOLD=METRIC::EMG::EMG2_MED+3*METRIC::EMG::EMG2_SD
-! /FRAME_WINDOW=8
-! /FRAME_OFFSET=0
-/ASCENDING=TRUE
-! /DESCENDING=FALSE
-! /ENSURE_RANGE_FRAMES_BEFORE_THRESHOLD_CROSSING=FALSE
-/ENSURE_RANGE_FRAMES_AFTER_THRESHOLD_CROSSING=TRUE
-! /START_AT_EVENT=
-! /END_AT_EVENT=
-/EVENT_INSTANCE=1
 ;
 </code>
-==== Pipeline for computing the number of frames the signal is on ====
-To calculate how long a signal is on or how long a muscle is active, a signal will be created to be set as zero when below the threshold, and one when the signal is above the threshold.
-<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>
-===== EMG Onset based on TKEO =====
-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 [[Teager_Kaiser_Energy|Teager-Kaiser Energy Operator]] and [[Event_Threshold|Event_Threshold]]. This is not a definitive definition, but simply a starting point for the user. The following figure shows an EMG signal before and after the Teager Kaiser Energy Operator is implemented.
-[[File:Teager_Kaiser.jpg|{{/images/7/73/Teager_Kaiser.jpg?560x353}}]]
-==== Pipeline ====
-The following pipeline follows the same general format as above, yet at the beginning there is a Highpass Filter and a Teager Kaiser Energy Operator instead of a Lowpass filter and a HighPass Filter.
-<code>
-!Apply a high pass and low pass filter assuming that the surface emg signal contains information
-!from 20 to 500 Hz
-Highpass_Filter
-/SIGNAL_TYPES=ANALOG
-/SIGNAL_NAMES=EMG_C
-! /SIGNAL_FOLDER=ORIGINAL
-! /RESULT_SUFFIX=
-/RESULT_FOLDER=PROCESSED
-! /FILTER_CLASS=BUTTERWORTH
-/FREQUENCY_CUTOFF=20
-/NUM_REFLECTED=0
-/TOTAL_BUFFER_SIZE=100
-! /NUM_BIDIRECTIONAL_PASSES=1
-;
-!Apply the Teager-Kaiser Operator
-Teager_Kaiser_Energy
-/SIGNAL_TYPES=ANALOG
-/SIGNAL_NAMES=EMG_C
-/SIGNAL_FOLDER=PROCESSED
-! /SIGNAL_COMPONENTS=
-! /RESULT_FOLDER=PROCESSED
-! /EVENT_SEQUENCE=
-! /EXCLUDE_EVENTS=
-;
-!Create events that define the resting baseline of the muscle
-Event_Explicit
-/EVENT_NAME=REST1
-/FRAME=1
-! /TIME=
-;
-Event_Explicit
-/EVENT_NAME=REST2
-/FRAME=150
-! /TIME=
-;
-!Compute the median and standard deviation of the TKO envelope during rest
-Metric_Median
-/RESULT_METRIC_NAME=_MED
-/APPLY_AS_SUFFIX_TO_SIGNAL_NAME=TRUE
-/RESULT_METRIC_FOLDER=EMG
-/SIGNAL_TYPES=ANALOG
-/SIGNAL_NAMES=EMG_C
-/SIGNAL_FOLDER=PROCESSED
-! /SIGNAL_COMPONENTS=ALL_COMPONENTS
-/EVENT_SEQUENCE=REST1+REST2
-/EXCLUDE_EVENTS=
-/GENERATE_MEAN_AND_STDDEV=FALSE
-! /APPEND_TO_EXISTING_VALUES=FALSE
-;
-Metric_StdDev
-/RESULT_METRIC_NAME=_SD
-/APPLY_AS_SUFFIX_TO_SIGNAL_NAME=TRUE
-/RESULT_METRIC_FOLDER=EMG
-/SIGNAL_TYPES=ANALOG
-/SIGNAL_NAMES=EMG_C
-/SIGNAL_FOLDER=PROCESSED
-! /SIGNAL_COMPONENTS=ALL_COMPONENTS
-/EVENT_SEQUENCE=REST1+REST2
-/EXCLUDE_EVENTS=
-/GENERATE_MEAN_AND_STDDEV=FALSE
-! /APPEND_TO_EXISTING_VALUES=FALSE
-;
-!Identify the event EMG_C_ON based on a threshold crossing
-!of the median + 3 * standard deviation
-Event_Threshold
-/SIGNAL_TYPES=ANALOG
-/SIGNAL_NAMES=EMG_C
-/SIGNAL_FOLDER=PROCESED
-/EVENT_NAME=EMG_C_ON
-/SELECT_X=TRUE
-! /SELECT_Y=FALSE
-! /SELECT_Z=FALSE
-! /SELECT_RESIDUAL=FALSE
-/THRESHOLD=METRIC::EMG::EMG2_MED+3*METRIC::EMG::EMG2_SD
-! /FRAME_WINDOW=8
-! /FRAME_OFFSET=0
-/ASCENDING=TRUE
-! /DESCENDING=FALSE
-! /ENSURE_RANGE_FRAMES_BEFORE_THRESHOLD_CROSSING=FALSE
-/ENSURE_RANGE_FRAMES_AFTER_THRESHOLD_CROSSING=TRUE
-! /START_AT_EVENT=
-! /END_AT_EVENT=
-/EVENT_INSTANCE=1
-;
-</code>
-===== EMG Onset based on Automatic Detection =====
-The following meta-command may be used to process EMG analog signals to find "resting" thresholds automatically, without using a reference file. Processing motion files with multiple EMG signals can become tedious when manually extracting events and identifying points of “rest”. We can manually detect sections of minimal activity within a motion trial to set as a reference threshold. The pipeline to compute this automated detection has been provided as a meta-command, and you may copy and paste it into your **Visual3D x64 > Plugins > Meta-Commands** folder. For more information on Meta-Commands and how they work, see [[Pipeline_Commands:Meta_Commands|Pipeline Commands:Meta Commands]]. For a detailed tutorial on how to use Automatic_EMG_Events, head to the tutorial page:
-**Meta-Command**:
-**Tutorial page**: [[Visual3D:Tutorials:EMG:Automatic_EMG_Events|Tutorial: Automatic EMG Events]]
-Retrieved from ""