FR3131191A1 - Dispositif de traitement de signaux intracardiaques - Google Patents

Dispositif de traitement de signaux intracardiaques Download PDF

Info

Publication number
FR3131191A1
FR3131191A1 FR2114503A FR2114503A FR3131191A1 FR 3131191 A1 FR3131191 A1 FR 3131191A1 FR 2114503 A FR2114503 A FR 2114503A FR 2114503 A FR2114503 A FR 2114503A FR 3131191 A1 FR3131191 A1 FR 3131191A1
Authority
FR
France
Prior art keywords
wavelet
coefficient
qrs
instant
signal
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
FR2114503A
Other languages
English (en)
Inventor
Thomas Boudou
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Substrate HD SAS
Original Assignee
Substrate HD SAS
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Substrate HD SAS filed Critical Substrate HD SAS
Priority to FR2114503A priority Critical patent/FR3131191A1/fr
Priority to PCT/FR2022/052497 priority patent/WO2023118769A1/fr
Publication of FR3131191A1 publication Critical patent/FR3131191A1/fr
Pending legal-status Critical Current

Links

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/24Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/25Bioelectric electrodes therefor
    • A61B5/279Bioelectric electrodes therefor specially adapted for particular uses
    • A61B5/28Bioelectric electrodes therefor specially adapted for particular uses for electrocardiography [ECG]
    • A61B5/283Invasive
    • A61B5/287Holders for multiple electrodes, e.g. electrode catheters for electrophysiological study [EPS]
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/24Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/316Modalities, i.e. specific diagnostic methods
    • A61B5/318Heart-related electrical modalities, e.g. electrocardiography [ECG]
    • A61B5/346Analysis of electrocardiograms
    • A61B5/349Detecting specific parameters of the electrocardiograph cycle
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7203Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7235Details of waveform analysis
    • A61B5/7253Details of waveform analysis characterised by using transforms
    • A61B5/726Details of waveform analysis characterised by using transforms using Wavelet transforms

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Medical Informatics (AREA)
  • Surgery (AREA)
  • Biophysics (AREA)
  • Pathology (AREA)
  • Cardiology (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Veterinary Medicine (AREA)
  • Molecular Biology (AREA)
  • Physics & Mathematics (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Physiology (AREA)
  • Signal Processing (AREA)
  • Artificial Intelligence (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Psychiatry (AREA)
  • Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)

Abstract

Dispositif de traitement de signaux intracardiaques Un dispositif de traitement de signaux intracardiaques comprend une mémoire (4) agencée pour recevoir des données d’électrocardiogramme et des données d’électrogrammes synchronisées, un détecteur (6) agencé pour analyser les données d’électrocardiogramme et pour y détecter des instants d’onde QRS, un analyseur (8) agencé pour réaliser une transformation en ondelettes des données d’électrogrammes, un extracteur (10) agencé pour prélever dans la transformation en ondelettes des coefficients chacun associé à un instant d’onde QRS détecté par le détecteur (6), et pour les stocker dans un tampon (14), et un composeur (12) agencé pour extraire du tampon (14) un signal d’empreinte QRS et le soustraire dans la transformation en ondelettes aux instants d’onde QRS, et pour produire en sortie des données d’électrogrammes débruitées par transformation en ondelettes inverse du signal résultant. Fig.1

Description

Dispositif de traitement de signaux intracardiaques
L’invention concerne le domaine du traitement des signaux intracardiaques et leur débruitage.
La mesure des électrogrammes auriculaires à proximité des ventricules (par exemple dans la veine du sinus coronaire - le cathéter de référence utilisé pendant les procédures d'ablation), peut être perturbée par l'activité ventriculaire. Celle-ci génère des formes d'onde sur les électrogrammes intracardiaques appelées champ lointain (FF ou «far field» en anglais). Une telle distorsion par les ventricules rend l'analyse du signal des électrogrammes auriculaires plus difficile et ces informations non liées aux oreillettes y sont incluses à tort (par exemple, dans l’estimation de la longueur du cycle).
En particulier, l’invention vise à annuler la contribution ventriculaire en champ lointain dans les pistes d’électrogrammes issues de sondes intracardiaques, considéré comme du bruit, tout en conservant l'activité en champ proche («near-field» en anglais), même lorsqu’il y a superposition de ceux-ci.
Dans l’exercice clinique, c’est-à-dire dans les salles d’opérations, l’influence du champ lointain est traitée par effacement brutal du signal là où celui-ci est détecté. Cela est peu efficace, voire dégrade le signal.
La littérature scientifique offre plusieurs axes :
- la TMS (pour « Template matching and substraction » en anglais), qui consiste à moyenne un complexe QRS sur la durée d’un enregistrement, et à soustraire cette moyenne à chaque complexe QRS rencontré (voir par exemple Rieta, J. J., et al. "Atrial activity extraction based on blind source separation as an alternative to QRST cancellation for atrial fibrillation analysis." Computers in Cardiology 2000. Vol. 27 (Cat. 00CH37163), IEEE, 2000),
- l’analyse en composantes indépendantes (ou « ICA » pour «Independent Component Analysis» en anglais), qui consiste à trouver un ensemble de composantes qui minimise l'information mutuelle existant dans le signal ECG à différents intervalles. Les composantes indépendantes peuvent être regroupées dans les sous-espaces d’activité ventriculaire et d’activité atriale, ce qui permet de reconstruire l'activité atriale à chaque point d'observation à partir du sous-espace activité atriale (voir par exemple l’article de F. Castells, et. al. «Multidimensional ICA for the Separation of Atrial and Ventricular Activities from Single Lead ECGs in Paroxysmal Atrial Fibrillation Episodes», Lecture Notes in Computer science, vol. 3195, p.1229-1236, 2004),
- l’annulation ventriculaire adaptative (« AVC » pour «Adaptative Ventricular Cancellation» en anglais), qui consiste à utiliser un filtre à réponse impulsionnelle finie sur le canal de référence afin d’estimer les interférences, puis à retirer ces dernières du canal de référence (voir par exemple l’article de Widrow B., et al «Adaptive noise cancelling: principles and applications», Proc. IEEE 63, 1692–716, 1975),
- le couplage de l’ICA avec une décomposition en ondelettes, qui consiste à appliquer l’ICA sur les décompositions en ondelettes des électrogrammes et de l’électrocardiogramme (voir par exemple l’article de Simanto Saha, et.al. «A Ventricular Far-field Artefact Filtering Technique for Atrial Electrograms», 2019 Computing in Cardiology Conference, 2019, et
- l’analyse en composantes principales (« PCA » pour «Principal Component Analysis» en anglais), qui consiste à transformer des variables liées entre elles (dites « corrélées » en statistique) en nouvelles variables décorrélées les unes des autres (voir par exemple l’article de Christopher Schilling, «Analysis of Atrial Electrograms», Vol. 17 Karlsruhe Transactions on Biomedical Engineering, 2012.
Toutes ces méthodes présentent des désavantages considérables. Ainsi, la TMS déforme les électrogrammes atriaux, ce qui rend inusable toute méthode qui l’utilise comme base, l’ICA perd une grande partie de son efficacité dès lors que les signaux deviennent irréguliers et désorganisés, l’AVC est moins performante que l’ICA et l’efficacité de l’annulation dépend comme elle fortement du signal de référence, et la PCA vient améliorer la TMS, mais avec des défauts similaires.
Aucune méthode ne permet donc à ce jour de débruiter efficacement les signaux d’électrogrammes pour permettre leur exploitation dans un cadre d’arythmie des oreillettes.
L’invention vient améliorer la situation. À cet effet, elle propose un dispositif de traitement de signaux intracardiaques, comprenant une mémoire agencée pour recevoir des données d’électrocardiogramme et des données d’électrogrammes synchronisées, un détecteur agencé pour analyser les données d’électrocardiogramme et pour y détecter des instants d’onde QRS, un analyseur agencé pour réaliser une transformation en ondelettes des données d’électrogrammes, un extracteur agencé pour prélever dans la transformation en ondelettes des coefficients chacun associé à un instant d’onde QRS détecté par le détecteur, et pour les stocker dans un tampon, et un composeur agencé pour extraire du tampon un signal d’empreinte QRS et le soustraire dans la transformation en ondelettes aux instants d’onde QRS, et pour produire en sortie des données d’électrogrammes débruitées par transformation en ondelettes inverse du signal résultant.
Ce dispositif est particulièrement avantageux car il permet non seulement d’annuler l'activité du champ lointain, mais également de reconstruire le champ proche là où il se superpose avec le champ lointain, ce qui permet d’augmenter fortement le rapport signal sur bruit.
Selon divers modes de réalisation, l’invention peut présenter une ou plusieurs des caractéristiques suivantes :
- la mémoire est agencée pour recevoir des données d’électrogrammes qui correspondent à des pistes distinctes, le détecteur, l’analyseur, l’extracteur et le composeur étant agencés pour traiter indépendamment les données d’électrogrammes associées à des pistes distinctes,
- l’extracteur est agencé pour prélever des coefficients de tel sorte que, pour un coefficient choisi correspondant à un instant donné dans un niveau de transformée en ondelette donné, le signal d’ondelette qui correspond au niveau d’ondelette donné et qui est centré sur l’instant donné présente un recouvrement avec une fenêtre extraite des données d’électrogrammes dont est tiré le coefficient choisi, laquelle fenêtre est centrée sur l’instant d’onde QRS auquel est associé chaque coefficient respectif,
- l’extracteur est agencé pour pondérer les coefficients stockés dans le tampon en fonction du recouvrement temporel entre le signal d’ondelette qui correspond au niveau d’ondelette de chaque coefficient respectif et qui est centré sur l’instant correspondant à chaque coefficient respectif et la fenêtre centrée sur l’instant d’onde QRS auquel est associé chaque coefficient respectif,
- le composeur est agencé pour définir un signal d’empreinte QRS pour chaque niveau d’ondelette de la transformée en ondelettes, chacun à partir d’une fonction des coefficients prélevés par l’extracteur pour un niveau d’ondelette respectif,
- le composeur est agencé pour appliquer une fonction choisie dans le groupe comprenant la médiane géométrique, la PCA ou l’ICA,
- le composeur est agencé pour utiliser un nombre choisi de coefficients les plus récents dans le tampon pour chaque niveau d’ondelette, et
-l’extracteur est agencé pour réaliser une transformation en ondelettes de type SWT.
L’invention concerne également un procédé de traitement de signaux intracardiaques comprenant :
a) recevoir des données d’électrocardiogramme et des données d’électrogrammes synchronisées,
b) analyser les données d’électrocardiogramme pour y détecter des instants d’onde QRS,
c) réaliser une transformation en ondelettes des données d’électrogrammes,
d) prélever dans la transformation en ondelettes des coefficients, chacun associé à un instant d’onde QRS détecté à l’opération b), et les stocker dans un tampon,
e) extraire du tampon un signal d’empreinte QRS et le soustraire dans la transformation en ondelettes de l’opération c) aux instants d’onde QRS, et
f) réaliser une transformation en ondelettes inverse du signal de l’opération e), et retourner en sortie les données d’électrogrammes débruitées correspondantes.
Selon divers modes de réalisation, ce procédé peut présenter une ou plusieurs des caractéristiques suivantes :
- l’opération d) comprend prélever des coefficients de la transformation en ondelettes de l’opération c) tels que, pour un coefficient choisi correspondant à un instant donné d’un niveau donné de la transformée en ondelettes, le signal d’ondelette qui correspond au niveau donné de la transformée en ondelettes et qui est centré sur l’instant donné présente un recouvrement avec une fenêtre extraite des données d’électrogrammes dont est tiré le coefficient choisi centrée sur l’instant d’onde QRS auquel est associé le coefficient choisi,
- l’opération d) comprend pondérer chaque coefficient prélevé avant de le stocker dans le tampon, en fonction du recouvrement temporel entre le signal d’ondelette qui correspond au niveau de la transformation en ondelettes de chaque coefficient prélevé respectif et qui est centré sur l’instant correspondant à chaque coefficient prélevé respectif et la fenêtre centrée sur l’instant d’onde QRS associée à chaque coefficient prélevé respectif,
- l’opération e) comprend définir un signal d’empreinte QRS respectif pour chaque niveau d’ondelette de la transformée en ondelettes de l’opération b), chacun à partir d’une fonction des coefficients de l’opération d) correspondant à un niveau d’ondelette respectif,
- la fonction est choisie dans le groupe comprenant la médiane géométrique, la PCA ou l’ICA,
- la fonction utilise un nombre choisi de coefficients les plus récents dans le tampon pour chaque signal d’empreinte QRS, et
- l’opération c) réalise une transformation en ondelettes de type SWT
L’invention concerne également un programme informatique comprenant des instructions pour exécuter le procédé selon l'invention, un support de stockage de données sur lequel est enregistré un tel programme informatique et un système informatique comprenant un processeur couplé à une mémoire, la mémoire ayant enregistré un tel programme informatique.
D’autres caractéristiques et avantages de l’invention apparaîtront mieux à la lecture de la description qui suit, tirée d’exemples donnés à titre illustratif et non limitatif, tirés des dessins sur lesquels :
- la représente un schéma générique d’un dispositif selon l’invention,
- la représente un exemple de réalisation d’une fonction mise en œuvre par le dispositif de la ,
- la représente un exemple de signal d’électrogramme avant débruitage, après débruitage, la superposition ainsi que la différence entre ces signaux, et
- la représente un autre exemple de signal d’électrogramme avant débruitrage, après débruitage, la superposition ainsi que la différence entre ces signaux.
Les dessins et la description ci-après contiennent, pour l'essentiel, des éléments de caractère certain. Ils pourront donc non seulement servir à mieux faire comprendre la présente invention, mais aussi contribuer à sa définition, le cas échéant.
La représente un schéma générique d’un dispositif de traitement de signaux intracardiaques 2 selon l’invention. Le dispositif 2 comprend une mémoire 4, un détecteur 6, un analyseur 8, un extracteur 10, un composeur 12 et un tampon 14.
La mémoire 4 peut être tout type de stockage de données propre à recevoir des données numériques : disque dur, disque dur à mémoire flash, mémoire flash sous toute forme, mémoire vive, disque magnétique, stockage distribué localement ou dans le cloud, etc. Les données calculées par le dispositif peuvent être stockées sur tout type de mémoire similaire à la mémoire 4, ou sur celle-ci. Ces données peuvent être effacées après que le dispositif a effectué ses tâches ou conservées.
La mémoire 4 reçoit divers types de données :
- des données d’électrogrammes, qui représentent le signal mesuré sur une ou plusieurs pistes d’un cathéter cardiaque. Lorsque ces données proviennent de pistes distinctes, elles sont agencées de sorte que les données associées à chaque piste distincte peuvent être regroupée entre elles. Par exemple, il est classique que 5 pistes distinctes soient reçues,
- des données d’électrocardiogramme,
- des données d’instant d’onde QRS, dont la détermination sera explicitée plus bas,
- des données d’empreinte QRS, dont la détermination sera explicitée plus bas, et
- des données d’électrogrammes débruitées.
Le tampon 14 peut être réalisé comme un sous-ensemble de la mémoire 4, ou être séparé de celle-ci. Il pourra être mis en œuvre avec les mêmes moyens que ceux décrits ci-dessus en référence à la mémoire 4.
Le dispositif 2 exploite la diversité et la complémentarité des signaux disponibles dans une procédure d’électrophysiologie.
Comme on le verra plus bas, le détecteur 6 détermine d’abord la position des instants d’onde QRS en champ lointain sur l’électrocardiogramme. En effet, ceux-ci sont alignés avec le bruit de l'activité ventriculaire dans les dérivations intracardiaques. Ensuite, ce bruit est annulé des autres électrogrammes sans supprimer l'activité en champ proche, y compris lorsqu’il y a superposition du champ lointain et du champ proche, grâce à une analyse temps-fréquence.
L’analyseur 8 et l’extracteur 10 mettent en œuvre une approche temps-fréquence pour séparer les composantes de champ lointain, d'origine ventriculaire, des électrogrammes bipolaires enregistrés localement, c’est-à-dire le champ proche.
Des enregistrements simultanés d’électrocardiogrammes de surface et de signaux endocavitaires ont été utilisés pour segmenter le signal endocavitaire en plusieurs composantes temps-fréquence.
L’extracteur 10 collecte et stocke chaque activité ventriculaire de champ lointain, spécifique à chaque piste. Ces activités sont détectées grâce à leur synchronicité avec les instants d’onde QRS qui peuvent être détectés dans les données d’électrocardiogramme par une méthode de seuillage dans un domaine d'ondelettes. Elles sont stockées dans un tampon 14. En variante, l’algorithme de Pan-Tompkins pourrait être mis en œuvre pour détecter les instants d’onde QRS.
Au début, le tampon 14 est vide. Ensuite, il se remplit rapidement pour que l'on puisse calculer des empreintes respectives de chaque piste permettant d’annuler les activités de champs lointain.
Ensuite, le composeur 12 soustrait ces empreintes dans le domaine des ondelettes, puis reconstitue le signal débruité en effectuant une transformation en ondelettes inverse.
L'idée sous-jacente est que les activités auriculaires et ventriculaires peuvent être considérées comme statistiquement indépendantes et provenant de deux sources distinctes. Dans ce cas, les dérivations intracardiaques contiennent le mélange des contenus auriculaires et ventriculaires que le domaine des ondelettes permet de traiter efficacement.
Le détecteur 6, l’analyseur 8, l’extracteur 10, et le composeur 12 accèdent directement ou indirectement à la mémoire 4. Ils peuvent être réalisés sous la forme d’un code informatique approprié exécuté sur un ou plusieurs processeurs. Par processeurs, il doit être compris tout processeur adapté aux calculs décrits plus bas. Un tel processeur peut être réalisé de toute manière connue, sous la forme d’un microprocesseur pour ordinateur personnel, portable, tablette ou smartphone, d’une puce dédiée de type FPGA ou SoC, d’une ressource de calcul sur une grille ou dans le cloud, d’une grappe de processeurs graphiques (GPUs), d’un microcontrôleur, ou de toute autre forme propre à fournir la puissance de calcul nécessaire à la réalisation décrite plus bas. Un ou plusieurs de ces éléments peuvent également être réalisés sous la forme de circuits électroniques spécialisés tel un ASIC. Une combinaison de processeurs et de circuits électroniques peut également être envisagée.
La représente un exemple de mise en œuvre d’une fonction de débruitage mise en œuvre par le dispositif 2.
Dans une première opération 200, les données d’électrogrammes et les données d’électrocardiogramme synchronisées entre elles sont reçues. Typiquement, ces données sont reçues par tranches de 2s chacune.
Ensuite, dans une opération 210, le détecteur 6 exécute une fonction QRS() qui reçoit les données d’électrocardiogramme comme arguments, et retourne des instants d’onde QRS qui y sont détectés comme décrit plus haut.
Parallèlement, ou séquentiellement, l’analyseur 8 exécute une opération 220 dans laquelle une fonction SWT() reçoit les données d’électrogrammes comme argument et retourne une pluralité de décompositions en ondelettes, chacune correspondant à une piste distincte d’où sont issues les données d’électrogrammes. Dans l’exemple décrit ici, la fonction SWT() met en œuvre un algorithme du type algorithme à trous (« SWT » pour «Stationary Wavelet Transform» en anglais), qui présente l’avantage de combler l’absence d’invariance par translation de la transformée en ondelette discrète (« DWT » pour «Discrete Wavelet Transform» en anglais), et éviter de diviser la taille du signal par deux à chaque nouveau niveau de décomposition. En variante, une transformée en ondelette discrète pourrait être utilisée, ou une transformée en ondelette continue (« CWT » pour «Continuous Wavelet Transform» en anglais).
Une fois les transformées en ondelettes et les instants d’onde QRS obtenus, l’extracteur 10 exécute une fonction Buf() dans une opération 230. La fonction Buf() reçoit les transformées en ondelettes et les instants d’onde QRS et stocke dans le tampon 14 des extraits de ces transformées en ondelettes.
Chaque niveau d’ondelette de chaque piste est stocké séparément dans le tampon 14, de sorte qu’il est possible, pour chaque extrait dans le tampon 14, de déterminer à quelle piste et à quel niveau d’ondelette il correspond.
Cela est important, car la Demanderesse a découvert qu’il est intéressant de stocker les coefficients d’ondelettes qui correspondent aux instants d’onde QRS de chaque piste pour pouvoir retirer du signal la composante correspondant au champ lointain.
Plus précisément, une fenêtre QRS est définie autour de l’instant d’onde QRS qui définit la période pendant laquelle il est considéré que le signal d’électrogramme comprend du signal de champ lointain. Dans l’exemple décrit ici, cette fenêtre est centrée sur chaque instant d’onde QRS et dure 120ms, ce qui correspond à une durée commune des ondes QRS dans les connaissances médicales. En variante, la fenêtre pourrait être définie différemment, et par exemple être définie par la fonction QRS().
Chaque niveau d’ondelette est associé à une fonction d’ondelette qui présente une durée liée à la transformée en ondelettes utilisée ainsi qu’au niveau d’ondelette. Chaque coefficient d’ondelette de la transformée en ondelettes est donc associé à une fenêtre temporelle qui est centrée sur l’instant auquel ce coefficient est associé et qui présente la largeur de la fonction d’ondelette de ce niveau d’ondelette.
Ainsi, pour un niveau d'ondelette et une piste donnés, la fonction Buf() vient, à chaque instant d’onde QRS, sélectionner les coefficients de transformée en ondelettes dont la fenêtre temporelle présente un recouvrement avec la fenêtre QRS.
De manière optionnelle, la Demanderesse a découvert qu’il était avantageux de pondérer les coefficients d’ondelette stockés dans le tampon 14 en fonction de la quantité de recouvrement entre la fenêtre temporelle d’un coefficient donné et la fenêtre QRS.
En effet, pour les instants les plus à l’extrême de la fenêtre QRS, la fenêtre temporelle présente un recouvrement plus faible, puisqu’elle est centrée sur l’instant en question. Or, comme le signal correspondant à ce coefficient a vocation a être supprimé, il y a un risque de supprimer du signal qui ne correspond pas au champ lointain. De même, les ondelettes de niveau important peuvent présenter une fenêtre temporelle plus large que la fenêtre QRS, ce qui pose évidemment des problèmes similaires. La pondération correspondant à la quantité relative de recouvrement permet de limiter ces effets de bord.
Toujours en variante, un compromis peut être de limiter le niveau d’ondelette de sorte que la fenêtre temporelle la plus large présente une taille de l’ordre de celle de la fenêtre QRS.
Une fois le tampon 14 suffisamment rempli, par exemple lorsque le tampon 14 contient au moins 5 extraits de l’opération 230 pour chaque piste, le composeur 12 peut exécuter une fonction Out() dans une opération 240. Dans la fonction Out(), le composeur 12 détermine une médiane pour chaque piste et niveau en ondelette des extraits qui la concerne, et soustrait cette médiane de la transformation en ondelettes des données d’électrogrammes de la piste et du niveau en ondelette correspondant. Par médiane, il doit être entendu toute méthode permettant de déterminer une valeur représentative de l’ensemble des extraits, qu’elle soit obtenue par le biais d’une médiane géométrique, d’une PCA ou d’une ICA. Enfin, la fonction Out() exécute une transformation en ondelettes inverse du signal dans laquelle la médiane a été soustraite, et retourne le signal débruité correspondant.
La et la permettent de montrer le gain apporté par le dispositif 2, en particulier dans les zones dans lesquelles le champ lointain et le champ proche se superposent. Sur ces figures, la Demanderesse a représenté en haut le signal des données d’électrogrammes (c’est-à-dire non débruité), le signal débruité, la superposition et la différence entre ces signaux qui mettent en valeur l’effet de la soustraction des signaux d’empreintes QRS dans la transformée en ondelettes inversée.
Ces figures montrent l’efficacité du dispositif 2, qui vient complètement recréer le signal en champ proche là où il y avait superposition avec le champ lointain, et sans dégrader le signal en dehors de ces zones.
En outre, la Demanderesse a mené une étude afin de quantifier les gains obtenus grâce au dispositif 2. Pour cela, la Demanderesse a comparé trois méthodes pour analyser leurs pertinences, à partir d’électrogrammes dont la durée de cycle est connue. Pour cela, la Demanderesse a pris les signaux d’électrogrammes correspondants et en a créé trois copies : une première copie inchangée, une deuxième copie dans laquelle le signal aux instants d’onde QRS ventriculaire est supprimé, c'est à dire remplacé par une ligne isoélectrique gardant la continuité du signal, et une troisième copie débruitée avec le dispositif 2. Ensuite, les signaux résultants ont été utilisés pour déterminer la durée du cycle, puis comparés avec la durée réelle connue du cycle.
Cette étude a été réalisée sur une cohorte de 52 patients qui pouvaient être en rythme sinusal, en tachycardie atriale ou en fibrillation atriale. Le score qui a été calculé pour chaque copie est la racine carrée de l’erreur moyenne, selon la formule
Où c est l’indice de la copie considérée, yc(i) est la durée de cycle déterminée à partir du i-ème signal d’électrogramme de la copie d’indice c, et y(i) est la durée de cycle réelle.
Les scores sont les suivants : 13,36% pour la première copie, 8,45% pour la deuxième copie, et 4,06% pour la troisième copie.
Cela démontre le gain de rapport signal sur bruit obtenu grâce au dispositif 2.

