特開 2024-135382 再構成方法、再構成プログラム及び生体磁気計測装置

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

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

(11)【公開番号】P2024135382

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

(54)【発明の名称】再構成方法、再構成プログラム及び生体磁気計測装置

(51)【国際特許分類】

   A61B 5/242 20210101AFI20240927BHJP

   A61B 5/243 20210101ALI20240927BHJP

   A61B 5/245 20210101ALI20240927BHJP

   A61B 6/46 20240101ALI20240927BHJP

   A61B 5/055 20060101ALI20240927BHJP

   A61B 8/14 20060101ALI20240927BHJP

【ＦＩ】

A61B5/242

A61B5/243

A61B5/245

A61B6/00 360B

A61B5/055 380

A61B5/055 390

A61B8/14

【審査請求】未請求

【請求項の数】12

【出願形態】ＯＬ

(21)【出願番号】P 2023046034

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

(71)【出願人】

【識別番号】000006747

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

(71)【出願人】

【識別番号】504179255

【氏名又は名称】国立大学法人 東京医科歯科大学

(74)【代理人】

【識別番号】100107766

【弁理士】

【氏名又は名称】伊東 忠重

(74)【代理人】

【識別番号】100070150

【弁理士】

【氏名又は名称】伊東 忠彦

(74)【代理人】

【識別番号】100107515

【弁理士】

【氏名又は名称】廣田 浩一

(72)【発明者】

【氏名】渡部 泰士

(72)【発明者】

【氏名】川端 茂▲徳▼

(72)【発明者】

【氏名】横田 隆徳

(72)【発明者】

【氏名】関原 謙介

(72)【発明者】

【氏名】赤座 実穂

(72)【発明者】

【氏名】大谷 泰

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

4C093

4C096

4C127

4C601

【Ｆターム（参考）】

4C093AA01

4C093AA22

4C093FF35

4C096AA18

4C096AD14

4C096DC33

4C127AA10

4C127BB05

4C127HH13

4C601EE11

4C601LL33

(57)【要約】

【課題】神経電気活動のモデルを想定してLeadfield行列を作成することで、神経機能評価を簡易に行うことを可能にする。
【解決手段】生体電気活動により発生する生体磁気に基づいて生体内の生体電気活動を再構成する再構成方法であって、生体組織において生体電気活動が伝播する伝播方向上の点を挟んで、前記伝播方向上で対向する２つの電流双極子又は前記伝播方向の垂直方向で対向する少なくとも２つの電流双極子を電流源として設定し、設定した２つの電流双極子が発生する電流の広がりを電流保存則に基づいて計算することでLeadfield行列を生成し、生成したLeadfield行列を用いて生体電気活動を再構成する。
【選択図】図１

【特許請求の範囲】

【請求項1】

生体電気活動により発生する生体磁気に基づいて生体内の生体電気活動を再構成する再構成方法であって、
生体組織において生体電気活動が伝播する伝播方向上の点を挟んで、前記伝播方向上で対向する２つの電流双極子又は前記伝播方向の垂直方向で対向する少なくとも２つの電流双極子を電流源として設定し、
設定した２つの電流双極子が発生する電流の広がりを電流保存則に基づいて計算することでLeadfield行列を生成し、
生成したLeadfield行列を用いて生体電気活動を再構成することを特徴とする再構成方法。

【請求項2】

計測部により計測された生体磁気に対応する電流分布を、前記生成したLeadfield行列を用いて生体電気活動に再構成することを特徴とする請求項１に記載の再構成方法。

【請求項3】

前記２つの電流双極子が発生する電流が無限平面を流れるとして前記電流の広がりを計算することで、Leadfield行列を生成することを特徴とする請求項１又は請求項２に記載の再構成方法。

【請求項4】

前記２つの電流双極子の間隔は、前記生体組織の前記伝播方向に対して垂直方向の幅に設定されることを特徴とする請求項１又は請求項２に記載の再構成方法。

【請求項5】

前記伝播方向を、生体の形態情報に基づいて決定することを特徴とする請求項１又は請求項２に記載の再構成方法。

【請求項6】

前記生体の形態情報は、超音波画像、単純Ｘ線画像又はＭＲ画像であることを特徴とする請求項４に記載の再構成方法。

【請求項7】

前記生成したLeadfield行列を用いて再構成した前記生体電気活動を示す画像を、前記超音波画像、前記単純Ｘ線画像又は前記ＭＲ画像に重畳して表示部に表示することを特徴とする請求項６に記載の再構成方法。

【請求項8】

時間の経過とともに変化する前記生体電気活動の時間毎の画像を前記表示部に表示することを特徴とする請求項７に記載の再構成方法。

【請求項9】

前記生成したLeadfield行列を用いて空間フィルター法を用いて生体電気活動を再構成することを特徴とする請求項１又は請求項２に記載の再構成方法。

【請求項10】

生体電気活動により発生する生体磁気に基づいて生体内の生体電気活動を再構成する再構成プログラムであって、
生体組織において生体電気活動が伝播する伝播方向上の点を挟んで、前記伝播方向上で対向する２つの電流双極子又は前記伝播方向の垂直方向で対向する少なくとも２つの電流双極子を電流源として設定する設定処理と、
設定した２つの電流双極子が発生する電流の広がりを電流保存則に基づいて計算することでLeadfield行列を生成する生成処理と、
生成したLeadfield行列を用いて生体電気活動を再構成する再構成処理と、
をコンピュータに実行させることを特徴とする再構成プログラム。

【請求項11】

生体電気活動により発生する生体磁気を計測する計測部と、
生体組織において生体電気活動が伝播する伝播方向上の点を挟んで、前記伝播方向上で対向する２つの電流双極子又は前記伝播方向の垂直方向で対向する少なくとも２つの電流双極子を電流源として設定する設定部と、
設定した２つの電流双極子が発生する電流の広がりを電流保存則に基づいて計算することでLeadfield行列を生成する生成部と、
前記計測部が計測した生体磁気に対応する電流分布を、前記生成部が生成したLeadfield行列を用いて再構成する再構成部と、
前記再構成部により再構成された電流分布に基づいて生体電気活動を再構成する再構成部と、
を有することを特徴とする生体磁気計測装置。

【請求項12】

