0% ont trouvé ce document utile (0 vote)
126 vues63 pages

Traitement des Signaux Biomédicaux

Transféré par

Belalia
Copyright
© © All Rights Reserved
Nous prenons très au sérieux les droits relatifs au contenu. Si vous pensez qu’il s’agit de votre contenu, signalez une atteinte au droit d’auteur ici.
Formats disponibles
Téléchargez aux formats PDF, TXT ou lisez en ligne sur Scribd
0% ont trouvé ce document utile (0 vote)
126 vues63 pages

Traitement des Signaux Biomédicaux

Transféré par

Belalia
Copyright
© © All Rights Reserved
Nous prenons très au sérieux les droits relatifs au contenu. Si vous pensez qu’il s’agit de votre contenu, signalez une atteinte au droit d’auteur ici.
Formats disponibles
Téléchargez aux formats PDF, TXT ou lisez en ligne sur Scribd

SYS867 SUJETS SPÉCIAUX I:

TRAITEMENT DES SIGNAUX BIOMÉDICAUX

Mesures du mouvement humain


Classification descriptive des
habiletés motrices
• Mouvement discret
– Début et fin observables
– Durée variable
– Exemples: lancer, signer son nom

• Mouvement continu
– Début et fin non-définis, le mouvement s’arrête à un moment
arbitraire
– Durée habituellement plus longue (vs. mouv. discrets)
– Poursuite, poursuite compensatoire, poursuite par saut
– Exemple: Nager, courir, conduire une automobile

• Mouvement sériel
– Discret ou continu.
– Série de mouvement assemblés ensemble pour faire un tout
– L’ordre des mouvements est important
– Exemple: allumer un feu, jouer du piano
Classification basée sur la
prédictibilité de l’environnement
• Habileté ouverte
– Environnement changeant
– Planification entière du mouvement impossible
– Succès de la tâche dépend de la capacité d’adaptation.
– Exemple: réaction du gardien de but lors d’une échappée

• Habileté fermée
– Environnement presque parfaitement stable, prévisible
– Exemple: tir à l’arc
– Environnement peut être changeant mais les changements
sont prévisibles (dû à la pratique par exemple)
– Exemple: jongler
Classification des
habiletés motrices: continuums
Mouvement Mouvement Mouvement
discret sériel continu

Début et fin Assemblage de Début et fin


observables sous-mouvements non-définis

Habileté Habileté
fermée ouverte

Environnement Environnement Environnement


prévisible semi-prévisible imprévisible
(conduire une automobile)
Objectivité
• Une mesure est dite objective si deux évaluateurs
indépendants en arrivent au même résultat.

• Importance de l’objectivité dans l’étude du mouvement? Dans


le domaine de la santé? En laboratoire?

• Importance de la sensibilité (précision) de la mesure


– Exemple: saut en hauteur
Fiabilité
• Mesure reproductible, répétable

• Causes pour un manque de fiabilité


– Variabilité dans le matériel de mesure (ruban à mesurer qui s’étire…)
– Variabilité intra-sujet (plusieurs facteurs modifie la performance: fatigue,
motivation, médication,…). Exemple de l’UPDRS.

• Moyens pour augmenter la fiabilité


– Mesures en laboratoire (toutefois, ces mesures sont moins écologiques)
– Utiliser les mêmes instruments de mesure
– Éliminer toutes les sources de distraction. Environnement contrôlé
– Répéter exactement les mêmes procédures

• Malgré tout, les systèmes biologiques sont de nature instable. Pour réduire
la variabilité, il faut rendre plusieurs mesures du même mouvement et en
faire une moyenne
Validité
• Est-ce que le test mesure vraiment ce qu’il dit mesurer?

• Exemple: javelot
– Est-ce que la force brut d’un sujet permet d’évaluer l’habileté
d’un sujet à lancer du javelot?
– Vs. distance parcourue du javelot (face-valid test)
Mesures du mouvement: trois
niveaux d’analyse
• Niveau I: résultats du mouvement (en fonction du but). Est-ce que
la cible est atteinte?

• Niveau II: quantification de l’exécution du mouvement (cinétique,


cinématique, biomécanique)

