特開 2024-132728 電流源の推定時刻決定方法、情報処理装置およびプログラム

(19)【発行国】日本国特許庁(JP)

(12)【公報種別】公開特許公報(A)

(11)【公開番号】P2024132728

(43)【公開日】2024-10-01

(54)【発明の名称】電流源の推定時刻決定方法、情報処理装置およびプログラム

(51)【国際特許分類】

   A61B 5/245 20210101AFI20240920BHJP

   A61B 5/372 20210101ALI20240920BHJP

【ＦＩ】

A61B5/245

A61B5/372

【審査請求】未請求

【請求項の数】9

【出願形態】ＯＬ

(21)【出願番号】P 2023043620

(22)【出願日】2023-03-17

(71)【出願人】

【識別番号】000006747

【氏名又は名称】株式会社リコー

(71)【出願人】

【識別番号】504157024

【氏名又は名称】国立大学法人東北大学

(74)【代理人】

【識別番号】110002147

【氏名又は名称】弁理士法人酒井国際特許事務所

(72)【発明者】

【氏名】朝井 都

(72)【発明者】

【氏名】中里 信和

(72)【発明者】

【氏名】菅野 彰剛

【テーマコード（参考）】

4C127

【Ｆターム（参考）】

4C127AA03

4C127AA10

4C127GG05

(57)【要約】

【課題】医師の解析方法により即した形で特徴的な波形情報（ＩＥＤ）の時刻を決定する。
【解決手段】電流源の推定時刻決定方法は、センサにより取得した１つ以上の波形データについて、特徴的な波形が得られる解析区間を決定する解析区間決定ステップと、前記解析区間決定ステップで決定した前記解析区間の信号に対して電流源の推定を実施する電流源推定ステップと、前記電流源推定ステップで推定した前記電流源の座標を用いて、前記解析区間において当該電流源が連続移動する連続移動区間を検出する連続移動区間検出ステップと、前記連続移動区間検出ステップで検出された前記連続移動区間から１つの電流源推定時刻を決定する推定時刻決定ステップと、を含む。
【選択図】図１０

【特許請求の範囲】

【請求項1】

センサにより取得した１つ以上の波形データについて、特徴的な波形が得られる解析区間を決定する解析区間決定ステップと、
前記解析区間決定ステップで決定した前記解析区間の信号に対して電流源の推定を実施する電流源推定ステップと、
前記電流源推定ステップで推定した前記電流源の座標を用いて、前記解析区間において当該電流源が連続移動する連続移動区間を検出する連続移動区間検出ステップと、
前記連続移動区間検出ステップで検出された前記連続移動区間から１つの電流源推定時刻を決定する推定時刻決定ステップと、
を含むことを特徴とする電流源の推定時刻決定方法。

【請求項2】

前記解析区間決定ステップは、特徴的な波形情報が得られる区間として時間幅を設定し、設定した前記時間幅を前記解析区間とする、
ことを特徴とする請求項１に記載の電流源の推定時刻決定方法。

【請求項3】

前記推定時刻決定ステップは、前記電流源の評価指標を考慮して前記電流源推定時刻の決定を行う、
ことを特徴とする請求項１に記載の電流源の推定時刻決定方法。

【請求項4】

前記推定時刻決定ステップは、前記電流源の評価指標として、ＧｏＦ(Goodness of Fit)、correlation、信頼体積や、前記電流源の電流モーメント（intensity）の少なくとも何れか１つを含む、
ことを特徴とする請求項３に記載の電流源の推定時刻決定方法。

【請求項5】

前記推定時刻決定ステップは、前記電流源の評価指標の他、因果性解析結果（Causality）、周波数カップリングの少なくとも何れか１つを利用する、
ことを特徴とする請求項３に記載の電流源の推定時刻決定方法。

【請求項6】

前記連続移動区間検出ステップは、前記電流源の移動量が所定の条件を満たす場合に、前記解析区間の終端の時刻を採用して前記連続移動区間を検出する、
ことを特徴とする請求項１ないし５の何れか一項に記載の電流源の推定時刻決定方法。

【請求項7】

前記連続移動区間検出ステップは、
前記所定の条件として、
前記解析区間における前記電流源の座標の累積移動和が閾値を超えないこと、
前記解析区間の基準時刻の前記電流源の座標を基準として最大移動距離が閾値を超えないこと、
前記解析区間内の前記電流源が同じ脳回内に存在すること、
の少なくとも何れか１つを含む、
ことを特徴とする請求項６に記載の電流源の推定時刻決定方法。

【請求項8】

電流源の推定時刻の検出を行う推定時刻検出部を備えた情報処理装置において、
前記推定時刻検出部は、
センサにより取得した１つ以上の波形データについて、特徴的な波形が得られる解析区間を決定し、
決定した前記解析区間の信号に対して電流源の推定を実施し、
推定した前記電流源の座標を用いて、前記解析区間において当該電流源が連続移動する連続移動区間を検出し、
検出された前記連続移動区間から１つの電流源推定時刻を決定する、
ことを特徴とする情報処理装置。

【請求項9】

センサにより取得した１つ以上の波形データについて、特徴的な波形が得られる解析区間を決定する解析区間決定ステップと、
前記解析区間決定ステップで決定した前記解析区間の信号に対して電流源の推定を実施する電流源推定ステップと、
前記電流源推定ステップで推定した前記電流源の座標を用いて、前記解析区間において当該電流源が連続移動する連続移動区間を検出する連続移動区間検出ステップと、
前記連続移動区間検出ステップで検出された前記連続移動区間から１つの電流源推定時刻を決定する推定時刻決定ステップと、
をコンピュータに実行させるためのプログラム。

【発明の詳細な説明】

【技術分野】

【0001】

本発明は、電流源の推定時刻決定方法、情報処理装置およびプログラムに関する。

【背景技術】

【0002】