生体電気活動により発生する生体磁気を計測する計測部と、
生体組織において生体電気活動が伝播する伝播方向上の点を挟んで、前記伝播方向上で対向する２つの電流双極子又は前記伝播方向の垂直方向で対向する、電流源として設定された少なくとも２つの電流双極子が発生する電流の広がりを電流保存則に基づいて計算することで生成されたLeadfield行列が記憶された記憶部と、
前記計測部が計測した生体磁気に対応する電流分布を、前記記憶部に記憶されたLeadfield行列を用いて再構成する再構成部と、
前記再構成部により再構成された電流分布に基づいて生体電気活動を再構成する再構成部と、
を有することを特徴とする生体磁気計測装置。

【発明の詳細な説明】

【技術分野】

【0001】

本発明は、再構成方法、再構成プログラム及び生体磁気計測装置に関する。

【背景技術】

【0002】

被検体の磁気の計測データに基づいて、マトリックス状に配置された複数個所で再構成した電流成分と、被検体の形態画像上で指定された神経経路に沿って抽出した複数個所での電流成分とを表示する生体情報表示装置が知られている（例えば、特許文献１、非特許文献１参照）。

【発明の概要】

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

【0003】

この種の生体情報表示装置を使用して神経機能を評価する場合、評価者は、形態画像に重畳して表示される電流の分布の変化に基づいて、神経経路上を伝播する脱分極点の位置と、脱分極点での電流強度の時間変化とを推測していた。なお、神経活動電流の強度は、脱分極点が障害部位を通過する際に低下することが分かっている。

【0004】

このため、推定された脱分極点の位置が正しくない場合、神経機能を正しく評価することができなかった。換言すれば、電流の分布から脱分極点の位置と脱分極点での電流強度の時間変化とを推定する手法では、神経機能を直感的に評価することが難しく、神経機能を正しく評価するためには習熟を要していた。

【0005】

開示の技術は、上記の課題に鑑みてなされたものであり、神経軸索の電気活動のモデルを想定してLeadfield行列を作成することで、神経機能を簡易に評価可能にすることを目的とする。

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

【0006】

上記技術的課題を解決するため、本発明の一形態の再構成方法は、生体電気活動により発生する生体磁気に基づいて生体内の生体電気活動を再構成する再構成方法であって、生体組織において生体電気活動が伝播する伝播方向上の点を挟んで、前記伝播方向上で対向する２つの電流双極子又は前記伝播方向の垂直方向で対向する少なくとも２つの電流双極子を電流源として設定し、設定した２つの電流双極子が発生する電流の広がりを電流保存則に基づいて計算することでLeadfield行列を生成し、生成したLeadfield行列を用いて生体電気活動を再構成することを特徴とする。

【発明の効果】

【0007】

神経電気活動のモデルを想定してLeadfield行列を作成することで、神経機能を簡易に評価可能にすることができる。

【図面の簡単な説明】

【0008】

【図1】本発明の一実施形態に係る生体磁気計測装置の一例を示すブロック図である。

【図2】神経活動電流のモデルの一例を示す説明図である。

【図3】生体の形態画像に重畳して表示される仮想電極と生体電流の一例を示す説明図である。

【図4】図１の生体磁気計測装置により電流強度及び電流波形を表示する処理の一例を示すフロー図である。

【図5】図２に円で示した２つの電流双極子が作る磁気を定式化する場合の電流双極子の位置関係を示す説明図である。

【図6】Leadfield行列を生成する場合に仮定する電流源（電流双極子）と、電流源から発生する電流分布の例を示す説明図である。

【図7】シミュレーションにより発生させた磁気を仮想の磁気センサアレイで計測して得られる磁気分布の例を示す説明図である。

【図8】信号源を再構成するシミュレーションにおいて、電流四極子を仮定してLeadfield行列を生成した場合に得られる等電流分布の例を示す説明図である。

【図9】電流四極子を仮定して生成したLeadfield行列を使用して、実際に計測した磁気データから電流を再構成した場合の電流強度分布の変化の例を示す説明図である。

【図10】信号源を再構成するシミュレーションにおいて、電流四極子を仮定せずにLeadfield行列を生成した場合に得られる等電流分布の例を示す説明図である。

【図11】電流四極子を仮定せずに生成したLeadfield行列を使用して、実際に計測した磁気データから電流を再構成した場合の電流強度分布の変化の例を示す説明図である。

【図12】図１のデータ処理装置のハードウェア構成の一例を示すブロック図である。

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

【0009】

以下、図面を参照して実施の形態の説明を行う。なお、各図面において、同一構成部分には同一符号を付し、重複した説明を省略する場合がある。

【0010】

図１は、本発明の一実施形態に係る生体磁気計測装置の一例を示すブロック図である。例えば、図１に示す生体磁気計測装置１００は、磁気センサ部１０、信号取得部２０、データ処理装置３０、入力装置８０及び表示装置９０を有する。

【0011】

磁気センサ部１０及び信号取得部２０は、計測部の一例である。データ処理装置３０は、ＰＣ（Personal Computer）やサーバ等のコンピュータである。表示装置９０は、単純Ｘ線画像、ＭＲ画像等の被検体の形態画像と生体信号の電流強度及び電流波形等とを表示可能である。

【0012】

信号取得部２０は、ＦＬＬ（Flux Locked Loop）回路２１、アナログ信号処理部２２、ＡＤ（Analog-to-Digital）変換部２３及びＦＰＧＡ（Field-Programmable Gate Array）２４を有する。例えば、磁気センサ部１０及び信号取得部２０は、磁気を遮蔽するシールドルーム内に設置され、データ処理装置３０、入力装置８０及び表示装置９０は、シールドルーム外に設置される。

【0013】