• Niveau III: analyse de l’activité neuronale avant et pendant le


mouvement.
Mesures de précision:
• Sensibilité de la mesure: dichotomique (réussi vs. non réussi) →
infini.

• Type de précision
– Cinétique (p/r à une force)
– Cinématique (p/r à une vitesse)
– Spatiale (p/r à un emplacement)
– Temporel (p/r à un temps de mouvement)
Mesures de précision:
Erreur constante (EC)
• Erreur moyenne de la précision d’un participant

EC= Σ(xi - T)/n

xi : position finale
T : cible

• Interprétation:
– Représentation de la position finale d’un mouvement par
rapport à une cible
– Indique la déviation (biais) du mouvement. i.e. la direction de
l’erreur
– N’indique pas la variabilité ou la dispersion des réponses (par
exemple: -8, -15, 8, 15 = 0).
Mesures de précision:
Erreur cartésienne vs. vectorielle
Erreur cartésienne Erreur vectorielle

Target: X1, Y1 X1 Target: X1, Y1

Dc
V2-V1
Ec

V1
Y1
Endpoint: X2, Y2
Endpoint: X2, Y2

Y2 Y2
V2
θ

Starting base (0,0) Starting base (0,0)


X2 X2

Ec: Cartesian extent error (Y2-Y1) Ev: Vectorial extent error: V2-V1
Dc: Cartesian direction error: (X2-Y1) Dv: Vectorial direction error: θ: tan –1 X2/Y2
Mesures de précision:
Erreur variable (EV)
• Mesure de l’inconstance d’un participant. Variabilité d’un
participant. Variabilité en fonction de sa moyenne.

EV= Σ (xi – M)2 / n

xi : position finale
M: moyenne

• Interprétation:
– Différence en fonction de la moyenne et non de la cible. La précision
des réponses n’est pas importante.
– Reflète l’étalement des réponses autour de la moyenne du sujet.
– Utilisée pour évaluer la qualité de la mémorisation de la position
d’une cible.
Mesures de précision:
Ellipse
Mesures de précision:
Ellipse

Messier & Kalaska (1997)


Erreur constante vs. Erreur variable
A B
Mesures de précision:
Variabilité totale (E)
• Appelée aussi RMSE (root-mean square error)

• Mesure d’erreur totale. Mesure globale du succès de la


réponse.

E2 = EC2 + EV2
Ou
E= Σ (xi – T)2 / n

• Interprétation
– Combinaison exacte de EC et EV
– Représentation de la dispersion des réponses autour de la cible
(plutôt que de la moyenne; EV)
– Lorsque les réponses sont près de la cible, E=EV.
Mesures de précision:
Erreur absolue (EA)
• Mesure globale de précision. Sans égard à la direction de
l’erreur

EA= Σ |xi – T | / n

• Interprétation
– Similaire à E (numérateur similaire, élimination du signe)
– Contribution de EC et EV inconnue.
– E est maintenant préférée à EA
Mesures de précision:
Erreur constante absolue (ECA)
• Valeur absolue de EC

ECA= |EC |

• Le signe est enlevé après avoir fait la moyenne des erreurs.

• Interprétation
– Utilisée lorsque l’erreur constante est biaisée par la
distribution des réponses (la moitié des participants a des
valeurs positives, l’autre moitié a des valeurs négatives)
– Mesure indépendante du biais.
Exercice 1
• Calculer E, CE, VE, AE, |EC | pour les participants 1 et 2.

• Position de la cible: 100

• Participant 1: 80, 90, 85, 82, 87

• Participant 2: 95, 105, 100, 97, 102


Exercice 1
Participant 1: 80, 90, 85, 82, 87

Participant 2 (+15): 95, 105, 100, 97, 102

E EC EV EA |EC |

Participant 1 15.61 -15.2 3.54 15.2 15.2

Participant 2 3.55 -0.2 3.54 3.0 0.2

• Importance de présenter plusieurs mesures d’erreur:


-Quand le biais est grand, E représente surtout le biais.
-Quand le biais est faible, E représente surtout la variabilité.
Exercice 1
Table de calcul (participant 1)
Essais xi xi - T xi - M (xi – M)2 (xi-T)2 |xi - T|