従来、脳磁計や脳波計の解析において、特徴的な波形を解析することは非常に重要である。ただし、技術の進歩によって計測時のサンプリング周波数やセンサの数は増加傾向にあり、結果として目視での波形の探索にかかる時間も増加する傾向がある。

【0003】

例えば、臨床における脳磁計を用いたてんかんの診断では、等価電流双極子法と呼ばれる手法を用いててんかん病変部の局在性を評価する。等価電流双極子法では、頭皮上で計測された磁場を生み出す電流源（ダイポール）を推定する。ダイポールを推定するためには、複数センサの時系列の中から、特徴的な波形情報（ＩＥＤ：Interictal Epileptiform Discharge）が発生した時刻（起始部）とその波形情報が出現しているセンサの絞り込みが必要となる。

【0004】

現状では、医師は手動でＩＥＤの探索・起始部の決定を行っているが、脳磁計のデータは膨大であるため、個別のＩＥＤについて手動で正確にセンサおよびＩＥＤの時刻の抽出を行うのは困難である。

【0005】

また、医師が解析を行う際には、ＩＥＤの起始部を等価電流双極子法で解析することが理想ではあるが、起始部はＳ／Ｎが低いことが知られている。そのため、医師が解析を行う際には、実際には起始部からＩＥＤのピークまでの中でＩＥＤの時刻を決定している。

【0006】

しかしながら、電流源を推定する時刻の決定方法は、学会が出している指針（非特許文献８）やガイドライン（非特許文献９）においても厳密には定められておらず、医師がそれぞれの判断によって手動で決定しているため、ばらつきが生じてしまっている。

【0007】

そこで、特許文献１には、特徴的な波形情報（ＩＥＤ）が現れている時刻の決定とセンサの抽出をより正確に行う目的で、ＩＥＤの確率マップを用いて、等価電流双極子法に必要なＩＥＤの時刻の検出とセンサの選択を行うシステムが開示されている。

【発明の概要】

【発明が解決しようとする課題】

【0008】

しかしながら、特許文献１に開示の技術によれば、電流源を推定する時刻の決定方法が必ずしも医師の判断基準に即したものにならない、という問題があった。

【0009】

本発明は、上記に鑑みてなされたものであって、医師の解析方法により即した形で特徴的な波形情報（ＩＥＤ）の時刻を決定することを目的とする。

【課題を解決するための手段】

【0010】

上述した課題を解決し、目的を達成するために、本発明の電流源の推定時刻決定方法は、センサにより取得した１つ以上の波形データについて、特徴的な波形が得られる解析区間を決定する解析区間決定ステップと、前記解析区間決定ステップで決定した前記解析区間の信号に対して電流源の推定を実施する電流源推定ステップと、前記電流源推定ステップで推定した前記電流源の座標を用いて、前記解析区間において当該電流源が連続移動する連続移動区間を検出する連続移動区間検出ステップと、前記連続移動区間検出ステップで検出された前記連続移動区間から１つの電流源推定時刻を決定する推定時刻決定ステップと、を含むことを特徴とする。

【発明の効果】

【0011】

本発明によれば、医師の解析方法により即した形で特徴的な波形情報（ＩＥＤ）の時刻を決定することができる。

【図面の簡単な説明】

【0012】

【図1】図１は、第１の実施形態に係る生体信号計測システムの構成を示す概略図である。

【図2】図２は、サーバの機能構成を概略的に示すブロック図である。

【図3】図３は、情報処理装置のハードウェア構成の一例を示す図である。

【図4】図４は、情報処理装置の機能ブロック構成の一例を示す図である。

【図5】図５は、推定時刻検出処理の流れを示すフローチャートである。

【図6】図６は、機械学習を用いたＩＥＤの確率マップの生成の際の学習処理の流れを示すフローチャートである。

【図7】図７は、ＩＥＤの確率マップの一例を示す図である。

【図8】図８は、閾値処理後のＩＥＤ確率マップの一例を示す図である。

【図9】図９は、センサ数の拡張方法の一例を示す概念図である。

【図10】図１０は、推定時刻検出処理の流れを示すフローチャートである。

【図11】図１１は、解析区間におけるダイポールの軌跡の一例を示す図である。

【図12】図１２は、解析区間におけるダイポールの軌跡を脳スライス画像上に重畳した表示の一例を示す図である。

【図13】図１３は、第２の実施形態に係る閾値設定入力ＵＩの一例を示す図である。

【発明を実施するための形態】

【0013】

以下に添付図面を参照して、電流源の推定時刻決定方法、情報処理装置およびプログラムの実施の形態を詳細に説明する。また、以下の実施形態によって本発明が限定されるものではなく、以下の実施形態における構成要素には、当業者が容易に想到できるもの、実質的に同一のもの、およびいわゆる均等の範囲のものが含まれる。さらに、以下の実施の形態の要旨を逸脱しない範囲で構成要素の種々の省略、置換、変更および組み合わせを行うことができる。

【0014】

本実施形態は、電流源推定時刻の決定に際して、以下の特徴を有する。要するに、特徴的な波形を有する区間を検出し、検出区間の中で電流源（ダイポール）が連続的な動きをしている区間を探索する。これにより、特徴的な波形の活動を捉えた区間を探索することができ、この区間内で始まり付近のダイポールを検出することで、特徴的な波形の活動について活動の起始部部分の時刻を自動で検出することができることが特徴になっている。

【0015】

医師が解析を行う際には、ＩＥＤの起始部を等価電流双極子法で解析することが理想ではあるが、起始部はＳ／Ｎが低いことが知られているため、実際には起始部からＩＥＤのピークまでの中でＩＥＤの時刻を決定している。本実施形態では起始部に限定しないで解析を行う。

【0016】

［第１の実施形態］
（生体信号計測システムの概略）
図１は、第１の実施形態に係る生体信号計測システムの構成を示す概略図である。図１を参照しながら、本実施形態に係る生体信号計測システム１の概略について説明する。

