JP3067728B2 - Biomagnetic field measurement device - Google Patents

Description

【発明の詳細な説明】DETAILED DESCRIPTION OF THE INVENTION

【０００１】[0001]

【発明の属する技術分野】本発明は、生体の脳の神経活
動、心臓の心筋活動等により発生する生体磁場を、高感
度な量子干渉素子（ＳＱＵＩＤ:superconducting quant
um interferencedevice）からなる複数の磁束計を用い
て計測する生体磁場計測方法及び生体磁場計測装置に関
する。BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to a highly sensitive quantum interference device (SQUID: superconducting quant
The present invention relates to a biomagnetic field measurement method and a biomagnetic field measurement apparatus for performing measurement using a plurality of magnetometers composed of an um interference device.

【０００２】[0002]

【従来の技術】本発明は、生体の脳の神経活動、心臓の
心筋活動等により発生する生体磁場を、高感度な量子干
渉素子（ＳＱＵＩＤ：ｓｕｐｅｒｃｏｎｄｕｃｔｉｎｇ
ｑｕａｎｔｕｍ ｉｎｔｅｒｆｅｒｅｎｃｅ ｄｅｖ
ｉｃｅ）からなる複数の磁束計を用いて計測する生体磁
場計測方法及び生体磁場計測装置に関する。2. Description of the Related Art The present invention relates to a highly sensitive quantum interference device (SQUID: superconducting element) for generating a biomagnetic field generated by neural activity of the brain of a living body or cardiac muscle activity of the heart.
quantum interference dev
The present invention relates to a biomagnetic field measurement method and a biomagnetic field measurement apparatus for performing measurement using a plurality of magnetometers composed of a plurality of pieces of ice.

【０００３】生体磁場としては、電流ダイポールが作り
出す磁場の他、生体内を流れる体積電流による磁場があ
る。生体磁場の法線成分（Ｂ_z（直交座標系でのＺ成
分）又はＢ_r（極座標系での動径成分））の計測は、体
積電流の影響を受けにくいと考えられている。従来技術
では、生体表面に対してＳＱＵＩＤに接続した検出コイ
ルの面を平行に配置して、生体表面に垂直な法線成分で
あるＢ_z又はＢ_rを計測していた。生体磁場計測の結果
は、測定された磁場成分の時間変化波形、測定された磁
場成分の任意の時点での強度の等しい点を結ぶ等磁場線
図（コンターマップ）により表示されていた。また、得
られた等磁場曲線から、生体磁場を発生している磁場源
を解析する種々の解析方法が提案されているが、代表的
な解析方法では磁場源を電流ダイポールに置き換えて解
析を行なっていた。As a biomagnetic field, there is a magnetic field generated by a volume current flowing in a living body, in addition to a magnetic field generated by a current dipole. It is considered that the measurement of the normal component (B _z (Z component in the orthogonal coordinate system) or _Br (radial component in the polar coordinate system)) of the biomagnetic field is hardly affected by the volume current. In the prior art, the surface of the detection coil connected to the SQUID to a living body surface is disposed in parallel, was measured B _z or B _r is a perpendicular normal component to the surface of the living body. The result of the biomagnetic field measurement was displayed by a time-varying waveform of the measured magnetic field component and an isomagnetic field map (contour map) connecting points of the measured magnetic field component having the same intensity at an arbitrary point in time. In addition, various analysis methods have been proposed to analyze the magnetic field source that is generating the biomagnetic field from the obtained isomagnetic field curves.A typical analysis method is to perform analysis by replacing the magnetic field source with a current dipole. I was

【０００４】電流ダイポールが作る磁場の法線成分（Ｂ
_z又はＢ_r）の等磁場線図は、磁場源（電流ダイポール）
を中心として分離した位置に磁場の沸き出し極、磁場の
吸い込み極を持つパターンとなる。この２つの極での磁
場強度、２つの極の間の距離により、磁場源（電流ダイ
ポール）の大きさ、位置、方向等が解析されている。The normal component (B) of the magnetic field created by the current dipole
_z or B _r ) is a magnetic field source (current dipole)
Is a pattern having a magnetic field evacuation pole and a magnetic field suction pole at positions separated from the center. The magnitude, position, direction, and the like of the magnetic field source (current dipole) are analyzed based on the magnetic field strength at the two poles and the distance between the two poles.

【０００５】第１の従来技術（Ｈ．Ｈｏｓａｋａ ａｎ
ｄ Ｄ．Ｃｏｈｅｎ、Ｊ．Ｅｌｅｃｔｒｏｃａｒｄｉｏ
ｌ．、９（４）、４２６−４３２（１９７６））では、
心筋内の電流の方向や強さを見易くするため、計測され
た法線成分Ｂ_zの等磁場線図を用いて、心筋に分布する
電流源を表示する方法として、（数１）で定義される電
流ベクトル〈Ｊ（ｘ、ｙ）〉を各計測点上に矢印で表現
するアローマップが考案されている。なお以下の説明で
は、括弧〈 〉は〈 〉内の記号がベクトルであること
を示し、例えば、〈Ｊ〉はＪがベクトルであることを表
わす。[0005] The first prior art (H. Hosaka an
d D. Cohen, J .; Electrocardio
l. , 9 (4), 426-432 (1976))
For clarity the direction and intensity of the current in the myocardium, with the isomagnetic field of the measured normal component B _z, as a method for displaying a current source distributed in the myocardium, is defined by (Equation 1) An arrow map has been devised in which a current vector <J (x, y)> is expressed by an arrow on each measurement point. In the following description, parentheses <> indicate that the symbol in <> is a vector, for example, <J> indicates that J is a vector.

【０００６】 〈Ｊ（ｘ、ｙ）〉 ＝（∂Ｂ_z（ｘ、ｙ）／∂ｙ）〈ｅ_x〉−（∂Ｂ_z（ｘ、ｙ）／∂ｘ）〈ｅ_y〉 …（数１） （数１）に於いて、〈ｅ_x〉はｘ方向の単位ベクトル、
〈ｅ_y〉はｙ方向の単位ベクトルである。しかし、複数
の電流源が存在する時には、法線成分Ｂ_zの等磁場線図
から個々の電流源を判別しにくいという問題があった。[0006] <J (x, y)> = (∂B z (x, y) / ∂y) <e x> - (∂B z (x, y) / ∂x) <e y> ... ( number 1) (in the several 1), <e _x> is the x-direction of the unit vector,
<E _y > is a unit vector in the y direction. However, when a plurality of current sources are present, there is a problem that the isomagnetic field of the normal component B _z difficult to determine the individual current sources.

【０００７】第２の従来技術（Ｋ．Ｔｕｋａｄａ ｅｔ
ａｌ．、Ｒｅｖｅｉｗ ｏｆ ｔｈｅ Ｓｃｉｅｎｔ
ｉｆｉｃ Ｉｎｓｔｒｕｍｅｎｔｓ、６６（１０）５０
８５−５０９１（１９９５））では、分布する複数の電
流源を可視化するために、法線成分（Ｂ_z又はＢ_r）を計
測するのではなく、検出コイルの面を生体表面に対して
垂直に配置して、接線成分Ｂ_x及びＢ_yを計測している。
計測された接線成分Ｂ_x、Ｂ_yを各成分毎に等磁場線図と
して表示している。従来技術２で計測された接線成分Ｂ
_x、Ｂ_yは体積電流の影響が考えられるものの、（数２）
に従って、時刻ｔに於いて計測されたＢ_x及びＢ_yを合成
した２次元ベクトル強度Ｂ_xyの等磁場線図では、常に電
流ダイポールの直上にピークが得られることから、複数
の電流ダイポールが存在する場合でも、各電流ダイポー
ルを分離して可視化できる。A second prior art (K. Tukada et. Al.)
al. , Review of the Scientific
ific Instruments, 66 (10) 50
In 85-5091 (1995)), in order to visualize a plurality of current sources distributed, rather than measure the normal component (B _z or B _r), perpendicular surfaces of the detection coil with respect to the living body surface arranged to, measures the tangential components B _x and B _y.
It measured tangential components B _x, and a B _y displayed as isomagnetic field for each component. Tangent component B measured in prior art 2
_x, although B _y is considered the influence of the volume current, (Equation 2)
Accordingly the isomagnetic field of a two-dimensional vector magnitude B _xy obtained by combining the measured B _x and B _y at time t, always from a peak immediately above the current dipole can be obtained, there are a plurality of current dipoles In this case, each current dipole can be separated and visualized.

【０００８】 │Ｂ_xy（ｘ、ｙ、ｔ）│ ＝√｛（Ｂ_x（ｘ、ｙ、ｔ））²＋（Ｂ_y（ｘ、ｙ、ｔ））²｝ …（数２） 第３の従来技術（Ｙ．Ｙｏｓｈｉｄａ ｅｔ ａｌ．、
Ｔｅｎｔｈ Ｉｎｔｅｒｎａｔｉｏｎａｌ Ｃｏｎｆｅ
ｒｅｎｃｅ ｏｎ Ｂｉｏｍａｇｎｅｔｉｓｍ、Ｓａｎ
ｔａｎａ Ｆｅ、Ｎｅｗ Ｍｅｘｉｃｏ、Ｆｅｂ．２０
（１９９６））では、コイルの面がそれぞれ直交した３
つの検出コイルからなるベクトル磁場センサを用いて生
体磁場の法線成分と２つの接線成分を検出し、磁場成分
の検出結果を直交座標系に変換して、直交座標系の成分
Ｂ_x、Ｂ_y、Ｂ_zを求め、法線成分Ｂ_z及び２次元ベクトル
強度Ｂ_xyの等磁場線図をそれぞれ表示している。| B _xy (x, y, t) | = {(B _x (x, y, t)) ² + (B _y (x, y, t)) ² } (Equation 2) Prior art (Y. Yoshida et al.,
Tenth International Confes
rence on Biomagnetism, San
tana Fe, New Mexico, Feb. 20
In (1996)), three planes of the coil were orthogonal to each other.
One of using vector magnetic field sensor comprising a detecting coil to detect the normal component and two tangential components of the biomagnetic field, converts the detection result of the magnetic field components in orthogonal coordinate system, an orthogonal coordinate system of the component B _x, B _y , B _z are obtained, and isomagnetic field diagrams of the normal component B _z and the two-dimensional vector intensity B _xy are displayed, respectively.

【０００９】第４の従来技術（Ｋ．Ｔｓｕｋａｄａ、ｅ
ｔ．ａｌ．、Ｔｅｎｔｈ Ｉｎｔｅｒｎａｔｉｏｎａｌ
Ｃｏｎｆｅｒｅｎｃｅ ｏｎ Ｂｉｏｍａｇｎｅｔｉ
ｓｍ、Ｓａｎｔａｎａ Ｆｅ、Ｎｅｗ Ｍｅｘｉｃｏ、
Ｆｅｂ．２０（１９９６））では、生体磁場の２つの接
線成分Ｂ_x、Ｂ_yを検出し、｜Ｂ_xy｜＝｜Ｂ_x＋Ｂ_y｜に基
づく等磁場線図と法線成分Ｂ_zに基づく等磁場線図との
比較を行なっている。A fourth prior art (K. Tsukada, e.
t. al. , Tenth International
Conference on Biomagneti
sm, Santana Fe, New Mexico,
Feb. 20 (1996)), the two tangential components B _x biomagnetic field to detect _{B y, | B xy | =} | B x + B y | in based like magnetic field diagrams and the like magnetic field based on the normal component B _z Comparison with the diagram is performed.

【００１０】生体内の電気的生理学現象の計測結果を表
す図として、脳波計により計測して得る脳波図（ＭＥ
Ｇ、magnetoencephalogram）、心電計により計測して得
る心電図（ＥＣＧ、electrocardiogram）がある。心電
図の計測に於いて、複数の電極を用いて心電図形をマッ
ピングする体表面心電図(body surface potential map)
は周知の技術である。これらの脳波図、又は体表面心電
図は、等しい電位点を結ぶ等電位線図として表示されて
いた。An electroencephalogram (ME) obtained by measuring with an electroencephalograph as a diagram showing the measurement result of an electrophysiological phenomenon in a living body.
G, magnetoencephalogram) and an electrocardiogram (ECG, electrocardiogram) obtained by measuring with an electrocardiograph. A body surface potential map that maps ECGs using multiple electrodes in ECG measurement
Is a well-known technique. These electroencephalograms or body surface electrocardiograms were displayed as equipotential maps connecting equal potential points.

【００１１】第５の従来技術（Ｔ．Ｊ．Ｍｏｎｔａｇｕ
ｅ ｅｔ ａｌ．、Ｃｉｒｃｕｌａｔｉｏｎ ６３、Ｎ
ｏ．５、ｐｐ１１６６−１１７２（１９８１））では、
複数の電極の各電極の出力の時間変化波形を任意の時間
区間で積分した等積分図(isointegral map)を、体表面
心電図として表示している。A fifth prior art (TJ Montagu)
e et al. , Circulation 63, N
o. 5, pp1166-1172 (1981))
An isointegral map obtained by integrating the time-varying waveform of the output of each of the plurality of electrodes in an arbitrary time interval is displayed as a body surface electrocardiogram.

【００１２】[0012]

【発明が解決しようとする課題】以下の説明では、「生
体磁場」は「生体磁場から発する磁場」を意味し、「心
磁場計測」は、「心臓から発する磁場の計測」を意味
し、「心磁波形」は、「心磁場計測により得た心磁図
（ＭＣＧ、Magnetocardiogram）が表わす波形」を意味
する。また、「脳磁場計測」は、「脳から発する磁場の
計測」を意味し、「脳磁波形」は、「脳磁場計測により
得た脳磁図（ＭＥＧ、Magnetoencephalogram）が表わす
波形」を意味する。In the following description, "biomagnetic field" means "magnetic field generated from biomagnetic field", "cardiac magnetic field measurement" means "measurement of magnetic field generated from heart", The “magnetocardiogram waveform” means “a waveform represented by a magnetocardiogram (MCG, Magnetocardiogram) obtained by measuring a magnetocardiogram”. “Measurement of a brain magnetic field” means “measurement of a magnetic field generated from the brain”, and “magnetoencephalogram waveform” means “a waveform represented by a magnetoencephalogram (MEG) obtained by measurement of a brain magnetic field”.

【００１３】従来技術に於ける各成分毎の等磁場線図は
それぞれ特徴があり、単一電流ダイポールが存在する時
には、法線成分Ｂ_zの等磁場線図では、電流源の位置、
大きさ、方向等が容易に解析できる。一方、接線成分Ｂ
_x、Ｂ_yの計測結果から得る２次元ベクトル強度Ｂ_xyの等
磁場線図では、複数の電流ダイポールが存在する時で
も、容易に各電流ダイポールを判別できる特徴がある。
しかし、磁場を検出するコイルの数はｘ、ｙ方向それぞ
れに必要であるため、法線成分Ｂ_zのみの検出に比べて
コイル数が２倍になる。また、Ｂ_x、Ｂ_y、Ｂ_zの全ての
成分を計測するベクトル計測では、法線成分Ｂ_zのみの
検出に比べて３倍の数のコイルが必要となる。このた
め、検出コイルとＳＱＵＩＤからなる磁場センサの数は
増加し、更に、信号処理回路等も増加し、生体磁場計測
システムは高価なシステムとなってしまうという問題が
あった。また、第１の従来技術では、各計測点上にアロ
ーを表示するだけであり、電流源の詳細な分布状態が識
別しにくいという問題があった。[0013] isomagnetic field for each in the prior art component is characterized respectively, when a single current dipole exists, the isomagnetic field of the normal component B _z, of the current source position,
Size, direction, etc. can be easily analyzed. On the other hand, the tangent component B
_x, in isomagnetic field of a two-dimensional vector magnitude B _xy obtained from the measurement results of B _y, even when the plurality of current dipoles are present, is characterized readily determine each current dipole.
However, the number of coils for detecting a magnetic field x, because it is required in the y-direction, respectively, the number of coils is doubled as compared with the detection of the normal component B _z only. Also, B _x, B _y, the vector measurement for measuring all the components B _z, it is necessary to 3 times the number of coils in comparison to the detection of the normal component B _z only. For this reason, the number of magnetic field sensors including the detection coil and the SQUID increases, and further, the number of signal processing circuits and the like increases, and there is a problem that the biomagnetic field measurement system becomes an expensive system. Further, in the first conventional technique, there is a problem that only an arrow is displayed on each measurement point, and it is difficult to identify a detailed distribution state of the current source.

【００１４】生体磁場成分で表わした等磁場線図によ
り、任意の時点での生体内の電流源の位置、大きさ、方
向等を解析でき詳細な電流源の位置、大きさ、方向等の
情報の変化を知ることができる。従来技術では、装置に
表示、又は出力された多数の図を用いて各種情報のダイ
ナミックな変化をとらえ疾患等の診断を行っていた。し
かし、従来技術では、診断のために各種情報を表す多数
の図を必要とし、各種情報の変化の異常を経験的に行っ
ていた。この様に従来技術では、どの生体部位でどのよ
うな大きさの電流が流れたか、又は異常な生体電流が流
れている領域がどこであるか等を表わす総合的な情報を
１つの図として表示するための処理は実行されてなかっ
た。また、体表面心電図では、任意の時間間隔（Ｑ、
Ｒ、Ｓの各波の発生する期間、Ｓ波からＴ波の発生する
期間等）での積分値の等しい点を示す等積分図では、連
続する各時刻での等電位線図を複数必要とせず、１つの
心電図形で心臓の情報を得ることができる。しかし、等
電位線図では心臓内の電流源を１つの電流ダイポールと
仮定しておくと、電流ダイポールの直上ではなく電流ダ
イポールの直上から離れた位置に陽極のピークと陰極の
ピークが存在する図形となってしまうという問題があ
る。更に、電流ダイポールの位置が変化せず電流ダイポ
ールの方向が変化すると陽極及び陰極のピーク位置が変
化してしまい、電位を積分する時に電流源と積分値のピ
ークとが対応しなくなるという問題があった。また、生
体磁場計測により得る生体磁場の成分を単に積分して
も、心電図の場合と同様に、生体磁場成分のピーク位置
と電流源の位置が対応しないという問題があった。ま
た、心電図から得る等積分図のみでは、臓器の位置、大
きさ等の個人差があり単純に等積分図から疾患等の異常
を正確に判断することが困難であるという問題があっ
た。The position, size, direction, etc. of the current source in the living body at any time can be analyzed from the isomagnetic field diagram represented by the biomagnetic field component, and detailed information on the position, size, direction, etc. of the current source can be obtained. You can know the change. In the related art, diagnosis of a disease or the like is performed by capturing dynamic changes in various types of information using a large number of figures displayed or output on an apparatus. However, in the related art, a large number of diagrams representing various types of information are required for diagnosis, and abnormal changes in various types of information are empirically performed. As described above, in the related art, comprehensive information indicating what magnitude of current has flowed in which biological part, or where an abnormal biological current is flowing, is displayed as one diagram. Has not been executed. In the body surface electrocardiogram, an arbitrary time interval (Q,
In the isointegral diagram showing the points where the integrated values are equal in the period in which each of the R and S waves is generated, the period in which the S wave is generated in the T wave, and the like, a plurality of equipotential diagrams at successive times are required. Instead, heart information can be obtained with one electrocardiogram. However, assuming that the current source in the heart is a single current dipole in the equipotential diagram, a figure in which the anode peak and the cathode peak exist not directly above the current dipole but at a position away from directly above the current dipole There is a problem that it becomes. Further, if the position of the current dipole does not change and the direction of the current dipole changes, the peak positions of the anode and the cathode change, and there is a problem that the peak of the current source does not correspond to the peak of the integrated value when integrating the potential. Was. Further, even if the components of the biomagnetic field obtained by the biomagnetic field measurement are simply integrated, there is a problem that the peak position of the biomagnetic component does not correspond to the position of the current source as in the case of the electrocardiogram. In addition, there is a problem that it is difficult to accurately determine an abnormality such as a disease simply from the isointegral diagram due to individual differences in the position, size, etc. of the organs using only the isointegral diagram obtained from the electrocardiogram.

【００１５】本発明の目的は、従来技術で必要としてい
た図（マップ）の数よりもはるかに少数の図（マップ）
を用いて、生体部位の全体の状態を把握できる生体磁場
計測方法及び生体磁場計測装置を提供することにある。It is an object of the present invention to provide a much smaller number of maps (maps) than required in the prior art.
The present invention provides a biomagnetic field measurement method and a biomagnetic field measurement device that can grasp the entire state of a living body part using the method.

【００１６】本発明の他の目的は、検出コイルの数を増
加させることなく、生体磁場の垂直成分Ｂ_zを計測して
磁場源の解析を可能とする生体磁場計測方法及び生体磁
場計測装置を提供することにある。Another object of the present invention is to provide a biomagnetic field measuring method and a biomagnetic field measuring apparatus which can measure a vertical component _Bz of a biomagnetic field and analyze a magnetic field source without increasing the number of detection coils. To provide.

【００１７】[0017]

【課題を解決するための手段】本発明の生体磁場計測方
法では、（１）量子干渉素子（ＳＱＵＩＤ）からなり、
生体の外部に配置される複数の磁束計を用いて、生体か
ら発する生体磁場の生体の面に垂直な第１方向の磁場成
分の時間変化を計測する第１の工程と、第１方向と交叉
する第２方向及び第３方向に於ける第１方向の磁場成分
の変化率の２乗和の平方根に比例する値の時間変化を求
める第２の工程と、この第２の工程で得る値の時間変化
を所定の期間で積分し積分値を求める第３の工程と、こ
の第３の工程の工程で得る積分値を表示する第４の工程
とを有することに特徴があり、更に、（２）量子干渉素
子（ＳＱＵＩＤ）からなり、生体の外部に配置される複
数の磁束計を用いて、生体から発する生体磁場の生体の
面に平行な第１、第２方向の磁場成分の時間変化を計測
する第１の工程と、第１、第２方向の磁場成分の２乗和
の平方根に比例する値の時間変化を求める第２の工程
と、この第２の工程で得る値の時間変化を所定の期間で
積分し積分値を求める第３の工程と、この第３の工程の
工程で得る積分値を表示する第４の工程とを有すること
に特徴がある。また上記（１）、（２）の特徴を有する
生体磁場計測方法に於いて、上記の積分値を用いて、内
挿、外挿により、上記の第４の工程で積分値が等しい点
を結ぶ等積分図を表示すること、上記の第３の工程に於
いて、上記の第２の工程で得る上記の値の時間変化を所
定の期間で積分し積分値を求めることを、複数の所定の
期間で行ない積分値を複数個求め、この複数個の積分値
の間での、比、等加重を含む和又は差の何れかを求める
演算を行なうことにも特徴がある。なお、直交座標系
（ｘ、ｙ、ｚ）に於いて、生体表面に垂直な方向をｚ軸
とし、第１方向をｚ方向、第２方向をｘ方向、第３方向
をｙ方向とする。また、極座標系（ｒ、θ、φ）におい
て、生体表面に垂直な方向をｒ軸とし、第１方向をｒ方
向、第２方向をθ方向、第３方向をφ方向とする。According to the biomagnetic field measuring method of the present invention, (1) a quantum interference device (SQUID)
A first step of measuring a temporal change of a magnetic field component of a biomagnetic field emitted from the living body in a first direction perpendicular to the surface of the living body, using a plurality of magnetometers arranged outside the living body, and intersecting the first direction A second step of obtaining a time change of a value proportional to the square root of the sum of squares of the change rate of the magnetic field component in the first direction in the second direction and the third direction; It is characterized in that it has a third step of integrating the time change over a predetermined period to obtain an integrated value, and a fourth step of displaying the integrated value obtained in the third step. ) Using a plurality of magnetometers, which are composed of quantum interference devices (SQUIDs) and are arranged outside the living body, the time change of the magnetic field components of the biomagnetic field emitted from the living body in the first and second directions parallel to the surface of the living body is measured. The first step of measurement and the square root of the sum of squares of the magnetic field components in the first and second directions A second step of obtaining a time change of the value, a third step of integrating the time change of the value obtained in the second step in a predetermined period to obtain an integrated value, and an integration obtained in the step of the third step And a fourth step of displaying a value. Further, in the biomagnetic field measuring method having the features of (1) and (2), points having the same integral value in the fourth step are connected by interpolation and extrapolation using the integral value. Displaying an isointegral diagram and, in the third step, integrating a time change of the value obtained in the second step in a predetermined period to obtain an integral value are performed by a plurality of predetermined operations. It is also characterized in that a plurality of integrated values are calculated in a period, and a calculation is performed to obtain any of a sum or a difference including ratios and equal weights among the plurality of integrated values. In the orthogonal coordinate system (x, y, z), the direction perpendicular to the surface of the living body is set as the z-axis, the first direction is set as the z direction, the second direction is set as the x direction, and the third direction is set as the y direction. In the polar coordinate system (r, θ, φ), a direction perpendicular to the surface of the living body is defined as an r axis, a first direction is defined as an r direction, a second direction is defined as a θ direction, and a third direction is defined as a φ direction.