1 80 -20 -4.8 23.04 400 20


2 90 -10 -5.2 27.04 100 10
3 85 -15 0.2 0.04 225 15
4 82 -18 -2.8 7.84 324 18
5 87 -13 2.2 4.84 169 13
Somme 424 -76 62.8 1218 76
Moyenne 84.8 -15.2 12.56 243.6 15.2
Racine carrée 3.54 15.6

EC. |EC| = 15.2 EV E EA


Exercice 1
Table de calcul (participant 2)
Essais xi xi - T xi - M (xi – M)2 (xi-T)2 |xi - T|

1 95 -5 -4.8 23.04 25 5
2 105 5 5.2 27.04 25 5
3 100 0 0.2 0.04 0 0
4 97 -3 -2.8 7.84 9 3
5 102 2 2.2 4.84 4 2
Somme 413 -1 62.8 63 15
Moyenne 99.8 -0.2 12.56 12.6 3
Racine carrée 3.54 3.55

EC. |EC| = 0.2 EV E EA


Mesures temporelles:
Temps de réaction

• Temps de réaction (TR): temps minimum nécessaire pour


déclencher une activité motrice en réponse à une stimulation
sensorielle.

• Akinésie: retard dans l’initiation d’un mouvement. Observable


chez des patients atteints de la maladie de Parkinson et
d’Huntington.
Mesures temporelles:
Temps de réaction
Signal Début de Fin de la
d’avertissement Stimulus la réponse réponse

Période préparatoire Temps de réponse

Temps de réaction Temps de mouvement

TR TR
prémoteur moteur
(40-80 ms)

EMG
Mesures temporelles:
Temps de réaction

• Comment éviter l’anticipation?


– période préparatoire variable (entre 1s et 5s.)
– « catch trials »: essais pour lesquels le stimulus n’est pas
présenté (environ 15% des essais)
– réponse variable en fonction du stimulus (TRc): évite que le
participant ne planifie à l’avance son mouvement. Implique que
la direction ET l’amplitude varient d’un essai à l’autre.
Facteurs affectant le TR
• Nombre d’alternatives S-R
• Compatibilité spatiale de la relation S-R
• Demande du mouvement
• Intention
• Pas ou peu de pratique
• Irrégularité de la période préparatoire
• Motivation, éveil, fatigue, attention Éveil
• Prédictibilité (répéter le même mouvement)
• Saillance du stimulus
• Sexe
• Âge
• Dominance manuelle
• Consommation d’alcool / drogue / médication (stimulants)
• Cycle de respiration
• Tremblements des doigts
• Vision centrale vs. périphérique
• Maladie
Mesures temporelles:
Temps de mouvement
• Temps de mouvement (TM): temps séparant l’initiation (fin du TR)
et la fin du mouvement.

• TM ≠ TR. Mesurent différents processus. Important de les


mesurer de façon séparée.

• TR + TM = Temps de réponse (par exemple, le « brake reaction


time »)

• TM ≠ Vitesse

• Importance des compromis vitesse-précision lors de


l’interprétation de TM.

• Bradykinésie: lenteur des mouvements volontaires.


Mesures temporelles:
Temps de mouvement
• Comment déterminer le début et la fin d’un mouvement?
– Interrupteur
– Analyse du signal de vitesse (algorithme de Teasdale et al.
(1993))
– Fin du mouvement: Atteinte de la cible: avec ou sans
immobilisation
• Exemples d’immobilisation
– Vitesse inférieure à 3 cm/s pendant 1000 ms
– Curseur demeure à l’intérieur d’une zone de 2 mm
pendant 500 ms
Algorithme de Teasdale et al. (1993)
1. Filtrer les données.

2. Déterminer le premier échantillon


supérieur à 10% de la vitesse
maximale (Vmax/10). Vmax

3. De ce point, reculer jusqu’au premier


échantillon (S) inférieur ou égal à
(Vmax/10)-(Vmax/100)

4. Calculer l’écart-type entre le premier


échantillon et S (SD).

5. Reculer à partir de S et trouver le


