JP2022169141A - 信号ベクトル導出装置、方法、プログラム、記録媒体 - Google Patents
信号ベクトル導出装置、方法、プログラム、記録媒体 Download PDFInfo
- Publication number
- JP2022169141A JP2022169141A JP2021074991A JP2021074991A JP2022169141A JP 2022169141 A JP2022169141 A JP 2022169141A JP 2021074991 A JP2021074991 A JP 2021074991A JP 2021074991 A JP2021074991 A JP 2021074991A JP 2022169141 A JP2022169141 A JP 2022169141A
- Authority
- JP
- Japan
- Prior art keywords
- signal
- vector
- spectrum
- coefficient
- derivation
- 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.)
- Pending
Links
- 239000013598 vector Substances 0.000 title claims abstract description 165
- 238000009795 derivation Methods 0.000 title claims abstract description 132
- 238000000034 method Methods 0.000 title claims description 26
- 238000001228 spectrum Methods 0.000 claims abstract description 78
- 238000005259 measurement Methods 0.000 claims abstract description 53
- 238000012546 transfer Methods 0.000 claims description 26
- 230000005484 gravity Effects 0.000 claims description 10
- 230000006870 function Effects 0.000 description 28
- 239000011159 matrix material Substances 0.000 description 16
- 230000003595 spectral effect Effects 0.000 description 4
- 230000004907 flux Effects 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 241000238366 Cephalopoda Species 0.000 description 1
- 230000004807 localization Effects 0.000 description 1
- 238000010801 machine learning Methods 0.000 description 1
- 238000013519 translation Methods 0.000 description 1
- 239000002023 wood Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/02—Measuring direction or magnitude of magnetic fields or magnetic flux
- G01R33/0206—Three-component magnetometers
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
- A61B5/242—Detecting biomagnetic fields, e.g. magnetic fields produced by bioelectric currents
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/02—Measuring direction or magnitude of magnetic fields or magnetic flux
Landscapes
- Physics & Mathematics (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- General Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Engineering & Computer Science (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Measuring Magnetic Variables (AREA)
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
Abstract
【課題】磁場などの信号の測定精度を向上する。【解決手段】信号ベクトル導出装置1は、所定の方向を有するベクトルm(磁気双極子モーメント)により表される信号を受け、互いに直交する3軸X、Y、Zの成分Bx、By、Bzを測定する複数のセンサMS1~MS64の測定結果を受けてベクトルmの向きを導出する。信号ベクトル導出装置1は、センサの測定結果Bx、By、Bzと、第一係数vk1、vk2、vk3に第二係数ak、bk、ckを乗じて得られた値の和とに基づき得られるスペクトルであって、信号を出力する信号源S1、S2の存在するボクセルVで極大値をとるスペクトルを導出するスペクトル導出部16と、スペクトルを得るために用いられた第二係数ak、bk、ckに基づき、ベクトルの向きを導出する向き導出部18とを備える。【選択図】図2
Description
本発明は、磁場などの信号の測定に関する。
従来より、磁場の測定結果から信号源の位置を推定する方法(例えば、部分空間法(MUSIC法, SF, WSF等))が知られている(特許文献1、非特許文献1および非特許文献2を参照)。また、生体に関する磁場を測定することも知られている(特許文献2、特許文献3および特許文献4を参照)。
本発明は、磁場などの信号の測定精度の向上を課題とする。
本発明にかかる信号ベクトル導出装置は、所定の方向を有するベクトルにより表される信号を受け、互いに直交する3軸の成分を測定する複数のセンサの測定結果を受けて前記ベクトルの向きを導出するものであり、前記センサの測定結果が、前記ベクトルの前記3軸の成分の各々に第一係数を乗じたものの和に比例し、前記センサの測定結果と、前記第一係数に第二係数を乗じて得られた値の和とに基づき得られるスペクトルであって、前記信号を出力する信号源の存在するボクセルで極大値をとるスペクトルを導出するスペクトル導出部と、前記スペクトルを得るために用いられた前記第二係数に基づき、前記ベクトルの向きを導出する向き導出部とを備えるように構成される。
上記のように構成された信号ベクトル導出装置は、所定の方向を有するベクトルにより表される信号を受け、互いに直交する3軸の成分を測定する複数のセンサの測定結果を受けて前記ベクトルの向きを導出するものである。前記センサの測定結果が、前記ベクトルの前記3軸の成分の各々に第一係数を乗じたものの和に比例する。スペクトル導出部が、前記センサの測定結果と、前記第一係数に第二係数を乗じて得られた値の和とに基づき得られるスペクトルであって、前記信号を出力する信号源の存在するボクセルで極大値をとるスペクトルを導出する。向き導出部が、前記スペクトルを得るために用いられた前記第二係数に基づき、前記ベクトルの向きを導出する。
なお、本発明にかかる信号ベクトル導出装置は、前記ベクトルが、磁気双極子モーメントまたは電気双極子モーメントであるようにしてもよい。
なお、本発明にかかる信号ベクトル導出装置は、前記ベクトルが、電気双極子モーメントであり、前記測定結果の成分と同じ方向の前記ベクトルの成分が0であるようにしてもよい。
なお、本発明にかかる信号ベクトル導出装置は、前記第一係数が、前記ボクセルの各々と前記センサの各々との間の位置関係に基づき定められるようにしてもよい。
なお、本発明にかかる信号ベクトル導出装置は、前記第二係数のいずれか一つまたは二つが0であるようにしてもよい。
なお、本発明にかかる信号ベクトル導出装置は、前記スペクトルが、MUSIC法に則して求められるようにしてもよい。
なお、本発明にかかる信号ベクトル導出装置は、前記スペクトル導出部が、前記センサの測定結果から求められたノイズ部分空間の固有ベクトルに基づき前記スペクトルを導出するようにしてもよい。
なお、本発明にかかる信号ベクトル導出装置は、前記第一係数に第二係数を乗じて得られた値の和が、MUSIC法における伝達関数であるようにしてもよい。
なお、本発明にかかる信号ベクトル導出装置は、前記極大値が2つ以上存在するようにしてもよい。
なお、本発明にかかる信号ベクトル導出装置は、前記スペクトルに基づき、前記信号源の存在するボクセルの位置を導出する位置導出部を備えるようにしてもよい。
なお、本発明にかかる信号ベクトル導出装置は、前記位置導出部が、各ボクセルにおいて前記スペクトルの各々がとる値の最大値に基づき、前記信号源の存在するボクセルの位置を導出するようにしてもよい。
なお、本発明にかかる信号ベクトル導出装置は、前記位置導出部が、前記最大値の内でも最大の値から所定の範囲内の値をとるボクセルの重心を、前記所定の範囲を増やしながら、前記重心が所定量を超えて変化する回数に1を加えた値が前記信号源の個数になるまで、求め、前記重心を求めた前記ボクセルを、前記信号源の個数だけクラスタリングし、クラスタリングされた前記ボクセルの内で、前記スペクトルが最大の値をとるものを前記信号源の存在するボクセルの位置とするようにしてもよい。
なお、本発明にかかる信号ベクトル導出装置は、前記位置導出部が既に導出した前記信号源の存在するボクセルの位置に基づき、さらに、ボクセルの大きさを小さくして、前記信号源の存在するボクセルの位置を導出するようにしてもよい。
本発明は、所定の方向を有するベクトルにより表される信号を受け、互いに直交する3軸の成分を測定する複数のセンサの測定結果を受けて前記ベクトルの向きを導出する信号ベクトル導出方法であり、前記センサの測定結果が、前記ベクトルの前記3軸の成分の各々に第一係数を乗じたものの和に比例する信号ベクトル導出方法であって、前記センサの測定結果と、前記第一係数に第二係数を乗じて得られた値の和とに基づき得られるスペクトルであって、前記信号を出力する信号源の存在するボクセルで極大値をとるスペクトルを導出するスペクトル導出工程と、前記スペクトルを得るために用いられた前記第二係数に基づき、前記ベクトルの向きを導出する向き導出工程とを備えた信号ベクトル導出方法である。
本発明は、所定の方向を有するベクトルにより表される信号を受け、互いに直交する3軸の成分を測定する複数のセンサの測定結果を受けて前記ベクトルの向きを導出する信号ベクトル導出処理をコンピュータに実行させるためのプログラムであり、前記センサの測定結果が、前記ベクトルの前記3軸の成分の各々に第一係数を乗じたものの和に比例するプログラムであって、前記信号ベクトル導出処理が、前記センサの測定結果と、前記第一係数に第二係数を乗じて得られた値の和とに基づき得られるスペクトルであって、前記信号を出力する信号源の存在するボクセルで極大値をとるスペクトルを導出するスペクトル導出工程と、前記スペクトルを得るために用いられた前記第二係数に基づき、前記ベクトルの向きを導出する向き導出工程とを備えたプログラムである。
本発明は、所定の方向を有するベクトルにより表される信号を受け、互いに直交する3軸の成分を測定する複数のセンサの測定結果を受けて前記ベクトルの向きを導出する信号ベクトル導出処理をコンピュータに実行させるためのプログラムを記録したコンピュータによって読み取り可能な記録媒体であり、前記センサの測定結果が、前記ベクトルの前記3軸の成分の各々に第一係数を乗じたものの和に比例する記録媒体であって、前記信号ベクトル導出処理が、前記センサの測定結果と、前記第一係数に第二係数を乗じて得られた値の和とに基づき得られるスペクトルであって、前記信号を出力する信号源の存在するボクセルで極大値をとるスペクトルを導出するスペクトル導出工程と、前記スペクトルを得るために用いられた前記第二係数に基づき、前記ベクトルの向きを導出する向き導出工程とを備えた記録媒体である。
以下、本発明の実施形態を図面を参照しながら説明する。
図1は、本発明の実施形態にかかるボクセルVおよび磁気センサMSの斜視図である。図2は、本発明の実施形態にかかる信号ベクトル導出装置1の構成を示す機能ブロック図である。
図1を参照して、信号源S1および信号源S2が信号を出力する。信号は、所定の方向を有するベクトルmにより表される。ベクトルmは、例えば、磁気双極子モーメントである。なお、信号源の個数は、例えば2個であるが、磁気センサMSの個数未満であれば3個以上であってもよい。ただし、各信号源の出力する信号は、互いに周波数または位相が異なっているものとする。
また、信号源S1および信号源S2が存在する空間内の位置は、ボクセルV(例えば、10×10×10=1000個のボクセル)により表される。信号源S1および信号源S2は、それぞれ異なるボクセルVに位置している。なお、1000個のボクセルVの各々を、V1~V1000と表記する。
複数(例えば、8行8列の64個)の磁気センサMSは、信号(例えば、磁気双極子モーメント)を受け、互いに直交する3軸X、Y、Zの成分Bx、By、Bzを測定する。なお、64個の磁気センサMSの各々を、MS1~MS64と表記する。
ここで、ベクトルrを、信号源(磁気双極子)から磁気センサMSまでの方向ベクトルとすると、磁気センサMSにより測定される磁束密度B(ベクトルrの関数)は、ビオサバールの法則により、式(1)のように表される。ただし、μ0は、磁気定数である。また、ベクトルrは、ボクセルVの各々(V1~V1000)と、磁気センサMSの各々MS1~MS64との間の位置関係といえる。
式(1)より、Byは、以下の式(3)のように表される。
式(1)より、Bzは、以下の式(4)のように表される。
なお、vx1、vx2、vx3、vy1、vy2、vy3、vz1、vz2、vz3(第一係数)は、式(2)~(4)および式(2’)~(4’)を参照して、ベクトルrに基づき定められている。
図2を参照して、本発明の実施形態にかかる信号ベクトル導出装置1は、相対位置記録部11、第一係数導出部12、伝達関数導出部13、雑音固有ベクトル導出部14、スペクトル導出部16、向き導出部18、位置導出部19を備える。
信号ベクトル導出装置1は、複数のセンサMS1~MS64の測定結果を受けてベクトルmの向きを導出する。
相対位置記録部11は、ボクセルVの各々1000個と、磁気センサMSの各々MS1~MS64との間の相対位置であるベクトルrを記録する。
第一係数導出部12は、相対位置記録部11からベクトルrを読み出し、第一係数vx1、vx2、vx3、vy1、vy2、vy3、vz1、vz2、vz3を導出する(式(2)~(4)および式(2’)~(4’)を参照)。
例えば、第一係数vx1は、以下の式(5)のように表される。
他の第一係数vx2、vx3、vy1、vy2、vy3、vz1、vz2、vz3も同様に、1000×64通りの値をとる。
なお、第一係数vx1は、そのまま用いることも考えられるが、以下の式(6)のように正規化して、以後の処理に使用する。
他の第一係数vx2、vx3、vy1、vy2、vy3、vz1、vz2、vz3も同様に、正規化する。
第一係数導出部12は、上記のようにして正規化した第一係数を出力する。
雑音固有ベクトル導出部14は、磁気センサMSの測定結果Bx、By、Bzから、MUSIC法に則して、ノイズ部分空間の固有ベクトルを求める。
まず、磁気センサMSの測定結果Bxから以下の式(7)のように、X(t)xを求める。ただし、tは測定を行った時間である。Tは転置を意味する。
X(t)xを用いて、以下の式(8)のように相関行列を求める。
なお、ノイズ部分空間の固有ベクトルeyも同様に求めることができる。まず、式(7)のBxをByに置き換え、式(7)および式(8)のX(t)xをX(t)yに置き換えて、式(8)により相関行列を求める。あとは、同様に、小さな値の固有値に対応する固有ベクトルを求め、ノイズ部分空間の固有ベクトルeyとする。ノイズ部分空間の固有ベクトルeyは、64行1列のベクトルである。ノイズ部分空間の固有ベクトルeyは、小さな値の固有値に対応して62個存在する。
また、ノイズ部分空間の固有ベクトルezも同様に求めることができる。まず、式(7)のBxをBzに置き換え、式(7)および式(8)のX(t)xをX(t)zに置き換えて、式(8)により相関行列を求める。あとは、同様に、小さな値の固有値に対応する固有ベクトルを求め、ノイズ部分空間の固有ベクトルezとする。ノイズ部分空間の固有ベクトルezは、64行1列のベクトルである。ノイズ部分空間の固有ベクトルezは、小さな値の固有値に対応して62個存在する。
伝達関数導出部13は、以下の式(9)、式(10)および式(11)のように、伝達関数vx、vyおよびvzを導出する。第一係数vx1、vx2、vx3に、それぞれ第二係数ax、bx、cxを乗じて得られた値の和を導出する(式(9)参照)。導出結果を、伝達関数vxという。第一係数vy1、vy2、vy3に、それぞれ第二係数ay、by、cyを乗じて得られた値の和を導出する(式(10)参照)。導出結果を、伝達関数vyという。第一係数vz1、vz2、vz3に、それぞれ第二係数az、bz、czを乗じて得られた値の和を導出する(式(11)参照)。導出結果を、伝達関数vzという。
また、第二係数は、いずれも0以外の値であってもよい。例えば、ax=bx=cx=1(すなわち、vx=vx1+vx2+vx3)でもよいし、ay=by=1、cy=-1(すなわち、vy=vy1+vy2-vy3)でもよいし、az=1、bz=-1、cz=1(すなわち、vz=vz1-vz2+vz3)でもよい。
ただし、第二係数のいずれか一つまたは二つが0であってもよい。例えば、ax=1、bx=cx=0(すなわち、vx=vx1)でもよいし、ax=bx=1、cx=0(すなわち、vx=vx1+vx2)でもよい。
ここで、(ak、bk、ck)(ただし、k=x、y、z)が、(1,0,0)、(0,1,0)、(0,0,1)、(1,1,0)、(1,-1,0)、(0,1,1)、(0,1,-1)、(1,0,1)、(-1,0,1)、(1,1,1)、(-1,1,1)、(1,-1,1)、(1,1,-1)の13種類あるものとする。これらに対応するvkをそれぞれ、vk1、vk2、…、vk13とする。
例えば、(ax、bx、cx)が、(1,0,0)である場合、vx=vx1である。(ax、bx、cx)が、(1,1,0)である場合、vx=vx4=vx1+vx2である。(ax、bx、cx)が、(1,1,-1)である場合、vx=vx13=vx1+vx2-vx3である。
例えば、(ay、by、cy)が、(1,0,0)である場合、vy=vy1である。(ay、by、cy)が、(1,1,0)である場合、vy=vy4=vy1+vy2である。(ay、by、cy)が、(1,1,-1)である場合、vy=vy13=vy1+vy2-vy3である。
例えば、(az、bz、cz)が、(1,0,0)である場合、vz=vz1である。(az、bz、cz)が、(1,1,0)である場合、vz=vz4=vz1+vz2である。(az、bz、cz)が、(1,1,-1)である場合、vz=vz13=vz1+vz2-vz3である。
スペクトル導出部16は、信号源S1、S2の存在するボクセルVで極大値をとるスペクトルを導出する。このようなスペクトルは、MUSIC法に則して求められる。スペクトルの極大値は、信号源の個数に対応して、2つある。なお、信号源の個数が3つ以上あれば、スペクトルの極大値も3つ以上ある。
スペクトルは、磁気センサMSの測定結果Bx、By、Bz(から求められたノイズ部分空間の固有ベクトルex、ey、ez)と、第一係数に第二係数を乗じて得られた値の和(すなわち、伝達関数vx、vy、vz)(式(9)、(10)、(11))とに基づき、スペクトル導出部16により導出される。スペクトル導出部16は、伝達関数導出部13の出力する伝達関数vx、vy、vzと、雑音固有ベクトル導出部14の出力するノイズ部分空間の固有ベクトルex、ey、ezとに基づきスペクトルを導出する。
スペクトル導出部16は、以下のようにして、スペクトルPx1を導出する。
(1)ノイズ部分空間の固有ベクトルex(64行1列のベクトル)は、62個存在するが、これらのベクトルexを62列並べて、exを64行62列の行列とする。
(2)行列exを転置し、伝達関数vx1(64行1000列の行列)に乗じる。すなわち、exTvx1を求める。これは、62行1000列の行列となる。
(3)(2)で得た行列の各成分を2乗する。
(4)(3)で得た行列の成分を列ごとに合計し、1行に並べて、1行1000列の行列を得る。例えば、(3)で得た行列の1行Q列目+2行Q列目+…+62行Q列目が、(4)で得られる1行1000列の行列の1行Q列目の成分となる(ただし、Qは1~1000の整数)。
(5)(4)で得た行列の各成分の逆数をとると、スペクトルPx1(1行1000列の行列)が得られる。
(1)ノイズ部分空間の固有ベクトルex(64行1列のベクトル)は、62個存在するが、これらのベクトルexを62列並べて、exを64行62列の行列とする。
(2)行列exを転置し、伝達関数vx1(64行1000列の行列)に乗じる。すなわち、exTvx1を求める。これは、62行1000列の行列となる。
(3)(2)で得た行列の各成分を2乗する。
(4)(3)で得た行列の成分を列ごとに合計し、1行に並べて、1行1000列の行列を得る。例えば、(3)で得た行列の1行Q列目+2行Q列目+…+62行Q列目が、(4)で得られる1行1000列の行列の1行Q列目の成分となる(ただし、Qは1~1000の整数)。
(5)(4)で得た行列の各成分の逆数をとると、スペクトルPx1(1行1000列の行列)が得られる。
なお、スペクトルPx1の各列は、ボクセルV1~V1000に対応する。他のスペクトルも同様である。
スペクトル導出部16は、スペクトルPx2、Px3、…、Px13も導出する。スペクトルPx2、Px3、…、Px13は、上記(2)における伝達関数vx1を、vx2、vx3、…、vx13に置き換えれば、導出することができる。
スペクトル導出部16は、スペクトルPy1、Py2、Py3、…、Py13を導出する。上記(1)におけるノイズ部分空間の固有ベクトルexをeyに置き換え、上記(2)における伝達関数vx1を、vy1、vy2、vy3、…、vy13に置き換えれば、スペクトルPy1、Py2、Py3、…、Py13を導出することができる。
スペクトル導出部16は、スペクトルPz1、Pz2、Pz3、…、Pz13を導出する。上記(1)におけるノイズ部分空間の固有ベクトルexをezに置き換え、上記(2)における伝達関数vx1を、vz1、vz2、vz3、…、vz13に置き換えれば、スペクトルPz1、Pz2、Pz3、…、Pz13を導出することができる。
位置導出部19は、スペクトルPx1、Px2、Px3、…、Px13、Py1、Py2、Py3、…、Py13、Pz1、Pz2、Pz3、…、Pz13に基づき、信号源S1、S2の存在するボクセルVの位置を導出する。
スペクトル導出部16の出力するスペクトルは、以下の式(12)のように表される。
図3は、最大値Pのグラフの一例である。図3において、縦軸はスペクトルの値、横軸はボクセル(V1~V1000)を示す。
図3を参照して、ボクセルV750において最大値P(のスペクトル)の値がSP1(極大値)、ボクセルV250において最大値P(のスペクトル)の値がSP2(極大値)をとるものとする。ただし、SP1がSP2よりも大きいものとする。
位置導出部19は、最大値Pの内でも最大の値SP1から所定の範囲内(例えば、最大値Pの値が0.95SP1以上)の値をとるボクセルの重心を、所定の範囲を増やしながら(例えば、最大値Pの値が0.95SP1以上→0.90SP1以上→0.85SP1以上→…というように0.05SP1ずつ所定の範囲を広げる)、重心が所定量を超えて変化する回数に1を加えた値が信号源の個数(2個)になるまで、求める。
最大値Pの内でも最大の値SP1から近傍におけるボクセルの重心は、おおむねボクセルV750である。しかし、最大の値SP1から所定の範囲が広がって、SP2を含むようになると、ボクセルの重心がボクセルV750から、かなり小さくなる。すると、重心が所定量を超えて変化するので、その回数(1回)に1を加えた値が信号源の個数(2個)になるので、ボクセルの重心を求めることを終了する。
次に、重心を求めたボクセルを、信号源の個数(2個)だけクラスタリングする。例えば、教師なし機械学習のKmeansクラスタリングを行い、信号源の数だけラベリングする。
最後に、クラスタリングされたボクセルの内で、スペクトルが最大の値をとるものを信号源の存在するボクセルの位置とする。
なお、位置導出部19は、このようにして既に導出した信号源の存在するボクセルの位置に基づき、さらに、ボクセルの大きさを小さくして、信号源の存在するボクセルの位置を導出するようにしてもよい。このようにして、信号源の存在するボクセルの位置を、高精度かつ高速に計算することができる。
向き導出部18は、位置導出部19から信号源S1、S2に対応する(すなわち、Pにおいて極大値をとる)Pkj(ただし、k=x、y、zかつj=1、2、3、…)を受ける。さらに、向き導出部18は、信号源S1、S2に対応するPkjを得るために用いられた第二係数に基づき、ベクトルmの向きを導出する。
例えば、向き導出部18に、位置導出部19から信号源S1およびS2に対応するスペクトルとして、750列目(ボクセルV750)のPx13(Py13またはPz13)および250列目(ボクセルV250)のPx4(Px4またはPz4)が与えられたとする。
すると、向き導出部18は、ボクセルV750の信号源S1におけるベクトルmの向きとして、Px13(Py13またはPz13)を得るために用いられた第二係数(ak、bk、ck)(ただし、k=x、y、z)は、(1,1,-1)である。よって、向き導出部18は、ボクセルV750の信号源S1におけるベクトルmの向きを、ベクトル(1,1,-1)に平行であると導出する。ただし、ベクトル(1,1,-1)は、X成分1、Y成分1、Z成分-1のベクトルである。
さらに、向き導出部18は、ボクセルV250の信号源S2におけるベクトルmの向きとして、Px4(Px4またはPz4)を得るために用いられた第二係数(ak、bk、ck)(ただし、k=x、y、z)は、(1,1,0)である。よって、向き導出部18は、ボクセルV250の信号源S2におけるベクトルmの向きを、ベクトル(1,1,0)に平行であると導出する。ただし、ベクトル(1,1,0)は、X成分1、Y成分1、Z成分0のベクトルである。
次に、本発明の実施形態の動作を説明する。
第一係数導出部12により、相対位置記録部11からベクトルrが読み出され、第一係数vx1、vx2、vx3、vy1、vy2、vy3、vz1、vz2、vz3が導出される(式(2)~(4)および式(2’)~(4’)を参照)。
なお、第一係数は、1000×64通りの値をとり(式(5)参照)、正規化されて(式(6)参照)、伝達関数導出部13に与えられる。
伝達関数導出部13によって、第一係数および第二係数ax、bx、cx、ay、by、cy、az、bz、czに基づき、伝達関数vx、vy、vzが導出される(式(9)、式(10)および式(11)参照)。
雑音固有ベクトル導出部14によって、磁気センサMSの測定結果Bx、By、Bzから、MUSIC法に則して、ノイズ部分空間の固有ベクトルex、ey、ezが導出される。
スペクトル導出部16により、伝達関数vx、vy、vzおよびノイズ部分空間の固有ベクトルex、ey、ezに基づき、スペクトルPx1、Px2、Px3、…、Px13、Py1、Py2、Py3、…、Py13、Pz1、Pz2、Pz3、…、Pz13が導出される(式(12)参照)。
位置導出部19により、各ボクセルにおいてスペクトルの各々がとる値の最大値P(すなわち、式(12)の各列における最大値)が求められる(式(13)および図3を参照)。最大値Pに基づき、信号源S1、S2の存在するボクセル250、750が導出される。
向き導出部18により、信号源S1、S2に対応するPkjを得るために用いられた第二係数に基づき、ベクトルmの向きが導出される。
本発明の実施形態によれば、磁場などの信号の測定精度が向上する。
例えば、伝達関数vkが、vk1、vk2、vk3のみである場合、ベクトルmの向きが、X方向、Y方向またはZ方向に平行な場合しか、ベクトルmの向きを測定できない。ベクトルmの向きが、それ以外の向き、例えば、ベクトル(1,1,0)(すなわち、X成分1、Y成分1、Z成分0のベクトル)に平行である場合は、ベクトルmの向きを測定できない。
しかし、本発明の実施形態によれば、伝達関数vkが、vk1、vk2、vk3、…、vk13と多くの種類があるので、ベクトルmの向きが、X方向、Y方向およびZ方向に平行ではない場合であっても、ベクトルmの向きを測定できる。
なお、本発明の実施形態においては、信号のベクトルを磁気双極子モーメントであるとしてきたが、信号のベクトルは磁気双極子モーメントに限定されない。信号のベクトルは、例えば、電気双極子モーメント(ベクトルp)であってもよい。
磁気センサMSにより測定される磁束密度B(ベクトルrの関数)は、式(14)のように表される。
式(14)より、Byは、以下の式(16)のように表される。
式(14)より、Bzは、以下の式(17)のように表される。
信号ベクトル導出装置1の構成および動作は、信号のベクトルが磁気双極子モーメント(ベクトルm)の場合と同様であり、説明を省略する。
また、上記の実施形態は、以下のようにして実現できる。CPU、ハードディスク、メディア(USBメモリ、CD-ROMなど)読み取り装置を備えたコンピュータに、上記の各部分、例えば相対位置記録部11、第一係数導出部12、伝達関数導出部13、雑音固有ベクトル導出部14、スペクトル導出部16、向き導出部18および位置導出部19を実現するプログラムを記録したメディアを読み取らせて、ハードディスクにインストールする。このような方法でも、上記の機能を実現できる。
1 信号ベクトル導出装置
11 相対位置記録部
12 第一係数導出部
13 伝達関数導出部
14 雑音固有ベクトル導出部
16 スペクトル導出部
18 向き導出部
19 位置導出部
MS 磁気センサ
V ボクセル
B 磁束密度
vx1、vx2、vx3、vy1、vy2、vy3、vz1、vz2、vz3 第一係数
ak、bk、ck 第二係数
vk1、vk2、…、vk13 伝達関数
S1、S2 信号源
m ベクトル(磁気双極子モーメント)
ex、ey、ez ノイズ部分空間の固有ベクトル
P 最大値
11 相対位置記録部
12 第一係数導出部
13 伝達関数導出部
14 雑音固有ベクトル導出部
16 スペクトル導出部
18 向き導出部
19 位置導出部
MS 磁気センサ
V ボクセル
B 磁束密度
vx1、vx2、vx3、vy1、vy2、vy3、vz1、vz2、vz3 第一係数
ak、bk、ck 第二係数
vk1、vk2、…、vk13 伝達関数
S1、S2 信号源
m ベクトル(磁気双極子モーメント)
ex、ey、ez ノイズ部分空間の固有ベクトル
P 最大値
Claims (16)
- 所定の方向を有するベクトルにより表される信号を受け、互いに直交する3軸の成分を測定する複数のセンサの測定結果を受けて前記ベクトルの向きを導出する信号ベクトル導出装置であり、前記センサの測定結果が、前記ベクトルの前記3軸の成分の各々に第一係数を乗じたものの和に比例する信号ベクトル導出装置であって、
前記センサの測定結果と、前記第一係数に第二係数を乗じて得られた値の和とに基づき得られるスペクトルであって、前記信号を出力する信号源の存在するボクセルで極大値をとるスペクトルを導出するスペクトル導出部と、
前記スペクトルを得るために用いられた前記第二係数に基づき、前記ベクトルの向きを導出する向き導出部と、
を備えた信号ベクトル導出装置。 - 請求項1に記載の信号ベクトル導出装置であって、
前記ベクトルが、磁気双極子モーメントまたは電気双極子モーメントである、
信号ベクトル導出装置。 - 請求項1に記載の信号ベクトル導出装置であって、
前記ベクトルが、電気双極子モーメントであり、
前記測定結果の成分と同じ方向の前記ベクトルの成分が0である、
信号ベクトル導出装置。 - 請求項1に記載の信号ベクトル導出装置であって、
前記第一係数が、前記ボクセルの各々と前記センサの各々との間の位置関係に基づき定められる、
信号ベクトル導出装置。 - 請求項1に記載の信号ベクトル導出装置であって、
前記第二係数のいずれか一つまたは二つが0である、
信号ベクトル導出装置。 - 請求項1に記載の信号ベクトル導出装置であって、
前記スペクトルが、MUSIC法に則して求められる、
信号ベクトル導出装置。 - 請求項6に記載の信号ベクトル導出装置であって、
前記スペクトル導出部が、前記センサの測定結果から求められたノイズ部分空間の固有ベクトルに基づき前記スペクトルを導出する、
信号ベクトル導出装置。 - 請求項6に記載の信号ベクトル導出装置であって、
前記第一係数に第二係数を乗じて得られた値の和が、MUSIC法における伝達関数である、
信号ベクトル導出装置。 - 請求項1に記載の信号ベクトル導出装置であって、
前記極大値が2つ以上存在する、
信号ベクトル導出装置。 - 請求項1に記載の信号ベクトル導出装置であって、
前記スペクトルに基づき、前記信号源の存在するボクセルの位置を導出する位置導出部を備えた信号ベクトル導出装置。 - 請求項10に記載の信号ベクトル導出装置であって、
前記位置導出部が、各ボクセルにおいて前記スペクトルの各々がとる値の最大値に基づき、前記信号源の存在するボクセルの位置を導出する信号ベクトル導出装置。 - 請求項11に記載の信号ベクトル導出装置であって、
前記位置導出部が、
前記最大値の内でも最大の値から所定の範囲内の値をとるボクセルの重心を、前記所定の範囲を増やしながら、前記重心が所定量を超えて変化する回数に1を加えた値が前記信号源の個数になるまで、求め、
前記重心を求めた前記ボクセルを、前記信号源の個数だけクラスタリングし、
クラスタリングされた前記ボクセルの内で、前記スペクトルが最大の値をとるものを前記信号源の存在するボクセルの位置とする、
信号ベクトル導出装置。 - 請求項10に記載の信号ベクトル導出装置であって、
前記位置導出部が既に導出した前記信号源の存在するボクセルの位置に基づき、さらに、ボクセルの大きさを小さくして、前記信号源の存在するボクセルの位置を導出する、
信号ベクトル導出装置。 - 所定の方向を有するベクトルにより表される信号を受け、互いに直交する3軸の成分を測定する複数のセンサの測定結果を受けて前記ベクトルの向きを導出する信号ベクトル導出方法であり、前記センサの測定結果が、前記ベクトルの前記3軸の成分の各々に第一係数を乗じたものの和に比例する信号ベクトル導出方法であって、
前記センサの測定結果と、前記第一係数に第二係数を乗じて得られた値の和とに基づき得られるスペクトルであって、前記信号を出力する信号源の存在するボクセルで極大値をとるスペクトルを導出するスペクトル導出工程と、
前記スペクトルを得るために用いられた前記第二係数に基づき、前記ベクトルの向きを導出する向き導出工程と、
を備えた信号ベクトル導出方法。 - 所定の方向を有するベクトルにより表される信号を受け、互いに直交する3軸の成分を測定する複数のセンサの測定結果を受けて前記ベクトルの向きを導出する信号ベクトル導出処理をコンピュータに実行させるためのプログラムであり、前記センサの測定結果が、前記ベクトルの前記3軸の成分の各々に第一係数を乗じたものの和に比例するプログラムであって、
前記信号ベクトル導出処理が、
前記センサの測定結果と、前記第一係数に第二係数を乗じて得られた値の和とに基づき得られるスペクトルであって、前記信号を出力する信号源の存在するボクセルで極大値をとるスペクトルを導出するスペクトル導出工程と、
前記スペクトルを得るために用いられた前記第二係数に基づき、前記ベクトルの向きを導出する向き導出工程と、
を備えたプログラム。 - 所定の方向を有するベクトルにより表される信号を受け、互いに直交する3軸の成分を測定する複数のセンサの測定結果を受けて前記ベクトルの向きを導出する信号ベクトル導出処理をコンピュータに実行させるためのプログラムを記録したコンピュータによって読み取り可能な記録媒体であり、前記センサの測定結果が、前記ベクトルの前記3軸の成分の各々に第一係数を乗じたものの和に比例する記録媒体であって、
前記信号ベクトル導出処理が、
前記センサの測定結果と、前記第一係数に第二係数を乗じて得られた値の和とに基づき得られるスペクトルであって、前記信号を出力する信号源の存在するボクセルで極大値をとるスペクトルを導出するスペクトル導出工程と、
前記スペクトルを得るために用いられた前記第二係数に基づき、前記ベクトルの向きを導出する向き導出工程と、
を備えた記録媒体。
Priority Applications (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2021074991A JP2022169141A (ja) | 2021-04-27 | 2021-04-27 | 信号ベクトル導出装置、方法、プログラム、記録媒体 |
US18/272,888 US20240094311A1 (en) | 2021-04-27 | 2022-03-07 | Signal vector derivation apparatus, method, program, and recording medium |
PCT/JP2022/009613 WO2022230365A1 (ja) | 2021-04-27 | 2022-03-07 | 信号ベクトル導出装置、方法、プログラム、記録媒体 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2021074991A JP2022169141A (ja) | 2021-04-27 | 2021-04-27 | 信号ベクトル導出装置、方法、プログラム、記録媒体 |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2022169141A true JP2022169141A (ja) | 2022-11-09 |
Family
ID=83848295
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2021074991A Pending JP2022169141A (ja) | 2021-04-27 | 2021-04-27 | 信号ベクトル導出装置、方法、プログラム、記録媒体 |
Country Status (3)
Country | Link |
---|---|
US (1) | US20240094311A1 (ja) |
JP (1) | JP2022169141A (ja) |
WO (1) | WO2022230365A1 (ja) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2023105885A1 (ja) * | 2021-12-09 | 2023-06-15 | 株式会社アドバンテスト | 信号源特定装置、方法、プログラム、記録媒体 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2003038455A (ja) * | 2001-05-22 | 2003-02-12 | Shimadzu Corp | 生体活動解析装置 |
WO2005025416A2 (en) * | 2003-09-11 | 2005-03-24 | Regents Of The University Of Minnesota | Localizing neural sources in a brain |
FR2893434B1 (fr) * | 2005-11-17 | 2008-05-09 | Univ Rennes I Etablissement Pu | Procede et dispositif d'identification de parametre multidimensionnels : application a la localisation et la reconstruction d'activites electriques de profondeur au moyen d'observations de surfaces |
CN106970348B (zh) * | 2017-02-22 | 2019-08-06 | 西安电子科技大学 | 电磁矢量传感器阵列解相干二维music参数估计方法 |
-
2021
- 2021-04-27 JP JP2021074991A patent/JP2022169141A/ja active Pending
-
2022
- 2022-03-07 US US18/272,888 patent/US20240094311A1/en active Pending
- 2022-03-07 WO PCT/JP2022/009613 patent/WO2022230365A1/ja active Application Filing
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2023105885A1 (ja) * | 2021-12-09 | 2023-06-15 | 株式会社アドバンテスト | 信号源特定装置、方法、プログラム、記録媒体 |
Also Published As
Publication number | Publication date |
---|---|
US20240094311A1 (en) | 2024-03-21 |
WO2022230365A1 (ja) | 2022-11-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US7463024B2 (en) | Method and device for processing a multi-channel measurement of magnetic fields | |
Taulu et al. | Presentation of electromagnetic multichannel data: the signal space separation method | |
JP7402768B2 (ja) | 磁場計測装置、磁場計測方法、磁場計測プログラム | |
WO2022230365A1 (ja) | 信号ベクトル導出装置、方法、プログラム、記録媒体 | |
WO2020138170A1 (ja) | 磁場計測装置 | |
Nolte et al. | The effect of artifact rejection by signal-space projection on source localization accuracy in MEG measurements | |
US20170059318A1 (en) | Method and apparatus of correcting output value of geomagnetic sensor | |
JP2018004286A (ja) | 信号処理装置、信号処理方法、信号処理プログラム、及び磁場計測システム | |
Ma et al. | Improved SENSE imaging using accurate coil sensitivity maps generated by a global magnitude‐phase fitting method | |
JPH04117944A (ja) | 核磁気共鳴装置 | |
Bhayani et al. | Partially calibrated semi-generalized pose from hybrid point correspondences | |
Wijnholds et al. | Fundamental imaging limits of radio telescope arrays | |
Nolte et al. | Partial signal space projection for artefact removal in MEG measurements: a theoretical analysis | |
WO2023203820A1 (ja) | 画像出力装置、方法、プログラム、記録媒体 | |
Costa et al. | Using the axis of rotation of polar navigator echoes to rapidly measure 3D rigid‐body motion | |
CN114492082B (zh) | 光栅投影成像***的光栅相位提取方法 | |
JP2023147067A (ja) | 磁場計測装置、磁場計測方法、及び磁場計測プログラム | |
WO2023105885A1 (ja) | 信号源特定装置、方法、プログラム、記録媒体 | |
US10997462B2 (en) | Devices, systems, and methods for clustering reference images for non-destructive testing | |
Di Rienzo et al. | Theoretical lower error bound for comparative evaluation of sensor arrays in magnetostatic linear inverse problems | |
Kamondetdacha et al. | Multipole characterization of a magnetic source using a truncated SVD | |
US11268805B2 (en) | Measurement method | |
JP4775632B2 (ja) | Squid顕微鏡による画像の解析方法、及びsquid顕微鏡による画像の解析システム | |
US20160084630A1 (en) | Calculation method, storage medium, information processing apparatus, and measurement apparatus | |
EP4109044A1 (en) | System and a method for analyzing the motions of a mechanical structure |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20240209 |