Claims (15)

  1. Dispositif de traitement de signaux intracardiaques, comprenant une mémoire (4) agencée pour recevoir des données d’électrocardiogramme et des données d’électrogrammes synchronisées, un détecteur (6) agencé pour analyser les données d’électrocardiogramme et pour y détecter des instants d’onde QRS, un analyseur (8) agencé pour réaliser une transformation en ondelettes des données d’électrogrammes, un extracteur (10) agencé pour prélever dans la transformation en ondelettes des coefficients chacun associé à un instant d’onde QRS détecté par le détecteur (6), et pour les stocker dans un tampon (14), et un composeur (12) agencé pour extraire du tampon (14) un signal d’empreinte QRS et le soustraire dans la transformation en ondelettes aux instants d’onde QRS, et pour produire en sortie des données d’électrogrammes débruitées par transformation en ondelettes inverse du signal résultant.
  2. Dispositif selon la revendication 1, dans lequel la mémoire (4) est agencée pour recevoir des données d’électrogrammes qui correspondent à des pistes distinctes, le détecteur (6), l’analyseur (8), l’extracteur (10) et le composeur (12) étant agencés pour traiter indépendamment les données d’électrogrammes associées à des pistes distinctes.
  3. Dispositif selon la revendication 1 ou 2, dans lequel l’extracteur (10) est agencé pour prélever des coefficients de tel sorte que, pour un coefficient choisi correspondant à un instant donné dans un niveau de transformée en ondelette donné, le signal d’ondelette qui correspond au niveau d’ondelette donné et qui est centré sur l’instant donné présente un recouvrement avec une fenêtre extraite des données d’électrogrammes dont est tiré le coefficient choisi, laquelle fenêtre est centrée sur l’instant d’onde QRS auquel est associé chaque coefficient respectif.
  4. Dispositif selon la revendication 3, dans lequel l’extracteur (10) est agencé pour pondérer les coefficients stockés dans le tampon (14) en fonction du recouvrement temporel entre le signal d’ondelette qui correspond au niveau d’ondelette de chaque coefficient respectif et qui est centré sur l’instant correspondant à chaque coefficient respectif et la fenêtre centrée sur l’instant d’onde QRS auquel est associé chaque coefficient respectif.
  5. Dispositif selon l’une des revendications précédentes, dans lequel le composeur (12) est agencé pour définir un signal d’empreinte QRS pour chaque niveau d’ondelette de la transformée en ondelettes, chacun à partir d’une fonction des coefficients prélevés par l’extracteur (10) pour un niveau d’ondelette respectif.
  6. Dispositif selon la revendication 5, dans lequel le composeur (12) est agencé pour appliquer une fonction choisie dans le groupe comprenant la médiane géométrique, la PCA ou l’ICA.
  7. Dispositif selon la revendication 5 ou 6, dans lequel le composeur (12) est agencé pour utiliser un nombre choisi de coefficients les plus récents dans le tampon (14) pour chaque niveau d’ondelette.
  8. Dispositif selon l’une des revendications précédentes, dans lequel l’extracteur (8) est agencé pour réaliser une transformation en ondelettes de type SWT.
  9. Procédé de traitement de signaux intracardiaques comprenant :
    a) recevoir des données d’électrocardiogramme et des données d’électrogrammes synchronisées,
    b) analyser les données d’électrocardiogramme pour y détecter des instants d’onde QRS,
    c) réaliser une transformation en ondelettes des données d’électrogrammes,
    d) prélever dans la transformation en ondelettes des coefficients, chacun associé à un instant d’onde QRS détecté à l’opération b), et les stocker dans un tampon (14),
    e) extraire du tampon (14) un signal d’empreinte QRS et le soustraire dans la transformation en ondelettes de l’opération c) aux instants d’onde QRS, et
    f) réaliser une transformation en ondelettes inverse du signal de l’opération e), et retourner en sortie les données d’électrogrammes débruitées correspondantes.
  10. Procédé selon la revendication 9, dans lequel l’opération d) comprend prélever des coefficients de la transformation en ondelettes de l’opération c) tels que, pour un coefficient choisi correspondant à un instant donné d’un niveau donné de la transformée en ondelettes, le signal d’ondelette qui correspond au niveau donné de la transformée en ondelettes et qui est centré sur l’instant donné présente un recouvrement avec une fenêtre extraite des données d’électrogrammes dont est tiré le coefficient choisi centrée sur l’instant d’onde QRS auquel est associé le coefficient choisi.
  11. Procédé selon la revendication 10, dans lequel l’opération d) comprend pondérer chaque coefficient prélevé avant de le stocker dans le tampon (14), en fonction du recouvrement temporel entre le signal d’ondelette qui correspond au niveau de la transformation en ondelettes de chaque coefficient prélevé respectif et qui est centré sur l’instant correspondant à chaque coefficient prélevé respectif et la fenêtre centrée sur l’instant d’onde QRS associée à chaque coefficient prélevé respectif.
  12. Procédé selon l’une des revendications 9 à 11, dans lequel l’opération e) comprend définir un signal d’empreinte QRS respectif pour chaque niveau d’ondelette de la transformée en ondelettes de l’opération b), chacun à partir d’une fonction des coefficients de l’opération d) correspondant à un niveau d’ondelette respectif.
  13. Procédé selon la revendication 12, dans lequel la fonction est choisie dans le groupe comprenant la médiane géométrique, la PCA ou l’ICA.
  14. Programme informatique comprenant des instructions pour exécuter le procédé selon l’une des revendications 9 à 13 lorsque ledit programme est exécuté sur un ordinateur.
  15. Support de stockage de données sur lequel est enregistré le programme informatique selon la revendication 14.