premier échantillon inférieur ou égale à Vmax/10
S-SD.
SD S
6. Ce dernier point correspond au début
du mouvement.
Début
7. Le même algorithme est utilisé pour
déterminer la fin du mouvement.
Mesures cinématiques

Novak et al., 2000


Mesures cinématiques: vitesse
• Dérivé du déplacement (taux de changement du
déplacement).

• Profile en forme de cloche, habituellement symétrique


(mouv. simple)

• Mesures de vitesse
– Vitesse moyenne et vitesse maximale
– Moment d’atteinte de la vitesse maximale
– Nombre de pics de vitesse
– Durée précédent et suivant l’atteinte de la vitesse
maximale
Mesures cinématiques: vitesse
Méthodes de dérivation
A. Classique: (x i+1) – (x i) /∆t
B. « Central finite technique »: (x i+1) – (x i-1) /2∆t

P A B
22 --- ---
24 0.1 0.18

29 0.25 0.38
Taux d’échantillonnage: 50 hz
(20 ms entre chaque mesure)
39 0.5 0.58

52 0.65 0.73

68 0.8 ---
Mesures cinématiques: vitesse

Desmurget et al., 2004


Mesures cinématiques:
accélération
• Dérivée de la vitesse (taux de changement de la vitesse). Peut
aussi être mesurée directement (accéléromètre).

• Mesures d’accélération
– Accélération minimale, maximale
– Nombre de « zero-crossing », de transitions accélération-
décélération
– Durée des impulsions

• Le profil d’accélération permet de segmenter le mouvement en


sous-mouvements (nécessaire au calcul des coefficients de
variation).
Correction du mouvement
• L’atteinte d’une cible requière souvent
Cible
d’effectuer une correction (sous-
mouvement, deuxième impulsion). Vitesse

• Procédure. Une correction est présente


Distance
si l’une de ces 3 situations est
observée:
– Si le profil de vélocité passe du
positif au négatif (indique un
changement de direction)
– Si le profil d’accélération traverse la
ligne “zéro” (zero-crossing) une
deuxième fois (indique un
allongement du mouvement).
– Si une discontinuité dans le profil
de décélération entraîne un
changement de signe au niveau du
profil de jerk.
Coefficient de variation
CV: EV/[emplacement du marqueur – emplacement de la position de départ]

• Ce coefficient est habituellement calculé à différents marqueurs tels que:


– Au pic d’accélération
– Au pic de vitesse
– Au pic de décélération
– À la fin de la première impulsion
– À la fin du mouvement

• Interprétation
– Si l’exécution du mouvement reflète uniquement le processus de
planification, la variabilité devrait augmenter de façon linéaire jusqu’à
l’atteinte de la cible et le coefficient devrait demeurer stable d’un
marqueur à l’autre.
– Si le mouvement est corrigé on-line, une diminution du coefficient de
variation sera observée.
Coefficient de variation
0.36 Context-context
No Context- Context
0.34
Context- No Context
0.32 No Context- No Context

0.30
0.28
Coefficient of variation

0.26
0.24
0.22
0.20
0.18
0.16
0.14
0.12
0.10
Peak velocity Peak deceleration End 1st impulse End

Kinematic landmark
Figure 3
Mesures cinématiques: jerk

• Dérivée de l’accélération

• Utilisée comme mesure de l’irrégularité, de la « douceur », la


« fluidité » du mouvement. Permet de déterminer l'optimalité d'un
mouvement et de quantifier le degré d'irrégularité dans des
pathologies motrices.
– Indique les « secousses »
– Irrégularités dues à des commandes correctives ou des
blocages biomécaniques (par exemple, frapper un obstacle,
dispositif de pointage en butée).
Loi du jerk minimal
(coût cinématique)

Prémisse: l’exécution d’un mouvement optimal requière de


minimiser le jerk, de minimiser les « coûts »

Jerk minimal:

Le jerk minimisé (coût) est déterminé pour chaque axe (x, y,


z) de façon indépendante
Loi du jerk minimal
Mouvement rapide 2D, sujet en santé

CC==9665.25
9665.25cm
cm3/s
3 33
/s

Fimbel, 2005
Loi du jerk minimal
Mouvement rapide 2D, sujet parkinsonien

