[go: up one dir, main page]

JP7252109B2 - 診断支援システムおよび制御プログラム - Google Patents

診断支援システムおよび制御プログラム Download PDF

Info

Publication number
JP7252109B2
JP7252109B2 JP2019184841A JP2019184841A JP7252109B2 JP 7252109 B2 JP7252109 B2 JP 7252109B2 JP 2019184841 A JP2019184841 A JP 2019184841A JP 2019184841 A JP2019184841 A JP 2019184841A JP 7252109 B2 JP7252109 B2 JP 7252109B2
Authority
JP
Japan
Prior art keywords
atrial
subject
potential
unit
degree
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.)
Active
Application number
JP2019184841A
Other languages
English (en)
Other versions
JP2021058423A (ja
Inventor
嘉之 山海
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.)
Cyberdyne Inc
Original Assignee
Cyberdyne Inc
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 Cyberdyne Inc filed Critical Cyberdyne Inc
Priority to JP2019184841A priority Critical patent/JP7252109B2/ja
Priority to US17/764,103 priority patent/US12251230B2/en
Priority to EP20875276.6A priority patent/EP4042943A4/en
Priority to PCT/JP2020/036893 priority patent/WO2021070678A1/ja
Publication of JP2021058423A publication Critical patent/JP2021058423A/ja
Application granted granted Critical
Publication of JP7252109B2 publication Critical patent/JP7252109B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7235Details of waveform analysis
    • A61B5/7264Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems
    • A61B5/7267Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems involving training the classification device
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/24Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/316Modalities, i.e. specific diagnostic methods
    • A61B5/318Heart-related electrical modalities, e.g. electrocardiography [ECG]
    • A61B5/346Analysis of electrocardiograms
    • A61B5/349Detecting specific parameters of the electrocardiograph cycle
    • A61B5/361Detecting fibrillation
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/24Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/25Bioelectric electrodes therefor
    • A61B5/279Bioelectric electrodes therefor specially adapted for particular uses
    • A61B5/28Bioelectric electrodes therefor specially adapted for particular uses for electrocardiography [ECG]
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/24Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/316Modalities, i.e. specific diagnostic methods
    • A61B5/318Heart-related electrical modalities, e.g. electrocardiography [ECG]
    • A61B5/339Displays specially adapted therefor
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/24Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/316Modalities, i.e. specific diagnostic methods
    • A61B5/318Heart-related electrical modalities, e.g. electrocardiography [ECG]
    • A61B5/346Analysis of electrocardiograms
    • A61B5/349Detecting specific parameters of the electrocardiograph cycle
    • A61B5/352Detecting R peaks, e.g. for synchronising diagnostic apparatus; Estimating R-R interval
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/48Other medical applications
    • A61B5/4842Monitoring progression or stage of a disease
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16HHEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H40/00ICT specially adapted for the management or administration of healthcare resources or facilities; ICT specially adapted for the management or operation of medical equipment or devices
    • G16H40/60ICT specially adapted for the management or administration of healthcare resources or facilities; ICT specially adapted for the management or operation of medical equipment or devices for the operation of medical equipment or devices
    • G16H40/63ICT specially adapted for the management or administration of healthcare resources or facilities; ICT specially adapted for the management or operation of medical equipment or devices for the operation of medical equipment or devices for local operation
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16HHEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H50/00ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/20ICT 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
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16HHEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H50/00ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/30ICT 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
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16HHEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H50/00ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/70ICT 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)
  • Medical Informatics (AREA)
  • Public Health (AREA)
  • Biomedical Technology (AREA)
  • General Health & Medical Sciences (AREA)
  • Cardiology (AREA)
  • Pathology (AREA)
  • Physics & Mathematics (AREA)
  • Veterinary Medicine (AREA)
  • Surgery (AREA)
  • Molecular Biology (AREA)
  • Animal Behavior & Ethology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Biophysics (AREA)
  • Epidemiology (AREA)
  • Data Mining & Analysis (AREA)
  • Primary Health Care (AREA)
  • Artificial Intelligence (AREA)
  • Databases & Information Systems (AREA)
  • Evolutionary Computation (AREA)
  • Business, Economics & Management (AREA)
  • General Business, Economics & Management (AREA)
  • Fuzzy Systems (AREA)
  • Mathematical Physics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Physiology (AREA)
  • Psychiatry (AREA)
  • Signal Processing (AREA)
  • Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
  • Medical Treatment And Welfare Office Work (AREA)

Description

