DEMUSEtool – user manual
Contents
1. Preferences age Setup Print Print Selection 1 E grams 4boutDialog m 2 plotMUAP_3DMap_black m 3 E rograms MUAPmap2D m 4E rograms MUAPmap3D m howe Up hove Down Mowe to Bottom Exit MATLAB Ctrl Q Remove Save Close Revert Default Help Figure 2 Setting the path in matlab programme environment in this case DEMUSE tool was copied to the directory c DEMUSEtool Page 5 of 34 DEMUSE tool User Manual DEMUSE 4 Using DEMUSE tool 4 1 Starting DEMUSE tool To start the DEMUSEtool type the following command to matlab command window gt gt DEMUSEtool The main DEMUSEtool window appears Figure 3 This window comprises the following four groups of commands for e loading band pass filtering and visualization of the acquired sEMG signals e decomposition of acquired SEMG signals e graphical plots and animations of the decomposition results e saving and reloading of the decomposition results Channel selection frame loads and visualizes DEMUSE toal a E i 1 TULEE rites 2051 500 gt 5 ce ok Ron decomposton sEMG signals Rucni 7 displays _Plot MU discharges _ Plot discharge rates _ __seMGedior 8 graphs and Plot MUAPs Plot MUAP trains Mut alf animations Animate 2D MUAP map Animate 3D MUAP map Frames s 10 saves and About bewuse ff reloads the ee eae results Figure2. Time s Figure 8 Matlab figure with selected sEMG channels Displayed figures can be manipulated by using standard matlab graphical tools for zooming in out for saving and printing the figure Figure 9 Zoomed in version of Figure 8 is depicted in Figure 10 Figure No 1 o File Edik View Insert Tools Window Help JDSUS RL AAS PED Figure 9 Matlab figure toolbar with tools for zooming in out on a figure and for saving and printing the figure Page 9 of 34 DEMUSE tool User Manual DEMUSE O x lt Figure No 1 7 aa File Edit View Insert Tools Window Help De eHSG KAAS PPD Q r 2 8 3 Time s Figure 10 Zoomed in version of Figure 8 Figures can be closed by clicking on a corresponding buttons in the top right corner of each figure Figure 11 Figure 11 Buttons for minimization maximization and closing of the figure Page 10 of 34 DEMUSE tool User Manual DEMUSE 4 4 SEMG signal decomposition Decomposition of SEMG signals is fully automatic The user specifies only the number of decomposition runs Figure 12 and optionally the initial SEMG signal offset Figure 13 DEMUSEtool uses the gradient Convolution Kernel Compensation gCKC decomposition technique 3 gCKC is a sequential motor unit MU identification method and requires one iteration run per each reconstructed MU The user can predefine the number of iterations by moving th
3. ARA ES A A MU 1 MU 1 Amplitude mV Amplitude mV Figure 37 Pseudocolor animation of a MUAP generation propagation and attenuation process MUAP amplitudes on different sEMG channels black amp white circles specify the colours of corresponding points on the 2D map Colours of intermediate map points are calculated by the bilinear interpolation of the MUAP amplitudes in four adjacent SEMG channels Page 24 of 34 DEMUSE tool User Manual DEMUSE 4 6 5 3D MUAP map animation The second animation provided by the DEMUSEtool includes a 3D plot of MAUP amplitude in time By analogy with the animation of 2D MUAP map MUAP templates are first estimated by a spike triggered averaging of SEMG channels SEMG channels are then Spatially organized into a discrete 2D map reflecting the relative position of pick up electrodes The amplitudes of MUAP templates at a given time instant specify the height on this 2D map of channels Missing intermediate points on the map are calculated by bilinear interpolation of MUAP amplitudes in four adjacent SEMG channels In the next animation frame the time is moved forward by one signal sample and the 3D map is recalculated To start the 3D animation select the MU and the animation frame rate Figure 38 right panel After clicking on the Animate 3D MUAP map button Figure 38 left panel the animation window opens Figure 39 and the animation automatically starts The animation begins
4. DEMUSE l ee m en ats Bom I egaa Piot MUAP trains MUI al Frames fs 10 ELE Preep W Filler z0 0 HE LEE 8 2 8 E 8 8 fe a8 8 8 amp amp 8 8 4 2 amp T 8 4 8 2 E a 8 8 8 8 DEMUSEtool user manual Version 2 3 02 10 2007 Page 1 of 34 DEMUSE tool User Manual CONTENTS C EF OCUIC HON e E EE E EEE E E EE EE T E EE E EEE E E E E E EET 3 2 DEMU S EDOCONPON SA eaa E E E E E EE E 3 So DEMUSEHoNNStaNalO i sescctccvuccccsncsenansncdscvsncgcoepenasebamenearsenssenceeneoaoebepebeseeonseeadeensesaeeceenasreonieses 4 3 1 REGUETON E T A T 4 32 DEMUSE files and 6 6 2k ennn EAAS ee eee eee eee 5 tU DEMUS TOO eeine reenter nn ne ene ee ee ee eee 6 Ts an DEMOS E 0 l0 ee E nee ee 6 42 SENG SIQMALMOAGING ncorccnecassnsvesssenensdeeccewancanesnageeonsmageeeneoasebapanagenspenepeb epepeneesbenebebenanaseeenseans T AS SEMG signal visualization vesavetsdecocncdencdonsdenseebobanebancdancdcnodenddenstnnchovelunebsnetunsberewsnabcontuachensdes 8 4 4 SEMG signal CECOMPOSITION 0 ccc cccc cece eeeceeeceeeceeeseeseeseeeceeeseeeseeeseecseeseeeseeeseeeaeesanss 11 4 5 Visualisation and editing of the decomposition FESUItS ccccecceeceeeeeeeeeseeeeeeeeeeeeeeeeees 13 4 5 1 FOUNO SEMO SINA en ee ee 14 4 5 2 Displaying instantaneous MU discharge rate ccccccecceecseeceeeeeeseeeeeeseeeeeeseeseeenes 15 4 5 3 Ad
5. a 2D map heights of intermediate map points are calculated by the bilinear interpolation of the MUAP amplitudes in four adjacent SEMG channels Note Animation window cannot be closed during the run of animation Stop the f animation by clicking on or button before you close the figure During the 3D animation the axes of the 3D plot can be freely rotated To rotate a 3 D axes click on the axes and drag the cursor in the direction you want to rotate When you release the mouse button DEMUSEtool redraws the axes in the new orientation Figure 40 Page 26 of 34 DEMUSE tool User Manual DEMUSE MUAPmap3D a S E Q ga z L FaN lt MUAPmap3D F Amplitude mV Figure 40 3D rotation of the axes the axes of the 3D plot can be freely rotated by clicking on the axes and dragging the cursor in the direction you want to rotate the plot Page 27 of 34 DEMUSE tool User Manual DEMUSE 4 6 6 Plots of reconstructed MUAP trains DEMUSEtool provides tool for plotting the sum of reconstructed MUAP trains Superimposed to the original SEMG signals This proves beneficial when evaluating the efficiency of the decomposition process In surface EMG there are many small and deep MUs which cannot be recognized They contribute the background physiological noise The second source of noise is so called instrumentation or thermal noise which originates from the instrumentation s parasite capacities line in
6. approx 5 ms before the actual generation of the MUAP and ends approx 5 ms after the MUAP attenuation During the animation the propagation of MUAP along the muscle fibers can be observed Figure 39 Buttons on the top of the animation window Figure 39 enable the following actions e gt button re plays the animation e button pauses the animation e button stops the animation e lt button animates the previous animation frame i e step backward e s animates the next animation frame i e step forward The current animation frame is displayed in the top right corner of the animation window Figure 39 The animation window is a regular matlab figure and can be freely manipulated by all available matlab graphical tools e g Zooming coping printing etc _ annn eee een MU1 Ji ai map Animate 3D MUAP map Framesis 10 gt 1 F gt 1 Figure 38 Animate 3D MUAP map button left and MU selection and frame rate sliders right Page 25 of 34 DEMUSE tool User Manual DEMUSE E E Pmap3D a ol x Eee E E EARTE MU 1 MU 1 Amplitude mV Amplitude mV a ap3D a aS aS eae LSS Seeks see MU 1 MU 1 So o o xo xo St 500 A lt L 0 lt L Figure 39 3D animation of MUAP generation propagation and attenuation MUAP amplitudes on different SEMG channels red circles specify the height of corresponding points on
7. channel residual signal yellow line is clearly visible when compared to Figure 27 Note Several MUAP occurrences can be simultaneously selected This saves the f user s time and energy in the case of intensive manual editing Page 18 of 34 DEMUSE tool User Manual DEMUSE 4 5 5 Selecting sEMG channel spatial filter and cut off frequencies The sEMG channel to be displayed can be selected by clicking on a radio button representing the corresponding channel in the left bottom corner of the sEMG editor window Figure 29 Electrode rows and columns are denoted by red numbers displayed at left bottom of the channel selection panel Text labels right above the sEMG channel selection panel display cut off frequencies of the 4th order Butterworth band pass filter Cut off frequencies can be modified by typing into the aforementioned text labels Figure 29 The sEMG signals and corresponding MUAP templates are automatically recalculated when cut off frequencies are changed Default cut off frequencies values are set to 20 and 500 Hz respectively Ga eee eee eee eee Electrode configuration z ggde L OATIQUrAaION S ilter 20 510 Hz Filter 2i 500 Hz LEE SE ER ELER REESE SESE OE LER ER ER ELSE SESE SESE OE 2o e e e e a e e ESL SE ses E Ei fe es es g oe X 2 S S se E E so 8 se g ses fe bta hi es es one Figure 29 Command fr
8. of a MUAP occurrence must be explicitly confirmed by clicking on add selected MUAP button Page 17 of 34 DEMUSE tool User Manual DEMUSE 4 5 4 Moving and deleting MUAP occurrences Each MUAP occurrence displayed in the top panel can be manually moved to the left moved to the right and deleted respectively The user must first click on the displayed MUAP occurrence in the top panel in order to select it The selected MUAP occurrence is denoted by a thick red line Figure 27 Afterwards the selected MUAP occurrence can be moved or deleted by clicking on a lt lt move selected MUAPs move selected MUAPs gt gt and delete selected MUAPs button respectively Figure 27 Time 5 0 65 0 7 0 75 0 6 delete selected MU move left selected move right selected discharge MUAP MU discharge MUAP MU discharge MUAP Figure 27 Manipulation of selected MUAP occurrence depicted by red thick line The selected MUAP occurrence can be moved left right or deleted Each user action results in an immediate update of a signal residual and MUAP templates Figure 28 Selected MUAP occurrence can be deselected by a single click 7 2 4 3 i 2 4 3 7 5 5 8 6 i add selected MUAPs delete selected MUAPs lt lt move selected MUAPs move selected MUAPs gt gt ASIANA TOn SEATA O oo Ea a ARS Figure 28 Same as in Figure 27 with the selected MUAP occurrence moved to the left Increase of the sEMG
9. 3 Main DEMUSEtool window with explanations of different command groups Each group of commands is depicted in different colour Figure 3 Channel selection frame displayed on the right hand site in Figure 3 allows selection of different rows or columns of acquired SEMG channels Its detailed description is provided in Subsections 4 2 4 3 and 4 6 6 Page 6 of 34 DEMUSE tool User Manual DEMUSE 4 2 SEMG signal loading To load acquired sEMG signals into DEMUSEtool click on Load SIG file button Figure 4 left panel Load SIG file dialog window appears Figure 4 right panel Chose the SIG file and click on Open button Load 5IG file Look in 6 OOMaGaLi0705151518 2 d OOMaGaLi070515151813 sig OO0MasaLli070515151814 sig OOMasaLli070515151815 sig Files of type SIG Cancel Figure 4 Load SIG file button left panel and Load SIG file dialog window right panel DEMUSE tool automatically loads all the required technical information e g sampling frequency dimensions of the acquisition system electrode configurations etc from the corresponding measurement session abstract file 2 The number of acquired channels and their relative spatial configuration is displayed in Channel selection frame Figure 5 Figure 5 Channel selection frame the number of acquired channels and their relative spatial organization is displayed grey rectangles Red numbers
10. 31 49 Ll Channel MINA a bn ihi 2 8 Time s Figure 43 Matlab figure with selected sSEMG channels and corresponding MUAP trains short signal segment from the signal shown in Figure 42 Page 30 of 34 DEMUSE tool User Manual DEMUSE 4 7 Saving and reloading of the decomposition results Decomposition results can be saved by clicking on the Save results button Figure 44 The results are automatically saved into the directory containing the currently loaded SIG file The following file naming convention is used NameOfTheSIGfile offsetN runsM mat the name of the Initial Number of currently loaded signal decomposition SIG file offset iterations where NameOfTheSIGfile stands for the name of the currently loaded SIG file N is the initial signal offset in seconds see Section 4 4 for details and M is the number of decomposition runs see Section 4 4 The decomposition results of a SIG file Subject1 SIG with initial signal offset set equal to 0 and number of decomposition iterations set equal to 30 is saved in the following matlab file pUDTeCEL OffSelL Luns3sU cmat Saved results can be reloaded by clicking on the Load results button Figure 44 Load results dialog window opens Figure 45 Choose the mat file and click on Open button Once reloaded into the DEMUSEtool results can be freely edited and displayed graphical representations and animations of the reloade
11. P templates of different MUs as estimated by spike triggering averaging of the selected sEMG channel Figure 22 The identified MU discharges are taken as trigger The length of the averaging window is set equal to 10 ms addselected MUAPs delete selected MUAPs MUAP template of MU 1 MUAP template of MU 8 estimated from the estimated from the selected sSEMG channel selected sEMG channel Figure 22 The central panel of sEMG editor window with MUAP templates MUAP templates are estimated by spike triggered averaging of the selected SEMG channel The central panel allows the user to select a particular MU by clicking on corresponding MUAP template The MUAP occurrences of selected MUs are automatically displayed in the top panel while the bottom panel displays the instantaneous discharge rate plot Figure 24 Several MUs can be simultaneously selected Figure 25 Each circle in the bottom panel corresponds to a single MU discharge Figure 23 The horizontal position of the circle denotes the MUAP occurrence time whereas its vertical position reflects instantaneous MU discharge rate calculated as the quotient between the sampling frequency and the inter pulse interval preceding the selected MU discharge The position and the length of the signal window displayed in the top panel is also depicted grey rectangle User can move to the time of a particular MU discharge by simply clicking on the correspondin
12. ame of sEMG editor with sEMG channel panel Butterworth band pass filter cut off frequencies and Electrode configuration panel for single differential electrode configuration right and Laplacian electrode configuration left Above the filter cut off frequencies there is the Electrode configuration panel Figure 29 It allows the user to choose among different spatial filters applied to acquired SEMG signals The following four spatial filters are currently supported monopolar longitudinal single differential default longitudinal double differential and Laplacian sEMG channel selection panel SEMG channel displayed in the top panel and corresponding MUAP templates are automatically updated whenever a new spatial filter is selected Figure 29 Page 19 of 34 DEMUSE tool User Manual DEMUSE 4 6 Graphical results DEMUSEtool includes several tools for graphical representation of the decomposition results The user can plot the discharge patterns of reconstructed MUs instantaneous MU discharge rates multichannel MUAPs and reconstructed MUAP trains In addition MUAP generation propagation and attenuation can be animated for each identified MUs All the graphical results are depicted in matlab figures and can be easily manipulated by standard matlab s editing tools In the sequel description of each aforementioned graphical representation is provided 4 6 1 MU discharge patterns plot MU discharg
13. at the left hand site and at the bottom denote the corresponding rows and columns of electrodes respectively Page 7 of 34 DEMUSE tool User Manual DEMUSE DEMUSEtool uses 4 order Butterworth band pass filter to filter the raw EMG signals Filter s cut off frequencies can be controlled by typing new values into the text labels shown in Figure 6 Default cut off frequencies are set to 20 and 500 Hz respectively Fitter 20 FJ 500 Figure 6 Text labels controlling the cut off frequencies of built in Butterworth band pass filter 4 3 SEMG signal visualization To display loaded and band pass filtered SEMG signals select first the corresponding electrode row or electrode column Figure 7 left panel and click on Plot signals button Figure 7 right panel Matlab figure with selected sEMG channels appears Figure 8 Note Due to the large number of acquired SEMG channels only selected row column of f SEMG channels can be displayed in one figure The number of figures however is g not limited You can display all the sEMG channels by consecutively selecting the different electrode columns for example 1 3 Figure 7 Channel selection frame left panel and plot signals button right panel 6 electrode row is selected Page 8 of 34 DEMUSE tool User Manual DEMUSE J Figure No 1 a a File Edit Yiew Insert Tools Window Help DeHBkRAAS PEA i E Channel ERR FAP 1 2 3 4 5 6
14. constructed in the 13 iteration of the gCKC decomposition technique Each pulse corresponds to a single MU discharge Average discharge rate and coefficient of variability of inter pulse interval are outside the expected range of values Thus the selected MU discharges are depicted by red circles while the reconstructed discharge pattern is discarded ee Sr Se a Page 12 of 34 DEMUSE tool User Manual DEMUSE When the decomposition ends the total processing time is displayed Figure 17 RUN Ge COM poOSTion 7 ne Te n Se y A F i k m gf F rer Ira y ie rn eB Figure 17 Text label displaying the total processing time 4 5 Visualisation and editing of the decomposition results Decomposition results can be visually inspected and edited by sEMG editor Figure 19 To open sEMG editor click on sEMG editor button Figure 18 t i ic SLMG editor ac Figure 18 sEMG editor button lt SEMGeditor 1 0 x File Edit View Insert Tools Window Help Time s 0 26 0 28 0 3 0 52 0 34 0 36 0 38 0 4 0 42 0 44 0 46 ain a 4 2 1 2 5 6 lt lt CILCICICICICICICICIrCICIrI CICIKCICILCICICICICIrITI ele ee eee ee ee ele ee ee 2 ee Figure 19 sEMG editor window Page 13 of 34 DEMUSE tool User Manual DEMUSE 4 5 1 Plotting sEMG signals SEMG editor window comprises three panels Figure 19 The top panel
15. d results are fully supported Figure 44 Save results button left and Load results button right Look ir E signals do fe Subjecti_offsetO_rums30 mati Files of type MAT files rat Cancel Page 31 of 34 DEMUSE tool User Manual DEMUSE Figure 45 Load results dialog window 4 8 DEMUSEtool acknowledgements Information about the DEMUSEtool version copyrights and authors acknowledgement are displayed by clicking on the About DEMUSE button Figure 46 About DEMUSE dialog window opens Figure 47 13 ij I About DEMUSE Figure 46 About DEMUSE button AboutDialog MARIE CURIE ACTIONS Figure 47 About DEMUSE dialog window Page 32 of 34 DEMUSE tool User Manual DEMUSE 5 Technical support DEMUSEtool is copyrighted by the System Sofware Laboratory SSL from University of Maribor Slovenia and Laboratory of Engineering of Neuromuscular System and Motor Rehabilitation LISIN from Politectnico di Torino Italy Its development was supported by a Marie Curie Intra European Fellowship within the 6th European Community Framework Programme DE MUSE Contract No 023537 For further technical assistance and support please contact ales holobar delen polito it or ales holobar uni mb si System Software Laboratory SSL 2007 University of Maribor Slovenia Laboratory of Engineering of Neuromuscular System and Motor R
16. ding anew MUAP occurrence ccccceccescceeceeeceeceeceeeceeseeeceecuecauecaesaeecseseesaueeas 17 4 5 4 Moving and deleting MUAP OCCUITENCES ccccccecceecceeeeeeeeeceeeseeseeeseeseeeeeeseesaeeaes 18 4 5 5 Selecting SEMG channel spatial filter and cut off frequencies cccceeee eee eeees 19 Ae OPCIE UNS e e EEEE EENEN EEE AE EENE E EAE EE 20 4 6 1 MU discharge DANCING PlOleoicssreripen ene EEE 20 4 6 2 Instantaneous discharge rate PlOts cccccccecceecceeceeceeeceeceeeceecueceeeseeseeeseesessaeeeas 21 4 6 3 MOUCHEL MUAF POLS aeon e E E E E E 22 4 6 4 PAB MUAP Map animalON eee EE 23 4 6 5 OL MUAF Map SE AUTO iser eren EE 25 4 6 6 Plots of reconstructed MUAP trains cccccccceccceeeseeeeeeeeeece esse eeseeeaeesseeseeeseeeseeaes 28 4 7 Saving and reloading of the decomposition rESUITS cece c cece ceeeeeeeeeeeeeeseeeseeeseeeaes 31 4 8 DEMUSEtool acknowledgements ccccccccsecceeeceeeeeeceeeceeeceeeseeesaeeseesegeseueseeeseeeseeeaes 32 DS LOCMNICAl SUDDOT Es acccezccrencsencnsnancennsnneanunensceanespeesespemeasnaeemaneaansacgabaaeneseasesmeaeespennesbapemaneeentacetbaneer 33 polka o 0 o gt ee ene eee ee ee een ee ne ee ee ete eee ee errr 34 Page 2 of 34 DEMUSE tool User Manual DEMUSE 1 Introduction DEMUSEtool is a matlab 1 program for visualization and decomposition of multichannel surface electromyograms SEMG acquired by EMG acquisition softwar
17. displays original SEMG channel blue line and the residual after subtraction of identified MUAPs yellow line The ID of MU discharging at a particular time moment is also depicted Figure 20 ID of MU discharging at this particular time instant Time s og 0 2 03 0 32 0 34 0 36 0 38 0 4 0 42 0 44 O46 residual after subtraction of identified MUAPs original SEMG channel Figure 20 The top panel of sEMG editor window with original SEMG channel blue line and residual after subtraction of the identified MUAPs yellow line Buttons at the left and right hand side of the top panel determine the position and the size of the time window being depicted Buttons lt lt and gt gt move the displayed window left and right respectively Figure 21 Button zooms in out on a displayed signals decreases increases the length of the time window move left move right Tiny s 06 u7 O08 09 1 11 1 2 l 13 l i i Pro E S ba he _ delete sBtested MUAPs lt lt move selected MUAPs mgyers lected MUAPs gt Electrode configuration zoom in zoom out Figure 21 The top panel of sEMG editor window with buttons for moving left right and zooming in out on displayed signals Page 14 of 34 DEMUSE tool User Manual DEMUSE 4 5 2 Displaying instantaneous MU discharge rate The central panel of the sEMG editor window displays the MUA
18. e 37 The MU and the animation frame rate are selected by the sliders shown in Figure 36 right panel After clicking on the Animate 2D MUAP map button Figure 36 left panel animation window opens Figure 37 and the animation starts automatically The animation begins approx 5 ms before the actual generation of the multi channel MUAP and ends approx 5 ms after the MUAP attenuation During the animation the propagation of MUAP along the muscle fibers can be observed Figure 37 Buttons on the top of the animation window Figure 37 enable the following actions e gt button re plays the animation e button pauses the animation e button stops the animation e lt button animates the previous animation frame i e step backward e lt animates the next animation frame i e step forward Current animation frame is displayed in the top right corner of the animation window Figure 37 eee ES hs 11 Animate 2D MUAP map E Frames s 10 E p T Figure 36 Animate 2D MUAP map button left and sliders for selection of MU and frame rate right Note Animation window cannot be closed while the animation is running Stop the f animation by clicking on or button before you close the figure Page 23 of 34 DEMUSE tool User Manual DEMUSE ymuarmen2 gt z gt MuAPmapzD E A MU 1 MU 1 gt E V o ga ga Z Z rot rot X
19. e installed on the system in order to prevent extensive swapping of the memory space 3 1 Requirements Minimal hardware configuration e 1 GHz CPU e 20 MB disk e 512 MB RAM Recommended hardware configuration 2 GHz CPU or higher 100 Mb disk 2 MB L2 cache or more 1 GB RAM To run DEMUSE tool following software should be properly installed e matlab 1 version 6 5 or higher Page 4 of 34 DEMUSE tool User Manual DEMUSE 3 2 DEMUSE files and folders DEMUSEtool comprises several matlab s m and fig files which are located in the directory DEMUSEtool programs Program documentation is located in the directory DEMUSEtool documentation To install DEMUSEtool copy both directories to your hard disk e g to c DEMUSEtool directory and set the path in matlab environment to DEMUSEtool programs o x File Edit view Web Window Help 7 New a Current Directory CADEMUSEtoo M zix PE 5 x Open Ctrl O All changes take effect immediately Import Data Save Workspace As MATLAB search path Add Folder Add with Subfolders hove to Top 4 CADEMUSEtoolhprograms Gg CiMatlab6 Sitoolboxlocal 3 CiMatlabs Sitoolbox matlabigeneral Gg UiMatlabs Sitoolboximatlablops Gg CiMatlabs Sitoolboximatlabvand eg CiMatlabb6 Sitoolboximatlabvelmat Gg CiMatlabs Sitoolboximatlableltun eg CiMatlab6 Sitoolbox matlab spectun eg CiMatlabs Sitoolboximatlabimatfun m
20. e patterns are plotted by a click ona Plot MU discharges button Figure 30 A matlab figure opens with a plot of all reconstructed MU discharge patterns Figure 31 Each circle in the figure corresponds to a single MU discharge The horizontal position of the circle denotes the time of MU discharge whereas its vertical position reflects instantaneous MU discharge rate calculated as a quotient between the sampling frequency and the inter pulse interval preceding the given MU discharge Discharge patterns of different MU are depicted one above the other Plot MU discharges F Figure 30 Plot MU discharge button i 2 File Edit View Insert Tools Window Help lDB HSlKAAT SE2 oi i Q E A instantaneous discharge rate pps w 5 Time s Figure 31 Plots of reconstructed MU discharge patterns Vertical axis on the left displays MU IDs vertical axis on the right displays the instantaneous discharge rates in pulses per second pps The tick lines on the right denote the discharge rates of 5 10 and 15 pps respectively Discharge patterns of different MUs are depicted one above the other Page 20 of 34 DEMUSE tool User Manual DEMUSE g 2 Instantaneous discharge rate plots Ta IN OQ MU instantaneous discharge rate plots are plotted by a click on a Plot discharge rates button Figure 32 A matlab figure opens with a different colour lines depicting the ins
21. e slider Decomp runs Figure 12 As a general rule the number of iterations should be larger or equal to the number of expected MUs excluding the small and deep MUs which contribute the background noise only As the exact number of MUs is difficult to estimate the number of decomposition runs should be large default value is set to 30 DEMUSEtool automatically tests the reconstructed MU discharge patterns against the predefined ranges of physiological variables i e discharge rate variability of inter pulse interval etc and discards all the outliers As a result from 5 to 15 most reliably reconstructed MUs are taken into consideration Reconstructed MUs are additionally sorted with respect to the aforementioned degree of decomposition reliability the first MU being the most reliable one Figure 12 Slider for selection of decomposition runs For practical reasons the user can also specify the initial signal offset Figure 13 This allows him to discard the initial signal portions where for example the contraction level is not yet stabilized etc The initial signal offset is measured in seconds Figure 13 Slider for selection of initial SEMG signal offset The decomposition starts by clicking on Run decomposition button Figure 14 ae Run decomposition Processing time Figure 14 Run decomposition button Page 11 of 34 DEMUSE tool User Manual DEMUSE During the decompositi
22. e v1 32 2 It runs on a standard PC and enables the user to e load and visualize the multichannel surface electromyograms e decompose the sEMG signals into contributions of individual motor units MUs e inspect and edit the results attained by automatic decomposition e display graphs of decomposition results including plots of the MU discharge patterns instantaneous discharge rate motor unit action potentials MUAPs and their 2D and 3D animations e compare the original SEMG signals to the reconstructed MUAP trains e save and reload the decomposition results All the graphs are displayed as regular matlab figures and can be freely manipulated by standard matlab graphic tools i e figure resizing Zooming rotating printing etc User is referred to matlab documentation for further details on the use of matlab graphic user interface For further information on EMG acquisition tool see 2 Note The current version of the DEMUSEtoo l supports decomposition of isometric sEMG f signals only i e the signals acquired during isometric muscle contraction Intensive 8 work on decomposition of dynamic sEMG signals is currently in progress Support for dynamic conditions will be built in the future versions of the DEMUSEtool 2 DEMUSEtool components DEMUSEtool represents the third layer in three tier system architecture Figure 1 The first two layers comprise the 128 channel EMG USB electromyographic signal amplifier 4 and t
23. ehabilitation LISIN 2007 Politecnico di Torino Italy i L DEMUSE MARIE CURIE ACTIONS This work was supported by a Marie Curie Intra European Fellowship within the 6th European Community Framework Programme DE MUSE Contract No 023537 Page 33 of 34 DEMUSE tool User Manual DEMUSE References 1 Matlab the language of technical computing MathWorks Inc web address http matworks com 2 Acquisition Software User Manual v1 62 OT Bioelettronica SIRIO Automazione srl and LISIN Bioengineering Center Politecnico di Torino February 2007 3 A Holobar D Zazula Gradient Convolution Kernel Compensation Applied to Surface Electromyograms ICA 2007 LNCS 4666 pp 617 624 2007 4 EMG USB electromyographic signal amplifier User Manual v 1 32 OT Bioelettronica SIRIO Automazione srl and LISiIN Bioengineering Center Politecnico di Torino September 2006 Page 34 of 34
24. f gt pa LJIIAL TMT Figure 34 Plot MU discharge button left and MU selection slider right Figure No 1 File Edit Yiew Insert Tools Window Help DaeHarR AAS PPD 3 lt 2 ie a gt n lt 5 gt f OO OO OO OO OO OO OO OO OO OO OO O 0 25 Time ms oO Figure 35 Multi channel MUAPs estimated by spike triggered averaging of sEMG signals sEMG signals were recorded with a grid of 61 electrodes arranged in 5 columns and 13 rows The location of the innervation zone tendon regions and propagation of motor unit action potentials are visible Page 22 of 34 DEMUSE tool User Manual DEMUSE 4 6 4 2D MUAP map animation DEMUSEtool offers two animations of MUAP generation propagation and attenuation process The first one so called 2D MUAP map is a pseudocolor plot of the estimated MUAP amplitude in a given time instant First MUAP templates are estimated by a spike triggered averaging of the sEMG channels sEMG channels are then spatially organized into a discrete 2D map reflecting the relative position of pick up electrodes The amplitudes of MUAP templates at a given time instant specify the colour on this 2D map of channels The missing intermediate points on the map are calculated by bilinear interpolation of the MUAP amplitudes in four adjacent sSEMG channels In the next animation frame the animation time is advanced by one sample and the 2D MUAP map is recalculated Figur
25. g circle position and length of the time interval currently displayed in the top panel men an E amp 2 z il a S oo x rT a oO vertical scale denotes the instantaneous discharge each circle depicts rate of selected MUs in pulses per second pps single MU discharge Figure 23 Bottom panel of the sEMG editor window with the instantaneous discharge rate plot of a selected MUs Page 15 of 34 DEMUSE tool User Manual DEMUSE lt SEMGeditor m 10 x Fie Edit View Insert Tools Window Help Time s 0 65 add selected MUAPs delete selected MUAPs lt lt move selected MUAPs move selected MUAPs gt gt Electrode configuration C Monopolar Single differential Double differential Laplace Filter 20 SEMG channel e ry oN OAR WON e ee e e ee ee Discharge rate pps CU 9 Figure 24 MU 1 is selected by clicking on its MUAP template in the central panel MUAP occurrences and instantaneous discharge rates of selected MU are automatically depicted in the top and the bottom panel respectively lt SEMGeditor e 10 x File Edit view Insert Tools Window Help Time s 6 C Monopolar Single differential Double differential Laplace Filter 20 500 Hz SEMG channel e ry oN OT wWN e e ee ee ee e ee ee 2 e 2 CICICICICIC ICICI IITs Discharge ra
26. he SEMG acquisition software 2 respectively SEMG signals under investigation are first acquired with a 2D matrix of electrodes The matrix is put on the surface of the skin above the investigated muscle and connected to the EMG USB signal amplifier 4 which amplifies and band pass filters the signals and send them to the sEMG acquisition software 2 The acquisition software acquires the sEMG signals displays them on the screen for immediate visual inspection of the signal quality and saves them into so called SIG files 2 Information about the measurement session including subject data and technical specifications i e number of channels sampling frequency gain etc is saved into purposely designed abstract file 2 Finally both SIG and abstract files are loaded into DEMUSEtool and processed off line Figure 1 Page 3 of 34 DEMUSE tool User Manual DEMUSE Tier 3 DEMUSE decomposition tool EM G system E M G USB R Figure 1 Three tier architecture with indicated data flows of the sEMG decomposition system DEMUSEtool installation DEMUSEtool v2 3 is still a prototype and runs in the matlab programme environment 1 As such it is supported by several personal computer PC platforms including Windows Linux and Mac OS It does not require any special hardware configuration However in the case of large number of channels and long sEMG signals a sufficient amount of RAM 512 MB or more should b
27. on the reconstructed MU discharge patterns are displayed Figure 15 Figure 16 and the instantaneous MU discharge rate and inter pulse interval variability are automatically calculated The reconstructed discharge pattern is put into the set of reconstructed MUs if and only if the calculated values fall within the expected range of values i e discharge rate between 6 and 45 pulses per second Coefficient of Variability of inter pulse interval smaller than 20 The selected MU discharges of accepted discharge pattern are depicted by green circles Figure 15 The reconstructed discharge pattern with the calculated values outside the expected range of values is discarded The selected MU discharge pulses of discarded pattern are depicted by red circles Figure 16 Figure No 1 SS xl Decomp run 12 cp pO gt e a p 2 De aF L SS Laa auaa Iai lMM iiII Iiei IMs gt 7 ro Time sl Figure 15 MU discharge pattern reconstructed in the 12 iteration of the gCKC decomposition technique Each pulse corresponds to a single MU discharge Average discharge rate and coefficient of variability of inter pulse interval are within the expected range of values Thus the selected MU discharges are depicted by green circles lol x Decomp run 13 m N i ary D b di al Vo o D g YY G ey 0 ip o K S p p iT A o cy iit 5 Time sl Figure 16 MU discharge pattern re
28. tantaneous discharge rates of different MUs one line per each MU The thick grey line depicts the exerted muscle force when measured during the acquisition of SEMG signals Instantaneous discharge rates are calculated by low pass filtering of the reconstructed discharge patterns 4 order Butterworth filter with cut off frequency set to 50 Hz s pete tee Figure 32 Plot MU discharge button E lol x File Edit View Insert Tools Window Help DaeWe kK AAS PPO mu WW MU 5 a Q Q Q _ o O O A Ee n 5 Q A E a MUG MUS Figure 33 Plot of instantaneous discharge rates coloured thin lines Vertical axis depicts the instantaneous discharge rates in pulses per second pps Thick grey line depicts the measured muscle force Page 21 of 34 DEMUSE tool User Manual DEMUSE 4 6 3 Multichannel MUAP plots Multichannel MUAP plots so called MU fingerprints can be plotted by clicking on a Plot MUAPs button Figure 34 left panel MU to be depicted is selected by a MU selection slider Figure 34 right panel A matlab figure opens Figure 35 with MUAP shapes as estimated by a spike triggered averaging of each acquired sSEMG channel Displayed MUAPs are spatially organized in rows and columns reflecting the relative position of pick up electrodes in longitudinal single differential configuration Plot MUAPs E MUI aj l o
29. te pps Figure 25 The same as in Figure 24 with all MUs selected Page 16 of 34 DEMUSE tool User Manual DEMUSE 4 5 3 Adding a new MUAP occurrence A double click on a MUAP template in the central panel followed by a click on the add selected MUAP button adds the MUAP of the selected MU to the top panel The optimal MUAP position is automatically determined by the minimum squared error between the residual signal Figure 26 yellow line and the selected MUAP template In the central panel the selected MUAP template is denoted by a red rectangle Figure 26 central panel A further double click on the already selected MUAP template e g MUAP template of MU 1 in Figure 26 deselects the corresponding MU SEMGeditor File Edt View Insert Too Window Help Time s 0 65 0 7 DAVIDA ARID AD ASV VA add selected MUAPs delete selected MUAPs lt lt move selected MUAPs move selected MUAPs gt gt yo Age pee y Fi al a E a E ai li mra A ne 2 3 4 5 6 7 8 Figure 26 Double click on a MUAP template of MU 1 followed by a click on a add selected MUAP button adds the MUAP to the top panel The MUAP template selected by double click is depicted by a red rectangle Note Double click on a MUAP template could directly result in adding of the new MUAP occurrence to the top panel However this would cause errors in the case of unintentional double clicks on the central panel Therefore adding
30. ted SEMG channels and corresponding MUAP trains appears Figure 42 Note Due to the large number of acquired sSEMG channels only selected f row column of SEMG channels can be displayed in one figure The number of figures however is not limited You can display the MUAP trains on all the sEMG channels by consecutively selecting the different electrode columns for example Page 28 of 34 DEMUSE tool User Manual DEMUSE s f E PlotMUAPtrains e r o A T e J 6c m a a E 8c CE a D E E E iuc s E E E E E E J EERTE Figure 41 Channels selection frame left and Plot MUAP trains button right ix File Edit View Insert Tools Window Help lO eH SIRAASL BSE2 Row 5 SIR 18 54 Time s Figure 42 Matlab figure of selected SEMG channels and corresponding reconstructed MUAP trains The range of SIRs of the depicted channels is displayed on the top of the figure Page 29 of 34 DEMUSE tool User Manual DEMUSE Plots of reconstructed MUAP trains are displayed as matlab figures and can be freely manipulated by matlab figure editing tools i e figure resizing Zooming rotating printing etc Zoomed in portion of Figure 42 is depicted in Figure 43 The user is referred to matlab documentation for further details on the use of the matlab s graphic user interface lt Figure No 1 a File Edit view Insert Tools Window Help DOaeHaG rR AAS PED Row 6 SIR
31. terference etc All together these sources add to the measurement noise and affect the efficiency of the SEMG decomposition By comparing the sum of the reconstructed MUAP trains to the original SEMG signal one can estimate the relative proportion of recognized SEMG components and hence speculate on the power and even on the nature of the measurement noise In the DEMUSEtool the reconstructed MUAP trains are calculated as follow Firstly the MUAP shapes are estimated by spike triggered averaging of the acquired SEMG channel using the identified MU discharge instants as trigger The estimated MUAP shapes are then convolved with the identified MU discharge patterns and summed together The sum of MUAP trains is subtracted from the original SEMG signals and the following signal to interference ratio SIR between the original SEMG signals and the residue after the subtraction is calculated E C n dz o SIR i 1 2 a E a n 100 where x n denotes the i th sEMG measurement z 7 stands for the j th MU s MUAP train reconstructed from the i th sEMG measurement and E stands for sample mean Finally the range of SIRs of depicted sEMG channels is displayed together with the reconstructed MUAP trains Figure 42 To display reconstructed MUAP trains select the corresponding electrode row or electrode column Figure 41 left panel and click on Plot MUAP trains button Figure 41 right panel Matlab figure with selec
Download Pdf Manuals
Related Search