CC==11575,50
11575,50cm
cm3/s
3 33
/s

Fimbel, 2005
Mesures cinématiques:
Jerk
• Les mesures globales de jerk (par exemple, jerk minimal) sont
associées au degré d’optimalité du mouvement.

• Contre-exemple: Patients en réadaptation (ACV)

Début de la pratique Fin de la pratique

2 sous-mouvements séparés Les 2 sous-mouvements s’unissent


par une période de latence Mesure globale de jerk: élevé.
Mesure globale de jerk: faible.

Roehrer & al. 2003


Mesures cinématiques:
Jerk normalisé
• Jerk normalisé: sqrt(1/2 · j2(t) dt · d5/l2)
j: 3ième dérivé du déplacement
d: durée
l: amplitude du mouvement

• Procédure
– Filtrer (filtre lowpass)
– Estimer la vélocité [cm/s] et le jerk [cm/s**3]
– Segmenter en sous-mouvements (d’une vélocité nulle à une vélocité nulle)
– Estimer l’amplitude du sous-mouvement [cm]
– Estimer la durée du sous-mouvement [s]
– Mettre les données du jerk au carré [cm**2/s**6]
– Intégrer par rapport au temps [cm**2/s**5]
– Diviser par la distance **2
– Multiplier par la durée **5
– Racine carrée et diviser par 2
Mesures cinématiques:
Jerk normalisé
• Normalisé pour l’amplitude et la durée du mouvement

• S’applique au mouvement entier (ne peut pas être calculé sur les
corrections seulement).

• Jerk sans dimension (« unit free ») : améliore les possibilités de


comparaisons inter-sujets et inter-mouvements
– Temps d'exécution différent
– Distance parcourue différente

• Jerk normalisé pour mouvement optimal (sans secousse)


– Mouvement rectiligne: 7.75 (pi**3/2**2)
– Cercle: 10.96
Mesures cinématiques:
Jerk normalisé

Teulings et al., 1997


Irrégularité et
maladie d’Huntington

Smith et al. (2000)


Irrégularité et
maladie d’Huntington

• Mouvements de pointage
rapides sont plus irréguliers
chez les patients MH, mais
seulement en fin de
mouvement.

• Cette irrégularité serait


associée à un problème de
traitement du feedback
pendant la phase de
correction du mouvement.

Smith et al. (2000)


Irrégularité et
maladie d’Huntington
• Une impulsion de 6, 12 ou
18 N était appliqué par un
bras haptique.

• Les perturbations causent


chez les patients:
– Une augmentation de la
longueur de tracé.
– Une augmentation du
niveau de jerk.

Smith et al. (2000)


Autres mesures d’irrégularité
• Snap: dérivé du jerk (quatrième dérivée du déplacement).

• Compter le nombre de pics de vitesse ou de transitions accélération-


décélération.

• Compter le nombre de singularités de jerk (par exemple, inversion du


jerk). Approche qualitative. Voir Fimbel et al. (2003)

Fimbel, 2005
Trajectoire
• Index de linéarité (Atkeson & Cible
Hollerbach, 1985)*. Représente le
degré de courbure de la trajectoire

IL: d / a d
d: distance maximale d’une ligne
séparant un vecteur joignant la
base de départ et la cible et la
a
trajectoire
a: amplitude du mouvement

• L’exécution d’un mouvement avec


le doigts entraîne une trajectoire
plus courbée que les mouvements Base de départ
effectués avec un objet (stylet,
souris)

*Voir aussi Contreras-Vidal, J.L., Teulings, H.L., & Stelmach, G.E. (1998).
Variabilité de la trajectoire

• La variabilité de la trajectoire initiale est utilisée afin d’évaluer


l’efficacité de la phase de planification.

• Procédure:
– L’erreur constante angulaire est
mesurée pour tous les essais vers une
même cible xi

– EV= Σ (xi – M)2 / n


où xi: = angle
M= moyenne
Préhension: atteinte de l’objet
2 phases interdépendantes
• Transport (atteinte)
– Contrôlé par des muscles proximaux
• Manipulation (ouverture des mains)
– Contrôlé par des muscles distaux