【0017】

生体信号計測システム１は、特定の発信源（生体部位）からの被検者の複数種類の生体信号（例えば、脳磁（ＭＥＧ：Magneto-encephalography）信号、および脳波（ＥＥＧ：Electro-encephalography）信号を計測し、表示するシステムである。なお、本発明において、測定対象となる生体信号は、脳磁信号および脳波信号に限られるものではなく、例えば、心臓の活動に応じて発生する電気信号（心電図として表現可能な電気信号）などであってもよい。

【0018】

図１に示すように、生体信号計測システム１は、被検者の１以上の生体信号を測定する測定装置３と、測定装置３で測定された１種類以上の生体信号を記録するサーバ４０と、サーバ４０に記録された１種類以上の生体信号を解析する生体信号表示装置である情報処理装置５０と、を含む。測定装置３は、例えば脳磁場や刺激を与えたタイミング等の計測値を収集する脳磁計である。なお、図１では、サーバ４０と情報処理装置５０とが別々に記載されているが、例えば、サーバ４０が有する機能の少なくとも一部が情報処理装置５０に組み込まれる形態であってもよい。

【0019】

図１の例では、被検者（被測定者）は、頭に脳波測定用の電極（またはセンサ）を付けた状態で測定テーブル４に仰向けで横たわり、測定装置３のデュワ３１の窪み３２に頭部を入れる。デュワ３１は、液体ヘリウムを用いた極低温環境の保持容器であり、デュワ３１の窪み３２の内側には脳磁測定用の多数の磁気センサが配置されている。測定装置３は、電極からの脳波信号と、磁気センサからの脳磁信号とを収集し、収集した脳波信号および脳磁信号を含むデータ（以下、「測定データ」と称する場合がある）をサーバ４０へ出力する。サーバ４０へ出力された測定データは、情報処理装置５０に読み出されて表示され、解析される。一般的に、磁気センサを内蔵するデュワ３１および測定テーブル４は、磁気シールドルーム内に配置されているが、図１では便宜上、磁気シールドルームの図示を省略している。

【0020】

情報処理装置５０は、複数の磁気センサからの脳磁信号の波形データと、複数の電極からの脳波信号の波形データを、同じ時間軸上に同期させて表示する装置である。脳波信号とは、神経細胞の電気的な活動（シナプス伝達の際のニューロンの樹状突起で起きるイオン電荷の流れ）を電極間の電圧値として表される信号である。脳磁信号とは、脳の電気活動により生じた微小な電場変動を表す信号である。脳磁場は、高感度の超伝導量子干渉計（ＳＱＵＩＤ：Superconducting Quantum Interference Device）センサで検知される。これらの脳波信号および脳磁信号は、「生体信号」の一例である。

【0021】

ここで、図２はサーバ４０の機能構成を概略的に示すブロック図である。図２に示すように、サーバ４０は、データ取得部４１と、データ記憶部４２とを有する。

【0022】

データ取得部４１は、測定装置３から測定データを定期的に取得する。ここで、測定データは、測定装置３のデュワ３１の複数の磁気センサで計測した個々の波形データである。

【0023】

データ記憶部４２は、測定装置３から取得した測定データを記憶する。

【0024】

（情報処理装置のハードウェア構成）
図３は、情報処理装置５０のハードウェア構成の一例を示す図である。図３を参照しながら、情報処理装置５０のハードウェア構成について説明する。

【0025】

図３に示すように、情報処理装置５０は、ＣＰＵ（Central Processing Unit）１０１と、ＲＡＭ（Random Access Memory）１０２と、ＲＯＭ（Read Only Memory）１０３と、補助記憶装置１０４と、ネットワークＩ／Ｆ１０５と、入力装置１０６と、表示装置１０７と、を有し、これらがバス１０８で相互に接続されている。

【0026】

ＣＰＵ１０１は、情報処理装置５０の全体の動作を制御し、各種の情報処理を行う演算装置である。ＣＰＵ１０１は、ＲＯＭ１０３または補助記憶装置１０４に記憶された情報表示プログラムを実行して、測定収集画面および解析画面の表示動作を制御する。

【0027】

ＲＡＭ１０２は、ＣＰＵ１０１のワークエリアとして用いられ、主要な制御パラメータおよび情報を記憶する揮発性の記憶装置である。ＲＯＭ１０３は、基本入出力プログラム等を記憶する不揮発性の記憶装置である。例えば、上述の情報表示プログラムがＲＯＭ１０３に記憶されているものとしてもよい。

【0028】

補助記憶装置１０４は、ＨＤＤ（Hard Disk Drive）またはＳＳＤ（Soliｄ State Drive）等の記憶装置である。補助記憶装置１０４は、例えば、情報処理装置５０の動作を制御する制御プログラム、ならびに、情報処理装置５０の動作に必要な各種のデータおよびファイル等を記憶する。

【0029】

ネットワークＩ／Ｆ１０５は、サーバ４０等のネットワーク上の機器と通信を行うための通信インターフェースである。ネットワークＩ／Ｆ１０５は、例えば、ＴＣＰ（Transmission Control Protocol）／ＩＰ（Internet Protocol）に準拠したＮＩＣ（Network Interface Card）等によって実現される。

【0030】

入力装置１０６は、タッチパネルの入力機能、キーボード、マウスおよび操作ボタン等のユーザインターフェース等である。表示装置１０７は、各種の情報を表示するディスプレイ装置である。表示装置１０７は、例えば、タッチパネルの表示機能、液晶ディスプレイ（ＬＣＤ：Liquid Crystal Display）または有機ＥＬ（Electro－Luminescence）等によって実現される。表示装置１０７は、測定収集画面および解析画面を表示し、入力装置１０６を介した入出力操作に応じて画面が更新される。

【0031】

