EP4642334A1 - Dispositif de traitement de signaux de dispersion spatiotemporelle d'electrogrammes cardiaques - Google Patents

Dispositif de traitement de signaux de dispersion spatiotemporelle d'electrogrammes cardiaques

Info

Publication number
EP4642334A1
EP4642334A1 EP23841291.0A EP23841291A EP4642334A1 EP 4642334 A1 EP4642334 A1 EP 4642334A1 EP 23841291 A EP23841291 A EP 23841291A EP 4642334 A1 EP4642334 A1 EP 4642334A1
Authority
EP
European Patent Office
Prior art keywords
value
dispersion
signal
duration
spatiotemporal
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
EP23841291.0A
Other languages
German (de)
English (en)
Inventor
Thomas DEMARCY
Léandre GENIN
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
Publication of EP4642334A1 publication Critical patent/EP4642334A1/fr
Pending legal-status Critical Current

Links

Classifications

    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/316—Modalities, i.e. specific diagnostic methods
    • A61B5/318—Heart-related electrical modalities, e.g. electrocardiography [ECG]
    • A61B5/346—Analysis of electrocardiograms
    • A61B5/349—Detecting specific parameters of the electrocardiograph cycle
    • A61B5/361—Detecting fibrillation
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/25—Bioelectric electrodes therefor
    • A61B5/279—Bioelectric electrodes therefor specially adapted for particular uses
    • A61B5/28—Bioelectric electrodes therefor specially adapted for particular uses for electrocardiography [ECG]
    • A61B5/283—Invasive
    • A61B5/287—Holders for multiple electrodes, e.g. electrode catheters for electrophysiological study [EPS]
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/316—Modalities, i.e. specific diagnostic methods
    • A61B5/318—Heart-related electrical modalities, e.g. electrocardiography [ECG]
    • A61B5/367—Electrophysiological study [EPS], e.g. electrical activation mapping or electro-anatomical mapping
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7235—Details of waveform analysis
    • A61B5/7264—Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems

Definitions

  • the invention relates to the field of analysis of signals taken from cardiac electrograms.
  • This processing is based on the detection of a quantity of cardiac electrogram signals called dispersion by the Applicant. More precisely, the Applicant has discovered that the measurement of the spatiotemporal dispersion of cardiac electrograms, that is to say based both on the temporal evolution of the cardiac electrogram signals from each electrode, but also taking into account cardiac electrogram signals from neighboring electrodes, is particularly effective in determining the areas of the heart which are the cause of atrial fibrillation. The measurement of this quantity has been the subject of a patent granted in many countries around the world, and published in Europe under the number EP 3 236 843.
  • the Applicant has integrated the dispersion into a machine learning tool, in order to determine in near real time the presence of a dispersion within cardiac electrograms during a procedure, which makes it possible to alert the practitioner, accelerate the ablation procedure as well as the quality of the results obtained while reducing operational risks.
  • the reliability of the procedure is increased and makes it possible to limit the number of patients having to undergo another procedure for the same causes.
  • VX1 This work was integrated into software called VX1 which received approval from the FDA ("Food and Drug Administration") and CE marking.
  • the VX1 software without offering a diagnosis, assists the practitioner in carrying out the diagnosis: based on the dispersion measurements presented in the software, the practitioner can choose to mark areas of the heart to be ablated by catheter.
  • the VX1 software uses a database of annotated cardiac electrograms to train a machine learning engine which returns a value indicating the chance that an area for which a cardiac electrogram signal has been analyzed or the subject of atrial fibrillation. It is therefore not a measurement of dispersion strictly speaking, nor information constituting a medical diagnosis, but an indication that there is a more or less significant chance that a cardiac electrogram will demonstrate the occurrence of atrial fibrillation.
  • the value returned by the VX1 software is a value between 0 and 1 (0 indicating that there is almost zero probability that the area being measured participates in fibrillation, and 1 that this probability is a virtual certainty).
  • the software is configured such that a value of 0.5 is used as the threshold at which the practitioner's attention is drawn to the potential presence of dispersion.
  • the invention improves the situation.
  • a device for processing spatiotemporal dispersion signals of cardiac electrograms comprising a memory arranged to receive spatiotemporal dispersion signals of cardiac electrograms associated on the one hand with a time marker and on the other part of a cardiac electrogram track, a calculator arranged to receive as input an electrogram track identifier and a time marker, to analyze the spatiotemporal dispersion signal of cardiac electrograms associated with this electrogram track identifier, analyze a signal extract between the time marker and the first previous time marker whose spatiotemporal dispersion signal value of cardiac electrograms indicates an absence of dispersion, extract from this signal extract on the one hand a flatness value of the signal, and on the other hand a duration value taken from the duration of the signal extract, and return a track priority value calculated from the flatness value and the duration value, and a monitor arranged to receive the spatiotemporal dispersion signals d cardiac electrograms associated with each cardiac electrogram track, and,
  • This device is particularly advantageous because it makes it possible to offer additional information about the dispersion value which makes it possible to indicate to the practitioner a degree of priority for participation in the fibrillation of the area of the heart at the origin of the value. dispersion returned.
  • the track priority value is information which allows the practitioner to classify cardiac zones for ablation.
  • the track priority value will tend to increase when a practitioner stops in a given area and this is the location of dispersion.
  • the track priority value allows the practitioner to carry out his gesture more efficiently: when a dispersion is detected, he knows that it is desirable to stop on the area concerned, and to wait to see if the value track priority increases. If this is not the case, he can have confidence that it is desirable to continue his exploration.
  • the invention may have one or more of the following characteristics:
  • the spatiotemporal dispersion signals of cardiac electrograms are sequences of values taken from electrogram signals each indicating a degree of confidence in the fact that dispersion takes place for the electrogram track and the time marker concerned, and in which the monitor and/or the calculator are arranged to determine whether a spatiotemporal dispersion signal of cardiac electrograms indicates a relevant dispersion by comparing the value of the spatiotemporal dispersion signal of cardiac electrograms to a threshold value,
  • the calculator is arranged to calculate the flatness value from at least one value among the standard deviation of the signal extract, the total variation of the signal extract, the entropy of the signal extract or a value taken from one or more derivatives of order greater than or equal to one of the signal extract,
  • the calculator is arranged to calculate the duration value by comparing the duration of the signal extract to a minimum value and/or to a maximum value, and by returning the value 0 if the duration of the signal extract is less to the minimum value, re- turning the value 1 if the duration of the signal extract is greater than the maximum value, and returning a value between 0 and 1 otherwise,
  • the calculator is arranged to determine the value between 0 and 1 by applying to the duration of the signal extract a function of projection of the interval between the minimum value and the maximum value towards the interval between 0 and 1, which projection function being chosen from the group comprising affine functions, exponential functions, polynomials and threshold functions, and
  • the calculator receive spatiotemporal dispersion signals of cardiac electrograms associated on the one hand with a time marker and on the other hand with a cardiac electrogram track, b) determine whether a value of a spatiotemporal dispersion signal d cardiac electrograms indicate dispersion, c) if operation b) is negative, repeat it with a dispersion value having a later time marker, d) if operation b) is positive, analyze the spatiotemporal dispersion signal d cardiac electrograms associated with the corresponding electrogram track identifier by analyzing a signal extract between the time marker and the first preceding time marker whose spatio-temporal dispersion signal value of cardiac electrograms indicates an absence of dispersion, and by drawing from this signal extract on the one hand a flatness value of the signal, and on the other hand a duration value taken from the duration of the signal extract, and e) returning a priority value of track calculated from the flatness value and the duration value of operation d).
  • the process may have one or more of the following characteristics:
  • the spatiotemporal dispersion signals of cardiac electrograms are sequences of values taken from electrogram signals each indicating a degree of confidence in the fact that dispersion takes place for the electrogram track and the time marker concerned, and in which operation b) comprises the comparison of a value of the spatiotemporal dispersion signal of cardiac electrograms to a threshold value,
  • - operation d) comprises calculating the flatness value from at least one value among the standard deviation of the signal extract, the total variation of the signal extract, the entropy of the 'signal extract or a value taken from one or more derivatives of order greater than or equal to one of the signal extract,
  • - operation d) comprises calculating the duration value by comparing the duration of the signal extract to a minimum value and/or to a maximum value, and by returning the value 0 if the duration of the signal extract is lower than the minimum value, in re- turning the value 1 if the duration of the signal extract is greater than the maximum value, and returning a value between 0 and 1 otherwise,
  • - operation d) includes determining the value between 0 and 1 by applying to the duration of the signal extract a projection function from the interval between the minimum value and the maximum value towards the interval between 0 and 1, which projection function being chosen from the group comprising affine functions, exponential functions, polynomials and threshold functions, and
  • - operation e) comprises calculating the track priority value by producing a weighted harmonic, harmonic, weighted or arithmetic average of the flatness value and the duration value.
  • the invention also relates to a computer program comprising instructions for executing the method according to the invention, a data storage medium on which such a computer program is recorded and a computer system comprising a processor coupled to a memory, the memory having recorded such a computer program.
  • FIG.l represents a block diagram of a device according to the invention.
  • FIG.2 represents an example of implementation of an operating loop of the device of [Fig.l],
  • FIG.l represents a schematic example of a device 2 according to the invention.
  • the signals that are used by the device are based on cardiac electrograms measured by pairs of electrodes from a catheter in a patient's heart.
  • Device 2 includes a memory 4, a monitor 6 and a computer 8.
  • Memory 4 is arranged to receive all the data from device 2, whether input or output, of a global or local nature.
  • Memory 4 can be any type of data storage suitable for receiving digital data: hard disk, hard disk with flash memory, flash memory in any form, random access memory, magnetic disk, storage distributed locally or in the cloud, etc.
  • the memory 4 receives all the data which concerns the device 2, that is to say the programs and software instantiating the monitor 6 and the computer 8, the parameters and possible hyperparameters of those -here, the weights of any neural networks, the outputs and intermediate data of the neural networks, the spatiotemporal dispersion signal data of cardiac electrograms received as input (if applicable), the signal flatness and duration values signal, buffered data, and output track priority value data.
  • the data calculated by the device can be stored on any type of memory similar to memory 4, or on it. This data may be erased after the device has performed its tasks or retained.
  • the flatness value of the signal and the duration value are two values used to qualify whether the dispersion value signal is associated with a cardiac zone whose ablation has priority. These values are combined to produce a track priority value that indicates whether the dispersion value that is determined is associated with dispersion that has significant potential to be associated with an area that is the source of atrial fibrillation. Even more than in the other patent applications filed by the Applicant, the track priority value does not constitute a diagnosis but an indication which allows a doctor to make a decision, as blood pressure would be in another context.
  • the monitor 6 and the computer 8 directly or indirectly access the memory 4. They can be produced in the form of an appropriate computer code executed on one or more processors.
  • processors we must understand any processor adapted to the calculations described below.
  • Such a processor can be produced in any known manner, in the form of a microprocessor for a personal computer, laptop, tablet or smartphone, of a dedicated chip of the FPGA or SoC type, of a calculation resource on a grid or in the cloud, a cluster of graphics processors (GPUs), a microcontroller, or any other form capable of providing computing power necessary for the achievement described below.
  • GPUs graphics processors
  • microcontroller or any other form capable of providing computing power necessary for the achievement described below.
  • One or more of these elements can also be produced in the form of specialized electronic circuits such as an ASIC.
  • a combination of processor and electronic circuits can also be considered.
  • Processors dedicated to machine learning could also be considered.
  • the function of the monitor 6 is to analyze the data stream of spatiotemporal dispersion signals of cardiac electrograms received as input, and to detect the fact that a dispersion value indicates that treatment is necessary.
  • a dispersion value greater than or equal to 0.5 is significant.
  • the detection of this value will depend on the values taken by the dispersion signal received as input. For example, this could be generated in the opposite way to the VX1 software (for example 1 - Value from VX1), in which case, it would rather be a value less than or equal to 0.5 which would be significant. This determination could also be made differently, based on a value taken from the derivative of the dispersion value signal, or in any other relevant way.
  • the monitor 6 could for example measure the derivative of the input dispersion value signal and consider that a detection is positive if the derivative is positive, etc. Generally speaking, the monitor 6 could rely on several tests in order to qualify the detection of a relevant dispersion value.
  • monitor 6 The role of monitor 6 is therefore of an “interruptive” nature. Indeed, in the absence of a relevant dispersion value, it is useless to calculate a track priority value, since dispersion is not detected. On the other hand, as soon as a relevant dispersion value is detected, the monitor 6 calls the computer 8 to calculate the track priority value. Thus, the operation of the device 2 could appear as a loop of detections by the monitor 6 for each track, with the execution of the calculator 8 each time a relevant dispersion value is detected. Of course other implementations could be considered.
  • the role of the calculator 8 is to calculate the track priority value for a track whose dispersion value has been considered relevant by the monitor 6.
  • the dispersion signal is a signal which associates a dispersion marker time to a dispersion value.
  • This dispersion value is itself taken from an analysis of several cardiac electrogram signal values.
  • the dispersion values are updated approximately every 300ms, based on extracts of cardiac electrogram signals lasting approximately 1.5s. Alternatively, this update can take place every 100ms, every 500ms or other.
  • the determination of a track priority value is based on an extract of dispersion values which can have a continuous duration of 1.5s, or approximately 5 dispersion values, up to several tens of seconds, or around a hundred dispersion values.
  • the calculator 8 analyzes each time a signal extract of dispersion values which ends with the dispersion value which has just been determined by the monitor 6 as relevant.
  • This extract contains exclusively dispersion values including the time markers which are successive, associated with the same track, and which are considered relevant by the monitor 6. It will appear that this extract can be obtained in numerous ways:
  • the monitor 6 can generate extracts during its operation, by adding a current dispersion value detected as relevant to a current extract if the immediately preceding dispersion value has also been detected as relevant, or create a new extract in the case opposite,
  • the calculator 8 upon receipt of a dispersion value associated with a given time marker, can analyze a buffer of past dispersion values, and stop at the oldest value considered relevant by the monitor 6, or Again,
  • the calculator 8 can recover a buffer of past dispersion values from the time marker associated with a dispersion value received as input, and properly determine an extract from this buffer which it considers relevant.
  • the Applicant's work has demonstrated that the temporal continuity of relevant dispersion values is all the more useful as the dispersion value signal is precise. Indeed, with a dispersion considered "noisy", one could be tempted to ignore an irrelevant dispersion value to have more data allowing the calculation of the track priority value.
  • the Applicant's work has demonstrated that the combination of the signal flatness value and the duration value makes it possible to avoid artificially enlarging the size of the extracts and to obtain better results.
  • the calculator 8 operates by carrying out two measurements on the extract defined above, a measurement of the flatness value of the signal, and a measurement of the duration value. In both cases, the aim is to determine whether the dispersion value presents a form of stability over time. Indeed, the Applicant's work revealed that the signals of stable dispersion values were associated with priority zones in terms of ablation with a view to suppressing atrial fibrillation.
  • the flatness value of the signal is taken from the standard deviation of the extract data.
  • the Applicant's work has revealed that the standard deviation is the measurement providing the best results. Nevertheless, the Applicant has determined that other types of measurements could be retained, such as variance, total variation of the extract, entropy of the extract, a value taken from one or more higher order derivatives or equal to one of the extract, range, interquartile range or other similar measure.
  • the calculator 8 also determines a duration value which makes it possible to indicate the extent to which the current extract is long in relation to a time interval which is considered to indicate that the cardiac zone associated with this dispersion is a priority in terms of ablation.
  • the calculator 8 projects the duration of the extract onto a standard interval comprised between a minimum duration and a maximum duration.
  • the calculator 8 determines the duration of the extract and projects its value onto the interval [minimum duration; maximum duration] in order to determine a value between 0 and 1.
  • the projection can be of any type: linear, polynomial, exponential, threshold, etc. This involves indicating, for a given extract duration, whether this duration is characteristic of a priority cardiac zone or not.
  • the minimum duration has an obvious usefulness as a floor value.
  • the maximum duration has an important operational utility: when the first dispersion value of an extract longer than the maximum duration is detected, the output priority value will necessarily be low because the duration value will be low. As the clip grows, the Track Priority value will increase with the Duration value. When an operator sees that the track priority value is no longer moving, he or she can determine that it is because the duration value cannot increase any further, and that it is time to move the catheter.
  • the calculator 8 could determine the duration value differently, independently of the range [minimum duration; maximum duration]. Still alternatively, the minimum duration and the maximum duration could be variable during the procedure or be personalized for each patient.
  • the flatness value of the signal and the duration value can be determined in parallel with each other. Alternatively, one can be calculated before the other.
  • the value between 0 and 1 is chosen because of the way in which the calculator 8 determines the track priority value from the signal flatness value and the duration value. Indeed, the computer 8 operates in the example described by operating a weighted harmonic average. Alternatively, this average can be harmonic, weighted or arithmetic of the flatness value and the duration value.
  • FIG.2 represents an example of an operating loop of device 2.
  • monitor 6 is called with the current dispersion value for a given track.
  • an operation 210 is triggered in which the extract is determined from the time marker of the dispersion value of the operation 200, the associated track identifier , as well as dispersion value signals already received for this track identifier.
  • the calculator 8 can determine the flatness value and the duration value in an operation 220, then the track priority value and return this in an operation 230.
  • the device 2 will determine a track priority value by track for which a dispersion value signal is received as input. This makes it possible to enrich the information transmitted to the practitioner to allow him to carry out the diagnosis determining whether an area associated with a given track must be the subject of an ablation or not.

