JP2017176464A - 磁界計測装置、磁界計測方法、磁界計測プログラム - Google Patents
磁界計測装置、磁界計測方法、磁界計測プログラム Download PDFInfo
- Publication number
- JP2017176464A JP2017176464A JP2016068429A JP2016068429A JP2017176464A JP 2017176464 A JP2017176464 A JP 2017176464A JP 2016068429 A JP2016068429 A JP 2016068429A JP 2016068429 A JP2016068429 A JP 2016068429A JP 2017176464 A JP2017176464 A JP 2017176464A
- Authority
- JP
- Japan
- Prior art keywords
- magnetic field
- current
- elliptical
- elliptic
- numerical parameter
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 19
- 230000000694 effects Effects 0.000 claims abstract description 41
- 238000005259 measurement Methods 0.000 claims description 35
- 238000004364 calculation method Methods 0.000 claims description 16
- 230000010354 integration Effects 0.000 claims description 8
- 238000001228 spectrum Methods 0.000 claims description 8
- 238000004458 analytical method Methods 0.000 claims description 5
- 230000008569 process Effects 0.000 claims description 3
- 238000000691 measurement method Methods 0.000 claims description 2
- 230000005389 magnetism Effects 0.000 claims 4
- 230000002159 abnormal effect Effects 0.000 abstract description 7
- 238000004088 simulation Methods 0.000 description 19
- 238000010586 diagram Methods 0.000 description 10
- 238000012545 processing Methods 0.000 description 7
- 206010003662 Atrial flutter Diseases 0.000 description 5
- 210000004165 myocardium Anatomy 0.000 description 5
- 230000006870 function Effects 0.000 description 4
- 230000008859 change Effects 0.000 description 3
- 230000003595 spectral effect Effects 0.000 description 3
- 210000000779 thoracic wall Anatomy 0.000 description 3
- 210000000591 tricuspid valve Anatomy 0.000 description 3
- 238000012795 verification Methods 0.000 description 3
- 240000004050 Pentaglottis sempervirens Species 0.000 description 2
- 235000004522 Pentaglottis sempervirens Nutrition 0.000 description 2
- 230000000747 cardiac effect Effects 0.000 description 2
- 238000002591 computed tomography Methods 0.000 description 2
- 238000002595 magnetic resonance imaging Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 239000013589 supplement Substances 0.000 description 2
- 238000011282 treatment Methods 0.000 description 2
- 206010003658 Atrial Fibrillation Diseases 0.000 description 1
- 108010038083 amyloid fibril protein AS-SAM Proteins 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 210000000988 bone and bone Anatomy 0.000 description 1
- 238000013153 catheter ablation Methods 0.000 description 1
- 210000000038 chest Anatomy 0.000 description 1
- 239000012141 concentrate Substances 0.000 description 1
- 239000004020 conductor Substances 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000009472 formulation Methods 0.000 description 1
- 230000014509 gene expression Effects 0.000 description 1
- 210000002837 heart atrium Anatomy 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 210000004072 lung Anatomy 0.000 description 1
- 238000005404 magnetometry Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 230000000877 morphologic effect Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 210000000056 organ Anatomy 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 238000005086 pumping Methods 0.000 description 1
- 210000005245 right atrium Anatomy 0.000 description 1
- 210000005241 right ventricle Anatomy 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Images
Landscapes
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
- Measuring Magnetic Variables (AREA)
Abstract
Description
図1は、本発明の実施形態1に係る磁界計測装置(心磁計)1の全体構成を示す概略図である。磁界計測装置1の構成要素は、磁気シールドルーム2の内部と外部とに分かれて配置される。
磁界計測装置1は、磁気センサ11を用いて、被験者の心臓から発生する磁場を計測する。演算装置8は、駆動回路6およびアンプフィルタユニット7を介してその計測結果を取得し、これを磁界時系列データとして取り込む。磁界時系列データは、計測時刻と磁界値のペアとして取得される。
演算装置8は、ステップS101において取得した磁界時系列データに対して周波数解析を実施することにより、各磁気センサ11の計測結果について周波数スペクトルを取得する。本ステップの具体例については、後述する図7〜図8において、検証シミュレーションの一部として改めて説明する。本ステップにおける周波数解析法としては、周知の高速フーリエ変換法やピリオドグラム法などがある。これらによりフーリエスペクトルやパワースペクトル密度を計算することができる。
演算装置8は、ステップS102における周波数解析結果に基づき、旋回する電気的活動の周期を算出する。周波数スペクトルは磁界時系列データ内に含まれる周波数成分を表しているので、支配的な周波数成分を特定することにより、電気的活動の周期を特定することができる。本ステップの具体例については、後述する図7〜図8において、検証シミュレーションの一部として改めて説明する。
演算装置8は、磁界時系列データが表す磁界値を、ステップS103において求めた旋回する電気的活動の1周期分、時間積分する。演算装置8は、この時間積分を磁気センサ11毎に実施する。磁気センサ11が各時刻において計測した磁界値を1周期分加算することにより、その1周期の電気的活動によって瞬時的に発生した磁界を求めることができる。この磁界を用いて、楕円電流の数値パラメータを以下のステップにより特定する。本ステップにおいて1周期分加算する磁界値には、心室の電気的活動を反映するQRS波やT波が含まれていないことが望ましい。
演算装置8は、楕円電流を表す数値パラメータを様々に変化させながら楕円電流が発生させる磁界値を算出することを繰り返すことにより、1周期分の楕円電流が発生させる磁界と、ステップS104において求めた磁界との間の誤差を最小化する。誤差が最小になったとき、ステップS104において求めた磁界を発生させる電気的活動が、その楕円電流によって最もよくモデル化されたと推定することができる。
楕円電流を表す数値パラメータとしては、(a)楕円の中心座標、(b)楕円の方位角、(c)楕円の仰俯角、(d)楕円の径、(e)電流値、(f)電流が流れる方向、が挙げられる。中心座標は、図2におけるXYZ座標によって表される。方位角は、XY平面においてX軸に対して楕円の長軸(または短軸)が形成する角度である。仰俯角は、XY平面に対して楕円平面が形成する角度である。楕円が真円である場合、楕円の径は半径または直径である。楕円が真円でない場合、楕円の径は長径と短径である。電流値と電流方向は、磁界値とその方向から求めることができる。これらのパラメータをノイマンの公式に対して適用することにより、楕円電流が発生させる磁界を計算することができる。
本ステップにおいて両磁界値間の誤差を最小化する手順としては、例えば滑降シンプレックス法や準ニュートン法などによって両磁界値間の2乗誤差を最小化する手法を用いることができる。
表示装置9は、演算装置8が算出した楕円電流の数値パラメータを画面表示する。これらパラメータに加えて、楕円電流の形状を画面表示することもできる。
楕円電流によって旋回する電気的活動をモデル化することの有効性を検証するため、心房粗動を模擬した旋回する電気的活動が発生させる磁場をシミュレーション計算した。このシミュレーションは電流双極子が楕円旋回することにより発生させる磁場を計算することにより実施した。さらに、そのシミュレーションにより得られた磁場を磁場時系列データの代わりに用いて図3のフローチャートを実施することにより、楕円電流が発生させる磁場をシミュレーション計算した。これらシミュレーション結果を比較することにより、楕円電流モデルの妥当性を検証した。以下シミュレーションの詳細を説明する。
心臓内を旋回する電気的活動を楕円電流によってどの程度の精度でモデル化することができたかについての評価結果を提示することにより、磁界計測装置1のユーザは数値パラメータの有効性を判断することができる。そこで本発明の実施形態2では、楕円電流モデルの適合度を計算する方法を説明する。
図11は、表示装置9が表示する画面900の1例を示す図である。画面900は、位置表示欄910とパラメータ表示欄920を有する。位置表示欄910は、磁気センサ11と楕円電流の位置関係を表示する。パラメータ表示欄920は、推定した楕円電流の数値パラメータを表示する。
図12は、本発明の実施形態4に係る磁界計測装置1の構成図である。本実施形態4に係る磁界計測装置1は実施形態1と概ね同様の構成を備えるが、第2演算装置20と第2表示装置21をさらに備える点が異なる。以下では主にこれら差異点について説明する。
本発明は上記した実施例に限定されるものではなく、様々な変形例が含まれる。例えば、上記した実施例は本発明を分かりやすく説明するために詳細に説明したものであり、必ずしも説明した全ての構成を備えるものに限定されるものではない。また、ある実施例の構成の一部を他の実施例の構成に置き換えることが可能であり、また、ある実施例の構成に他の実施例の構成を加えることも可能である。また、各実施例の構成の一部について、他の構成の追加・削除・置換をすることが可能である。
2:磁気シールドルーム
3:クライオスタット
4:ガントリ
5:ベッド
6:駆動回路
7:アンプフィルタユニット
8:演算装置
9:表示装置
10:胸壁
11:磁気センサ
13:剣状突起
20:第2演算装置
21:第2表示装置
Claims (14)
- 被験者の複数の位置における磁気をそれぞれ計測する複数の磁気センサ、
前記複数の磁気センサから計測結果として磁界時系列データを取得するデータ取得部、
前記磁界時系列データから前記被験者の旋回する生体的な電気的活動の旋回周期を算出するとともに1周期分の前記電気的活動と等価な楕円電流を特定する数値パラメータを算出する演算装置、
前記演算装置が算出した前記数値パラメータを出力する出力部、
を備えることを特徴とする磁界計測装置。 - 前記演算装置は、各前記磁気センサから取得した前記磁界時系列データに対して周波数解析を実施することにより前記磁気センサごとに前記磁界時系列データの周波数スペクトルを算出するとともに、各前記磁気センサの前記周波数スペクトルを合算し、
前記演算装置は、前記合算によって得られた周波数成分のうち直流成分を除くもののなかでパワーが最も高いものを、前記旋回周期として特定する
ことを特徴とする請求項1記載の磁界計測装置。 - 前記演算装置は、前記磁界時系列データが記述している磁場を前記旋回周期の1周期分積算するとともに、前記積算の結果と前記楕円電流が発生させる磁場とを比較することにより、前記数値パラメータを算出する
ことを特徴とする請求項2記載の磁界計測装置。 - 前記演算装置は、前記楕円電流を数学的に表す変数を変化させながら前記楕円電流が発生させる磁場を計算することを繰り返すことにより、前記楕円電流が発生させる磁場と前記積算の結果との間の差分を最小化し、
前記演算装置は、前記最小化によって得られた前記変数を前記数値パラメータとして特定する
ことを特徴とする請求項3記載の磁界計測装置。 - 前記数値パラメータは、前記楕円電流の中心座標、前記楕円電流の方位角、前記楕円電流の仰俯角、前記楕円電流の径を少なくとも含む
ことを特徴とする請求項1記載の磁界計測装置。 - 前記楕円電流は、真円経路を流れる電流または楕円経路を流れる電流であり、
前記演算装置は、前記楕円電流が真円経路を流れる電流である場合は、前記楕円電流の径として前記真円の半径または直径を算出し、
前記演算装置は、前記楕円電流が楕円経路を流れる電流である場合は、前記楕円電流の径として前記楕円の短径と長径を算出する
ことを特徴とする請求項1記載の磁界計測装置。 - 前記数値パラメータは、前記楕円電流が流れる方向、前記楕円電流の電流値、を少なくとも含む
ことを特徴とする請求項1記載の磁界計測装置。 - 前記演算装置は、前記楕円電流が発生させる磁場を計算するとともに、その計算結果と前記磁界時系列データが表す磁場との間の合致率を算出し、
前記出力部は、前記演算装置が算出した前記合致率を出力する
ことを特徴とする請求項1記載の磁界計測装置。 - 前記数値パラメータは、前記楕円電流の中心座標を少なくとも含み、
前記出力部は、前記数値パラメータとして前記中心座標を出力する
ことを特徴とする請求項1記載の磁界計測装置。 - 前記数値パラメータは、前記楕円電流の方位角、前記楕円電流の仰俯角、前記楕円電流の径を少なくとも含み、
前記出力部は、前記数値パラメータとして、前記方位角、前記仰俯角、前記径を出力する
ことを特徴とする請求項1記載の磁界計測装置。 - 前記演算装置は、前記数値パラメータとして、前記楕円電流の位置および形状を特定する変数、前記楕円電流が流れる方向を算出し、
前記出力部は、前記演算装置による算出結果に基づき、前記楕円電流の位置および形状と前記楕円電流が流れる方向を出力する
ことを特徴とする請求項1記載の磁界計測装置。 - 前記データ取得部は、各前記磁気センサが磁気を計測する位置と関連付けて前記磁界時系列データを取得し、
前記出力部は、各前記磁気センサの位置と前記楕円電流の位置および形状を同一座標系の下で出力する
ことを特徴とする請求項1記載の磁界計測装置。 - 複数の磁気センサを用いて被験者の複数の位置における磁気をそれぞれ計測した結果を表す磁界時系列データを取得するデータ取得ステップ、
前記磁界時系列データから前記被験者の旋回する生体的な電気的活動の旋回周期を算出するとともに1周期分の前記電気的活動と等価な楕円電流を特定する数値パラメータを算出する演算ステップ、
前記演算ステップにおいて算出した前記数値パラメータを出力する出力ステップ、
を有することを特徴とする磁界計測方法。 - 磁界を計測する処理をコンピュータに実行させる磁界計測プログラムであって、前記コンピュータに、
複数の磁気センサを用いて被験者の複数の位置における磁気をそれぞれ計測した結果を表す磁界時系列データを取得するデータ取得ステップ、
前記磁界時系列データから前記被験者の旋回する生体的な電気的活動の旋回周期を算出するとともに1周期分の前記電気的活動と等価な楕円電流を特定する数値パラメータを算出する演算ステップ、
前記演算ステップにおいて算出した前記数値パラメータを出力する出力ステップ、
を実行させることを特徴とする磁界計測プログラム。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2016068429A JP6532425B2 (ja) | 2016-03-30 | 2016-03-30 | 磁界計測装置、磁界計測方法、磁界計測プログラム |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2016068429A JP6532425B2 (ja) | 2016-03-30 | 2016-03-30 | 磁界計測装置、磁界計測方法、磁界計測プログラム |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2017176464A true JP2017176464A (ja) | 2017-10-05 |
JP6532425B2 JP6532425B2 (ja) | 2019-06-19 |
Family
ID=60003149
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2016068429A Expired - Fee Related JP6532425B2 (ja) | 2016-03-30 | 2016-03-30 | 磁界計測装置、磁界計測方法、磁界計測プログラム |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP6532425B2 (ja) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2020101451A (ja) * | 2018-12-21 | 2020-07-02 | 日置電機株式会社 | 電流分布解析装置、電流分布解析方法、及びプログラム |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH10305019A (ja) * | 1997-03-07 | 1998-11-17 | Hitachi Ltd | 生体磁場計測方法及び生体磁場計測装置 |
JPH11221198A (ja) * | 1998-02-06 | 1999-08-17 | Shimadzu Corp | 生体磁気計測装置 |
JP2002320596A (ja) * | 2001-04-27 | 2002-11-05 | Hitachi Ltd | 生体磁場計測装置 |
JP2003061924A (ja) * | 2001-08-30 | 2003-03-04 | Katsu Yamatani | 生体活動解析装置 |
US20110081058A1 (en) * | 2009-10-06 | 2011-04-07 | Chenyu Wu | Reconstructing Three-Dimensional Current Sources from Magnetic Sensor Data |
US20130324832A1 (en) * | 2012-05-31 | 2013-12-05 | Chenyu Wu | 2D Dipole Localization Using Absolute Value of MCG Measurements |
-
2016
- 2016-03-30 JP JP2016068429A patent/JP6532425B2/ja not_active Expired - Fee Related
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH10305019A (ja) * | 1997-03-07 | 1998-11-17 | Hitachi Ltd | 生体磁場計測方法及び生体磁場計測装置 |
JPH11221198A (ja) * | 1998-02-06 | 1999-08-17 | Shimadzu Corp | 生体磁気計測装置 |
JP2002320596A (ja) * | 2001-04-27 | 2002-11-05 | Hitachi Ltd | 生体磁場計測装置 |
JP2003061924A (ja) * | 2001-08-30 | 2003-03-04 | Katsu Yamatani | 生体活動解析装置 |
US20110081058A1 (en) * | 2009-10-06 | 2011-04-07 | Chenyu Wu | Reconstructing Three-Dimensional Current Sources from Magnetic Sensor Data |
US20130324832A1 (en) * | 2012-05-31 | 2013-12-05 | Chenyu Wu | 2D Dipole Localization Using Absolute Value of MCG Measurements |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2020101451A (ja) * | 2018-12-21 | 2020-07-02 | 日置電機株式会社 | 電流分布解析装置、電流分布解析方法、及びプログラム |
Also Published As
Publication number | Publication date |
---|---|
JP6532425B2 (ja) | 2019-06-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US20220265161A1 (en) | Systems and methods for reconstruction of intra-body electrical readings to anatomical structure | |
Samuelsson et al. | Spatial fidelity of MEG/EEG source estimates: A general evaluation approach | |
Pfeiffer et al. | Localizing on-scalp MEG sensors using an array of magnetic dipole coils | |
RU2622371C2 (ru) | Система обеспечения карты электрической активности | |
Huang et al. | A novel integrated MEG and EEG analysis method for dipolar sources | |
EP3139826B1 (en) | Estimating distribution, fluctuation and/or movement of electrical activity through a heart tissue | |
KR101310747B1 (ko) | 비침습적 심근 전기활동 매핑 방법 | |
JP2012179352A (ja) | 電流双極子を構築するシステムおよび電流双極子を構築する方法 | |
JP2015504704A (ja) | 電磁的トラッキング補正、検出及びエラー訂正のための歪みフィンガープリンティング | |
JP2012157696A (ja) | 心磁図システムおよび心磁図画像を作成する方法 | |
Andersen et al. | Similarities and differences between on-scalp and conventional in-helmet magnetoencephalography recordings | |
Alday et al. | Comparison of electric-and magnetic-cardiograms produced by myocardial ischemia in models of the human ventricle and torso | |
Nurminen et al. | Improving MEG performance with additional tangential sensors | |
JP2008142154A (ja) | 生体磁場計測装置および生体モデルへの平行投影方法 | |
JP6532425B2 (ja) | 磁界計測装置、磁界計測方法、磁界計測プログラム | |
JP3067728B2 (ja) | 生体磁場計測装置 | |
WO2016193185A1 (en) | Determination of positions of objects, such as brachytherapy seeds | |
PT1005677E (pt) | Sistema e método para a tomografia da corrente eléctrica primária do cérebro e do coração | |
Wierenga | MEG, EEG and the integration with magnetic resonance images | |
Bradshaw et al. | Surface current density mapping for identification of gastric slow wave propagation | |
Allouch et al. | Methods used to estimate EEG source-space networks: A comparative simulation-based study | |
Nagarajan et al. | Magnetoencephalographic imaging | |
JP4791797B2 (ja) | 生体磁場計測装置 | |
Yadan et al. | An expert review of the inverse problem in electrocardiographic imaging for the non-invasive identification of atrial fibrillation drivers | |
Primin et al. | Magnetocardiographic technology for human heart investigation |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20180316 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20181218 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20181219 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20190131 |
|
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: 20190508 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20190521 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 6532425 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
S531 | Written request for registration of change of domicile |
Free format text: JAPANESE INTERMEDIATE CODE: R313531 |
|
S533 | Written request for registration of change of name |
Free format text: JAPANESE INTERMEDIATE CODE: R313533 |
|
R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
LAPS | Cancellation because of no payment of annual fees |