なお、図３に示す情報処理装置５０のハードウェア構成は一例であり、これ以外の装置が備えられるものとしてもよい。また、図３に示す情報処理装置５０は、例えば、ＰＣ（Personal Computer）を想定したハードウェア構成であるが、これに限定されるものではなく、タブレット等のモバイル端末であってもよい。この場合、ネットワークＩ／Ｆ１０５は、無線通信機能を有する通信インターフェースであればよい。

【0032】

（情報処理装置の機能ブロック構成）
図４は、情報処理装置５０の機能ブロック構成の一例を示す図である。図４を参照しながら、情報処理装置５０の機能ブロック構成について説明する。

【0033】

図４に示すように、情報処理装置５０は、前処理部５０１と、ＩＥＤ確率マップ算出部５０２と、閾値処理部５０３と、後処理部５０４と、ピーク検出部５０５と、推定時刻検出部５０６と、を備える。

【0034】

前処理部５０１は、センサの抽出・拡張、ダウンサンプリングや周波数フィルタの適用、アーチファクト除外、不良チャネル処理、時間窓切り出し、磁場データの標準化などの前処理を行う。

【0035】

センサの抽出・拡張は、後述のＩＥＤ確率マップ算出の際に、全センサを使って算出する方法以外にも、予めグルーピングした一定数センサ群のみを使用して算出することも可能である。センサのグルーピングには、側頭葉や前頭葉など解剖学的な基準に即して設定したものや、単純に近くにある任意の数のセンサをグルーピングしたものなどが考えられる。また、後述の学習過程において、使用したセンサ数よりも、算出の際のセンサ数が少ない場合には、架空のセンサを拡充することも可能である。

【0036】

ダウンサンプリングは、学習時に使用したサンプリング周波数に合わせる目的で適用される。周波数フィルタも同様で学習時に適用したフィルタリングと同じものを適用する。よく使われるフィルタは、３５Ｈｚのローパスフィルタ、３Ｈｚ－３５Ｈｚのバンドパスフィルタなどである。

【0037】

アーチファクト除外は、心拍性のアーチファクトや瞬き・体動によるアーチファクトを除外する目的で、ＩＣＡ（非特許文献１を参照）やＤＳＳＰ（非特許文献２を参照）などが適用される。

【0038】

不良チャネル処理は、予め設定した閾値を超えるような磁場変化が観測されたセンサを除外したり、周囲のセンサ値を使用して補間したりする処理を指す。

【0039】

時間窓の切り出し方は、時間窓の長さ分だけずらしていくオーバーラップなしのものや、時間窓の半分の長さをオーバーラップさせる方法、１／４の長さをオーバーラップさせる方法などがある。オーバーラップさせる場合には、後述のＩＥＤ確率マップ算出の際にオーバーラップ分の加算平均を行う。

【0040】

磁場データの標準化は、切り出した時間窓の中で平均が０、分散が１になるような標準化が適用される。標準化の他に、予め設定した磁場の範囲内を－１～１となるように正規化する手法などを用いてもよい。

【0041】

ＩＥＤ確率マップ算出部５０２は、特徴的な波形（ＩＥＤ：Interictal Epileptiform Discharge）の確率マップを算出する。

【0042】

ＩＥＤとは、てんかん病変部に特徴的なSpike Wave、Spike and Wave、Poly Spike and Wave、Sharp Waveなどの波形を指す。本実施形態では、ＩＥＤ確率マップ算出部５０２は、事前に学習した機械学習のモデルを用いて確率マップを算出する。機械学習のモデルの作成方法については後述する。

【0043】

閾値処理部５０３は、閾値用いて、ＩＥＤ確率マップ算出部５０２において得られたＩＥＤの確率マップからＩＥＤの確率が高い時間とセンサの領域の絞り込みを行う。

【0044】

後処理部５０４は、閾値処理部５０３で閾値処理したＩＥＤの確率マップからダイポール推定に使用するためのサンプリング時点およびセンサの抽出を行う後処理を実行する。また、後処理部５０４は、ピーク検出手法やマップが複数のピークを持つ場合（病変部が２つ存在する場合などにみられるケース）には、マップを分離するための後処理を行う。

【0045】

ピーク検出部５０５は、ピーク検出部５０５以前の処理で得られたサンプリング地点を用いてピーク検出を行う。

【0046】

最後に、推定時刻検出部５０６は、ピーク検出部５０５で得られたピーク時刻を用いて、電流源推定時刻の検出を行う。

【0047】

次に、推定時刻検出処理の流れについて説明する。

【0048】

ここで、図５は推定時刻検出処理の流れを示すフローチャートである。処理に先立ち、測定装置３は、脳磁計計測を実施し、デュワ３１の複数の磁気センサで計測した個々の波形データをサーバ４０に出力する。サーバ４０は、測定装置３のデュワ３１の複数の磁気センサで計測した個々の波形データをデータ記憶部４２に記憶する。

【0049】

そして、図５に示すように、まず、前処理部５０１は、測定装置３のデュワ３１の複数の磁気センサで計測した個々の波形データをサーバ４０のデータ記憶部４２から取得する（ステップＳ１）。

【0050】

次に、前処理部５０１は、取得したデータに対し、センサの抽出・拡張、ダウンサンプリングや周波数フィルタの適用、アーチファクト除外、不良チャネル処理、時間窓切り出し、磁場データの標準化などの前処理を行う（ステップＳ２）。

【0051】

次に、ＩＥＤ確率マップ算出部５０２は、特徴的な波形（ＩＥＤ）の確率マップを算出する（ステップＳ３）。

【0052】

ＩＥＤの確率マップの算出については、機械学習を用いて算出したモデルを適用する方法のほかにも、従来のスパイク位置の検出アルゴリズム（非特許文献４を参照）を用いたものでもよい。機械学習を用いる場合には、医師が手動で検出したＩＥＤの時点を正解として学習させるため、従来のスパイク位置を検出する手法と比較して、医師が解析した場合と類似した結果が得られるという効果がある。従来の手法では、ＩＥＤの起始部もしくはＩＥＤのピークを検出することを目的として開発されているが、実際の解析では、医師は波形や所見などからダイポール推定を行う時刻を起始部からピークまでの間で調整しながら解析している。

