ホーム > 特許ランキング > 日本電信電話株式会社
(19)【発行国】日本国特許庁(JP)
(12)【公報種別】特許公報(B2)
(11)【特許番号】
(24)【登録日】2023-12-19
(45)【発行日】2023-12-27
(54)【発明の名称】心起電力推定方法、心起電力推定装置及びプログラム
(51)【国際特許分類】
A61B 5/346 20210101AFI20231220BHJP
A61B 5/353 20210101ALI20231220BHJP
A61B 5/355 20210101ALI20231220BHJP
A61B 5/352 20210101ALI20231220BHJP
【ＦＩ】
A61B5/346
A61B5/353
A61B5/355
A61B5/352
【請求項の数】
10
(21)【出願番号】P 2022507004
(86)(22)【出願日】2020-03-09
(86)【国際出願番号】 JP2020010025
(87)【国際公開番号】W WO2021181465
(87)【国際公開日】2021-09-16
【審査請求日】2022-06-24
(73)【特許権者】
【識別番号】000004226
【氏名又は名称】日本電信電話株式会社
(74)【代理人】
【識別番号】110001634
【氏名又は名称】弁理士法人志賀国際特許事務所
(72)【発明者】
【氏名】塚田 信吾
【審査官】永田 浩司
(56)【参考文献】
【文献】特開２０１７－１６４１６４（ＪＰ，Ａ）
【文献】真島三郎，心電図学の基礎に関する諸問題，心電図，1981年，Vol.1 No. 2，117-125
【文献】堤健，1.心室内各部の活動電位持続時間の差異とＴ波，JPN. J. ELECTROCARDIOLOGY，1989年，Vol.9 No.1，3-16
【文献】真島三郎，Ventricular Gradient の概念とその応用，心臓，1969年，1巻10号，1003-1012
(58)【調査した分野】(Int.Cl.，ＤＢ名)
Ａ６１Ｂ ５／２４
(57)【特許請求の範囲】
【請求項1】
心起電力の時系列の推定対象の心電位の時系列を、心筋の１回の拍動によって生じる心電位を表す時系列が互いに分離されるように区分する区分ステップと、前記区分ステップにおいて区分された各時系列を単位心電位時系列として、単位心電位時系列ごとに単位心電位時系列の各データを脱分極期に属するデータと不応期に属するデータと再分極期に属するデータと休止期に属するデータとに分類することで、前記単位心電位時系列ごとに、脱分極期に属するデータの集合である脱分極心電位時系列と不応期に属するデータの集合である不応心電位時系列と再分極期に属するデータの集合である再分極心電位時系列とを取得する状態時系列取得ステップと、
前記脱分極心電位時系列に基づき前記脱分極期の各時点における心起電力を示す脱分極心起電力時系列を推定する脱分極心起電力推定ステップと、
前記不応心電位時系列と前記脱分極期の終了時点における心起電力の推定値とに基づき前記不応期の各時点における心起電力を示す不応心起電力時系列を推定する不応心起電力推定ステップと、
前記再分極心電位時系列と前記不応期の終了時点における心起電力の推定値とに基づき前記再分極期の各時点における心起電力を示す再分極心起電力時系列を推定する再分極心起電力推定ステップと、
を有し、
前記心電位の時系列は心起電力によって生じる前記推定対象の体表面の心電位の時系列である、
心起電力推定方法。
【請求項2】
前記脱分極心起電力推定ステップにおいて前記脱分極心起電力時系列に含まれる心起電力を推定する場合、前記脱分極期の開始時点から、推定される心起電力に対応する時点まで、の隣接する心電位の差の絶対値の総和を前記開始時点における心電位に足し算した値を心起電力として取得する、請求項１に記載の心起電力推定方法。
【請求項3】
前記脱分極心起電力推定ステップでは、予め取得済みの脱分極心電位時系列と脱分極心起電力時系列との間の対応関係を示す全単射の写像を用いて前記推定対象の脱分極心起電力時系列を推定する、請求項１に記載の心起電力推定方法。
【請求項4】
前記再分極心起電力推定ステップにおいて前記再分極心起電力時系列に含まれる心起電力を推定する場合、前記再分極心電位時系列における前記不応期の終了時点から、推定される心起電力に対応する時点まで、の隣接する心電位の差に基づく値の総和を前記終了時点における心起電力の推定値から引き算した値を心起電力として取得する、請求項１から３のいずれか一項に記載の心起電力推定方法。
【請求項5】
前記再分極心起電力推定ステップでは、予め取得済みの再分極心電位時系列と再分極心起電力時系列との間の対応関係を示す全単射の写像を用いて前記推定対象の再分極心起電力時系列を推定する、請求項１から３のいずれか一項に記載の心起電力推定方法。
【請求項6】
前記不応心起電力推定ステップにおいて前記不応心起電力時系列に含まれる心起電力を推定する場合、前記脱分極期の終了時点における心起電力の推定値から前記不応期の対応する時点における心電位に関する値を引き算した値を心起電力として取得する、請求項１から５のいずれか一項に記載の心起電力推定方法。
【請求項7】
前記不応心起電力推定ステップにおいて推定する各時点の心起電力は、予め機械学習によって不応心電位時系列と不応心起電力時系列との関係を学習済みの学習済みモデルを用いて推定された値である、請求項１から５のいずれか一項に記載の心起電力推定方法。
【請求項8】
前記不応心起電力推定ステップにおいて前記不応心起電力時系列に含まれる心起電力を推定する場合、前記不応心電位時系列が示す心電位に定数を乗算した値を前記脱分極期の終了時点における心起電力の推定値の値に足し算した値を心起電力として取得する、請求項１から５のいずれか一項に記載の心起電力推定方法。
【請求項9】
心起電力の時系列の推定対象の心電位の時系列である心電位時系列を取得する取得部と、前記心電位時系列を、心筋の１回の拍動によって生じる心電位を表す時系列が互いに分離されるように区分する区分処理と、前記区分処理によって区分された各時系列を単位心電位時系列として、単位心電位時系列ごとに単位心電位時系列の各データを脱分極期に属するデータと不応期に属するデータと再分極期に属するデータと休止期に属するデータとに分類することで、前記単位心電位時系列ごとに、脱分極期に属するデータの集合である脱分極心電位時系列と不応期に属するデータの集合である不応心電位時系列と再分極期に属するデータの集合である再分極心電位時系列とを取得する状態時系列取得処理と、前記脱分極心電位時系列に基づき前記脱分極期の各時点における心起電力を示す脱分極心起電力時系列を推定する脱分極心起電力推定処理と、前記不応心電位時系列と前記脱分極期の終了時点における心起電力の推定値とに基づき前記不応期の各時点における心起電力を示す不応心起電力時系列を推定する不応心起電力推定処理と、前記再分極心電位時系列と前記不応期の終了時点における心起電力の推定値とに基づき前記再分極期の各時点における心起電力を示す再分極心起電力時系列を推定する再分極心起電力推定処理と、を実行する制御部と、
を備え、
前記心電位の時系列は心起電力によって生じる前記推定対象の体表面の心電位の時系列である、
心起電力推定装置。
【請求項10】
請求項９に記載の心起電力推定装置としてコンピュータを機能させるためのプログラム。
【発明の詳細な説明】
【技術分野】
【0001】
本発明は、心起電力推定方法、心起電力推定装置及びプログラムに関する。
【背景技術】
【0002】
心筋梗塞や狭心症に代表される心筋虚血による心筋の異常は心筋の膜電位の集合である心起電力の異常として現れることが知られている。また心筋や刺激伝導系の障害、心不全などの心臓への負荷も心起電力に影響を及ぼす。しかしながら、心起電力を直接測定することは難しい。そこで、医療現場では多くの場合心起電力によって生じる体表面の心電位の時系列である心電図を見て心筋の異常を推定する。具体的には、心電位の時系列のＰ、ＱＲＳ、ＳＴ、Ｔの各成分部とその間隔の変化を見て心筋の異常を推定する。しかし原信号の心起電力と比較して体表面の心電位は直流成分を中心に著しく減衰するため、心筋の異常によって生じる変化は小さく、そのため異常の検出が困難な場合がある。そこで、心電位の時系列を再現する生物物理、数理モデルを用いたシミュレーションの実行により心起電力を含む心機能を推定する方法が提案されている。
【先行技術文献】
【非特許文献】
【0003】
【文献】大江 透 「心臓電気生理学的検査（単相性心筋活動電位記録）」 日本循環器学会専門医誌 循環器専門医第２８巻 ｐ．１２９－１３５ ２０１９年８月
【文献】K.H.W.J.ten Tusscher, D.Noble, P.J.Noble, and A.V.Panfilov, “A model for human ventricular tissue”, Am J Physiol Heart Circ Physiol. 286:H1573-H1589, 2004
【文献】田中 博 「心電図逆問題における方法論」 医用電子と生体工学 １９８５年 第２３巻 第３号 ｐ．１４７－１５８
【文献】MICHAEL R. FRANZ, “Long-Term Recording of Monophasic Action Potentials From Human Endocardium” Am J Cardiol. 1983 Jun; 51, p.1629-1634
【発明の概要】
【発明が解決しようとする課題】
【0004】
しかしながら、生物物理、数理モデルを用いた非線形有限要素法等のシミュレーションは計算負荷が高く、実測値と整合する際に経験と高度な技術が求められるという問題がある。そのため、高い性能の計算機を用意する必要ということや、結果が得られるまでの時間がかかるなどの問題があった。
【0005】
また、心筋梗塞や致死的不整脈などの重症な心臓疾患は、異常の発生から死に至る時間が短く、迅速な異常の検知と評価が求められる。しかしながら、従来技術では心起電力を迅速に算出することは困難であったため、体表面の心電位が明確な異常を呈するまで判定することができず、病状が深刻になるという問題があった。
【0006】
上記事情に鑑み、本発明は、心起電力を推定する際の計算負荷を軽減する技術を提供することを目的としている。
【課題を解決するための手段】
【0007】
本発明の一態様は、心起電力の時系列の推定対象の心電位時系列を、心筋の１回の拍動によって生じる心電位を表す時系列が互いに分離されるように区分する区分ステップと、前記区分ステップにおいて区分された各時系列を単位心電位時系列として、単位心電位時系列ごとに単位心電位時系列の各データを脱分極期に属するデータと不応期に属するデータと再分極期に属するデータと休止期に属するデータとに分類することで、前記単位心電位時系列ごとに、脱分極期に属するデータの集合である脱分極心電位時系列と不応期に属するデータの集合である不応心電位時系列と再分極期に属するデータの集合である再分極心電位時系列とを取得する状態時系列取得ステップと、前記脱分極心電位時系列に基づき前記脱分極期の各時点における心起電力を示す脱分極心起電力時系列を推定する脱分極心起電力推定ステップと、前記不応心電位時系列と前記脱分極期の終了時点における心起電力の推定値とに基づき前記不応期の各時点における心起電力を示す不応心起電力時系列を推定する不応心起電力推定ステップと、前記再分極心電位時系列と前記不応期の終了時点における心起電力の推定値とに基づき前記再分極期の各時点における心起電力を示す再分極心起電力時系列を推定する再分極心起電力推定ステップと、を有する、心起電力推定方法である。
【発明の効果】
【0008】
本発明により、心起電力を推定する際の計算負荷を軽減することが可能となる。
【図面の簡単な説明】
【0009】
【図1】実施形態の心起電力推定方法の概略を説明する説明図。
【図2】実施形態の心起電力推定方法において実行される処理の流れの一例を示すフローチャート。
【図3】実施形態の主単位心起電力時系列推定処理の流れの一例を示すフローチャート。
【図4】実施形態の心起電力推定方法による推定結果の確からしさを、実測された心電位時系列を用いて説明する第１の説明図。
【図5】実施形態の心起電力推定方法による推定結果の確からしさを、実測された心電位時系列を用いて説明する第２の説明図。
【図6】実施形態の心起電力推定方法による推定結果の確からしさを、実測された心電位時系列を用いて説明する第３の説明図。
【図7】実施形態の心起電力推定方法による推定結果の確からしさを、実測された心電位時系列を用いて説明する第４の説明図。
【図8】実施形態の心起電力推定方法による推定結果の確からしさを、実測された心電位時系列を用いて説明する第５の説明図。
【図9】実施形態の心起電力推定方法による推定結果の確からしさを、実測された心電位時系列を用いて説明する第６の説明図。
【図10】実施形態の心起電力推定方法による推定結果の確からしさを、実測された心電位時系列を用いて説明する第７の説明図。
【図11】実施形態の心起電力推定方法による推定結果の確からしさを、実測された心電位時系列を用いて説明する第８の説明図。
【図12】実施形態の心起電力推定方法による推定結果の確からしさを、実測された心電位時系列を用いて説明する第９の説明図。
【図13】実施形態の心起電力推定方法による推定結果の確からしさを、実測された心電位時系列を用いて説明する第１０の説明図。
【図14】実施形態の心起電力推定方法による推定結果の確からしさを、実測された心電位時系列を用いて説明する第１１の説明図。
【図15】実施形態の心起電力推定方法による推定結果の確からしさを、実測された心電位時系列を用いて説明する第１２の説明図。
【図16】実施形態の心起電力推定方法による推定結果の確からしさを、実測された心電位時系列を用いて説明する第１３の説明図。
【図17】実施形態の心起電力推定方法による推定結果の確からしさを、実測された心電位時系列を用いて説明する第１４の説明図。
【図18】実施形態の心起電力推定方法による推定結果の確からしさを、実測された心電位時系列を用いて説明する第１５の説明図。
【図19】実施形態の心起電力推定方法による推定結果の確からしさを、実測された心電位時系列を用いて説明する第１６の説明図。
【図20】実施形態の心起電力推定方法による推定結果の確からしさを、実測された心電位時系列を用いて説明する第１７の説明図。
【図21】実施形態の心起電力推定方法による推定結果の確からしさを、実測された心電位時系列を用いて説明する第１８の説明図。
【図22】実施形態の心起電力推定方法による推定結果の確からしさを、実測された心電位時系列を用いて説明する第１９の説明図。
【図23】実施形態の心起電力推定方法による推定結果の確からしさを、実測された心電位時系列を用いて説明する第２０の説明図。
【図24】実施形態の心起電力推定方法による推定結果の確からしさを、実測された心電位時系列を用いて説明する第２１の説明図。
【図25】実施形態の心起電力推定方法による推定結果の確からしさを、実測された心電位時系列を用いて説明する第２２の説明図。
【図26】実施形態の心起電力推定方法による推定結果の確からしさを、実測された心電位時系列を用いて説明する第２３の説明図。
【図27】実施形態の心起電力推定方法による推定結果の確からしさを、実測された心電位時系列を用いて説明する第２４の説明図。
【図28】実施形態の心起電力推定方法による推定結果の確からしさを、実測された心電位時系列を用いて説明する第２５の説明図。
【図29】実施形態の心起電力推定方法による推定結果の確からしさを、実測された心電位時系列を用いて説明する第２６の説明図。
【図30】実施形態の心起電力推定方法による推定結果の確からしさを、実測された心電位時系列を用いて説明する第２７の説明図。
【図31】実施形態の心起電力推定方法による推定結果の確からしさを、実測された心電位時系列を用いて説明する第２８の説明図。
【図32】実施形態の心起電力推定方法による推定結果の確からしさを、実測された心電位時系列を用いて説明する第２９の説明図。
【図33】実施形態の心起電力推定方法による推定結果の確からしさを、実測された心電位時系列を用いて説明する第３０の説明図。
【図34】実施形態の心起電力推定方法による推定結果の確からしさを、実測された心電位時系列を用いて説明する第３１の説明図。
【図35】実施形態の心起電力推定方法による推定結果の確からしさを、実測された心電位時系列を用いて説明する第３２の説明図。
【図36】実施形態の心起電力推定方法による推定結果の確からしさを、実測された心電位時系列を用いて説明する第３３の説明図。
【図37】実施形態の心起電力推定方法を実行する心起電力推定システム１００の構成の一例を示す図。
【図38】変形例における一致処理が奏する効果を説明する説明図。
【発明を実施するための形態】
【0010】
以下、実施形態の心起電力推定方法を、図面を参照して説明する。なお、本発明における「脱分極期」は、心電図におけるＱＲＳと定義した。脱分極期は、心室の心筋細胞の集団における脱分極相及びスパイク（Ｐｈａｓｅ０及びＰｈａｓｅ１）に相当する。
【0011】
また、本発明における「再分極期」は、心電図のＴと定義した。再分極期は、心室の心筋細胞の集団における再分極相（Ｐｈａｓｅ３）に相当する。また、本発明における「不応期」は、心電図のＳＴセグメントと定義した。不応期は、心室の心筋細胞の集団におけるプラトー相（Ｐｈａｓｅ２）に相当する。本発明における「休止期」は、脱分極期、不応期又は再分極期のいずれでもない期間と定義した。休止期は、心室の心筋細胞の集団における静止電位の期間（Ｐｈａｓｅ４）に相当する。
【0012】
なお、心室の心筋細胞の集団における脱分極期のスパイク（Ｐｈａｓｅ１）については、脱分極期から省略し、代わりに、不応期に含めても良い。
【0013】
図１は、実施形態の心起電力推定方法の概略を説明する説明図である。心起電力推定方法は、心電位の時系列（以下「心電位時系列」という。）に基づき、推定対象の心電位時系列が示す期間における心起電力の時系列（以下「心起電力時系列」という。）を推定する方法である。推定対象は人であってもよいし動物であってもよい。以下、説明の簡単のため推定対象が人である場合を例に心起電力推定方法を説明する。図１において、時系列Ａ１００は、心電位時系列の一例である。図１において心電位の単位はマイクロボルトである。図１において時刻の単位はミリ秒である。
【0014】
心起電力推定方法では、まず心電位時系列を取得する。心起電力推定法では、取得した心電位時系列を、心筋の１回の拍動によって生じる心電位を表す時系列が互いに分離されるように区分する。心筋の１回の拍動とは、心筋に脱分極が生じた時点から、脱分極期と、不応期と、再分極期とを経て休止期に入るまでの期間における心筋の一連の動作である。
【0015】
図１において時系列Ａ２００は、心筋の１回の拍動によって生じる心電位を表す時系列を含む時系列であって区分後の１つの時系列（以下「単位心電位時系列」という。）の一例である。単位心電位時系列は、例えば、始まりがＰ波の開始の時点であり終わりが次のＰ波の開始直前の時点である期間における心電位時系列である。単位心電位時系列は、例えば、Ｐ波の検出が困難である場合、基線からＲ波Ｔ波を経て基線に戻るまでという条件を満たす期間における心電位時系列であってもよい。
【0016】
単位心電位時系列は、脱分極心電位時系列、不応心電位時系列、再分極心電位時系列及び休止心電位時系列を有する。脱分極心電位時系列は、脱分極期に生じた心電位を示す時系列である。不応心電位時系列は、不応期に生じた心電位を示す時系列である。再分極心電位時系列は、再分極期に生じた心電位を示す時系列である。休止心電位時系列は、休止期に生じた心電位を示す時系列である。すなわち休止心電位時系列は、単位心電位時系列が有する各データのうち、脱分極心電位時系列、不応心電位時系列又は再分極心電位時系列のいずれにも属さないデータの集合である。
【0017】
図１において、時刻Ｔ１から時刻Ｔ２までの期間における時系列は脱分極心電位時系列の一例である。図１において、時刻Ｔ２から時刻Ｔ３までの期間における時系列は不応心電位時系列の一例である。図１において、時刻Ｔ３から時刻Ｔ４までの期間における時系列は再分極心電位時系列の一例である。図１において、時刻Ｔ０から時刻Ｔ１までの期間における時系列と時刻Ｔ４以降の時系列とは休止心電位時系列の一例である。
【0018】
【0019】
心起電力推定法では、単位心電位時系列ごとに各単位心電位時系列が示す期間における心起電力の時系列（以下「単位心起電力時系列」という。）を推定する。具体的には、まず、単位心電位時系列が示す各データを、脱分極心電位時系列に属するデータと、不応心電位時系列に属するデータと、再分極心電位時系列に属するデータと、休止心電位時系列に属するデータとに分類する。
【0020】
以下、脱分極心電位時系列に属するデータを脱分極心電位時系列データという。以下、不応心電位時系列に属するデータを不応心電位時系列データという。以下、再分極心電位時系列に属するデータを再分極心電位時系列データという。以下、休止心電位時系列に属するデータを休止心電位時系列データという。
【0021】
次に、脱分極期、不応期及び再分極期の各期間に応じた規則にしたがって心起電力時系列を推定する。以下、図２及び図３を用いて、単位心電位時系列ごとに心起電力の時系列を推定する具体的な推定方法の説明とともに心起電力推定方法を説明する。
【0022】
図２は、実施形態の心起電力推定方法において実行される処理の流れの一例を示すフローチャートである。心起電力推定方法は、例えばコンピュータが実行する。
【0023】
心電位時系列を取得する（ステップＳ１００）。次に、取得した心電位時系列を、単位心電位時系列ごとに区分する（ステップＳ２００）。ステップＳ２００の処理は、心電位時系列を単位心電位時系列で区分する処理であるので単位心電位時系列を生成する処理である。
【0024】
次に、主心起電力時系列推定処理を実行する（ステップＳ３００）。主心起電力時系列推定処理は、ステップＳ２００で生成された全ての単位心電位時系列について各単位心電位時系列が示す期間における単位心起電力時系列を推定する処理である。
【0025】
主心起電力時系列推定処理は、具体的には、副心起電力時系列推定処理を全ての単位心電位時系列に対して実行する処理である。副心起電力時系列推定処理は、１つの単位心電位時系列に対して実行される処理であって実行対象の単位心電位時系列が示す期間における単位心起電力時系列を推定する処理である。以下、副心起電力時系列推定処理の実行対象の単位心電位時系列を対象単位心電位時系列という。
【0026】
図３は、実施形態の主単位心起電力時系列推定処理の流れの一例を示すフローチャートである。所定の規則に従い、ステップＳ２００で生成された単位心電位時系列の１つを対象単位心電位時系列に決定する（ステップＳ３０１）。所定の規則は、例えば、ステップＳ２００で取得された単位心電位時系列のうち最も早い期間の単位心電位時系列を対象単位心電位時系列に決定する、という規則である。
【0027】
以下で説明するステップＳ３０２からステップＳ３０５の処理が副心起電力時系列推定処理の一例である。ステップＳ３０１の次に、対象単位心電位時系列の各データを、脱分極心電位時系列データ、不応心電位時系列データ、再分極心電位時系列データ又は休止心電位時系列データのいずれか１つに分類する（ステップＳ３０２）。
【0028】
ステップＳ３０２の処理によって脱分極心電位時系列データに分類されたデータの集合が脱分極心電位時系列である。ステップＳ３０２の処理によって不応心電位時系列データに分類されたデータの集合が不応心電位時系列である。ステップＳ３０２の処理によって再分極心電位時系列データに分類されたデータの集合が再分極心電位時系列である。ステップＳ３０２の処理によって休止心電位時系列データに分類されたデータの集合が休止心電位時系列である。
【0029】
このように、ステップＳ３０２の処理は、脱分極心電位時系列、不応心電位時系列、再分極心電位時系列及び休止心電位時系列を取得する処理である。
【0030】
対象単位心電位時系列分類方法は、対象単位心電位時系列の各データを、脱分極心電位時系列データ、不応心電位時系列データ、再分極心電位時系列データ又は休止心電位時系列データのいずれか１つに分類できればどのような方法であってもよい。対象単位心電位時系列分類方法は、対象単位心電位時系列の各データを、脱分極心電位時系列データ、不応心電位時系列データ、再分極心電位時系列データ又は休止心電位時系列データのいずれか１つに分類する方法である。
【0031】
対象単位心電位時系列分類方法において脱分極心電位時系列データを判定する方法は例えば、脱分極期の開始時点と脱分極期の終了時点とを判定した後、データの示す時点が脱分極期か否かを判定する方法である。期間の開始時点は期間の開始の時点を意味する。期間の終了時点は期間の終了の時点を意味する。
【0032】
脱分極期の開始時点は例えば、対象単位心電位時系列の示す心電位の単位時間当たりの変化量が所定の量以上という条件を最初に満たした時点である。脱分極期の開始時点は、例えば、２００Ｈｚサンプリング周波数の場合、単位時間５ｍｓｅｃ当たりの心電位の変化量の絶対値が４０μボルト以上という条件であってもよい。閾値はこの値に限らず心電位の大きさや皮膚の性状、生体電極、計測機器の利得などの計測条件によって調整された値であってもよい。また、脱分極期の開始時点は、例えば、心電位の微分等によって求めた変曲点から算出された時点であってもよい。
【0033】
脱分極期の終了時点は例えば、脱分極の開始時点以降の時点であって、対象単位心電位時系列の示す心電位の時間軸正方向の傾きの符号が少なくとも１回反転した後に対象単位心電位時系列の示す心電位の変化量が初めて所定の閾値を下回った時点である。脱分極期の終了時点は、例えば、単位時間５ｍｓｅｃ当たりの心電位の変化量の絶対値が３０μボルト以下という条件であってもよい。閾値はこの値に限らず心電位の大きさや皮膚の性状、雑音、生体電極、計測機器の利得などの計測条件によって調整された値であってもよい。また、脱分極期の終了時点は、例えば、心電位の微分等によって求めた変曲点から算出された時点であってもよい。
【0034】
対象単位心電位時系列分類方法において不応心電位時系列データを判定する方法は例えば、不応期の開始時点と不応期の終了時点とを判定した後、データの示す時点が不応期か否かを判定する方法である。
【0035】
不応期の開始時点は、脱分極期の終了時点の直後の時点である。不応期の終了時点は例えば、対象単位心電位時系列の示す心電位の単位時間当たりの負方向の変化量が所定の量以上という条件を不応期の開始時点以降に最初に満たした時点である。不応期の終了時点は、例えば、心電位の変化量の絶対値が１０μボルト以上という条件であってもよい。閾値はこの値に限らず心電位の大きさや皮膚の性状、雑音、生体電極、計測機器の利得などの計測条件によって調整された値であってもよい。また、不応期の終了時点は、例えば、心電位の微分等によって求めた変曲点から算出された時点であってもよい。
【0036】
対象単位心電位時系列分類方法において再分極心電位時系列データを判定する方法は例えば、再分極期の開始時点と再分極期の終了時点とを判定した後、データの示す時点が再分極期か否かを判定する方法である。
【0037】
再分極期の開始時点は、不応期の終了時点の直後の時点である。再分極期の終了時点は例えば、対象単位心電位時系列の示す心電位の単位時間当たりの変化量の絶対値が所定の量未満という条件を再分極期の開始時点以降に最初に満たした時点である。再分極期の終了時点は、心電位の変化量の絶対値が７μボルト未満であってもよい。閾値はこの値に限らず心電位の大きさや皮膚の性状、雑音、生体電極、計測機器の利得などの計測条件によって調整された値であってもよい。また、再分極期の終了時点は、例えば、心電位の微分等によって求めた変曲点から算出された時点であってもよい。
【0038】
ステップＳ３０２の次に、ステップＳ３０２で取得された脱分極心電位時系列に基づき、脱分極期における心起電力の時系列（以下「脱分極心起電力時系列」という。）を推定する（ステップＳ３０３）。脱分極心起電力時系列は、脱分極心電位時系列の時点ごとに脱分極心起電力推定値を有するデータである。各時点の脱分極心起電力推定値は、脱分極期の各時点における心起電力の推定値である。
【0039】
脱分極心起電力推定値は具体的には、ステップＳ３０２で取得された脱分極心電位時系列が示す脱分極期の開始時点の心電位に、ステップＳ３０２で取得された脱分極心電位時系列が示す過去の脱分極差分電位の総和を足し算した値である。
【0040】
脱分極差分電位は、脱分極期における隣接する２つの時点間の心電位の差の絶対値である。なお、過去の脱分極差分電位は、異なる単位心電位時系列における脱分極差分電位は含まない。そのため、過去の脱分極差分電位とは、脱分極期の開始時点以降の過去の脱分極差分電位である。
【0041】
このように脱分極心起電力時系列を推定する処理は、脱分極期の時点ごとに各時点における心起電力を推定する処理である。具体的には、脱分極心起電力時系列を推定する処理は、脱分極期の開始時点から推定対象の心起電力に対応する時点までの隣接する心電位の差の絶対値の総和を脱分極期の開始時点における心電位に足し算した値を心起電力として取得する処理である。
【0042】
ステップＳ３０３の次に、ステップＳ３０２で取得された不応心電位時系列とステップＳ３０３で取得された終了時点脱分極心起電力推定値とに基づき不応期における心起電力の時系列（以下「不応心起電力時系列」という。）を推定する（ステップＳ３０４）。終了時点脱分極心起電力推定値は、脱分極期の終了時点における脱分極心起電力推定値である。
【0043】
不応心起電力時系列は、不応心電位時系列の時点ごとに不応心起電力推定値を有するデータである。各時点の不応心起電力推定値は、不応期の各時点における心起電力の推定値である。各時点の不応心起電力推定値は具体的には、終了時点脱分極心起電力推定値の値から不応心電位時系列の対応する時点における心電位に関する値を引き算した値である。
【0044】
引き算に用いる不応心電位時系列の心電位に関する値は、例えば、終了時点脱分極心起電力推定値の値から不応心電位時系列の対応する時点における心電位を引き算した値である。引き算に用いる不応心電位時系列の心電位に関する値は、心電位の実測値または、心電位の絶対値であってもよい。心電位のパターンによっては、基線との位置関係により、反転や折り返しや屈曲などの不自然な変形を生じる場合がある。そのような場合、引き算に用いる不応心電位時系列の心電位に関する値は、不応期心電位に変形処理を加えてから引き算した値であっても良い。この変形では、ゲイン（倍数）やバイアス（定数）を加える、累積を用いるなどの調整が行われてもよい。なお、不応心起電力推定値の定義における、対応する時点とは、不応心電位時系列データが示す時点であって、不応心起電力推定値を示す不応心起電力時系列データが示す時点に同一の時点である。
【0045】
このように不応心起電力時系列を推定する処理は、不応期の時点ごとに各時点における心起電力を推定する処理である。具体的には、脱分極期の終了時点における心起電力の推定値から不応期の対応する時点における心電位に関する値を引き算した値を心起電力として取得する処理である。
【0046】
不応心起電力時系列を推定する処理は、例えば、不応期の時点ごとに各時点における心起電力を推定する処理である。具体的には、脱分極期の終了時点における心起電力の推定値から不応期の対応する時点における心電位を引き算した値を心起電力として取得する処理である。不応心起電力時系列を推定する処理は、例えば、脱分極期の終了時点における心起電力の推定値から不応期の対応する時点における変形した心電位を引き算した値を心起電力として取得する処理であってもよい。
【0047】
以下、説明の簡単のため、引き算に用いる不応心電位時系列の心電位に関する値が終了時点脱分極心起電力推定値の値から不応心電位時系列の対応する時点における心電位を引き算した値である場合を例に、心起電力推定方法を説明する。このような場合、不応心起電力時系列を推定する処理は、不応心起電力時系列を推定する処理は、例えば、不応期の時点ごとに各時点における心起電力を推定する処理である。
【0048】
ステップＳ３０４の次に、ステップＳ３０２で取得された再分極心電位時系列とステップＳ３０４で取得された終了時点不応心起電力推定値に基づき再分極期における心起電力の時系列（以下「再分極心起電力時系列」という。）を推定する（ステップＳ３０５）。終了時点不応心起電力推定値は、不応期の終了時点における不応心起電力推定値である。終了時点不応心起電力推定値は再分極期の開始時点における再分極心起電力推定値に一致する。
【0049】
再分極心起電力時系列は、再分極心電位時系列の時点ごとに再分極心起電力推定値を有するデータである。各時点における再分極心起電力推定値は、再分極期の各時点における心起電力の推定値である。各時点における再分極心起電力推定値は具体的には、終了時点不応心起電力推定値から、ステップＳ３０２で取得された再分極心電位時系列が示す過去の再分極差分電位の総和を引き算した値である。
【0050】
再分極差分電位は、再分極期における隣接する２つの時点間の心電位の差に基づく値である。再分極差分電位は、例えば、再分極期における隣接する２つの時点間の心電位の差の絶対値である。なお、過去の再分極差分電位は、異なる単位心電位時系列における再分極差分電位は含まない。そのため、過去の再分極差分電位とは、不応期の終了時点以降の過去の再分極差分電位である。この引き算に用いる再分極期心電位は、例えば基線と交差する場合などに、上述の心電位の絶対値ではなく、上述の心電位の絶対値にバイアスやゲインを加える、累積を用いるなどの調整が行われた値であってもよい。以下、説明の簡単のため、再分極差分電位が再分極期における隣接する２つの時点間の心電位の差の絶対値である場合を例に、心起電力推定方法を説明する。
【0051】
このように再分極心起電力時系列を推定する処理は、再分極期の時点ごとに各時点における心起電力を推定する処理である。具体的には、再分極心電位時系列における不応期の終了時点から、推定される心起電力に対応する時点まで、の隣接する心電位の差の絶対値の総和を不応期の終了時点における心起電力の推定値から引き算した値を心起電力として取得する処理である。
【0052】
ステップＳ３０２からステップＳ３０５の一連の処理の１回の実行によって推定された脱分極心起電力時系列、不応心起電力時系列及び再分極心起電力時系列の集合が、対象単位心電位時系列が示す期間における推定された単位心起電力時系列である。
【0053】
ステップＳ３０５の次に、ステップＳ２００で生成された全ての単位心電位時系列について単位心起電力時系列を推定したか否かを判定する（ステップＳ３０６）。すなわち、未推定単位心電位時系列が無いか否かを判定する。未推定単位心電位時系列は、副心起電力時系列推定処理が実行されていない単位心電位時系列である。
【0054】
未推定単位心電位時系列がある場合（ステップＳ３０６：ＹＥＳ）、所定の規則に則り、未推定単位心電位時系列のうちの１つを対象単位心電位時系列に決定する（ステップＳ３０７）。所定の規則は、例えば、未推定単位心電位時系列のうちの最も早い期間の未推定単位心電位時系列を対象単位心電位時系列に決定する、という規則である。ステップＳ３０７の次に、ステップＳ３０２に戻る。一方、未推定単位心電位時系列が無い場合（ステップＳ３０６：ＮＯ）、処理を終了する。
【0055】
図１の説明に戻る。ステップＳ３００の処理によって推定された各単位心起電力時系列の集合を、ステップＳ１００で取得された心電位時系列が示す期間における心起電力時系列の推定結果として取得する（ステップＳ４００）。
【0056】
以下、図４～図３４を用いて実施形態の心起電力推定方法による推定結果の確からしさの検証結果を、実測された心電位時系列を用いて説明する。図４～図３４の説明においては説明の簡単のため、心電位時系列が心筋の１回の拍動によって生じる心電位の時系列である場合を例に心起電力推定方法を説明する。そのため、推定結果の心起電力系列は、１つの単位心電位時系列である。
【0057】
【0058】
【0059】
図４は、実施形態の心起電力推定方法による推定結果の確からしさを、実測された心電位時系列を用いて説明する第１の説明図である。より具体的には、図４は心筋の動作が正常な人であって運動負荷試験が行われる前の人の心電位時系列の一例とその心電位時系列を用いて心起電力推定方法によって推定した推定結果とを示す。
【0060】
図４の時系列Ｄ１０１は、実測された心電位時系列を示す。図４の時系列Ｄ１０２は、時系列Ｄ１０１に基づき心起電力推定方法によって推定された心起電力時系列を示す。推定結果の時系列Ｄ１０２は、心筋の動作が正常な人の運動負荷試験実行前の心起電力時系列が有する特徴として医学的に認知されている特徴と同様の特徴を備える。具体的には、推定結果の時系列Ｄ１０２は、脱分極電位の持続によるプラトーを有する心筋の活動電位の特徴を備え、単相性心筋活動電位記録等で観測される心筋の波形の特徴を具備している（非特許文献１及び４参照）。そのため、図４は、実施形態の心起電力推定方法について、心筋の動作が正常な人の運動負荷試験実行前の心起電力時系列を高い精度で推定する方法であることを示す。
【0061】
図５は、実施形態の心起電力推定方法による推定結果の確からしさを、実測された心電位時系列を用いて説明する第２の説明図である。より具体的には、図５は心筋に一過性の虚血を呈した人の心電位時系列の一例とその心電位時系列を用いて心起電力推定方法によって推定した推定結果とを示す。
【0062】
図５の時系列Ｄ１０３は、実測された心電位時系列を示す。図５の時系列Ｄ１０４は、時系列Ｄ１０３に基づき心起電力推定方法によって推定された心起電力時系列を示す。推定結果の時系列Ｄ１０４は、心筋の虚血時の心起電力時系列が有する特徴として医学的に認知されている特徴と同様の特徴を備える。具体的には、脱分極の遅延、プラトーの上昇や下降、屈曲などの変化、活動電位持続時間の延長の特徴を備える。そのため、図５は、心起電力推定方法について、心筋の虚血時の心起電力時系列を高い精度で推定する方法であることを示す。
【0063】
図６は、実施形態の心起電力推定方法による推定結果の確からしさを、実測された心電位時系列を用いて説明する第３の説明図である。より具体的には、図６は心筋の虚血時の心電位時系列の一例とその心電位時系列を用いて心起電力推定方法によって推定した推定結果とを示す。
【0064】
図６の時系列Ｄ１０５は、実測された心電位時系列を示す。図６の時系列Ｄ１０６は、時系列Ｄ１０５に基づき心起電力推定方法によって推定された心起電力時系列を示す。推定結果の時系列Ｄ１０６は、心筋の虚血時の心起電力時系列が有する特徴として医学的に認知されている特徴と同様の特徴を備える。具体的には、脱分極の遅延、プラトーの上昇や下降、屈曲などの変化、活動電位持続時間の延長特徴を備える。そのため、図６は、心起電力推定方法について、心筋の虚血時の心起電力時系列を高い精度で推定する方法であることを示す。
【0065】
図７は、実施形態の心起電力推定方法による推定結果の確からしさを、実測された心電位時系列を用いて説明する第４の説明図である。より具体的には、図７は心筋の虚血時の心電位時系列の一例とその心電位時系列を用いて心起電力推定方法によって推定した推定結果とを示す。
【0066】
図７の時系列Ｄ１０７は、実測された心電位時系列を示す。図７の時系列Ｄ１０８は、時系列Ｄ１０７に基づき心起電力推定方法によって推定された心起電力時系列を示す。推定結果の時系列Ｄ１０８は、心筋の虚血時の心起電力時系列が有する特徴として医学的に認知されている特徴と同様の特徴を備える。具体的には、脱分極の遅延、プラトーの上昇や下降、屈曲などの変化、活動電位持続時間の延長の特徴を備える。そのため、図７は、心起電力推定方法について、心筋の虚血時の心起電力時系列を高い精度で推定する方法であることを示す。
【0067】
図８は、実施形態の心起電力推定方法による推定結果の確からしさを、実測された心電位時系列を用いて説明する第５の説明図である。より具体的には、図８は心筋の虚血時の心電位時系列の一例とその心電位時系列を用いて心起電力推定方法によって推定した推定結果とを示す。
【0068】
図８の時系列Ｄ１０９は、実測された心電位時系列を示す。図８の時系列Ｄ１１０は、時系列Ｄ１０９に基づき心起電力推定方法によって推定された心起電力時系列を示す。推定結果の時系列Ｄ１１０は、心筋の虚血時の心起電力時系列が有する特徴として医学的に認知されている特徴と同様の特徴を備える。具体的には、脱分極の遅延、プラトーの上昇や下降、屈曲などの変化、活動電位持続時間の延長の特徴を備える。そのため、図８は、心起電力推定方法について、心筋の虚血時の心起電力時系列を高い精度で推定する方法であることを示す。
【0069】
図９は、実施形態の心起電力推定方法による推定結果の確からしさを、実測された心電位時系列を用いて説明する第６の説明図である。より具体的には、図９は、図５、図６、図７及び図８に示す時系列を重ねて表示した結果を示す。
【0070】
図９は、縦軸が電位を示し横軸が時刻を示す。縦軸が示す電位は、時系列Ｄ１０３、Ｄ１０５、Ｄ１０７及びＤ１０９については心電位を示し、時系列Ｄ１０４、Ｄ１０６、Ｄ１０８及びＤ１１０については心起電力を示す。
【0071】
図９は、心起電力推定方法によって推定された２つの心起電力時系列が示す各形状には、心電位時系列の違いを反映した違いがあることを示す。具体的には、図９は、心起電力時系列間には不応期における電位差があり、心電位時系列にも不応期における電位差があることを示す。また、図９の不応期における違いは、心起電力時系列間の違いの方が心電位時系列間における違いよりも大きい。そのため、図９は、心電位時系列間の違いを増幅した違いが心起電力時系列間の違いとして現れることを示す。
【0072】
図１０～図３４を用いて、実施形態の心起電力推定方法による推定結果の確からしさを説明する。図１０～図３４は、それぞれ実施形態の心起電力推定方法による推定結果の確からしさを、実測された心電位時系列を用いて説明する第７～第３１の説明図である。より具体的には、図１０～図３４の各結果はそれぞれ、心筋の動作が正常な人の心電位時系列とその心電位時系列を用いて心起電力推定方法によって推定した推定結果とを示す例である。
【0073】
図１０～図３４は、結果Ｅ１～結果Ｅ２５を示す。結果Ｅ１～結果Ｅ２５はそれぞれ、実測された心電位時系列と心起電力推定方法による推定結果の心起電力時系列とを示す。結果Ｅ１～結果Ｅ２５の各グラフの縦軸は電位を示す。電位の単位はマイクロボルトである。結果Ｅ１～結果Ｅ２５の各グラフの横軸は時刻を示す。時刻の単位はミリ秒である。
【0074】
図１０～図３４において、時系列Ｄ２０１、Ｄ２０３、Ｄ２０５、Ｄ２０７、Ｄ２０９、Ｄ２１１、Ｄ２１３、Ｄ２１５、Ｄ２１７、Ｄ２１９、Ｄ２２１、Ｄ２２３、Ｄ２２５、Ｄ２２７、Ｄ２２９、Ｄ２３１、Ｄ２３３、Ｄ２３５、Ｄ２３７、Ｄ２３９、Ｄ２４１、Ｄ２４３、Ｄ２４５、Ｄ２４７及びＤ２４９はそれぞれ実測された心電位時系列である。
【0075】
図１０～図３４において、時系列Ｄ２０２、Ｄ２０４、Ｄ２０６、Ｄ２０８、Ｄ２１０、Ｄ２１２、Ｄ２１４、Ｄ２１６、Ｄ２１８、Ｄ２２０、Ｄ２２２、Ｄ２２４、Ｄ２２６、Ｄ２２８、Ｄ２３０、Ｄ２３２、Ｄ２３４、Ｄ２３６、Ｄ２３８、Ｄ２４０、Ｄ２４２、Ｄ２４４、Ｄ２４６、Ｄ２４８及びＤ２５０はそれぞれ心起電力推定方法による推測結果の心起電力時系列である。
【0076】
時系列Ｄ２０２、Ｄ２０４、Ｄ２０６、Ｄ２０８、Ｄ２１０、Ｄ２１２、Ｄ２１４、Ｄ２１６、Ｄ２１８、Ｄ２２０、Ｄ２２２、Ｄ２２４、Ｄ２２６、Ｄ２２８、Ｄ２３０、Ｄ２３２、Ｄ２３４、Ｄ２３６、Ｄ２３８、Ｄ２４０、Ｄ２４２、Ｄ２４４、Ｄ２４６、Ｄ２４８及びＤ２５０はいずれも心筋の動作が正常な人の心起電力時系列が有する特徴として医学的に認知されている特徴と同様の特徴を備える。そのため、図１０～図３４の結果Ｅ１～結果Ｅ２５は、実施形態の心起電力推定方法について、心筋の動作が正常な人の心起電力時系列を高い精度で推定する方法であることを示す。
【0077】
図３５は、実施形態の心起電力推定方法による推定結果の確からしさを、実測された心電位時系列を用いて説明する第３３の説明図である。図３５は時系列Ｄ３０１及び時系列Ｄ３０２を示す。時系列Ｄ３０１は実測された心電位時系列である。時系列Ｄ３０２は時系列Ｄ３０１に基づき心起電力推定方法によって推定された心起電力時系列である。
【0078】
時系列Ｄ３０１は、時刻Ｔ５から時刻Ｔ６までの期間に心筋の動作に期外収縮が起きていることを示す。心筋の動作に期外収縮が起きているとは、心筋の動作が異常であることを意味する。時系列Ｄ３０１は、それ以外の期間では心筋の動作が洞調律であることを示す。心筋の動作が洞調律であるとは、心筋の動作が正常であることを意味する。
【0079】
時系列Ｄ３０２は、期外収縮時の心起電力の波形に医学的に知られた特徴が生じていることを示す。時系列Ｄ３０２において期外収縮時の波形は、時刻Ｔ５から時刻Ｔ６までの期間の波形である。期外収縮時の心起電力の波形の特徴は、具体的には、プラトー相（すなわち不応期）の長さが正常の場合に比較して短いという特徴と、再分極相（すなわち再分極期）の開始が正常の場合に比較して早いという特徴とである。
【0080】
このように図３５は、実施形態の心起電力推定方法について、期外収縮を生じている人の心起電力時系列を高い精度で推定する方法であることを示す。
【0081】
図３６は、実施形態の心起電力推定方法による推定結果の確からしさを、実測された心電位時系列を用いて説明する第３４の説明図である。図３６は時系列Ｄ４０１及び時系列Ｄ４０２を示す。時系列Ｄ４０１は実測された心電位時系列である。時系列Ｄ４０２は時系列Ｄ４０１に基づき心起電力推定方法によって推定された心起電力時系列である。
【0082】
時系列Ｄ４０１は、時刻Ｔ７から時刻Ｔ８までの期間に心筋の動作に期外収縮が起きていることを示す。時系列Ｄ４０１は、時刻Ｔ９から時刻Ｔ１０までの期間に心筋の動作に期外収縮が起きていることを示す。
【0083】
時系列Ｄ４０２は、時刻Ｔ７から時刻Ｔ８までの期間の波形が、プラトー相の長さが正常の場合に比較して短いという特徴と、再分極相の開始が正常の場合に比較して早いという特徴とを有することを示す。また、時系列Ｄ４０２の波形は時刻Ｔ７から時刻Ｔ８までの期間の波形が、房室ブロックに伴う異常であることも示す。房室ブロックに伴う異常は、異所性の刺激伝導路によるＱＲＳ波やＳＴの異形性という特徴として波形に現れている。
【0084】
時系列Ｄ４０２は、時刻Ｔ９から時刻Ｔ１０までの期間の波形が、プラトー相の長さが正常の場合に比較して短いという特徴と、再分極相の開始が正常の場合に比較して早いという特徴とを有することを示す。また、時系列Ｄ４０２の波形は時刻Ｔ７から時刻Ｔ８までの期間の波形が、房室ブロックに伴う異常であることも示す。房室ブロックに伴う異常は、異所性の刺激伝導路によるＱＲＳ波やＳＴの異形性という特徴として波形に現れている。
【0085】
このように、時系列Ｄ４０２は、時刻Ｔ７から時刻Ｔ８までの期間の波形に、房室ブロックに伴う期外収縮時の心起電力の波形の医学的に知られた特徴が生じていることを示す。また、時系列Ｄ４０２は、時刻Ｔ９から時刻Ｔ１０までの期間の波形に、房室ブロックに伴う期外収縮時の心起電力の波形の医学的に知られた特徴が生じていることを示す。
【0086】
このように図３６は、実施形態の心起電力推定方法について、期外収縮を生じている人の心起電力時系列を高い精度で推定する方法であることを示す。
【0087】
図３７は、実施形態の心起電力推定方法を実行する心起電力推定システム１００の構成の一例を示す図である。心起電力推定システム１００は、心起電力推定装置１及び心電位測定装置２を備える。心起電力推定装置１は、心電位測定装置２が測定した心電位時系列を取得し、取得した心電位時系列に基づき取得した心電位時系列が示す期間における心起電力時系列を推定する。
【0088】
心電位測定装置２は、心起電力時系列の推定対象の心電位を測定する。心電位測定装置２は、測定結果の心電位の時系列である心電位時系列を心起電力推定装置１に出力する。心電位測定装置２は、例えば、心電計である。
【0089】
心起電力推定装置１は、バスで接続されたＣＰＵ（Central Processing Unit）等のプロセッサ９１とメモリ９２とを備える制御部１０を備え、プログラムを実行する。心起電力推定装置１は、プログラムの実行によって制御部１０、通信部１１、記憶部１２及びユーザインタフェース１３を備える装置として機能する。
【0090】
より具体的には、心起電力推定装置１は、プロセッサ９１が記憶部１２に記憶されているプログラムを読み出し、読み出したプログラムをメモリ９２に記憶させる。プロセッサ９１が、メモリ９２に記憶させたプログラムを実行することによって、心起電力推定装置１は、制御部１０、通信部１１、記憶部１２及びユーザインタフェース１３を備える装置として機能する。
【0091】
制御部１０は、心起電力推定装置１が備える各機能部の動作を制御する。制御部１０は例えば通信部１１の動作を制御する。制御部１０は、例えば通信部１１の動作を制御して、心電位測定装置２から心電位時系列を取得する。制御部１０は取得した心電位時系列を記憶部１２に記録する。
【0092】
制御部１０は例えば、心起電力推定方法を用いて、心起電力時系列を推定する。より具体的には、制御部１０は、通信部１１又はユーザインタフェース１３を介して入力された心電位時系列を取得し、取得した心電位時系列に基づき心起電力推定方法によって心電位時系列が示す期間における心起電力時系列を推定する。制御部１０は取得した心起電力時系列を記憶部１２に記録する。
【0093】
制御部１０は例えば、出力部１３２の動作を制御して、出力部１３２に推定結果の心起電力時系列を出力させる。
【0094】
通信部１１は、心電位測定装置２を含む１又は複数の外部装置に心起電力推定装置１を接続するための通信インタフェースを含んで構成される。外部装置の１つは心電位測定装置２である。通信部１１は、心電位測定装置２から心電位時系列を取得する。外部装置の１つは、例えば、推定結果の心起電力時系列を出力するプリンタである。このような場合、通信部１１は外部装置の１つであるプリンタに推定結果の心起電力時系列を送信する。
【0095】
記憶部１２は、磁気ハードディスク装置や半導体記憶装置などの記憶装置を用いて構成される。記憶部１２は心起電力推定装置１に関する各種情報を記憶する。記憶部１２は例えば、心起電力推定装置１が備える各機能部の動作を制御するプログラムを予め記憶する。記憶部１２は例えば、予め心起電力推定方法を実行するプログラムを記憶する。記憶部１２は例えば、推定結果の心電位時系列を記憶する。記憶部１２は例えば、推定結果の心起電力時系列を記憶する。
【0096】
ユーザインタフェース１３は、心起電力推定装置１に対する入力を受け付ける入力部１３１と心起電力推定装置１に関する各種情報を表示する出力部１３２とを備える。ユーザインタフェース１３は、例えば、タッチパネルである。入力部１３１は、自装置に対する入力を受け付ける。入力部１３１は、例えばマウスやキーボード、タッチパネル等の入力端末である。入力部１３１は、例えば、これらの入力端末を自装置に接続するインタフェースとして構成されてもよい。入力部１３１が受け付ける入力は、例えば、心起電力の推定の開始を指示する入力である。入力部１３１が受け付ける入力は、例えば、心電位時系列であってもよい。
【0097】
出力部１３２は、例えば液晶ディスプレイ、有機ＥＬ（Electro Luminescence）ディスプレイ等の表示装置である。出力部１３２は、例えば、これらの表示装置を自装置に接続するインタフェースとして構成されてもよい。出力部１３２は、例えばスピーカー等の音声出力装置であってもよい。出力部１３２が出力する情報は、例えば、推定結果の心起電力時系列である。出力部１３２の動作は制御部１０によって制御される。そのため、出力部１３２は制御部１０の指示にしたがって推定結果の心起電力時系列を出力する。
【0098】
このように構成された心起電力推定方法は、心電位時系列の各データを拍動ごとに脱分極心電位時系列データ、不応心電位時系列データ、再分極心電位時系列データ及び休止心電位時系列データの４種類のいずれか１つに分類する。
【0099】
心起電力推定方法は、分類の結果得られる脱分極心電位時系列に基づき、脱分極心起電力時系列を推定する。上述したように、脱分極心起電力時系列は、脱分極心電位時系列の時点ごとに脱分極心起電力推定値を有するデータである。上述したように、脱分極心起電力推定値は、脱分極心電位時系列が示す脱分極期の開始時点の心電位に、脱分極心電位時系列が示す過去の脱分極差分電位の総和を足し算した値である。
【0100】
そのため心起電力推定方法は、推定シミュレーション法よりも計算負荷の小さい演算によって脱分極心起電力時系列を推定することができる。推定シミュレーションは、心電位時系列を再現する数理モデルを用いたシミュレーションにより逆問題を解くことで心起電力を推定する方法である。
【0101】
心起電力推定方法は、分類の結果得られる不応心電位時系列に基づき、不応心起電力時系列を推定する。上述したように、不応心起電力時系列は、不応心電位時系列の時点ごとに不応心起電力推定値を有するデータである。上述したように、各時点の不応心起電力推定値は、終了時点脱分極心起電力推定値の値から不応心電位時系列の対応する時点における心電位の絶対値を引き算した値である。
【0102】
そのため心起電力推定方法は、推定シミュレーション法よりも負荷の小さな演算によって不応心起電力時系列を推定することができる。
【0103】
心起電力推定方法は、分類の結果得られる再分極心電位時系列に基づき、再分極心起電力時系列を推定する。上述したように、再分極心起電力時系列は、再分極心電位時系列の時点ごとに再分極心起電力推定値を有するデータである。上述したように、再分極心起電力推定値は、終了時点不応心起電力推定値から、再分極心電位時系列が示す過去の再分極差分電位の総和を引き算した値である。
【0104】
そのため心起電力推定方法は、推定シミュレーション法よりも負荷の小さな演算によって再分極心起電力時系列を推定することができる。
【0105】
このように、心起電力推定方法は心電位時系列を脱分極期、不応期、再分極期及び休止期に分類し、分類したデータごとに心起電力を推定する。ところで、心電位時系列を再現する数理モデルを用いたシミュレーションにより逆問題を解くことで心起電力を推定する方法は心電位時系列を再現する必要がある。そのためこのような数理モデルを用いたシミュレーションでは、１つの時点の心起電力を推定するために１つの拍動で生じた心電位時系列の全データを必要とする。
【0106】
一方、心起電力推定方法であれば、１つの時点の心起電力を推定するためには必ずしも心電位時系列の全てを必要としない。例えば脱分極期の心起電力の推定は、脱分極期の心起電力のデータさえあれば可能である。このように、心起電力推定方法は、数理モデルを用いたシミュレーションよりも負荷の小さな演算によって心起電力時系列を推定することができる。そのため、心起電力推定方法は心起電力を推定する際の計算負荷を軽減することができる。
【0107】
（変形例）
心起電力推定方法では、心起電力時系列が示す各単位心起電力時系列に対して一致処理を実行してもよい。一致処理は、再分極期の終了時点以降の期間におけるデータが示す心起電力の平均値が脱分極期の開始時点以前のデータが示す心起電力の平均値に略一致であるように、再分極期の終了時点以降の期間におけるデータが示す心起電力を定数倍する処理である。
心起電力推定方法では、心起電力時系列が示す各単位心起電力時系列に対して一致処理を実行してもよい。一致処理は、再分極期の終了時点以降の期間におけるデータが示す心起電力の平均値が脱分極期の開始時点以前のデータが示す心起電力の平均値に略一致であるように、再分極期の終了時点以降の期間におけるデータが示す心起電力を定数倍する処理である。
【0108】
【0109】
図３８は、上段のグラフに心電位時系列Ｄ５０１と心起電力時系列Ｄ５０２とを示す。心起電力時系列Ｄ５０２は、心電位時系列Ｄ５０１に基づき心起電力推定方法によって推定された心起電力時系列であって一致処理が実行されていない心起電力時系列である。図３８の上段のグラフにおける横軸は時刻を表す。図３８の上段のグラフにおける縦軸は電位を表す。図３８の上段のグラフにおける縦軸の単位はミリボルトである。
【0110】
図３８は、下段のグラフに心起電力時系列Ｄ５０３を示す。心起電力時系列Ｄ５０３は、心電位時系列Ｄ５０１に基づき心起電力推定方法によって推定された心起電力時系列であって一致処理の実行後の心起電力時系列である。図３８の下段のグラフにおける横軸は時刻を表す。図３８の下段のグラフにおける縦軸は電位を表す。図３８の下段のグラフにおける縦軸の単位はミリボルトである。
【0111】
図３８の上段のグラフは、心起電力時系列Ｄ５０２における電位が負のピークの電位と基線が示す電位との電位差を示す。具体的には、図３８の上段のグラフにおける電位差Ｗ１、電位差Ｗ２、電位差Ｗ３及び電位差Ｗ４が、心起電力時系列Ｄ５０２における電位が負のピークの電位と基線が示す電位との電位差である。
【0112】
図３８の下段のグラフが示す心起電力時系列Ｄ５０３は、心電位時系列Ｄ５０２に比べて電位が負のピークの電位と基線が示す電位との電位差が小さいことを示す。
【0113】
このように、一致処理が行われない場合、心起電力時系列における不応期の電位の平均値は必ずしも略一定ではない。一方、一致処理の実行により不応期の電位の平均値が略一定の心起電力時系列が取得される。
【0114】
また、図３８は、心起電力時系列Ｄ５０３の基線は心起電力時系列Ｄ５０２の基線よりも心電位時系列Ｄ５０１の基線に対する平行の度合が高いことを示す。
【0115】
このように、一致処理は、心起電力時系列における心電位時系列の基線に対する平行の度合を一致処理が行われない場合よりも高くすることができる。
【0116】
このように一致処理の実行により不応期の電位の平均値が略一定になるため心起電力時系列のグラフはユーザにとって見やすいグラフになる。そのため、ユーザは、一致処理が行われない場合よりも少ない労力で心起電力時系列が示す内容を取得することができる。
【0117】
また、一致処理の実行により、心起電力時系列における心電位時系列の基線に対する平行の度合が高まるため、心起電力時系列のグラフはユーザにとって見やすいグラフになる。そのため、ユーザは、一致処理が行われない場合よりも少ない労力で心起電力時系列が示す内容を取得することができる。
【0118】
なお不応期心起電力時系列を取得する処理は、本明細書に記載の処理に限らない。例えば、不応期心起電力時系列を取得する処理は、不応心電位時系列の各データに対応する心起電力として変形不応心起電力推定値を取得する処理（以下「不応心起電力推定変形処理」という。）であってもよい。例えば、変形不応心起電力推定値は、不応心電位時系列が示す心電位に定数を乗算した値を終了時点脱分極心起電力推定値の値に足し算した値であってもよい。定数は再分極心電位時系列の高さに応じて定まる値である。定数は例えば、Ｔ波最大電位の１００分の１から２００分の１である。定数は、例えば、心電位の大きさによって調整されてもよい。
【0119】
不応心起電力推定変形処理を有する心起電力推定方法によって推定された心起電力時系列が示す不応期における波形が実施形態の心起電力推定方法によって推定された心起電力時系列が示す不応期における波形よりも複雑である。波形が複雑であることは、情報量が多いことを意味する。そのため、不応心起電力推定変形処理を有する心起電力推定方法によって推定された心起電力時系列は、実施形態の心起電力推定方法によって推定された心起電力時系列よりも多くの心筋の動作に関する情報を示す。
【0120】
なお、不応心起電力時系列を推定する方法は、予め機械学習によって不応心電位時系列と不応心起電力時系列との関係を学習済みの学習済みモデルを用いて、不応心起電力時系列を推定する方法であってもよい。このような場合、学習済みモデルの説明変数は不応心電位時系列であり、従属変数は不応心起電力時系列である。機械学習の方法は、再帰型ニューラルネットワークであっても、深層学習であってもよい。
【0121】
なお、脱分極心起電力時系列を推定する方法は、本明細書に記載の方法に限らない。例えば、脱分極心起電力時系列を推定する方法は、予め取得済みの脱分極心電位時系列と脱分極心起電力時系列との間の対応関係を示す全単射の写像を用いて推定する方法であってもよい。
【0122】
なお、再分極期心起電力時系列を推定する方法は、本明細書に記載の方法に限らない。例えば、再分極期心起電力時系列を推定する方法は、予め取得済みの脱分極心電位時系列と脱分極心起電力時系列との間の対応関係を示す全単射の写像を用いて推定する方法であってもよい。
【0123】
なお、ステップＳ２００の処理は、区分ステップの一例である。ステップＳ３０２の処理は、状態時系列取得ステップの一例である。ステップＳ３０３の処理は、脱分極心起電力推定ステップの一例である。ステップＳ３０４の処理は、不応心起電力推定ステップの一例である。ステップＳ３０５の処理は、再分極心起電力推定ステップの一例である。
【0124】
通信部１１又は入力部１３１はそれぞれ取得部の一例である。ステップＳ２００の処理は、区分処理の一例である。ステップＳ３０２の処理は、状態時系列取得処理の一例である。ステップＳ３０３の処理は、脱分極心起電力推定処理の一例である。ステップＳ３０４の処理は、不応心起電力推定処理の一例である。ステップＳ３０５の処理は、再分極心起電力推定処理の一例である。
【0125】
心起電力推定装置１は、ネットワークを介して通信可能に接続された複数台の情報処理装置を用いて実装されてもよい。この場合、心起電力推定装置１が備える各機能部は、複数の情報処理装置に分散して実装されてもよい。
【0126】
なお、心起電力推定装置１の各機能の全て又は一部は、ＡＳＩＣ（Application Specific Integrated Circuit）やＰＬＤ（Programmable Logic Device）やＦＰＧＡ（Field Programmable Gate Array）、ＰＲＣ（Physical Reservoir Computing）等のハードウェアを用いて実現されてもよい。プログラムは、コンピュータ読み取り可能な記録媒体に記録されてもよい。コンピュータ読み取り可能な記録媒体とは、例えばフレキシブルディスク、光磁気ディスク、ＲＯＭ、ＣＤ－ＲＯＭ等の可搬媒体、コンピュータシステムに内蔵されるハードディスク等の記憶装置である。プログラムは、電気通信回線を介して送信されてもよい。
【0127】
以上、この発明の実施形態について図面を参照して詳述してきたが、具体的な構成はこの実施形態に限られるものではなく、この発明の要旨を逸脱しない範囲の設計等も含まれる。
【符号の説明】
【0128】
１００…心起電力推定システム、１…心起電力推定装置、２…心電位測定装置、１０…制御部、１１…通信部、１２…記憶部、１３…ユーザインタフェース
【図1】
【図2】
【図3】
【図4】
【図5】
【図6】
【図7】
【図8】
【図9】
【図10】
【図11】
【図12】
【図13】
【図14】
【図15】
【図16】
【図17】
【図18】
【図19】
【図20】
【図21】
【図22】
【図23】
【図24】
【図25】
【図26】
【図27】
【図28】
【図29】
【図30】
【図31】
【図32】
【図33】
【図34】
【図35】
【図36】
【図37】
【図38】