本発明は、診断支援システムおよび診断支援方法に関し、特に心臓の不整脈をモニタリングしながら持続性の心房細動を判定する診断支援システムおよび診断支援方法に適用して好適なものである。
心房細動は、我が国における死因の上位を占める心不全や脳梗塞を起こす重篤な疾患である。心房細動の有効な治療として、経皮的カテーテル心筋焼灼術がある。この経皮的カテーテル心筋焼灼術は心臓内部で計測した心電図を元に医師が焼灼部位を判定し、高周波通電を用いて、電気的に細胞を壊死させることにより心房細動を治療する手術である。
経皮的カテーテル心筋焼灼術の焼灼部位として、心房複雑電位(CFAEs:Complex Fractionated Atrial Electrograms)の計測される部位が重要な標的であることが従来から知られている。
心房複雑電位とは、120[ms]未満の短い間隔で発生する心房波と連続的に電位が発生する心房電位である連続電位という2つの特徴が1つ以上存在する心電図信号であり、心房細動を持続させることに寄与する。
従来から経皮的カテーテル心筋焼灼術における心房複雑電位の判定は、視覚により確認する手法以外にも、周波数解析や閾値判定による自動判定手法が提案されている。これらの周波数解析や閾値判定を用いた手法により、心房複雑電位の特徴の一つである短い間隔の心房波を判定することが可能となった。
しかし、連続的に電位が発生する心房電位である連続電位の計測される部位への焼灼も、経皮的カテーテル心筋焼灼術における焼灼標的として有効であるものの、自動判定を行う手法の確立には至っていないのが現状である。
従来から心房細動の検出方法としては、複数の脈波のうち全ての発生間隔を平均した平均発生間隔を中心とする予め設定された発生間隔判定範囲内にない発生間隔を有する不規則脈波を抽出し、複数の脈波の全ての発生間隔の数に対する不規則脈波の数の割合が予め設定された異常発生間隔発生割合判定値よりも大きい場合に、生体の心房細動が発生していると判定するようにした心房細動判定装置が提案されている(特許文献1参照)。
またユーザの心臓活動に関連する生体信号であるそれぞれの拍動と直前の拍動との間の拍動間隔を連続する拍動間隔シーケンスとして作成し、当該拍動間隔シーケンスを、連続する拍動間隔からなる既定の時間的長さの複数の有限のシーケンスセグメントに分割しておき、各シーケンスセグメントにおいて、拍動間隔シーケンスに関連するデータの分散を表す散布図を作成してシーケンスセグメントを暗号化する。
そして散布図の各特定の領域にプロットされたプロットポイントの量を表すヒストグラムを作成して当該散布図を暗号化し、散布図のプロットポイントの合計で各ビン(ヒストグラムに含まれる既定幅かつ既定数のビン)のビンカウントを除算してヒストグラムを正規化することにより、シーケンスセグメントが心房細動のリズムまたは非心房細動のリズムである確率を決定するようになされたユーザの心房細動の兆候を解析および同定するための方法も提案されている(特許文献2参照)。
さらに被検者の心臓の拍動間隔を測定して得られる拍動間隔データから当該被検者の心臓の拍動に由来しない不正な拍動間隔を除外した後、その除外後データを用いて心房細動が発生したか否かを判定して表示手段に表示することにより、被検者の心房細動の有無を検出する心房細動検出システムが提案されている(特許文献3参照)。
特開2018-153487号公報 特表2019-502437号公報 特開2019-129954号公報
ところが、特許文献1では、複数の心拍同期波の数に対する異常発生間隔心拍同期波の数の割合に応じて心房細動か否かを判定する手法であるが、予め設定した基準範囲を閾値とするため、判定精度の向上には限界がある。
また特許文献2では、連続する拍動間隔に基づくシーケンスセグメントが心房細動のリズムまたは非心房細動のリズムである確率を人工ニューラルネットワークを用いて決定するものであるが、医師による実際の判断結果を含むことなく、心房細動のリズム検出のアルゴリズムの精度向上に努めているにすぎない。
さらに特許文献3では、心臓の拍動間隔から不正な拍動間隔を除外して信号細動の有無を検出するにとどまり、実際の医師による判断結果とは無関係である。
実際に心房細動の症状については、患者本人が自覚している場合も少なく、医師が判断するには1ヶ月程度の期間分の拍動間隔を測定しておく必要があるため、患者の判断と医師の判断とが曖昧になるケースが多い。
本発明は以上の点を考慮してなされたもので、持続性の心房細動を医師が正確に判断できるように支援する診断支援システムおよび制御プログラムを提案するものである。
かかる課題を解決するため本発明においては、診断支援システムについて、対象者の心電図を計測する心電図計測部と、心電図計測部により計測された心電図データに基づいて、時系列的に隣り合うR波のピーク間隔である心拍間隔における心房細動の特徴を表す心房複雑電位を判定する心房複雑電位判定部と、心房複雑電位判定部による判定結果を対象者に固有の判定データとして時系列的に記憶するデータ記憶部と、データ記憶部から読み出した対象者における所定期間にわたって収集された判定データに基づいて、当該期間内の発生頻度に応じた頻度ヒストグラムを生成する頻度ヒストグラム生成部と、頻度ヒストグラム生成部による対象者における頻度ヒストグラムに基づいて、当該対象者が患う心房細動の進行度合いを判断するための判断学習を行う判断学習部とを備えるようにした。
この診断支援システムでは、対象者の心電図データから得られる所定期間分の心房複雑電位の判定結果に基づく発生頻度に応じた頻度ヒストグラムに基づいて、対象者が患う心房細動の進行度合いを判断して学習するようにしたことにより、持続性の心房細動を対象者自身が認識できない状況下でも医師が正確に診断することができるように支援することが可能となる。
また本発明においては、心房複雑電位判定部により判定される心房複雑電位は、心房波の発生頻度および連続的に発生する心房電位の組合せである。この結果、診断支援システムでは、経皮的カテーテル心筋焼灼術の焼灼部位を求めるための計測指標となる心房複雑電位(CFAEs)を自動的に判定することができる。
さらに本発明においては、頻度ヒストグラム生成部による対象者における頻度ヒストグラムを表示する表示部と、表示部に表示されている頻度ヒストグラムに対して、対応する対象者が患う心房細動の進行度合いを医師が診断結果として入力する入力部と、対象者について、判断学習部による判断学習の結果と、入力部による医師の診断結果とを参照し、心房細動の進行度合いを判断するための基準を示す基準関数を生成する基準関数生成部とを備えるようにした。
この結果、診断支援システムでは、医師が表示部に表示されている対象者の頻度ヒストグラムを目視確認しながら、所定期間にわたる対象者の心房細動の進行度合いを医師が入力部を介して入力することにより、当該対象者が患う心房細動の進行度合いを判断するための判断学習の精度を格段と向上させることができる。
さらに本発明においては、診断支援システムに、対象者の心電図を計測して得られる心電図データに基づいて、時系列的に隣り合うR波のピーク間隔である心拍間隔における心房細動の特徴を表す心房複雑電位を判定する第1ステップと、第1ステップによる判定結果を対象者に固有の判定データとして時系列的にデータ記憶部に記憶しておく第2ステップと、データ記憶部から読み出した対象者における所定期間にわたって収集された判定データに基づいて、当該期間内の発生頻度に応じた頻度ヒストグラムを生成する第3ステップと、第3ステップによる対象者における頻度ヒストグラムに基づいて、当該対象者が患う心房細動の進行度合いを判断するための判断学習を行う第4ステップとからなる一連の処理を制御プログラムに実行させるようにした。
この結果、診断支援システムの制御プログラムでは、対象者の心電図データから得られる所定期間分の心房複雑電位の判定結果に基づく発生頻度に応じた頻度ヒストグラムに基づいて、対象者が患う心房細動の進行度合いを判断して学習するようにしたことにより、持続性の心房細動を対象者自身が認識できない状況下でも医師が正確に診断し得るように支援することが可能となる。
以上のように本発明によれば、対象者に身体的負荷をかけることなく医師が心房細動の進行度合いを正確に診断し得るように支援することができる診断支援システムおよび制御プログラムを実現できる。
本発明の実施形態に係る診断支援システムの全体構成を示すブロック図である。 図1に示す情報管理装置における制御装置の構成を示すブロック図である。 持続性の心房細動を患う対象者の心電図を示すグラフである。 心房波の除去前後の対象者の心電図を示すグラフである。 心房波の判定結果に対する閾値超えの説明に供するグラフである。 心房波を除去した後の心房複雑電位を示すグラフである。 心房波を除去した後の2乗の電位や閾値の算出結果を示すグラフである。 心房複雑電位の判定処理の説明に供するフローチャートである。 発生頻度値の説明に供するグラフである。
以下図面について、本発明の一実施の形態を詳述する。
(1)本発明による診断支援システムの構成
図1は本実施の形態における診断支援システム1は、対象者の左胸部に非侵襲に装着される複数の電極からなる心電図センサ(心電図計測部)2と、当該対象者の手首に装着される通信端末部3とを有する。
心電図センサ2は、対象者の心筋の興奮による生じる活動電位を示す心電図を連続的に計測しながら、当該計測結果である心電図データを無線または有線にて腕時計型の通信端末部3に送信する。通信端末部3は、心電図データを内部メモリ(図示せず)に記憶するとともに、インターネット等の通信ネットワーク4を介して通信可能な情報管理装置10に送信する。
情報管理装置10は、装置全体を制御するCPU(Central Processing Unit)からなる制御装置11と、通信端末部か3ら通信ネットワーク4を介して通信インタフェース12にて受信される心電図データを記憶するデータ記憶部13とを有する。
制御装置11は、心房複雑電位判定部20と、頻度ヒストグラム生成部21と、判断学習部22とを有する。心房複雑電位判定部20は、通信インタフェース12を介して受信した心電図データに基づいて、時系列的に隣り合うR波のピーク間隔である心拍間隔における心房細動の特徴を表す心房複雑電位を判定し、当該判定結果を対象者に固有の判定データとしてデータ記憶部13に記憶する。
頻度ヒストグラム生成部21は、データ記憶部13から読み出した対象者における所定期間にわたって収集された判定データに基づいて、当該期間内の発生頻度に応じた頻度ヒストグラムを生成する。
判断学習部22は、頻度ヒストグラム生成部21による対象者における頻度ヒストグラムに基づいて、当該対象者が患う心房細動の進行度合いを判断するための判断学習を行う。
さらに診断支援システム1においては、情報管理装置10に加えて、医師が対象者の心房細動の進行度合いを診断するための医師側端末部30が通信ネットワーク4を介してワイヤレス接続されている。この医師側端末部30は、タッチパネル型の無線通信端末であり、通信インタフェース31を介して受信した各種情報に基づく映像等を表示する表示部32と、指先を用いた接触により入力操作が可能なタッチパネル式の入力部33と、本端末全体を統括制御する端末制御部34とを有する。
(2)制御装置における診断支援処理
ここで、心電図センサ2により測定した心電図を図2(A)~(C)に示す。心房波が120[ms]以上の間隔で計測される特徴がある、心房細動が生じていない心房複雑電位の心電図を図2(A)に示し、短い間隔の心房波と連続電位の二つの特徴が確認できる心房複雑電位の心電図をそれぞれ図2(B)および(C)に示す。
経皮的カテーテル心筋焼灼術において、心房複雑電位の判定を行う心房から計測される心電図は、心房電位と同時に心室波が計測される場合がある。この心室波は持続性心房細動に寄与しないため、体表で計測した心室波を用いて、心房で計測される心電図信号から心室波を除去する。
体表で計測される心室波の判定には、後述する心房波間隔判定アルゴリズムを用いる。除去する心室波の興奮時間は60[ms]から100[ms]である。
制御装置11において、心房複雑電位判定部20は、心房波間隔判定アルゴリズムによって判定した体表の各心室波に対し、最大値を求め、心房で計測した電位に対して、各心室波の最大値が判定された時間の前後30[ms]に計測された電位を0[mV]にする。
心房で計測された心電図と心室波除去後の心電図をそれぞれ図3(A)および(B)に示す。図3の(B)は、図3(A)から、心室波を除去した心房電位である。
心房複雑電位判定部20は、図3(B)のように、計測された心電図データから心室波除去処理を施した心電図データを用いて、心房複雑電位の特徴である短い間隔の心房波と連続電位を判定する。
心房複雑電位判定部20は、心房波間隔判定アルゴリズムを用いて、心房複雑電位の特徴の一つである短い間隔の心房波を判定する。具体的には、図2(B)に示すような、心房波の間隔が120[ms]未満である部分を判定する。
まず心房複雑電位判定部20は、心電図データ中の心房波と連続電位との差を強調し、電位を正値にするため、電位e(t)を2乗する。e(t)に対して心房波間隔の判定を行う。心房波の電圧値は、計測部位、心臓の筋肉量等により大きく変動するため、波形毎に閾値を生成する必要がある。そこで、2乗した各電位e(t)に対して、電位e(t)の分散を用いて閾値eth1を生成し、当該閾値eth1を次式(1)にて定義する。
Figure 0007252109000001
図2(B)に対する2乗した電位e(t)および閾値eth1を図4に示す。この心房波間隔判定アルゴリズムを用いて、200[ms]間隔で2乗した電位e(t)が閾値eth1を超える回数が求められる。図2(B)に対して、閾値eth1を超える回数を求めた結果が図5に示される。
ここで、閾値を超える回数が2回以上のとき、心房波の間隔が120[ms]未満となり、短い間隔の心房波と判定される。図5の網掛け部に示すような閾値を超える回数が2回以上である部分の閾値を超えた回数を積分した値Iを算出する。値Iの値が大きい程、心房複雑電位の特徴である短い間隔の心房波が存在することを示す。
心房波の最小値を決定するパラメータKは、3人の持続性の心房細動患者(男性2名、女性1名、平均年齢59±6.4)から計測した計30の心電図を用いた判定結果から設定した。医師が判定した心房波に対して、この心房波間隔判定アルゴリズムにより検出した心房波の正診率が100[%]となるパラメータK=3に設定した。
続いて心房複雑電位判定部20は、連続電位判定アルゴリズムを用いて、図2(B)および(C)に示す連続電位を判定する。連続電位は、連続的に計測される心房電位であるが、ここでは振れ幅の大きい心房電位(心房波)とは区別する必要がある。そのため、連続電位判定アルゴリズムにより判定された心房波を電位から除去し、当該アルゴリズムで心房波を判定することを避ける。図6に連続電位判定アルゴリズムにより判定された心房波を除去した電位e(t)を実線、除去された心房波を点線でそれぞれ示す。
連続電位は連続した振動する心房電位であり、その電圧値は、計測部位、心臓の筋肉量等により変動する。そこで心房波の最大値から閾値を生成し、閾値を超えた部分を連続電位として判定する。心房波の最大値emaximumから、連続電位を判定するための閾値eth2を次式(2)にて定義する。
Figure 0007252109000002
設定された閾値に対して、閾値を超えた回数の積分値を値Iとして計算する。値Iは連続電位の継続時間を表しており、値Iが大きい程、電位が持続的で安定した心房複雑電位であることを示す。
連続電位の最小値を決定するパラメータKは、パラメータKと同様の計30の心電図を用いた判定から設定した。心電図に対して、医師が連続電位の見られる電位と連続電位の見られない電位に判定した結果と連続電位判定アルゴリズムにより判定された連続電位の正診率が最も高い93[%]となるK=1/40に設定した。図6の心房波を除去した電位e(t)に対して2乗の電位e(t)と閾値eth2を求めた結果を図7に示す。
心房複雑電位判定部20は、心房波間隔判定アルゴリズムにより求めた値Iと、連続電位判定アルゴリズムにより求めた値Iを用いて、心房複雑電位の特徴である短い間隔の心房波と連続電位の2項目から心房複雑電位を判定する。値Iは短い間隔の心房波の発生頻度を反映しており、値Iは連続電位の持続時間を反映している。
続いて心房複雑電位判定部20は、心房複雑電位の特徴である短い間隔の心房波と連続電位の2項目を反映した値I、値Iと計測時間tから、心房複雑電位の二つの特徴が発生する頻度を表す値(発生頻度値)Qを次式(3)により算出する。
Figure 0007252109000003
本手法では、短い間隔の心房波と連続電位を反映した値Qの値に閾値QThresholdを設定し、心房複雑電位と非心房複雑電位の自動判定を行う。Q≦QThresholdの場合、その電位を非心房複雑電位とし、Q>QThresholdの場合、心房複雑電位として判定する。
非心房複雑電位および心房複雑電位の閾値QThresholdは、パラメータK、Kの設定で対象とした3人とは別の持続性心房細動患者7人(男性6名、女性1名、平均年齢56±7.8)から計測した53の電位(心房複雑電位:20、非心房複雑電位:33)に対して値Qを求め、心房複雑電位と非心房複雑電位の判定の正診率が最も高い94[%]となる閾値QThreshold=4と設定した。53の電位(心房複雑電位:20、非心房複雑電位:33)に対し、本手法により計算された値Qの結果を図9(A)に示す。
上述した心房複雑電位判定部20による心房複雑電位または非心房複雑電位を判定するフローチャートを図8に示す。情報管理装置10の制御装置11における心房複雑電位判定部20は、対象者が装着する通信端末部3から心電図データを受け取ると、図8に示す心房複雑電位判定処理手順RT1をステップSP0から開始する。
心房複雑電位判定部20は、対象者の表面心電図(体表で計測した心電図)を用いた心房波間隔判定アルゴリズムにより心房波を検出した後(ステップSP1)、心房波をフィルタリングにより除去する(ステップSP2)。
続いて心房複雑電位判定部20は、心房波間隔判定アルゴリズムにより心房波の発生頻度を反映する値Iを算出した後(ステップSP3)、連続電位判定アルゴリズムにより連続電位の持続時間を反映する値Iを算出する(ステップSP4)。
そして心房複雑電位判定部20は、値Iおよび値Iを用いて発生頻度値Qを算出し(ステップSP5)、当該発生頻度値Qを自動判定基準となる所定の閾値QThresholdと比較する(ステップSP6)。
心房複雑電位判定部20は、このステップSP6において、発生頻度値Qが閾値より大きい場合には、心房複雑電位として判定する一方(ステップSP7)発生頻度値Qが閾値以下の場合には、非心房複雑電位として判定し(ステップSP8)、当該処理手順RT1を終了する(ステップSP9)。
実際に以下に示すような実験を実行して本発明手法の有効性を検証する。すなわち、パラメータK、K、閾値QThresholdの設定で対象とした10人とは別の8人の持続性心房細動患者(男性6名、女性2名、平均年齢60±10)の心臓内部から計測された心電図に対する自動判定を行う実験を通して、本発明の有効性を検証した。
心電図は4秒から32秒間計測された計100の心電図を用いた。それぞれの電位に対して、心房波の間隔が120[ms]未満である、連続電位が存在するという心房複雑電位の定義に従い、医師が心房複雑電位と非心房複雑電位の判定を行った。その結果、計100の電位は心房複雑電位が62、非心房複雑電位が38であった。この100の電位を用いて、心房複雑電位と非心房複雑電位の自動判定を行った。
心房複雑電位の特徴である短い間隔の心房波と連続電位の2項目から心房複雑電位の判定を行う本手法を用いて、心房複雑電位と非心房複雑電位の自動判定行い、自動判定結果の正診率から、心房細動における心房複雑電位の自動判定に対する本発明手法の有効性を検証した。
従来技術である高速フーリエ変換(FFT:Fast Fourier Transform)による判定手法によって解析を行った。20[Hz]以下でパワースペクトルが最大となる周波数DF(Dominant Frequency)について、心房複雑電位の特徴である心房波の間隔が120[ms]未満であることから8.3[Hz]をDFの閾値とする。DFが8.3[Hz]未満である電位を非心房複雑電位、DFが8.3[Hz]以上である電位を心房複雑電位として判定を行った。本提案手法と従来手法の正診率から有効性を検証する。
本発明手法を用いて、心房複雑電位と非心房複雑電位を自動判定した結果を示す。8人の持続性心房細動患者の心房から計測された計100の心電図信号(心房複雑電位:62、非心房複雑電位:38)に対して、本手法により算出したQの値を図9(B)に示す。これにより、計100の電位(心房複雑電位:62、非心房複雑電位:38)に対して、陽性適中率(positive predictive value)は97[%]、陰性適中率(negative predictive value)は88[%]、正診率(accuracy)は93[%]で判定することができた。
同様の8人の持続性心房細動患者の心房から計測された計100の電位(心房複雑電位:62、非心房複雑電位:38)に対して、FFTを用いた解析手法によりDFを算出した結果を図9(C)に示す。FFTを用いた解析手法により算出したDFから自動判定した結果、計100の電位(心房複雑電位:62、非心房複雑電位:38)に対して、陽性適中率は71[%]、陰性適中率は42[%]、正診率は50[%]で判定することができた。
本発明手法と従来手法を比較すると、陽性適中率、陰性適中率、正診率のすべてにおいて、本発明手法はより高い精度で心房複雑電位と非心房複雑電位を自動判定することができた。
また図9(C)に示すように、従来手法の心房複雑電位に対する自動判定の適中率を示す感度は32[%]であったのに対して、図9(B)に示すように、提案手法の心房複雑電位に対する自動判定の結果は感度92[%]であった。この結果は、心房波間隔を判定する従来手法と比較して、短い間隔の心房波と連続電位の2項目から心房複雑電位を判定する本発明手法が、従来手法で定義された心房複雑電位より高い精度で判定可能であることを示す。
(3)本実施の形態による心房細動の進行度合いの判断学習
上述のように心房複雑電位判定部20は、心房複雑電位および非心房複雑電位を判定し、当該判定結果を表す判定データをデータ記憶部13に記憶する。その後、頻度ヒストグラム生成部21は、データ記憶部13から読み出した対象者における所定期間にわたって収集された判定データに基づいて、当該期間内の発生頻度に応じた頻度ヒストグラムを生成する。
判断学習部22は、頻度ヒストグラム生成部21による対象者における頻度ヒストグラムに基づいて、当該対象者が患う心房細動の進行度合いを判断するための判断学習を行う。すなわち、経皮的カテーテル心筋焼灼術の焼灼部位を求めるための計測指標となる心房複雑電位の判定結果を所定期間にわたって収集して得られる頻度ヒストグラムに基づいて、判断学習部22は、対象者の心房細動の進行度合いを判断して学習する。
この結果、持続性の心房細動を対象者自身が認識できない状況下でも医師が正確に診断することができるように支援することが可能となる。
(4)医師による診断方法
診断支援システム1においては、医師が医師側端末部30を用いて、対象者の心房細動の進行度合いを診断する。医師側端末部30では、表示部32に頻度ヒストグラム生成部21による対象者における頻度ヒストグラムが表示された状態で、医師が当該頻度ヒストグラムを目視確認しながら、入力部33を用いて、対象者が患う心房細動の進行度合いを診断結果として入力する。
端末制御部(基準関数生成部)34は、対象者について、判断学習部22による判断学習の結果と、入力部33による医師の診断結果とを参照し、心房細動の進行度合いを判断するための基準を示す基準関数を生成する。
この結果、診断支援システム1では、医師が表示部32に表示されている対象者の頻度ヒストグラムを目視確認しながら、所定期間にわたる対象者の心房細動の進行度合いを医師が入力部33を介して入力することにより、当該対象者が患う心房細動の進行度合いを判断するための判断学習の精度を格段と向上させることができる。
(5)他の実施の形態
なお上述のように本実施の形態においては、診断支援システム1を、対象者の胸部に心電図センサ2を非接触で装着するとともに当該対象者の手首に通信端末部3を装着するようにした場合について述べたが、本発明はこれに限らず、対象者の体表からの心電図を非侵襲で計測することができれば、これ以外の種々の構成に広く適用するようにしてもよい。ただし、所定期間(4週間)程度の間継続して計測することができるように、対象者への身体的負荷が極力かからない構成にするのが望ましい。
また本実施の形態においては、頻度ヒストグラム生成部21は、頻度ヒストグラムを生成するための対象者における判定データを収集する所定期間として、4週間を設定した場合について述べたが、本発明はこれに限らず、1週間以上であれば医師との相談の上で最適なデータ収集期間を所定期間として設定するようにしてもよい。
1…診断支援システム、2…心電図センサ(心電図計測部)、3…通信端末部、4…通信ネットワーク、10…情報管理装置、11…制御装置、12、31…通信インタフェース、13…データ記憶部、20…心房複雑電位判定部、21…頻度ヒストグラム生成部、22…判断学習部、30…医師側端末部、32…表示部、33…入力部、34…端末制御部、RT1…心房複雑電位判定処理手順。