Landscapes

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

Abstract

Un dispositif de traitement de signaux de dispersion spatiotemporelle d'électrogrammes cardiaques comprend une mémoire (4) agencée pour recevoir des signaux de dispersion spatio- temporelle d'électrogrammes cardiaques associés d'une part à un marqueur de temps et d'autre part à une piste d'électrogramme cardiaque, un calculateur (8) agencé pour recevoir en entrée un identifiant de piste d'électrogramme et un marqueur de temps, pour analyser le signal de dispersion spatiotemporelle d'électrogrammes cardiaques associés à cet identifiant piste d'électrogramme, analyser un extrait de signal compris entre le marqueur de temps et le premier marqueur de temps précédent dont la valeur de signal de dispersion spatiotemporelle d'électrogrammes cardiaques indique une absence de dispersion, tirer de cet extrait de signal d'une part une valeur de planéité du signal, et d'autre part une valeur de durée tirée de la durée de l'extrait de signal, et retourner une valeur de priorité de piste calculée à partir de la valeur de planéité et de la valeur de durée, et un moniteur (6) agencé pour recevoir les signaux de dispersion spatiotemporelle d'électrogrammes cardiaques associés à chaque piste d'électrogramme cardiaque, et, lorsque la valeur d'un signal de dispersion spatiotemporelle d'électrogrammes cardiaques indique une dispersion pertinente, pour appeler le calculateur (8) avec le marqueur de temps et l'identifiant de piste d'électrogramme correspondants.

Description