FR2114503A 2021-12-24 2021-12-24 Dispositif de traitement de signaux intracardiaques Pending FR3131191A1 (fr)

Priority Applications (2)

Application Number Priority Date Filing Date Title
FR2114503A FR3131191A1 (fr) 2021-12-24 2021-12-24 Dispositif de traitement de signaux intracardiaques
PCT/FR2022/052497 WO2023118769A1 (fr) 2021-12-24 2022-12-23 Dispositif de traitement de signaux intracardiaques

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
FR2114503A FR3131191A1 (fr) 2021-12-24 2021-12-24 Dispositif de traitement de signaux intracardiaques
FR2114503 2021-12-24

Publications (1)

Publication Number Publication Date
FR3131191A1 true FR3131191A1 (fr) 2023-06-30

Family

ID=80933404

Family Applications (1)

Application Number Title Priority Date Filing Date
FR2114503A Pending FR3131191A1 (fr) 2021-12-24 2021-12-24 Dispositif de traitement de signaux intracardiaques

Country Status (2)

Country Link
FR (1) FR3131191A1 (fr)
WO (1) WO2023118769A1 (fr)

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20040111121A1 (en) * 2002-12-04 2004-06-10 Brown Mark L. Methods and apparatus for discriminating polymorphic tachyarrhythmias from monomorphic tachyarrhythmias facilitating detection of fibrillation

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20040111121A1 (en) * 2002-12-04 2004-06-10 Brown Mark L. Methods and apparatus for discriminating polymorphic tachyarrhythmias from monomorphic tachyarrhythmias facilitating detection of fibrillation