Claims (10)

  1. 対象者の心電図を計測する心電図計測部と、
    前記心電図計測部により計測された心電図データに基づいて、時系列的に隣り合うR波のピーク間隔である心拍間隔における心房細動の特徴を表す心房複雑電位を判定する心房複雑電位判定部と、
    前記心房複雑電位判定部による判定結果を前記対象者に固有の判定データとして時系列的に記憶するデータ記憶部と、
    前記データ記憶部から読み出した前記対象者における所定期間にわたって収集された前記判定データに基づいて、当該期間内の発生頻度に応じた頻度ヒストグラムを生成する頻度ヒストグラム生成部と、
    前記頻度ヒストグラム生成部による前記対象者における前記頻度ヒストグラムに基づいて、当該対象者が患う心房細動の進行度合いを判断するための判断学習を行う判断学習部と
    を備えることを特徴とする診断支援システム。
  2. 前記心房複雑電位判定部により判定される心房複雑電位は、心房波の発生頻度および連続的に発生する心房電位の組合せである
    ことを特徴とする請求項1に記載の診断支援システム。
  3. 前記頻度ヒストグラム生成部による前記対象者における前記頻度ヒストグラムを表示する表示部と、
    前記表示部に表示されている前記頻度ヒストグラムに対して、対応する前記対象者が患う心房細動の進行度合いを医師が診断結果として入力する入力部と、
    前記対象者について、前記判断学習部による判断学習の結果と、前記入力部による前記医師の診断結果とを参照し、前記心房細動の進行度合いを判断するための基準を示す基準関数を生成する基準関数生成部と
    を備えることを特徴とする請求項1または2に記載の診断支援システム。
  4. 前記基準関数生成部により生成された前記基準関数により、前記頻度ヒストグラム生成部による前記対象者における前記頻度ヒストグラムから当該対象者が患う心房細動の進行度合いを判断する進行度合判断部
    を備えることを特徴とする請求項3に記載の診断支援システム。
  5. 前記心房細動の進行度合いは、前記心電図データに基づく心電図波形の基線が細かく動揺しており、かつP波が存在せず、前記心拍間隔が不規則である状態が進行する度合いである
    ことを特徴とする請求項3または4に記載の診断支援システム。
  6. 診断支援システムに、
    対象者の心電図を計測して得られる心電図データに基づいて、時系列的に隣り合うR波のピーク間隔である心拍間隔における心房細動の特徴を表す心房複雑電位を判定する第1ステップと、
    前記第1ステップによる判定結果を前記対象者に固有の判定データとして時系列的にデータ記憶部に記憶しておく第2ステップと、
    前記データ記憶部から読み出した前記対象者における所定期間にわたって収集された前記判定データに基づいて、当該期間内の発生頻度に応じた頻度ヒストグラムを生成する第3ステップと、
    前記第3ステップによる前記対象者における前記頻度ヒストグラムに基づいて、当該対象者が患う心房細動の進行度合いを判断するための判断学習を行う第4ステップと
    からなる一連の処理を実行させるための制御プログラム
  7. 前記第1ステップにより判定される心房複雑電位は、心房波の発生頻度および連続的に発生する心房電位の組合せである
    ことを特徴とする請求項6に記載の制御プログラム
  8. 請求項6または7に記載の診断支援システムに、
    前記第3ステップによる前記対象者における前記頻度ヒストグラムを表示部に表示しておき、当該表示部に表示されている前記頻度ヒストグラムに対して、対応する前記対象者が患う心房細動の進行度合いが診断結果として外部入力された際、前記対象者について、前記第4ステップによる判断学習の結果と、当該診断結果とを参照し、前記心房細動の進行度合いを判断するための基準を示す基準関数を生成する第5ステップ
    からなる処理を続けて実行させることを特徴とする制御プログラム
  9. 請求項8に記載の診断支援システムに、
    前記第5ステップにより生成された前記基準関数により、前記第3ステップによる前記対象者における前記頻度ヒストグラムから当該対象者が患う心房細動の進行度合いを判断する第6ステップ
    からなる処理を続けて実行させることを特徴とする制御プログラム
  10. 前記心房細動の進行度合いは、前記心電図データに基づく心電図波形の基線が細かく動揺しており、かつP波が存在せず、前記心拍間隔が不規則である状態が進行する度合いである
    ことを特徴とする請求項6から9までのいずれか一項に記載の制御プログラム