Description
Titre de l’invention : Dispositif de traitement de signaux de dispersion spatiotemporelle d’électrogrammes cardiaques
[0001] L’invention concerne le domaine de l'analyse de signaux tirés d'électrogrammes cardiaques.
[0002] Le domaine du traitement des fibrillations atriales a connu des progrès importants dans la dernière décennie. Pour traiter les fibrillations atriales, les praticiens opèrent en introduisant des cathéters munis d'une pluralité d'électrodes. Ces électrodes sont déplacées au sein du cœur afin de mesurer le signal électrique qui se propage dans celui-ci. Les signaux obtenus sont appelés électrogrammes cardiaques. Ces électrogrammes cardiaques sont traités afin d'aider les praticiens à détecter la ou les zones du cœur qui sont à l'origine de la fibrillation atriale. Une fois ces zones identifiées, le praticien les brûle afin de les rendre inactives, ce qui a pour conséquence de rétablir le fonctionnement normal du cœur et de supprimer les fibrillations atriales.
[0003] La plupart des solutions existantes sont basées sur l'analyse des CFAE pour " Complex Fractionated Atrial Electrograms” (ou "Électrogrammes atriaux fractionnés complexes" en français). Le principe est d'essayer de trouver des emplacements des oreillettes où les électrogrammes perdent leur continuité, c’est-à-dire deviennent fractionnés.
[0004] La Demanderesse a développé un traitement des électrogrammes qui pique l'intérêt de la communauté scientifique depuis la publication d'études cliniques très positives. Celles-ci ont mené à la publication de plusieurs articles dont celui de Seitz et al. "AF Ablation Guided by Spatiotemporal Electrogram Dispersion Without Pulmonary Vein Isolation: A Wholly Patient-Tailored Approach” , Journal of the American College of Cardiology, Volume 69, Issue 3, 24 January 2017, pages 303-321.
[0005] Ce traitement repose sur la détection d'une grandeur des signaux d'électrogrammes cardiaques appelée dispersion par la Demanderesse. Plus précisément, la Demanderesse a découvert que la mesure de la dispersion spatiotemporelle d'électrogrammes cardiaques, c’est-à-dire basée à la fois sur l'évolution temporelle des signaux d'électrogrammes cardiaques de chaque électrode, mais en tenant compte également des signaux d'électrogrammes cardiaques d'électrodes voisines, est particulièrement efficace en vue de déterminer les zones du cœur qui sont à l'origine de la fibrillation atriale. La mesure de cette grandeur a fait l'objet d'un brevet délivré dans de nombreux pays du monde, et publié en Europe sous le numéro EP 3 236 843.
[0006] En prolongeant ses travaux, la Demanderesse a intégré la dispersion dans un outil d'apprentissage automatique, afin de déterminer en temps quasi-réel la présence d'une dispersion au sein d'électrogrammes cardiaques lors d'une procédure, ce qui permet d'alerter le praticien, d'accélérer la procédure d'ablation ainsi que la qualité des résultants obtenus tout en diminuant les risques opératoires. De plus, la fiabilité de la procédure est augmentée et permet de limiter le nombre de patients devant subir une autre procédure pour les mêmes causes.
[0007] Ces travaux ont été intégrés dans un logiciel dénommé VX1 qui a reçu une approbation de la FDA ("Food and Drug Administration"') et un marquage CE. Le logiciel VX1, sans offrir de diagnostic, assiste le praticien dans la réalisation du diagnostic : sur la base des mesures de dispersion présentées dans le logiciel, le praticien peut choisir de marquer des zones du cœur devant faire l'objet d'une ablation par cathéter.
[0008] Dans sa version actuelle, le logiciel VX1 utilise une base de données d'électrogrammes cardiaques annotés pour entraîner un moteur d'apprentissage automatique qui retourne une valeur indiquant la chance qu'une zone pour laquelle un signal d'électrogramme cardiaque a été analysé soit l'objet d'une fibrillation atriale. Il ne s'agit donc pas d'une mesure de la dispersion à proprement parler, ni d'une information constituant un diagnostic médical, mais d'une indication qu'il existe une chance plus ou moins importante qu'un électrogramme cardiaque témoigne de l'occurrence d'une fibrillation atriale. Ainsi, la valeur retournée par le logiciel VX1 est une valeur comprise entre 0 et 1 (0 indiquant qu'il n'y a qu'une probabilité quasi-nulle que la zone faisant l'objet de la mesure participe à la fibrillation, et 1 que cette probabilité est une quasi-certitude). Le logiciel est configuré de telle sorte qu'une valeur de 0,5 est utilisée comme seuil à partir duquel l'attention du praticien est attirée sur la présence potentielle d'une dispersion.
[0009] Au cours du développement du logiciel VX1, la Demanderesse a constaté que, dans certains cas, de nombreuses zones du cœur retournent une valeur supérieure à 0,5, ce qui complique le travail du praticien, puisque ce dernier doit procéder au moins d'ablation possible. De plus, lorsqu'un nombre trop important de zones à l'origine de dispersion sont identifiées, le praticien a tendance à perdre confiance dans la qualité de mesure du logiciel, puisque celui retourne une information hautement bruitée.
[0010] L'invention vient améliorer la situation. A cet effet, elle propose un dispositif de traitement de signaux de dispersion spatiotemporelle d'électrogrammes cardiaques comprenant une mémoire agencée pour recevoir des signaux de dispersion spatio- temporelle d'électrogrammes cardiaques associés d'une part à un marqueur de temps et d'autre part à une piste d'électrogramme cardiaque, un calculateur agencé pour recevoir en entrée un identifiant de piste d'électrogramme et un marqueur de temps, pour analyser le signal de dispersion spatiotemporelle d'électrogrammes cardiaques associés à cet identifiant piste d'électrogramme, analyser un extrait de signal compris entre le marqueur de temps et le premier marqueur de temps précédent dont la valeur de signal de dispersion spatiotemporelle d'électrogrammes cardiaques indique une absence de dispersion, tirer de cet extrait de signal d'une part une valeur de planéité du signal, et d'autre part une valeur de durée tirée de la durée de l'extrait de signal, et retourner une valeur de priorité de piste calculée à partir de la valeur de planéité et de la valeur de durée, et un moniteur agencé pour recevoir les signaux de dispersion spatiotemporelle d'électrogrammes cardiaques associés à chaque piste d'électrogramme cardiaque, et, lorsque la valeur d'un signal de dispersion spatiotemporelle d'électrogrammes cardiaques indique une dispersion pertinente, pour appeler le calculateur avec le marqueur de temps et l'identifiant de piste d'électrogramme correspondants.
[0011] Ce dispositif est particulièrement avantageux car il permet d'offrir une information complémentaire de la valeur de dispersion qui permet d'indiquer au praticien un degré de priorité de participation à la fibrillation de la zone du cœur à l'origine de la valeur de dispersion retournée. Ainsi, la valeur de priorité de piste est une information qui permet au praticien de classifier les zones cardiaques en vue de l'ablation.
[0012] De plus, la valeur de priorité de piste aura tendance à croître lorsqu'un praticien s'arrête sur une zone donnée et que celle-ci est le lieu d'une dispersion. Ainsi, la valeur de priorité de piste permet au praticien de réaliser son geste plus efficacement : lorsqu'une dispersion est détectée, il sait qu'il est souhaitable de s'arrêter sur la zone concernée, et d'attendre pour voir si la valeur de priorité de piste augmente. Si ce n'est pas le cas, il peut avoir confiance dans le fait qu'il est souhaitable de continuer son exploration.
[0013] Selon divers modes de réalisation, l’invention peut présenter une ou plusieurs des caractéristiques suivantes :
- les signaux de dispersion spatiotemporelle d'électrogrammes cardiaques sont des suites de valeurs tirées de signaux d'électrogramme indiquant chacune un degré de confiance dans le fait qu'une dispersion a lieu pour la piste d'électrogramme et le marqueur de temps concernés, et dans lequel le moniteur et/ou le calculateur sont agencés pour déterminer si un signal de dispersion spatiotemporelle d'électrogrammes cardiaques indique une dispersion pertinente en comparant la valeur du signal de dispersion spatiotemporelle d'électrogrammes cardiaques à une valeur seuil,
- le calculateur est agencé pour calculer la valeur de planéité à partir d'au moins une valeur parmi l'écart-type de l'extrait de signal, de la variation totale de l'extrait de signal, de l'entropie de l'extrait de signal ou d'une valeur tirée d'une ou plusieurs dérivées d'ordre supérieur ou égal à un de l'extrait de signal,
- le calculateur est agencé pour calculer la valeur de durée en comparant la durée de l'extrait de signal à une valeur minimale et/ou à une valeur maximale, et en retournant la valeur 0 si la durée de l'extrait de signal est inférieure à la valeur minimale, en re- tournant la valeur 1 si la durée de l'extrait de signal est supérieure à la valeur maximale, et en retournant une valeur comprise entre 0 et 1 sinon,
- le calculateur est agencé pour déterminer la valeur comprise entre 0 et 1 en appliquant à la durée de l'extrait de signal une fonction de projection de l'intervalle compris entre la valeur minimale et la valeur maximale vers l'intervalle compris entre 0 et 1, laquelle fonction de projection étant choisie parmi le groupe comprenant les fonctions affines, les fonctions exponentielles, les polynômes et les fonctions à seuil, et
- le calculateur recevoir des signaux de dispersion spatiotemporelle d'électrogrammes cardiaques associés d'une part à un marqueur de temps et d'autre part à une piste d'électrogramme cardiaque, b) déterminer si une valeur d'un signal de dispersion spatiotemporelle d'électrogrammes cardiaques indique une dispersion, c) si l'opération b) est négative, la répéter avec une valeur de dispersion présentant un marqueur de temps ultérieur, d) si l'opération b) est positive, analyser le signal de dispersion spatiotemporelle d'électrogrammes cardiaques associés à l'identifiant piste d'électrogramme correspondant en analysant un extrait de signal compris entre le marqueur de temps et le premier marqueur de temps précédent dont la valeur de signal de dispersion spatio- temporelle d'électrogrammes cardiaques indique une absence de dispersion, et en tirant de cet extrait de signal d'une part une valeur de planéité du signal, et d'autre part une valeur de durée tirée de la durée de l'extrait de signal, et e) retourner une valeur de priorité de piste calculée à partir de la valeur de planéité et de la valeur de durée de l'opération d).
[0014] Selon divers variantes, le procédé peut présenter une ou plusieurs des caractéristiques suivantes :
[0015] - les signaux de dispersion spatiotemporelle d'électrogrammes cardiaques sont des suites de valeurs tirées de signaux d'électrogramme indiquant chacune un degré de confiance dans le fait qu'une dispersion a lieu pour la piste d'électrogramme et le marqueur de temps concernés, et dans lequel l'opération b) comprend la comparaison d'une valeur du signal de dispersion spatiotemporelle d'électrogrammes cardiaques à une valeur seuil,
- l'opération d) comprend calculer la valeur de planéité à partir d'au moins une valeur parmi l'écart-type de l'extrait de signal, de la variation totale de l'extrait de signal, de l'entropie de l'extrait de signal ou d'une valeur tirée d'une ou plusieurs dérivées d'ordre supérieur ou égal à un de l'extrait de signal,
- l'opération d) comprend calculer la valeur de durée en comparant la durée de l'extrait de signal à une valeur minimale et/ou à une valeur maximale, et en retournant la valeur 0 si la durée de l'extrait de signal est inférieure à la valeur minimale, en re- tournant la valeur 1 si la durée de l'extrait de signal est supérieure à la valeur maximale, et en retournant une valeur comprise entre 0 et 1 sinon,
- l'opération d) comprend déterminer la valeur comprise entre 0 et 1 en appliquant à la durée de l'extrait de signal une fonction de projection de l'intervalle compris entre la valeur minimale et la valeur maximale vers l'intervalle compris entre 0 et 1, laquelle fonction de projection étant choisie parmi le groupe comprenant les fonctions affines, les fonctions exponentielles, les polynômes et les fonctions à seuil, et
- l'opération e) comprend calculer la valeur de priorité de piste en réalisant une moyenne harmonique pondérée, harmonique, pondérée ou arithmétique de la valeur de planéité et de la valeur de durée.
[0016] 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.
[0017] 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 [Fig.l] représente un schéma de principe d'un dispositif selon l'invention, et
- la [Fig.2] représente un exemple de mise en œuvre d'une boucle de fonctionnement du dispositif de la [Fig.l],
[0018] 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.
[0019] La [Fig.l] représente un exemple schématique d'un dispositif 2 selon l'invention. Comme cela a été précisé dans l'introduction, les signaux qui sont utilisés par le dispositif sont basés sur des électrogrammes cardiaques mesurés par des paires d'électrodes d'un cathéter dans le cœur d'un patient.
[0020] Cependant, dans le cas particulier de l'invention, ce ne sont pas ces signaux qui sont traités, mais la mesure de dispersion qui en est tirée. Comme indiqué en introduction, l'article de Seitz et al. "AF Ablation Guided by Spatiotemporal Electrogram Dispersion Without Pulmonary Vein Isolation: A Wholly Patient-Tailored Approach" , Journal of the American College of Cardiology, Volume 69, Issue 3, 24 January 2017, pages 303-321 et le brevet publié sous le numéro EP 3 236 843 permettent de mieux comprendre ce qu'est la dispersion tant en ce qui concernent les points communs que les différences avec les CFAE.
[0021] Compte tenu du domaine considéré et du fait que la dispersion et les CFAE représentent des phénomènes assez distincts dans la mesure où les CFAE ignorent tout aspect spatial, ces derniers ne présentent pas d'intérêt majeur dans le cadre de l'invention. En effet, comme on le verra plus bas, l'invention vise d'une manière générale à qualifier la stabilité du signal formé par les valeurs de dispersion. Cette même analyse n'aurait que peu voire pas de sens dans le cas des CFAE.
[0022] Le dispositif 2 comprend une mémoire 4, un moniteur 6 et un calculateur 8.
[0023] La mémoire 4 est agencée pour recevoir toutes les données du dispositif 2, qu'elles soient d'entrée ou de sortie, de nature globale ou locale. 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.
[0024] Dans l'exemple décrit ici, la mémoire 4 reçoit toutes les données qui concernent le dispositif 2, c’est-à-dire les programmes et logiciels instanciant le moniteur 6 et le calculateur 8, les paramètres et éventuels hyperparamètres de ceux-ci, les poids des éventuels réseaux de neurones, les sorties et données intermédiaires des réseaux de neurones, les données de signaux de dispersion spatiotemporelle d'électrogrammes cardiaques reçues en entrée (le cas échéant), les valeurs de planéité du signal et de durée de signal, les données stockées en mémoire tampon, ainsi que les données de valeur de priorité de piste en sortie. 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.
[0025] Comme on le verra plus bas, la valeur de planéité du signal et la valeur de durée sont deux valeurs utilisées pour qualifier si le signal de valeurs de dispersion est associé à une zone cardiaque dont l'ablation est prioritaire. Ces valeurs sont combinées afin de produire une valeur de priorité de piste qui indique si la valeur de dispersion qui est déterminée est associée à une dispersion qui présente un potentiel important d'être associée à une zone qui est la source d'une fibrillation atriale. Encore plus que dans les autres demandes de brevet déposées par la Demanderesse, la valeur de priorité de piste ne constitue pas un diagnostic mais une indication qui permet à un médecin de prendre une décision, comme le serait la pression sanguine dans un autre contexte.
[0026] Le moniteur 6 et le calculateur 8 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 processeur et de circuits électroniques peut également être envisagée. Des processeurs dédiés à l’apprentissage automatique pourront aussi être envisagés.
[0027] Le moniteur 6 a pour fonction d'analyser le flux de données de signaux de dispersion spatiotemporelle d'électrogrammes cardiaques reçus en entrée, et d'y détecter le fait qu'une valeur de dispersion indique qu'un traitement est nécessaire. Comme indiqué en introduction, dans le cadre du logiciel VX1, une valeur de dispersion supérieure ou égale à 0,5 est significative. Bien sûr, la détection de cette valeur dépendra des valeurs prises par le signal de dispersion reçu en entrée. Par exemple, celui-ci pourrait être généré de manière inverse au logiciel VX1 (par exemple 1 - Valeur issue de VX1), auquel cas, ce serait plutôt une valeur inférieure ou égale à 0,5 qui serait significative. Cette détermination pourrait également être réalisée différemment, sur la base d'une valeur tirée de la dérivée du signal de valeur de dispersion, ou de toute autre manière pertinente.
[0028] Du fait de la nature continue du traitement par le dispositif 2 qui apparaîtra mieux plus bas, une fois qu'une valeur de dispersion a été détectée par le moniteur 6, la détection pour les valeurs suivantes (mais concernant la même piste bien sûr) pourra être différent ou simplifié. Ainsi, dans le cas décrit plus haut, plutôt que comparer la valeur de dispersion courante à un seuil, le moniteur 6 pourra par exemple mesurer la dérivée du signal de valeur de dispersion en entrée et considérer qu'une détection est positive si la dérivée est positive, etc. D'une manière générale, le moniteur 6 pourrait reposer sur plusieurs tests afin de qualifier la détection d'une valeur de dispersion pertinente.
[0029] Le rôle du moniteur 6 est donc de nature "interruptive". En effet, en l'absence d'une valeur de dispersion pertinente, il est inutile de calculer une valeur de priorité de piste, puisqu'une dispersion n'est pas détectée. En revanche, dès qu'une valeur de dispersion pertinente est détectée, le moniteur 6 appelle le calculateur 8 pour calculer la valeur de priorité de piste. Ainsi, le fonctionnement du dispositif 2 pourra apparaître comme une boucle de détections par le moniteur 6 pour chaque piste, avec l'exécution du calculateur 8 à chaque fois qu'une valeur de dispersion pertinente est détectée. Bien sûr d'autres mises en œuvre pourront être envisagées.
[0030] Le calculateur 8 a pour rôle de calculer la valeur de priorité de piste pour une piste dont une valeur de dispersion a été considérée comme pertinente par le moniteur 6. Pour rappel, le signal de dispersion est un signal qui associe un marqueur de temps à une valeur de dispersion. Cette valeur de dispersion est elle-même tirée d'une analyse de plusieurs valeurs de signal d'électrogramme cardiaque. Dans le cadre du logiciel VX1, les valeurs de dispersions sont mises à jour environ toutes les 300ms, sur la base d'extraits de signaux d'électrogrammes cardiaques d'une durée d'environ 1,5s. En variante, cette mise à jour peut avoir lieu toutes les 100ms, toutes les 500ms ou autre.
[0031] Comme on le verra plus bas, la détermination d'une valeur de priorité de piste se base sur un extrait de valeurs de dispersion pouvant présenter une durée continue de 1,5s, soit environ 5 valeurs de dispersion, jusqu'à plusieurs dizaines de secondes, soit une centaine de valeurs de dispersion.
[0032] Plus précisément, le calculateur 8 analyse à chaque fois un extrait de signal de valeurs de dispersion qui se termine par la valeur de dispersion qui vient d'être déterminée par le moniteur 6 comme pertinente. Sur cet extrait contient exclusivement des valeurs de dispersion dont les marqueurs de temps qui sont successifs, associés à la même piste, et qui sont considérées comme pertinentes par le moniteur 6. Il apparaîtra que cet extrait peut être obtenu de nombreuses manières :
- le moniteur 6 peut générer des extraits au cours de son fonctionnement, en ajoutant une valeur de dispersion courante détectée comme pertinente à un extrait courant si la valeur de dispersion immédiatement précédente a également été détectée comme pertinente, ou créer un nouvel extrait dans le cas contraire,
- le calculateur 8, à réception d'une valeur de dispersion associée à un marqueur de temps donné, peut analyser un tampon de valeurs de dispersion passées, et s'arrêter à la valeur la plus ancienne considérée comme pertinente par le moniteur 6, ou encore,
- le calculateur 8 peut récupérer un tampon de valeurs de dispersions passées à partir du marqueur de temps associé à une valeur de dispersion reçue en entrée, et déterminer de manière propre un extrait de ce tampon qu'il considère comme pertinent.
[0033] Les travaux de la Demanderesse ont démontré que la continuité temporelle de valeurs de dispersion pertinentes est d'autant plus utile que le signal de valeur de dispersion est précis. En effet, avec une dispersion considérée comme "bruitée", on pourrait être tenté d'ignorer une valeur de dispersion non pertinente pour avoir plus de données permettant de calculer la valeur de priorité de piste. Les travaux de la Demanderesse ont démontré que la combinaison de la valeur de planéité de signal et de la valeur de durée permet de ne pas agrandir artificiellement la taille des extraits et d'obtenir de meilleurs résultats.
[0034] Le calculateur 8 opère en réalisant deux mesures sur l'extrait défini plus haut une mesure de valeur de planéité du signal, et une mesure de valeur de durée. Dans les deux cas, il s'agit de chercher à déterminer si la valeur de dispersion présente une forme de stabilité dans le temps. En effet, les travaux de la Demanderesse ont révélé que les signaux de valeurs de dispersion stables étaient associés à des zones prioritaires en termes d'ablation en vue de supprimer la fibrillation atriale.
[0035] Dans l'exemple décrit ici, la valeur de planéité du signal est tirée de l'écart-type des données de l'extrait. Les travaux de la Demanderesse ont révélé que l'écart-type est la mesure fournissant les meilleurs résultats. Néanmoins, la Demanderesse a déterminé que d'autres types de mesures pouvaient être retenues, comme la variance, la variation totale de l'extrait, l'entropie de l'extrait, une valeur tirée d'une ou plusieurs dérivées d'ordre supérieur ou égal à un de l'extrait, l'étendue, l'écart interquartile ou une autre mesure semblable.
[0036] En parallèle, le calculateur 8 détermine également une valeur de durée qui perme d'indiquer à quel point l'extrait courant est long par rapport à un intervalle de temps dont il est considéré comme indiquant que la zone cardiaque associée à cette dispersion est prioritaire en termes d'ablation. Dans l'exemple décrit ici, le calculateur 8 vient projeter la durée de l'extrait sur un intervalle étalon compris entre une durée minimale et une durée maximale. Ces deux valeurs, estimées de manière empirique par la Demanderesse, indiquent respectivement la durée minimale que doit avoir un extrait pour désigner une zone cardiaque prioritaire, et la durée maximale compte tenu du fait que le praticien ne peut pas se permettre de rester trop longtemps sur chaque zone s'il veut conduire son geste dans des délais raisonnables et minimiser le risque opératoire.
[0037] Dans l'exemple décrit ici, le calculateur 8 détermine la durée de l'extrait et projette sa valeur sur l'intervalle [durée minimale ; durée maximale] afin de déterminer une valeur comprise entre 0 et 1. La projection peut être de tout type : linéaire, polynomiale, exponentielle, à seuil, etc. Il s'agit ici d'indiquer, pour une durée d'extrait donné, si cette durée est caractéristique d'une zone cardiaque prioritaire ou pas.
[0038] La durée minimale présente une utilité évidente de valeur plancher. La durée maximale présente une utilité opératoire importante : lorsque la première valeur de dispersion d'un extrait plus long que la durée maximale est détectée, la valeur de priorité en sortie sera nécessairement basse car la valeur de durée sera faible. Au fur et à mesure que l'extrait grandit, la valeur de priorité de piste va s'accroître avec la valeur de durée. Lorsqu'un opérateur voit que la valeur de priorité de piste ne bouge plus, il peut déterminer que c'est parce que la valeur de durée ne peut plus augmenter, et qu'il est temps de déplacer le cathéter.
[0039] En variante, le calculateur 8 pourrait déterminer la valeur de durée différemment, indépendamment de la plage [durée minimale ; durée maximale]. Toujours en variante, la durée minimale et la durée maximale pourrait être variable au cours de la procédure ou être personnalisées pour chaque patient.
[0040] La valeur de planéité du signal et la valeur de durée peuvent être déterminées en parallèle l'une de l'autre. En variante, l'une peut être calculée avant l'autre.
[0041] La valeur comprise entre 0 et 1 est choisie en raison de la manière dont le calculateur 8 déterminer la valeur de priorité de piste à partir de la valeur de planéité de signal et de la valeur de durée. En effet, le calculateur 8 opère dans l'exemple décrit en opérant une moyenne harmonique pondérée. En variante, cette moyenne peut être harmonique, pondérée ou arithmétique de la valeur de planéité et de la valeur de durée.
[0042] Dans le cas où la valeur de planéité du signal et la valeur de durée ne sont pas comprises dans des plages de valeurs identiques, un redressement relatif peut être opéré, ou une autre formule peut être retenue pour calculer la valeur de priorité de piste.
[0043] La [Fig.2] représente un exemple de boucle de fonctionnement du dispositif 2. Dans une opération 200, le moniteur 6 est appelé avec la valeur de dispersion courante pour une piste donnée. Lorsque le moniteur 6 détermine qu'une valeur de dispersion est pertinente, une opération 210 est déclenchée dans laquelle l'extrait est déterminé à partir du marqueur de temps de la valeur de dispersion de l'opération 200, de l'identifiant de piste associé, ainsi que des signaux de valeurs de dispersion déjà reçus pour cet identifiant de piste. Une fois l'extrait déterminé, le calculateur 8 peut déterminer la valeur de planéité et la valeur de durée dans une opération 220, puis la valeur de priorité de piste et retourner celle-ci dans une opération 230.
[0044] Il apparaît que, le cas échéant, le dispositif 2 déterminera une valeur de priorité de piste par piste pour laquelle un signal de valeur de dispersion est reçu en entrée. Cela permet d'enrichir l'information transmise au praticien pour lui laisser la réalisation du diagnostic déterminant si une zone associée à une piste donnée doit faire l'objet d'une ablation ou pas.

