Brevet FR3115124A1 — publié le 20220415

La présente invention concerne une méthode de calibration en ligne d’une interface neuronale directe mettant en œuvre une régression REW-NPLS entre un tenseur de calibration de sortie et un tenseur de calibration d’entrée. La régression REW-NPLS comprend une décomposition itérative PARAFAC du tenseur de covariance croisée entre le tenseur de calibration d’entrée et le tenseur de calibration de sortie, chaque itération PARAFAC comprenant une séquence de étapes élémentaires (240 1 , 240 1 ,.. 240 M ) de minimisation d’une métrique selon la méthode des moindres carrés alternés, chaque étape élémentaire de minimisation portant sur un projecteur et considérant les autres comme constants, ladite métrique comportant un terme de pénalisation fonction de la norme de ce projecteur, les éléments de ce projecteur ne faisant pas l’objet d’une pénalisation lors d’une itération PARAFAC n’étant pas pénalisables lors des itérations PARAFAC suivantes. Ladite méthode de calibration permet d’obtenir un modèle prédictif dont les coefficients non nuls sont rares par blocs. Fig. 2 MÉTHODE DE CALIBRATION EN LIGNE ET À FACTEUR D’OUBLI D’UNE INTERFACE NEURONALE DIRECTE À RÉGRESSION MULTIVARIÉE PÉNALISÉE La présente invention concerne le domaine des interfaces neuronales directes encore dénommées BCI ( Brain Computer Interface ) ou BMI ( Brain Machine Interface ). Elle trouve notamment à s’appliquer à la commande neuronale directe d’une machine, telle qu’un exosquelette ou un ordinateur. ETAT DE LA TECHNIQUE ANTÉRIEURE Les interfaces neuronales directes utilisent les signaux électro-physiologiques émis par le cortex cérébral pour élaborer un signal de commande. Ces interfaces neuronales ont fait l’objet de nombreuses recherches notamment dans le but de restaurer une fonction motrice à un sujet paraplégique ou tétraplégique à l’aide d’une prothèse ou d’une orthèse motorisée. Les interfaces neuronales peuvent être de nature invasive ou non invasive. Les interfaces neuronales invasives utilisent des électrodes intracorticales (c’est-à-dire implantées dans le cortex) ou des électrodes corticales (disposées à la surface du cortex) recueillant dans ce dernier cas des signaux électrocorticographiques (ECoG). Les interfaces neuronales non invasives utilisent des électrodes placées sur le cuir chevelu pour recueillir des signaux électroencéphalographiques (EEG). D’autres types de capteurs ont été également envisagés comme des capteurs magnétiques mesurant les champs magnétiques induits par l'activité électrique des neurones du cerveau. On parle alors de signaux magnétoencéphalographiques (MEG). Les interfaces neuronales directes utilisent avantageusement des signaux de type ECoG, présentant l’avantage d’un bon compromis entre biocompatibilité (matrice d’électrodes implantées à la surface du cortex) et qualité des signaux recueillis. Les signaux ECoG ainsi mesurés doivent être traités afin d’estimer la trajectoire du mouvement désiré par le sujet et en déduire les signaux de commande de l’ordinateur ou de la machine. Par exemple, lorsqu’il s’agit de commander un exosquelette, l’interface BCI estime la trajectoire du mouvement désiré à partir des signaux électro-physiologiques mesurés et en déduit les signaux de contrôle permettant à l’exosquelette de reproduire la trajectoire en question. De manière similaire, lorsqu’il s’agit de commander un ordinateur, l’interface BCI estime par exemple la trajectoire souhaitée d’un pointeur ou d’un curseur à partir des signaux électro-physiologiques et en déduit les signaux de commande du curseur/pointeur. L’estimation de trajectoire, et plus précisément celle des paramètres cinématiques (position, vitesse, accélération), est encore dénommée décodage neuronal dans la littérature. Le décodage neuronal permet notamment de commander un mouvement (d’une prothèse ou d’un curseur) à partir de signaux ECoG. L’estimation de trajectoire et le calcul des signaux de contrôle de l’exosquelette ou de l’effecteur, nécessite généralement une phase d’apprentissage ou de calibration préalable, dite off-line. Lors de cette phase, le sujet imagine, observe ou réalise un mouvement selon une trajectoire déterminée pendant un intervalle de calibration donné. Les signaux électrophysiologiques mesurés pendant cet intervalle sont exploités en relation avec cette trajectoire pour construire un modèle prédictif et plus précisément pour calculer les paramètres de ce modèle. La validité du modèle prédictif est cependant limitée dans le temps en raison de la non-stationnarité des signaux neuronaux. Pour cette raison, il est nécessaire de procéder à une calibration en ligne ( on-line ) du modèle prédictif, c’est-à-dire au fur et à mesure que les signaux neuronaux sont observés et la commande appliquée. Différentes méthodes de calibration en ligne d’une interface neuronale ont été décrites dans l’état de la technique. Toutefois, ces méthodes de calibration en ligne peuvent nécessiter la manipulation d’une quantité considérable de données résultant de la juxtaposition de données d’étapes de calibration antérieures et celles d’une nouvelle étape de calibration. Une méthode de calibration on-line d’une interface BCI a été décrite dans l’article de A. Eliseyev et al. intitulé « Recursive exponentially weighted N-way Partial Least Squares regression with recursive validation of hyper-parameters in Brain-Computer Interface applications » publié dans Scientific Reports, vol. 7, n° 1, p. 16281, nov. 2017 ainsi que dans la demande de brevet FR-A-3 061 318. Cette méthode sera désignée par la suite sous l’acronyme REW-NPLS ( Recursive Exponentially Weighted N-way Partial Least Squares ). D’autre part, différentes variantes de méthode de calibration off-line ont été proposées dans la littérature, dans le but d’améliorer la robustesse de la prédiction et d’accroître la rareté ( sparsity ) des coefficients non nuls dans le tenseur du modèle prédictif. En particulier, la rareté des coefficients non nuls dans le mode spatial de ce tenseur permet de ne sélectionner que certaines électrodes pertinentes pour la prédiction de la trajectoire et le calcul de la commande. Une telle variante utilisant une régression multivoie pénalisée dans un tel but a été décrite dans l’article de A. Eliseyev et al. intitulé « L1-penalized N-way PLS for subset of electrodes selection in BCI experiments » publié dans J. Neural Eng. vol. 9, n° 4, p. 045010 ainsi que dans la demande de brevet FR-A-3 046 471. Toutefois, une telle méthode de calibration utilisant une régression multivoie pénalisée ne peut être mise en œuvre en ligne en raison de la quantité importante de calculs qu’elle requiert. L’objet de la présente invention est par conséquent de proposer une méthode de calibration en ligne d’une interface neuronale directe qui puisse également accroître la rareté des coefficients non nuls du tenseur de prédiction, notamment de manière à pouvoir ne sélectionner à chaque fenêtre d’observation qu’un sous-ensemble pertinent d’électrodes. PRÉSENTATION DE L’INVENTION La présente invention est définie par une méthode de calibration en ligne d’une interface neuronale directe destinée à recevoir une pluralité de signaux électrophysiologiques pendant une pluralité de fenêtres d’observation associées à des instants d’observation, pour former un tenseur d’entrée et en déduire, au moyen d’un modèle prédictif, un tenseur de sortie fournissant des signaux de commande destinés à commander un ou plusieurs effecteurs, ladite calibration est effectuée lors d’une pluralité d’étapes de calibration, chaque étape de calibration étant réalisée à partir de données de calibration d’entrée, représentées par un tenseur de calibration d’entrée, , et de données de calibration de sortie représentées par un tenseur de calibration de sortie , ladite étape de calibration mettant en œuvre une régression REW-NPLS comportant : le calcul d’un tenseur de covariance à partir d’un tenseur de covariance obtenu lors d’une étape de calibration précédente et du tenseur , ainsi que le calcul d’un tenseur de covariance croisée à partir d’un tenseur de covariance croisée obtenu lors de ladite étape de calibration précédente et des tenseurs et ; une décomposition PARAFAC itérative du tenseur de covariance croisée , chaque itération PARAFAC étant associée à une dimension de l’espace des variables latentes sur lequel est effectuée ladite régression, chaque itération fournissant une pluralité de projecteurs de tailles respectives et un modèle prédictif défini par un tenseur de coefficients de prédiction et un tenseur de biais . Ladite méthode de calibration est originale en ce que chaque itération PARAFAC comprend une séquence de étapes élémentaires de minimisation d’une métrique selon la méthode des moindres carrés alternés, chaque étape élémentaire de minimisation portant sur un projecteur et considérant les autres comme constants, ladite métrique comportant un terme de pénalisation fonction de la norme de ce projecteur, les éléments de ce projecteur ne faisant pas l’objet d’une pénalisation lors d’une itération PARAFAC n’étant pas pénalisables lors des itérations PARAFAC suivantes. Avantageusement, une étape élémentaire de minimisation portant sur un projecteur vise à minimiser la métrique pénalisée où est un projecteur non normalisé et , sont des projecteurs normalisés, est une matrice de taille obtenue par déploiement du tenseur de covariance croisée en prenant le mode comme référence, est le produit de Kronecker, est un coefficient strictement positif, , est une norme définie par , , où est le élément du projecteur , est le symbole de Kronecker, et est l’ensemble des indices correspondant aux éléments pénalisables du projecteur à l’itération PARAFAC . Alternativement, une étape élémentaire de minimisation portant sur un projecteur vise à minimiser les métriques pénalisées , où est le kème vecteur colonne de la matrice de taille obtenue par déploiement du tenseur de covariance croisée en prenant le mode comme référence, est le produit de Kronecker, est un coefficient strictement positif, et Les fonctions , sont définies par : si et sinon si et sinon si et sinon où est l’ensemble des indices correspondant aux éléments pénalisables du projecteur à l’itération PARAFAC . Dans un cas particulier de réalisation, et les éléments résultant de la minimisation sont donnés par : si et sinon où représente la solution de minimisation au sens des moindres carrés : et que le seuil est donné par . Dans un second cas particulier de réalisation, et les éléments résultant de la minimisation sont donnés par : si et si et sinon où et où représente la solution de minimisation au sens des moindres carrés : et le seuil est donné par . Dans un troisième cas particulier de réalisation, et les éléments résultant de la minimisation sont donnés par : si et si et sinon où représente la solution de minimisation au sens des moindres carrés : et le seuil est donné par Typiquement, ladite méthode de calibration en ligne peut s’appliquer à des signaux électrophysiologiques constitués de signaux ECoG. BRÈVE DESCRIPTION DES FIGURES D’autres caractéristiques et avantages de l’invention apparaîtront à la lecture d’un mode de réalisation préférentiel de l’invention, décrit en référence aux figures jointes parmi lesquelles : représente de manière schématique la décomposition d’un tenseur par analyse factorielle parallèle (PARAFAC) ; représente de manière schématique le développement d’un tenseur par rapport à un mode ; représente de manière schématique une décomposition utilisant une méthode des moindres carrés alternés (ALS) ; représente de manière schématique l’ordinogramme d’une méthode de calibration récursive d’une interface neuronale directe à régression multivariée pénalisée, selon un mode de réalisation de la présente invention ; représente de manière schématique le mécanisme de mise à jour du modèle prédictif dans la méthode de calibration de la Fig. 2. On considérera dans la suite une interface neuronale directe à décodage continu. Les signaux électrophysiologiques issus des différentes électrodes sont échantillonnés et rassemblés par blocs de données, chaque bloc correspondant à une fenêtre glissante d’observation de largeur . Chaque fenêtre d’observation est définie par un instant d’observation ( epoch ) auquel démarre la fenêtre en question. Les signaux électrophysiologiques peuvent faire l’objet d’un prétraitement. Ce prétraitement peut notamment comporter une suppression de la moyenne prise sur l’ensemble des électrodes, puis une analyse temps-fréquence est effectuée sur chacune des fenêtres d’observation. L’analyse temps-fréquence est basée sur une décomposition en ondelettes, notamment en ondelettes de Morlet. L’homme du métier comprendra toutefois que d’autres type d’analyse temps-fréquence pourront être envisagés par l’homme du métier. Ces résultats de l’analyse temps-fréquence peuvent en outre faire l’objet d’un lissage fréquentiel ou d’une décimation. Ainsi, à chaque fenêtre d’observation, ou instant d’observation , est associé un tenseur d’ordre 3 de données d’observation, d’où la génération d’un tenseur d’entrée d’ordre 4: le premier mode correspond aux fenêtres d’observation successives, le second mode correspond à l’espace, autrement dit aux capteurs, le troisième mode correspond au temps au sein d’une fenêtre d’observation, autrement dit aux positions des ondelettes, et le quatrième mode correspond à la fréquence, autrement dit au nombre de bandes de fréquence utilisées pour la décomposition en ondelettes sur une fenêtre d’observation. De manière plus générale, le tenseur d’entrée (ou tenseur d’observation) sera d’ordre , le premier mode étant dans tous les cas celui relatif aux instants d’observation ( epochs ). Le tenseur d’entrée (ou tenseur d’observation) est noté et est de dimension . De la même façon, la trajectoire du mouvement imaginé, observé ou réalisé est décrite par un tenseur de sortie (ou tenseur de commande) d’ordre , noté , de dimension , dont le premier mode correspond aux instants successifs auxquels s’appliqueront les commandes (en règle générale, ce premier mode correspond également aux fenêtres d’observation), les autres modes correspondant aux commandes de différents effecteurs ou aux différents degrés de liberté d’un robot multi-axe. Plus précisément, le tenseur de sortie fournit des blocs de données de commande, chacun des blocs permettant de générer les signaux de commande relatifs aux différents effecteurs ou degrés de liberté. Ainsi, on comprendra que la dimension de chaque bloc de données pourra dépendre du cas d’usage envisagé et notamment du nombre de degrés de liberté de l’effecteur. On notera dans la suite le tenseur d’observation à l’instant . Ce tenseur est par conséquent d’ordre et de dimension . Il prend ses valeurs dans un espace où est l’ensemble de réels. De manière similaire, on notera le tenseur de commande à l’instant . Ce tenseur de sortie est d’ordre et de dimension . Il prend ses valeurs dans un espace . La calibration ou, de manière équivalente, le mécanisme de mise à jour du modèle prédictif selon l’algorithme REW-NPLS ( Recursive Exponentially Weigthed N-way Partial Least Squares ) est rappelé ci-après. Si l’on note les blocs de données de calibration ( training data set ) lors de l’itération , l’idée de base de l’algorithme REW-NPLS est de mettre à jour le modèle prédictif à partir des tenseurs de covariance et de covariance croisée obtenus à l’étape de calibration précédente, respectivement notés et . Ces deux tenseurs permettent d’obtenir une représentation condensée du modèle prédictif à l’étape précédente, sans qu’il soit besoin de stocker l’historique des blocs de données de calibration précédents : . Plus précisément, les tenseurs de covariance et de covariance croisée sont mis à jour au moyen des relations de récurrence : (1-1) (1-2) où désigne le produit de tenseur de mode 1 et est un coefficient d’oubli avec . On comprend ainsi qu’une période de calibration courante, réalisée lors de l’itération tient compte d’une période de calibration précédente réalisée lors de l’itération et des tenseurs des données d’observation de la période de calibration courante, l’historique étant pondéré par le coefficient d’oubli avec . Cette opération de mise à jour des tenseurs de covariance est effectuée de manière itérative sur les périodes de calibration. Pour chaque mise à jour, le modèle prédictif est ensuite obtenu de manière itérative au moyen d’une décomposition PARAFAC (généralisation au cas tensoriel d’une décomposition en valeurs singulières) du tenseur de covariance croisée, . On rappelle que, d’une manière générale, une décomposition PARAFAC d’un tenseur d’ordre , soit permet de représenter ce tenseur sous la forme d’une combinaison linéaire de produits extérieurs de vecteurs (tenseurs d’ordre 1) à un tenseur résiduel près : (2) avec , où est le produit extérieur et est un tenseur résiduel et est l’indice de l’itération PARAFAC et est un entier donné. Une décomposition PARFAC d’un tenseur a été représenté schématiquement en Fig. 1A. Dans le cas d’une étape de calibration REW-NPLS, chaque itération PARAFAC ne porte que sur une décomposition élémentaire, c’est-à-dire . Elle permet d’extraire du tenseur un ensemble de vecteurs de projection, encore dénommés projecteurs, , de tailles respectives où le numéro d’itération donne la dimension d’un espace de variables latentes sur lequel est projet le tenseur d’entrée . Afin de simplifier la présentation et sans préjudice de généralisation à un tenseur d’ordre quelconque, nous supposerons dans la suite que l’ordre du tenseur est égal à 3 ( ), c’est-à-dire avec . Par exemple, le premier mode correspond à l’époque, le second à la fréquence et le troisième à l’espace (c’est-à-dire aux électrodes). A l’itération de la décomposition PARAFAC, on recherche les vecteurs de projection qui vérifient : (3) avec et Ce problème de minimisation peut être résolu par une méthode des moindres carrés alternés dans laquelle on effectue séquentiellement une minimisation pour chacun des vecteurs de la décomposition, soit : (4-1) (4-2) (4-3) où est le produit de Kronecker et est une matrice obtenue par déploiement du tenseur en prenant le mode comme mode de référence ( tensor unfolding or flattening along mode ). Par exemple, la matrice est définie par où les sont les colonnes de la matrice . Les opérations de minimisation (4-1), (4-2) et (4-3) sont répétées jusqu’à la convergence. Les solutions de minimisation au sens des moindres carrés (LS) pour chacune de ces opérations sont données par les projecteurs : (5-1) (5-2) (5-3) Chaque itération de la décomposition PARAFAC fournit une pluralité de projecteurs et un tenseur de coefficients de prédiction, ou modèle prédictif de l’interface neuronale directe, défini par un tenseur de coefficients de prédiction et un tenseur de biais , permettant prédire le tenseur de commande à partir du tenseur d’observation. On obtient ainsi une pluralité de modèles, , correspondant à différentes dimensions successives de l’espace des variables latentes. On détermine ensuite, à l’étape de calibration suivante, , parmi , , le modèle prédictif correspondant à l’erreur de prédiction la plus faible où l’on a noté l’indice correspondant à cette plus faible erreur de prédiction. Cette détermination du modèle prédictif à l’étape de calibration suivante est désignée dans la littérature relative sous le terme de validation récursive. La Fig. 1C représente de manière schématique une décomposition utilisant une méthode des moindres carrés alternés (ALS). La décomposition par les moindres carrés alternés a été illustrée pour une étape de calibration REW-NPLS et une itération PARAFAC . A l’étape 110, on estime le vecteur (correspondant au mode époque) à partir de la matrice et des vecteurs et courants. Cette opération peut être considérée comme une projection des vecteurs colonnes de la matrice sur l’espace tensoriel généré par les modes 2 et 3. A l’étape 120, on estime le vecteur (correspondant au mode fréquence ) à partir de la matrice et des vecteurs et courants. Cette opération peut être considérée comme une projection des vecteurs colonnes de la matrice sur l’espace tensoriel généré par les modes 1 et 3. A l’étape 130, on estime le vecteur (correspondant au mode spatial) à partir de la matrice et des vecteurs et courants. Cette opération peut être considérée comme une projection des vecteurs colonnes de la matrice sur l’espace tensoriel généré par les modes 1 et 2. Les étapes 110, 120, 130 sont itérées jusqu’à ce qu’un critère de convergence soit rempli, par exemple lorsque la somme des écarts entre les vecteurs de deux itérations ALS successives est inférieure à une valeur de seuil prédéterminée. L’idée à la base de la présente invention est d’effectuer une pénalisation par tranche à chaque itération de l’algorithme PARAFAC, une pénalisation par tranche se traduisant par une élimination des éléments relatifs à une valeur d’indice dans le tenseur du modèle de prédiction. Les éléments ne faisant pas l’objet d’une pénalisation lors d’une itération ne sont plus pénalisables pour les itérations suivantes, . En revanche les éléments pénalisés lors d’une itération sont pris en compte lors de la décomposition PARAFAC de l’itération suivante. On promeut ainsi une rareté par blocs ( blockwise sparsity ) des éléments non nuls dans le tenseur du modèle de prédiction, ce qui simplifie considérablement le calcul de la commande. Plus précisément, lors d’une itération de la décomposition PARAFAC, la minimisation des écarts quadratiques selon (4-1), (4-2) et (4-3) dans la méthode des moindres carrés alternés est remplacée par une minimisation prenant en compte un terme de pénalisation favorisant cette rareté, à savoir : (6-1) (6-2) (6-3) où , et sont des projecteurs non normalisés et , , sont les projecteurs normalisés correspondants et sont des coefficients strictement positifs et inférieurs à 1. On notera que du fait de la présence des termes de pénalisation dans les expressions (6-1), (6-2) et (6-3), la minimisation ne peut plus être réalisée sur des vecteurs de projection normalisés comme dans la méthode ALS classique. Les normes sont définies par : (7-1) (7-2) (7-3) où est le élément du vecteur de projection , est le symbole de Kronecker, et est l’ensemble des indices correspondant aux éléments pénalisables, étant compris que puisque les éléments non-pénalisés d’un projecteur à une itération ne sont pas pénalisables à l’itération . On comprendra notamment que la norme de l’expression (7-1) donne la somme des éléments non nuls dont les indices appartiennent à . Ainsi par exemple, dans ce cas, la pénalisation sera d’autant plus forte que le nombre d’éléments dont les indices appartiennent à est plus élevé. Afin de simplifier le calcul de minimisation dans les expressions (6-1), (6-2), (6-3), on pourra avantageusement opérer élément par élément, à savoir : , (8-1) , (8-2) , (8-3) où les vecteurs , , sont les vecteurs colonnes des matrices , , obtenues par déploiement du tenseur de covariance croisée selon les modes 1, 2 et 3, respectivement. Les fonctions , s’appliquent à des valeurs scalaires et sont définies par : si et sinon (9-1) si et sinon (9-2) si et sinon (9-3) On peut montrer que le résultat de l’opération de minimisation pénalisée selon (8-1), (8-2), (8-3) pour la norme (c’est-à-dire )donne : si et (10-1) sinon (10-2) où représente la solution de minimisation au sens des moindres carrés (5-1), (5-2), (5-3), c’est-à-dire : (11-1) (11-2) (11-3) et les seuils , sont donnés par : (12-1) (12-2) (12-3) L’opération de minimisation pénalisée revient par conséquent à effectuer un seuillage sur les éléments des projecteurs obtenus à chaque étape de la méthode ALS. De manière similaire, on peut montrer que l’opération de minimisation pénalisée selon (8-1), (8-2), (8-3) pour la norme (c’est-à-dire )a pour solution : si et (13-1) si et (13-2) sinon (13-3) où et les seuils , sont donnés par : (14-1) (14-2) (14-3) Les coefficients , sont définis comme les solutions respectives des cubiques , , avec : (15-1) (15-2) (15-3) Enfin, de manière similaire, on peut montrer que l’opération de minimisation pénalisée selon (8-1), (8-2), (8-3) pour la norme (c’est-à-dire )a pour solution : si et (16-1) si et (16-2) sinon (16-3) où les seuils , sont donnés par : (17-1) (17-2) (17-3) On comprend que, quelle que soit la norme utilisée, la minimalisation pénalisée se réduit avantageusement à un simple seuillage des éléments d’un projecteur à chaque étape de la méthode ALS. Cette minimalisation pénalisée permet de promouvoir une rareté par blocs des coefficients non nuls dans le modèle prédictif ( blockwise sparsity ) et par conséquent de simplifier le calcul. Il est important de noter que celle-ci peut être réalisée en ligne dans le but de réaliser une calibration incrémentale du modèle prédictif. La Fig. 2 représente de manière schématique l’ordinogramme d’une méthode de calibration récursive d’une interface neuronale directe à régression multivariée pénalisée, selon un mode de réalisation de la présente invention. A l’étape 200, on prend en compte de nouvelles données de calibration d’entrée et de sortie, respectivement sous la forme d’un tenseur d’entrée et d’un tenseur de sortie . Les données de calibration d’entrée résultent de la mesure des signaux électrophysiologiques. Les données de calibration de sortie résultent par exemple de la mesure de paramètres cinématiques de la trajectoire ou de signaux de commande d’effecteurs. Les tenseurs d’entrée et de sortie peuvent faire le cas échéant l’objet d’un prétraitement visant à les centrer et les normaliser, de manière connue en soi. A l’étape 210, on détermine parmi l’ensemble des modèles prédictifs , obtenus à l’étape précédente, le modèle prédictif correspondant à l’erreur de prédiction la plus faible. Ce modèle est utilisé pour calculer le tenseur de commande (tenseur de sortie) à partir du tenseur d’observation (tenseur d’entrée). A l’étape 220, on met à jour les tenseurs de covariance et de covariance croisée, , selon les relations (1-1) et (1-2), en tenant compte du facteur d’oubli. On entre ensuite dans une première boucle itérative NPLS, dans laquelle on itère sur la dimension de l’espace des variables latentes, avec où est une dimension maximale. Cette boucle itérative permet de générer successivement une pluralité de modèles prédictifs de l’interface neuronale directe, chaque modèle prédictif étant défini par un tenseur de coefficients de prédiction et un tenseur de biais . A l’étape 230, on initialise avec et les ensembles , avec pour . Pour chaque valeur de , on effectue une décomposition PARAFAC du tenseur de covariance croisée à l’aide d’une méthode des moindres carrés alternés. A la différence d’une décomposition PARAFAC conventionnelle (cf. expressions (4-1), (4-2), (4-3)) utilisée dans la méthode REW-NPLS, on utilise ici une décomposition PARAFAC régularisée au moyen d’une métrique pénalisée, le terme de pénalisation visant à promouvoir une rareté par blocs des coefficients non nuls au sein du tenseur représentatif du modèle prédictif. Plus précisément, en 240, pour chaque valeur de , on recherche le -uplet de projecteurs , de dimensions respectives permettant de minimiser ladite métrique pénalisée, comme exprimé en (6-1),(6-2),(6-3) dans le cas . Cette minimisation se fait séquentiellement par ALS, projecteur par projecteur en 240 1 ,…,240 M , les autres projecteurs étant considérés comme constants. Avantageusement, à chacune de ces étapes, la minimisation porte sur les éléments du projecteur en question, comme décrit en relation avec les expressions (8-1) à (8-3) et (9-1) à (9-3). Plus précisément, à chaque étape, les éléments du projecteur permettant de minimiser la métrique pénalisée sont obtenus au moyen d’un seuillage des éléments du projecteur minimisant la métrique non pénalisée, comme décrit en relation avec les expressions (10-1),(10-2) ; (13-1) à (13-3) ; (16-1) à (16-3). Il est important de noter qu’à chaque itération , seuls les éléments pénalisables sont pris en compte pour la minimisation de la métrique. A cet effet, on élimine de l’ensemble des indices pénalisables , les indices pour lesquels les éléments n’ont pas été pénalisée à l’itération , pour générer l’ensemble des indices pénalisables relatif à l’itération suivante, . Le tenseur du modèle prédictif, , est obtenu à partir des projecteurs , , comme dans la méthode de régression REW-NPLS conventionnelle. On trouvera le détail du calcul de à partir des projecteurs en question dans l’annexe A de l’article de A. Eliseyev et al. publié en 2017, cité dans la partie introductive. Les éléments du modèle prédictif correspondant aux indices pour lesquels les éléments ont été pénalisés, sont mis à zéro en 250. Après chaque itération , on élimine en 260 de et l’information redondante déjà prise en compte dans la régression sur l’espace de variables latentes, aux itérations précédentes . Le tenseur de covariance croisée dans lequel cette redondance a été éliminée est utilisée pour la nouvelle itération de décomposition PARAFAC. L’étape 260 est connue sous le nom de déflation dans la régression NPLS. Au terme des itérations, en 270, on dispose par conséquent de modèles prédictifs, , . Après l’acquisition de nouveaux blocs de données de calibration d’entrée, et de données de calibration de sortie, à l’étape de calibration suivante, 280, on détermine en 290, le modèle prédictif optimal, . Ce modèle prédictif est utilisé pour calculer le tenseur de commande à partir du tenseur d’observation, jusqu’à l’étape de calibration suivante, . La Fig. 3 représente de manière schématique le mécanisme de mise à jour du modèle prédictif dans la méthode de calibration de la Fig. 2. On a représenté en 310 la prise en compte d’un nouveau bloc de données de calibration d’entrée, représenté par le tenseur et d’un nouveau bloc de calibration de sortie, représenté par le tenseur . Les tenseurs de covariance et de covariance croisée, , en sont déduits en 320. Pour chaque valeur de , une décomposition PARAFAC du tenseur de covariance croisée est effectuée à l’aide d’une méthode des moindres carrés alternés (ALS) en 330, chaque étape de la décomposition PARAFAC portant sur un projecteur, cette décomposition étant régularisée au moyen d’un terme de pénalisation portant sur ce projecteur. La régularisation se traduit par un seuillage de certains éléments du projecteur obtenu par la méthode des moindres carrés conventionnelle, c’est-à-dire non pénalisée. A chaque itération PARAFAC représentée par la boucle 335, on en déduit un modèle prédictif en 340. Pour chaque modèle prédictif, on élimine en 350 l’information redondante dans le tenseur de covariance croisée (déflation) et l’on effectue une nouvelle itération PARAFAC sur le tenseur de covariance et le nouveau tenseur de covariance croisée (noté sur la figure) déflaté de ladite redondance, 360. La boucle de validation récursive est représentée en 380. Elle permet de déterminer parmi les modèles prédictifs , , obtenus à l’étape de calibration précédente, celui qui minimise l’erreur de prédiction sur les nouveaux blocs de données de calibration d’entrée, , et de sortie, . Méthode de calibration en ligne d’une interface neuronale directe destinée à recevoir une pluralité de signaux électrophysiologiques pendant une pluralité de fenêtres d’observation associées à des instants d’observation, pour former un tenseur d’entrée et en déduire, au moyen d’un modèle prédictif, un tenseur de sortie fournissant des signaux de commande destinés à commander un ou plusieurs effecteurs, ladite calibration est effectuée lors d’une pluralité d’étapes de calibration, chaque étape de calibration étant réalisée à partir de données de calibration d’entrée, représentées par un tenseur de calibration d’entrée, , et de données de calibration de sortie représentées par un tenseur de calibration de sortie , ladite étape de calibration mettant en œuvre une régression REW-NPLS comportant : le calcul (220) d’un tenseur de covariance à partir d’un tenseur de covariance obtenu lors d’une étape de calibration précédente et du tenseur , ainsi que le calcul d’un tenseur de covariance croisée à partir d’un tenseur de covariance croisée obtenu lors de ladite étape de calibration précédente et des tenseurs et ; une décomposition PARAFAC itérative du tenseur de covariance croisée , chaque itération PARAFAC étant associée à une dimension de l’espace des variables latentes sur lequel est effectuée ladite régression, chaque itération fournissant une pluralité de projecteurs de tailles respectives et un modèle prédictif défini par un tenseur de coefficients de prédiction et un tenseur de biais ; caractérisée en ce que chaque itération PARAFAC comprend une séquence de étapes élémentaires (240 1 , 240 1 ,.. 240 M ) de minimisation d’une métrique selon la méthode des moindres carrés alternés, chaque étape élémentaire de minimisation portant sur un projecteur et considérant les autres comme constants, ladite métrique comportant un terme de pénalisation fonction de la norme de ce projecteur, les éléments de ce projecteur ne faisant pas l’objet d’une pénalisation lors d’une itération PARAFAC n’étant pas pénalisables lors des itérations PARAFAC suivantes. Méthode de calibration en ligne d’une interface neuronale directe selon la revendication 1, caractérisée en ce qu’une étape élémentaire de minimisation portant sur un projecteur vise à minimiser la métrique pénalisée où est un projecteur non normalisé et , sont des projecteurs normalisés, est une matrice de taille obtenue par déploiement du tenseur de covariance croisée en prenant le mode comme référence, est le produit de Kronecker, est un coefficient strictement positif, , est une norme définie par , , où est le élément du projecteur , est le symbole de Kronecker, et est l’ensemble des indices correspondant aux éléments pénalisables du projecteur à l’itération PARAFAC Méthode de calibration en ligne d’une interface neuronale directe selon la revendication 1, caractérisée en ce qu’une étape élémentaire de minimisation portant sur un projecteur vise à minimiser les métriques pénalisées , où est le kème vecteur colonne de la matrice de taille obtenue par déploiement du tenseur de covariance croisée en prenant le mode comme référence, est le produit de Kronecker, est un coefficient strictement positif, et les fonctions , sont définies par : si et sinon si et sinon si et sinon où est l’ensemble des indices correspondant aux éléments pénalisables du projecteur à l’itération PARAFAC . Méthode de calibration en ligne d’une interface neuronale directe selon la revendication 3, caractérisée en ce que et que les éléments résultant de la minimisation sont donnés par : si et sinon où représente la solution de minimisation au sens des moindres carrés : et que le seuil est donné par . Méthode de calibration en ligne d’une interface neuronale directe selon la revendication 3, caractérisée en ce que et que les éléments résultant de la minimisation sont donnés par : si et si et sinon où et où représente la solution de minimisation au sens des moindres carrés : et le seuil est donné par . Méthode de calibration en ligne d’une interface neuronale directe selon la revendication 3, caractérisée en ce que et que les éléments résultant de la minimisation sont donnés par : si et si et sinon où représente la solution de minimisation au sens des moindres carrés : et le seuil est donné par Méthode de calibration en ligne d’une interface neuronale directe selon l’une des revendications précédentes, caractérisée en ce que les signaux électrophysiologiques sont des signaux ECoG.