【００１８】本発明の生体磁場計測装置では、（１）量
子干渉素子（ＳＱＵＩＤ）からなり生体から発する生体
磁場を信号として検出する、生体の外部に配置される複
数の磁束計と、信号の演算処理を行なう演算処理手段
と、演算処理結果を表示する表示手段とを有し、生体磁
場分布を計測する生体磁場計測装置に於いて、磁束計
は、生体磁場の生体の面に垂直な第１方向の磁場成分の
時間変化を検出し、演算処理手段は、第１方向と交叉す
る第２方向及び第３方向に於ける第１方向の磁場成分の
変化率の２乗和の平方根に比例する値の時間変化を求め
る演算と、この値の時間変化を所定の期間で積分し積分
値を求める演算とを行ない、表示手段に積分値を表示す
ることに特徴があり、更に、（２）同上の生体磁場計測
装置に於いて、磁束計は、生体磁場の生体の面に平行な
第１、第２方向の磁場成分の時間変化を検出し、演算処
理手段は、第１、第２方向の磁場成分の２乗和の平方根
に比例する値の時間変化を求める演算と、この値の時間
変化を所定の期間で積分し積分値を求める演算とを行な
い、表示手段に積分値を表示することに特徴がある。ま
た、上記（１）、（２）の特徴を有する生体磁場計測装
置に於いて、表示手段に、内挿、外挿により積分値の等
しい点を結ぶ等積分図が表示されること、演算処理手段
は、上記値の時間変化を所定の期間で積分し積分値を求
めることを、複数の所定の期間で行ない積分値を複数個
求め、この複数個の積分値の間での、比、等加重を含む
和又は差の何れかを求める演算を行なうこと、複数の磁
束計が、生体の面に等間隔に配置されることにも特徴が
ある。本発明の生体磁場計測装置では、心臓から発する
磁場の、胸面に対する法線（垂直）成分、接線（平行）
成分の同時表示が可能である。なお、直交座標系（ｘ、
ｙ、ｚ）に於いて、生体表面に垂直な方向をｚ軸とし、
第１方向をｚ方向、第２方向をｘ方向、第３方向をｙ方
向とする。また、極座標系（ｒ、θ、φ）において、生
体表面に垂直な方向をｒ軸とし、第１方向をｒ方向、第
２方向をθ方向、第３方向をφ方向とする。According to the biomagnetic field measuring apparatus of the present invention, (1) a plurality of magnetometers, which are composed of quantum interference elements (SQUIDs) and detect a biomagnetic field emitted from a living body as a signal, and are arranged outside the living body, In a biomagnetic field measuring apparatus for measuring a biomagnetic field distribution, the magnetometer has a first magnetic field perpendicular to the surface of the living body of the biomagnetic field. The temporal change of the magnetic field component in the direction is detected, and the arithmetic processing means is proportional to the square root of the sum of squares of the rate of change of the magnetic field component in the first direction in the second direction and the third direction crossing the first direction. It is characterized in that a calculation for obtaining a time change of the value and a calculation for integrating the time change of the value over a predetermined period to obtain an integrated value are performed, and the integrated value is displayed on the display means. Magnetic flux meter Detecting the time change of the magnetic field components of the biomagnetic field in the first and second directions parallel to the surface of the living body, and calculating a value proportional to the square root of the sum of squares of the magnetic field components in the first and second directions. And a calculation for obtaining an integrated value by integrating the time change of this value over a predetermined period, and displaying the integrated value on the display means. Further, in the biomagnetic field measuring apparatus having the features of (1) and (2), the display means displays an isointegral diagram connecting points having equal integral values by interpolation and extrapolation. The means integrates a time change of the value over a predetermined period to obtain an integral value in a plurality of predetermined periods, obtains a plurality of integral values, and calculates a ratio between the plurality of integral values. It is also characterized in that the calculation for obtaining either the sum or the difference including the weight is performed, and that a plurality of magnetometers are arranged at equal intervals on the surface of the living body. In the biomagnetic field measuring apparatus of the present invention, the normal (vertical) component and the tangent (parallel) of the magnetic field generated from the heart with respect to the chest surface
Simultaneous display of components is possible. Note that a rectangular coordinate system (x,
y, z), the direction perpendicular to the living body surface is defined as the z axis,
The first direction is the z direction, the second direction is the x direction, and the third direction is the y direction. In the polar coordinate system (r, θ, φ), a direction perpendicular to the surface of the living body is defined as an r axis, a first direction is defined as an r direction, a second direction is defined as a θ direction, and a third direction is defined as a φ direction.

【００１９】本発明の本質的な特徴は、生体表面に垂直
な方向を直交座標（ｘ、ｙ、ｚ）のｚ軸とし、生体表面
に平行な面を（ｘ、ｙ）平面とする時、生体磁場の生体
表面に垂直な法線成分Ｂ_z(ｘ、ｙ）を検出し、生体磁場
の生体表面に平行な接線成分Ｂ_x、Ｂ_yをそれぞれ、法線
成分Ｂ_zのｘ方向、ｙ方向に於ける変化率から推定する
ことに特徴がある。The essential feature of the present invention is that when the direction perpendicular to the surface of the living body is the z-axis of the rectangular coordinates (x, y, z) and the plane parallel to the surface of the living body is the (x, y) plane, perpendicular normal component to the living body surface of the biomagnetic field B _z (x, y) is detected and tangential component parallel to the biological surface of the biomagnetic field B _x, B _y, respectively, x-direction of the normal component B _z, y The feature is that it is estimated from the rate of change in the direction.

【００２０】本発明によれば、接線成分Ｂ_x、Ｂ_yを測定
する検出コイルを必要とすることなく、生体の磁場分布
を２次元（ｘ、ｙ）平面に投影した等磁場線図を得るこ
とができ、等磁場線図のピークパターンから生体内の電
流源を判別でき、複数の電流ダイポールの（ｘ、ｙ）座
標での位置を知ることができる。According to the invention, obtained tangential components B _x, without the need for a detection coil for measuring the B _y, the magnetic field distribution of the living body 2D (x, y) the isomagnetic field projected on a plane The current source in the living body can be determined from the peak pattern of the isomagnetic field map, and the positions of the plurality of current dipoles at the (x, y) coordinates can be known.

【００２１】以下、本発明に於ける演算処理手段（複数
の磁束計により計測された信号を収集し、信号に対して
以下の演算処理を行なうパソコン等の計算機、又は専用
的にハードウエア化され演算処理を行なう電子回路）に
て行なう演算処理の内容に付いて説明する。Hereinafter, the arithmetic processing means of the present invention (a computer such as a personal computer which collects signals measured by a plurality of magnetometers and performs the following arithmetic processing on the signals, or a hardware dedicated to the arithmetic processing means) The contents of the arithmetic processing performed by the electronic circuit for performing the arithmetic processing) will be described.

【００２２】量子干渉素子（ＳＱＵＩＤ）からなる複数
の磁束計を用いて、生体表面の位置（ｘ、ｙ）に於いて
生体から発する磁場の接線成分（生体の面に平行な成
分）Ｂ_x（ｘ、ｙ、ｔ）、Ｂ_y（ｘ、ｙ、ｔ）を計測する
場合には（但し、直交座標系（ｘ、ｙ、ｚ）に於いて生
体の面に平行な面をｘｙ面、生体の面に垂直な軸をｚと
する）、接線成分Ｂ_x（ｘ、ｙ、ｔ）とＢ_y（ｘ、ｙ、
ｔ）の２乗和の平方根から２次元ベクトル強度│Ｂ
_xy（ｘ、ｙ）│（以下、│ │は絶対値を表わす）を
（数３）により求める。Using a plurality of magnetometers composed of quantum interference devices (SQUIDs), a tangential component (a component parallel to the surface of the living body) of a magnetic field emitted from the living body at a position (x, y) on the surface of the living body, B _x ( x, y, t), B y (x, y, in the case of measuring the t) (where the orthogonal coordinate system (x, y, z) xy plane parallel to the plane of the living body at the biological of the z axis perpendicular to the plane), the tangential component _{B x (x, y, t} ) and B _y (x, y,
t) From the square root of the sum of squares of two-dimensional vector intensity | B
_xy (x, y) | (hereinafter, | | represents an absolute value) is obtained by (Equation 3).

【００２３】 │Ｂ_xy（ｘ、ｙ、ｔ）│＝√｛（Ｂ_x（ｘ、ｙ、ｔ））² ＋（Ｂ_y（ｘ、ｙ、ｔ））²｝ …（数３） 次いで、各点（ｘ、ｙ）について任意の期間での波形│
Ｂ_xy（ｘ、ｙ、ｔ）│の積分値Ｉ₁（ｘ、ｙ）を（数
４）により求め、内挿、外挿により各点（ｘ、ｙ）での
積分値Ｉ₁（ｘ、ｙ）が同じ値の点を結ぶ等積分図を求
めて、等積分図を表示画面に表示する。| B _xy (x, y, t) | = {(B _x (x, y, t)) ² + (B _y (x, y, t)) ² } (Equation 3) Waveform at an arbitrary period for each point (x, y) |
_{B xy (x, y, t} ) │ integrated value I ₁ (x, y) and calculated by equation (4), the interpolation, the integral value I ₁ at each point by extrapolation (x, y) (x, y) obtains an isointegral diagram connecting points having the same value, and displays the isointegral diagram on the display screen.

【００２４】 Ｉ₁（ｘ、ｙ）＝∫│Ｂ_xy（ｘ、ｙ、ｔ）│ｄｔ …（数４） 以下、計測された生体の面に垂直な磁場成分Ｂ_z（ｘ、
ｙ、ｔ）（法線成分）から、接線成分Ｂ_x、Ｂ_yを推定す
ること説明する。I ₁ (x, y) = ∫│B _xy (x, y, t) │dt (Equation 4) Hereinafter, the measured magnetic field component B _z (x,
The estimation of the tangent components B _x and B _y from y, t) (normal component) will be described.

【００２５】生体磁場の体表面に平行な接線成分は、体
表面直下を流れる電流を最もよく反映していることを利
用すると、電流の流れる向きと磁場の向きの関係から、
計測された磁場の接線ベクトル（Ｂ_x、Ｂ_y）を反時計回
りに９０°回転させることにより、生体内の電流分布を
生体表面に平行な２次元平面に投影して概観できる。即
ち、〈ｅ_x〉、〈ｅ_y〉をそれぞれｘ軸方向、ｙ軸方向の
単位ベクトルとして、各計測点に於ける接線成分Ｂ_x、
Ｂ_yから、（数５）に示す電流ベクトク〈Ｊ〉を求め、
各計測点（ｘ、ｙ）に於ける電流ベクトル場の分布（ア
ローマップ）として表現することができる。Taking advantage of the fact that the tangential component of the biomagnetic field parallel to the body surface best reflects the current flowing just below the body surface, the relationship between the current flowing direction and the magnetic field direction is given by
By rotating the tangent vectors (B _x , B _y ) of the measured magnetic field by 90 ° counterclockwise, the current distribution in the living body can be projected on a two-dimensional plane parallel to the surface of the living body to give an overview. _{That, <e x>, <e} y> the x-axis direction, a unit vector in the y-axis direction, in the tangential component B _x in the respective measurement points,
From B _y, seek Bekutoku <J> current shown in (5),
It can be expressed as a distribution (arrow map) of a current vector field at each measurement point (x, y).

【００２６】 〈Ｊ〉＝−Ｂ_y〈ｅ_x〉＋Ｂ_x〈ｅ_y〉 …（数５） 一方、磁場の生体表面に垂直な法線成分Ｂ_zを計測する
場合、（数１）により表現される電流ベクトルを用いた
アローマップが定義されている（第１の従来技術：Ｈ、
Ｈｏｓａｋａ ａｎｄ Ｄ．Ｃｏｈｅｎ（１９７
６））。[0026] _{<J> = - B y <} e x> + B x <e y> ... ( 5) On the other hand, when measuring the perpendicular normal component B _z to a biological surface of the magnetic field, expressed by equation (1) An arrow map using a current vector is defined (first conventional technique: H,
Hosaka and D.S. Cohen (197
6)).

【００２７】 〈Ｊ〉＝（∂Ｂ_z／∂ｙ）〈ｅ_x〉−（∂Ｂ_z／∂ｘ）〈ｅ_y〉 …（数１） 本願発明の発明者らは、（数１）と（数５）との比較か
ら、（数６）及び（数７）が成立する可能性、即ち、計
測された磁場の法線成分Ｂ_zから接線成分Ｂ_x及びＢ_yを
導出できる可能性があることを見い出し、種々の検討を
行なった。以下、検討の結果を詳細に説明する。[0027] _{<J> = (∂B z /} ∂y) <e x> - (∂B z / ∂x) <e y> ... ( Equation 1) The inventors of the present invention, the (number 1) comparison with equation (5), (6) and possibly (7) is satisfied, that is, the possibility of deriving the tangential components B _x and B _y of the normal component B _z of the measured magnetic field We found something and made various studies. Hereinafter, the result of the study will be described in detail.

【００２８】 Ｂ_x＝−（∂Ｂ_z／∂ｘ） …（数６） Ｂ_y＝−（∂Ｂ_z／∂ｙ） …（数７） 図１は、心臓の活動による磁場（心磁場）の発生を、無
限平面導体中の電流ダイポールから発生する磁場により
モデル化して解析するための図である。図１に於いて、
Ｐは直交座標系（ｘ、ｙ、ｚ）のｘｙ面に表面を持つ無
限平面導体、〈Ｑ〉は位置ベクトル〈ｒ₀（ｘ₀、ｙ₀、
ｚ₀）〉で示される位置に存在する電流ダイポールのモ
ーメント、〈ｒ（ｘ、ｙ、ｚ）〉は磁束密度（磁場）
〈Ｂ（ｒ）〉を計測する計測点の位置ベクトルを示す。
図１に示すモデルに於いて、無限平面導体Ｐの外部に生
じる磁場〈Ｂ（ｒ）〉は、Ｓａｒｖａｓ（文献：Ｐｈｙ
ｓ．Ｍｅｄ．Ｂｉｏｌ．、Ｖｏｌ．３２、Ｎｏ．１、１
１−２２（１９８７））により定式化されており、（数
８）により表現される。[0028] _{B x = - (∂B z /} ∂x) ... ( number _{6) B y = - (∂B} z / ∂y) ... ( 7) FIG. 1, the magnetic field due to activity of the heart (cardiac magnetic field) FIG. 5 is a diagram for modeling and analyzing the occurrence of に よ り with a magnetic field generated from a current dipole in an infinite plane conductor. In FIG.
P is an infinite plane conductor having a surface on the xy plane of a rectangular coordinate system (x, y, z), and <Q> is a position vector <r ₀ (x ₀ , y ₀ ,
z ₀ )> the moment of the current dipole present at the position indicated by <r (x, y, z)> is the magnetic flux density (magnetic field)
It shows the position vector of the measurement point for measuring <B (r)>.
In the model shown in FIG. 1, the magnetic field <B (r)> generated outside the infinite planar conductor P is Sarvas (reference: Phy).
s. Med. Biol. Vol. 32, no. 1, 1
1-22 (1987)), and is expressed by (Equation 8).

【００２９】 〈Ｂ（ｒ）〉＝｛μ₀／（４πＫ²）｝｛〈Ｑ〉×〈ａ〉・〈ｅ_z〉∇Ｋ −Ｋ〈ｅ_z〉×〈Ｑ〉｝ …（数８） （数８）に於いて、μ₀は真空の透磁率、〈ｅ_z〉はｚ軸
方向の単位ベクトル、×はベクトル積、・はスカラ積、
∇はｇｒａｄ＝（∂／∂ｘ、∂／∂ｙ、∂／∂ｚ）を表
わし、〈ａ〉は（数９）、ａは（数１０）、Ｋは（数１
１）、∇Ｋは（数１２）により示される。｜ ｜は絶対
値を示す。[0029] <B (r)> = { μ 0 / (4πK 2)} {<Q> × <a> · <e z> ∇K -K <e z> × <Q>} ... ( 8) In (Equation 8), μ ₀ is the magnetic permeability of vacuum, <e _z > is a unit vector in the z-axis direction, × is a vector product, • is a scalar product,
∇ represents grad = (∂ / ∂x, ∂ / ∂y, ∂ / ∂z), <a> is (Formula 9), a is (Formula 10), and K is (Formula 1)
1) and ΔK are represented by (Equation 12). || indicates an absolute value.