【0053】

ここで、ＩＥＤ確率マップ算出部５０２が、ＩＥＤの確率マップの生成において機械学習を用いる場合について説明する。

【0054】

図６は、機械学習を用いたＩＥＤの確率マップの生成の際の学習処理の流れを示すフローチャートである。図６に示すように、ＩＥＤ確率マップ算出部５０２は、前処理部５０１で前処理を実行した前処理済み脳磁図データを取得すると（ステップＳ１１）、前処理済み脳磁図データに対して、ＩＥＤの確率マップを出力するための学習を行う（ステップＳ１２）。なお、学習の際に用いる正解データには、ダイポール推定を行った際に記録されているサンプリング時点とその時使用したセンサの情報を用いる。

【0055】

具体的には、非特許文献５に記載のＵ－Ｎｅｔに代表されるSemantic Segmentationというタスクでよく使用されているネットワークを応用することで実現される。Semantic Segmentationは、入力となる全配列の各要素についてラベリングを行うタスクのことで、医療画像においては、脳腫瘍や癌の領域を自動で推定する際に用いられる。ここでは、医師がダイポール推定を行った際に選択されたセンサと推定を行った時刻から作成したマスクデータを正解として、入力の脳磁図データがマスクデータとなるように学習を行う。このとき、汎化性能を高める目的で、センサの順番をランダムに並び替えたり、特定のセンサのある範囲の値を０にしたりといったデータ拡張手法を適用してもよい。

【0056】

また、Semantic Segmentation以外のタスクであるObject DetectionやInstance Segmentationといったタスクで使われるネットワークでも同様に実現可能である。Object Detectionは、入力となる配列中のどこに目的となるobjectが存在しているかを矩形などで推定する手法で、この手法を使用する場合には、ＩＥＤ確率マップは算出できないが、検出された矩形からダイポール推定を行う時刻（例えば矩形の中点）とセンサを直接取得することができる。Instance Segmentationは、Semantic SegmentationとObject Detectionを混ぜたようなタスクで、検出したObjectの分類と個数に加えてそのObjectの領域も算出することが可能である。Instance Segmentationでは、ＩＥＤの確率マップを算出することができるので、そのまま応用可能である。

【0057】

次いで、ＩＥＤ確率マップ算出部５０２は、ステップＳ２の学習が終わると、学習済みモデルを作成する（ステップＳ１３）。作成された学習済みモデルは、機械学習の推論、つまり未知データに対してＩＥＤを検出する時に用いられる。

【0058】

ここで、ＩＥＤ確率マップ算出部５０２が算出するＩＥＤの確率マップについて説明する。

【0059】

図７は、ＩＥＤの確率マップの一例を示す図である。図７においては、横軸が時間、縦軸はセンサのインデックスを表している。図７に示す例の場合、時間のインデックスが１００、センサのインデックスが４０～８０あたりにＩＥＤが出現したと推定される。なお、図７においては、ＩＥＤの確率に応じて白黒の濃淡で表現したが、ＩＥＤの確率に応じて色が変化するようなカラー表示でもよい。

【0060】

図５に戻り、次に、閾値処理部５０３は、ＩＥＤ確率マップ算出部５０２において得られたＩＥＤの確率マップに対して閾値処理を行う（ステップＳ４）。具体的には、閾値処理部５０３は、閾値用いて、ＩＥＤ確率マップ算出部５０２において得られたＩＥＤの確率マップからＩＥＤの確率が高い時間とセンサの領域の絞り込みを行う。

【0061】

なお、閾値処理部５０３においては、閾値を高くすればするほど、よりＩＥＤらしい波形のみを選択できるが、ダイポール推定を行うＩＥＤそのものの数が少なくなってしまう。一方、閾値処理部５０３においては、閾値を低くした場合、誤検出は増えるが、より多くのＩＥＤを検出することができる。さらに、脳磁図のデータは個人差の影響が強いため、閾値を下げなければ、ＩＥＤがほとんど検出できない可能性もある。

【0062】

そこで、閾値処理部５０３においては、ＵＩ（User Interface）と連携させてインタラクティブに検出するＩＥＤの数を制御する構成としてもよい。

【0063】

閾値処理部５０３は、閾値として、基本的には「０．８」などのプリセット値を用いる。

【0064】

ここで、図８は閾値処理後のＩＥＤ確率マップの一例を示す図である。図８に示す例においては、閾値処理部５０３は、閾値として「０．５」を用いた。図８に示す２次元マップであるＩＥＤ確率マップ６０５では、５０％以上の確率でＩＥＤが存在している箇所を抽出したことに相当する。

【0065】

図５に戻り、次に、後処理部５０４は、閾値処理部５０３で閾値処理したＩＥＤ確率マップ６０５からダイポール推定に使用するためのサンプリング時点およびセンサの抽出を行う後処理を実行する（ステップＳ５）。

【0066】

具体的には、閾値処理部５０３を適用すると、所定の確率以上の値を持つ時刻とセンサの領域とが抽出される。ピーク検出、もしくはＩＥＤの検出時点を明確にする目的で、抽出されたセンサの領域に対して、ガウシアンフィルタ等のフィルタをピーク検出の前に適用するようにしてもよい。抽出されたセンサの領域に対して、ガウシアンフィルタ等のフィルタをピーク検出の前に適用した後にピーク検出を行うことで、ＩＥＤの時点が決定される。ピーク検出を行う際には、抽出されたセンサの領域を含む確率マップをセンサ方向に加算平均した１次元の確率マップに対して行ってもよい。加算平均を行わない場合は、センサごとにピーク検出を行い、その時点を平均することでＩＥＤの時刻を決定するようにしてもよい。