Claims

Revendications
[Revendication 1] Dispositif de traitement de signaux de dispersion spatiotemporelle d'électrogrammes cardiaques comprenant une mémoire (4) agencée pour recevoir des signaux de dispersion spatiotemporelle d'électrogrammes cardiaques associés d'une part à un marqueur de temps et d'autre part à une piste d'électrogramme cardiaque, un calculateur (8) agencé pour recevoir en entrée un identifiant de piste d'électrogramme et un marqueur de temps, pour analyser le signal de dispersion spatio- temporelle d'électrogrammes cardiaques associés à cet identifiant piste d'électrogramme, analyser un extrait de signal compris entre le marqueur de temps et le premier marqueur de temps précédent dont la valeur de signal de dispersion spatiotemporelle d'électrogrammes cardiaques indique une absence de dispersion, tirer de cet extrait de signal d'une part une valeur de planéité du signal, et d'autre part une valeur de durée tirée de la durée de l'extrait de signal, et retourner une valeur de priorité de piste calculée à partir de la valeur de planéité et de la valeur de durée, et un moniteur (6) agencé pour recevoir les signaux de dispersion spatiotemporelle d'électrogrammes cardiaques associés à chaque piste d'électrogramme cardiaque, et, lorsque la valeur d'un signal de dispersion spatiotemporelle d'électrogrammes cardiaques indique une dispersion pertinente, pour appeler le calculateur (8) avec le marqueur de temps et l'identifiant de piste d'électrogramme correspondants.
[Revendication 2] Dispositif selon la revendication 1, dans lequel les signaux de dispersion spatiotemporelle d'électrogrammes cardiaques sont des suites de valeurs tirées de signaux d'électrogramme indiquant chacune un degré de confiance dans le fait qu'une dispersion a lieu pour la piste d'électrogramme et le marqueur de temps concernés, et dans lequel le moniteur (6) et/ou le calculateur (8) sont agencés pour déterminer si un signal de dispersion spatiotemporelle d'électrogrammes cardiaques indique une dispersion pertinente en comparant la valeur du signal de dispersion spatiotemporelle d'électrogrammes cardiaques à une valeur seuil.
[Revendication 3] Dispositif selon la revendication 1 ou 2, dans lequel le calculateur (8) est agencé pour calculer la valeur de planéité à partir d'au moins une valeur parmi l'écart-type de l'extrait de signal, de la variation totale de l'extrait de signal, de l'entropie de l'extrait de signal ou d'une valeur tirée d'une ou plusieurs dérivées d'ordre supérieur ou égal à un de l'extrait de signal.
[Revendication 4] Dispositif selon l'une des revendications précédentes, dans lequel le calculateur (8) est agencé pour calculer la valeur de durée en comparant la durée de l'extrait de signal à une valeur minimale et/ou à une valeur maximale, et en retournant la valeur 0 si la durée de l'extrait de signal est inférieure à la valeur minimale, en retournant la valeur 1 si la durée de l'extrait de signal est supérieure à la valeur maximale, et en retournant une valeur comprise entre 0 et 1 sinon.
[Revendication 5] Dispositif selon la revendication 4, dans lequel le calculateur (8) est agencé pour déterminer la valeur comprise entre 0 et 1 en appliquant à la durée de l'extrait de signal une fonction de projection de l'intervalle compris entre la valeur minimale et la valeur maximale vers l'intervalle compris entre 0 et 1, laquelle fonction de projection étant choisie parmi le groupe comprenant les fonctions affines, les fonctions exponentielles, les polynômes et les fonctions à seuil.
[Revendication 6] Dispositif selon l'une des revendications précédentes, dans lequel le calculateur (8) est agencé pour calculer la valeur de priorité de piste en réalisant une moyenne harmonique pondérée, harmonique, pondérée ou arithmétique de la valeur de planéité et de la valeur de durée.
[Revendication 7] Procédé de détermination d'une valeur de priorité de piste de signaux de dispersion spatiotemporelle d'électrogrammes cardiaques comprenant les opérations suivantes : a) recevoir des signaux de dispersion spatiotemporelle d'électrogrammes cardiaques associés d'une part à un marqueur de temps et d'autre part à une piste d'électrogramme cardiaque, b) déterminer si une valeur d'un signal de dispersion spatiotemporelle d'électrogrammes cardiaques indique une dispersion, c) si l'opération b) est négative, la répéter avec une valeur de dispersion présentant un marqueur de temps ultérieur, d) si l'opération b) est positive, analyser le signal de dispersion spatio- temporelle d'électrogrammes cardiaques associés à l'identifiant piste d'électrogramme correspondant en analysant un extrait de signal compris entre le marqueur de temps et le premier marqueur de temps précédent dont la valeur de signal de dispersion spatiotemporelle d'électrogrammes cardiaques indique une absence de dispersion, et en tirant de cet extrait de signal d'une part une valeur de planéité du signal, et d'autre part une valeur de durée tirée de la durée de l'extrait de signal, et e) retourner une valeur de priorité de piste calculée à partir de la valeur de planéité et de la valeur de durée de l'opération d).
[Revendication 8] Procédé selon la revendication 7, dans lequel les signaux de dispersion spatiotemporelle d'électrogrammes cardiaques sont des suites de valeurs tirées de signaux d'électrogramme indiquant chacune un degré de confiance dans le fait qu'une dispersion a lieu pour la piste d'électrogramme et le marqueur de temps concernés, et dans lequel l'opération b) comprend la comparaison d'une valeur du signal de dispersion spatiotemporelle d'électrogrammes cardiaques à une valeur seuil.
[Revendication 9] Procédé selon la revendication 7 ou 8, dans lequel l'opération d) comprend calculer la valeur de planéité à partir d'au moins une valeur parmi l'écart-type de l'extrait de signal, de la variation totale de l'extrait de signal, de l'entropie de l'extrait de signal ou d'une valeur tirée d'une ou plusieurs dérivées d'ordre supérieur ou égal à un de l'extrait de signal.
[Revendication 10] Procédé selon l'une des revendications 7 à 9, dans lequel l'opération d) comprend calculer la valeur de durée en comparant la durée de l'extrait de signal à une valeur minimale et/ou à une valeur maximale, et en retournant la valeur 0 si la durée de l'extrait de signal est inférieure à la valeur minimale, en retournant la valeur 1 si la durée de l'extrait de signal est supérieure à la valeur maximale, et en retournant une valeur comprise entre 0 et 1 sinon.
[Revendication 11] Procédé selon la revendication 10, dans lequel l'opération d) comprend déterminer la valeur comprise entre 0 et 1 en appliquant à la durée de l'extrait de signal une fonction de projection de l'intervalle compris entre la valeur minimale et la valeur maximale vers l'intervalle compris entre 0 et 1, laquelle fonction de projection étant choisie parmi le groupe comprenant les fonctions affines, les fonctions exponentielles, les polynômes et les fonctions à seuil.
[Revendication 12] Procédé selon l'une des revendications 7 à 11, dans lequel l'opération e) comprend calculer la valeur de priorité de piste en réalisant une moyenne harmonique pondérée, harmonique, pondérée ou arithmétique de la valeur de planéité et de la valeur de durée.
[Revendication 13] Programme informatique comprenant des instructions pour exécuter le procédé selon l’une des revendications 7 à 12 lorsqu'il est exécuté par ordinateur.
[Revendication 14] Support de stockage de données sur lequel est enregistré le programme informatique selon la revendication 13.
EP23841291.0A 2022-12-27 2023-12-19 Dispositif de traitement de signaux de dispersion spatiotemporelle d'electrogrammes cardiaques Pending EP4642334A1 (fr)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
FR2214567A FR3143969B1 (fr) 2022-12-27 2022-12-27 Dispositif de traitement de signaux de dispersion spatiotemporelle d'électrogrammes cardiaques
PCT/FR2023/052063 WO2024141729A1 (fr) 2022-12-27 2023-12-19 Dispositif de traitement de signaux de dispersion spatiotemporelle d'electrogrammes cardiaques