例えば、磁気センサ部１０は、複数のＳＱＵＩＤ（Superconducting QUantum Interference Device；超伝導量子干渉素子）磁気センサを含んでいる。なお、磁気センサ部１０は、ＳＱＵＩＤ磁気センサに代えて、ＭＲ(Magneto Resistive）センサ又はＯＰＡＭ（Optically Pumped Atomic Magnetometer）センサ等の他の方式の磁気センサを含んでもよい。

【0014】

磁気センサ部１０により計測される生体信号は、生体の骨格筋由来の信号、心筋由来の信号、平滑筋由来の信号、又は神経由来の信号である。骨格筋、心筋、平滑筋及び神経は、生体組織の一例である。例えば、骨格筋は、手掌、前腕、上腕、大腿、下腿等の筋である。例えば、心筋は、心臓を構成する筋肉である。例えば、平滑筋は、胃の収縮や腸の蠕動運動を生み出す筋肉である。例えば、神経由来の信号は、体表から刺激可能な神経から発生する。以下で説明する各実施形態は、ＳＱＵＩＤ磁気センサにより計測された骨格筋由来、心筋由来、平滑筋由来又は神経由来の生体信号から生成される電流波形について適用可能である。

【0015】

データ処理装置３０は、入力制御部４０、表示制御部５０、動作制御部６０及び記憶部７０を有する。例えば、動作制御部６０の機能の一部は、データ処理装置３０に搭載されるＣＰＵ（Central Processing Unit）等のプロセッサが実行する生体電気活動の信号源の再構成プログラムがハードウェアと協働して信号源の再構成方法を実施することで実現される。

【0016】

生体磁気計測装置１００は、脳磁計（ＭＥＧ：Magnetoencephalograph）、心磁計（ＭＣＧ:Magnetocardiograph）又は脊磁計（ＭＳＧ：Magnetospinograph）などとして使用される。なお、生体磁気計測装置１００は、脊髄以外の神経磁気又は筋肉磁気の計測に適用されてもよい。

【0017】

磁気センサ部１０は、被検体が発生する磁気を計測し、計測した磁気を電圧として出力する。磁気センサ部１０は、例えば、ベッドに横たわる被検体の磁気の計測部位に対向して設置される複数のＳＱＵＩＤ磁気センサを有する。ＦＬＬ回路２１は、複数のＳＱＵＩＤ磁気センサにより計測された非線形の磁気－電圧特性をそれぞれ線形化することで、ダイナミックレンジを向上させる。

【0018】

アナログ信号処理部２２は、ＦＬＬ回路２１から出力される線形化されたアナログ信号を増幅し、増幅したアナログ信号のフィルタ処理等を実施する。ＡＤ変換部２３は、フィルタ処理されたアナログ信号（磁気信号）をデジタル値に変換する。ＦＰＧＡ２４は、ＡＤ変換部２３によりデジタル化された磁気データのさらなるフィルタ処理や間引き処理等を実施してデータ処理装置３０に転送する。

【0019】

データ処理装置３０において、入力制御部４０は、位置入力部４１及び波形領域指定部４２を有し、マウスやキーボード等の入力装置８０を介してデータ処理装置３０の操作者から入力される各種情報の入力処理を実施する。データ処理装置３０の操作者は、後述する生体信号の電流波形の評価者でもよい。位置入力部４１は、例えば、液晶ディルプレイ等の表示装置９０に表示された被検体の形態画像上での仮想電極ＶＥ（図３）の位置を受け付ける。

【0020】

波形領域指定部４２は、例えば、入力装置８０から受け付けたウィンドウ設定コマンドに基づいて、生体信号の電流波形を表示するウィンドウの座標等を示す情報を波形表示部５２に指示する。波形表示部５２は、表示装置９０の表示画面において、指示された座標で囲まれる領域に生体信号の電流波形等を表示する。

【0021】

表示制御部５０は、画像表示部５１を有し、表示装置９０に単純Ｘ線画像、電流波形及びノイズレベル等を表示する制御を実施する。画像表示部５１の機能は、図３等で説明される。なお、入力装置８０及び表示装置９０は、データ処理装置３０に含まれてもよい。また、プリンタ等の出力装置が、データ処理装置３０に接続されてもよい。

【0022】

動作制御部６０は、電流双極子設定部６１を含むLeadfield行列生成部６２と、信号源再構成実行部６３とを有する。電流双極子設定部６１を含むLeadfield行列生成部６２と、信号源再構成実行部６３の機能は、図３等で説明される。

【0023】

記憶部７０は、例えば、ＨＤＤ（Hard Disk Drive）等の記憶装置により実現され、生体磁気データ７１、形態画像データ７２及び各種の解析パラメータ７３を記憶する領域を有する。生体磁気データ７１は、磁気センサ部１０で計測されて信号取得部２０で処理された磁気の計測データを含む。形態画像データ７２は、図示しないＸ線撮影装置で撮影された単純Ｘ線画像データ、磁気共鳴断層撮影装置で撮影されたＭＲ（Magnetic Resonance）画像データ及び超音波画像の１つまたは複数を含む。以下では、単純Ｘ線画像データから生成される被検体のＸ線形態画像をＸ線画像と称し、ＭＲ画像データから生成される被検体の断面形態画像をＭＲ画像と称する。

【0024】

解析パラメータ７３は、信号取得部２０に設けられるフィルタ（ハイパスフィルタ、ロウパスフィルタ）の設定値や、各種信号処理を実施する時間範囲、要素数等の生体磁気の計測に必要なパラメータと、生体情報表示方法の実施に必要なパラメータとを含む。また、解析パラメータ７３は、Leadfield行列生成部６２により生成されたLeadfield行列を含む。

【0025】

図２は、神経活動電流のモデルの一例を示す説明図である。図２は、図の上下方向に直線状に走行する神経の活動により電流が発生する様子を示しており、図の下側が末梢側であり、図の上側が中枢側である。神経の活動は、生体電気活動の一例である。例えば、末梢神経に電気刺激を与えることで、神経活動電流が神経軸索を下側から上側に向けて伝播する。

【0026】

このとき、図の上側（順方向）に流れる軸索内電流及び図の下側（逆方向）に流れる軸索内電流と、神経軸索外を流れ、脱分極点に帰ってくる電流成分である体積電流とが発生する。図の上側に向けて流れる軸索内電流をリーディング成分と称され、図の下側に向けて流れる軸索内電流をトレイリング成分と称される。

【0027】

図２の脱分極の位置において、神経の幅方向の両側に実線で示す２つ円は、図４及び図５で説明するシミュレーション時に仮想的に設定される２つの電流双極子（電流源）を示す。なお、２つの電流双極子は、図２では、神経活動が伝播する神経軸索方向上の点を挟んで垂直方向で対向する２つの電流双極子が設定される例が示されるが、神経軸索方向上の点を挟んで神経軸索方向上で対向する２つの電流双極子が設定されてもよい。神経軸索方向は、生体電気活動が伝播する伝播方向の一例である。

【0028】

神経機能を詳細に評価するためには、神経軸索に沿って流れる軸索内電流と、体積電流によって神経軸索に垂直に流入する電流とを抽出して、生体電流波形として表示装置９０の画面等に視覚的に表示することが好ましい。また、神経軸索の一部に異常がある場合、異常個所での生体電流の伝播が、正常な個所での生体電流の伝播より弱くなる場合がある。このため、神経を伝播する生体電流の強度の時間的な変化を表示装置９０の画面等に視覚的に表示することが好ましい。

【0029】

本実施形態において、時間tに置いてセンサアレイで計測される磁場データは式（１）で表現される。

【数1】

以下は簡単のため、単一時間点の場合を考える。

【0030】

空間フィルター法などの信号の線形性を前提とした再構成手法は、信号源の存在する領域をボクセルの集合体として定義し、各ボクセルにおける単位信号源の発生する信号を計算することによりLeadfield行列を生成し、得られたLeadfield行列にもとづいて信号源の分布を算出する。

【0031】

位置rにおけるLeadfield行列は、式（２）のように表現できる。

【数2】

【0032】

ｌ_Ｍ（ｒ）は、位置rに単位信号源を置いたとき、M番目のセンサで得られる出力を示す。このとき、計測された磁場データyは、式（３）により表現される。

【数3】

【0033】

ここで、ｓ（ｒ_１），...，ｓ（ｒ_Ｎ）は、位置ｒ_１，...，ｒ_Ｎにおける電流強度である。換言すれば、計測された磁場データは、各ボクセルに置かれた単位信号源が発生する磁場に各ボクセルの信号源強度をかけ合わせたものを、全て足し合わせたものであるといえる。

【0034】

ここから、再構成手法の一種である空間フィルター法は式（４）により各ボクセルの強度を推定する。

【数4】

ここで、

【数5】

は、推定された各ボクセルの強度を並べたベクトルであり、ｗ（ｒ）は、位置rにおける重み行列である。空間フィルター法は重み行列を種々変更することにより、様々な特性を持つフィルターを構成する。

【0035】

例えば、最小二乗法の考え方に基づいて定義されるMinimum-normフィルターの場合、式（５）となる。

【数6】

【0036】

本実施形態は、上記再構成手続きの中で、特にLeadfield行列の定義に関する。電流双極子設定部６１は、例えば、操作者が、表示装置９０に表示された形態情報を見ながら指定した神経軸索の軸方向に対し平行または垂直に、各ボクセル点を中心として一定の間隔を空けて対向するよう電流双極子を定義する。

【0037】

ここで、神経軸索の軸方向の対しての垂直とは、図２における紙面を貫く方向を含めてもよい。換言すれば、神経の軸索に対して垂直に、３６０°に対向する電流双極子を定義してもよい。

【0038】

Leadfield行列生成部６２は、電流双極子設定部６１が設定した少なくとも２つの電流双極子が発生する電流の広がりを考慮してセンサアレイで得られる磁場を計算することで、Leadfield行列を生成する処理を行う。Leadfield行列生成部６２は、生成したLeadfield行列を解析パラメータ７３として記憶部７０に記憶させる。Leadfield行列生成部６２は、Leadfield行列を生成する生成部の一例である。

【0039】

信号源再構成実行部６３は、解析パラメータとして保存されたLeadfield行列を使用して、磁気センサ部１０で計測された被検体の磁気データから、各ボクセルにおける信号源強度を再構成する。

【0040】

信号源再構成実行部６３は、ボクセル毎の電流の強度に基づいて、画像上で同じ電流強度を有する位置を抽出する。信号源再構成実行部６３は、抽出した位置を結ぶ線を電流強度分布線として記憶部７０に記憶させ、画像表示部５１に指示して、電流強度分布線を形態画像上に重畳表示させる。

【0041】

図４は、図１の生体磁気計測装置１００により電流強度及び電流波形を表示する処理の一例を示すフロー図である。図４に示す例では、被検体の計測対象部位の形態画像を取得するため、被検体の磁気信号の計測とは別に、ステップＳ１０において、Ｘ線撮影装置を使用して被検体のＸ線画像が取得される。また、ステップＳ１１において、磁気共鳴断層撮影装置を使用して被検体のＭＲ画像データが取得され、ステップＳ１２において、計測対象にするスライス（ＭＲ画像）が、ＭＲ画像データから抽出される。なお、ステップＳ１１では、ＭＲ画像に代えて、ＣＴ画像が取得されてもよい。また、ステップＳ１２は、データ処理装置３０により実施されてもよい。

【0042】

なお、Ｘ線撮影装置をＳＱＵＩＤ磁気センサとともにシールドルーム内に設置することで、例えば、被検体をシールドルーム内のベッドに横たわらせた状態でＸ線撮影と生体磁気の計測との両方を実施することができる。これにより、Ｘ線撮影と生体磁気の計測とを別々の場所で実施する場合に比べて、Ｘ線画像上での磁気の計測位置を正確に得ることができる。

【0043】

データ処理装置３０により実施されるステップＳ２０からステップＳ２９までの処理は、データ処理装置３０が実行するデータ処理プログラムに基づいて実施されるデータ処理方法により実現される。また、ステップＳ２１，Ｓ２２，Ｓ２３，Ｓ２４，Ｓ２５に示す処理は、データ処理プログラム中に含まれる生体電気活動を再構成する再構成プログラムに基づいて実施される、生体電気活動を再構成する再構成方法により実現される。

【0044】

ステップＳ２０において、データ処理装置３０を使用して、Ｘ線画像とスライスしたＭＲ画像とを重ね合わせて新たな形態画像を生成する処理が実施される。例えば、画像の重ね合わせは、データ処理装置３０の操作者が、表示装置９０に表示された２つの画像の一方をマウス等の入力装置８０を使用して画面上で移動、回転させることで実施される。

【0045】

この場合、入力制御部４０が操作者の操作を受け付け、受け付けた操作にしたがって画像表示部５１が、表示装置９０に表示されている画像を移動又は回転させる。重ね合わせられた画像は、形態画像データ７２として記憶部７０に記憶される。画像表示部５１は、記憶部７０に記憶された形態画像データ７２が示す形態画像を表示装置９０に表示する。

【0046】

なお、Ｘ線画像とＭＲ画像とを重ね合わせる処理は、画像認識技術を使用して自動的に実施されてもよい。この際、ディープラーニング等の機械学習の手法が利用されてもよい。画像を自動的に重ね合わせる場合、ステップＳ２０の処理は、例えば、動作制御部６０に設けられる図示しない重ね合わせ処理部により実施される。

【0047】

次に、ステップＳ２１において、電流双極子設定部６１は、例えば、操作者が、表示装置９０に表示された形態画像を見ながら設定した神経の軸方向の定義を受け付ける。次に、ステップＳ２２において、電流双極子設定部６１は、受け付けた神経の軸方向の定義で示される形態画像上での神経経路の両側に２つの電流双極子の座標を設定する。

【0048】

次に、ステップＳ２３において、Leadfield行列生成部６２は、電流双極子設定部６１により設定された２つの電流双極子が発生する電流の広がりを電流保存則に基づいて計算することでLeadfield行列を生成する。Leadfield行列生成部６２は、生成したLeadfield行列を解析パラメータ７３として記憶部７０に記憶させる。

【0049】

一方、ステップＳ３０において、生体磁気計測装置１００は、例えば、被検体の末梢神経に与える電気パルス刺激に同期して、ＳＱＵＩＤ磁気センサにより被検体の生体磁気を計測する。次に、ステップＳ３１において、信号取得部２０は、計測した磁気に対応してＳＱＵＩＤ磁気センサが出力する電圧に基づいて、デジタルの磁気データを生成し、生成した磁気データを生体磁気データ７１として記憶部７０に記憶させる。なお、ステップＳ３０、Ｓ３１は、繰り返し実施される。ここまでにより、被検体の計測対象部位の生体電流の強度と電流波形とを取得する準備が完了する。

【0050】

ステップＳ２３の後、ステップＳ２４において、信号源再構成実行部６３は、ステップＳ２３で生成されて記憶部７０に記憶されたLeadfield行列を使用して、計測時間毎の磁気データに基づいて電流成分を再構成することで信号源を推定する。信号源再構成実行部６３は、再構成により得た電流の向き、強度、座標等を示す電流情報を、形態画像データ７２として記憶部７０に記憶させる。

【0051】

次に、ステップＳ２５において、信号源再構成実行部６３は、画像表示部５１に指示して、再構成により得たボクセル毎の電流の向きと強度とを示す矢印を、ステップＳ２０で重ね合わせた画像上に表示させる。また、信号源再構成実行部６３は、ボクセル毎の電流の強度に基づいて、画像上で同じ電流強度を有する位置を抽出し、電流強度分布線として画像に重ね合わせて表示させる。信号源再構成実行部６３は、電流強度分布線を示す情報を形態画像データ７２として記憶部７０に記憶させる。これにより、例えば、図３に示したように、形態画像上と、計測時間毎に再構成した電流の分布（電流の矢印と電流強度分布線）とが、画面上に重ね合わせて表示される。

【0052】

図５は、図２に円で示した２つの電流双極子（すなわち、４極子）が作る磁気を定式化する場合の電流双極子の位置関係を示す説明図である。図５に示すように、中心がｒ'で、距離ｄ離れ、互いに逆方向の２つの電流双極子Ｑが作る磁気分布を考える。まず、電流源（primary current）である４極子の電流が作る磁気分布のみを考え、その後、体積電流が作る磁気分布も含めて考える。なお、２つの電流双極子が発生する電流強度は、偏りがないものとする。

【0053】

以下の式において、記号Ｑは、電流双極子の強度ベクトルを示し、記号ｒ'は、２つの電流双極子の中間位置を示す。記号ｄは、２つの電流双極子間の距離ベクトルを示し、記号ｒは、センサの位置を示す。記号μ_０は、透磁率を示し、記号Ｂ（ｒ）は、位置ｒで得られる磁気を示す。

【0054】

記号Ｅは、電流双極子の方向の単位ベクトルを示す。記号Ｕ（ｒ）は、スカラポテンシャルを示し、記号ｔは、Ｚ方向無限大から位置ｒまでの経路（積分経路）を表すパラメータを示す。記号ｅ_ｚは、図５のＸＹ平面に対して垂直方向の単位ベクトルを示す。

【0055】

互いに対向する２つの電流双極子Ｑ（４極子）により位置ｒに作られる磁気Ｂ_０（ｒ）は、式（６）で示される。

【数7】

【0056】

式（６）において、ａ＝ｒ－ｒ'とおき、｜ｄ｜＜＜｜ａ｜としてｄ／／Ｑに注意して計算すると、式（７）が得られる。

【数8】

【0057】

また、Ｅを式（７）中のＱおよびｄの方向の単位ベクトル（Ｅ＝Ｑ／｜Ｑ｜＝ｄ／｜ｄ｜）とすると、式（８）が得られる。

【数9】

【0058】

次に、均一な伝導率σの半無限平面Ｇ（図５のＸＹ平面に平行な面）内に上記の逆向きで等価の強度の電流源がある場合に、その電流源に伴う体積電流による磁気も含めて考える。平面に対して垂直成分の電流源とそれに伴う体積電流は周囲に磁気を作らないことが知られており（Sarvas 1987）、Ｅのｚ成分は、０（Ｅ・ｅ_Ｚ＝０）とする。

【0059】

この磁気の無限半平面に対して垂直成分Ｂ_Ｚ（ｒ）は、Geselowitzの式を用いると、式（９）となり、無限半平面に対して垂直成分は電流源のみによる磁場の垂直成分と一致する。

【数10】

【0060】

なお、式（９）において、Ｖ（ｒ'）は、半無限平面表面の電位を表し、σは、半無限平面の伝導率を表し、ｄＳは、半無限平面表面の微小面積を表す。

【数11】

は、半無限平面表面での積分を表す。

【0061】

被検体の体外には電流は存在しないため、∇×Ｂ＝０が成り立つことから、磁気Ｂ（ｒ）は、スカラポテンシャルＵを用いて式（１０）で表され、その勾配を調べることで平面に対して平行な磁気の成分Ｂｘ，Ｂｙが分かる。

【数12】

【0062】

式（１０）中のスカラポテンシャルＵ（ｒ）を求めるために、式（１１）に示される無限遠から無限平面に垂直方向に沿った積分を考える。

【数13】

【0063】

式（１１）に式（１２）を代入して、式（１３）が得られる。

【数14】

【数15】

【0064】

ｅ_Ｚ・Ｅ＝０、ｅ_Ｚ×Ｅ・ｅ_Ｚ＝０なので、ａ＝ｒ－ｒ_０として、式（１４）が求まる。

【数16】

【0065】

式（１４）の積分部分を式（１５）に示すようにＭとおく。ここで、ａ＝｜ａ｜、ｂ＝ａ・ｅ_Ｚとした。

【数17】

【0066】

式（１５）の右項は、式（１６）となる。

【数18】

【0067】

そして、式（１４）の積部分をＭとおいて式（１７）の勾配を求めることで、式（１８）に示す目的の磁気を得ることができる。すなわち、互いに対向する電流双極子を電流源とする場合に電流源（primary current）と体積電流とにより発生する磁気を定式化することができる。

【数19】

【数20】

【0068】

上記の計算に、▽ａ・Ｅ＝Ｅ、▽［（ａ×Ｅ）・ｅ_Ｚ］＝Ｅ×ｅ_Ｚを用いた。但し、ａ＝｜ａ｜、ｂ＝ａ・ｅ_Ｚとして、Ｍは、式（１９）であり、▽Ｍは、式（２０）である。

【数21】

【数22】

【0069】

以上の一連の計算により、互いに対向する２つの電流双極子（４極子）が発生する位置ｒでの磁気を計算することが可能になる。互いに対向する２つの電流双極子を電流源としてLeadfield行列を生成するアイデアは、従来提案されていない。例えば、式（６）、式（８）、式（１８）は、本発明者らにより導かれた固有の式である。

【0070】

信号源を表すベクトルＳ（＝｜ｄ｜｜Ｑ｜Ｅ）と磁気分布Ｂの関係について、ａ・Ｅ及びａ×Ｅ（なお、Ｅ＝Ｓ／｜Ｓ｜）の項があるために、Ｓの各成分とＢの間に線形関係がない。このため、例えば、Ｓ＝（１，０，０）における磁気分布とＳ＝（０，１，０）における磁気分布を加算してもＳ＝（１，１，０）における磁気分布とはならない。

【0071】

そこで、形態学的情報などによって単位ベクトルＥの方向を指定することで線形性を持たせる。これにより、例えば、空間フィルター法のような信号源の線形性を前提としたアルゴリズムへの応用が可能となる。

【0072】

以上より、互いに対向する２つの電流双極子を電流源として発せられる磁気を計算する上記の式の手続きを用いてLeadfield行列を生成し、生成されたLeadfield行列に基づいて空間フィルター法を用いて計測データから体内の生体活動の強度を再構成することができる。すなわち、時間とともに移動する神経軸索上の電流源に着目して生体活動の強度を再構成できるため、より直接的に生体活動を評価することが可能となる。

【0073】

これに対して、従来手法では、マトリックス状に配置したボクセルでの電流分布を被検体の体内での電流分布として再構成していたため、時間とともに移動する神経軸索上の電流源での電流強度は、位置が固定された複数個所での電流強度の変化に基づいて評価者が推測するしかなかった。

【0074】

図６は、Leadfield行列を生成する場合に仮定する電流源（電流双四極子）と、電流源から発生する電流分布の例を示す説明図である。図６において、神経軸索内の生体電流の伝播経路は、Ｙ方向で示され、ＸＹ平面での電流分布が算出されるとする。２つの電流双極子Ｑは、神経軸索内の生体電流の伝播経路の幅方向（Ｘ方向）の両側に設定され、２つの電流双極子Ｑの間隔は、伝播経路の幅とほぼ等しく設定される。例えば、２つの電流双極子Ｑの間隔は、１ｃｍ程度以下であることが好ましい。

【0075】

図２と同様に、電流源である電流双極子ＱからＹ＋方向（順方向）とＹ－方向（逆方向）とに軸索内電流がそれぞれ流れる。また、電流双極子Ｑから発生し、神経軸索外を流れる二次電流は、体積電流（内向き電流）として電流双極子Ｑに帰ってくる。太い矢印及び細い矢印で示す電流で示される電流分布は、信号源を示す。

【0076】

従来手法では、図１の動作制御部６０が実施する生体電流の再構成により、二次電流を含むＸＹ平面上の全ての電流成分が抽出されており、神経機能の評価に有用でない電流成分も表示装置９０の画面に表示されていた。このため、神経機能の評価者は、画面に表示される電流成分から神経内を伝播する電流の強度及び時間変化を推測することで、神経機能を評価していた。したがって、従来手法では、神経機能を直感的に評価することが難しく、神経機能を評価するためには習熟を要していた。

【0077】

本発明者等は、以下の条件を設定して、電流四極子を仮定してLeadfield行列を生成した場合と、電流四極子を仮定せずにLeadfield行列を生成した場合（従来手法）とのシミュレーションをそれぞれ行った。ここで、電流双極子を仮定しないとは、体積電流を考慮しないことを示し、電流双極子を仮定するとは、体積電流を考慮することを示す。

【0078】

（条件１）Ｙ方向に走行する十分に細い神経のＸ方向の５ｍｍの領域が脱分極していることを想定して、神経の軸からＸ＋方向とＸ－方向とにそれぞれ２．５ｍｍずつ離れた位置に、互いに逆向きで同じ強度（例えば、１Ａｍ）の電流双極子（電流源）をそれぞれ設定した。
（条件２）電流源と磁気センサアレイとのＺ方向の距離は３０ｍｍとした。
（条件３）Sarvas式（無限平面モデル）を用いて電流双極子による磁気と、電流双極子に基づいて発生する体積電流とによる磁気とを発生させ、各磁気センサにより計測される値を計算した。これにより、仮想的に神経軸索活動が発生する磁気データを得ることができる。
（条件４）電流源と同じ平面上に関心領域（ＲＯＩ：Region of Interest）をボクセル間隔＝２．５ｍｍで設定し、（条件３）で得た磁気データに空間フィルター法（UGRENS法）を適用し信号源の推定を行った。

【0079】

図７は、シミュレーションにより発生させた磁気を仮想の磁気センサアレイで計測して得られる磁気分布の例を示す説明図である。図７のＸ＝０、Ｙ＝０は、図５において２つ電流双極子Ｑが設定される位置での神経軸索の中心を示している。

【0080】

分布中に示す矢印は、ＸＹ方向の磁気の向きを示している。色が比較的薄い領域は、磁気のＺ方向への吹き出しを示し、色が比較的濃い領域は、磁気のＺ方向からの吸い込みを示している。

【0081】

図８は、信号源を再構成するシミュレーションにおいて、電流四極子を仮定してLeadfield行列を生成した場合に得られる等電流分布の例を示す説明図である。電流分布において、等電流を示す線の間隔が小さいほど、電流の変化が大きいことを示し、等電流を示す多重の円のうち、内側の円ほど電流強度が高いことを示している。

【0082】

電流四極子を仮定してLeadfield行列を生成した場合の等電流分布では、神経活動のうち、体積電流が神経軸索に流れ込む、着目したい主要部の電流の電流強度をはっきりと認識することができる。したがって、時間とともに変化する電流強度と、時間とともに変化する電流強度が強い位置とを観察することで、従来に比べて神経機能を簡易に評価できることが分かる。

【0083】

図９は、電流四極子を仮定して生成したLeadfield行列を使用して、実際に計測した磁気データから電流を再構成した場合の電流強度分布の変化の例を示す説明図である。図９は、表示装置９０の形態画像用の表示ウィンドウに表示される計測時刻毎に変化する電流強度分布の画像を切り出して並べている。実際の形態画像用の表示ウィンドウに表示される電流強度分布の画像は、計測時刻をマウス等で変更することで、表示が切り替えられる。

【0084】

例えば、図９では、被検体の足首で腓腹神経を電気刺激し、電気刺激により誘発された神経磁気を下腿部で計測し、得られた磁気データを使用して図１の信号源再構成実行部６３により電流成分を再構成した。

【0085】

関心領域ＲＯＩは、神経が存在していると想定される曲面を設定した。なお、図９の形態画像内に表示されている５つの白い円と、円につながる白い線は、体表面の電位を計測するために体表面に貼付された皿電極と信号線である。

【0086】

図９では、１つの強い電流強度が時間とともに神経軸索を移動することを認識できる。したがって、着目したい神経活動がどこにあり、どのように伝播しているのかを画面から直感的に認識することができ、習熟していなくても神経機能を簡易に評価することができる。

【0087】

図１０は、信号源を再構成するシミュレーションにおいて、電流四極子を仮定せずにLeadfield行列を生成した場合に得られる等電流分布の例を示す説明図である。図８と同様に、電流分布において、等電流を示す線の間隔が小さいほど、電流の変化が大きいことを示し、等電流を示す多重の円のうち、内側の円ほど電流強度が高いことを示している。

【0088】

電流四極子を仮定せずにLeadfield行列を生成した場合の等電流分布では、神経軸索のＹ＋方向とＹ－方向とにそれぞれ流れる電流と、Ｘ＋方向とＸ－方向とから神経軸索に流れ込む体積電流とが発生していることを認識できる。しかしながら、着目したい主要部の電流（体積電流が神経軸索に流れ込む電流）の電流強度は、はっきりとは認識することができない。

【0089】

したがって、神経機能の評価者は、画面に表示される電流強度分布から神経軸索内を伝播する電流の強度及び時間変化を推定することで、神経機能を評価しなくてはならない。この結果、神経機能を直感的に評価できず、神経機能を評価するためには習熟を要する。

【0090】

図１１は、電流四極子を仮定せずに生成したLeadfield行列を使用して、実際に計測した磁気データから電流を再構成した場合の電流強度分布の変化の例を示す説明図である。図１１においても、図９と同様に、表示装置９０の形態画像用の表示ウィンドウに表示される計測時刻毎に変化する電流強度分布の画像を切り出して並べている。

【0091】

図１１では、図１０と同様に、神経軸索のＹ＋方向とＹ－方向とにそれぞれ流れる電流と、Ｘ＋方向とＸ－方向とから神経軸索に流れ込む体積電流とを示す電流強度分布の時間変化を認識することができる。しかしながら、着目したい神経活動がどこにあるのかを直感的に判別することが困難であり、神経機能を評価するためには習熟を要する。

【0092】

図１２は、図１のデータ処理装置３０のハードウェア構成の一例を示すブロック図である。データ処理装置３０は、ＣＰＵ３０１とＲＯＭ（Read Only Memory）３０２とＲＡＭ（Random Access Memory）３０３と外部記憶装置３０４とを有する。また、データ処理装置３０は、入力インタフェース部３０５と出力インタフェース部３０６と入出力インタフェース部３０７と通信インタフェース部３０８とを有する。例えば、ＣＰＵ３０１とＲＯＭ３０２とＲＡＭ３０３と外部記憶装置３０４と入力インタフェース部３０５と出力インタフェース部３０６と入出力インタフェース部３０７と通信インタフェース部３０８とは、バスＢＵＳを介して相互に接続される。

【0093】

ＣＰＵ３０１は、ＯＳ及びアプリケーション等の各種プログラムを実行し、データ処理装置３０の全体の動作を制御する。ＲＯＭ３０２は、各種プログラムをＣＰＵ３０１により実行可能にするための基本プログラムや各種パラメータ等を保持する。ＲＡＭ３０３は、ＣＰＵ３０１により実行される各種プログラムや、プログラムで使用するデータを記憶する。各種プログラムは、２つの電流双極子を設定する設定処理と、Leadfield行列を生成する生成処理と、Leadfield行列を用いて生体電気活動を再構成する推定処理とを含む推定プログラムと、図３に示した形態画像を表示する形態画像表示プログラムとを含む。

【0094】

外部記憶装置３０４は、ＨＤＤ（Hard Disk Drive）又はＳＳＤ（Solid State Drive）等であり、ＲＡＭ３０３に展開する各種プログラムを記憶する。各種プログラムには、上述した推定プログラムと形態画像表示プログラムとが含まれてもよい。

【0095】

入力インタフェース部３０５には、データ処理装置３０を操作する操作者等からの入力を受け付けるキーボード、マウスやタブレット等の入力装置８０が接続される。出力インタフェース部３０６には、ＣＰＵ３０１が実行する各種プログラムにより生成される表示画面等を表示する表示装置（例えば、図１の表示装置９０）やプリンタ等の出力装置９２が接続される。

【0096】

入出力インタフェース部３０７には、ＵＳＢ（Universal Serial Bus）メモリ等の記録媒体４００が接続される。例えば、記録媒体４００には、上述した形態画像表示プログラム及び電流波形表示プログラム等の各種プログラムが格納されてもよい。この場合、プログラムは、入出力インタフェース部３０７を介して記録媒体４００からＲＡＭ３０３に転送される。なお、記録媒体４００は、ＣＤ－ＲＯＭやＤＶＤ（Digital Versatile Disc：登録商標）等でもよく、この場合、入出力インタフェース部３０７は、接続する記録媒体４００に対応するインタフェースを有する。通信インタフェース部３０８は、データ処理装置３０をネットワーク等に接続する。

【0097】

以上、この実施形態では、形態画像において生体軸索上に２つの電流双極子を電流源として設定してLeadfield行列を生成し、電流成分を再構成することで、神経軸索に流れ込む体積電流（内向き電流）を形態画像に重畳して表示することができる。これにより、評価者は、時間とともに神経軸索を移動する生体電流と、生体電流の強度の変化とを認識することができる。

【0098】

例えば、図９に示すように、Leadfield行列を用いて再構成した電流強度分布を示す画像は、形態画像に重畳されて表示装置９０の表示部に表示される。また、時間の経過とともに変化する生体電流の強度分布の時間毎の画像は、形態画像に重畳して表示部に表示可能である。

【0099】

これにより、評価者は、着目したい神経活動がどこにあり、どのように伝播しているのかを画面から直感的に認識することができ、習熟していなくても神経機能を簡易に評価することができる。すなわち、神経電気活動のモデルを想定してLeadfield行列を作成することで、神経機能を簡易に評価可能にすることができる。

【0100】

また、電流双極子設定部６１は、神経経路が写るＸ線画像、ＭＲ画像等の形態画像上で指定された点に基づいて、点を挟んで神経経路上で対向する２つの電流双極子、又は、点を挟んで神経経路の垂直方向で対向する２つの電流双極子を設定する。これにより、形態画像を使用せずに点を設定する場合に比べて、神経軸索に流れ込む体積電流の電流強度の再構成の精度を高くすることができる。

【0101】

本発明の態様は、例えば、以下の通りである。
＜１＞
生体電気活動により発生する生体磁気に基づいて生体内の生体電気活動を再構成する再構成方法であって、
生体組織において生体電気活動が伝播する伝播方向上の点を挟んで、前記伝播方向上で対向する２つの電流双極子又は前記伝播方向の垂直方向で対向する少なくとも２つの電流双極子を電流源として設定し、
設定した２つの電流双極子が発生する電流の広がりを電流保存則に基づいて計算することでLeadfield行列を生成し、
生成したLeadfield行列を用いて生体電気活動を再構成することを特徴とする再構成方法。
＜２＞
計測部により計測された生体磁気に対応する電流分布を、前記生成したLeadfield行列を用いて生体電気活動に再構成することを特徴とする＜１＞に記載の再構成方法。
＜３＞
前記２つの電流双極子が発生する電流が無限平面を流れるとして前記電流の広がりを計算することで、Leadfield行列を生成することを特徴とする＜１＞又は＜２＞に記載の再構成方法。
＜４＞
前記２つの電流双極子の間隔は、前記生体組織の前記伝播方向に対して垂直方向の幅に設定されることを特徴とする＜１＞ないし＜３＞のいずれか１項に記載の再構成方法。
＜５＞
前記伝播方向を、生体の形態情報に基づいて決定することを特徴とする＜１＞ないし＜４＞のいずれか１項に記載の再構成方法。
＜６＞
前記生体の形態情報は、超音波画像、単純Ｘ線画像又はＭＲ画像であることを特徴とする＜４＞又は＜５＞に記載の再構成方法。
＜７＞
前記生成したLeadfield行列を用いて再構成した前記生体電気活動を示す画像を、前記超音波画像、前記単純Ｘ線画像又は前記ＭＲ画像に重畳して表示部に表示することを特徴とする＜６＞に記載の再構成方法。
＜８＞
時間の経過とともに変化する前記生体電気活動の時間毎の画像を前記表示部に表示することを特徴とする＜７＞に記載の再構成方法。
＜９＞
前記生成したLeadfield行列を用いて空間フィルター法を用いて生体電気活動を再構成することを特徴とする＜１＞ないし＜８＞のいずれか１項に記載の再構成方法。
＜１０＞
生体電気活動により発生する生体磁気に基づいて生体内の生体電気活動を再構成する再構成プログラムであって、
生体組織において生体電気活動が伝播する伝播方向上の点を挟んで、前記伝播方向上で対向する２つの電流双極子又は前記伝播方向の垂直方向で対向する少なくとも２つの電流双極子を電流源として設定する設定処理と、
設定した２つの電流双極子が発生する電流の広がりを電流保存則に基づいて計算することでLeadfield行列を生成する生成処理と、
生成したLeadfield行列を用いて生体電気活動を再構成する再構成処理と、
をコンピュータに実行させることを特徴とする再構成プログラム。
＜１１＞
生体電気活動により発生する生体磁気を計測する計測部と、
生体組織において生体電気活動が伝播する伝播方向上の点を挟んで、前記伝播方向上で対向する２つの電流双極子又は前記伝播方向の垂直方向で対向する少なくとも２つの電流双極子を電流源として設定する設定部と、
設定した２つの電流双極子が発生する電流の広がりを電流保存則に基づいて計算することでLeadfield行列を生成する生成部と、
前記計測部が計測した生体磁気に対応する電流分布を、前記生成部が生成したLeadfield行列を用いて再構成する再構成部と、
前記再構成部により再構成された電流分布に基づいて生体電気活動を再構成する再構成部と、
を有することを特徴とする生体磁気計測装置。
＜１２＞
生体電気活動により発生する生体磁気を計測する計測部と、
生体組織において生体電気活動が伝播する伝播方向上の点を挟んで、前記伝播方向上で対向する２つの電流双極子又は前記伝播方向の垂直方向で対向する、電流源として設定された少なくとも２つの電流双極子が発生する電流の広がりを電流保存則に基づいて計算することで生成されたLeadfield行列が記憶された記憶部と、
前記計測部が計測した生体磁気に対応する電流分布を、前記記憶部に記憶されたLeadfield行列を用いて再構成する再構成部と、
前記再構成部により再構成された電流分布に基づいて生体電気活動を再構成する再構成部と、
を有することを特徴とする生体磁気計測装置。

【0102】

以上、各実施形態に基づき本発明の説明を行ってきたが、上記実施形態に示した要件に本発明が限定されるものではない。これらの点に関しては、本発明の主旨をそこなわない範囲で変更することができ、その応用形態に応じて適切に定めることができる。

【符号の説明】

【0103】

１０ 磁気センサ部
２０ 信号取得部
２１ ＦＬＬ回路
２２ アナログ信号処理部
２３ ＡＤ変換部
２４ ＦＰＧＡ
３０ データ処理装置
４０ 入力制御部
４１ 位置入力部
４２ 波形領域指定部
５０ 表示制御部
５１ 画像表示部
６０ 動作制御部
６１ 電流双極子設定部
６２ Leadfield行列生成部
６３ 信号源再構成実行部
７０ 記憶部
７１ 生体磁気データ
７２ 形態画像データ
７３ 解析パラメータ
８０ 入力装置
９０ 表示装置
１００ 生体磁気計測装置
３０１ ＣＰＵ
３０２ ＲＯＭ
３０３ ＲＡＭ
３０４ 外部記憶装置
３０５ 入力インタフェース部
３０６ 出力インタフェース部
３０７ 入出力インタフェース部
３０８ 通信インタフェース部
４００ 記録媒体
ＶＥｃ，ＶＥｌ，ＶＥｒ 仮想電極

【先行技術文献】

【特許文献】

【0104】

【特許文献1】特開２０２１－０８３４７９号公報

【非特許文献】

【0105】

【非特許文献1】Andrei Irimia et al., Novel functional imaging technique for the brachial plexus based on magnetoneurography", Clinical Neurophysiology, Volume 130, Issue 11, November 2019, Pages 2114-2123

【図1】

【図2】

【図3】

【図4】

【図5】

【図6】

【図7】

【図8】

【図9】

【図10】

【図11】

【図12】