【0067】

閾値処理部５０３は、検出されたピーク時点における所定の確率以上の値をもつセンサをダイポール推定に使用するセンサとして抽出する。このとき、ダイポール推定を行う際のセンサの数が少なくなると、ダイポール推定解の安定性が低くなることが知られている。そのためここでは、選択されたセンサ数が規定の数に満たない場合には、そのＩＥＤ時点ではダイポール推定を行わないようにすることも可能である。

【0068】

また、別のセンサ抽出方法として、予めセンサのグループを定義し、抽出されたセンサをそのセンサが属するグループのすべてのセンサへと拡張することでダイポール推定解の安定性を高めることもできる。ここでグループとは、側頭葉付近や前頭葉付近といった解剖学的な基準と照らして設定することを基本とするが、単純に任意の数を含む近傍領域を一つのグループとしてもよい。

【0069】

ここで、図９はセンサ拡張方法の一例を示す概念図である。図９（ａ）は、被検者の頭部２００と、頭部２００を覆うように配置された脳磁図のセンサ２１０と、予め設定されたセンサ２１０のグループ２２０と、を表している。図９（ｂ）において黒丸で表されたセンサ２１０ａは、閾値処理された段階で抽出されたセンサを表している。これらの閾値処理された段階で抽出されたセンサ２１０ａのみを用いてダイポール推定を行うことも可能だが、前述の通り、ダイポール推定解の安定性を高めるために、センサの拡張を実施することも可能である。図９（ｃ）は、センサ数の拡張を行った結果である。すなわち、図９（ｃ）に示すように、閾値処理された段階で抽出されたセンサ２１０ａをそのセンサ２１０ａが属するグループ２２０のすべてのセンサ２１０へと拡張する。これにより、ダイポール推定解の安定性を高めることができる。

【0070】

また、後処理部５０４は、ＩＥＤ確率マップが複数のピークを持つ場合（病変部が２つ存在する場合などにみられるケース）には、マップを分離するための後述の後処理を行うこともできる。

【0071】

近傍領域の設定の別の例として、選択されたセンサ座標から凸包を構成し、次に凸包を所定の閾値でモルフォロジー変換を行って空間を拡張し、拡張された空間に含まれるセンサを一つのグループとするなどが挙げられる。

【0072】

また、選択センサデータの磁場の最大値を持つセンサと（湧き出し）、最小値を持つセンサ（吸い込み）の空間座標の中点を基準に最近傍の点から順に取得してもよい。

【0073】

また、センサデータから取得した磁場情報から、所定の磁場閾値で空間的な等磁場線図を再構成し、選択されたセンサの磁場情報から、既に選択されたセンサの磁場情報の近傍でかつ同じ等磁場線に属するセンサを選択しても良い。

【0074】

また、近傍領域の磁場の空間分布を用いた決定手段と解剖的な基準との組み合わせによるグループを構成することによって、空間的かつ解剖的な視点でのセンサグループを構成しても良い。

【0075】

ピーク検出部５０５は、以前の処理で得られたサンプリング地点を用いてピーク検出を行う（ステップＳ６）。例えば、ピーク検出部５０５は、ノイズの中に埋もれている特定のパターンの波形を、機械学習を用いたパターンマッチングなどによりピーク候補を検出する機能を有している。ピーク検出部５０５は、ピーク候補を検出する際に、決められた手順条件によって定量化したスコアが用いてピークを決定する機能を有している。本実施形態では、「スコア」は、パターンマッチングの強さ、２つの波形の適合率といった概念を定量化したものとして利用する。

【0076】

最後に、推定時刻検出部５０６は、ピーク検出部５０５で得られたピーク時刻を用いて、電流源推定時刻の検出を行う（ステップＳ７）。

【0077】

次に、推定時刻検出部５０６が推定時刻を検出する流れを説明する。ここで、図１０は推定時刻検出処理の流れを示すフローチャートである。

【0078】

図１０に示すように、推定時刻検出部５０６による推定時刻検出の処理フローは、解析区間決定ステップＳ２１と、ダイポール推定ステップＳ２２と、連続移動区間検出ステップＳ２３と、推定時刻決定ステップＳ２４と、で構成される。

【0079】

まず、解析区間決定ステップＳ２１では、推定時刻検出部５０６は、ピーク検出部５０５で検出した波形のピーク時刻を終点として、特徴的な波形が得られる区間として固定幅の時間を設定し、解析区間とする。解析区間の時間幅は、スパイク波形の起始部から頂点までの持続時間の統計的情報（非特許文献１０）などを考慮して設定することが望ましい。その他にも、医師や技師が所望の時間幅を設定できるようにしてもよい。

【0080】

次に、ダイポール推定ステップ（電流源推定ステップ）Ｓ２２では、推定時刻検出部５０６は、解析区間決定ステップＳ２１で決定した区間の信号に対して、全てのサンプリング点でダイポール推定（非特許文献３を参照）を行う。ここでは、ダイポール推定以外にも、最小ノルム法（非特許文献６を参照）やLCMV Beamformer法（非特許文献７を参照）といった空間フィルタ法による電流源解析法で代替することが可能である。空間フィルタ法は、予め脳内に数万程度のダイポールを配置し、各ダイポールのおける電流の経時変化を求めることができる手法である。

【0081】

次に、連続移動区間検出ステップＳ２３では、推定時刻検出部５０６は、ダイポール推定ステップＳ２２で推定したダイポール（電流源）の座標を用いて、解析区間においてピーク時刻を基準としてダイポールの移動量が小さい区間（連続して移動する区間）を抽出する。

【0082】

ここで、図１１は解析区間におけるダイポールの軌跡の一例を示す図、図１２は解析区間におけるダイポールの軌跡を脳スライス画像上に重畳した表示の一例を示す図である。例えば図１１では、解析区間の全サンプリング点でダイポール推定を実施したダイポールの座標例を示している。また、図１１のダイポールを脳スライス画像上に重畳した例が図１２である。

