JP5450269B2 - 磁場解析方法、プログラムおよび磁場解析装置 - Google Patents
磁場解析方法、プログラムおよび磁場解析装置 Download PDFInfo
- Publication number
- JP5450269B2 JP5450269B2 JP2010130972A JP2010130972A JP5450269B2 JP 5450269 B2 JP5450269 B2 JP 5450269B2 JP 2010130972 A JP2010130972 A JP 2010130972A JP 2010130972 A JP2010130972 A JP 2010130972A JP 5450269 B2 JP5450269 B2 JP 5450269B2
- Authority
- JP
- Japan
- Prior art keywords
- magnetic field
- time
- current
- change amount
- change
- 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.)
- Expired - Fee Related
Links
Landscapes
- Measuring Magnetic Variables (AREA)
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
Description
この心磁信号を用いて、心臓の2次元、あるいは、3次元の電流分布を再構成する方法が提案されている(例えば非特許文献3,4参照)。心臓の電流分布を再構成することによって、電流強度を等高線で表し、あわせて電流ベクトルを矢印で表すことができる。このような電流分布図を電流アロー図ともいう。このように心臓の電流分布を可視化することによって、心臓の興奮伝導を直感的に理解することができ、臨床上、多くの有効性が報告されている(例えば非特許文献5、6参照)。
その他の解決手段については、明細書中で適宜説明する。
まず、図1〜図8を参照して、本発明に係る第1実施形態について説明する。
第1実施形態では、電流ベクトルの角度(方向)変化の大きい時刻を検出することを目的とする。すなわち第1実施形態では、磁場計測装置において記録された心磁信号から2次元の電流分布を計算し、電流ベクトルの平均角度の時間的な変化量を計算し、変化量の波形を表示することで心臓における興奮部位の変化を視覚的に分かりやすく表示する。さらに、電流ベクトルの平均角度の時間的な変化量の波形の極大点を抽出し、極大点の時刻を興奮伝導時間(具体的には、電流ベクトルの角度に大きな変化が生じた時刻)として自動検出する。
図1は、本実施形態に係る生体磁場計測システムの全体構成例を示す図である。
生体磁場計測システム10は、磁気シールドルーム5の内部と外部とに分かれて配置される。
磁気シールドルーム5の内部には、複数のSQUID(Superconducting Quantum Interference Device)磁束計を内部に配置して極低温に保持しているクライオスタット2と、クライオスタット2を保持しているガントリ3と、被験者(図示せず)が横になるベッド4が配置されている。ベッド4は、短軸方向(y方向)の移動、長軸方向(x方向)での移動、上下方向(z方向)での移動が可能であり、クライオスタット2の計測領域の位置合わせを容易に行うことができる。クライオスタット2、ガントリ3をあわせて磁場計測装置6と称することがある。
なお、本実施形態では、演算装置1にその特徴を有している。そのため、本実施形態では、磁場計測装置6として、例えば、脳磁計測装置、肺磁計測装置、筋磁計測装置などを適宜適用することができる。
図2は、クライオスタットの構成例を示す図である。
クライオスタット2の底部には、図2の上段に示すように、複数のSQUID磁束計21が並んで配置されている。個々のSQUID磁束計21は、ボビンに超伝導線材を巻きつけたコイルと、コイル上部において、コイルと接続されているSQUIDセンサとを有している。
個々のSQUID磁束計21は、その底部がクライオスタット2の底部と平行になるように設置されている。そして、個々のSQUID磁束計21は、z方向の磁場成分Bzを経時的に計測する。なお、磁気の距離変化量を的確に捉えられるように、クライオスタット2内には複数のSQUID磁束計21がx方向およびy方向には等間隔に配列されている。本実施形態では、個々のSQUID磁束計21間の距離が0.025mであって、その計測面22が0.175m×0.175m、SQUID磁束計21の数が8×8のアレー状に配置した64チャンネルとなっているものとするがこのような配置に限定されるものではない。
また、本実施形態では、煩雑になるのを避けるため、図2下段のように、SQUID磁束計21を8×8のマス目(計測面22)で表現することとする。
図3は、測定時におけるSQUID磁束計の配置を示す図である。
計測時では、胸壁33に対してSQUID磁束計の計測面22が並行となるようSQUID磁束計を配置し、例えば、符号31で示す8行4列目に対応するSQUID磁束計が胸部の剣状突起32の真上に位置するように、計測面22の位置合わせを行う。このとき、例えば、符号34で示す1行8列目のSQUID磁束計を座標系の原点Oとする。そして、計測面22が磁場計測領域となる。
なお、SQUID磁束計21(図2)のコイル面(計測面22)の向きを変えることによって、x軸方向の磁場成分Bxや、y軸方向の磁場成分Byを検出することもできる。
図4は、第1実施形態に係る演算装置の構成例を示す図である。
演算装置1は、情報を処理する処理部100と、情報を記憶するHDD(Hard Disk Drive)などの記憶部120と、情報を表示するディスプレイなどの表示部110(表示装置9)とを有している。
処理部100は、電流分布計算部101、平均角度計算部(電流ベクトル角度計算部)102、平均角度変化量計算部(電流ベクトル角度計算部)103、表示処理部104および極大点抽出部(最大変化点検出部)105とを有してなる。
電流分布計算部101は、クライオスタット2から送信された磁場信号を基に、SQUID磁束計の計測面における電流の大きさを等高線で示し、電流の方向をベクトル(電流ベクトル)で示した電流アロー図(詳細は後記)を各時刻毎に算出する機能を有している。
平均角度計算部102は、各時刻における電流ベクトルの角度(方向)の平均値を計算する機能を有している。
平均角度変化量計算部103は、平均角度計算部102が計算した各時刻毎の電流ベクトルの平均角度を基に、各時刻間における電流ベクトルの平均角度の変化量を計算する機能を有する。
表示処理部104は、平均角度変化量計算部103が算出した結果を表示部110に表示する機能を有する。
極大点抽出部105は、電流ベクトルの平均角度の変化量の極大点を抽出する機能を有する。
次に、図4を参照しつつ、図5に沿って、第1実施形態に係る心磁信号解析処理を説明する。
図5は、第1実施形態に係る心磁信号解析処理の手順を示すフローチャートである。
なお、解析対象となっている各時刻毎に、ステップS101〜S105の各処理を行う。また、各処理の詳細は後記することとする。
心磁信号解析処理のプログラムが起動されることによって、処理が開始されると、まず、電流分布計算部101が、記憶部120に格納されている各時刻における心磁信号を取得し、心磁信号から各時刻毎における電流分布を計算する(S101)ことで、後記する電流アロー図を生成する。
次に、平均角度計算部102が、電流分布から各時刻毎における電流ベクトルの平均角度を計算する(S102)。
そして、平均角度変化量計算部103が、ステップS102で計算された電流ベクトルの平均角度を基に、各時刻間における電流ベクトルの平均角度の時間的な変化量を計算する(S103)。
続いて、表示処理部104が、ステップS103で算出した変化量を時間的変化をグラフ(変化量グラフ)の形式で表示する(S104)。
そして、極大点抽出部105が、電流ベクトルの平均角度の変化量から極大点の時刻を抽出し(S105)、処理を終了する。
図6は記憶部に格納されている心磁信号の例を示す図である。
図6は、健常者(典型例)の心室脱分極過程(QRS波)の心磁波形を示すグラフである。図6の心磁波形は、64個のSQUID磁束計で記録された心磁波形を重ね合わせて表示したものである。図6から、健常者の心磁波形は3つの成分(初期成分、中期成分、後期成分)で構成されている。
ステップS101における心磁信号から電流分布を計算する方法として電流アロー図を用いることとする。電流アロー図とは、図7に示すように胸壁に垂直な法線(z)方向の心磁(Bz)から解析的に接線(xおよびy)方向の心磁を算出し、この接線方向の心磁を擬似的な電流ベクトルとして計測平面上に投影し、表示したものである。よって、電流アロー図法は計測点と同数の電流ベクトルを再構成でき、電流ベクトルの大きさを矢印の長さで、電流ベクトルの方向を矢印の向きで表示する。
Iy,i ∝ −∂Bz,i/∂x ・・・ (2)
θi = (180/π)arctan(Iy,i/Ix,i) ・・・ (4)
次に、ステップS102の詳細を説明する。
次に、ステップS102において、平均角度計算部102は、ある時刻tにおける平均電流ベクトルIaveを求め、その平均電流ベクトルの角度を電流ベクトルの平均角度とする。ここで、平均角度計算部102は、平均電流ベクトルIaveのx成分Iavexおよびy成分Iaveyを式(5)および式(6)から計算する。
次に、ステップS103の詳細を説明する。
ステップS103において、平均角度変化量計算部103は、時間的に連続した平均電流ベクトルの角度の変化量を求めるため、ある時刻tを中心とする時間範囲(t―Δt〜t+Δt)内の平均電流ベクトルの角度の変化量Δθavetを式(8)から求める。
ステップS104において、表示処理部104は、式(8)から得られる計算範囲(T1〜T2)の平均電流ベクトルの角度の変化量Δθavetの波形をグラフ(変化量グラフ)として表示部110に表示する。
図8は、第1実施形態に係る変化量グラフの例を示す図である。
図8の上段は、記憶部に格納されている心磁波形のグラフ81であり、図6のグラフと同一のものである。
図8の下段は、変化量グラフ82の例を示す図である。変化量グラフ82は、横軸が時間であり、縦軸が平均電流ベクトルの平均角度の変化量(単位°)である。
そして、ステップS105で、極大点抽出部105が、式(8)の計算範囲(T1〜T2)内の極大点の時刻を検出する。ここで、同じ値を有する連続した時間の平均電流ベクトルの角度の変化量Δθtaの極大点がある場合、時刻の早い方を極大点として抽出する。
なお、心磁信号の立ち上がり時刻から、所定時間後(例えば、90msec後)後のピークは、心臓が安静状態へ移行したことに伴う変化によるピークであるため、極大点の候補から除外する。
次に、第1実施形態の手法を健常例(典型例)の心磁信号に適用した例を説明する。
この例では、健常者1例(男性、31歳)を対象とした。
以下、これまでの心磁信号の解析方法を説明する。心磁波形の初期成分(図6参照)に対応する210ms、220msの電流アロー図は電流強度から強い電流ベクトルは中央領域に位置し、その方向は左へ向かっていることが分かる。次に、心磁波形の初期成分に対応する230msの電流アロー図は、心磁波形の初期成分に対応する210ms、220msの電流アロー図とは異なり、電流強度の強い電流ベクトルは2か所(左下領域と中央やや右下の領域)に認められ、その電流ベクトルの方向は左下と左上を向いている。
以上より、230msと260msの付近で、興奮部位が大きく変化していることが分かり、より詳細に調べると、228msと260msであることが目視で分かる。
以上が、これまでの心磁波形の解析手順である。
前記したように、図8の上段は、健常者の心室脱分極区間(QRS波)の心磁波形を示す図であり、図8の下段は、心磁波形と同じ区間の平均電流ベクトルの角度の変化量の波形を示す図である。
図8の下段は、図5のステップS104で表示部110に表示される変化量グラフ82である。
変化量グラフ82から、平均電流ベクトルの角度の変化量の波形が大きく変化する時刻が2つ(230ms付近と250ms付近)あることが目視で認識できる。この時間帯は、前記した図7の電流アロー図から目視で確認された興奮部位が変化する時刻(228msと260ms)とほぼ同じである。よって、図7のような複数の電流アロー図を、1つ1つ目視で解析するような手間をかけることなく、平均電流ベクトルの角度の変化量を図8の下段のような変化量グラフ82とすることで、興奮部位の変化の時刻(具体的には、電流ベクトルの方向が大きく変化する時刻)をユーザが容易に認識することができる。
なお、本実施形態では、電流ベクトルの時刻毎の平均角度を評価値として用いているが、これに限らず、例えば、各時刻において長さが最も長い電流ベクトルの角度や各時刻において長さがある値以上の電流ベクトルの平均角度を評価値として用いてもよい。
次に、図9〜図11を参照して、本発明に係る第2実施形態について説明する。
第2実施形態では、生体磁場計測システム10において記録された心磁信号から2次元の電流分布を計算し、電流ベクトルの強度の強い座標の時間的な変化量を計算し、変化量の波形を表示することで心臓における興奮部位の変化を視覚的に分かりやすく表示する。さらに、電流ベクトルの強度の強い座標の時間的な変化量の波形の極大点を抽出し、それらの極大点の時刻を興奮伝導時間として自動検出する。
つまり、第2実施形態では電流強度の強い箇所(図7において、色が最も濃い箇所)が大きく移動する時刻を検出することを目的とする。
第2実施形態において、生体磁場計測システム10の全体構成は、演算装置1が、図9に示す演算装置1aと置き換わったこと以外は、図1と同様であるため図示および説明を省略する。
図9は、第2実施形態に係る演算装置の構成例を示す図である。
演算装置1aが、図4の演算装置1と異なる点は、処理部100aにおいて、平均角度計算部102の代わりに最大電流座標計算部106を有し、平均角度変化量計算部103の代わりに最大電流座標変化量計算部107を有している点である。
最大電流座標計算部106は、電流アロー図(図7)の各時刻において、電流強度が最も大きい座標(最大電流座標)を計算する機能を有する。
最大電流座標変化量計算部107は、各時刻間における最大電流座標の変化量を算出する機能を有する。
次に、図9を参照しつつ、図10に沿って、第2実施形態に係る心磁信号解析処理を説明する。
図10は、第2実施形態に係る心磁信号解析処理の手順を示すフローチャートである。
ステップS201は、図5のステップS101と同様の処理であるため、説明を省略する。
ステップS202において、最大電流座標計算部106が、ある時刻tの電流アロー図における最大の電流強度を持つ電流ベクトルの座標(最大電流座標)(xm,t,ym,t)を計算する。
そして、ステップS203において、最大電流座標変化量計算部107が、時間的に連続した最大電流座標の時間的な変化量を算出するため、ある時刻tを中心とする時間範囲(t―Δt〜t+Δt)内の最大電流座標の変化量Δrtを式(9)から求める。
そして、ステップS204において、表示処理部104が、計算範囲(T1〜T2)における最大電流座標の変化量Δrtの波形を変化量グラフとして表示部110に表示する。
さらに、ステップS205において、極大点抽出部105が、式(9)の計算範囲(T1〜T2)内における最大電流座標の変化量Δrtの極大点の時刻を抽出する。ここで、連続した時間における最大電流座標の変化量Δrtが極大点として同じ値を有している場合、時刻の早い方を極大点として用いる。
次に、第2実施形態の手法を第1実施形態と同一の健常例(典型例)の心磁信号に適用した例を説明する。
図11の上段は、健常者(第一の実施形態で用いた被験者と同じ)の心室脱分極区間(QRS波)の心磁波形を示す図である。図11の上段におけるグラフ1101は、図6および図8の上段と同一であるため、説明を省略する。
図11の下段は、図11の上段における心磁波形と同じ区間の最大電流座標の変化量の波形を示す変化量グラフ1102である。変化量グラフ1102から、最大電流座標の変化量の波形が急激に変化する時間帯が2つ(230ms付近と260ms付近)存在することがみてとれる。この時間帯は、第1実施形態で説明した健常者の電流アロー図から目視によって確認された興奮部位が変化する時間帯(228msと260ms)とほぼ同じである。よって、最大電流座標の変化量グラフ1102によって、複数の電流アロー図を比較検討することなく、興奮部位の変化(具体的には、最大電流座標が大きく変化する時刻)を容易に認識することができることがわかる。
次に、図12〜図14を参照して、本発明に係る第3実施形態について説明する。
第3実施形態では、生体磁場計測システム10において記録された心磁信号から2次元の電流分布を計算し、電流ベクトルの平均角度の変化量と電流ベクトルの強度の強い座標の時間的な変化量とを合成した合成変化量を計算し、変化量の波形を表示することで心臓における興奮部位の変化を視覚的に分かりやすく表示する。さらに、合成変化量の波形の極大点を抽出し、それらの極大点の時刻を興奮伝導時間として自動検出する。
つまり、第3実施形態では、第1実施形態の平均電流ベクトルの角度の変化量と、第2実施形態の最大電流座標の変化量とをあわせて考慮することによって、両者が大きく変化した時刻を検出することを目的とする。
なお、第2実施形態において、生体磁場計測システム10の全体構成は、演算装置1が、図12に示す演算装置1bと置き換わったこと以外は、図1と同様であるため図示および説明を省略する。
図12は、第3実施形態に係る演算装置の構成例を示す図である。
演算装置1bが、図4の演算装置1と異なる点は、処理部100bにおいて、図4の演算装置1の構成に加えて、図9の最大電流座標計算部106および最大電流座標変化量計算部107を有し、さらに平均電流ベクトルの角度の変化量と、第2実施形態の最大電流座標の変化量とをあわせた合成変化量を計算する合成変化量計算部108を有する点である。
次に、図12を参照しつつ、図13に沿って、第3実施形態に係る心磁信号解析処理を説明する。
図13は、第3実施形態に係る心磁信号解析処理の手順を示すフローチャートである。
まず、ステップS301〜S303は、図5のステップS101〜S103と同様の処理であるため、説明を省略する。また、ステップS304,S305は、図10のステップS202,S203と同様であるため説明を省略する。
ステップS306において、合成変化量計算部108が、ステップS303で計算された平均電流ベクトルの角度の変化量Δθavet(式(8)の結果)と、ステップS305で計算された最大電流座標の変化量Δrt(式(9)の結果)を合成した合成変化量ΔFtを式(10)から計算する。
なお、合成変化量は、平均電流ベクトルの角度の変化量および最大電流座標の変化量をあわせて考慮したものであればよく、式(10)に限らない。例えば、式(11)でもよいし、平均電流ベクトルの角度の変化量および最大電流座標の変化量の関数などで表現してもよい。
次に、ステップS307において、表示処理部104が、計算範囲(T1〜T2)の合成変化量の変化量ΔFの波形を変化量グラフとして表示部110に表示する。
さらに、ステップS308において、極大点抽出部105が、計算範囲(T1〜T2)内における合成変化量ΔFtの極大点の時刻を抽出する。ここで、連続した時間における合成変化量ΔFが極大点として同じ値を有している場合、時刻の早い方を極大点として用いる。
図14の上段は、健常者(第一の実施形態で用いた被験者と同じ)の心室脱分極区間(QRS波)の心磁波形を示す図である。図14の上段におけるグラフ1401は、図6、図8の上段および図11の上段と同一であるため説明を省略する。
図14の下段は、図14の上段における心磁波形と同じ区間の合成変化量の波形を示す変化量グラフ1402を示す図である。変化量グラフ1402から、最大電流座標の変化量の波形が急激に変化する時間帯が2つ(230ms付近と260ms付近)存在することが分かる。この時間帯は、健常者の電流アロー図から目視によって確認された興奮部位が変化する時間帯(228msと260ms)とほぼ同じである。よって、合成変化量のグラフによって、複数の電流アロー図を比較検討することなく、興奮部位の変化(具体的には、電流ベクトルの方向と、最大電流座標の双方が大きく変化した時刻)を容易に認識することができることがわかる。
次に、図15を参照して、本実施形態の画面例を説明する。
図15は、表示処理部104によって表示部110に表示される画面例を示す図である。
表示画面1500には、磁場波形表示エリア1501と、変化量表示エリア1502と、電流アロー図表示エリア1503と、時間設定バー1504と、波形表示設定エリア1505と、電流アロー図表示設定エリア1506とを有する。
磁場波形表示エリア1501は、演算装置1(図1)が磁場計測装置6(図1)から取得した各SQUID磁束計における電流値波形を示すグラフを表示するエリアである。ここで、表示されるグラフは、図6、図8の上段、図11の上段および図14の上段などで示されるグラフと同一のものである。
そして、磁場波形表示エリア1501においても、極大点を示す時刻が細実線1521,1522および細破線1523,1524で表示されている。ちなみに、磁場波形表示エリア1501では、平均電流ベクトル角度変化量波形から計算される極大点の時刻が細実線1521,1522で示され、最大電流座標変化量波形から計算される極大点の時刻が細破線1523,1524で示されている。
時間設定エリアは、磁場波形表示エリアおよび変化量表示エリアにおける時間幅を設定する時間幅入力窓1551と、時間幅のオフセットを設定する時間幅オフセット入力窓1552とを有する。
磁場波形設定エリアは、磁場波形における振幅の幅(縦軸の長さ)を設定する振幅入力窓1553と、振幅のオフセットを設定する振幅オフセット入力窓1554とを有する。
平均電流ベクトルの角度の変化量設定エリアは、平均電流ベクトルの角度の変化量波形の振幅の幅を設定する振幅入力窓1555と、振幅の幅のオフセットを設定する振幅オフセット入力窓1556とを有する。
最大電流座標の変化量設定エリアは、電流座標の変化量波形の振幅の幅を設定する振幅入力窓1557と、振幅のオフセットを設定する振幅オフセット入力窓1558とを有する。
そして、ユーザが各窓1551〜1558の入力を終え、実行ボタン1559を選択入力すると、各窓1551〜1558の設定値が磁場波形表示エリアおよび変化量表示エリアに反映される。
なお、各表示波形の時間幅のオフセットは、時間設定バー1504を操作することによっても変更されてもよい。また、図15の波形表示設定エリア1505において、合成変化量の波形設定のための入力窓が設けられてもよい。
本実施形態によれば、電流分布から得られる平均電流ベクトルの角度の変化量の波形や、電流強度の強い電流ベクトルの座標の変化量の波形や、平均電流ベクトルの角度の変化量と電流強度の強い電流ベクトルの座標の変化量との合成変化量の波形をグラフとして表示することによって、時間的に連続した複数枚の電流アロー図を表示し、各々の電流アロー図を比較するという手間をかけることなく、心臓の興奮部位の変化を容易に理解することができる。
また、それらの変化量の波形の極大点の時刻を検出することで、心臓の興奮伝導時間に関する情報を自動検出することができる。
2 クライオスタット
3 ガントリ
4 ベッド
5 磁気シールドルーム
6 磁場計測装置
7 駆動回路
8 アンプフィルタユニット
9 表示装置
10 磁場計測システム
21 SQUID磁束計
22 計測面(磁場計測領域)
100,100a,100b 処理部
101 電流分布計算部
102 平均角度計算部(電流ベクトル角度計算部)
103 平均角度変化量計算部(電流ベクトル角度計算部)
104 表示処理部
105 極大点抽出部(最大変化点検出部)
106 最大電流座標計算部
107 最大電流座標変化量計算部
108 合成変化量計算部
110 表示部
120 記憶部
Claims (13)
- 磁場を計測する磁場計測装置から、磁場信号を取得し、前記取得した磁場信号を基に、磁場の解析を行う磁場解析装置による磁場解析方法であって、
前記磁場解析装置は、
前記磁場計測装置から取得された磁場信号を基に、各時刻における電流ベクトルを算出し、
前記電流ベクトルの角度の時間毎の変化量を算出し、
前記変化量が最も大きい時刻を検出する
ことを特徴とする磁場解析方法。 - 前記磁場解析装置は、
前記磁場信号を、時刻毎に複数取得し、
前記複数の磁場信号から、磁場計測領域における各位置毎の電流ベクトルを算出し、
前記算出された各位置毎の電流ベクトルの角度から、電流ベクトルの角度の平均値を算出する
ことを特徴とする請求項1に記載の磁場解析方法。 - 前記磁場解析装置は、
情報を表示する表示部を、さらに有し、
前記変化量の時間変化と、前記検出した時刻とを、前記表示部に表示する
ことを特徴とする請求項1に記載の磁場解析方法。 - 磁場を計測する磁場計測装置から、磁場信号を取得し、前記取得した磁場信号を基に、磁場の解析を行う磁場解析装置による磁場解析方法であって、
前記磁場解析装置は、
前記磁場計測装置から取得された各時刻における複数の磁場信号から、磁場計測領域における各位置毎の電流強度を算出し、
前記算出した電流強度のうち、最も大きい電流強度を有する、前記磁場計測領域における座標を時刻毎に算出し、
前記座標を基に、最も大きい電流強度の座標の時間毎の移動距離を算出し、
前記算出した移動距離のうち、最も大きい移動距離を示す時刻を検出する
ことを特徴とする磁場解析方法。 - 前記磁場解析装置は、
情報を表示する表示部を、さらに有し、
前記変化量の時間変化と、前記検出した時刻とを、前記表示部に表示する
ことを特徴とする請求項4に記載の磁場解析方法。 - 磁場を計測する磁場計測装置から、磁場信号を取得し、前記取得した磁場信号を基に、磁場の解析を行う磁場解析装置による磁場解析方法であって、
前記磁場解析装置は、
前記磁場計測装置から取得された磁場信号を基に、各時刻における電流ベクトルを算出し、
前記電流ベクトルの角度の時間毎の変化量を算出し、
前記算出した電流強度のうち、最も大きい電流強度を有する、前記磁場計測領域における座標を時刻毎に算出し、
前記座標を基に、最も大きい電流強度の座標の時間毎の移動距離を算出し、
前記変化量と、前記移動距離と、に関する合成変化量を算出し、
前記合成変化量が最も大きくなる時刻を検出する
ことを特徴とする磁場解析方法。 - 前記磁場解析装置は、
各時刻における前記変化量と、すべての変化量のうち、最も大きい変化量との比である変化量比を各時刻において算出し、
各時刻における前記移動距離と、すべての移動距離のうち、最も大きい移動距離との比である移動距離比を各時刻において算出し、
前記合成変化量として、前記変化量比と、前記移動距離比と、を加算した値を算出する
ことを特徴とする請求項6に記載の磁場解析方法。 - 前記磁場解析装置は、
前記複数の磁場信号から、磁場計測領域における各位置毎の電流ベクトルを算出し、
前記算出された各位置毎の電流ベクトルの角度から、電流ベクトルの角度の平均値を算出する
ことを特徴とする請求項6に記載の磁場解析方法。 - 前記磁場解析装置は、
情報を表示する表示部を、さらに有し、
前記変化量の時間変化と、前記検出した時刻とを、前記表示部に表示する
ことを特徴とする請求項6に記載の磁場解析方法。 - 請求項1から請求項9のいずれか一項に記載の磁場解析方法をコンピュータに実行させるためのプログラム。
- 磁場を計測する磁場計測装置から、磁場信号を取得し、前記取得した磁場信号を基に、磁場の解析を行う磁場解析装置であって、
前記磁場計測装置から取得された磁場信号を基に、各時刻における電流ベクトルを算出する電流分布計算部と、
前記電流ベクトルの角度の時間毎の変化量を算出する電流ベクトル角度計算部と、
前記変化量が最も大きい時刻を検出する最大変化点検出部と、
を有することを特徴とする磁場解析装置。 - 磁場を計測する磁場計測装置から、磁場信号を取得し、前記取得した磁場信号を基に、磁場の解析を行う磁場解析装置であって、
前記磁場計測装置から取得された各時刻における複数の磁場信号から、磁場計測領域における各位置毎の電流強度を算出する電流分布計算部と、
前記算出した電流強度のうち、最も大きい電流強度を有する、前記磁場計測領域における座標を時刻毎に算出する最大電流座標計算部と、
前記座標を基に、最も大きい電流強度の座標の時間毎の移動距離を算出する最大電流座標変化量計算部と、
前記算出した移動距離のうち、最も大きい移動距離を示す時刻を検出する最大変化点検出部と、
を有することを特徴とする磁場解析装置。 - 磁場を計測する磁場計測装置から、磁場信号を取得し、前記取得した磁場信号を基に、磁場の解析を行う磁場解析装置であって、
前記磁場計測装置から取得された磁場信号を基に、各時刻における電流ベクトルを算出する電流分布計算部と、
前記電流ベクトルの角度の時間毎の変化量を算出電流ベクトル角度計算部と、
前記算出した電流強度のうち、最も大きい電流強度を有する、前記磁場計測領域における座標を時刻毎に算出最大電流座標計算部と、
前記座標を基に、最も大きい電流強度の座標の時間毎の移動距離を算出最大電流座標変化量計算部と、
前記変化量と、前記移動距離と、に関する合成変化量を算出する合成変化量計算部と、
前記合成変化量が最も大きくなる時刻を検出する最大変化点検出部と、
を有することを特徴とする磁場解析装置。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2010130972A JP5450269B2 (ja) | 2010-06-08 | 2010-06-08 | 磁場解析方法、プログラムおよび磁場解析装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2010130972A JP5450269B2 (ja) | 2010-06-08 | 2010-06-08 | 磁場解析方法、プログラムおよび磁場解析装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2011254944A JP2011254944A (ja) | 2011-12-22 |
JP5450269B2 true JP5450269B2 (ja) | 2014-03-26 |
Family
ID=45471825
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2010130972A Expired - Fee Related JP5450269B2 (ja) | 2010-06-08 | 2010-06-08 | 磁場解析方法、プログラムおよび磁場解析装置 |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP5450269B2 (ja) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP6438814B2 (ja) * | 2015-03-23 | 2018-12-19 | 株式会社日立ハイテクノロジーズ | 心磁場解析装置、心磁場解析システム及び心磁場解析方法 |
JP6821368B2 (ja) * | 2016-09-29 | 2021-01-27 | 日置電機株式会社 | 検出センサおよび検査装置 |
JP6962188B2 (ja) * | 2017-12-28 | 2021-11-05 | 株式会社リコー | 情報処理装置、生体信号計測システム、及びプログラム |
-
2010
- 2010-06-08 JP JP2010130972A patent/JP5450269B2/ja not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
JP2011254944A (ja) | 2011-12-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Fenici et al. | Clinical application of magnetocardiography | |
Kwong et al. | Diagnostic value of magnetocardiography in coronary artery disease and cardiac arrhythmias: A review of clinical data | |
Hänninen et al. | Detection of exercise‐induced myocardial ischemia by multichannel magnetocardiography in single vessel coronary artery disease | |
Yamada et al. | Magnetocardiograms in clinical medicine: unique information on cardiac ischemia, arrhythmias, and fetal diagnosis | |
Correa et al. | Novel technique for ST-T interval characterization in patients with acute myocardial ischemia | |
JP4027867B2 (ja) | 生体磁場計測装置 | |
JP4027746B2 (ja) | 生体磁場計測装置 | |
ter Haar et al. | Difference vectors to describe dynamics of the ST segment and the ventricular gradient in acute ischemia | |
Shin et al. | Incremental diagnostic value of combined quantitative and qualitative parameters of magnetocardiography to detect coronary artery disease | |
Li et al. | Diagnostic outcomes of magnetocardiography in patients with coronary artery disease | |
Stinstra | The reliability of the fetal magnetocardiogram | |
JP5450269B2 (ja) | 磁場解析方法、プログラムおよび磁場解析装置 | |
JP2004507305A (ja) | Mcgにおける虚血の識別、定量、および部分的位置判定 | |
CN108697356A (zh) | 矢量心磁图方法及矢量心磁图系统 | |
Chen et al. | Age and sex dependent variations in the normal magnetocardiogram compared with changes associated with ischemia | |
Correa et al. | Acute myocardial ischemia monitoring before and during angioplasty by a novel vectorcardiographic parameter set | |
Gapelyuk et al. | Cardiac magnetic field mapping quantified by Kullback–Leibler entropy detects patients with coronary artery disease | |
JP3427811B2 (ja) | 生体磁場計測装置 | |
JPH10305019A (ja) | 生体磁場計測方法及び生体磁場計測装置 | |
Kandori et al. | Space‐time database for standardization of adult magnetocardiogram‐making standard MCG parameters | |
Van Leeuwen et al. | Quantification of cardiac magnetic field orientation during ventricular de-and repolarization | |
Shiono et al. | Detection of repolarization abnormalities in patients with cardiomyopathy using current vector mapping technique on magnetocardiogram | |
JP4791797B2 (ja) | 生体磁場計測装置 | |
Ogata et al. | Visualization of three-dimensional cardiac electrical excitation using standard heart model and anterior and posterior magnetocardiogram | |
Watanabe et al. | Magnetocardiography in early detection of electromagnetic abnormality in ischemic heart disease |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20130214 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20131118 |
|
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: 20131203 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20131225 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 5450269 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
LAPS | Cancellation because of no payment of annual fees |