Publications (1)

Publication Number Publication Date
EP4642334A1 true EP4642334A1 (fr) 2025-11-05

Family

ID=85726658

Family Applications (1)

Application Number Title Priority Date Filing Date
EP23841291.0A Pending EP4642334A1 (fr) 2022-12-27 2023-12-19 Dispositif de traitement de signaux de dispersion spatiotemporelle d'electrogrammes cardiaques

Country Status (5)

Country Link
EP (1) EP4642334A1 (fr)
JP (1) JP2025542106A (fr)
FR (2) FR3143969B1 (fr)
IL (1) IL321142A (fr)
WO (1) WO2024141729A1 (fr)

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2206690A5 (fr) 1972-11-15 1974-06-07 Jeunet Jean Pierre
FR3030217B1 (fr) 2014-12-23 2021-09-03 Clement Bars Cartographie regionale haute densite du substrat de la fibrillation atriale
FR3119537B1 (fr) 2021-02-09 2023-04-07 Substrate Hd Dispositif informatique pour analyse en temps réel d’électrogrammes

Also Published As

Publication number Publication date
FR3143969A1 (fr) 2024-06-28
FR3143969B1 (fr) 2024-11-29
FR3143968B1 (fr) 2024-11-29
IL321142A (en) 2025-07-01
FR3143968A1 (fr) 2024-06-28
WO2024141729A1 (fr) 2024-07-04
JP2025542106A (ja) 2025-12-25