【００３０】 〈ａ〉＝〈ｒ（ｘ、ｙ、ｚ）〉−〈ｒ₀（ｘ₀、ｙ₀、ｚ₀）〉 …（数９） ａ＝｜〈ａ〉｜ …（数１０） Ｋ＝ａ（ａ＋〈ａ〉・〈ｅ_z〉） …（数１１） ∇Ｋ＝（２＋ａ^-1〈ａ〉・〈ｅ_z〉）〈ａ〉＋ａ〈ｅ_z〉 …（数１２） （数８）により示される〈Ｂ〉（ｒ）の無限平面導体Ｐ
に平行な接線成分Ｂ_x及びＢ_yと、無限平面導体Ｐに垂直
なな法線成分Ｂ_zは、それぞれ（数１３）、（数１
４）、（数１５）により与えられる。<a> = <r (x, y, z)> − <r ₀ (x ₀ , y ₀ , z ₀ )> (Equation 9) a = | <a> | (Equation 10) K = a (a + <a> · < e z>) ... ( number 11) ∇K = (2 + a -1 <a> · <e z>) <a> + a <e z> ... ( Equation 12) (Equation 8 )), The infinite planar conductor P of <B> (r)
Parallel tangential components B _x and B _y, the infinite plane conductor perpendicular Do normal component B _z to P, respectively (equation 13), (Equation 1
4), given by (Equation 15).

【００３１】 Ｂ_x＝｛μ₀／（４πＫ²）｝ ×［｛Ｑ_x（ｙ−ｙ₀）−Ｑ_y（ｘ−ｘ₀）｝（∇Ｋ）_x＋ＫＱ_y］…（数１３） Ｂ_y＝｛μ₀／（４πＫ²）｝ ×［｛Ｑ_y（ｙ−ｙ₀）−Ｑ_x（ｘ−ｘ₀）｝（∇Ｋ）_y＋ＫＱ_x］…（数１４） Ｂ_Z＝｛μ₀／（４πＫ²）｝ ×［｛Ｑ_x（ｙ−ｙ₀）−Ｑ_y（ｘ−ｘ₀）｝（∇Ｋ）_z］ …（数１５） 一方、（数１３）により示される法線成分Ｂ_Zのｘ方向
に於ける微分は（数１６）により表わされる。B _x = {μ ₀ / (4πK ² )} × [{Q _x (y−y ₀ ) −Q _y (x−x ₀ )} (∇K) _x + KQ _y ] (Equation 13) B _y = {μ ₀ / (4πK ² )} × [{Q _y (y−y ₀ ) −Q _x (x−x ₀ )} (∇K) _y + KQ _x ] (Equation 14) B _Z = ｛μ ₀ / (4πK ² )｝ × [｛Q _x (y-y ₀ ) -Q _y (x-x ₀ )｝ (∇K) _z ] (Equation 15) On the other hand, the normal represented by (Equation 13) The derivative of the component _{BZ in} the x direction is represented by (Equation 16).

【００３２】 ∂Ｂ_Z／∂ｘ＝｛μ₀／（４πＫ²）｝×［｛Ｑ_x（ｙ−ｙ₀）−Ｑ_y（ｘ−ｘ₀）｝ ｛−２（∇Ｋ）_z（∇Ｋ）_x／Ｋ−ａ^-3（ｘ−ｘ₀）（ｚ−ｚ₀）²＋ａ^-1（ｘ−ｘ₀ ）｝−（∇Ｋ）_zＱ_y］ …（数１６） 同様に、法線成分Ｂ_Zのｙ方向に於ける微分は（数１
７）により表わされる。{B _Z / ∂x = {μ ₀ / (4πK ² )} × [{Q _x (y−y ₀ ) −Q _y (x−x ₀ )}} − 2 (∇K) _z (∇ K) _x / K−a ⁻³ (x−x ₀ ) (z−z ₀ ) ² + a ⁻¹ (x−x ₀ )｝ − (∇K) _z Q _y ] (Equation 16) The derivative of the line component B _{Z in} the y direction is (Equation 1)
7).

【００３３】 ∂Ｂ_Z／∂ｙ＝−｛μ₀／（４πＫ²）｝×［｛Ｑ_x（ｙ−ｙ₀）−Ｑ_y（ｘ−ｘ₀） ｝｛２（∇Ｋ）_z（∇Ｋ）_y／Ｋ＋ａ^-3（ｙ−ｙ₀）（ｚ−ｚ₀）²−ａ^-1（ｙ−ｙ₀ ）｝＋（∇Ｋ）_zＱ_x］ （数１７） （数１６）、（数１７）に於いて、 α＝（∇Ｋ）_z／Ｋ …（数１８） β_x＝−ａ^-3（ｘ−ｘ₀）（ｚ−ｚ₀）²＋ａ^-1（ｘ−ｘ₀） …（数１９） β_y＝−ａ^-3（ｙ−ｙ₀）（ｚ−ｚ₀）²＋ａ^-1（ｙ−ｙ₀） …（数２０） と置く時、（数１６）、（数１７）はそれぞれ（数２
１）、（数２２）により表わされる。{B _Z / ∂y = − ｛μ ₀ / (4πK ² )} × [｛Q _x (y−y ₀ ) −Q _y (x−x ₀ )｝ ｛2 (∇K) _z (∇ K) _y / K + a ⁻³ (y−y ₀ ) (z−z ₀ ) ² −a ⁻¹ (y−y ₀ )｝ + (∇K) _z Q _x ] (Equation 17) (Equation 16), (Equation 16) In Equation 17, α = (∇K) _z / K (Equation 18) β _x = −a ⁻³ (x−x ₀ ) (z−z ₀ ) ² + a ⁻¹ (x−x ₀ ) ... (Equation 19) β _y = −a ⁻³ (y−y ₀ ) (z−z ₀ ) ² + a ⁻¹ (y−y ₀ ) (Equation 20) 17) is (Equation 2)
1) and (Expression 22).

【００３４】 ∂Ｂ_Z／∂ｘ＝−｛μ₀／（４πＫ²）｝×［｛Ｑ_x（ｙ−ｙ₀）−Ｑ_y（ｘ−ｘ₀） ｝｛２α（∇Ｋ）_x−β_x｝＋αＫＱ_y］ …（数２１） ∂Ｂ_Z／∂ｙ＝−｛μ₀／（４πＫ²）｝×［｛Ｑ_x（ｙ−ｙ₀）−Ｑ_y（ｘ−ｘ₀） ｝｛２α（∇Ｋ）_y−β_y｝＋αＫＱ_x］ …（数２２） 簡単のために、（数１３）、（数２１）、（数１４、
（数２２）を共通因子である｛μ₀／（４πＫ²）｝によ
り規格化して変形を行ない、（数１３’）、（数２
１’）、（数１４’）、（数２２’）を得る。∂B _Z / ∂x =-｛μ ₀ / (4πK ² )｝ × [｛Q _x (y-y ₀ ) -Q _y (x-x ₀ )｝ ｛2α (∇K) _x -β _x ｝ + αKQ _y ] (Equation 21) {B _Z / ∂y = − ｛μ ₀ / (4πK ² )} × [｛Q _x (y−y ₀ ) −Q _y (x−x ₀ )） 2α (∇K) _y −β _y ｝ + αKQ _x ] (Equation 22) For simplicity, (Equation 13), (Equation 21), (Equation 14)
(Equation 22) is normalized by {μ ₀ / (4πK ² )}, which is a common factor, and is transformed to obtain (Equation 13 ′) and (Equation ² ).
1 ′), (Equation 14 ′) and (Equation 22 ′) are obtained.