Non-Patent Citations (7)

* Cited by examiner, † Cited by third party
Title
CHRISTOPHER SCHILLING: "Analysis of Atrial Electrograms", KARLSRUHE TRANSACTIONS ON BIOMÉDICAL ENGINEERING, vol. 17, 2012
EL BOUNY LAHCEN ET AL: "A Wavelet Denoising and Teager Energy Operator-Based Method for Automatic QRS Complex Detection in ECG Signal", CIRCUITS, SYSTEMS AND SIGNAL PROCESSING, CAMBRIDGE, MS, US, vol. 39, no. 10, 7 April 2020 (2020-04-07), pages 4943 - 4979, XP037221692, ISSN: 0278-081X, [retrieved on 20200407], DOI: 10.1007/S00034-020-01397-8 *
F. CASTELLS: "Multidimensional ICA for the Se-paration ofatrial and Ventricular Activities from Single Lead ECGs in Paroxysmal Atrial Fibrillation Episodes", LECTURE NOTES IN COMPUTER SCIENCE, vol. 3195, 2004, pages 1229 - 1236
RIETA, J. J. ET AL.: "Computers in Cardiology", vol. 27, 2000, IEEE, article "Atrial activity extraction based on blind source séparation as an alternative to QRST can-cellation for atrial fibrillation analysis", pages: 2000
SAHA SIMANTO ET AL: "A Ventricular Far-Field Artefact Filtering Technique for Atrial Electrograms", 2019 COMPUTING IN CARDIOLOGY (CINC), CREATIVE COMMONS, 8 September 2019 (2019-09-08), XP033720823, DOI: 10.23919/CINC49843.2019.9005813 *
SIMANTO SAHA: "A Ventricular Far-field Artefact Filtering Technique for Atrial Electrograms", COMPUTING IN CARDIOLOGY CONFÉRENCE, 2019, pages 2019
WIDROW B. ET AL.: "Adaptive noise cancelling: principles and applications", PROC. IEEE, vol. 63, 1975, pages 1692 - 716, XP000762555