Similar Documents

Publication Publication Date Title
CN110393519B (zh) 心电信号的分析方法、装置、存储介质和处理器
CA2971775C (fr) Cartographie regionale haute densite du substrat de la fibrillation atriale
Masdjedi et al. Validation of novel 3‐dimensional quantitative coronary angiography based software to calculate fractional flow reserve post stenting
FR3119537A1 (fr) Dispositif informatique pour analyse en temps réel d’électrogrammes
EP4547105A1 (fr) Dispositif de détermination d'une longueur de cycle d'activation cardiaque
FR2831698A1 (fr) Station d'imagerie medicale a fonction d'extraction de trajectoire au sein d'un objet ramifie
WO2015121503A1 (fr) Procédé et système de surveillance du système nerveux autonome d'un sujet
EP3773161B1 (fr) Dispositif informatique de détection de troubles du rythme cardiaque
FR2807851A1 (fr) Procede et dispositif de detection de la fibrillation auriculaire cardiaque par la methode des ondelettes
EP4642334A1 (fr) Dispositif de traitement de signaux de dispersion spatiotemporelle d'electrogrammes cardiaques
EP3365868B1 (fr) Analyse numerique d'une image numerique representant une plaie pour sa caracterisation automatique
WO2023247734A1 (fr) Systeme et methode d'assistance a la realisation de mouvements physiques
EP4642333A1 (fr) Dispositif de traitement de signaux de dispersion spatiotemporelle d'electrogrammes cardiaques
Mohedano et al. Exploring EEG for object detection and retrieval
CN118395387B (zh) 特征融合的情感识别方法、装置、电子设备及存储介质
CH627358A5 (fr) Procede et appareillage d'analyse rapide d'un enregistrement electrocardiographique sur bande.
FR2983055A1 (fr) Detection et estimation du complexe qrs pour le suivi d'une activite cardiaque et pulmonaire
CN114190942B (zh) 由计算机执行的基于音频分析的抑郁症检测的方法
FR3139274A1 (fr) Dispositif de traitement de signaux cardiaques
EP3068282A1 (fr) Procédé et système de détermination d'endophénotypes caractéristiques d'une maladie, telle la schizophrénie
FR3066097A1 (fr) Dispositif et methode d'evaluation neuropsychologique
WO2025133123A1 (fr) Procede de reconstruction d'un ecg en temps reel a partir d'un signal ppg
WO2021084165A1 (fr) Procede et systeme de caracterisation d'un mouvement d'une entite en mouvement
CN120655604A (zh) 一种基于多阶段语义分割的主动脉斑块检测方法、设备及系统
CN109033167A (zh) 电影分类方法和系统

Legal Events

Date Code Title Description
STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: UNKNOWN

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: THE INTERNATIONAL PUBLICATION HAS BEEN MADE

PUAI Public reference made under article 153(3) epc to a published international application that has entered the european phase

Free format text: ORIGINAL CODE: 0009012

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: REQUEST FOR EXAMINATION WAS MADE

17P Request for examination filed

Effective date: 20250709

AK Designated contracting states

Kind code of ref document: A1

Designated state(s): AL AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HR HU IE IS IT LI LT LU LV MC ME MK MT NL NO PL PT RO RS SE SI SK SM TR

DAV Request for validation of the european patent (deleted)
DAX Request for extension of the european patent (deleted)