JP2024506633A - 電位図をリアルタイム解析するコンピュータ装置 - Google Patents
電位図をリアルタイム解析するコンピュータ装置 Download PDFInfo
- Publication number
- JP2024506633A JP2024506633A JP2023548220A JP2023548220A JP2024506633A JP 2024506633 A JP2024506633 A JP 2024506633A JP 2023548220 A JP2023548220 A JP 2023548220A JP 2023548220 A JP2023548220 A JP 2023548220A JP 2024506633 A JP2024506633 A JP 2024506633A
- Authority
- JP
- Japan
- Prior art keywords
- electrogram
- evaluator
- probability
- probability array
- real
- 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
Links
- 238000010801 machine learning Methods 0.000 claims abstract description 20
- 239000006185 dispersion Substances 0.000 claims abstract description 15
- 230000000877 morphologic effect Effects 0.000 claims abstract description 5
- 238000013527 convolutional neural network Methods 0.000 claims description 12
- 238000004590 computer program Methods 0.000 claims description 6
- 238000000034 method Methods 0.000 claims description 6
- 238000003491 array Methods 0.000 claims description 4
- 238000013500 data storage Methods 0.000 claims description 2
- 238000010586 diagram Methods 0.000 abstract description 5
- 238000012731 temporal analysis Methods 0.000 abstract description 4
- 206010003658 Atrial Fibrillation Diseases 0.000 description 19
- 239000013598 vector Substances 0.000 description 13
- 238000004458 analytical method Methods 0.000 description 9
- 239000000284 extract Substances 0.000 description 7
- 238000002679 ablation Methods 0.000 description 5
- 230000004913 activation Effects 0.000 description 5
- 238000001994 activation Methods 0.000 description 5
- 238000004422 calculation algorithm Methods 0.000 description 5
- 238000013459 approach Methods 0.000 description 3
- 206010003119 arrhythmia Diseases 0.000 description 3
- 210000002837 heart atrium Anatomy 0.000 description 3
- 210000002569 neuron Anatomy 0.000 description 3
- 238000001356 surgical procedure Methods 0.000 description 3
- 238000012549 training Methods 0.000 description 3
- 238000013528 artificial neural network Methods 0.000 description 2
- 230000001746 atrial effect Effects 0.000 description 2
- 230000000747 cardiac effect Effects 0.000 description 2
- 238000013153 catheter ablation Methods 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000002955 isolation Methods 0.000 description 2
- 230000003902 lesion Effects 0.000 description 2
- 238000013507 mapping Methods 0.000 description 2
- 210000003492 pulmonary vein Anatomy 0.000 description 2
- 230000008685 targeting Effects 0.000 description 2
- ORILYTVJVMAKLC-UHFFFAOYSA-N Adamantane Natural products C1C(C2)CC3CC1CC2C3 ORILYTVJVMAKLC-UHFFFAOYSA-N 0.000 description 1
- 230000005856 abnormality Effects 0.000 description 1
- 230000006793 arrhythmia Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 206010061592 cardiac fibrillation Diseases 0.000 description 1
- 238000013145 classification model Methods 0.000 description 1
- 229940079593 drug Drugs 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 238000002565 electrocardiography Methods 0.000 description 1
- 230000007831 electrophysiology Effects 0.000 description 1
- 238000002001 electrophysiology Methods 0.000 description 1
- 230000002600 fibrillogenic effect Effects 0.000 description 1
- 230000001939 inductive effect Effects 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000011176 pooling Methods 0.000 description 1
- 238000004393 prognosis Methods 0.000 description 1
- 238000011897 real-time detection Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 230000000717 retained effect Effects 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 230000001225 therapeutic effect Effects 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Images
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
-
- 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/339—Displays specially adapted therefor
-
- 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/35—Detecting specific parameters of the electrocardiograph cycle by template matching
-
- 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/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/7253—Details of waveform analysis characterised by using transforms
- A61B5/7257—Details of waveform analysis characterised by using transforms using Fourier transforms
-
- 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
- A61B5/7267—Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems involving training the classification device
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
- G16H50/20—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for computer-aided diagnosis, e.g. based on medical expert systems
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
- G16H50/30—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for calculating health indices; for individual health risk assessment
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
- G16H50/70—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for mining of medical data, e.g. analysing previous cases of other patients
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Public Health (AREA)
- Medical Informatics (AREA)
- Biomedical Technology (AREA)
- Cardiology (AREA)
- General Health & Medical Sciences (AREA)
- Pathology (AREA)
- Physics & Mathematics (AREA)
- Heart & Thoracic Surgery (AREA)
- Biophysics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Veterinary Medicine (AREA)
- Data Mining & Analysis (AREA)
- Physiology (AREA)
- Databases & Information Systems (AREA)
- Epidemiology (AREA)
- Primary Health Care (AREA)
- Artificial Intelligence (AREA)
- Mathematical Physics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Psychiatry (AREA)
- Signal Processing (AREA)
- Evolutionary Computation (AREA)
- Fuzzy Systems (AREA)
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
- Debugging And Monitoring (AREA)
- Complex Calculations (AREA)
Abstract
電位図をリアルタイムで解析するためのコンピュータ装置に関する。該装置は、複数の電極からそれぞれ発信されるリアルタイム電位図信号を受信するように配置されたメモリと、抽出器、および勾配ブースティングに基づく機械学習モジュールを含む第1の評価器と、を備える。抽出器は、電位図信号のセット内の各電位図信号から、少なくとも1つの時間的解析特徴および少なくとも1つの形態学的特徴を含む特徴のセットを抽出するように、および得られた特徴のセットを前記勾配ブースティングに基づく機械学習モジュールに供給するように配置される。機械学習モジュールは、関連する電位図信号が分散を示すかどうかを示す値でラベル付けされた特徴のセットを含むデータに基づいて訓練され、前記電位図信号のセットごとに出力を行うように配置される。【選択図】図1
Description
本発明は、電位図をリアルタイム解析するコンピュータ装置に関する。
心房細動の分野は、この5年の間に急速に発展してきた。
心房細動を促進する心臓部位の検出には、ToperaやCardiolnsightなどのソフトウェアを使用するオフラインのコンピューティング装置によるものと、CARTOなどのソフトウェアを使用するリアルタイムのコンピューティング装置によるものが含まれる。
Toperaソフトウェアは、心臓の心房の電気的活性化を再構築することを目的としている。信号の取得は、「バスケットカテーテル」(心房全体上に展開されるカテーテル)を用いて行われる。解析は、取得開始から2分以上経過した遅延時間で行われる。このタイプのカテーテルは挿入が難しく、電極の接触が保証されないので、この解決策は不便である。解析時間が非常に長く、再構築では高精細なマップが得られず、複雑な電気的活性化(心房細動の70%を占める)の場合には単純すぎる再構築のために失敗することもある。
Cardiolnsightソフトウェアは、患者の皮膚に装着する多電極心電図測定ベストを使用して、心臓の電気活動を再構築することを目的としている。患者は手術前にベストを着用し、データは手術中に解析され、且つアクセス可能である。解析は、取得開始後、遅延時間(15分以上)で行われる。この解決策には計算時間が非常に長いという欠点があり、延期時間を生じさせてセットアップが困難になり(手術前にデータを抽出できるように、患者は数日前に来てベストを着用する必要がある)、コストがかなりかかる。また、再構築では、測定のセットアップが離れすぎているため高精細なマップが得られず、複雑な電気的活性化の場合には単純すぎる再構築のために失敗する。
Biosense&Webster社によるCARTOソフトウェアは、心臓不整脈を検出するアルゴリズムを実現することができる。
CFAE(複雑性分裂心房電位)アルゴリズムは、信号中の偏位点の数(導関数の符号の変化)を計算し、リアルタイムでカラーマップを作製する。医師は、これらの複雑なマップを解釈して、対象となる部位を決定する。このアルゴリズムはあまり特殊ではなく、比較的単純である。
Abbott Laboratories社によるEnsiteソフトウェアは、心臓不整脈を検出するアルゴリズムを実現することができる。
Rippleアルゴリズム(Rippleマッピング)は、CARTOソフトウェアのモジュールであり、最初のカテーテル通過後の心房の電気活動を再現することができる。このモジュールにより、電気波の振幅と伝播の両方が可視化される。その結果、医師が理解するには非常に複雑なマップが作製される。これは単純なケースでは有効かもしれないが、複雑な電気的活性化の場合にはマップの解釈が難しすぎるため、解析は失敗する。
本出願人は、改良されたシステムを発明し、心房細動を促進する心臓部位をリアルタイムで検出するために、連動する2つの方法を使用することを記載したフランス共和国特許出願第1852850号を提出した。より詳細には、この発明には、電位図データに2つの分類モデルを適用するように配置された分類器の使用が記載されている。一方のモデルは他方のモデルよりも適用が速く、精度が低いので、これは、心臓領域の潜在的な対象部位を医師に警告するために使用される。
本出願人のさらなる研究により、心房細動を促進する心臓部位をリアルタイムで検出するための改良されたコンピュータ装置を決定することが可能になった。この装置は、複数の電極からそれぞれ発信されるリアルタイム電位図信号を受信するように配置されたメモリと、抽出器、および勾配ブースティングに基づく機械学習モジュールを含む第1の評価器と、を備える。抽出器は、電位図信号のセット内の各電位図信号から、少なくとも1つの時間的解析特徴および少なくとも1つの形態学的特徴を含む特徴のセットを抽出するように、および得られた特徴のセットを勾配ブースティングに基づく機械学習モジュールに供給するように配置される。機械学習モジュールは、関連する電位図信号が分散を示すかどうかを示す値でラベル付けされた特徴のセットを含むデータに基づいて訓練され、電位図信号のセットごとに、電位図信号のセットのそれぞれの電位図信号が分散を示すかどうかを示す第1の確率配列を出力するように配置される。また、該装置は、第2の評価器を備える。第2の評価器は、畳み込みニューラルネットワークを含む。畳み込みニューラルネットワークは、リアルタイム電位図信号のセットを受信し、入力されたリアルタイム電位図信号が分散を示す第2の確率配列を出力する。畳み込みニューラルネットワークは、それぞれの電位図信号が分散を示すかどうかを示す値でそれぞれラベル付けされた電位図信号のデータベースに基づいて訓練される。また、該装置は、予測器を備える。予測器は、電位図信号の所与のセットに対して決定された第1の確率配列および第2の確率配列に基づいて、第1の確率配列および第2の確率配列の値の加重平均に少なくとも部分的に基づく第3の確率配列を返す。
このコンピュータ装置が特に興味深いのは、本出願人がリアルタイムで患者の分散マップを構築し、患者らの心房細動の終結率が88%に達するようにこれらのマップを使用することができたからである。これは、分散マップを構築して心房細動を治療するために、リアルタイムで医師を支援するように本出願人の発明を確実に使用することができることを意味する。
様々な実施形態において、コンピュータ装置は、以下の特徴のうちの1つまたは複数を含んでもよい。
・ 抽出器は、第1の周期長推定値、第2の周期長推定値、および最大振幅を有する電位図信号の高速フーリエ変換内の周波数からなるグループの中から少なくとも1つの時間的解析特徴を抽出するように配置され、
・ 抽出器は、電位図信号のユークリッドノルム、および電位図信号の積分された絶対導関数からなるグループの中から少なくとも1つの形態学的特徴を抽出するように配置され、
・ さらに、該装置は、リアルタイム電位図信号を、選択された持続時間を有する一連の電位図信号に分割するように配置され、
・ さらに、該装置は、同じ選択された持続時間を有する電位図信号のセットを提供するように配置され、
・ さらに、該装置は、リアルタイム電位図信号を、選択された持続時間を有する一連の電位図信号に分割するように配置され、
・ さらに、予測器は、同じ電位図信号に対応する、第1の確率配列中の確率と第2の確率配列中の確率との差の絶対値がしきい値を超えたときに、第3の確率配列において第1の確率配列中の確率を使用するように配置され、
・ さらに、該装置は、第3の配列中の値に関連する色を決定するように配置されており、電極ごとに、第3の確率配列における対応する電位図信号について決定された確率に関連する色を出力するように配置されたディスプレイをさらに備える。
・ 抽出器は、第1の周期長推定値、第2の周期長推定値、および最大振幅を有する電位図信号の高速フーリエ変換内の周波数からなるグループの中から少なくとも1つの時間的解析特徴を抽出するように配置され、
・ 抽出器は、電位図信号のユークリッドノルム、および電位図信号の積分された絶対導関数からなるグループの中から少なくとも1つの形態学的特徴を抽出するように配置され、
・ さらに、該装置は、リアルタイム電位図信号を、選択された持続時間を有する一連の電位図信号に分割するように配置され、
・ さらに、該装置は、同じ選択された持続時間を有する電位図信号のセットを提供するように配置され、
・ さらに、該装置は、リアルタイム電位図信号を、選択された持続時間を有する一連の電位図信号に分割するように配置され、
・ さらに、予測器は、同じ電位図信号に対応する、第1の確率配列中の確率と第2の確率配列中の確率との差の絶対値がしきい値を超えたときに、第3の確率配列において第1の確率配列中の確率を使用するように配置され、
・ さらに、該装置は、第3の配列中の値に関連する色を決定するように配置されており、電極ごとに、第3の確率配列における対応する電位図信号について決定された確率に関連する色を出力するように配置されたディスプレイをさらに備える。
また、本発明は、第1の評価器、第2の評価器、および予測器を実現するための命令を含むコンピュータプログラムと、このコンピュータプログラムを記録したデータ記憶媒体と、リアルタイム電位図信号を受信し、本発明による第1の評価器、第2の評価器、および予測器を実行し、第3の確率配列を返す、コンピュータで実現される方法に関する。
以下、本発明の他の特徴および利点を、本発明による例示的な実施形態を示す添付の図面を参照して説明する。
本発明によるコンピュータ装置の一実施形態に関する全体図である。
図1の実施形態で使用される畳み込みニューラルネットワークの模式図である。
添付の図面および以下の説明には、主に、本発明の有効且つ良好に定義された特徴が含まれる。したがって、これらは、本発明をよりよく理解するのに役立つだけでなく、その定義に適宜貢献ためにも使用され得る。
持続性心房細動は、依然として治療上の課題を有する。カテーテルアブレーションの適応となる患者の数はますます増加しているが、持続性心房細動に対する最適なアブレーション戦略はいまだ不明である。STAR-AF2(Verma A,Jiang C-y,Betts TRら、「Approaches to catheter ablation for persistent atrial fibrillation」、New England Journal of Medicine、2015;372:1812-1822)は、肺静脈隔離術(PVI)を超えるアブレーション病変の追加処置は、長期予後の改善にはつながらないと結論付けている。
その一方で、いくつかの観察研究(Jadidi AS,Lehrmann H,Keyl Cら、「Ablation of persistent atrial fibrillation targeting low-voltage areas with selective activation characteristics」、Circulation:Arrhythmia and Electrophysiology、2016;9:e002962、およびSeitz J,Bars C,Theodore Gら、「AF ablation guided by spatiotemporal electrogram dispersion without pulmonary vein isolation:a wholly patient-tailored approach」、Journal of the American College of Cardiology、2017;69:303-321参照)は、電位図異常を有する領域を標的とする非PVI病変が持続性心房細動患者にとって有益であることを示唆している。
特に、後者の研究では、多極性電位図の分散を示す心房領域への高周波エネルギーの適用を支持する結果が示された。他の研究(Narayan SM,Krummen DE,Shivkumar K,Clopton P,Rappel W-J,Miller JM、「Treatment of atrial fibrillation by the ablation of localized sources:CONFIRM(Conventional Ablation for Atrial Fibrillation With or Without Focal Impulse and Rotor Modulation)trial」、Journal of the American College of Cardiology 2012;60:628-636)では、その代わりに高度な信号解析が行われ、心房細動のドライバーを誘導するためのメカニズムに基づいた表示が行われた。
電位図の異常や解析にかかわらず、同様の電位図に基づくマッピングアプローチを実現しようとする施設間では、経験や学習曲線の段階に大きな差がある。
図1は、本発明によるコンピュータ装置の一実施形態に関する全体図である。
コンピュータ装置2は、メモリ4と、抽出器8および勾配ブースティングに基づく機械学習モジュール10を含む第1の評価器6と、第2の評価器12と、予測器14と、を備える。
使用上の観点から、操作者は、患者の心臓に挿入された電極を使用する。電極は、参照符号16によって示され、コンピュータ装置2のディスプレイ上に表示される。図1に示すように、16個の様々なタイプの電極が設けられてもよい。
電極による感知結果は、メモリ4に供給されるリアルタイム電位図信号18のセットである。これらの電位図信号が解析され、予測器14が電極の修正図20を返す。これにより、医師は、どの電極が分散に関連する可能性のある領域を感知したかを特定することができる。その見返りとして、医師は、患者の心臓上の対応する位置に印をつける(または関連する位置を自動的に返す)ことができる。これにより、ドットで示されている心房細動の観点から治療すべき領域を有する心臓のマップ22が得られる。
メモリ4は、リアルタイム電位図信号(例えば、すべての電位図信号が同じ長さ、例えば3秒を有するように、一定時間間隔でデジタル化および分割された後)と、それらが分散を示すかどうか、すなわち、心房細動に関連するかどうかを示す値でラベル付けされた275000以上の電位図信号のデータベースとを格納する。同一の持続時間を有する電位図信号は、連続していてもよいし、互いに重なっていてもよい。また、メモリ4は、本発明を実行する過程で生成され得る任意の一時的データを格納してもよい。また、メモリ4は、予測器14の動作の結果として得られたマップを、場合によっては医師が作成した注釈およびカテーテル形状と組み合わせて格納してもよい。
本明細書に記載の実施例において、メモリ4は、ハードディスクドライブ、ソリッドステートドライブ、フラッシュメモリ、プロセッサに組み込まれたメモリ、またはクラウドでアクセス可能な遠隔ストレージなど、適切な方法で実現されてもよい。
本明細書に記載の実施例において、抽出器8および勾配ブースティングに基づく機械学習モジュール10を含む第1の評価器6、第2の評価器12、ならびに予測器14は、1つまたは複数のプロセッサで実行されるコンピュータプログラムである。このようなプロセッサは、CPU、GPU、CPUおよび/またはGPUグリッド、遠隔計算グリッド、特別に構成されたFPGA、特別に構成されたASIC、SOCまたはNOCなどの専用チップ、またはAI専用チップなど、自動計算を実行するために知られている任意の手段を含む。
本明細書に記載の実施例において、第1の評価器6は、機械学習モジュールである。機械学習モジュールは、電位図信号から特徴を抽出し、分散確率を決定するために勾配ブースティング手法にそれらを通すように配置される。
抽出器8は、リアルタイム電位図信号のそれぞれの特徴のセットを決定するように配置される。これは、機械学習手法を適用した場合に優れた分散検出を達成することができることを本出願人が特定した。これらの特徴は、2つの主要な特徴カテゴリに分類されてもよい。
第1のカテゴリは、電位図信号の時間的解析、すなわち、電位図信号の時間領域または周波数領域の解析に本質的に関連する特徴に関する。本明細書に記載の実施例において、抽出器8は、第1の周期長推定値、第2の周期長推定値、および/または最も高い振幅を有する電位図信号の高速フーリエ変換内の周波数を抽出してもよい。後者の場合、振幅も特徴として保持されてもよい。
第1の周期長推定値および第2の周期長推定値は、以下のように決定されてもよい。
リアルタイム電位図信号データは、ベクトルに格納され、その各要素は、そのサンプルに対応する。次に、自己相関パラメータTを使用して、最初のT個のベクトルサンプルを含む第1のベクトル、および最後のT個のベクトルサンプルを含む第2のベクトルを含む、サイズTを有する2つのベクトルを画定する。
本明細書に記載の実施例において、第1の推定値は、第1のベクトルと第2のベクトルのスカラー積を、第1のベクトルと第2のベクトルのユークリッドノルムの積で割って正規化自己相関を測定することによって算出される。Tを変化させることにより、第1の推定値の最大値と最小値が決定され、このT値で計算された第1の推定値が、第1の推定値の最大値から、第1の推定値の最大値と最小値の差の0.3倍を引いた値よりも大きくなる第1の局所極値をT値が提供する値として、第1の推定値に選択される値が選択される。
代替的に、例えば0.3の係数を変更することによって、あるいは経験的、網羅的、あるいはその他の方法で推定値を最適化するT値を探すことによって、第1の推定値を決定することもできる。
本明細書に記載の実施例において、第2の推定値は、各T値について、第1のベクトルと第2のベクトルとの差のユークリッドノルムの2乗を、第1のベクトルの絶対値が最大となるサンプルと第2のベクトルの絶対値が最大となるサンプルとの積で割った値を計算することによって算出される。第1の推定値と同様に、第2の推定値の最大値と最小値が決定され、このT値で計算された第2の推定値が、第2の推定値の最小値に第2の推定値の最大値と最小値の差の0.2倍を加算した値よりも小さくなる第1の局所極値をT値が提供する値として、第2の推定値に選択される値が選択される。代替的に、例えば0.2の係数を変更することによって、あるいは経験的、網羅的、あるいはその他の方法で推定値を最適化するT値を探すことによって、第1の推定値を決定することもできる。
第2のカテゴリは、信号の形態学的解析、すなわち、数学的形態学に基づく、またはそれに関連する解析に関する。本明細書に記載の実施例において、抽出器8は、電位図信号のユークリッドノルムおよび電位図信号の積分された絶対導関数を抽出してもよい。
上述したすべての特徴を使用することで、勾配ブースティングに基づく機械学習モジュール10は、リアルタイム電位図信号ごとに、最大6つの異なる特徴を使用することができる。言うまでもなく、他の実施形態において、抽出器8は、より多くの特徴を抽出してもよい。これらの特徴は、同一のタイムコードを共有することによって、異なる電極に関連する電位図信号の組み合わせまたは比較に関連することができる。一部の態様において、本出願人は、勾配ブースティングに基づく機械学習モジュール10で最大70個の特徴を利用した。
勾配ブースティングに基づく機械学習モジュール10の訓練は、分散を示すかどうかを示す値でラベル付けされた275000個の電位図信号を用いて行われる。まず、これらの電位図信号ごとに特徴が抽出され、次に、勾配ブースティングに基づく機械学習モジュール10が、本明細書に記載の実施例では対数損失を用いて訓練される。訓練後、電位図信号のセットから抽出された特徴のセットを入力として受信するため、および電位図信号のセットのうちの対応する1つの電位図信号が分散を示す確率をそれぞれ示す値の配列を返すために、勾配ブースティングに基づく機械学習モジュール10が使用され得る。任意選択で、訓練データは、カテーテル形状情報および/またはそれに関連する3D情報を含んでもよい。
本明細書に記載の実施例において、第2の評価器12は、畳み込みニューラルネットワークである。畳み込みニューラルネットワークは、電位図信号のセットを受信し、電位図信号のセットのうちの対応する1つの電位図信号が分散を示す確率をそれぞれ示す値の配列を出力する。
図2は、第2の評価器12の畳み込みニューラルネットワークのための例示的な構造を示す。
本明細書に記載の実施例において、ニューラルネットワークは、5つの畳み込み層を有する畳み込みニューラルネットワークである。したがって、電位図信号のセット100(例えば16個の電位図信号を含む)は、200個の特徴を抽出する第1の畳み込み層110によって処理され、次に、層110から3000個の特徴を抽出する第2の畳み込み層120、層120から20000個の特徴を抽出する第3の畳み込み層130、層130の最大プーリングから8000個の特徴を抽出する第4の畳み込み層140、および層140から1000個の特徴を抽出する第5の畳み込み層150によって処理される。
第5の畳み込み層150は、ニューラルネットワークの完全接続層160にリンクされている。完全接続層160は、それぞれ50ニューロンおよび10ニューロンを含む2層のニューロンの連鎖を含む。完全接続層160は、出力層に確率配列170を返す。
この畳み込みニューラルネットワークの訓練は、分散を示すかどうかを示す値でラベル付けされた275000個の電位図信号を用いて行われる。電位図信号は、16個ずつのグループごとに使用される。各電位図信号は、分散の有無を示す0または1の対応する値を有する。本明細書に記載の実施例において、訓練は、アダムオプティマイザとバイナリクロスエントロピー損失を用いて行われる。他のオプティマイザと損失を使用してもよく、オプティマイザをバイパスしてもよい。
予測器14は、第1の評価器6の予測値と第2の評価器12の予測値を整合させるために使用される。より詳細には、多くの場合、第2の評価器12は、第1の評価器6よりも過学習に対してより堅牢であり、その結果、第1の評価器6は、通常よりも正確である。そのため、予測器14は、第1の評価器6および第2の評価器12によってそれぞれ出力された確率配列の値の加重平均を実行することができる。一実施形態において、重みは、第1の評価器6に対して0.7、第2の評価器12に対して0.3とすることができる。畳み込みニューラルネットワークが勾配ブースティングに基づく機械学習手法よりもドリフトしやすいので、第1の評価器6と第2の評価器12による予測値の差が大きすぎる場合、予測器14は、任意選択で、第1の評価器6の値のみを保持するように構成されてもよい。したがって、予測器14の出力は、所与のリアルタイム電位図信号が分散を示す別の確率配列である。
さらに、コンピュータ14は、予測配列に基づいて色を決定してもよい。より詳細には、この決定は、連続するタイムコードに関連する予測配列に基づいて行われ得る。患者の心房細動が遅いと判定された場合、連続するタイムコードに対応する2つの予測配列を平均化することができる。この配列の値は、次のように関連付けることができる:値が0.35未満の場合は青色、値が0.35以上、0.65以下の場合はオレンジ色、それ以外の場合は赤色。患者の心房細動が速いと判定された場合、連続するタイムコードに対応する4つの予測配列を平均化することができる。この配列の値は、次のように関連付けることができる:値が0.5未満の場合は青色、値が0.5以上、0.8以下の場合はオレンジ色、それ以外の場合は赤色。心房細動が速いか遅いかの判断は医師が行い、医師が装置に入力する。これにより、(図1の参照符号20で示すように)ディスプレイ上に視覚的な反応を提供することができる。これは、図1の参照符号22のドットで示す分散領域をもたらす心臓領域を自動的に、手動で、または半手動でタグ付けするために医師が使用することができる。
以上、2つの評価器を使用する場合について説明したが、本発明による第1の評価器および第2の評価器を含む2つ以上の評価器を使用する装置または方法も、添付の特許請求の範囲に含まれるものとする。
Claims (10)
- 電位図をリアルタイムで解析するためのコンピュータ装置であって、
複数の電極からそれぞれ発信される複数のリアルタイム電位図信号を受信するように配置されたメモリ(4)と、
抽出器(8)、および勾配ブースティングに基づく機械学習モジュール(10)を含む第1の評価器(6)であって、前記抽出器(8)は、前記電位図信号のセット内の各電位図信号から、少なくとも1つの時間的解析特徴および少なくとも1つの形態学的特徴を含む特徴のセットを抽出するように、および得られた前記特徴のセットを前記勾配ブースティングに基づく機械学習モジュール(10)に供給するように配置され、前記機械学習モジュール(10)は、関連する前記電位図信号が分散を示すかどうかを示す値でラベル付けされた特徴のセットを含むデータに基づいて訓練され、前記電位図信号のセットごとに、前記電位図信号のセットのそれぞれの電位図信号が分散を示すかどうかを示す第1の確率配列を出力するように配置される、第1の評価器(6)と、
畳み込みニューラルネットワークを含む第2の評価器(12)であって、前記畳み込みニューラルネットワークは、リアルタイム電位図信号のセットを受信し、入力されたリアルタイム電位図信号が分散を示す第2の確率配列を出力し、前記畳み込みニューラルネットワークは、それぞれの前記電位図信号が分散を示すかどうかを示す値でそれぞれラベル付けされた前記電位図信号のデータベースに基づいて訓練される、第2の評価器(12)と、
予測器(14)であって、前記電位図信号の所与のセットに対して決定された前記第1の確率配列および前記第2の確率配列に基づいて、前記第1の確率配列および前記第2の確率配列の値の加重平均に少なくとも部分的に基づく第3の確率配列を返す、予測器(14)と、
を備えるコンピュータ装置。 - 前記抽出器(8)は、第1の周期長推定値、第2の周期長推定値、および最大振幅を有する前記電位図信号の高速フーリエ変換内の周波数からなるグループの中から少なくとも1つの時間的解析特徴を抽出するように配置される、請求項1に記載のコンピュータ装置。
- 前記抽出器(8)は、前記電位図信号のユークリッドノルム、および前記電位図信号の積分された絶対導関数からなるグループの中から少なくとも1つの形態学的特徴を抽出するように配置される、請求項1または2に記載のコンピュータ装置。
- リアルタイム電位図信号を、選択された持続時間を有する一連の電位図信号に分割するように配置される、請求項1~3のいずれか1項に記載のコンピュータ装置。
- 同じ前記選択された持続時間を有する電位図信号のセットを提供するように配置される、請求項4に記載のコンピュータ装置。
- 前記予測器(14)は、同じ前記電位図信号に対応する、前記第1の確率配列中の確率と前記第2の確率配列中の確率との差の絶対値がしきい値を超えたときに、前記第3の確率配列において前記第1の確率配列中の確率を使用するように配置される、請求項1~5のいずれか1項に記載のコンピュータ装置。
- 前記コンピュータ装置は、前記第3の確率配列中の値に関連する色を決定するように配置され、
前記電極ごとに、前記第3の確率配列における対応する前記電位図信号について決定された確率に関連する色を出力するように配置されたディスプレイをさらに備える、
請求項1~6のいずれか1項に記載のコンピュータ装置。 - 請求項1~7のいずれか1項に記載の第1の評価器(6)、第2の評価器(12)、および予測器(14)を実現するための命令を含むコンピュータプログラム。
- 請求項8に記載のコンピュータプログラムを記録するデータ記憶媒体。
- リアルタイム電位図信号を受信し、請求項1~7のいずれか1項に記載の第1の評価器(6)、第2の評価器(12)、および予測器(14)を実行し、第3の確率配列を返す、コンピュータで実現される方法。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
FR2101234A FR3119537B1 (fr) | 2021-02-09 | 2021-02-09 | Dispositif informatique pour analyse en temps réel d’électrogrammes |
FR2101234 | 2021-02-09 | ||
PCT/EP2022/053048 WO2022171638A1 (en) | 2021-02-09 | 2022-02-08 | Computer device for real-time analysis of electrograms |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2024506633A true JP2024506633A (ja) | 2024-02-14 |
Family
ID=77710778
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2023548220A Pending JP2024506633A (ja) | 2021-02-09 | 2022-02-08 | 電位図をリアルタイム解析するコンピュータ装置 |
Country Status (7)
Country | Link |
---|---|
US (1) | US20240041379A1 (ja) |
EP (1) | EP4291094A1 (ja) |
JP (1) | JP2024506633A (ja) |
CN (1) | CN116981405A (ja) |
FR (1) | FR3119537B1 (ja) |
IL (1) | IL304609A (ja) |
WO (1) | WO2022171638A1 (ja) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR3143969A1 (fr) | 2022-12-27 | 2024-06-28 | Substrate Hd | Dispositif de traitement de signaux de dispersion spatiotemporelle d'électrogrammes cardiaques |
WO2024141730A1 (fr) | 2022-12-27 | 2024-07-04 | Substrate Hd | Dispositif de traitement de signaux de dispersion spatiotemporelle d'electrogrammes cardiaques |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP3558101B1 (en) * | 2016-12-21 | 2022-06-08 | Emory University | Methods and systems for determining abnormal cardiac activity |
US11389116B2 (en) * | 2018-01-19 | 2022-07-19 | Biosense Webster (Israel) Ltd. | Apparatus and method for heartbeat classification based on time sequence and morphology of intracardiac and body surface electrocardiogram (ECG) signals |
FR3079405B1 (fr) * | 2018-03-30 | 2023-10-27 | Substrate Hd | Dispositif informatique de detection de troubles du rythme cardiaque |
EP3616605A1 (en) * | 2018-08-30 | 2020-03-04 | Tata Consultancy Services Limited | Non-invasive detection of coronary heart disease from short single-lead ecg |
US12048487B2 (en) * | 2019-05-06 | 2024-07-30 | Biosense Webster (Israel) Ltd. | Systems and methods for improving cardiac ablation procedures |
-
2021
- 2021-02-09 FR FR2101234A patent/FR3119537B1/fr active Active
-
2022
- 2022-02-08 EP EP22704908.7A patent/EP4291094A1/en active Pending
- 2022-02-08 WO PCT/EP2022/053048 patent/WO2022171638A1/en active Application Filing
- 2022-02-08 JP JP2023548220A patent/JP2024506633A/ja active Pending
- 2022-02-08 CN CN202280014252.2A patent/CN116981405A/zh active Pending
- 2022-02-08 US US18/264,564 patent/US20240041379A1/en active Pending
-
2023
- 2023-07-20 IL IL304609A patent/IL304609A/en unknown
Also Published As
Publication number | Publication date |
---|---|
IL304609A (en) | 2023-09-01 |
FR3119537B1 (fr) | 2023-04-07 |
WO2022171638A1 (en) | 2022-08-18 |
FR3119537A1 (fr) | 2022-08-12 |
CN116981405A (zh) | 2023-10-31 |
US20240041379A1 (en) | 2024-02-08 |
EP4291094A1 (en) | 2023-12-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Yang et al. | Localization of origins of premature ventricular contraction by means of convolutional neural network from 12-lead ECG | |
US20190090774A1 (en) | System and method for localization of origins of cardiac arrhythmia using electrocardiography and neural networks | |
EP3744282A2 (en) | Systems for improving cardiac ablation procedures | |
CN106419896B (zh) | 匹配并跟踪心脏激活的时间序列 | |
JP2024506633A (ja) | 電位図をリアルタイム解析するコンピュータ装置 | |
US20220079499A1 (en) | Identification of ablation gaps | |
CN107007279A (zh) | 一种基于堆栈式自编码器的无创心内异常激动点定位方法 | |
Baumert et al. | Quantitative-electrogram-based methods for guiding catheter ablation in atrial fibrillation | |
CN109077720A (zh) | 信号处理方法、装置、设备和存储介质 | |
Liao et al. | Deep learning classification of unipolar electrograms in human atrial fibrillation: application in focal source mapping | |
Ganesan et al. | Computer-aided clinical decision support systems for atrial fibrillation | |
Yang et al. | Automated diagnosis of atrial fibrillation using ECG component-aware transformer | |
US20220005198A1 (en) | Automatic contiguity estimation of wide area circumferential ablation points | |
EP3278269B1 (fr) | Procédé et système d'identification d'un isthme dans une cartographie tridimensionnelle | |
US20210391082A1 (en) | Detecting atrial fibrillation and atrial fibrillation termination | |
Chen et al. | Disentangling inter-subject variations: Automatic localization of ventricular tachycardia origin from 12-lead electrocardiograms | |
US20200054233A1 (en) | Method for determining a plurality of action potentials in the heart | |
Shi et al. | Universal 12-lead ECG representation for signal denoising and cardiovascular disease detection by fusing generative and contrastive learning | |
Lee et al. | Using beat score maps with successive segmentation for ECG classification without R-peak detection | |
US20230146716A1 (en) | Digital twin of atria for atrial fibrillation patients | |
Abdi | Atrial fibrillation fingerprinting | |
Siekierski et al. | Heart beats classification method using a multi-signal ECG spectrogram and convolutional neural network with residual blocks | |
EP3960083A1 (en) | Separating abnormal heart activities into different classes | |
Potekhin et al. | ANALYSIS OF ACUTE MYOCARDIAL INFARCTION USING ECG SIGNALS AND MACHINE LEARNING ALGORITHMS | |
US10932686B2 (en) | Identifying activations in an atrial fibrillation electrogram |