Also Published As

Publication number Publication date
WO2023118769A1 (fr) 2023-06-29

Similar Documents

Publication Publication Date Title
Slimane et al. QRS complex detection using empirical mode decomposition
Nagendra et al. Application of wavelet techniques in ECG signal processing: an overview
FR2847361A1 (fr) Dispositif d'analyse d'un signal, notamment d'un signal physiologique tel qu'un signal ecg
Satija et al. Noise‐aware dictionary‐learning‐based sparse representation framework for detection and removal of single and combined noises from ECG signal
CN111657936B (zh) 基于小波变换和全变差正则化的信号去噪方法
EP3552146B1 (fr) Détection de nerfs dans une série d'images échographiques
Islam et al. Wavelet based denoising algorithm of the ECG signal corrupted by WGN and Poisson noise
WO2001080042A1 (fr) Procede et dispositif de detection de la fibrillation auriculaire cardiaque par la methode des ondelettes
Yadav et al. Denoising and SNR improvement of ECG signals using wavelet based techniques
Li et al. Signal denoising based on the Schrödinger operator's eigenspectrum and a curvature constraint
FR3131191A1 (fr) Dispositif de traitement de signaux intracardiaques
Slim et al. Removal of ECG interference from surface respiratory electromyography
JP2008229307A (ja) 脳波測定方法及び脳波測定装置並びに記録媒体
El Bouny et al. ECG noise reduction based on stationary wavelet transform and zero-crossings interval thresholding
Abibullaev et al. Classification of Cardiac Arrhythmias using Biorthogonal Wavelets and Support Vector Machines.
Castells et al. Atrial activity extraction from atrial fibrillation episodes based on maximum likelihood source separation
Agrawal et al. Wavelet subband dependent thresholding for denoising of phonocardiographic signals
FR2983055A1 (fr) Detection et estimation du complexe qrs pour le suivi d'une activite cardiaque et pulmonaire
US8670821B2 (en) Electroencephalogram signal processing method
Qureshi et al. Noise reduction of electrocardiographic signals using wavelet transforms
Chauhan et al. Optimal choice of thresholding rule for denoising ECG using DWT
Linhares Arrhythmia detection coefficient based on wavelet shrinkage
WO2024003509A1 (fr) Dispositif de détermination d'une longueur de cycle d'activation cardiaque
CN109919878B (zh) 基于光学相干断层成像的降噪方法及其系统
Ramesh et al. ECG Signal Enhancement through Subband Adaptive Soft Thresholding and EMD for Efficient Cardiac Arrhythmia Analysis.

Legal Events

Date Code Title Description
PLFP Fee payment

Year of fee payment: 2

PLSC Publication of the preliminary search report

Effective date: 20230630

PLFP Fee payment

Year of fee payment: 3