【0083】

図１１および図１２について、ピークから起始部への時間方向に向けてダイポールの動き方に注目すると、ピーク時点からある一定区間は１サンプルあたりの移動量が小さく、ダイポールが連続的な動きをしている。その後、ダイポールの位置が大きく変化し、連続的な動きをしていることがわかる。この例では、推定時刻検出部５０６は、ピーク時点からダイポールの位置が大きく変化する直前の区間を、連続移動区間とする。

【0084】

連続的に移動する区間を抽出する方法としては、次の方法が例示できる。
（１）１つ前のサンプリング点とのダイポールの距離が閾値を超えない区間を抽出する方法
（２）ある一定サンプリング幅のダイポールの移動距離が閾値を超えない区間を抽出する方法
（３）ピーク時刻を起点とした累積移動和が閾値を超えない区間を抽出する方法

【0085】

なお、収録状況によりサンプリング周波数が異なることが考えられるため、上記の手法に関してサンプリング間隔を正規化し、ある一定の時刻あたりのダイポールの移動距離などを使用してもよい。

【0086】

最後に、推定時刻決定ステップＳ２４では、推定時刻検出部５０６は、連続移動区間の中から最終的な１つのダイポール推定時刻（電流源推定時刻）を決定する。

【0087】

連続区間の中からダイポール推定時刻を決める方法としては、次のようなものが例示できる。
（１）連続移動区間の先頭時刻を推定時刻とする方法
（２）連続推定区間の先頭から固定時間後ろの時刻を推定時刻とする方法
（３）連続推定区間の先頭から固定割合後ろ（例えば連続推定区間の１／３時点など）の時刻を推定時刻とする方法
（４）連続推定区間の末尾（ピーク時刻）から固定時間前の時刻を推定時刻とする方法
（５）連続推定区間の末尾から固定割合前（例えば連続推定区間の２／３時点など）の時刻を推定時刻とする方法

【0088】

以上のフローにより、推定時刻検出部５０６は、特徴的な波形の活動を捉えた区間を探索し、区間の中で開始時刻を選択することで、特徴的な波形の活動について、活動の起始部に相当する部分の時刻を自動で検出することができる。

【0089】

このように本実施形態によれば、センサにより取得した波形データから特徴的な波形情報（ＩＥＤ）を有する解析区間を検出し、解析区間の中で電流源（ダイポール）が連続的な動きをしている区間を探索する。これにより、特徴的な波形情報（ＩＥＤ）の起始部を考慮して電流源を推定する時刻を決定することになるので、医師の解析方法により即した形で特徴的な波形情報（ＩＥＤ）の時刻を決定することができる。

【0090】

［第２の実施形態］
次に、第２の実施形態について説明する。

【0091】

第２の実施形態に係る電流源の推定時刻決定方法について、第１の実施形態に係る電流源の推定時刻決定方法と相違する点を中心に説明する。以下、第２の実施の形態の説明では、第１の実施の形態と同一部分の説明については省略し、第１の実施の形態と異なる箇所について説明する。

【0092】

本実施形態では、第１の実施形態の推定時刻検出部５０６が実行する推定時刻決定ステップＳ２４において、推定したダイポールを評価指標に基づいて取捨選択する。

【0093】

より詳細には、本実施形態の推定時刻検出部５０６は、ダイポール推定の近似の良好さを示すＧｏＦ(Goodness of Fit)、correlation、信頼体積(Confidence Volume)、推定したダイポールの電流モーメント（intensity）の少なくとも何れか１つを含むダイポール推定の評価指標について閾値を設定する。本実施形態の推定時刻検出部５０６は、前段で取得した連続移動区間の中で、各時刻のダイポールの評価指標が閾値以内である時刻をさらに抽出し、抽出した時刻の中で最終的な推定時刻を選択する。

【0094】

推定時刻検出部５０６は、ダイポール推定の結果の妥当性を評価するために、ダイポール推定の評価指標の他、時系列データやネットワークの因果性解析結果（Causality）（非特許文献１１）や異なる周波数帯での同期活動を評価する周波数カップリング（非特許文献１２）などを利用してもよい。

【0095】

一般的に起始部はＳＮ比が低いことが知られており、ダイポール推定結果の信頼度が低い可能性がある。本実施形態の推定時刻検出部５０６によれば、ダイポール推定の評価指標の閾値を設けることにより、ダイポール推定の信頼度を考慮に入れることができる、という効果がある。

【0096】

ここで、図１３は第２の実施形態に係る閾値設定入力ＵＩの一例を示す図である。図１３に示すＧＵＩ（Graphical User Interface）画面１００には、ダイポール推定の評価指標の閾値として、ＧｏＦ、correlation、推定したダイポールの電流モーメント（intensity）、信頼体積(Confidence Volume)のそれぞれの閾値を設定することが可能になっている。このように図１３に示すＧＵＩ画面１００と連携することで、医師や技師がダイポール推定の評価指標の閾値として任意の数値を入力することで、入力に応じた結果を提示することができる。

【0097】

［第３の実施形態］
次に、第３の実施形態について説明する。

【0098】

第３の実施形態に係る電流源の推定時刻決定方法について、第１の実施形態から第２の実施形態に係る電流源の推定時刻決定方法と相違する点を中心に説明する。以下、第３の実施の形態の説明では、第１の実施の形態または第２の実施の形態と同一部分の説明については省略し、第１の実施の形態または第２の実施の形態と異なる箇所について説明する。

【0099】

本実施形態では、第１の実施形態または第２の実施形態の推定時刻検出部５０６における連続移動区間検出ステップＳ２３において、推定したダイポールの移動量が所定の条件を満たす場合（推定したダイポールの移動量が少ない場合）に、推定時刻にピーク検出部５０５で検出した解析区間の終端の時刻を採用して連続移動区間を検出する。

【0100】