JP2019184841A 2019-10-07 2019-10-07 診断支援システムおよび制御プログラム Active JP7252109B2 (ja)

Priority Applications (4)

Application Number Priority Date Filing Date Title
JP2019184841A JP7252109B2 (ja) 2019-10-07 2019-10-07 診断支援システムおよび制御プログラム
US17/764,103 US12251230B2 (en) 2019-10-07 2020-09-29 Diagnosis support system and diagnosis support method
EP20875276.6A EP4042943A4 (en) 2019-10-07 2020-09-29 Diagnosis assistance device and diagnosis assistance method
PCT/JP2020/036893 WO2021070678A1 (ja) 2019-10-07 2020-09-29 診断支援システムおよび診断支援方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2019184841A JP7252109B2 (ja) 2019-10-07 2019-10-07 診断支援システムおよび制御プログラム

Publications (2)

Publication Number Publication Date
JP2021058423A JP2021058423A (ja) 2021-04-15
JP7252109B2 true JP7252109B2 (ja) 2023-04-04

Family

ID=75380987

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2019184841A Active JP7252109B2 (ja) 2019-10-07 2019-10-07 診断支援システムおよび制御プログラム

Country Status (4)

Country Link
US (1) US12251230B2 (ja)
EP (1) EP4042943A4 (ja)
JP (1) JP7252109B2 (ja)
WO (1) WO2021070678A1 (ja)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2017006230A (ja) 2015-06-18 2017-01-12 セイコーエプソン株式会社 生体情報測定装置
US20180052957A1 (en) 2016-08-16 2018-02-22 Braemar Manufacturing, Llc End of service summary report for mobile cardiac outpatient telemetry
JP2019502437A (ja) 2015-12-07 2019-01-31 スマート ソリューションズ テクノロジーズ,エス.エル. 心房細動検出システムおよび使用方法
JP2019129954A (ja) 2018-01-30 2019-08-08 ユニオンツール株式会社 心房細動検出システム

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9456759B2 (en) * 2013-09-16 2016-10-04 Biotronik Se & Co. Kg Device for automatic mapping of complex fractionated atrial electrogram
JP6874258B2 (ja) 2017-03-17 2021-05-19 株式会社エー・アンド・デイ 生体の心房細動判定装置

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2017006230A (ja) 2015-06-18 2017-01-12 セイコーエプソン株式会社 生体情報測定装置
JP2019502437A (ja) 2015-12-07 2019-01-31 スマート ソリューションズ テクノロジーズ,エス.エル. 心房細動検出システムおよび使用方法
US20180052957A1 (en) 2016-08-16 2018-02-22 Braemar Manufacturing, Llc End of service summary report for mobile cardiac outpatient telemetry
JP2019129954A (ja) 2018-01-30 2019-08-08 ユニオンツール株式会社 心房細動検出システム