Jeannerod (1981, 1984)


Préhension: atteinte de l’objet
Préhension: atteinte de l’objet
• Variables d’intérêts:
– Phase de transport (habituellement mesurée à l’aide d’un
marqueur placé sur le poignet)
• temps de mouvement
• vitesse maximale
• durée (absolue et relative) de la phase d’accélération/décélération

– Phase de manipulation (habituellement mesurée à l’aide de


marqueurs placés sur le bout de l’index et du pouce)
• ouverture maximale des doigts,
• durée (absolue et relative) précédent et suivant l’ouverture
maximale des doigts.

– Corrélation entre les phases de transport et de préhension.


Préhension: atteinte de l’objet.
Rôle de la vision

Churchill et al., 2000


Préhension: saisie de l’objet
Mesure de force (tâches isométriques,
saisie d’un objet)
- jauge de contrainte (transducers)
- cybergloves(R).

M= mass
LF: force de soulèvement
GF: force de préhension
ACC: accélération
G: gravité
Préhension : saisie de l’objet
Lorsqu’un objet est soulevé:
• La force de préhension est légèrement supérieure au minimum nécessaire
pour éviter que l’objet glisse (marge de sécurité).
• La force de soulèvement augmente de façon proportionnelle à la force de
préhension. Ce ratio (Fs vs. Fp) constitue la variable de contrôle.
• L’absence de délai entre le soulèvement d’un objet et la modification de la
force de préhension suggère l’utilisation de modèle interne (prédictif).
Préhension : saisie de l’objet

L’objet est immobile Force minimale nécessaire pour éviter


d’échapper l’objet
L’objet se déplace verticalement

Blennerhassett et al., 2006


Conclusions
• Mesures de précision
– Erreur constante (cartésienne, vectorielle)
– Variabilité de la position finale (erreur variable, ellipse)
– Coefficient de variation

• Mesures temporelles
– Temps de réaction et de mouvement

• Mesures cinématiques et d’irrégularité


– Vitesse
– Accélération
– Jerk

• Évaluation de la trajectoire
– Index de linéarité
– Variabilité initiale

• Préhension et mesures isométriques


Références
Index de linéarité:
Atkeson CG, Hollerbach JM (1985) Kinematic features of unrestrained
vertical arm movements. J Neurosci 5: 2318-2330
Contreras-Vidal, J.L., Teulings, H.L., & Stelmach, G.E. (1998). Elderly
subjects are impaired in spatial coordination in fine motor control. Acta
Psychologica, 100, 25-35.

Vitesse:
M. Desmurget, V. Gaveau, P. Vindras, R. S. Turner, E. Broussolle and S.
Thobois (2004). On-line motor control in patients with Parkinson's
disease. Brain,127 (8), 1755-1773

Cinématique:
K.E. Novak, L.E. Miller, J.C. Houk (2000). Kinematic properties of rapid
hand movements in a knob turning task. Experimental Brain Research,
132 (4), 419 - 433
Références
Début et fin de mouvement:
Teasdale N, Bard C, Fleury M, Young DE, Proteau L. (1993). Determining
movement onsets from temporal series. J Mot Behav. 25(2), 97-106.

Jerk:
Fimbel E, Dubarry AS, Philibert M, Beuter A. (2003). Event identification in
movement recordings by means of qualitative patterns. Neuroinformatics.
1(3):239-57.
Teulings HL, Contreras-Vidal JL, Stelmach GE, Adler CH. (1997).
Parkinsonism reduces coordination of fingers, wrist, and arm in fine motor
control. Exp Neurol. 146(1),159-70.

Préhension:
Churchill, A., Hopkins, B., Ronnqvist, L., Vogt, S. (2000). Vision of the hand
and environmental context in human prehension. Exp. Brain Res, 134,
81-89
Liens
Début et fin de mouvement :
• [Link]
04/[Link]

Pour tester votre TR


• [Link]
Liens
Instruments de mesure:
• [Link]
• [Link]
• [Link]
• [Link]
[Link]/mpbt/haptics/hand_controllers/freedom/freedom.
html
• [Link]
• [Link]
• [Link]

Vous aimerez peut-être aussi