ダイポールの移動量が少ないと判断する具体的な方法としては、下記の方法が考えられる。
（１）解析区間におけるダイポール座標の累積移動和が閾値を超えない範囲であるか否かで判断する方法
（２）解析区間の基準時刻のダイポールの座標を基準として最大移動距離が閾値を超えない範囲であるか否かで判断する方法
（３）解析区間内のダイポールが同じ脳回内に存在するか否かで判断する方法

【0101】

本実施形態によれば、推定時刻検出部５０６における連続移動区間検出ステップＳ２３において、よりＳＮ比が良い時点でダイポール推定を行う方がダイポール推定結果の信頼性が高くなるため、起始部での推定位置とピークでの推定位置の差異が小さい場合は、ピーク時刻でダイポール推定を行うことでより信頼性が高い結果が得られる、という効果がある。

【0102】

なお、上述の各実施形態において、生体信号計測システム１の各機能部の少なくとも何れかがプログラムの実行によって実現される場合、そのプログラムは、ＲＯＭ等に予め組み込まれて提供される。また、上述の実施形態に係る生体信号計測システム１で実行されるプログラムは、インストール可能な形式または実行可能な形式のファイルでＣＤ－ＲＯＭ、フレキシブルディスク（ＦＤ）、ＣＤ－Ｒ（Compact Disc Recordable）、ＤＶＤ（Digital Versatile Disc）等のコンピュータで読み取り可能な記録媒体に記録して提供するように構成してもよい。

【0103】

また、上述の各実施形態の生体信号計測システム１で実行されるプログラムを、インターネット等のネットワークに接続されたコンピュータ上に格納し、ネットワーク経由でダウンロードさせることにより提供するように構成してもよい。

【0104】

また、上述の各実施形態の生体信号計測システム１で実行されるプログラムを、インターネット等のネットワーク経由で提供または配布するように構成してもよい。また、上述の実施形態の生体信号計測システム１で実行されるプログラムは、上述した各機能部のうち少なくとも何れかを含むモジュール構成となっており、実際のハードウェアとしてはＣＰＵがＲＯＭ等からプログラムを読み出して実行することにより、上述の各機能部が主記憶装置上にロードされて生成されるようになっている。

【符号の説明】

【0105】

５０ 情報処理装置
５０６ 推定時刻検出部

【先行技術文献】

【特許文献】

【0106】

【特許文献1】特開２０２１－０６９９２９号公報

【非特許文献】

【0107】

【非特許文献1】E. Javier, H. Roberto, A. Daniel, F. Alberto, and L. C. Miguel, “Artifact removal in magnetoencephalogram background activity with independent component analysis,” IEEE Trans Biomed Eng, vol. 54, no. 11, pp. 1965-1973, 2007.

【非特許文献2】K. Sekihara, Y. Kawabata, S. Ushio, S. Sumiya, S. Kawabata, Y. Adachi, and S. S. Nagarajan, “Dual signal subspace projection (DSSP): a novel algorithm for removing large interference in biomagnetic measurements,” Journal of Neural Engineering, vol. 13, no. 3, p. 036007, 2016.

【非特許文献3】M. Scherg, “Fundamentals of dipole source potential analysis,” in Auditory Evoked Magnetic Fields and Potentials, M. Hoke, F. Grandori, and G. L. Romani, Eds. Basel, Switzerland: Karger, 1989, vol. 6

【非特許文献4】A. Ossadtchi, S. Baillet, J. Mosher, D. Thyerlei, W. Sutherling, and R. Leahy, “Automated interictal spike detection and source localization in magnetoencephalography using independent components analysis and spatio-temporal clustering,” Clinical Neurophysiology, vol. 115, no. 3, pp. 508-522, 2004.

【非特許文献5】O. Ronneberger P. Fischer and T. Brox "U-net: Convolutional networks for biomedical image segmentation," Proc. Int. Conf. Medical Image Comput. Comput.-Assisted Intervention, pp. 234-241 2015.

【非特許文献6】K. Sekihara, M. Sahani, and S. S. Nagarajan, “Localization bias and spatial resolution of adaptive and non-adaptive spatial filters for MEG source reconstruction,” NeuroImage, vol. 25, no. 4, pp. 1056-1067, 2005.

【非特許文献7】B. V. Veen, W. V. Drongelen, M. Yuchtman, and A. Suzuki, “Localization of brain electrical activity via linearly constrained minimum variance spatial filtering,” IEEE Transactions on Biomedical Engineering, vol. 44, no. 9, pp. 867-880, 1997.

【非特許文献8】高橋勲, 柿木隆介, 白石秀明, 中里信和, 長峰隆, 渡辺裕貴: 臨床脳磁図検査解析指針, 臨床神経生理学, 2005年, 33巻2号, pp69-86

【非特許文献9】Bagic AI, Knowlton RC, Rose DF, Ebersole JS, et al: American Clinical Magnetoencephalography Society Clinical Practice Guideline 1: recording and analysis of spontaneous cerebral activity, J Clin Neurophysiol 2011, 28, 1.

【非特許文献10】Nowak, Rafal, et al. “Toward a Definition of Meg Spike: Parametric Description of Spikes Recorded Simultaneously by Meg and Depth Electrodes.” Seizure, vol. 18, no. 9, pp. 652-655, 2009.

【非特許文献11】Vaudano, Anna Elisabetta, et al. “Causality within the Epileptic Network: An EEG-fMRI Study Validated by Intracranial EEG.” Frontiers in Neurology, vol. 4, 2013.

【非特許文献12】Edakawa, Kohtaroh, et al. “Detection of Epileptic Seizures Using Phase-Amplitude Coupling in Intracranial Electroencephalography.” Scientific Reports, vol. 6, no. 1, 2016.

【図1】

【図2】

【図3】

【図4】

【図5】

【図6】

【図7】

【図8】

【図9】

【図10】

【図11】

【図12】

【図13】