Also Published As

Publication number Publication date
US12251230B2 (en) 2025-03-18
EP4042943A4 (en) 2024-02-28
WO2021070678A1 (ja) 2021-04-15
US20220378355A1 (en) 2022-12-01
JP2021058423A (ja) 2021-04-15
EP4042943A1 (en) 2022-08-17

Similar Documents

Publication Publication Date Title
US8666483B2 (en) System for cardiac medical condition detection and characterization
CN107582042B (zh) 多个电解剖图的自动创建
US8706201B2 (en) Apparatus and method for analysis of high frequency QRS complexes
US8929975B2 (en) System for heart monitoring, characterization and abnormality detection
US8380305B2 (en) System and method for predicting successful defibrillation for ventricular fibrillation cardiac arrest
US20090216144A1 (en) Hopping methods for the detection of QRS onset and offset
JP5996111B2 (ja) 電流密度変動に基づく心筋損傷の評価装置
US8684942B2 (en) System for cardiac impairment detection based on signal regularity
EP3139828B1 (en) Detection of pulmonary vein isolation
US8868168B2 (en) System for cardiac condition characterization using electrophysiological signal data
EP3721794B1 (en) Vcg vector loop bifurcation
US20150126884A1 (en) System for Heart Performance Characterization and Abnormality Detection
US11771360B2 (en) Cardiac monitoring system with normally conducted QRS complex identification
JP2020517334A5 (ja)
EP3845127B1 (en) Methods and systems for estimation of residual ecg noise level and adaptive noise threshold
JP7252109B2 (ja) 診断支援システムおよび制御プログラム
CN118415615A (zh) 一种创伤患者远程护理监测方法及系统
CN111096740B (zh) 心电信号分析方法及装置、信号记录仪和三维标测系统
EP4091547B1 (en) Medical apparatus for diagnostic and site determination of cardiac arrhythmias
Ittatirut et al. Detection of premature ventricular contraction for real-time applications
JP6173905B2 (ja) Twa測定装置及びtwa測定装置の作動方法
JP2021094381A (ja) ユニポーラ信号及びバイポーラ信号を使用した心室活動の検出
Khriji et al. New adaptive thresholding-based ECG R-peak detection technique
KR20150081763A (ko) 심전도 신호의 저전력 고효율 r파 검출 방법 및 시스템
JP2536410B2 (ja) 心電図解析装置

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20220727

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20230214

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20230306

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20230314

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20230323

R150 Certificate of patent or registration of utility model

Ref document number: 7252109

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150