【００３５】 Ｂ_x＝（∇Ｋ）_x｛Ｑ_x（ｙ−ｙ₀）−Ｑ_y（ｘ−ｘ₀）｝＋ＫＱ_y …（数１３’） ∂Ｂ_Z／∂ｘ＝ −２α（∇Ｋ）_x｛Ｑ_x（ｙ−ｙ₀）−Ｑ_y（ｘ−ｘ₀）｝−αＫＱ_y ＋β_x｛Ｑ_x（ｙ−ｙ₀）−Ｑ_y（ｘ−ｘ₀）｝＝ −２αＢ_x＋αＫＱ_y＋β_x｛Ｑ_x（ｙ−ｙ₀）−Ｑ_y（ｘ−ｘ₀）｝ …（数２１’） Ｂ_y＝（∇Ｋ）_y｛Ｑ_y（ｙ−ｙ₀）−Ｑ_x（ｘ−ｘ₀）｝＋ＫＱ_x …（数１４’） ∂Ｂ_Z／∂ｙ＝ −２α（∇Ｋ）_y｛Ｑ_x（ｙ−ｙ₀）−Ｑ_y（ｘ−ｘ₀）｝−αＫＱ_x］ ＋β_y｛Ｑ_x（ｙ−ｙ₀）−Ｑ_y（ｘ−ｘ₀）｝＝ −２αＢ_y＋αＫＱ_x＋β_y｛Ｑ_x（ｙ−ｙ₀）−Ｑ_y（ｘ−ｘ₀）｝ …（数２２’） （数１３’）と（数２１’）とから明らかなように、∂
Ｂ_Z／∂ｘの値は、接線成分Ｂ_xの−２α倍に等しい項
に、２つの付加項を加算した値に等しく、（数１４’）
と（数２２’）とから明らかなように、∂Ｂ_Z／∂ｙの
値は、接線成分Ｂ_yの−２α倍に等しい項に、２つの付
加項を加算した値に等しい。B _x = (∇K) _x ｛Q _x (y−y ₀ ) −Q _y (x−x ₀ )} + KQ _y (Equation 13 ′) ∂B _Z / ∂x = −2α (∇K ) _X {Q _x (y-y ₀ ) -Q _y (x-x ₀ )}-αKQ _y + β _x {Q _x (y-y ₀ ) -Q _y (x-x ₀ )} = − 2αB _x + αKQ _{y + β x {Q x (} y-y 0) -Q y (x-x 0)} ... ( number _{21 ') B y = (∇K} ) y {Q y (y-y 0) -Q x (x −x ₀ )｝ + KQ _x (Equation 14 ′) ∂B _Z / ∂y = −2α (∇K) _y ｛Q _x (y−y ₀ ) −Q _y (x−x ₀ )｝ − α KQ _x ] + Β _y {Q _x (y−y ₀ ) −Q _y (x x ₀ )} = − 2αB _y + α KQ _x + β _y {Q _x (y−y ₀ ) −Q _y (x−x ₀ )} ( (Equation 22 ') As is clear from (Equation 13') and (Equation 21 '), ∂
The value of B _Z / ∂x is equal to the sum of a term equal to −2α times the tangential component B _x and two additional terms,
And (Equation 22 ') and as apparent from the value of .differential.B _Z / ∂y is a term equal to -2α times the tangential component B _y, equal to a value obtained by adding the two additional terms.

【００３６】ここで、図２に概略位置を示すように、無
限平面導体Ｐの内部の点〈ｒ₀（０、０、−ｚ₀）〉、ｚ
₀＝０．０５［ｍ］に、電流ダイポールのモーメント
〈Ｑ〉＝（Ｑ_x、Ｑ_y、０）、Ｑ_x＝Ｑ_y＝５０［ｎＡｍ］
が存在する場合に、Ｂ_x（（数１３））と−∂Ｂ_Z／∂ｘ
（（数１６））を比較する。ｘ₀＝ｙ₀＝ｙ＝０、Ｑ_z＝
０を（数１３）、（数１６）に代入して、（数２３）、
（数２４）を得る。Here, as shown schematically in FIG. 2, the points <r ₀ (0, 0, −z ₀ )>, z inside the infinite planar conductor P
_{When 0} = 0.05 [m], the moment of the current dipole <Q> = (Q _x , Q _y , 0), Q _x = Q _y = 50 [nAm]
Exist, B _x ((Equation 13)) and −∂B _Z / ∂x
((Equation 16)) is compared. x ₀ = y ₀ = y = 0, Q _z =
By substituting 0 into (Equation 13) and (Equation 16), (Equation 23),
(Equation 24) is obtained.

【００３７】 Ｂ_x（ｘ、０） ＝｛μ₀／（４πＫ²）｝｛−（∇Ｋ）_xＱ_yｘ＋ＫＱ_y｝ …（数２３） ∂Ｂ_Z（ｘ、０）／∂ｘ＝ ｛μ₀／（４πＫ²）｝｛２α（∇Ｋ）_xＱ_yｘ−αＫＱ_y−β_xＱ_yｘ｝ …（数２４） 図３は、無限平面導体Ｐの上でのＢ_x（（数２３））及
び−∂Ｂ_Z／∂ｘ（（数２４））をそれぞれの最大値で
規格化した相対磁場強度曲線Ｃ₁、Ｃ₂で示す。B _x (x, 0) = {μ ₀ / (4πK ² )} − (∇K) _x Q _y x + KQ _y … (equation 23) ∂B _Z (x, 0) / ∂x = ｛ ^{_{μ 0 / (4πK 2)}}} {2α (∇K) x Q y x-αKQ y -β x Q y x} ... ( Equation 24) 3, B _x ((number of on an infinite plane conductor P 23)) and -∂B _Z / ∂x ((Equation 24)) are shown by relative magnetic field strength curves C ₁ and C ₂ normalized by their respective maximum values.

【００３８】即ち、曲線Ｃ₁はＢ_x（ｘ、０）／ｍａｘ｜
Ｂ_x（ｘ、０）｜を、曲線Ｃ₂は｛−∂Ｂ_Z（ｘ、０）／
∂ｘ｝／ｍａｘ｜∂Ｂ_Z（ｘ、０）／∂ｘ｜を表わす。
図３から明らかなように、Ｂ_x及び−∂Ｂ_Z／∂ｘの分布
は何れも電流ダイポールが存在する真上の原点（ｘ＝
０）にピークを持ち、何れも共に電流ダイポールが存在
する点の真上に計測点がある時に最大の信号を検出可能
であることを示している。また、曲線Ｃ₂の方が曲線Ｃ₁
よりも鋭いピークを与え、−∂Ｂ_Z／∂ｘ（（数１
６））による磁場分布はＢ_x（（数１３））による磁場
分布よりも空間分解能が高いことを示している。That is, the curve C ₁ is represented by B _x (x, 0) / max |
B _x (x, 0) | and the curve C ₂ is ｛−∂B _Z (x, 0) /
{X} / max | ∂B _Z (x, 0) / ∂x |
As is clear from FIG. 3, the distributions of B _x and −∂B _Z / ∂x are both the origin (x =
0) has a peak, and both indicate that the maximum signal can be detected when the measurement point is located immediately above the point where the current dipole exists. Also, the curve C ₂ is better than the curve C ₁
Giving a sharper peak than -∂B _Z / ∂x ((Equation 1
The magnetic field distribution according to 6)) has a higher spatial resolution than the magnetic field distribution according to B _x ((Equation 13)).

【００３９】図４に示す磁場強度曲線Ｃ₃、Ｃ₄、Ｃ₅は
それぞれ、−∂Ｂ_Z（ｘ、０）／∂ｘの第１項、第２
項、第３項を示す。図４に示す結果から、第３項は第１
項及び第２項に対して無視でき、−∂Ｂ_Z（ｘ、０）／
∂ｘの形状は第１項、第２項により決定されていると見
なせ、（数２４）は（数２４’）と近似できる。The magnetic field strength curves C ₃ , C ₄ and C ₅ shown in FIG. 4 are respectively the first term and the second term of −∂B _Z (x, 0) / ∂x.
Item and Item 3 are shown. From the results shown in FIG. 4, the third term is the first
Term and the second term are negligible, −∂B _Z (x, 0) /
It can be considered that the shape of て い る x is determined by the first and second terms, and (Formula 24) can be approximated to (Formula 24 ′).

【００４０】 ∂Ｂ_Z（ｘ、０）／∂ｘ＝ ｛μ₀／（４πＫ²）｝｛２α（∇Ｋ）_xＱ_yｘ−αＫＱ_y｝ …（数２４’） 図５は、（数１３）、（数１６）のそれぞれの第１項と
第２項を規格化の後に比較した相対磁場強度曲線を示す
図である。図５に於いて、曲線Ｃ₆は｛Ｂ_x（ｘ、０）の
第１項｝／ｍａｘ｜Ｂ_x（ｘ、０）｜、即ち、｛−（∇
Ｋ）_xＱ_yｘ｝／ｍａｘ｜Ｂ_x（ｘ、０）｜を表わし、曲
線Ｃ₇は｛−∂Ｂ_Z（ｘ、０）／∂ｘの第１項｝／ｍａｘ
｜∂Ｂ_Z（ｘ、０）／∂ｘ｜、即ち、｛−２α（∇Ｋ）_x
Ｑ_yｘ｝／ｍａｘ｜∂Ｂ_Z（ｘ、０）／∂ｘ｜を表わし、
曲線Ｃ₈は｛Ｂ_x（ｘ、０）の第２項｝／ｍａｘ｜Ｂ
_x（ｘ、０）｜、即ち、｛ＫＱ_y｝／ｍａｘ｜Ｂ_x（ｘ、
０）｜を表わし、曲線Ｃ₉は｛−∂Ｂ_Z（ｘ、０）／∂ｘ
の第２項｝／ｍａｘ｜∂Ｂ_Z（ｘ、０）／∂ｘ｜、即
ち、｛αＫＱ_y｝／ｍａｘ｜∂Ｂ_Z（ｘ、０）／∂ｘ｜を
表わす。∂B _Z (x, 0) / ∂x = ｛μ ₀ / (4πK ² )｝ ｛2α (∇K) _x Q _y x−α KQ _y … (Equation 24 ′) FIG. 13) is a diagram showing a relative magnetic field strength curve comparing the first term and the second term of (Equation 16) after normalization. In FIG. 5, curve C ₆ is the first term of {B _x (x, 0)} / max | B _x (x, 0) |, that is, ｛− (∇
K) _x Q _y x} / max | B _x (x, 0) |, and curve C ₇ is the first term of {− {B _Z (x, 0) / {x}} / max
| ∂B _Z (x, 0) / ∂x |, that is, ｛-2α (∇K) _x
Q _y x｝ / max | ∂B _Z (x, 0) / ∂x |
Curve C ₈ is the second term of {B _x (x, 0)} / max | B
_x (x, 0) |, that is, {KQ _y } / max | B _x (x,
0) |, and the curve C ₉ is represented by ｛−∂B _Z (x, 0) / ∂x
｝ / Max | ∂B _Z (x, 0) / ∂x |, that is, {αKQ _y } / max | ∂B _Z (x, 0) / ∂x |.

【００４１】図５に示す結果から、−∂Ｂ_Z（ｘ、０）
／∂ｘの第１項、第２項の分布は共にそれぞれ、Ｂ
_x（ｘ、０）の第１項、第２項の分布よりも鋭く、分布
の尖鋭度は（数１８）で定義されているα＝（∇Ｋ）_z
／Ｋにより規定されている。From the results shown in FIG. 5, -∂B _Z (x, 0)
The distribution of the first and second terms of / ∂x is B
_x (x, 0) is sharper than the distribution of the first and second terms, and the sharpness of the distribution is α = (∇K) _z defined by (Equation 18)
/ K.

【００４２】図６に於いて、磁場曲線Ｃ₁₀はα＝（∇
Ｋ）_z／Ｋを、磁場曲線Ｃ₁₁は−｛（数２４）の第１
項｝／｛（数２３）の第１項｝、即ち、２α（∇Ｋ）_x
Ｑ_yｘ／（∇Ｋ）_xＱ_yｘ＝２αを、磁場曲線Ｃ₁₂は−
｛（数２４）の第２項｝／｛（数２３）の第２項｝、即
ち、αＫＱ_y／ＫＱ_y＝αをそれぞれ示す。図６に示すよ
うに、α＝（∇Ｋ）_z／Ｋ（曲線Ｃ₁₀）は電流ダイポー
ルが存在する原点にピーク点を有し、ピーク値は２／
（ｚ−ｚ₀）である。−∂Ｂ_Z（ｘ、０）／∂ｘの大きさ
は、Ｂ_x（ｘ、０）の大きさとピーク点で２／（ｚ−
ｚ₀）だけ異なる。（ｚ−ｚ₀）は電流ダイポールの存在
する深さである。実際の磁場計測からは（ｚ−ｚ₀）を
決定することは困難である。（数２３）と（数２４’）
との比較から（数２５）を得る。[0042] In FIG. 6, the magnetic field curve C ₁₀ is alpha = (∇
The K) _z / K, the magnetic field curve C ₁₁ is - first {(number 24)
The first term の of the term｝ / ｛(Equation 23), that is, 2α (∇K) _x
Q _y x / (∇K) _x Q _y x = 2α, and the magnetic field curve C ₁₂ is −
{2nd term of (Expression 24)} / {2nd term of (Expression 23)}, that is, αKQ _y / KQ _y = α. As shown in FIG. 6, α = (∇K) _z / K (curve C ₁₀ ) has a peak point at the origin where the current dipole exists, and the peak value is 2 /
(Z−z ₀ ). The magnitude of −∂B _Z (x, 0) / ∂x is 2 / (z−) at the peak point with the magnitude of B _x (x, 0).
z ₀ ). (Z-z ₀₎ is the depth of the presence of current dipole. It is difficult to determine (z-z ₀ ) from actual magnetic field measurements. (Equation 23) and (Equation 24 ')
(Equation 25) is obtained from the comparison with

【００４３】 −∂Ｂ_Z（ｘ、０）／∂ｘ ＝｛μ₀／（４πＫ²）｝｛−２α（∇Ｋ）_xＱ_yｘ＋αＫＱ_y｝ ＝２αＢ_x（ｘ、０）−｛μ₀／（４πＫ）｝αＱ_y …（数２５） 即ち、（数２５）の第２項が第１項に対して小さい場合
には、近似的に（数２６）が成立すると見做せる。− {B _Z (x, 0) / ∂x = ｛μ ₀ / (4πK ² )} − 2α (∇K) _x Q _y x + α KQ _y ＝= 2αB _x (x, 0) − ｛μ ₀ / (4πK)｝ αQ _y (Equation 25) That is, when the second term of (Equation 25) is smaller than the first term, it can be considered that (Equation 26) is approximately established.

【００４４】 −∂Ｂ_Z（ｘ、０）／∂ｘ＝２αＢ_x（ｘ、０） …（数２６） 一般化して（数２１’）に於いて、−２αＢ_x以外の２
つの付加項が−２αＢ_xに対して小さい場合には、近似
的に（数２７）が成立すると見做せる。−∂B _Z (x, 0) / ∂x = 2αB _x (x, 0) (Equation 26) In generalization (Equation 21 ′), 2 other than −2αB _x
One of the cases additional term is small relative -2ArufaB _x is做causes seen as an approximation (Equation 27) is satisfied.

【００４５】 ∂Ｂ_Z／∂ｘ＝−２αＢ_x …（数２７） 以上では、−∂Ｂ_Z／∂ｘとＢ_xの関係について検討した
結果であるが、同様のことが−∂Ｂ_Z／∂ｙとＢ_yの関係
についても成立し、（数２２’）から近似的に（数２
８）成立すると見做せる。∂B _Z / ∂x = −2αB _x (Equation 27) The above is the result of examining the relationship between −∂B _Z / ∂x and B _x , but the same is true for −∂B _Z / also established the relationship between ∂y and B _y, approximately (number 2 from equation (22 ')
8) It can be considered that it holds.

【００４６】 ∂Ｂ_Z／∂ｙ＝−２αＢ_y …（数２８） 以下、（数２７）、（数２８）からそれぞれ、Ｂ_xは−
∂Ｂ_Z／∂ｘ、Ｂ_yは−∂Ｂ_Z／∂ｙに比例すると仮定し
て、計測された法線成分Ｂ_zから接線成分Ｂ_x、Ｂ_yを推
定して等磁場線図を求める手順を詳細に説明する。∂B _Z / ∂y = −2αB _y (Equation 28) From the following (Equation 27) and (Equation 28), B _x is −
.Differential.B _Z / ∂x, B _y is assumed to be proportional to -∂B _Z / ∂y, tangential components B _x from the measured normal component B _z, obtains the isomagnetic field by estimating a B _y The procedure will be described in detail.

【００４７】生体の面に垂直な磁場成分Ｂ_z（ｘ、ｙ、
ｔ）を計測した場合、Ｂ_z（ｘ、ｙ、ｔ）のｘ方向の変
化率∂Ｂ_z（ｘ、ｙ、ｔ）／∂ｘと、Ｂ_z（ｘ、ｙ、ｔ）
の方向の変化率∂Ｂ_z（ｘ、ｙ、ｔ）／∂ｙと求め、
（数３３）に示すように２乗和の平方根Ｓ（ｘ、ｙ、
ｔ）を求める。A magnetic field component B _z (x, y,
When t) is measured, the rate of change of B _z (x, y, t) in the x direction ∂B _z (x, y, t) / ∂x and B _z (x, y, t)
方向 B _z (x, y, t) / ∂y
As shown in (Equation 33), the square root S (x, y,
Find t).

【００４８】 Ｓ（ｘ、ｙ、ｔ）＝√［｛∂Ｂ_z（ｘ、ｙ、ｔ）／∂ｘ｝² ＋｛∂Ｂ_z（ｘ、ｙ、ｔ）／∂ｙ｝²］ …（数３３） 次いで、各点（ｘ、ｙ）について任意の期間での波形Ｓ
_t（ｔ、ｘ、ｙ）の積分値Ｉ₂（ｘ、ｙ）を（数３４）に
より求め、内挿、外挿により各点（ｘ、ｙ）での積分値
Ｉ₂（ｘ、ｙ）が同じ値の点を結ぶ等積分図を求めて、
等積分図を表示画面に表示する。S (x, y, t) = {[｛∂B _z (x, y, t) / ∂x｝ ² + ｛∂B _z (x, y, t) / ∂y｝ ² ] ( (Equation 33) Next, the waveform S for each point (x, y) in an arbitrary period
_t (t, x, y) integrated value I ₂ of (x, y) and determined by (number 34), the interpolation, the integral value at each point by extrapolation (x, y) I ₂ (x, y) Finds an isointegral diagram connecting points with the same value,
Display the equal integral diagram on the display screen.

【００４９】 Ｉ₂（ｘ、ｙ）＝∫│Ｓ_t（ｘ、ｙ、ｔ）│ｄｔ …（数３４） なお、（数４）、（数３４）の積分範囲としては、例え
ば、心臓を測定の対象とする時には、Ｑ、Ｒ、Ｓの各波
の発生する期間、Ｑ波からＳ波の発生するＱＲＳ波(QRS
complex)の期間、Ｔ波の発生する期間等をとる。更
に、（数４）、（数３４）の積分範囲として複数の積分
範囲をとり求めた複数の積分値の間での、等加重を含む
和又は差、比を求める等の演算を行ない、内挿、外挿に
より演算結果が同じ値の点を結ぶ等積分図を求めて、等
積分図を表示画面に表示する。例えば、第１の積分範囲
としてＱＲＳ波の発生する期間Ｔ₁、第２の積分範囲と
してＴ波の発生する期間Ｔ₂を設定し、（数４）に従っ
て、期間Ｔ₁に関する積分値Ｉ ₁、_T1（ｘ、ｙ）、Ｉ₂、
_T1（ｘ、ｙ）、（数３４）に従って、期間Ｔ₂に関する
積分値Ｉ₁、_T2（ｘ、ｙ）、Ｉ₂、_T2（ｘ、ｙ）をそれぞ
れを求め、積分値Ｉ₁、_T1（ｘ、ｙ）と積分値Ｉ₁、
_T2（ｘ、ｙ）との間、又は積分値Ｉ₂、_T1（ｘ、ｙ）と
積分値Ｉ₂、_T2（ｘ、ｙ）との間で、等加重を含む和Ｉ
_sum（ｘ、ｙ）、又は差Ｉ_dif（ｘ、ｙ）、比ｒ（ｘ、
ｙ）を、（数３５）〜（数３６）、（数３７）〜（数３
８）、（数３９）〜（数４０）に従って求める。I_Two(X, y) = ∫ | S_t(X, y, t) │dt (Equation 34) The integration ranges of (Equation 4) and (Equation 34) are, for example, as follows.
For example, when measuring the heart, the Q, R, and S waves
Is generated, the QRS wave (QRS
 complex), the period during which a T-wave is generated, and the like. Change
And a plurality of integrals as integral ranges of (Equation 4) and (Equation 34)
Including equal weights between multiple integrals over a range
Perform calculations such as sum or difference, ratio, etc., for interpolation and extrapolation.
Find an equal integral diagram that connects points of the same value
Display the integral diagram on the display screen. For example, the first integration range
The period T during which the QRS wave is generated₁, The second integration range and
The period T in which the T-wave occurs_TwoAnd follow (Equation 4)
And period T₁Integral I ₁,_T1(X, y), I_Two,
_T1According to (x, y), (Equation 34), the period T_TwoAbout
Integral value I₁,_T2(X, y), I_Two,_T2(X, y)
And the integral value I₁,_T1(X, y) and integral value I₁,
_T2(X, y) or integrated value I_Two,_T1(X, y) and
Integral value I_Two,_T2(X, y) and the sum I including equal weights
_sum(X, y) or difference I_dif(X, y), the ratio r (x,
y) is expressed by (Equation 35) to (Equation 36), (Equation 37) to (Equation 3)
8), calculated according to (Expression 39) to (Expression 40).

【００５０】 Ｉ_sum（ｘ、ｙ）＝ ｗ₁×Ｉ₁、_T1（ｘ、ｙ）＋ｗ₂×Ｉ₁、_T2（ｘ、ｙ） …（数３５） Ｉ_sum（ｘ、ｙ）＝ ｗ₁×Ｉ₂、_T1（ｘ、ｙ）＋ｗ₂×Ｉ₂、_T2（ｘ、ｙ） …（数３６） Ｉ_dif（ｘ、ｙ）＝ ｗ₂×Ｉ₁、_T2（ｘ、ｙ）−ｗ₁×Ｉ₁、_T1（ｘ、ｙ） …（数３７） Ｉ_dif（ｘ、ｙ） ＝ｗ₂×Ｉ₂、_T2（ｘ、ｙ）−ｗ₁×Ｉ₂、_T1（ｘ、ｙ） …（数３８） ｒ（ｘ、ｙ）＝Ｉ₁、_T1（ｘ、ｙ）／Ｉ₁、_T2（ｘ、ｙ） …（数３９） ｒ（ｘ、ｙ）＝Ｉ₂、_T1（ｘ、ｙ）／Ｉ₂、_T2（ｘ、ｙ） …（数４０） （数３５）〜（数３６）、（数３７）〜（数３８）、
（数３９）〜（数４０）の演算の結果、個人差による等
積分図のばらつきが改善され、疾患等による生体機能の
異常を検出できる。I _sum (x, y) = w ₁ × I ₁ , _{T 1} (x, y) + w ₂ × I ₁ , _{T 2} (x, y) ( _Equation 35) I _sum (x, y) = w ₁ × I ₂ , _{T 1} (x, y) + w ₂ × I ₂ , _{T 2} (x, y) ( _Expression 36) I _dif (x, y) = w ₂ × I ₁ , _{T 2} (x, y) −w ₁ × I ₁ , _{T 1} (x, y) ( _Expression 37) I _dif (x, y) = w ₂ × I ₂ , _{T 2} (x, y) −w ₁ × I ₂ , _{T 1} (x, y) ( (38) r (x, y) = I ₁ , _T1 (x, y) / I ₁ , _T2 (x, y) (expression 39) r (x, y) = I ₂ , _T1 (x, y) / I ₂ , _T2 (x, y) (Equation 40) (Equation 35) to (Equation 36), (Equation 37) to (Equation 38),
As a result of the calculations of (Equation 39) to (Equation 40), variations in the isointegral diagram due to individual differences are improved, and abnormalities in biological functions due to diseases and the like can be detected.

【００５１】本発明で得られる等積分図によれば、従来
技術で必要としていた生体部位の各時刻に於ける状態を
表わす多数の図（マップ）を用いて生体現象を解析する
ことなく、従来技術で必要としていた図（マップ）の数
よりもはるかに少数の図（マップ）を用いて、生体部位
の全体の状態を把握できる。また、生体磁場の接線成
分、又は法線成分を用いて得られる等積分図のピーク位
置と、生体内で電流が多く流れる部位が一致するので、
等積分図から任意の時間帯での生体内のどの部位で多く
電流が流れたかを判別できる。生体磁場分布は個人差が
大きいが、本発明では、生体磁場の各方向成分の時間変
化を表わす波形から得る任意の時間（期間）での積分値
を用いるので、より定量的な生体磁場分布を少数の図
（マップ）を用いて表示でき、個人毎の疾患、異常を客
観的、定量的に把握できる。According to the isointegral diagram obtained in the present invention, the biological phenomena can be analyzed without analyzing the biological phenomena using a large number of maps (maps) representing the state of the living body part at each time required in the prior art. The entire state of the living body part can be grasped using a much smaller number of maps (maps) than the number of maps (maps) required by the technology. In addition, since the tangent component of the biomagnetic field, or the peak position of the isointegral diagram obtained using the normal component, coincides with the portion where a large amount of current flows in the living body,
From the isointegral diagram, it is possible to determine at which part in the living body a large amount of current flows in an arbitrary time zone. Although the biomagnetic field distribution has a large individual difference, the present invention uses an integral value at an arbitrary time (period) obtained from a waveform representing the time change of each direction component of the biomagnetic field. It can be displayed using a small number of figures (maps), and it is possible to objectively and quantitatively grasp diseases and abnormalities for each individual.

【００５２】本発明では、生体の面に垂直な磁場成分Ｂ
_z（ｘ、ｙ、ｔ）を計測して、Ｂ_xをＢ_z（ｘ、ｙ、ｔ）
のｘ方向の変化率∂Ｂ_z（ｘ、ｙ、ｔ）／∂ｘから、Ｂ_y
をＢ_z（ｘ、ｙ、ｔ）の方向の変化率∂Ｂ_z（ｘ、ｙ、
ｔ）／∂ｙから推定して求めるので、隣接する各計測点
（ｘ、ｙ）に共通して存在する背景となる磁場（妨害磁
場）は、ｘ方向、及びｙ方向で各々キャンセルされるこ
ととなる。In the present invention, the magnetic field component B perpendicular to the surface of the living body
Measure _z (x, y, t) and convert B _x to B _z (x, y, t)
In the x direction of the change rate _{∂B z (x, y, t} ) / from ∂x, B _y
Is the rate of change in the direction of B _z (x, y, t) ∂B _z (x, y,
t) / ∂y, the background magnetic field (disturbing magnetic field) common to adjacent measurement points (x, y) is canceled in the x and y directions, respectively. Becomes

【００５３】[0053]

【発明の実施の形態】生体磁場計測に於ける座標系とし
て直交座標系（ｘ、ｙ、ｚ）（磁場成分をＢ_x、Ｂ_y、Ｂ
_zとする）や極直交座標系（ｒ、θ、φ）が用いられ
る。計測対象が心臓等である場合には、胸壁をｘｙ平面
とする直交座標系（ｘ、ｙ、ｚ）が用いられる。計測対
象が脳部等である場合には、頭部が球に近い形状である
ため極座標系（ｒ、θ、φ）（磁場成分をＢ_r、Ｂ_θ、
Ｂ_φとする）が用いられる。本実施例では、生体表面に
垂直な磁場成分（法線成分）はＢ_z、Ｂ_rで表わされ、生
体の面に平行な成分（接線成分）は、Ｂ_x、Ｂ_y、Ｂ_θ、
Ｂ_φで表わされる。以下、本実施例では、直交座標系
（ｘ、ｙ、ｚ）を用いて説明するが、極座標系（ｒ、
θ、φ）を用いる場合には、Ｂ_zをＢ_rに、Ｂ_xをＢ
_θに、Ｂ_yをＢ_φにそれぞれ読み替えれば良い。DESCRIPTION OF THE PREFERRED EMBODIMENTS As a coordinate system for measuring a biomagnetic field, a rectangular coordinate system (x, y, z) (magnetic field components are represented by B _x , B _y , B
_z ) and a polar orthogonal coordinate system (r, θ, φ) are used. When the measurement target is a heart or the like, an orthogonal coordinate system (x, y, z) in which the chest wall is an xy plane is used. When the measurement target is the brain or the like, since the head has a shape close to a sphere, the polar coordinate system (r, θ, φ) (magnetic field components are _Br , _Bθ ,
B _φ ) is used. In this embodiment, magnetic field component perpendicular to the biological surface (normal component) is represented by B _z, B _r, parallel component (tangential component) in a surface of the living _{body, B x, B y, B} θ,
Represented by B _φ. Hereinafter, the present embodiment will be described using the rectangular coordinate system (x, y, z), but the polar coordinate system (r,
θ, φ), B _z is B _r and B _x is B
to _theta, it may be read as respectively B _y in B _phi.

【００５４】図７は本発明が実施される生体磁場計測装
置の概略構成を示す。心磁場計測を行なう生体磁場計測
装置は、量子干渉素子（ＳＱＵＩＤ）からなる複数の磁
場センサを用いる。環境磁場雑音の影響を除去するため
に、心磁場計測は磁場シールドルーム１の内部で行なわ
れる。被検者２はベッド３に横たわり計測する（図１１
に示すように、ｘｙ面がベッドの面となるように直交座
標系（ｘ、ｙ、ｚ）を設定する）。被検者２の胸部の上
方に、ＳＱＵＩＤとそのＳＱＵＩＤに接続した検出コイ
ルとが一体化された磁場センサを複数個収納し、液体Ｈ
ｅを満たしたデュワ４が配置される。液体Ｈｅは磁場シ
ールドルーム１の外部の自動補給装置５により、連続的
に液体Ｈｅが補充されている。FIG. 7 shows a schematic configuration of a biomagnetic field measuring apparatus according to the present invention. A biomagnetic field measurement device that performs a cardiac magnetic field measurement uses a plurality of magnetic field sensors including a quantum interference device (SQUID). In order to remove the influence of the environmental magnetic field noise, the cardiac magnetic field measurement is performed inside the magnetic field shield room 1. The subject 2 lays down on the bed 3 and measures (FIG. 11)
(X, y, z) is set so that the xy plane is the bed plane as shown in FIG. Above the chest of the subject 2, a plurality of magnetic field sensors in which a SQUID and a detection coil connected to the SQUID are integrated are stored.
The dewar 4 satisfying e is arranged. The liquid He is continuously replenished with the liquid He by the automatic replenishing device 5 outside the magnetic field shield room 1.

【００５５】磁場センサからの出力は、検出コイルが検
出した磁場強度に比例する電圧を出力するＦＬＬ（Flux
Locked Loop）回路６に入力される。このＦＦＬ回路は
ＳＱＵＩＤの出力を一定に保つようＳＱＵＩＤに入力さ
れた生体磁場の変化を帰還コイルを介してキャンセルし
ている。この帰還コイルに流した電流を電圧に変換する
ことにより、生体磁場信号の変化に比例した電圧出力が
得られる。この電圧出力は、増幅器（図示せず）により
増幅され、フイルター回路７により周波数帯域が選択さ
れ、ＡＤ変換器で（図示せず）ＡＤ変換され、計算機８
に取り込まれる。計算機８では、各種の演算処理が実行
され、演算処理結果がデイスプレイに表示され、更に、
プリンタにより出力される。The output from the magnetic field sensor is FLL (Flux) which outputs a voltage proportional to the magnetic field strength detected by the detection coil.
Locked Loop) circuit 6. This FFL circuit cancels the change of the biomagnetic field input to the SQUID via the feedback coil so as to keep the output of the SQUID constant. By converting the current passed through the feedback coil into a voltage, a voltage output proportional to a change in the biomagnetic field signal can be obtained. This voltage output is amplified by an amplifier (not shown), a frequency band is selected by a filter circuit 7, and AD converted (not shown) by an AD converter (not shown).
It is taken in. In the calculator 8, various calculation processes are executed, the calculation process results are displayed on a display, and furthermore,
Output by printer.

【００５６】磁場の接線成分を検出する検出コイルとし
て、コイル面がｘ方向、及びｙ方向を向いた２つのコイ
ルを使用し、磁場の法線成分を検出する検出コイルとす
る。また磁場の法線成分を検出するコイルとしては、ｚ
方向を向いたコイルを使用する。これら磁場センサ（２
０−１、２０−２、〜、２０−８、２１−１、〜、２１
−８、２２−１、〜、２２−８、２３−２、〜、２３−
８、２４−１、〜、２４−８、２５−１、〜、２５−
８、２６−１、〜、２６−８、２７−１、〜、２７−
８）の配置図を図８に示す。磁場センサ９はデュワ内部
の底部から垂直の方向に設置し、また各センサ間は磁場
の距離変化量を正確に捕らえるようにｘ方向、ｙ方向に
等間隔になるようにした。ここで、センサ間距離は２５
ｍｍとし、センサ数は８×８の６４チャンネルとした。As the detection coil for detecting the tangential component of the magnetic field, two coils whose coil surfaces are oriented in the x direction and the y direction are used to detect the normal component of the magnetic field. As a coil for detecting the normal component of the magnetic field, z
Use oriented coils. These magnetic field sensors (2
0-1, 20-2, ..., 20-8, 21-1, ..., 21
-8, 22-1, ~ 22-8, 23-2, ~ 23
8, 24-1, ~, 24-8, 25-1, ~, 25-
8, 26-1, ..., 26-8, 27-1, ..., 27-
FIG. 8 shows the layout of 8). The magnetic field sensors 9 were installed in a direction perpendicular to the bottom inside the Dewar, and the sensors were arranged at equal intervals in the x and y directions so as to accurately capture the distance change of the magnetic field. Here, the distance between the sensors is 25
mm and the number of sensors was 64 channels of 8 × 8.

【００５７】この配列方法に従って、設置した磁場セン
サの１本の概略図を図９及び図１０に示す。図９の磁場
センサは生体表面に対して垂直な成分Ｂ_zを測定するセ
ンサで、超伝導線（ＮｂーＴｉ線）で作製したコイルの
面がｚ方向を向いている。このコイルは２つの逆向きの
コイルを組み合わせたもので生体に近い方を検出コイル
１０とし、遠い方のコイルを外部磁場雑音を除去する参
照コイル(reference coil)１１とし１次微分コイルを形
成している。ここでコイル径を２０ｍｍφ、コイル間の
ベースラインを５０ｍｍとした。外部磁場雑音は生体よ
り遠い信号源から生じており、これらは検出コイル及び
参照コイルで同じように検出される。一方、生体からの
信号はコイルに近いため検出コイル１０でより強く検出
される。このため、検出コイル１０では信号と雑音が検
出され、参照コイル１１では雑音のみが検出される。従
って、両者のコイルで捕らえた磁場の差を取ることによ
りＳ／Ｎの高い計測ができる。FIGS. 9 and 10 show schematic views of one magnetic field sensor installed according to this arrangement method. The magnetic field sensor shown in FIG. 9 is a sensor for measuring a component _Bz perpendicular to the surface of a living body, and the surface of a coil made of a superconducting wire (Nb-Ti wire) is oriented in the z direction. This coil is a combination of two opposite coils. The one closer to the living body is used as the detection coil 10, and the one farther away is used as a reference coil 11 for removing external magnetic field noise to form a primary differential coil. ing. Here, the coil diameter was 20 mmφ, and the baseline between the coils was 50 mm. External magnetic field noise originates from signal sources farther from the body, and these are similarly detected by the detection coil and the reference coil. On the other hand, since the signal from the living body is close to the coil, it is detected more strongly by the detection coil 10. Therefore, the detection coil 10 detects a signal and noise, and the reference coil 11 detects only noise. Therefore, a high S / N measurement can be performed by taking the difference between the magnetic fields captured by the two coils.

【００５８】１次微分コイルはＳＱＵＩＤ１２を実装し
た実装基板の超伝導配線を介してＳＱＵＩＤのインプッ
トコイルに接続し、コイルで検出した生体磁場をＳＱＵ
ＩＤに伝達する。生体磁場成分の接線成分Ｂ_x、Ｂ_yを検
出する磁場センサの概略図を図１０に示す。この磁場セ
ンサは平面型のコイルを使用しており、検出コイル１
０’、１０”及び参照コイル１１’、１１”が１つの平
面に並び、コイル径は２０ｍｍ×２０ｍｍ、ベースライ
ンは５０ｍｍとした。コイルは法線成分用と同様にＳＱ
ＵＩＤ１２’、１２”の実装基板に接続している。４角
柱の支持体の互いに直交する２面に、これらｘ成分検出
用磁場センサ１３とｙ成分検出用磁場センサ１４を張り
付けることにより、ｘ及びｙ成分を測定できる磁場セン
サを形成している。この４角柱は図８に示すようにアレ
イ状に配置した。The primary differential coil is connected to the input coil of the SQUID through the superconducting wiring of the mounting board on which the SQUID 12 is mounted, and the biomagnetic field detected by the coil is applied to the SQUID.
Communicate to ID. Tangential component B _x biomagnetic field components, a schematic diagram of a magnetic field sensor for detecting the B _y shown in FIG. 10. This magnetic field sensor uses a planar coil, and a detection coil 1
0 ′, 10 ″ and reference coils 11 ′, 11 ″ were arranged in one plane, the coil diameter was 20 mm × 20 mm, and the baseline was 50 mm. The coil is SQ like the normal component
It is connected to the mounting substrates of UIDs 12 'and 12 ". By attaching the x-component detecting magnetic field sensor 13 and the y-component detecting magnetic field sensor 14 to two mutually orthogonal surfaces of a quadrangular prism support, x And a magnetic field sensor capable of measuring the y component.The quadrangular prisms are arranged in an array as shown in FIG.

【００５９】磁場センサを内蔵したデュワは、ベットに
横たわった被験者の胸部上方に配置し心臓から発生する
磁場を計測する。ここで、体の横方向をｘ軸とし、体の
上下方向をｙ軸とする。磁場センサ（２０−１、〜、２
０−８、２１−１、〜、２１−８、２２−１、〜、２２
−８、２３−２、〜、２３−８、２４−１、〜、２４−
８、２５−１、〜、２５−８、２６−１、〜、２６−
８、２７−１、〜、２７−８）の配置と胸部３０との位
置関係を図１１に示す。この位置関係で計測した生体磁
場信号を図１２（ａ）、（ｂ）、（ｃ）に示す。The Dewar with a built-in magnetic field sensor is arranged above the chest of the subject lying on the bed and measures the magnetic field generated from the heart. Here, the horizontal direction of the body is defined as the x-axis, and the vertical direction of the body is defined as the y-axis. Magnetic field sensors (20-1, ..., 2
0-8, 21-1, ..., 21-8, 22-1, ..., 22
-8, 23-2, ~, 23-8, 24-1, ~, 24
8, 25-1, ..., 25-8, 26-1, ..., 26-
8, 27-1,..., 27-8) and the positional relationship with the chest 30 are shown in FIG. 12 (a), 12 (b) and 12 (c) show biomagnetic field signals measured in this positional relationship.

【００６０】図１２（ａ）、（ｂ）、（ｃ）は、各磁場
センサ（８×８のアレイ状に並んだ磁場センサ）によ
る、ある健常者の心臓から発する磁場の時間変化を示
し、各図の中の６４個の波形の横軸が時間軸、縦軸が検
出された磁場強度を示している。図１２（ａ）は接線成
分Ｂ_x、図１２（ｂ）は接線成分Ｂ_y、図１２（ｃ）は法
線成分Ｂ_z、の各成分の時間（横軸）の変化を、各磁場
成分毎に信号強度の最も大きいチャンネルの絶対値の最
大値で規格化して表示している。FIGS. 12 (a), 12 (b) and 12 (c) show the time change of the magnetic field emitted from the heart of a healthy person by each magnetic field sensor (magnetic field sensors arranged in an 8 × 8 array). The horizontal axis of the 64 waveforms in each figure is the time axis, and the vertical axis is the detected magnetic field strength. FIG. 12A shows the change in time (horizontal axis) of each component of the tangential component B _x , FIG. 12B shows the tangential component B _y , and FIG. 12C shows the normal component B _z . Each channel is standardized and displayed with the maximum value of the absolute value of the channel having the largest signal strength.

【００６１】図１３に示す点線、実線は、健常者につい
て計測された特定の２チャンネルに関する接線成分（Ｂ
_x）の時間波形を実線、点線で示している。心臓の心室
が脱分極したＱＲＳ波が出現する時間帯Ｔ₁でのＱ波、
Ｒ波、及びＳ波のピーク（極値）を与える時点を図１３
中にそれぞれｔ_Q、ｔ_R、ｔ_sで示した。また、心臓の再
分極過程であるＴ波の出現する時間帯Ｔ₂とし、ピーク
（極値）を与える時点をｔ_Tで示した。The dashed and solid lines shown in FIG. 13 indicate tangent components (B
The time waveform of _x ) is shown by a solid line and a dotted line. A Q wave in a time zone T ₁ at which a QRS wave in which the ventricle of the heart is depolarized appears;
FIG. 13 shows the points at which the peaks (extreme values) of the R wave and the S wave are given.
It is shown in t _Q, t _R, t _s, respectively in. Further, the time period T ₂ to the appearance of T-wave is repolarization of the heart, showing a point that gives a peak (extreme value) at t _T.

【００６２】図１３に於いて、Ｐ波は心房の興奮（脱分
極(depolarization)）を示し、Ｑ波、Ｒ波、及びＳ波か
らなるＱＲＳ波は心室の興奮（脱分極）を示し、Ｔ波は
ＱＲＳ波に続いて出現するゆくやかなふれであり、心筋
の再分極(repolarization)を示している。脱分極は、は
じめに筋の中を興奮が広がる過程であり、再分極は、興
奮した筋が静止状態に戻る過程である。In FIG. 13, a P wave indicates atrial excitation (depolarization), a QRS wave composed of Q, R, and S waves indicates ventricular excitation (depolarization). The wave is a gradual shake that appears following the QRS wave, indicating repolarization of the myocardium. Depolarization is the process by which the excitement initially spreads through the muscle, and repolarization is the process by which the excited muscle returns to a resting state.

【００６３】図１４（ａ）、（ｂ）、（ｃ）は、ｔ_Q、
ｔ_R、ｔ_sの時点での心磁場の強度の等しい点を線で結ん
だ等磁場線図を示す。図１４（ａ）、（ｂ）、（ｃ）
は、（数４）の│Ｂ_xy（ｘ、ｙ、ｔ）│で示され、６４
個所で計測された接線成分Ｂ_x、Ｂ_yを合成した２次元の
ベクトル強度分布を示している。更に、図１４（ａ）、
（ｂ）、（ｃ）中の矢印は、６４個所の各測定点での電
流源が各測定点での磁場を作っているものとして仮定し
た時の２次元の電流ベクトルを示している。この電流ベ
クトルにより心臓内での電流方向及び分布が推定でき
る。図１４（ａ）、（ｂ）、（ｃ）の各図の横軸ｘ、縦
軸ｙは磁場センサが配置されている座標を示す。図１４
（ａ）に示すように、Ｑ波のピーク時では、心臓内を流
れる電流は心室中隔で右下方向に流れ、図１４（ｂ）に
示すように、Ｒ波のピーク時では、左心室全体で斜め下
方向に電流が大きく流れ、図１４（ｃ）に示すように、
Ｓ波のピーク時では、心室基部の方向の左斜め上方向に
電流が流れ、心室の脱分極過程が終了することが分か
る。このように、図１４（ａ）、（ｂ）、（ｃ）の等磁
場線図により各時間での心臓内の活動部位及び電流方向
が可視化できることが分かる。FIGS. 14A, 14B and 14C show t _Q ,
The isomagnetic field diagram which connected the point where the intensity | strength of the magnetocardiogram at t _R and t _s is equal by a line is shown. 14 (a), (b), (c)
Is represented by | B _xy (x, y, t) | in (Equation 4), and 64
Tangential component B _x which is measured at the point, shows a two-dimensional vector intensity distribution obtained by combining the B _y. Further, FIG.
Arrows in (b) and (c) show two-dimensional current vectors when it is assumed that current sources at 64 measurement points generate a magnetic field at each measurement point. With this current vector, the current direction and distribution in the heart can be estimated. 14A, 14B, and 14C, the horizontal axis x and the vertical axis y indicate coordinates where the magnetic field sensor is arranged. FIG.
As shown in FIG. 14A, at the time of the peak of the Q wave, the current flowing in the heart flows in the lower right direction in the septum of the ventricle, and as shown in FIG. As a whole, a large current flows obliquely downward, and as shown in FIG.
At the time of the peak of the S wave, it can be seen that a current flows in a diagonally upper left direction in the direction of the ventricle base, and the ventricular depolarization process ends. Thus, it can be seen that the active site and the current direction in the heart at each time can be visualized by the isomagnetic field diagrams of FIGS. 14 (a), (b) and (c).

【００６４】図１５は、心磁波形のＱ波からＳ波までの
ＱＲＳ波が出現する時間帯Ｔ₁に於いて検出された２つ
の接線成分Ｂ_x、Ｂ_yから得た２次元ベクトル強度│Ｂ_xy
（ｘ、ｙ、ｔ）│を各点（ｘ、ｙ）について、（数４）
の積分を行ない、同じ積分値の点を結んだ等積分図であ
る。図１５のｘ軸、ｙ軸は、生体表面に配置された磁場
センサの座標を表し、等積分図の各曲線の黒丸の近傍に
示した数値はその曲線のもつ積分値を示す。図１５か
ら、ＱＲＳ波の時間帯に心筋に流れた電流の多くは心筋
の厚みが大きい左心室で流れたことが分かり、等積分図
でのピーク位置と心臓に流れる電流量の多い部位とがよ
く対応することが分かった。[0064] Figure 15, the two tangential components detected at a time zone T ₁ which QRS wave from the Q wave of the magnetocardiogram waveforms until S wave appears B _x, 2-dimensional vector magnitude obtained from B _y │ B _xy
(X, y, t) | for each point (x, y) (Equation 4)
Is an integral diagram obtained by integrating the points and connecting points having the same integral value. The x-axis and y-axis in FIG. 15 represent the coordinates of the magnetic field sensor disposed on the surface of the living body, and the numerical values shown near the black circles of the respective curves in the isometric diagram indicate the integrated values of the curves. From FIG. 15, it can be seen that most of the current flowing to the myocardium during the time zone of the QRS wave flows in the left ventricle where the thickness of the myocardium is large. It turns out that it corresponds well.

【００６５】図１６は、図１２（ａ）、（ｂ）、（ｃ）
から図１５のデータを求めたのと同一の健常者につい
て、法線線分Ｂ_zを各点（ｘ、ｙ）に於いて計測し、
（数３３）によりＳ（ｘ、y）を求め、ＱＲＳ波の時間
帯Ｔ₁について、（数３４）の積分を行ない同じ積分値
の点を結んだ等積分図である。以下、図１６から図２１
に於いて、ｘ軸、ｙ軸は、生体表面に配置された磁場セ
ンサの位置座標（単位はｍである）を表わす。図１６か
ら図２１の曲線の黒丸の近傍に示した数値はその曲線の
もつ積分値を示す。FIGS. 16A, 12B, 12C, and 12C are views of FIG.
For the same healthy subjects and were determined the data of FIG. 15, at the normal line B _z at each point (x, y) measured from,
FIG. 34 is an isointegral diagram in which S (x, y) is obtained by (Equation 33), the integration of (Equation 34) is performed for the time zone T ₁ of the QRS wave, and points of the same integration value are connected. Hereinafter, FIGS. 16 to 21
, The x-axis and the y-axis represent the position coordinates (the unit is m) of the magnetic field sensor arranged on the surface of the living body. The numerical values shown in the vicinity of the black circles of the curves in FIGS. 16 to 21 indicate the integral values of the curves.

【００６６】図１５に示す磁場の接線成分Ｂ_x、Ｂ_yから
求めた等積分図と、図１６に示す磁場の法線成分Ｂ_zか
ら求めた等積分図のパターンは一致することが判明し
た。この一致は、（数６）及び（数７）、又は（数２
７）及び（数２８）が実際の実験データでほぼ成立して
いることを意味している。[0066] the tangential component of the magnetic field shown in FIG. 15 B _x, and isointegral obtained from B _y, the pattern of isointegral obtained from normal component B _z of the magnetic field shown in FIG. 16 to be consistent found . This match is expressed by (Equation 6) and (Equation 7), or (Equation 2)
It means that 7) and (Equation 28) are almost satisfied with actual experimental data.

【００６７】図１７は、図１５を求めたのと同一の健常
者について、Ｔ波の時間帯Ｔ₂に於いて検出された２つ
の接線成分Ｂ_x、Ｂ_yから得た２次元ベクトル強度│Ｂ_xy
（ｘ、ｙ）│を各点（ｘ、ｙ）について、（数４）の積
分を行ない同じ積分値の点を結んだ等積分図である。図
１７に於いて、１ｅ＋００３は、１０００を示す。[0067] Figure 17, the same healthy subjects and were determined to 15, the two tangential components B _x detected at the time slot T ₂ of the T-wave, two-dimensional vector magnitude │ obtained from B _y B _xy
(X, y) | is an integral diagram obtained by performing integration of (Equation 4) for each point (x, y) and connecting points having the same integration value. In FIG. 17, 1e + 003 indicates 1000.

【００６８】図１８は、時間帯Ｔ₂についての（数４）
の積分値と、ＱＲＳ波が発生した期間帯Ｔ₁についての
（数４）の積分値との差（数３７）を表わす等高線図で
ある。即ち、図１８は図１７に示す等積分図から図１５
に示す等積分図を差し引いた図である。Ｔ波の時間帯Ｔ
₂の方が、ＱＲＳ波の時間帯Ｔ₁よりも長い。また、図１
７のパターンは、図１５に示すパターンと似ている。こ
のため、図１８に示す等高線図は全体が正の値となる。
図１７、図１８の曲線の黒丸の近傍に示した数値はその
曲線のもつ上記の積分値の差の値を示す。FIG. 18 shows (Equation 4) for the time zone T _2.
And integral value of a contour plot representative of the difference (number 37) of the integrated value of (number 4) for the period zone T ₁ which QRS wave occurs. That is, FIG. 18 is based on the isometric diagram shown in FIG.
3 is a diagram obtained by subtracting an equal integral diagram shown in FIG. Time zone T of T wave
_2/5 is longer than the time period T ₁ of the QRS wave. FIG.
The pattern of No. 7 is similar to the pattern shown in FIG. Therefore, the contour map shown in FIG. 18 has positive values as a whole.
Numerical values shown in the vicinity of the black circles in the curves in FIGS. 17 and 18 indicate the difference between the above-mentioned integral values of the curves.

【００６９】次ぎに、心筋梗塞の患者の心磁場計測に関
する結果を、図１９、図２０、図２１に示す。図１９
は、ＱＲＳ波の時間帯Ｔ₁について図１５と同様にして
求めた等積分図、図２０は、Ｔ波の時間帯Ｔ₂について
図１７と同様にして求めた等積分図、図２１は、Ｔ波の
時間帯Ｔ₂についての積分値（数４）と、ＱＲＳ波の時
間帯Ｔ₁についての積分値（数４）との差（数３５）を
表わし、図１８と同様にして求めた等高線図である。即
ち、図２１は、図２０に示す等積分図から図１９に示す
等積分図を差し引いた図である。図１９、図２０の曲線
の黒丸の近傍に示した数値はその曲線のもつ積分値を示
し、図２１の曲線の黒丸の近傍に示した数値はその曲線
の持つ上記の積分値の差の値を示す。Next, FIG. 19, FIG. 20, and FIG. 21 show the results regarding the measurement of the magnetocardiogram of the patient with myocardial infarction. FIG.
Is isointegral obtained in the same manner as in FIG. 15 for the time period T ₁ of the QRS wave, FIG. 20, isointegral obtained in the same manner as in FIG. 17 for the time period T ₂ of the T-wave, Fig. 21, integral values for the time slot T ₂ of the T-wave and (Equation 4), represents the difference (number 35) of the integral value for the time period T ₁ of the QRS wave (number 4) were determined in the same manner as FIG. 18 It is a contour map. That is, FIG. 21 is a diagram obtained by subtracting the isometric diagram shown in FIG. 19 from the isometric diagram shown in FIG. The numerical values shown in the vicinity of the black circles of the curves in FIGS. 19 and 20 indicate the integrated values of the curves, and the numerical values shown in the vicinity of the black circles of the curves in FIG. 21 are the values of the difference of the above integrated values of the curves. Is shown.

【００７０】図１９に示す時間帯Ｔ₁での等積分図は、
図１５及び図１６に示す等積分図とあまり差のないパタ
ーンであり、左心室に電流が多く流れたことが分かる。
しかし、図２０に示す時間帯Ｔ₂での等積分図は、図１
９に示す時間帯Ｔ₁での等積分図とは異なるパターンと
なり、心筋梗塞のために、時間帯Ｔ₁と時間帯Ｔ₂では心
臓に流れる電流量のパターンが大きく異なることが明確
に分かる。更に、図２１に示す等高線図は全体が負の値
をもち、全体が正の値をもつ図１８に示す健常者の等高
線図とは大きく異なり、心筋梗塞の患者では、時間帯Ｔ
₂で心臓に流れる電流が障害を受けていることが明確に
分かる。The isointegral diagram in the time zone T ₁ shown in FIG.
The patterns are not so different from the equal integration diagrams shown in FIGS. 15 and 16, and it can be seen that a large amount of current flows in the left ventricle.
However, isointegral in the time zone T ₂ shown in FIG. 20, FIG. 1
The pattern is different from the isointegral diagram in the time zone T ₁ shown in FIG. 9, and it is clearly seen that the pattern of the amount of current flowing through the heart differs greatly between the time zone T ₁ and the time zone T ₂ due to myocardial infarction. Further, the contour map shown in FIG. 21 has a negative value as a whole and is significantly different from the contour map of a healthy person shown in FIG. 18 having a positive value as a whole.
_In Figure ₂ , it can be clearly seen that the current flowing to the heart is damaged.

【００７１】以上説明したように、心臓の時間帯Ｔ₁と
時間帯Ｔ₂に於ける磁場強度を画像化するすることによ
り、患者に苦痛を与えることなく非侵襲的に、１分以下
の短時間で、健康な状態と異常な状態（例えば、心筋梗
塞の状態、虚血状態等）とを容易に判別できる。即ち、
逆問題を解くことな疾患部位の早期発見、推定が可能と
なる。As described above, by imaging the magnetic field strength in the time zone T ₁ and the time zone T ₂ of the heart, the patient can be invasively non-invasively without pain for less than one minute. With time, a healthy state and an abnormal state (for example, a state of myocardial infarction, an ischemic state, etc.) can be easily distinguished. That is,
This enables early detection and estimation of a disease site that does not solve the inverse problem.

【００７２】図２２には生体磁場計測装置のコンピュー
タの画面上での処理画像例を示す。マルチウィンド形式
になっており、各処理画像をそれぞれのウィンド上に表
示できる。また、先に説明した図１５から図２１では磁
場強度や積分値の高低がわかるように各曲線に数値を入
れたが、ディスプレイ上では等高線の高低によって色分
けをして３次元カラー表示している。同時に、図１３に
示すように時間波形（心磁図）、更には心電図も表示で
きるようになっており、心臓疾患に関する総合的な解析
ができるようにしている。FIG. 22 shows an example of a processed image on the screen of the computer of the biomagnetic field measuring apparatus. It is a multi-window format, and each processed image can be displayed on each window. In FIG. 15 to FIG. 21 described above, numerical values are entered in each curve so that the level of the magnetic field intensity and the integrated value can be understood. However, on the display, color is classified according to the level of the contour lines, and three-dimensional color display is performed. . At the same time, as shown in FIG. 13, a time waveform (magnetocardiogram) and an electrocardiogram can also be displayed, so that a comprehensive analysis relating to a heart disease can be performed.

【００７３】図２３は本発明の生体磁場計測装置のデス
プレイに表示された処理画像の一例を示す図である。図
２３に於いて、ＭＣＧは心磁図の例、ＱＲＳは積分範囲
をＱＲＳ波の発生する期間Ｔ₁とし（数３４）により得
られた第１の等磁場曲線図、Ｔは積分範囲をＴ波の発生
する期間Ｔ₂とし（数３４）により得られた第２の等磁
場曲線図、（Ｔ−ＱＲＳ）は第１及び第２の等磁場曲線
図の差の各例を示す。図２２、図２３に示すディスプレ
イ上の表示例では、等高線の高低によって色分けをして
３次元カラー表示している。FIG. 23 is a view showing an example of a processed image displayed on the display of the biomagnetic field measuring apparatus of the present invention. In FIG. 23, MCG Examples of MCG, first isomagnetic curve diagram obtained by the QRS is set to a period T ₁ of occurrence of QRS wave integration range (number 34), T is the integration range T wave isomagnetic curve diagram of the time period T ₂ for generating second obtained by (number 34), (T-QRS) indicate each example of the difference between the isomagnetic curves of the first and second. In the display examples on the display shown in FIGS. 22 and 23, three-dimensional color display is performed by color-coding according to the level of contour lines.

【００７４】なお、（数４）、（数３４）に於いて、積
分を行なわず簡便な方法により、Ｉ₁（ｘ、ｙ）、Ｉ
₂（ｘ、ｙ）を求めることもできる。即ち、以下の（数
４１）〜（数４４）からＩ₁（ｘ、ｙ）、Ｉ₂（ｘ、ｙ）
を求めて、更に、（数３５）〜（数４０）を適用する。
生体から発する磁場の接線成分（生体の面に平行な成
分）Ｂ_x（ｘ、ｙ、ｔ）、Ｂ_y（ｘ、ｙ、ｔ）を計測する
場合には（但し、直交座標系（、ｚ）に於いて生体の面
に平行な面をｘｙ面、生体の面に垂直な軸をｚとす
る）、接線成分Ｂ_xとＢ_yの２乗和の平方根から２次元ベ
クトル強度│Ｂ_xy（ｘ、ｙ）│（│ │は絶対値を表わ
す）を（数４１）により求める。In Equations (4) and (34), I ₁ (x, y) and I ₁ (x, y) are calculated by a simple method without integration.
₂ (x, y) can also be obtained. That is, from the following (Formula 41) to (Formula 44), I ₁ (x, y) and I ₂ (x, y)
And (Equation 35) to (Equation 40) are further applied.
When measuring the tangential component of the magnetic field emanating from the living body (a component parallel to the plane of the living _{body) B x (x, y,} t), B y (x, y, t) and is (however, orthogonal coordinate system (, z xy plane parallel to the plane of the living body in), and z axis perpendicular to the plane of the living body), 2-dimensional vector magnitude from the square root of the square sum of tangential components B _x and B _y │B _xy ( x, y) | (|| represents an absolute value) by (Equation 41).

【００７５】 │Ｂ_xy（ｘ、ｙ、ｔ₀）│ ＝√｛（Ｂ_x（ｘ、ｙ、ｔ₀））²＋（Ｂ_y（ｘ、ｙ、ｔ₀））²｝ …（数４１） 次いで、各点（ｘ、ｙ）について任意の時点での波形│
Ｂ_xy（ｘ、ｙ、ｔ₀）│の値Ｉ₁（ｘ、ｙ）を（数１４）
により求め、内挿、外挿により各点（ｘ、ｙ）でのＩ₁
（ｘ、ｙ）が同じ値の点を結ぶ等積分図を求めて、等積
分図を表示画面に表示する。| B _xy (x, y, t ₀ ) | = {(B _x (x, y, t ₀ )) ² + (B _y (x, y, t ₀ )) ² } (Equation 41) ) Next, the waveform at each point (x, y) at any time |
The value I ₁ (x, y) of B _xy (x, y, t ₀ ) |
I ₁ at each point (x, y) by interpolation and extrapolation.
An isointegral diagram connecting points of the same value (x, y) is obtained, and the isointegral diagram is displayed on the display screen.

【００７６】 Ｉ₁（ｘ、ｙ）＝│Ｂ_xy（ｘ、ｙ、ｔ₀）│ …（数４２） 生体の面に垂直な磁場成分Ｂ_z（ｘ、ｙ、ｔ）を計測す
る場合には、垂直な磁場成分Ｂ_z（ｘ、ｙ、ｔ₀）のｘ方
向の変化率∂Ｂ_z（ｘ、ｙ、ｔ₀）／∂ｘと、Ｂ_z（ｘ、
ｙ、ｔ₀）の方向の変化率∂Ｂ_z（ｘ、ｙ、ｔ₀）／∂ｙ
と求め、（数４３）に示すように２乗和の平方根Ｓ
（ｘ、ｙ、ｔ）を求める。I ₁ (x, y) = │B _xy (x, y, t ₀ ) │ (Expression 42) When measuring the magnetic field component B _z (x, y, t) perpendicular to the surface of the living body Is the rate of change 垂直 B _z (x, y, t ₀ ) / ∂x of the perpendicular magnetic field component B _z (x, y, t ₀ ) in the x direction, and B _z (x,
y, t ₀ ) direction of change ∂B _z (x, y, t ₀ ) / ∂y
And the square root S of the sum of squares as shown in (Equation 43)
(X, y, t) is obtained.

【００７７】 Ｓ（ｘ、ｙ、ｔ₀）＝√［｛∂Ｂ_z（ｘ、ｙ、ｔ₀）／∂ｘ｝² ＋｛∂Ｂ_z（ｘ、ｙ、ｔ₀）／∂ｙ｝²］ …（数４３） 次いで、各点（ｘ、ｙ）について任意の時点での波形Ｓ
_t0（ｔ₀、ｘ、ｙ）の値Ｉ₂（ｘ、ｙ）を（数４４）によ
り求め、内挿、外挿により各点（ｘ、ｙ）での値Ｉ
₂（ｘ、ｙ）が同じ値の点を結ぶ等積分図を求めて、等
積分図を表示画面に表示する。S (x, y, t ₀ ) = {[｛∂B _z (x, y, t ₀ ) / ∂x｝ ² + ｛∂B _z (x, y, t ₀ ) / ∂y｝ ² (Equation 43) Next, the waveform S at each point (x, y) at an arbitrary time point
_t0 (t _0, x, y) value I ₂ (x, y) of the determined by (number 44), the value I at each point interpolation, by extrapolation (x, y)
₂ Find an isointegral diagram connecting points with the same value of (x, y) and display the isointegral diagram on the display screen.

【００７８】 Ｉ₂（ｘ、ｙ）＝│Ｓ_t0（ｘ、ｙ、ｔ₀）│ｄｔ …（数４４） なお、（数４１）〜（数４４）に於いてｔ₀として、例
えば、心臓を測定の対象とする時には、心室が収縮した
時のＱ、Ｒ、Ｓの各波の極大値を与える時点をとる。更
に、（数４１）〜（数４４）に於いてｔ₀として、複数
のｔ₀をとり求めた複数の値の間での、等加重を含む和
又は差、比を求める等の演算を行ない、内挿、外挿によ
り演算結果が同じ値の点を結ぶ等積分図を求めて、等積
分図を表示画面に表示する。このような方法によって
も、上記で説明した（数４）、（数３４）を用いる方法
とほぼ同様な結果を得ることができる。[0078] _{I 2 (x, y) =} │S t0 (x, y, t 0) │dt ... ( number 44) As t ₀ In equation (41) through (44), for example, heart Is measured, the time point at which the maximum value of each of the Q, R, and S waves when the ventricle contracts is taken. Further, as (t ₀ ) in (Equation 41) to (Equation 44), calculations such as finding a sum or difference including equal weights and a ratio among a plurality of values obtained by taking a plurality of t ₀ are performed. , Interpolation and extrapolation to obtain an isointegral diagram connecting points having the same value of the operation result, and display the isointegral diagram on the display screen. Even with such a method, it is possible to obtain substantially the same results as the methods using (Equation 4) and (Equation 34) described above.

【００７９】従来方法により法線成分Ｂ_zを測定して得
た患者Ｘの心磁図のＱ波、Ｒ波、Ｓ波の極値が出現する
時点での等磁場線図を、図２４（ａ）、（ｂ）、（ｃ）
に示す。図２４（ａ）、（ｂ）、（ｃ）に於いて、点線
は吸い込まれる磁場の等磁場線図を示し、実線は沸き出
す磁場の等磁場線図を示し、白抜き矢印は電流ダイポー
ルの大きさ、方向を示している。図２４（ａ）、
（ｂ）、（ｃ）に示す等磁場線図には、心臓内に存在す
る電流源を１つと仮定した時の電流ダイポールの位置を
白抜き矢印により示して重ねて表示している。図２４
（ａ）に示すように、Ｑ波の極値が出現する時点では、
心室中隔で右下方向に電流が流れ、図２４（ｂ）に示す
ように、Ｒ波の極値が出現する時点では、左室全体で左
斜め下方向に電流が大きく流れる。また、図２４（ｃ）
に示すように、Ｓ波の極値が出現する時点では、心室基
部方向に右斜め上に電流が流れ、心室の脱分極過程が終
了するのが分かる。FIG. 24A shows an isomagnetic field diagram at the time when the extreme values of the Q, R, and S waves of the magnetocardiogram of the patient X obtained by measuring the normal component B _{z according} to the conventional method. ), (B), (c)
Shown in In FIGS. 24 (a), (b) and (c), the dotted line shows the isomagnetic field map of the magnetic field to be sucked, the solid line shows the isomagnetic field map of the boiling magnetic field, and the white arrow shows the current dipole. The size and direction are shown. FIG. 24 (a),
In the isomagnetic field diagrams shown in (b) and (c), the position of the current dipole assuming that there is one current source in the heart is indicated by a white arrow and superimposed. FIG.
As shown in (a), when the extreme value of the Q wave appears,
A current flows in the lower right direction in the interventricular septum, and as shown in FIG. 24 (b), at the time when the extreme value of the R wave appears, a large current flows diagonally to the lower left in the entire left ventricle. FIG. 24 (c)
As shown in the figure, when the extreme value of the S-wave appears, a current flows obliquely to the upper right in the direction of the base of the ventricle, and the depolarization process of the ventricle ends.

【００８０】上記患者Ｘの心臓から発する磁場の接線成
分Ｂ_x、Ｂ_yを測定し、Ｑ波、Ｒ波、Ｓ波の各極値が出現
する時点に於いて、接線成分を（数４１）、（数４２）
に基づいて合成した等磁場線図を、図２５（ａ）、
（ｂ）、（ｃ）に示す。[0080] the tangential component of the magnetic field generated from the heart of the patient X B _x, the B _y is measured, Q wave, R wave, at the time of the extreme value of the S wave appears, the tangential component (number 41) , (Equation 42)
FIG. 25 (a) shows an isomagnetic field diagram synthesized based on
(B) and (c).

【００８１】図２５（ａ）のパターンと図２４（ａ）の
パターン、図２５（ｂ）のパターンと図２４（ｂ）のパ
ターン、図２５（ｃ）のパターンと図２４（ｃ）のパタ
ーン、はそれぞれほぼ一致する。しかし、図２５（ｂ）
に示すＲ波の極値が出現する時点のパターンでは、心筋
は広い領域で活動しており、図２４（ｂ）のＲ波の極値
が出現する時点のパターンでは鮮明でなかった複数の電
流源が容易に判別でき、電流源の１つは左方向に存在
し、他の電流源は下方に存在することが分かる。The patterns shown in FIGS. 25 (a) and 24 (a), the patterns shown in FIGS. 25 (b) and 24 (b), the patterns shown in FIGS. 25 (c) and 24 (c) , Are almost the same. However, FIG.
In the pattern at the time when the extreme value of the R-wave appears as shown in FIG. 24, the myocardium is active in a wide area, and a plurality of currents that are not clear in the pattern at the time when the extreme value of the R-wave in FIG. The sources can be easily identified, and it can be seen that one of the current sources is to the left and the other current source is below.

【００８２】図２４（ａ）、（ｂ）、（ｃ）に示す、Ｑ
波、Ｒ波、Ｓ波の各極値が出現する時点での法線成分Ｂ
_zの等磁場線図データをそれぞれ用いて、（数４３）、
（数４４）に基づいて求めた、Ｑ波、Ｒ波、Ｓ波の各極
値が出現する時点の等磁場線図を、図２６（ａ）、
（ｂ）、（ｃ）に示す。図２６（ａ）、（ｂ）、（ｃ）
に示す結果から、図２４（ａ）、（ｂ）、（ｃ）に示す
法線成分Ｂ_zの等磁場線図や、（数１）に基づくアロー
マップでは判別しにくかった複数の電流源が判別でき
る。図２６（ａ）、（ｂ）、（ｃ）のパターンは、図２
５（ａ）、（ｂ）、（ｃ）に示すパターン（接線成分Ｂ
_x、Ｂ_y合成から得られるＢ_xyの等磁場線図）と同等であ
ることが分かる。このことは、（数６）及び（数７）、
又は（数２７）及び（数２８）が実際の実験データでほ
ぼ成立していることを意味している。The Q shown in FIGS. 24 (a), (b) and (c)
Normal component B at the time when each extreme value of wave, R wave and S wave appears
Using _z isomagnetic field data of each (number 43),
FIG. 26A shows the isomagnetic field diagram at the time when the extreme values of the Q wave, the R wave, and the S wave appear based on (Equation 44).
(B) and (c). FIGS. 26 (a), (b), (c)
From the results shown in FIG. 24 (a), (b), and isomagnetic field of normal component B _z of (c), the plurality of current source which was difficult to determine the arrow map that is based on equation (1) Can be determined. The patterns in FIGS. 26A, 26B and 26C are shown in FIG.
5 (a), (b) and (c) show patterns (tangent component B
_x, it is found that equivalent to isomagnetic field) of the resulting B _xy from B _y synthesis. This means (Equation 6) and (Equation 7),
Or, it means that (Equation 27) and (Equation 28) are almost satisfied with actual experimental data.

【００８３】なお、図２４（ａ）から図２６（ｃ）の各
図に於いて、横軸ｘ、縦軸ｙは、生体表面に配置された
磁場センサの位置座標を表わす。In each of FIGS. 24A to 26C, the horizontal axis x and the vertical axis y represent the position coordinates of the magnetic field sensor arranged on the living body surface.

【００８４】以上の説明では、心磁場計測に関する例を
とって本発明を説明したが、脳磁図（ＭＥＧ）を得る脳
磁場計測の場合にも本発明が適用できることは言うまで
もない。In the above description, the present invention has been described by taking an example relating to magnetocardiographic measurement. However, it goes without saying that the present invention can also be applied to cerebral magnetic field measurement for obtaining a magnetoencephalogram (MEG).

【００８５】図２７は脳磁場を計測する脳磁場計測シス
テムの脳磁場計測用デュワの内部構成の一部を示す断面
図である。図２７に示すように、脳磁場を計測する場合
には、胸部と異なり頭部は球状であるため、ＳＱＵＩＤ
磁束計１０３−１、１０３−２、…、１０３−Ｎを内蔵
する頭部計測用デュワ１０２の底面の形状を半球として
頭部１００を覆うようにする。ＳＱＵＩＤ磁束計１０３
−１、１０３−２、…、１０３−Ｎは頭部計測用デュワ
１０２の内側の面に沿って放射状に配置され、各ＳＱＵ
ＩＤ磁束計の先端面（磁場計測面）は半球面の接線面に
ほぼ平行となるように配置されている。半球の中心が頭
部の脳部のほぼ中心と一致するように脳部を球と仮定し
て半球の半径は設定され、この半径は成人でも測定でき
るよう約１０ｃｍとした。頭部計測用デュワ１０２の内
部には熱輻射シールド部材１０４が配置され頭部計測用
デュワの上部は上板１０５により密閉されている。ＳＱ
ＵＩＤ磁束計１０３−１、…、１０３−Ｎにより検出さ
れた信号は信号線１０６−１、…、１０６−Ｎを通して
頭部計測用デュワの外部に取り出される。FIG. 27 is a sectional view showing a part of the internal configuration of a dewar for brain magnetic field measurement of a brain magnetic field measurement system for measuring a brain magnetic field. As shown in FIG. 27, when the brain magnetic field is measured, since the head is spherical unlike the chest, the SQUID
.., 103-N, the bottom surface of the head measurement dewar 102 is formed as a hemisphere so as to cover the head 100. SQUID magnetometer 103
, 103-2,..., 103-N are radially arranged along the inner surface of the dewar 102 for head measurement.
The tip surface (magnetic field measurement surface) of the ID magnetometer is disposed so as to be substantially parallel to the tangent surface of the hemisphere. The radius of the hemisphere was set assuming that the brain is a sphere so that the center of the hemisphere substantially coincides with the center of the brain of the head, and this radius was set to about 10 cm so that even adults could measure it. A heat radiation shield member 104 is disposed inside the head measurement dewar 102, and an upper portion of the head measurement dewar is sealed by an upper plate 105. SQ
The signals detected by the UID magnetometers 103-1,..., 103-N are taken out of the head measuring dewar through the signal lines 106-1,.

【００８６】図２８は図２７に示す脳磁場計測システム
により計測可能な磁場成分と頭部の関係を説明する図で
ある。頭部の上方に放射状に複数の位置の１つＯ’配置
されたＱＵＩＤ磁束計により計測可能な脳磁場Ｂの成分
は、Ｏを原点とする極座標（ｒ、θ、φ）に於けるｒ方
向の成分Ｂ_r（法線成分）である。図２８に於いて、成
分Ｂ_θ、Ｂ_φは頭部表面に平行な接線成分を示し、原点
Ｏは脳部を球と仮定した時の球の中心である。体性感覚
として右中指に電気刺激を与え、図２７に示す脳磁場計
測システムにより法線成分Ｂ_rを検出し、電気刺激を与
えてから約１００ｍｓｅｃ後に出現する脳波が最大とな
る時点での等磁場曲線図を求める。図２９（ａ）、
（ｂ）は、図２７に示す脳磁場計測システムにより得ら
れる等磁場曲線図の一例を示す図であり、図２９（ａ）
は従来の方法による法線成分Ｂ_rの等磁場曲線図、図２
９（ｂ）は以下に示す本発明の（数４５）を使用して得
られる等磁場曲線図（地球儀に示された地図の如く、脳
部を近似する球面に表示された脳磁場の強度分布を示
す。）を示す。FIG. 28 is a view for explaining the relationship between a magnetic field component and the head which can be measured by the brain magnetic field measurement system shown in FIG. The component of the cerebral magnetic field B that can be measured by the QUID magnetometer radially arranged at one of a plurality of positions O ′ above the head is the r direction in polar coordinates (r, θ, φ) with O as the origin. Is a component _Br (normal component) of In FIG. 28, components B _θ and B _φ indicate tangential components parallel to the head surface, and the origin O is the center of the sphere when the brain is assumed to be a sphere. An electrical stimulus is given to the right middle finger as a somatic sensation, a normal component _Br is detected by the cerebral magnetic field measurement system shown in FIG. Obtain a magnetic field curve diagram. FIG. 29 (a),
FIG. 29B is a diagram showing an example of an isomagnetic field curve diagram obtained by the brain magnetic field measurement system shown in FIG. 27, and FIG.
Isomagnetic curves of normal component B _r by the conventional method, FIG. 2
9 (b) is an isomagnetic field curve diagram obtained by using the following (Equation 45) of the present invention (the intensity distribution of the brain magnetic field displayed on a spherical surface approximating the brain as in a map shown on a globe) Is shown.)

【００８７】 Ｓ（θ、φ、ｔ） ＝√｛（∂Ｂ_r（ｔ）／∂θ）²＋（∂Ｂ_r（ｔ）／∂φ）²｝ …（数４５） 図２９（ａ）に示す等磁場線図には、脳内に存在する電
流源を１つと仮定した時の電流ダイポールの位置を白抜
き矢印により示して重ねて表示している。図２９（ａ）
において、点線は吸い込まれる磁場の等磁場線図を示
し、実線は沸き出す磁場の等磁場線図を示し、白抜き矢
印は電流ダイポールの大きさ、方向を示している。図２
９（ａ）に示す法線成分Ｂ_rの等磁場曲線図で従来推定
していた電流源（白抜き矢印で示す電流ダイポール）
が、図２９（ｂ）に示す等磁場曲線図ではピーク位置Ａ
に対応して出現していることが容易に直視できる。な
お、図２７に図示しない脳磁場計測システムのその他の
構成は基本的に図７に示す生体磁場計測装置の構成と同
一である。S (θ, φ, t) = √ ｛(∂B _r (t) / ∂θ) ² + (∂B _r (t) / ∂φ) ² … (Formula 45) FIG. 29 (a) In the isomagnetic field diagram shown in FIG. 5, the position of the current dipole when the number of current sources existing in the brain is assumed to be one is indicated by a white arrow and superimposed. FIG. 29 (a)
In FIG. 7, a dotted line shows an iso-magnetic field diagram of a magnetic field to be absorbed, a solid line shows an iso-magnetic field diagram of a boiling magnetic field, and outline arrows show the size and direction of a current dipole. FIG.
9 the current source has been estimated conventionally with an equal magnetic field curves of normal component B _r shown in (a) (current dipole shown by a white arrow)
However, in the isomagnetic field curve diagram shown in FIG.
Can easily be seen directly. The other configuration of the brain magnetic field measurement system not shown in FIG. 27 is basically the same as the configuration of the biomagnetic field measurement device shown in FIG.

【００８８】以上説明した本発明による各種の方法によ
り得られる心磁場、脳磁場に関する等磁場線図を使っ
て、磁場源を解析する方法として、逆問題を解く様々の
アルゴリズムが考えられる。実際に多く使用されている
単純なアルゴリズムは、磁場源として単一あるいは２つ
程度の電流ダイポールを想定し、これら電流ダイポール
が存在する位置座標を初期条件として任意に仮定して、
各位置座標に存在する電流ダイポールが、ビオサバール
の式で表される磁場を作るものとして、実測した磁場の
計測点（ｘ、ｙ）での磁場を計算する。計算された磁場
〈Ｂ_c（ｘ、ｙ）〉と実測値の磁場〈Ｖ_m（ｘ、ｙ）〉
（ｍ＝１、２、…、Ｍ：Ｍは実測される磁場の計測点の
総数）との差で表される次の（数１７）に示す評価関数
を計算し、各電流ダイポールの位置座標を変化させて、
評価関数Ｌの最小値を解析的に求めていく。（数４６）
に於いて、Ｇは定数、〈ｎ_s〉は法線又はｚ方向の単位
ベクトルであり、加算記号Σは、ｍ＝１、２、…、Ｍに
関する加算を示す。Various algorithms for solving the inverse problem can be considered as a method of analyzing the magnetic field source using the isomagnetic maps relating to the cardiac magnetic field and the brain magnetic field obtained by the various methods according to the present invention described above. A simple algorithm that is often used in practice is to assume one or two current dipoles as the magnetic field source, and arbitrarily assuming the position coordinates where these current dipoles exist as initial conditions.
The magnetic field at the measurement point (x, y) of the actually measured magnetic field is calculated on the assumption that the current dipole existing at each position coordinate generates a magnetic field represented by the Biot-Savart equation. The calculated magnetic field <B _c (x, y)> and the actually measured magnetic field <V _m (x, y)>
(M = 1, 2,..., M: where M is the total number of actually measured magnetic field measurement points) and the following evaluation function expressed by (Expression 17) is calculated, and the position coordinates of each current dipole are calculated. To change
The minimum value of the evaluation function L is obtained analytically. (Equation 46)
, G is a constant, < _ns > is a normal vector or a unit vector in the z direction, and the addition symbol Σ indicates addition for m = 1, 2,..., M.

【００８９】 Ｌ＝Σ｛〈Ｖ_m（ｘ、ｙ）〉−Ｇ（［〈Ｂ_c（ｘ、ｙ）〉］・ｎ_s）｝² …（数４６） しかし、（数４６）に基づく方法では、磁場の広い測定
領域を解析する場合、最小値に収束しない場合も出てく
る。本発明では、評価関数Ｌを算出の際のダイポールの
位置と個数の初期条件を、（数４）、（数３４）、又は
（数４４）に基づく等磁場線図に於けるピーク位置をダ
イポールの位置とし、更に、等磁場線図に於けるピーク
の個数をダイポールの個数として予め決める。このよう
に初期条件を与え評価関数Ｌを解くことにより、磁場源
解析が必ず収束する。ディスプレイ上に表示される、
（数４）、（数３４）、又は（数４４）に基づく心磁
場、脳磁場に関する等磁場線図上での各ピーク位置を指
定することにより、自動的に各ピーク位置の座標とその
個数が上記の初期値として自動的に装置に入力され、評
価関数Ｌが解かれ、収束する磁場源解析結果が得られ
る。L = {<V _m (x, y)> − G ([<B _c (x, y)>] · n _s )} ² (Equation 46) However, the method based on (Equation 46) Then, when analyzing a measurement region with a wide magnetic field, there may be cases where the measurement value does not converge to the minimum value. In the present invention, the initial condition of the position and the number of dipoles at the time of calculating the evaluation function L is defined as the peak position in the isomagnetic field map based on (Equation 4), (Equation 34), or (Equation 44). And the number of peaks in the isomagnetic field diagram is determined in advance as the number of dipoles. By thus giving the initial conditions and solving the evaluation function L, the magnetic field source analysis always converges. Displayed on the display,
By specifying each peak position on the isomagnetic field map relating to the cardiac magnetic field and brain magnetic field based on (Equation 4), (Equation 34), or (Equation 44), the coordinates of each peak position and the number thereof are automatically set. Is automatically input to the apparatus as the above initial value, the evaluation function L is solved, and a converged magnetic field source analysis result is obtained.

【００９０】従って、従来技術のように試行錯誤的に初
期値を設定するのではなく、計測の結果得られる等磁場
線図のデータに基づいて、初期値設定をほぼ一義的にか
つ容易に可能ででき、効率よくより正確に逆問題を解く
ことが可能となる。Therefore, the initial value can be set almost unambiguously and easily based on the data of the isomagnetic field diagram obtained as a result of the measurement, instead of setting the initial value by trial and error unlike the prior art. It is possible to efficiently and more accurately solve the inverse problem.

【００９１】なお、以上の説明に於いて使用した等磁場
線図を表わす各図では、医療の分野で行なわれている通
例に従い、人体の右側を各図の左側に表示し、人体の左
側を各図の右側に表示している。In each of the figures showing the isomagnetic field map used in the above description, the right side of the human body is displayed on the left side of each figure and the left side of the human body is displayed in accordance with the usual practice in the medical field. It is displayed on the right side of each figure.

【００９２】[0092]

【発明の効果】本発明では、ベクトル計測により接線成
分Ｂ_x、Ｂ_yを計測することなく、法線成分Ｂ_zの計測の
みから、従来技術に於けるＢ_xyに基づく等磁場線図と等
価的な等磁場線図が得られる。従来技術の於ける法線成
分Ｂ_zから直接得る等磁場線図では、複数の電流源は判
別しにくかったが、本発明の等磁場線図では、従来技術
に於けるＢ_xyに基づく等磁場線図と同様に、電流源の直
上にピークパターンが出現するので、生体内の複数の電
流源を直読でき、複数の電流源の位置、大きさ等を解析
する逆問題が容易に解けるようになる。本発明の装置に
よれば、心筋梗塞、虚血等の発見、不整脈を生じている
位置の発見、心筋肥大の発見、術前術後の心筋状態の変
化の評価等の心臓に関する疾患の発見、状態の確認等が
容易にできる。In the present invention, the tangential component B _x by vector measurement, without measuring the B _y, only the measurement of the normal component B _z, such a magnetic field diagrams equivalent based on in B _xy in the prior art Is obtained. The isomagnetic field obtained directly from in normal component B _z of the prior art, a plurality of current sources is was difficult to determine, etc. In the field diagram, such a magnetic field based on at B _xy in the prior art of the present invention As in the diagram, a peak pattern appears directly above the current source, so that multiple current sources in the living body can be read directly, and the inverse problem of analyzing the position, size, etc. of multiple current sources can be easily solved. Become. According to the device of the present invention, myocardial infarction, discovery of ischemia, etc., discovery of the position where arrhythmia has occurred, discovery of myocardial hypertrophy, discovery of heart-related diseases such as evaluation of changes in myocardial condition before and after surgery, It is easy to check the status.

【図面の簡単な説明】[Brief description of the drawings]

【図１】本発明に於いて、心磁場の発生を、無限平面導
体中の電流ダイポールから発生する磁場によりモデル化
して解析するための図。FIG. 1 is a diagram for modeling and analyzing generation of a cardiac magnetic field using a magnetic field generated from a current dipole in an infinite planar conductor in the present invention.

【図２】本発明に於いて、無限平面導体の内部に存在す
る電流ダイポールのモーメントの概略位置を示す図。FIG. 2 is a diagram showing a schematic position of a moment of a current dipole existing inside an infinite planar conductor in the present invention.

【図３】本発明に於いて、無限平面導体の上でのＢ_x及
び−∂Ｂ_Z／∂ｘをそれぞれの最大値で規格化した相対
磁場強度曲線Ｃ₁、Ｃ₂を示す図。FIG. 3 is a diagram showing relative magnetic field strength curves C ₁ and C _{2 in} which B _x and −ΔB _Z / Δx on an infinite plane conductor are normalized by their respective maximum values in the present invention.

【図４】本発明に於いて、−∂Ｂ_Z（ｘ、０）／∂ｘの
第１項、第２項、第３項を示す磁場強度曲線Ｃ₃、Ｃ₄、
Ｃ₅を示す図。FIG. 4 is a graph showing magnetic field intensity curves C ₃ , C ₄ , C ₁ , C ₂ , and C _{3 of} −∂B _Z (x, 0) / ∂x according to the present invention.
Shows the C _5.

【図５】本発明に於いて、Ｂ_x、∂Ｂ_Z／∂ｘのそれぞれ
の第１項と第２項を規格化の後に比較した相対磁場強度
曲線Ｃ₆、Ｃ₇、Ｃ₈、Ｃ₉を示す図。FIG. 5 shows relative magnetic field intensity curves C ₆ , C ₇ , C ₈ , C ₈ in which the first and second terms of B _x , ∂B _Z / ∂x are compared after normalization. shows a _9.

【図６】本発明に於いて、α＝（∇Ｋ）_z／Ｋ、｛−∂
Ｂ_Z（ｘ、０）／∂ｘの第１項｝／｛Ｂ_x（ｘ、０）の第
１項｝、｛−∂Ｂ_Z（ｘ、０）／∂ｘの第２項｝／｛Ｂ_x
（ｘ、０）の第２項｝の各々の磁場強度曲線Ｃ₁₀、
Ｃ₁₁、Ｃ₁₂を示す図。FIG. 6: In the present invention, α = (∇K) _z / K, ｛−∂
B _Z (x, 0) / the first term of ∂x} / {first term of _{B x (x, 0)}} , {- ∂B Z (x, 0) / second term of ∂x} / { B _x
The magnetic field strength curve C ₁₀ of each of the second term の of (x, 0),
Shows the C _11, C _12.

【図７】本発明が実施される心磁場計測を行なう生体磁
場計測装置の概略構成を示す図。FIG. 7 is a diagram showing a schematic configuration of a biomagnetic field measuring apparatus for performing a cardiac magnetic field measurement according to the present invention.

【図８】本発明が実施される心磁場計測を行なう生体磁
場計測装置に於ける磁場センサの配置構成を示す図。FIG. 8 is a diagram showing an arrangement configuration of a magnetic field sensor in a biomagnetic field measuring apparatus that performs a cardiac magnetic field measurement according to the present invention.

【図９】本発明が実施される心磁場計測を行なう生体磁
場計測装置に於ける磁場の法線成分を検出する磁場セン
サ単体の構成を示す図。FIG. 9 is a diagram showing the configuration of a single magnetic field sensor for detecting a normal component of a magnetic field in a biomagnetic field measuring apparatus for performing a cardiac magnetic field measurement according to the present invention.

【図１０】本発明が実施される心磁場計測を行なう生体
磁場計測装置に於ける磁場の接線成分を検出する磁場セ
ンサ単体の構成を示す図。FIG. 10 is a diagram showing the configuration of a single magnetic field sensor for detecting a tangential component of a magnetic field in a biomagnetic field measurement apparatus for performing a cardiac magnetic field measurement according to the present invention.

【図１１】本発明が実施される心磁場計測を行なう生体
磁場計測装置に於ける磁場センサの配置と人体の胸部と
の位置関係を示す図。FIG. 11 is a diagram showing a positional relationship between an arrangement of a magnetic field sensor and a chest of a human body in a biomagnetic field measuring apparatus for performing a cardiac magnetic field measurement according to the present invention.

【図１２】本発明の実施例に於いて、各磁場センサ位置
に於いて計測した健常者の心臓から発する磁場の各方向
の成分に於ける時間波形図。FIG. 12 is a time waveform chart in each direction component of a magnetic field emitted from a healthy person's heart measured at each magnetic field sensor position in the embodiment of the present invention.

【図１３】本発明の実施例に於いて、健常者について計
測された特定の２チャンネルに関する接線成分（Ｂ_x）
の時間波形を示す図。FIG. 13 is a diagram illustrating a tangential component (B _x ) of two specific channels measured for a healthy person in an embodiment of the present invention.
FIG.

【図１４】本発明の実施例に於いて、磁場の接線成分Ｂ
_x、Ｂ_yを計測した健常者の心磁波形から得た、Ｑ波、Ｒ
波、Ｓ波の各波のピーク時に於ける等磁場線図。FIG. 14 shows a tangential component B of a magnetic field in the embodiment of the present invention.
_x, to obtain a B _y from healthy subjects of the magnetocardiogram waveform measured, Q wave, R
FIG. 3 is an isomagnetic field diagram at the peak of each of the wave and the S wave.

【図１５】本発明の実施例に於いて、健常者の心磁波形
のＱＲＳ波が出現する時間帯に於いて検出された２つの
接線成分から得た等積分図。FIG. 15 is an isometric diagram obtained from two tangent components detected in a time zone in which a QRS wave of a magnetocardiogram of a healthy person appears in the example of the present invention.

【図１６】本発明の実施例に於いて、健常者の心磁波形
のＱＲＳ波が出現する時間帯に於いて検出された法線線
分から得た等積分図。FIG. 16 is an isointegral graph obtained from a normal line segment detected in a time zone in which a QRS wave of a magnetocardiogram of a healthy person appears in the example of the present invention.

【図１７】本発明の実施例に於いて、健常者の心磁波形
のＴ波が出現する時間帯に於いて検出された２つの接線
成分から得た等積分図。FIG. 17 is an isointegral diagram obtained from two tangent components detected in a time zone in which a T wave of a magnetocardiogram of a healthy person appears in the example of the present invention.

【図１８】図１７に示す等積分図から図１５に示す等積
分図を差し引いた図。FIG. 18 is a diagram obtained by subtracting the isometric diagram shown in FIG. 15 from the isometric diagram shown in FIG. 17;

【図１９】本発明の実施例に於いて、心筋梗塞の患者の
心磁波形のＱＲＳ波が出現する時間帯に於いて検出され
た２つの接線成分から得た等積分図。FIG. 19 is an isointegral diagram obtained from two tangent components detected in a time zone during which a QRS wave of a magnetocardiographic waveform of a patient with myocardial infarction appears in an embodiment of the present invention.

【図２０】本発明の実施例に於いて、心筋梗塞の患者の
心磁波形のＴ波が出現する時間帯に於いて検出された２
つの接線成分から得た等積分図。FIG. 20 is a graph showing the detection of the T-wave of the magnetocardiographic waveform of a patient with myocardial infarction in the time slot in the embodiment of the present invention.
Isointegral diagram obtained from two tangent components.

【図２１】図２０に示す等積分図から図１９に示す等積
分図を差し引いた図。FIG. 21 is a diagram obtained by subtracting the isometric diagram shown in FIG. 19 from the isometric diagram shown in FIG. 20;

【図２２】本発明が実施される心磁場計測を行なう生体
磁場計測装置のパソコンでの出力画面の例を示す図。FIG. 22 is a diagram showing an example of an output screen on a personal computer of a biomagnetic field measurement apparatus for performing cardiac magnetic field measurement according to the present invention.

【図２３】本発明の生体磁場計測装置のデスプレイに表
示された処理画像の一例を示す図。FIG. 23 is a view showing an example of a processed image displayed on the display of the biomagnetic field measuring apparatus of the present invention.

【図２４】従来方法により法線成分Ｂ_zを測定して得
た、心磁図（ＭＣＧ）のＱ波、Ｒ波、Ｓ波の極値が出現
する時点での等磁場線図を示す図。FIG. 24 is a diagram showing an isomagnetic field diagram at the time when the extreme values of the Q wave, R wave, and S wave of the magnetocardiogram (MCG) appear, obtained by measuring the normal component B _z by the conventional method.

【図２５】本発明の実施例に於いてそれぞれ、心臓から
の磁場の接線成分Ｂ_x、Ｂ_yを測定し、Ｑ波、Ｒ波、Ｓ波
の極値が出現する時点に於いて、接線成分を合成したＢ
_xyの等磁場線図を示す図。[Figure 25], respectively In the embodiment of the present invention, the tangential component B _x of a magnetic field from the heart, the B _y is measured, at the time of the Q wave, R wave, the extreme values of the S-wave occurrence, tangent B that synthesized the components
_The figure which shows the _isomagnetic field map of _xy .

【図２６】本発明の実施例に於いて、図２４に示す、Ｑ
波、Ｒ波、Ｓ波の極値が出現する時点での法線Ｂ_zの等
磁場線図データをそれぞれ用いて、（数４３）、（数４
４）に基づいて求めた、各時点での等磁場線図を示す
図。FIG. 26 is a graph showing the operation of Q shown in FIG. 24 in the embodiment of the present invention;
Waves, using R-wave, the isomagnetic field data of the normal B _z at the time the extreme values of the S wave appears respectively, (number 43), (Expression 4
FIG. 4 is a view showing isomagnetic field diagrams at each time point obtained based on 4).

【図２７】本発明の実施例に於いて、脳磁場を計測する
脳磁場計測システムの脳磁場計測用デュワの内部構成の
一部を示す断面図。FIG. 27 is a cross-sectional view showing a part of an internal configuration of a dewar for measuring a cerebral magnetic field in a cerebral magnetic field measuring system for measuring a cerebral magnetic field in an embodiment of the present invention.

【図２８】図２７に示す脳磁場計測システムにより計測
可能な磁場成分と頭部の関係を説明する図。28 is a view for explaining the relationship between a magnetic field component and the head which can be measured by the brain magnetic field measurement system shown in FIG. 27.

【図２９】図２７に示す脳磁場計測システムにより得ら
れる等磁場曲線図の一例を示す図。FIG. 29 is a diagram showing an example of an isomagnetic field curve diagram obtained by the brain magnetic field measurement system shown in FIG. 27;

【符号の説明】[Explanation of symbols]

１…磁場シールドルーム、２…被検者、３…ベッド、４
…デュワ、５…自動補給装置、６…ＦＦＬ回路、７…フ
イルター回路、８…計算機、１０、１０’、１０”…検
出コイル、１１、１１’、１１”…参照コイル、１２、
１２’、１２”…ＳＱＵＩＤ、１３…ｘ成分検出用磁場
センサ、１４…ｙ成分検出用磁場センサ、２０−１、２
０−２、〜、２０−８、２１−１、〜、２１−８、２２
−１、〜、２２−８、２３−２、〜、２３−８、２４−
１、〜、２４−８、２５−１、〜、２５−８、２６−
１、〜、２６−８、２７−１、〜、２７−８…磁場セン
サ、３０…胸部、１０３−１、１０３−２、…、１０３
−Ｎ…ＳＱＵＩＤ磁束計、１００…頭部、１０２…頭部
計測用デュワ、１０４…熱輻射シールド部材、１０５…
上板、１０６−１、…、１０６−Ｎ…信号線。1 ... magnetic shield room, 2 ... subject, 3 ... bed, 4
... Dua, 5 ... Automatic replenishing device, 6 ... FFL circuit, 7 ... Filter circuit, 8 ... Computer, 10, 10 ', 10 "... Detection coil, 11, 11', 11" ... Reference coil, 12,
12 ', 12 "... SQUID, 13 ... Magnetic field sensor for detecting x component, 14 ... Magnetic field sensor for detecting y component, 20-1, 2
0-2, ~, 20-8, 21-1, ~, 21-8, 22
-1, ~ 22-8, 23-2, ~ 23-8, 24-
1, to 24-8, 25-1, to 25-8, 26-
1, ..., 26-8, 27-1, ..., 27-8 ... magnetic field sensor, 30 ... chest, 103-1, 103-2, ..., 103
-N: SQUID magnetometer, 100: head, 102: head dewar, 104: heat radiation shield member, 105:
.., 106-N... Signal lines.

───────────────────────────────────────────────────── フロントページの続き (72)発明者 鈴木 博之 茨城県ひたちなか市市毛882番地 株式 会社日立製作所計測器事業部内 (72)発明者 近藤 昭二 茨城県ひたちなか市市毛882番地 株式 会社日立製作所計測器事業部内 (72)発明者 小見山 泰明 茨城県ひたちなか市市毛882番地 株式 会社日立製作所計測器事業部内 (72)発明者 岡島 健一 東京都国分寺市東恋ケ窪一丁目280番地 株式会社日立製作所中央研究所内 (56)参考文献 特開 平２−249530（ＪＰ，Ａ) 特開 平10−248821（ＪＰ，Ａ) ”Ｔａｎｇｅｎｔｉａｌ（ｘ− ａｎ ｄ ｙ−）Ｃｏｍｐｏｎｅｎｔ Ｍｅａ ｓｕｒｅｍｅｎｔ ｏｆ Ｃａｒｄｉａ ｃ Ｍａｇｎｅｔｉｃ Ｆｉｅｌｄ ａ ｎｄ Ｃｏｍｐａｒｉｓｏｎ ｗｉｔｈ Ｃｏｎｖｅｎｔｉｏｎａｌ ｚ−Ｃｏ ｍｐｏｎｅｎｔ Ｍｅａｓｕｒｅｍｅｎ ｔ" Ｐｒｅｓｅｎｔｅｄ ａｔ ｔｈｅ 10ｔｈ Ｉｎｔｅｒｎａｔｉｏｎａｌ Ｃｏｎｆｅｒｅｎｃｅ ｏｎ Ｂｉｏｍ ａｇｎｅｔｉｓｍ，ＳａｎｔａＦｅ，Ｎ ｅｗＭｅｘｉｃｏ，Ｆｅｂｒｕａｒｙ 20，1996 (58)調査した分野(Int.Cl.⁷，ＤＢ名) A61B 5/05 ──────────────────────────────────────────────────の Continuing on the front page (72) Inventor Hiroyuki Suzuki 882 Ma, Hitachinaka-shi, Ibaraki Pref.Measurement Division, Hitachi, Ltd. (72) Inventor Shoji Kondo 882 Moo, Hitachinaka-shi, Ibaraki Pref. (72) Inventor Yasuaki Omiyama 882 Ma, Hitachinaka-shi, Ibaraki Pref.Hitachi, Ltd.Measurement Instruments Division (72) Inventor Kenichi Okajima 1-280, Higashi-Koikekubo, Kokubunji-shi, Tokyo Within the Hitachi, Ltd. 56) References JP-A-2-249530 (JP, A) JP-A-10-248821 (JP, A) "Tentential (x-andy-) Component Measuring of Cardiac Magnetic Field and Comparator" rison with Conventional z-Co mponent Measuremen t "Presented at the 10th International Conference on Biom agnetism, SantaFe, N ewMexico, February 20,1996 (58) investigated the field (Int.Cl. ^7, DB name) A61B 5/05

Claims (17)

(57)【特許請求の範囲】(57) [Claims] 【請求項１】生体の外部に配置され，前記生体から発す
る生体磁場を検出する量子干渉素子（ＳＱＵＩＤ）から
なり，前記生体磁場の前記生体の面に垂直な第１方向の
磁場成分の時間変化を検出する複数の磁束計と，前記第
１方向と交叉する第２方向及び第３方向に於ける前記第
１方向の磁場成分の変化率の２乗和の平方根に比例する
値の時間変化を求める演算と，前記値の時間変化を所定
の期間で積分し積分値を求める演算とを行なう演算処理
手段と，前記積分値を表示する表示手段とを有すること
を特徴とする生体磁場計測装置。1. A time-dependent change in a magnetic field component of a biomagnetic field in a first direction perpendicular to a surface of a living body, the quantum interference element being arranged outside the living body and detecting a biomagnetic field emitted from the living body. And a time change of a value proportional to the square root of the sum of squares of the change rate of the magnetic field component in the first direction in the second direction and the third direction crossing the first direction. A biomagnetic field measuring apparatus, comprising: an arithmetic processing unit for performing an operation to obtain the value, an operation of integrating the time change of the value over a predetermined period to obtain an integrated value, and a display unit for displaying the integrated value. 【請求項２】請求項１に記載の生体磁場計測装置に於い
て，前記表示手段に，前記積分値の等しい点を結ぶ等積
分図が表示されることを特徴とする生体磁場計測装置。2. The biomagnetic field measuring apparatus according to claim 1, wherein an isointegral diagram connecting points having the same integral value is displayed on the display means. 【請求項３】請求項１に記載の生体磁場計測装置に於い
て，前記演算処理手段は，前記積分値を求める演算を，
複数の前記所定の期間で行ない前記積分値を複数個求
め，該複数個の前記積分値の間での，比，等加重を含む
和または差の何れかを求める演算を行なうことを特徴と
する生体磁場計測装置。3. The biomagnetic field measuring apparatus according to claim 1, wherein the arithmetic processing means performs an arithmetic operation for calculating the integral value.
A plurality of the integrated values are obtained in the plurality of predetermined periods, and a calculation is performed to obtain any of a sum or a difference including a ratio and an equal weight between the plurality of the integrated values. Biomagnetic field measurement device. 【請求項４】請求項１に記載の生体磁場計測装置に於い
て，前記複数の磁束計が等間隔に配置されることを特徴
とする生体磁場計測装置。4. A biomagnetic field measuring apparatus according to claim 1, wherein said plurality of magnetometers are arranged at equal intervals . 【請求項５】生体の外部に配置され，前記生体から発す
る生体磁場を検出する量子干渉素子（ＳＱＵＩＤ）から
なり，前記生体磁場の前記生体の面に平行な第１，第２
方向の磁場成分の時間変化を検出する複数の磁束計と，
前記第１，第２方向の磁場成分の２乗和の平方根に比例
する値の時間変化を求める演算と，前記値の時間変化を
所定の期間で積分し積分値を求める演算とを行なう演算
処理手段と，前記積分値を表示する表示手段とを有する
ことを特徴とする生体磁場計測装置。5. A first and second quantum interference device (SQUID) disposed outside a living body for detecting a biomagnetic field emitted from the living body, the first and second quantum interference devices being parallel to a surface of the living body of the biomagnetic field.
Multiple magnetometers that detect the time change of the magnetic field component in the direction
Arithmetic processing for calculating the time change of a value proportional to the square root of the sum of squares of the magnetic field components in the first and second directions, and calculating the integrated value by integrating the time change of the value over a predetermined period. biomagnetic field measuring apparatus, characterized by chromatic means, and display means for displaying said integrated value. 【請求項６】請求項５に記載の生体磁場計測装置に於い
て，前記表示手段に，前記積分値の等しい点を結ぶ等積
分図が表示されることを特徴とする生体磁場計測装置。6. The biomagnetic field measuring apparatus according to claim 5, wherein an isointegral diagram connecting points having the same integral value is displayed on the display means. 【請求項７】請求項５に記載の生体磁場計測装置に於い
て，前記演算処理手段は，前記積分値を求める演算を，
複数の前記所定の期間で行ない前記積分値を複数個求
め，該複数個の前記積分値の間での，比，等加重を含む
和または差の何れかを求める演算を行なうことを特徴と
する生体磁場計測装置。7. The biomagnetic field measuring apparatus according to claim 5, wherein the arithmetic processing means performs an arithmetic operation for obtaining the integral value.
A plurality of the integrated values are obtained in the plurality of predetermined periods, and a calculation is performed to obtain any of a sum or a difference including a ratio and an equal weight between the plurality of the integrated values. Biomagnetic field measurement device. 【請求項８】請求項５に記載の生体磁場計測装置に於い
て，前記複数の磁束計が等間隔に配置されることを特徴
とする生体磁場計測装置。8. The biomagnetic field measuring apparatus according to claim 5, wherein said plurality of magnetometers are arranged at equal intervals . 【請求項９】生体表面に平行な面を直交座標のｘ，ｙ平
面，前記生体表面に垂直な方向を前記直交座標のｚ軸と
し，前記生体の外部に配置され，前記生体から発する生
体磁場を検出する量子干渉素子（ＳＱＵＩＤ）からな
り，前記生体表面に垂直な磁場成分Ｂ_z（ｘ，ｙ）を検
出する複数の磁束計と，前記磁場成分Ｂ_z（ｘ，ｙ）に
基づいて，ｘ方向の磁場成分Ｂ _x （ｘ，ｙ）及びｙ方向
の磁場成分Ｂ _y （ｘ，ｙ）を推定する演算処理手段と，
前記ｘ方向の磁場成分Ｂ _x （ｘ，ｙ）及び前記ｙ方向の
磁場成分Ｂ _y （ｘ，ｙ）から求めた任意の時点での等磁
場線図を表示する表示手段とを有することを特徴とする
生体磁場計測装置。9. A biomagnetic field disposed outside the living body and emitted from the living body, wherein a plane parallel to the living body surface is defined as an x, y plane of orthogonal coordinates, and a direction perpendicular to the living body surface is defined as a z-axis of the rectangular coordinates. consists quantum interference device for detecting the (SQUID), the magnetic field component perpendicular to the biological surface B _z (x, y) a plurality of magnetometers for detecting the magnetic field component B _z (x, y) to
Based on the magnetic field component B _x (x, y) in the x direction and the y direction
Arithmetic processing means for estimating the magnetic field component B _y (x, y) of
The x-direction magnetic field component B _x (x, y) and the y-direction
Equimagnetism at any point determined from magnetic field component B _y (x, y)
A biomagnetic field measuring apparatus, comprising: display means for displaying a field map . 【請求項１０】請求項９に記載の生体磁場計測装置に於
いて，演算処理手段は，前記ｘ方向の磁場成分Ｂ
_x （ｘ，ｙ）を前記磁場成分Ｂ _z （ｘ，ｙ）のｘ方向での
変化率から，前記ｙ方向の磁場成分Ｂ _y （ｘ，ｙ）を前
記磁場成分Ｂ _z （ｘ，ｙ）のｙ方向での変化率から推定
し，前記ｘ方向での変化率と前記ｙ方向での変化率の２
乗和の平方根に比例する値を求め，前記任意の時点での
前記値の等しい点を結ぶ前記等磁場線図が，前記表示手
段に表示されることを特徴とする生体磁場計測装置。10. The biomagnetic field measuring apparatus according to claim 9,
AndThe arithmetic processing means calculates the magnetic field component B in the x direction.
_x (X, y) is converted to the magnetic field component B _z (X, y) in x direction
From the rate of change, the magnetic field component B in the y direction _y Before (x, y)
Magnetic field component B _z Estimated from the change rate of (x, y) in the y direction
And the rate of change in the x direction and the rate of change in the y direction are 2
Find a value proportional to the square root of the sum of squares,
The isomagnetic field diagram connecting the points having the same value is displayed on the display means.
Displayed in columnsA biomagnetic field measurement apparatus, characterized in that: 【請求項１１】請求項１０に記載の生体磁場計測装置に
於いて，複数の前記任意の時点をとり求めた複数の前記
値の間での，比，等加重を含む和または差の何れかを求
める演算を行なうことを特徴とする生体磁場計測装置。11. The biomagnetic field measuring apparatus according to claim 10,
A plurality of said plurality of said arbitrary time points
Finds the sum or difference between values, including ratios, equal weights
A biomagnetic field measurement apparatus characterized by performing a calculation to obtain 【請求項１２】生体の外部に配置され，前記生体から発
する生体磁場を検出する量子干渉素子（ＳＱＵＩＤ）か
らなり，前記生体磁場の前記生体の面に垂直な第１方向
の磁場成分の時間変化を検出する複数の磁束計と，前記
第１方向と交叉する第２方向及び第３方向に於ける前記
第１方向の磁場成分の変化率の２乗和の平方根に比例す
る値の時間変化を求める演算と，前記値の時間変化を所
定の期間で積分し積分値を求める演算とを行なう演算処
理手段と，前記積分値を表示する表示手段とを有し，前
記演算処理手段は，前記積分値を求める演算を，複数の
前記所定の期間で行ない前記積分値を複数個求め，該複
数個の前記積分値の間での，比，等加重を含む和または
差の何れかを求める演算を行ない，前記表示手段は，前
記積分値の等しい点を結ぶ等積分図を等高線の高低によ
って色分けをして３次元カラー表示することを特徴とす
る生体磁場計測装置。12. A time change of a magnetic field component of a biomagnetic field in a first direction perpendicular to a surface of the living body, comprising a quantum interference device (SQUID) arranged outside the living body and detecting a biomagnetic field emitted from the living body. And a plurality of magnetometers for detecting a change in a magnetic field component in a second direction and a third direction intersecting with the first direction in proportion to a square root of a sum of squares of a change rate of a magnetic field component in the first direction .
Calculating the time change of the value, and calculating the time change of the value.
And processing means for performing a calculation for obtaining the integrating integrated values at a constant period, and a display means for displaying the integrated value possess, before
The arithmetic processing means performs an operation for obtaining the integral value by a plurality of operations.
A plurality of the integral values are obtained in the predetermined period, and
A sum between several integrals, including ratios, equal weights, or
An operation for obtaining any of the differences is performed.
The contour map connecting points with the same integral value is expressed by the height of the contour line.
A three-dimensional color display by color-coding the biomagnetic field. 【請求項１３】生体の外部に配置され，前記生体から発
する生体磁場を検出する量子干渉素子（ＳＱＵＩＤ）か
らなり，前記生体磁場の前記生体の面に垂直な第１方向
の磁場成分の時間変化を検出する複数の磁束計と，前記
第１方向と交叉する第２方向及び第３方向に於ける前記
第１方向の磁場成分の変化率の２乗和の平方根に比例す
る値の時間変化を求める演算と，前記値の時間変化を所
定の期間で積分し積分値を求める演算とを行なう演算処
理手段と，前記積分値を表示する表示手段とを有し，前
記複数の磁束計が等間隔に配置され，前記演算処理手段
は，前記積分値を求める演算を，複数の前記所定の期間
で行ない前記積分値を複数個求め，該複数個の前記積分
値の間での，比，等加重を含む和または差の何れかを求
める演算を行ない，前記表示手段は，前記積分値の等し
い点を結ぶ等積分図を等高線の高低によって色分けをし
て３次元カラー表示することを特徴とする生体磁場計測
装置。13. A temporal change in a magnetic field component of a biomagnetic field in a first direction perpendicular to a surface of the living body, the quantum interference device being arranged outside the living body and detecting a biomagnetic field emitted from the living body. And a plurality of magnetometers for detecting a change in a magnetic field component in a second direction and a third direction intersecting with the first direction in proportion to a square root of a sum of squares of a change rate of a magnetic field component in the first direction .
Calculating the time change of the value, and calculating the time change of the value.
An arithmetic process for performing integration over a fixed period to obtain an integrated value
And management means, and display means for displaying the integrated value possess, before
The plurality of magnetometers are arranged at equal intervals, and
Calculates the integral value by a plurality of predetermined periods
To obtain a plurality of the integral values,
Finds the sum or difference between values, including ratios, equal weights
And the display means calculates the equality of the integral value.
The contour map connecting the points is colored according to the level of the contour line.
A biomagnetic field measuring apparatus for displaying a three-dimensional color image . 【請求項１４】前記生体表面に平行な面を直交座標の
ｘ，ｙ平面，前記生体表面に垂直な方向を直交座標のｚ
軸とし，前記生体の外部に配置され，前記生体から発す
る生体磁場を検出する量子干渉素子（ＳＱＵＩＤ）から
なり，前記生体磁場のｚ方向成分Ｂ_z（ｘ，ｙ）を検出
する複数の磁束計と，前記ｚ方向成分Ｂ_z（ｘ，ｙ）か
ら， Ｂ_xy＝｛（∂Ｂ_z（ｘ，ｙ）／∂ｘ）²＋（∂Ｂ_z（ｘ，
ｙ）／∂ｙ）²｝ の平方根に比例する値の時間変化を求める演算と，前記
値の時間変化を所定の期間で積分し積分値を求める演算
とを行なう演算処理手段と，前記積分値の等しい点を結
ぶ等積分図を表示する表示手段とを有することを特徴と
する生体磁場計測装置。14. A plane parallel to the surface of the living body is defined by an x, y plane of rectangular coordinates, and a direction perpendicular to the surface of the living body is defined by z of orthogonal coordinates.
A plurality of magnetometers, each of which is a quantum interference device (SQUID) disposed outside of the living body and configured to detect a biomagnetic field emitted from the living body, and detects a z-direction component B _z (x, y) of the biomagnetic field; From the z-direction component B _z (x, y), B _xy = ｛(∂B _z (x, y) / ∂x) ² + (∂B _z (x,
y) / {y) An operation for determining a time change of a value proportional to the square root of ² ｝ ;
Calculation to integrate the time change of the value over a predetermined period to obtain the integral value
And a display means for displaying an isointegral diagram connecting points having the same integral value . 【請求項１５】前記生体表面に平行な面を直交座標の
ｘ，ｙ平面，前記生体表面に垂直な方向を直交座標のｚ
軸とし，前記生体の外部に配置され，前記生体から発す
る生体磁場を検出する量子干渉素子（ＳＱＵＩＤ）から
なり，前記生体磁場のｚ方向成分Ｂ_z（ｘ，ｙ）を検出
する複数の磁束計と，前記ｚ方向成分Ｂ_z（ｘ，ｙ）か
ら， Ｂ_xy＝｛（∂Ｂ_z（ｘ，ｙ）／∂ｘ）²＋（∂Ｂ_z（ｘ，
ｙ）／∂ｙ）²｝ の平方根に比例する値の時間変化を求める演算と，前記
値の時間変化を所定の期間で積分し積分値を求める演算
とを行なう演算処理手段と，前記積分値の等しい点を結
ぶ等積分図を等高線の高低によって色分けをして３次元
カラー表示する表示手段とを有することを特徴とする生
体磁場計測装置。15. A plane parallel to the surface of the living body is defined as an x, y plane in rectangular coordinates, and a direction perpendicular to the surface of the biological body is defined as z in orthogonal coordinates.
A plurality of magnetometers, each of which is a quantum interference device (SQUID) disposed outside of the living body and configured to detect a biomagnetic field emitted from the living body, and detects a z-direction component B _z (x, y) of the biomagnetic field; From the z-direction component B _z (x, y), B _xy = ｛(∂B _z (x, y) / ∂x) ² + (∂B _z (x,
y) / {y) An operation for determining a time change of a value proportional to the square root of ² ｝ ;
Calculation to integrate the time change of the value over a predetermined period to obtain the integral value
And a point that equals the integral value.
3D by color-coding the contour isointegral diagram according to the level of the contour line
A biomagnetic field measurement apparatus, comprising: display means for performing color display. 【請求項１６】生体の外部に配置され，前記生体から発
する生体磁場を検出する量子干渉素子（ＳＱＵＩＤ）か
らなり，前記生体磁場の前記生体の面に垂直な法線成分
Ｂ _z （ｘ，ｙ）を検出する複数の磁束計と，前記生体の
面に平行な接線成分Ｂ _x を前記法線成分Ｂ _z （ｘ，ｙ）の
ｘ方向の変化率から，接線成分Ｂ _y を前記法線成分Ｂ
_z （ｘ，ｙ）のｙ方向の変化率から推定する演算を行な
う演算処理手段と，演算処理結果が表示される表示手段
とを有することを特徴とする生体磁場計測装置。16.It is located outside the living body and emits from said living body.
Quantum Interference Device (SQUID) to Detect the Changing Biomagnetic Field
A normal component of the biomagnetic field perpendicular to the surface of the living body
B _z A plurality of magnetometers for detecting (x, y);
Tangent component B parallel to the plane _x To the normal component B _z (X, y)
From the rate of change in the x direction, the tangent component B _y To the normal component B
_z An operation for estimating the change rate of (x, y) in the y direction is performed.
Calculation processing means and display means for displaying the calculation processing result
HavingA biomagnetic field measurement apparatus, characterized in that: 【請求項１７】生体の外部に配置され，前記生体から発
する生体磁場を検出する量子干渉素子（ＳＱＵＩＤ）か
らなり，前記生体磁場の前記生体の面に垂直な法線成分
を検出する複数の磁束計と，前記法線成分に基づいて，
２つの接線成分を推定する演算処理手段と，前記２つの
接線成分から求めた任意の時点での等磁場線図を表示す
る表示手段とを有することを特徴とする生体磁場計測装
置。17. A plurality of magnetic fluxes arranged outside of a living body, comprising a quantum interference device (SQUID) for detecting a biomagnetic field emitted from the living body, and detecting a normal component of the biomagnetic field perpendicular to the surface of the living body. And the normal component ,
Arithmetic processing means for estimating two tangent components ;
A display unit for displaying an isomagnetic field diagram at an arbitrary time point obtained from a tangent component .