JP3944383B2 - 心臓磁場計測装置 - Google Patents
心臓磁場計測装置 Download PDFInfo
- Publication number
- JP3944383B2 JP3944383B2 JP2001350920A JP2001350920A JP3944383B2 JP 3944383 B2 JP3944383 B2 JP 3944383B2 JP 2001350920 A JP2001350920 A JP 2001350920A JP 2001350920 A JP2001350920 A JP 2001350920A JP 3944383 B2 JP3944383 B2 JP 3944383B2
- Authority
- JP
- Japan
- Prior art keywords
- magnetic field
- subject
- cardiac magnetic
- data
- information
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Lifetime
Links
- 0 C(C1)*2C11C2CCC1 Chemical compound C(C1)*2C11C2CCC1 0.000 description 1
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
- A61B5/242—Detecting biomagnetic fields, e.g. magnetic fields produced by bioelectric currents
- A61B5/243—Detecting biomagnetic fields, e.g. magnetic fields produced by bioelectric currents specially adapted for magnetocardiographic [MCG] signals
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Heart & Thoracic Surgery (AREA)
- Surgery (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Engineering & Computer Science (AREA)
- Biomedical Technology (AREA)
- Cardiology (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Physics & Mathematics (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
- Measuring Magnetic Variables (AREA)
- Medical Treatment And Welfare Office Work (AREA)
Description
【発明の属する技術分野】
本発明は、生体の脳の神経活動や心臓の心筋活動等、生体内の電流などが原因で発生する生体の磁場を計測する生体磁場計測装置に係り、特に医師が脳疾患や心疾患を診断する診断サポート機能を備えた生体磁場計測装置に関する。
【0002】
【従来の技術】
磁気センサである超伝導量子干渉素子(SQUID)を用いて、生体から発生する微弱な磁場の分布を測定し、その測定結果から、生体内部の活動電流の位置を推定し、その分布をイメージングする多チャンネルの生体磁気イメージング装置の開発が進められている。たとえば特許第3140731号,3140732号等にそのような生体磁気イメージング装置に関する技術が開示されている。
【0003】
【発明が解決しようとする課題】
上記のような従来技術は、生体磁場図表示装置に関する動作原理および計測データを効率的に分りやすく表示するための方法に関するものであり、時間波形表示,等磁場線図表示,磁場−時間積分図表示などについて言及されている。しかし、これらの表示だけを見て心疾患等を診断するには心機能や心疾患に関する高度の知識、および心臓磁気計測の動作原理やデータ再構成方法に関する知識を必要とするため、一般の医師では、装置を使いこなすのが難しい状況である。更に、計測によって得られるデータが膨大であるために診断に用いるべき表示条件を選択,設定するのに長時間を要す、等の問題があった。
【0004】
本発明の目的は生体磁場計測装置で得られた計測データの特徴を定量的にとらえることによって心疾患と思われる被検者を推定し、さらに該被検者の計測データに基づき疾患の候補を推定することによって、医師の診断をサポートでき、疾患の見落しを防止し、また診断時間の大幅短縮が可能な診断支援機能を備える生体磁場計測装置を提供することにある。
【0005】
【課題を解決するための手段】
上記目的を達成するための本発明の構成は以下の通りである。
【0006】
(1)被検者の生体から発せられる磁場を複数の位置で計測する生体磁場計測装置において、磁気センサで計測された信号情報を記憶し、記憶された信号情報から複数の特徴パラメータを算出し、算出された特徴パラメータを、主に健常者の磁場計測結果に基づいて作成された基準データベースと比較することによって被検者が健常であるかどうかを識別する生体磁場計測装置。
【0007】
基準データベースとは健常者のデータに基づいて構築されたマハラノビス空間,ニューラル・ネットワークの判定論理を表すネットワーク等種々のものが採用できる。
【0008】
(2)被検者の生体内から発せられる磁場を複数の位置で計測する生体磁場計測装置において、磁気センサで計測された信号情報を記憶し、記憶された信号情報から複数の特徴パラメータを算出し、算出された特徴パラメータを、予め定められたルールと比較対照することによって被検者が健常であるかどうかを識別する生体磁場計測装置。
【0009】
予め定められたルールとはif then else形式でルールを定めた集合体で形成されたエキスパートシステム等がある。
【0010】
(3)前記(1)または(2)において、計測する磁場は主に心臓から発せられる磁場であり、かつ抽出する特徴パラメータは、R波ピーク近傍の一定時間ごと電流方向である生体磁場計測装置。
【0011】
(4)前記(1)または(2)において、測定する磁場は主に心臓から発せられる磁場であり、かつ抽出する特徴パラメータは、計測した磁場強度のうちの最も振幅の大きな時刻を代表時刻とし、複数の磁気センサで測定されたその代表時刻での磁場を、磁場の等しい点で結んだ、いわゆる等磁場線図に基づいて抽出する生体磁場計測装置。
【0012】
(5)被検者の生体内から発せられる磁場を複数の位置で計測する生体磁場計測装置において、計測した信号情報を記憶し、記憶された信号情報から複数の特徴パラメータを算出し、算出された特徴パラメータに基づいて、本明細書中で後記するマハラノビス空間を構築し、そのマハラノビス空間における前記被検者の、複数の特徴パラメータに基づいたマハラノビスの距離を計算し、そのマハラノビスの距離の情報に基づいて被検者が健常であるかどうかを識別するための識別情報を算出する識別情報算出手段を備えた生体磁場計測装置。
【0013】
(6)(1),(2),(5)のいずれかにおいて、前記識別情報算出手段により算出された識別情報と、被検者が健常であるかどうかを識別するため予め記憶された識別基準情報とを比較して、被検者が健常であるかどうかを識別する識別手段を備えた生体磁場計測装置。
【0014】
(7)(6)において、識別基準情報を操作者が任意に設定可能である生体磁場計測装置。
【0015】
(8)(6)において、識別手段により被検者が健常でないと識別されたときには、測定された複数の特徴パラメータに基づいて、疾患の候補を表示する疾患推定機能を備えた生体磁場計測装置。
【0016】
(9)(5)〜(8)のいずれかにおいて、複数の被検者毎の特徴パラメータ,マハラノビスの距離を記憶する記憶手段を備え、選択した被検者の少なくとも1つの特徴パラメータ,マハラノビスの距離を表示する機能を備えた生体磁場計測装置。
【0017】
(10)(5)〜(8)のいずれかにおいて、複数の被検者毎の特徴パラメータ,マハラノビスの距離を記憶する記憶手段を備え、選択した被検者の少なくとも1つの特徴パラメータ,マハラノビスの距離と、構築されたマハラノビス空間におけるマハラノビスの距離及び特徴パラメータを同時に表示する機能を備えた生体磁場計測装置。
【0018】
(11)(9)または(10)において、選択した被検者のマハラノビス空間における平均値と標準偏差とを表示する機能を備えた生体磁場計測装置。
【0019】
(12)(9)または(10)において、選択した被検者の特徴パラメータと基準化した特徴パラメータを併記する機能を備えた生体磁場計測装置。
【0020】
(13)(9)または(10)において、マハラノビスの距離が、被検者が健常であるかどうかを識別するために、予め記憶された識別基準情報以上となるデータを有する被検者のデータを強調表示する機能を備えた生体磁場計測装置。
【0021】
(14)(9)または(10)において、マハラノビスの距離が、被検者が健常であるかどうかを識別するため予め記憶された識別基準情報以上となるデータの時に、基準化した特徴パラメータが予め定められたしきい値以上の特徴パラメータを強調表示する機能を備えた生体磁場計測装置。
【0022】
(15)(9)または(10)において、記憶された複数の被検者のデータを、マハラノビスの距離の大きさによって表示順を昇順あるいは降順に並び替える機能を備えた生体磁場計測装置。
【0023】
(16)(5)において、複数のマハラノビス空間データを記憶する記憶手段を備え、被検者情報に基づいて使用する該空間データを選択できる機能を備えた生体磁場計測装置。
【0024】
(17)(16)において、空間データごとに特徴パラメータを定義できる機能を備えた生体磁場計測装置。
【0025】
(18)(1)〜(5)のいずれかにおいて、被検者の生体内から発せられる磁場が主に心臓から発せられる心臓磁場であり、かつ心臓磁場信号情報のうちの、P波・QRS波・T波近傍の等磁線図を1画面上に表示する機能を備えた生体磁場計測装置。
【0026】
【発明の実施の形態】
以下に図を用いて本発明の実施例を説明する。
【0027】
図1は本発明が実施される生体磁気計測装置の一実施例の概略構成を示す。環境磁気雑音の影響を除去するために、生体磁気計測装置は磁気シールドルーム1内に設置される。生体からなる被検体である被検者2は、ベッド3に仰向け状態で計測が行われる。被検者の生体面(胸部の場合は一般に胸壁に平行面)はベッド3の面とほぼ平行であるとし、そしてこの面は直交座標系(x,y,z)のx−y平面と平行であるものとする。被検者の胸部は曲面であると共に傾いているが、説明を簡単にするためにほぼ平行とする。
【0028】
被検者2の胸部の上方には、冷媒である液体Heで満たされたデユワ4が配置され、該デユワは超伝導量子干渉素子(SQUID=Superconducting Quantum Interference Device )とそのSQUIDに接続された検出コイルとを含む複数個の磁気センサを収容している。液体Heは磁気シールド1の外部にある自動補給装置5から連続的に補給される。
【0029】
磁気センサの出力は、被検者2から発生して検出コイルにより検出される生体磁場の強度(磁束密度と考えることもできる)と特定の関係をもつ電圧を出力し、その出力がFLL(Flux Locked Loop)回路6に入力される。このFLL回路6は、SQUIDの出力を一定に保つように、SQUIDに入力された生体磁場(生体磁気)の変化を帰還コイルを介してキャンセルする(これを磁場ロックと呼ぶ)。その帰還コイルに流した電流を電圧に変換することにより、生体磁場信号の変化と特定の関係にある電圧出力を得ることができる。このように帰還コイルを介して検出する方式を取っているので、微弱の磁場を高感度に検出できる。
【0030】
上記出力電圧は増幅器・フィルタ・増幅器(AFA)7に入力され、その出力はサンプリングされて、A/D変換され、計算機8に取り込まれる。
【0031】
計算機8はパーソナルコンピュータからなり、8−1はそのディスプレイ部、8−2はキーボード、そして8−3はマウスを示す。マウス8−3は画面上でカーソルを移動させて処理対象を選択するのに用いられる。この操作はキーボードを操作することによっても行うことができる。AFA7の入力ゲイン(Igain)及び出力ゲイン(Ogain)は調整可能であり、また、AFA7は第1の基準周波数以下の周波数信号を通過させるローパスフィルタ(LPF),第1の基準周波数よりも低い第2の基準周波数以上の周波数信号を通過させるハイパスフィルタ(HPF)及び商用電源周波数をカットするノッチフィルタ(BEF)を含む。計算機8は各種の処理を行うことができ、その処理結果はディスプレイ部8−1に表示され得る。
【0032】
SQUIDとしては例えば1例として直流SQUIDが用いられる。SQUIDに外部磁場が与えられたときに、それに対応する電圧(V)が発生するようにSQUIDには直流バイアス電流(Ibias)が流される。その外部磁場を磁束Φで表すと、VのΦに対する特性曲線すなわちΦ−V特性曲線は周期関数で与えられる。計測に当っては、それに先立って、FLL回路6のオフセット電圧(VOFF)を調整してΦ−V特性曲線の直流電圧をゼロレベルにする操作が行われる。更に、AFA7の入力がゼロのときその出力がゼロとなるようにAFA7のオフセット電圧(AOFF)の調整が行われる。
【0033】
SQUIDに外部から大きな磁場がかかると、その磁場がSQUIDによってトラップされ、その正規の動作がなされなくなる。すなわち、計測信号の時間波形を見るに、信号のベースラインがステップ上に変動したり、ノイズが極端に大きくなる。またそのSQUIDセンサのΦ−V特性を見れば、正常時の特性曲線はΦの走査に対して位相が一定の周期曲線となるが、磁束トラップが起こると位相が変わり、また特性曲線の振幅が正常時と比較して小さくなるため、磁束トラップの発生を検知することができる。その場合は、SQUIDを加熱して一旦常伝導状態にし、そしてその後その加熱を止めて、トラップされた磁場を取り除くことができる。その場合の、SQUIDの加熱操作をヒートフラッシュと呼ぶ。
【0034】
図2は磁気センサの配置構成を示す。磁気センサの検出コイルには生体磁場の接線成分(生体面すなわちx−y平面にほぼ平行成分)を検出するコイルと生体磁場の法線成分(生体面すなわちx−y平面に直交する成分)を検出するコイルがある。生体磁場の接線成分を検出するコイルとしては、コイル面がx方向及びy方向をそれぞれ向いた2つのコイルが用いられ、また、生体磁場の法線成分を検出するコイルとしてはコイル面がz方向を向いたコイルが用いられる。複数個の磁気センサ20−1〜20−8,21−1〜21−8,22−1〜22−8,23−1〜23−8,24−1〜24−8,25−1〜25−8,26−1〜26−8及び27−1〜27−9は、図2に示されるように、生体面すなわちx−y平面とほぼ平行面上にマトリックス状に配置される。磁気センサの数は任意であってよいが、図2では、磁気センサのマトリックスは8行8列からなっているから、磁気センサの数は8×8=64である。各磁気センサは、図2に示されるように、その長手方向が生体面すなわちx−y平面に対して垂直な方向(z方向)と一致するように配置される。なお、この一実施例ではベッド面とセンサのX−Y面とを平行にしているが、測定精度を高めるには体に接近させる方が良く、傾けるようにすることができる。但し、被検者である人体は常に動いているので、人体に密着させるとこの動きが検出部を動かし、かえって高精度の検出が困難となる。
【0035】
図3は磁気センサの各々の、生体磁場の法線成分Bzを検出するセンサの構成を示す。同図において、超伝導線(Ni−Ti線)で作られたコイルはそのコイル面がz方向を向くように配置される。このコイルは互いに逆向きの2つのコイル10及び11の組み合わせからなり、被検者2に近い方のコイル10は検出コイルとされ、遠い方のコイル11は外部磁場雑音を検出する参照コイルとされる。外部磁場雑音は被検者よりも遠い信号源から生じており、したがって、その雑音信号は検出コイル10及び参照コイル11の両方によって検出される。一方、被検者からの磁場信号は微弱であり、したがって、その生体磁場信号は検出コイル10によって検出されるが、参照コイル11はその生体磁場信号にほとんど感応しない。このため、検出コイル10は生体磁場信号と外部磁場雑音信号を検出し、参照コイル11は外部磁場雑音信号を検出するから、両コイルで検出された信号の差をとることによりS/N比の高い生体磁場の計測が可能となる。これらのコイルはSQUID12を実装した実装基板の超伝導線を介してSQUIDの入力コイルに接続され、これによって、検出された生体磁場信号の法線方向の成分BzがSQUIDに伝達される。
【0036】
生体磁気の接線成分Bx,Byについては、検出コイル面と参照コイル面をx方向あるいはy方向に向けることによって検出することができるが(図示せず)、図3の磁気センサで得られた法線成分Bzをx,yについて偏微分することによって実際の接線成分Bx,Byと強い相関を持った信号を得ることができるため、法線成分Bzをxおよびyについて偏微分した値を接線成分としてもよい。この場合は一つの磁気センサで接線成分Bx,Byと法線成分Bzとの両方を検出し、測定することができる。
【0037】
図4は磁気センサと被検者2の被計測部である胸部30との位置関係を示す。示されている点は図2に示されるマトリックス上の行と列との交点すなわち被検者2の計測点すなわち計測位置を表す。これらの各計測位置をチャンネルとも呼ぶ。図からわかるように、この実施例では、被検者2の身長方向をy方向とし、被検者2の横方向をx方向としている。
【0038】
図5はある健常者について法線成分Bzの心磁波形を示す。
【0039】
心磁波形は心電波形のように、心房の興奮によるP波,心室の興奮過程によるQRS波、および心室興奮の回復過程によるT波からなる。P波とQ波の間、S波とT波の間には心臓の興奮による信号は略0(ゼロ)と考えられ、この間の信号レベルをベースラインと呼ぶ。心臓の心室が脱分極した時間帯T1つまり収縮期のQRS波でのそれぞれの波形ピークの時間がtQ ,tR 及びtS としてそれぞれ示されている。また、心臓の再分極過程(拡張期)であるT波の時間帯はT2として示されている。
【0040】
このような時間波形が64個のセンサで同時に計測されているため、ある時点における各センサの計測信号(すなわち磁場強度)を2次元で表した等磁場線図を作成することによってその時点における磁場分布を得ることができる。図中、各センサ位置に対応した位置に矢印を記しているが、この方向は計測された接線方向の磁場データあるいは法線方向の磁場データをXおよびY方向に偏微分して得られた擬似的な接線方向の磁場データから求めた磁場の方向を表しているが、簡単のために反時計方向に90度回転させることによって、その磁場を発生させる電流の方向に一致させている。また矢印の長さは磁場強度を表している。これによって接線方向の等磁場線図に矢印を合成したアローマップは、色の濃さや矢印の大きさによって磁場の分布を表し、矢印の方向によって電流方向を表している。
【0041】
図6は、典型的な健常者(37歳,男性)のP波,QRS波のピーク(R),T波の各事象での法線方向および接線方向の磁場によるアローマップを示している。(a)〜(b)はP波,QRS波,ST−T波の各事象での磁場分布の時間変化を調べた結果をそれぞれ示す。これら等磁場線図は同じ磁場強度のところを線で結んだ等高線を示している。また、図中の矢印はセンサ各点での磁場から等価的な電流に置き換えたものを示している。
【0042】
法線成分からは磁場の湧き出し、吸い込みの2方向の磁場があることが分かり、両極の谷間に電流ダイポールを想定することができる。その位置と方向から、電流ダイポールは右上部から真下に動いていることが分かる。これは右心房の洞結節から興奮が始まった現象を捉えていると考えられる。一方、法線成分から算出した接線成分を見ると、先ほど推定した興奮部位の場所や方向が一目で分かりやすい形で表現されている。
【0043】
QRS波では、初期には右下方の向きに電流ダイポールが推定される。すぐにこの向きは左下方に変わり大きな磁場が検出されている。この方向の磁場が最大値を示した後徐々に減衰している。後期では向きが右上方に変わり大きくなった後、減衰してQRSが終わっている。これらを接線成分で見ると興奮部位が画面中央で右下方向に電流が流れている。R波では急激に向きが変わり左下方向に電流が流れ始めている。時間とともに電流が大きくなりR波のピークに達した後、徐々に強度が減衰する。更に、S波になると右上の方向に電流が流れ消滅した。これらの興奮を心臓の部位で考えると、QRSの初期の電流は心室中隔の興奮が始まり右室に到達するまでの興奮に対応していると推定される。さらにR波では左室に興奮が移り心筋の大きさからも興奮が非常に大きく現れたと考えられる。左室の興奮が終わると流出路の方に興奮が移り心室での脱分極過程が終了した。
【0044】
ST−T波では、心室の再分極を表しており心尖部の方向に全体が興奮し続け、強度が増大したのち、徐々に減衰していった。その間ほとんど方向は変わらず、強度だけが変化した。この変化は法線成分および接線成分から同様に読み取ることができる。
【0045】
図7は57歳男性の陳旧性心筋梗塞患者のQRS波における等磁場線図である。初期ベクトルが左横に向いており強度が強いが、これは心電図の異常Q波に対応するもとの考えられる。また、T波においても方向性が健常者のものと異なり、特に興奮部位は2ヵ所に分離してマルチダイポール性が顕著である。これらから等磁場線図上での電流方向,時間変化及びマルチダイポール性で虚血性心疾患の診断が可能である。
【0046】
図8は22歳男性のWPW症候群患者の等磁場線図である。WPW症候群は正常の房室伝導路の他に異常な副伝導路があり、心房からの興奮を早く心室に伝える異常である。健常者の心磁図と比較すると、R波ピーク前の副伝導路による異常興奮があり、R波ピークにおける電流方向にも異常が見られる。等磁場線図を作成することにより、このようにして時間波形では分らない、心臓内における電気生理学的な情報を得ることができる。
【0047】
本発明では、等磁場線図や時間波形から疾患の特徴を表す特徴パラメータをパソコンによって自動的に抽出し、MTS技法(Mahalanobis-Taguchi System)によって被検者が健康か、あるいは何らかの異常があるかを判別するものである。MTS技法については日本規格協会:「技術開発のための品質工学」等の文献に詳細に述べられているが、ここで簡単に説明しておく。
【0048】
最初に既に健常者であることが分っている被検者に対して心磁計測を行い、それぞれのデータから特徴パラメータ(Y1 ,Y2 ,…,Yn )を抽出する。そして全ての特徴パラメータに対して平均値mと標準偏差σを求めて,項目ごとに次の基準化という変換を行う。
【0049】
【数1】
【0050】
基準化した特徴パラメータYi について次の式によって定義されるマハラノビス(Mahalanobis)の距離D2 を求める。
【0051】
【数2】
【0052】
ここでaijは相関行列の逆行列のaij成分の値、Y1 ,Y2 ,…,Yk は基準化した個人のk項目の値である。この場合、マハラノビス(Mahalanobis)の距離D2 は平均1.0となる。
【0053】
MTS技法は、例えば集団健康診断において得られた尿検査,血液検査等の結果に、更に性別,年齢,たばこ本数,喫煙年数,アルコール合数,飲酒年数,肥満度等の個人情報を加えた多次元情報に基づき基準空間を構築し、個々の被検者毎にマハラノビスの距離を算出する。この際に健常者の集団を特徴づけるものさしを作ると、マハラノビスの距離により、個々の被検者が健常であるか、何らかの異常をもつかが判別できるため、異常をもつ人を正常であると判断して早期治療の時期を逸したり、健常な人を異常であると判断して過剰な精密検査を行う誤りをなくすることができるのではないかと思われている。詳細は例えば品質工学Vol.7 No.2(1999年4月)MTSを用いた健康診断による総合判定の信頼性向上 長谷川良子他に記載されている。
【0054】
本発明は、このMTSを生体磁場計測方法に応用するべく検討している上でなされたものである。生体磁場計測は、複数の磁気センサで得られた磁場強度を示す信号が得られるだけである。得られた信号情報をどのように処理して診断に適用していくかは現在検討されている段階であり、まだ確固たる処理方法は提案されていない。本発明は、情報処理の一つの試みとして、複数の磁気センサの信号情報から複数の特徴パラメータを算出し、その特徴パラメータに基づき基準空間を構築してマハラノビスの距離を算出しようとするものである。マハラノビスの距離に基づいて、まず健常者か、何らかの異常を有する可能性をもった被検者かを判別し、異常を持っている可能性のある被検査者については、該被検者の特徴パラメータに基づいて疾患の推定を行う。疾患の推定については、いくつかの疾患候補を挙げる。それら疾患の可能性をそれぞれ定量的な数字であらわすこともできる。医師はその疾患候補に基づき、生体磁場計測以外の方法を用いて精密検査を行い、最終的な診断を下すことが可能になる。これにより、生体磁場計測に不慣れな医師でも、生体磁場計測装置を有効に使いこなすことができるようになる。また、被検者の生体磁場計測結果は、同一病院内での計測結果のみならず、広く市単位,県単位,国単位でのデータを共通のデータベースで管理することにより、より正確な診断が可能になる。すなわち、健常者かどうかを判別するための基準となるマハラノビスの距離を設定するに当たっては、より多くのデータを参照した上で決めることが望ましいが、上記データベースを用いることによりそれが可能となる。また生体磁場計測にあっては、性別,年齢、または被検者の居住地域等により、健常者かどうかを判別するための基準となるマハラノビスの距離が代わってくる可能性があるが、上記データベースのデータを性別,年齢等により細分したデータを健常者かどうかを判別するための基準となるマハラノビスの距離算出の基準とすることにより、より正確な診断が可能になる可能性がある。
【0055】
図9に本発明での特徴パラメータの1例を示す。特徴パラメータは31項目からなり、性別(sex),身長(height),体重(weight) ,P波のピークおよびR波のピークを与える時間の間隔(Tpr),R波ピークとT波ピークの間隔(Trt),P波ピークにおける磁場強度(P−pt(B)),電流方向(P−pt(Ap)),センサのX位置(P−pt(ixp))とY位置(P−pt(iyp))、同じくR波ピークにおける電流方向(R−pt(Ar)),センサのX位置(R−pt(ixr))とY位置(R−pt(iyr))、さらにT波ピークにおける磁場強度(T−pt(B)),電流方向(T−pt(At)),センサのX位置(P−pt(ixt))とY位置(P−pt(iyt)),R波ピーク時刻から10ミリ秒から90ミリ秒までの間で10ミリ秒ずつ間隔を置いて計測した電流方向AP2,…,Ap10,ST−T波における最大磁場強度(ST−Tmax)、QRS波における最大磁場強度(QRSmax)ST−T波とQRS波における最大磁場強度の差の最大値([ST−T−QRS]max)および最小値([ST−T−QRS]min)、QRS波時間積分値(QRSsum)、およびST−T波とQRS波との時間積分値の差([ST−T−QRS]sum)の項目を含んでいる。なお、上記の34個のパラメータにてマハラノビスの距離計算を行ったところ、R波ピーク10ms後の電流方向(Ap1)はR波における電流方向(R−pt(Ap))との相関、およびST−T波時間積分値(ST−Tsum)はST−T波における最大磁場強度(ST−Tmax)との相関が強すぎて相関行列の逆行列Aが計算できないため、使用しなかった。またMTS法によるパラメータの評価を行ったところ年齢を特徴パラメータとして使用した場合の望大特性SN比の改善が大きすぎ、基準空間を作成した年齢が20〜30代に偏っていたためと分り、本評価には使用しなかった。このため、31項目によるマハラノビス距離を診断の尺度とした。
【0056】
上記の特徴パラメータによって健常者48名から基準空間を生成し、虚血性心臓病患者42に対してマハラノビスの距離D2 を計算したところ、図10に示すように健常者と心臓病患者がしきい値D2 =2.0 で判別することができた。全ての健常者はマハラノビスの距離が2.0 以下であるのに対して、心臓病患者はすべてD2 >2.0 となった。この中にはマハラノビスの距離が10.0 以上で、グラフの表示範囲外となったデータが81%あり、このため図中に「級外81%」を記してある。健康なヒトが次第に心疾患に罹る場合には、マハラノビスの距離は健常者の範囲から心疾患が重くなるにつれて大きな値を示すようになる。このため、心疾患の初期の状態では健常者と区別するのは難しい場合もありうるが、適切な特徴パラメータとホモジニアス(均一な)健常者集団を用いればマハラノビスの距離D2は心疾患の診断に有効であると言える。
【0057】
図9に示した特徴パラメータの要因効果図を図25に示した。要因効果図の縦軸は望大特性のSN比[dB]をあらわしている。横軸は1から31番までのパラメータ番号を表しており、その上段の1,2という数値は第1水準と第2水準を表している。第1水準はその特徴パラメータを使用した場合、第2水準は使用しなかった場合である。ほとんどのパラメータは右下がりとなっているがこれはその特徴パラメータを使用することによって、健常者と患者をより精度良く識別できることを示しており、右下がりの度合いが大きいほど有効なパラメータと言うことができる。
【0058】
図11はマハラノビスの距離によって被検者が健常者であるか、心臓病患者であるかを判別するためのソフトウェアのブロック図を示しており、計算機8に実装されているプログラムの機能単位の関連を示したものである。
【0059】
被検者/データリスト表示部は計測やデータ処理の対象となる被検者やデータの一覧を表示し、オペレータの選択を受け付けるものである。選択した被検者情報あるいはデータ情報はデータ計測部やデータ解析・表示部に送られる。計測する被検者情報が心磁データベースに登録されていない場合には情報登録部によって被検者情報を登録する。
【0060】
データ計測部はFLL制御回路やAFA回路を制御して、被検者からの心磁信号の計測を行う。被検者情報は被検者/データリスト表示部から、計測条件は測定条件設定部から送られる。データ計測部で計測されたデータは心磁データベースに格納されると同時に、特徴パラメータ抽出部によって31項目の特徴パラメータが計算され、マハラノビス距離計算部によってマハラノビスの距離が計算される。マハラノビス距離計算部では基準空間データを参照しているが、これは基準空間登録部によって登録や更新が行われる。一般的に、基準空間は年齢や性別・職業・人種など、被検者の属性に関する情報も特徴パラメータとして計算することもできるが、基準空間を定義するための健常者データに偏りがあり、各パラメータに対するデータの均一性を確保しにくい場合には、これらの被検者の属性データによって健常者を層別し、それらに対応した基準空間を別々に登録してもよい。データ解析・表示部は被検者/データリスト表示部で選択された被検者情報およびデータ情報から心磁データベースに保存されている計測データを呼び出し、診断に用いる様々な解析図をGUIを通してインタラクティブに生成する。
【0061】
被検者の登録から、その登録された被検者のデータ計測を行って、その計測されたデータの解析を行うまでの一連の操作はディスプレイ8−1に表示される表示画面を見ながら行われる。このため、その一連の操作の説明に先立ってまずその表示画面のレイアウトを説明する。
【0062】
図12は、図1のディスプレイ8−1に表示される表示画面の基本的なレイアウトを示す。表示画面の上方部は上から順番に配置されたタイトルバー部801,メニューバー部802及びアイコンが配置されているツールバー部803によって占められる。上記各部は表示領域やエリアと考えることもできる。これらの配置は他の処理目的、たとえば、被検者の登録や読み出し,磁場の計測,計測データの解析のための処理、などにおいてもその表示画面で共通して表示される。これにより使用し易さが増し、計測や処理の時間が短縮できる。
【0063】
表示画面の中央部は左から右へ順番に配置された被検者情報部804−1および解析データについてのデータ情報部804−2,線図や波形のような解析データが表示される解析データ部805−1および参照波形部805−2及び操作領域部806によって占められている。
【0064】
被検者情報部804−1には、被検者リスト画面(図18)が表示されているときはその画面中の被検者リストにおいてカーソル91が置かれる被検者の情報が常時表示され、また、解析データ部に線図や波形のような解析データが表示されているときは(図21〜図24)、その表示されている解析データが得られた被検者の情報が常時表示される。これによって、表示されている解析データとその解析データが得られた被検者との関係を明確に知ることができる。同様にデータ情報部804−2には、被検者リスト画面(図18)が表示されているときはその画面中で選択された被検者についてのデータリストが表示されており、そのデータリスト中においてカーソル92が置かれるデータの情報が表示され、また、解析データ部に線図や波形のような解析データが表示されているときは(図21〜図24)、その表示されている解析データの情報が常時表示されるため、表示されている解析データに関する計測時刻や計測条件などの情報を明確に知ることができる。このように、このシステムの表示画面においては、メニューバー部802と同様に、常に、被検者情報部804−1およびデータ情報部804−2が表示画面の定位置(左側)に表示されているので、使用者は表示画面が変わるたびに被検者情報エリアを探す必要がなく、常に表示画面の所定位置(左側)を見れば知ることができる。
【0065】
図14にデータ情報表示部の例を示す。マハラノビスの距離D2 は計測データの属性であるから、データ情報部804−2に表示される。レーダーチャートは特徴パラメータの各項目の値を示している。これで異常の有無や、おおよその原因(特異な値を持つ特徴パラメータ)を知ることができるが、下の詳細ボタンを押すことによって図20の詳細なレーダーチャートが表示される。
【0066】
タイトルバー部にはフレームの名称、具体的には、「MultichannelMCGSystem」という名称が表示される。
【0067】
図13に操作メニューを示す。メニューバー部は操作メニューを選択する部分で、「ファイル(F)」,「被検者リスト(L)」,「データ計測(Q)」,「解析(A)」の各メニューを使用する事ができる。
【0068】
図13の操作メニューの内容は、これらのメニューの内容はそれぞれ対応するメニューボタンをクリックすることによってプルダウンメニューとして表示される。このため、操作メニューを必要としないときは、前記各メニューを呼び出すためのキーワードのみメニューバー部にコンパクトに表示しているので、解析データ部や操作領域部等の各作業で必要な表示エリアを広く設定することができる。そして、操作メニューが必要なときは、操作手順にしたがって配列された前記キーワードをメニューバー部から選択することによって表示して操作指示することができる。この際、前記キーワードは、文字の配列(左から右)に準じて配列されているので、自然な形で操作指示することができる。
【0069】
「ファイル(F)」のプルダウンメニューは、MultichannelMCGSystem を終了させる「心磁システムの終了(X)」という項目を含む。
【0070】
「被検者リスト(L)」のプルダウンメニューには「被検者リストを開く(O)」,「被検者登録(R)」,「被検者削除(D)」,「データ削除(E)」という項目を含む。「被検者リストを開く(O)」が選択されるとディスプレイ装置8−1に表示される画面が被検者リスト画面(図18)に切り替えられる。「被検者登録(R)」が選択されると被検者登録ダイアログ(図17)が表示されて、被検者ID,氏名,生年月日,身長,体重,性別、およびコメントなどの入力を受付ける。「被検者削除(D)」が選択されると被検者リスト(図18)の被検者情報カーソル91が置かれた被検者情報およびその被検者に関するすべてのデータが削除され、「データ削除(E)」が選択されると被検者リスト画面(図18)のデータリスト上のデータ情報カーソル92が置かれているデータが削除される。
【0071】
特権ユーザモードにおける「データ計測(Q)」のプルダウンメニューは「計測モニタ画面を開く(O)」と「計測開始(M)」という項目を含む。「計測モニタ画面を開く(O)」がクリックされると、ディスプレイ装置8−1上の画面表示を図21に示す計測モニタ画面に切り替える。「計測開始(M)」がクリックされると、センサ状態の自動調整を実行し、すべてのSQUIDセンサの磁場ロックを行い、指定された条件でデータの取り込みを行う。計測が終わると計測データは心磁データベースに保存されるが、同時に計測データから特徴パラメータが抽出されてマハラノビスの距離が計算され、計測したデータの情報として心磁データベースに保存される。
【0072】
「データ解析(A)」のプルダウンメニューは「時間波形表示(W)」,「等磁線図(B)」「時間積分図(T)」「サマリ表示(S)」「マハラノビス距離再計算(M)」「ベースライン補正(C)」「基準空間設定(N)」を含む。
【0073】
「時間波形表示(W)」「等磁線図(B)」「時間積分図(T)」「サマリ表示(S)」がクリックされるとそれぞれ時間波形表示画面(図21),等磁線図画面(図22),時間積分図(図23),サマリ表示(図24)が表示される。「マハラノビス距離再計算(M)」がクリックされると表示中のデータについてマハラノビスの距離を計算しなおして、データベースに登録されている値も更新する。「ベースライン補正(C)」がクリックされるとベースライン、すなわち心臓の電気生理学的な振舞いが休止している状態における心磁信号のレベルにオフセットを加えて0にする処理である(図示せず)。また、「基準空間設定(N)」を指定すると基準空間を定義してあるファイルと、それを適用する被検者の範囲を指定するためのダイアログが表示される。
【0074】
ツールバー部803には操作メニューのプルダウンメニューの項目の中で使用頻度の高いものと関係づけられたアイコンボタンを配置してある。
【0075】
次に、被検者の登録から、その登録された被検者のデータ計測を行って、その計測されたデータの解析を行うまでの一連の操作を図15〜図24を参照しながら説明する。
【0076】
図15は本実施例の生体磁気計測装置の全体に関する操作フローを示すもので、計算機8の電源がONにされると(S−1)、オペレーティングシステムが立ち上げられ、計算機で利用可能なプログラムの起動アイコンがディスプレイ部8−1に表示される(S−2)。そのアイコンの中からMultichannelMCGSystem のプログラムのアイコンがクリックされると、システム立ち上がりの初期画面として図18に示す被検者リストが表示される(S−3)。
【0077】
図18に示される被検者リスト画面について説明するに、その左上部は被検者情報部804で、左下部はデータ情報部805で占められる。また、その右側全体の上部には被検者リストが、下部にはデータリストが表示されるようになっている。被検者情報部には被検者リスト上でカーソル91が置かれた被検者の情報が表示され、カーソル91を移動すると表示の内容もそれに合わせて更新される。被検者リストの項目は、ID(被検者ID番号),氏名,登録年月日(データ登録された日),マハラノビスの距離D2 ,生年月日,年齢,身長,体重,コメント(被検者に関するコメント)等を含む。被検者リストについては、これを縦スクロールバーでスクロールすることができ、被検者リストの項目については、これを水平(横)スクロールバーでスクロールすることができる。選択された被検者の行は強調表示される。
【0078】
図18の被検者リスト画面の下半分には計測データリストが表示されており、データのIDやマハラノビスの距離D2 、および計測日付などのデータ属性を表示する。同様に等磁線図リストには計測データから等磁場線図を作成するための情報を表示する。
【0079】
被検者リストと計測データリストの両方にマハラノビスの距離D2 が表示されるが、前者は全計測データのマハラノビスの距離D2 の中で最大のものを表示している。図中、カーソルが置かれている上から2番目の患者(ID=0000000001,氏名=BBBBBBBB)はマハラノビスの距離D2 =3.3 を示しているが、これはその下の3つの計測データのマハラノビスの距離3.0 ,3.3,2.9のうち、最大の3.3 を表示している。
【0080】
図15の全体のフローに戻る。ステップS−4においては、被検者リスト画面上の被検者リストの中から所望の被検者の行が選択される。この後のフローはメニューにより4つに分岐される(S−5)。分岐の一つによれば、「ファイル(F)」というメニューの「心磁システム終了(X)」というサブメニューが選択され、この場合はウインドウを閉じる等の終了処理が行われ(S−8)、それによってシステムがシャットダウンされる(S−9)。その後、計算機8の電源がOFFにされ(S−10)、すべてが終了する。
【0081】
分岐の残りによれば、データ計測(S−6),データ解析(S−7)あるいは基準空間の設定(S−71)が行われる。データ計測は「データ計測(Q)」というメニューの「計測モニタ画面を開く(O)」というサブメニューを選択することにより実行可能である。更に、データ解析は「データ解析(A)」というメニューの「時間波形表示(W)」,「等磁線図(B)」「時間積分図(T)」「サマリ表示(S)」というサブメニューのいずれかを選択することにより実行可能である。ステップS−6,S−7及びS−71の終了のフローは、ステップS−3に戻って被検者リスト選択画面を表示する。ステップS−4の被検者リスト選択、ステップS−6のデータ計測及びステップS−7のデータ解析のフローについては、その詳細はそれぞれ図16〜図17に関連して以下に更に詳しく説明される。
【0082】
図16は図15のステップS−6におけるデータ計測の詳細なフローを示す。まず計測画面の初期状態として図21に示すグリッドマップによる時間波形を表示し(S−15−1)、波形モニタを起動する(ステップS−15−2)。チャンネルは8×8の64チャンネルからなり、「全チャンネル選択」ボタンをクリックするか、又はチャンネルマトリックスを対角線に沿って端から端までドラッグすることによって、全チャンネルを選択して表示することができる。
【0083】
波形モニタは、予め設定されている周期時間(例えば1秒間)で、データを取り込んで表示し、それを計測ボタンが押されるまで繰り返す。計測条件パラメータを設定・変更(ステップS−15−3)では、サンプリング時間やサンプリング間隔などの計測条件の設定および変更を行い、操作領域にある「計測」ボタンを押して計測を開始する(S−15−4)。サンプリングの時間(計測時間)及び間隔については、逆三角印のついた対応テキストボックスをクリックすれば、選択可能な数値のプルダウンメニューが開かれ、その中から所望の数字を選択することができる。その選択可能な数字は、時間については、たとえば1sec,5sec,10sec,30sec,1min及び2minであり、間隔についてはたとえば0.1msec,0.5msec,1.0msec,2.0msec,4.0msec,5.0msec 及び10.0msecである。時間は必要に応じて1sec程度から24h程度までの間で選ばれるようにしてもよい。「スケール」ボックス内の「時間」とあるのはmsec 単位の時間スケールすなわち水平方向のスケールを意味する。これらについても、サンプリングの時間及び間隔の選択と同様に、対応するテキストボックスをクリックすることによって開かれるプルダウンメニューの中から所望の数値が選択される。
【0084】
FLL制御回路6およびアンプ・フィルタ7に対して計測起動命令を発すると、これらは指定された計測が終了するまでデータ収集を行い、終了すると計算機8−1に割り込みをかけてデータ収集を通知する。データ計測制御はこのようにしてデータ収集終了を待ち合わせる(S−15−5)。計測されたデータは基準空間データを読み込み(S−15−6)、基準空間データで定義されている特徴パラメータを計測データから自動的に取り出し(S−15−7)、マハラノビスの距離を計算して(S−15−8)、被検者情報・計測条件・計測データおよびマハラノビスの距離を関係付けて心磁データベースに保存する。
【0085】
図17は図15のステップS−14におけるデータ解析のフローを示す。データ解析はいろいろな種類の波形や線図を表示して診断に必要な情報を得ようとするもので、図13のメニューを選択することによりいろいろな種類の波形や線図の画面を選択的に表示することができる。すなわち、「データ解析(A)」の「時間波形表示(W)」を選択すればステップS−14−4から時間波形画面が表示されるが、そのレイアウトは操作領域部を除けば図21に示す時間波形モニタと同じなのでここでは図示しない。「データ解析(A)」の「等磁線図(B)」を選択すればステップS−14−5から図22に示される等磁線図画面が表示される。「データ解析(A)」の「時間積分図(T)」を選択すればステップS−14−6から図23に示される時間積分図画面がそれぞれ表示される。そして「データ解析(A)」の「サマリ表示(S)」を選択すればステップS−14−7から図24に示されるサマリ表示画面がそれぞれ表示される。
【0086】
また、「ファイル(F)」の「心磁システムの終了(X)」を選択すればシステムが終了する。
【0087】
それぞれの画面において、ツールバーにあるアイコンボタン(808−1〜
808−15)をクリックすればそのクリックによって指定された波形又は線図の画面が代わって表示される。図17において、分岐の部分を「メニューで分岐」とせずに「メニュー又はアイコンボタンで分岐」(S−14−1)としたのはそのためである。したがって、この実施例によれば、前記図13のメニューを選択することなく、前記操作領域部にあるラジオボタンをクリックするだけで多様な解析データが得られるので、操作時間の短縮が図れるとともに、誤操作を軽減して操作性を向上することができる。
【0088】
データ解析を行う画面(図21〜図24)において、表示パラメータの一部を変更した場合、解析データを再計算して表示するのは好ましくない場合がある。例えば、等磁線図画面(図22)においてR波の近傍における等磁線図を表示しているが、続いてP波における磁場分布を調べるために等磁線図を表示する時には、表示時間を変更するだけでなく、カラーマップにおける配色を変更することが望ましい。何故ならば、R波における磁場強度の方がP波における磁場強度に比べて強いため、R波におけるカラーマップの配色をそのままP波における等磁線図に適用しようとすれば全体的に薄い色調でメリハリがなく、興奮部位を読み取りにくくなってしまう。また、所望するカラーマップに変更するまでパラメータを更新するたびに、画像の表示を更新するのは再描画の繰り返しによる応答時間の低下を招き、また意味のない画像を操作者に表示することになり、誤診に導く恐れがある。本実施例では、データ解析を行う画面に対して更新ボタン135およびキャンセルボタン134を設けることによってこれを解決している。すなわちデータ解析パラメータが変更され、データ解析パラメータが示している内容と解析データ部805−1に表示されている内容の対応がなくなった時に、更新ボタン135とキャンセルボタン134を活性状態とし、操作者が更新ボタン135を押した場合には更新されたデータ解析パラメータにしたがって再計算を行って解析データ部805−1の表示を更新する。またキャンセルボタン134が押された時にはデータ解析パラメータの変更を元に戻してデータ解析パラメータの内容が解析データ部805−1の内容と一致するようにする。すなわち、更新ボタン135とキャンセルボタン134は、データ解析パラメータの内容と解析データ部の表示内容との対応が保たれている間は不活性状態にあり、対応がなくなった時に活性化される。これにより、データ解析パラメータと解析データ部805−1の表示内容とが対応していない場合に、更新ボタンとキャンセルボタンの状態によって明確に区別できるため、解析データを誤って解釈する可能性が低下する。
【0089】
図22の等磁線図において、解析データ部の右端には縦に細長の磁場強度指標ボックス310が配置されている。その磁場強度指標ボックスは互いに色の異なる区画に区切られている。これは、等磁線図画面の各縞模様で示される磁場の強度範囲を色の種類で区別することで視覚的な(色彩的な)認識性を向上させるようにしたものである。すなわち、その磁場強度指標ボックス310の長手方向の中心位置311は磁場強度がゼロの位置で、その中心位置よりも上方の区画を中心位置に近い順番に第1〜第6区画とそれぞれ呼ぶことにすれば、たとえば第1区画は0〜2pTの磁場強度範囲に、第2区画は2〜4pTの磁場強度範囲に、第3区画は4〜6pTの磁場強度範囲に、第4区画は6〜8pTの磁場強度範囲に、第5区画は8〜10pTの磁場強度範囲に、そして第6区画は10〜12pTの磁場強度範囲にそれぞれ対応している。中心位置よりも下方の区画についてもまったく同じである。ただし、中心位置よりも上方の区画はプラス方向の磁場強度を、下方の区画はマイナス方向の磁場強度を表している。
【0090】
図22に示される等磁線図は、磁場強度指標ボックス310内の磁場強度範囲と色との所定の対応関係にしたがい、磁場強度に応じて色分け表示される。なお、色として、磁場強度のプラス側を暖色系、マイナス側は寒色系とし、中心部を黄色とするようにしてもよい。これにより、磁場の強弱を色彩的に認識することができるので視認性を向上することができる。しかも、この実施例によれば、解析データ部の近傍に磁場強度指標ボックス310を設けたので、比較対象の色、すなわち、マップに付された色と磁場強度指標ボックス310の所定の色とを目線移動を大きく移動させることなく比べながら確認することができるので、前記磁場の強弱のレベルと色との関係を明確に判断することができる。
【0091】
図22において、「再構成パラメータ」ボックス内の「マップ数」とあるのは表示される等磁線図の数を、「最大値」とあるのは磁場強度指標ボックス310の両端部に相当する磁場強度を、「間隔」とあるのは磁場強度指標ボックス310内の各区画の長さに対応する磁場範囲を意味する。その値については、これを対応するテキストボックスの三角又は逆三角ボタンをクリックして選択することができる。
【0092】
解析データ部の最下段には参照チャンネルの波形とマップ数として設定された16枚の等磁線図に対応した16本のカーソル線140が表示されている。カーソルの間隔は所定の時間で等間隔に設定されているが、間隔を設定するための手段を画面上に設けても良く、また複数の等磁場線図に対応したカーソル線を非等間隔で1本ずつ設定するための手段を設けても良い。カーソル線の位置についてはマウスでドラッグして左右に移動するが、テキストボックスから1枚目の等磁線図を作成する時刻を指定しても良い。図22では、表示されている等磁線図の数は16個であるが、これらの線図は波形上のカーソル線が位置する時点での線図であり、各マップについては、該マップがいつの時点のものであるかがわかるようにその時刻も表示される。
【0093】
これにより、図22で説明したと同様に、操作者は、現在解析データ部に表示されているマップが解析時間(参照波形の幅)の中でどのくらいの範囲を占めるかを示し、前記マップが示す範囲が解析時間の中でどこの範囲なのか等を一目で把握することができるので、視認性を向上させることができる。また、前記マップが示す範囲を2個のカーソルをマウス操作で簡単に移動させることで設定できるので操作が容易である。更に、各分割線の間隔を自由に設定するようにすれば、疑問のある部分を密にして他の部分を疎にする等、操作者に多様な解析環境を提供できる。
【0094】
図23の時間積分図において、データ表示領域には3つの等高線図が表示されているが、これはそれぞれ100−140ミリ秒における時間積分値による等高線図、180−240ミリ秒における時間積分値による等高線図、および両時間積分値の差による等高線図を表している。100−140ミリ秒はP波が発生している時間帯であり、180−240ミリ秒はQRS波が発生している時間帯であり、それぞれ参照チャンネルの波形において帯状カーソルによってその時間帯が示されている。通常はQRS波とT波の時間積分図および時間積分値の差を求めることが多い。
【0095】
図24のサマリ表示図は、特に心臓病の診断に重要と考えられるP波ピーク付近,QRS波,T波ピーク付近の等磁場線図を表示したものである。それぞれP波,R波,T波のピークとその前後1枚ずつの等磁場線図を表示している。これらは図7に示した特徴パラメータの中の、多くの情報を含んでいる。
【0096】
最後に、図26に疾患ごとの有効なパラメータを示す。他の検査によって心臓病と診断された患者を疾患名ごとに層別し、疾患ごとに要因効果図を作成して特異な値をとる特徴パラメータを明らかにしたものである。症例数は非常に少ないが、疾患ごとの有効なパラメータは顕著であった。これより、疾患と有効な特徴パラメータの関係をデータとして保持しておけば、被検者の計測データのマハラノビス距離が大きな値を示した場合に、特異な値を示す特徴パラメータから病名を推定できることを示す。
【0097】
図27はニューラル・ネットワークを用いた場合のブロック図を示している。マハラノビス距離の場合のブロック図(図11)と比べると、マハラノビス距離計算部が診断処理部,基準空間データがニューラル・ネットワーク,基準空間登録部がニューラル・ネットワーク登録部に置き換わっている。
【0098】
ニューラル・ネットワークは脳の情報処理方式を計算機で模倣したものである。情報処理単位であるニューロン(神経細胞)をモデル化したものを図27に、またニューロンの集合であるニューラル・ネットワークの例として階層型ニューラル・ネットワークを図28にそれぞれ1例として示す。
【0099】
図28のニューロンは、他のn個のニューロンからの入力Xj(j=1,…,n)があり、重み係数Wj によって重み付けされ、Xj*Wjとしてニューロンに入力されるものとする。ここで入力Xj は0から1までの間の値をとるものとする。ニューロンに入力された重み付けされた入力は最初にその総和を取り、単一のスカラー量に変換される(X1*W1+X2*W2+X3*W3+ … +Xn*Wn)。この総和を入力として出力関数fによって、そのニューロンの出力値が決定される。出力関数には次に示すしきい値関数やシグモイド関数,線形関数などが広く使われているが、何れも単調増加関数であり、重み付けされた総和の値域に対して0〜1の範囲の値を返す。また雑音を考慮して、これらの関数に対して確率的に入力値に変動を与えることもある。
しきい値関数
【0100】
【数3】
【0101】
シグモイド関数(1例としてボルツマン関数)
【0102】
【数4】
【0103】
線形関数
【0104】
【数5】
f(x)=ax
上述のニューロンを組合せて図29のようなニューラル・ネットワークによって心臓疾患の診断を行うことができる。入力層の各ニューロンには心臓磁気計測装置で計測された信号データから抽出した特徴パラメータ(例えば図9に示すもの)を割り当てればよい。ここで特徴パラメータは0〜1の間に規格化されているものとする。また出力層には考えられる診断結果が割り当てられているものとする。例えば健常者と心臓病患者とを区別するだけならば健常者のニューロンと心臓病患者のニューロンの2つがあるのみであり、診断処理の結果、2つのニューロンのうちのどちらかが発火する(すなわち1となる)。また、心臓病患者についてはその病名までを診断する必要があれば、出力層のニューロンは健常者と、心筋梗塞,狭心症,WPW症候群などの病名を表すニューロンからなる。中間層は診断処理の途中結果を保持するニューロンであり、必要に応じて最適な数の中間層を設けてもよい。
【0105】
【発明の効果】
本発明によれば、複数の計測位置の磁場強度の測定が容易で操作性の良好な生体の磁場の計測方法が提供される。
【図面の簡単な説明】
【図1】本発明が実施される生体磁場計測装置の一実施例の概略構成図。
【図2】図1の生体磁場計測装置に用いられる磁気センサの配置構成を示す斜視図。
【図3】図1の生体磁場計測装置において用いられる、磁場の法線成分を検出する磁気センサ単体の斜視図、および磁場の接線成分を検出する磁気センサ単体の斜視図。
【図4】図1の生体磁場計測装置における磁気センサと被検者の胸部との位置関係を示す図。
【図5】図1の生体磁場計測装置において、健常者について計測された特定の1チャンネルの心磁の接線成分の時間波形を示す図。
【図6】図1の生体磁場計測装置において、健常者について計測された心磁データから作成したP波,R波,T波のピーク近傍における等磁線図の例。
【図7】図1の生体磁場計測装置において、心筋梗塞患者について計測されたQRS波における等磁線図の例。
【図8】図1の生体磁場計測装置において、右脚ブロックについて計測されたQRS波における等磁線図の例。
【図9】心磁信号を定量評価するための特徴パラメータを表す表。
【図10】健常者および心疾患患者について得られたマハラノビスの距離の分布図。
【図11】図1の生体磁場計測装置において、計算機8−1で実行されるプログラムシステムの機能ブロック図。
【図12】図1の生体磁場計測装置において、ディスプレイ部に表示される表示画面の基本的なレイアウトを示す図。
【図13】図1の生体磁場計測装置において、特権ユーザモードでディスプレイ部に表示される表示画面のメニューバー部における操作メニューを示す図。
【図14】図1の生体磁場計測装置のディスプレイ部に表示される表示画面の基本的なレイアウトのうち、データ情報表示領域のレイアウト図。
【図15】図1の生体磁場計測装置において行われる全体の操作のフローを示す図。
【図16】図13の操作フロー中のデータ計測ステップにおけるデータ計測のフローを示す図。
【図17】図14の操作フロー中のデータ解析ステップにおけるデータ解析のフローを示す図。
【図18】図1の生体磁場計測装置のディスプレイ部に表示される被検者リスト画面を示す図。
【図19】図1の生体磁場計測装置において、ディスプレイ部に表示される表示画面中の操作メニューとして「被検者リスト(L)」−「被検者登録(R)」が選択された場合に開かれる被検者登録ダイアログボックスの内容を示す図。
【図20】図1の生体磁場計測装置において、ディスプレイ部に表示される表示画面中の操作メニューとして「被検者リスト(L)」−「被検者登録(R)」が選択された場合に開かれる被検者登録ダイアログボックスの内容を示す図。
【図21】図1の生体磁場計測装置のディスプレイ部に表示したデータ計測画面を示す図。
【図22】図1の生体磁場計測装置のディスプレイ部に表示される等磁線図表示画面を示す図。
【図23】図1の生体磁場計測装置のディスプレイ部に表示される時間積分図表示画面を示す図。
【図24】図1の生体磁場計測装置のディスプレイ部に表示されるサマリ表示画面を示す図。
【図25】特徴パラメータの要因効果図。
【図26】疾患ごとの有効なパラメータ。
【図27】ニューラル・ネットワークを用いた場合のブロック図。
【図28】情報処理単位であるニューロン(神経細胞)をモデル化したもの。
【図29】階層構造ニューラル・ネットワークの模式図。
【符号の説明】
1…磁気シールドルーム、2…被検者、3…ベッド、4…デユワ、5…自動補給装置、6…FLL回路、7…増幅器・フィルタ・増幅器、8…計算機、8−1…ディスプレイ部、8−2…キーボード、8−3…マウス、12…SQUIDセンサ、10…検出用コイル、11…参照用コイル、20−1〜20−8,21−1〜21−8,22−1〜22−8,23−1〜23−8,24−1〜24−8,25−1〜25−8,26−1〜26−8,27−1〜27−8…磁気センサ、10,10′,10″,11,11′,11″…コイル、12,12′,12″…SQUID、13,14…センサ、30…胸部、90…被検者情報並べ替えボタン、91…被検者情報カーソル、92…データ情報カーソル、93…データリスト選択タブ、100…計測/計測準備ボタン、101…解除ボタン、102…計測条件詳細設定ボタン、103…被検者情報詳細表示ボタン、121…アベレージング開始時刻カーソル、122…アベレージング基準時刻線、123…アベレージング終了時刻カーソル、124…時間波形の表示時間を表すスクロールボックス、126−1…時間波形の1/2頁送りボタン、126−2…頁送りボタン、127…しきい値カーソル、128…アベレージング処理開始ボタン、131…基線補正開始時刻カーソル、132…基線補正終了時刻カーソル、133…信号基線、134…表示更新キャンセルボタン、135…表示更新ボタン、136…チャンネル選択ボタン、140…1枚目の等磁線図等磁線図の時刻カーソル、142…センサ位置、261…スクロールボックス、262…スクロールバー、271…しきい値カーソル、273〜275…カーソル線、311,312,331〜334…カーソル、801…タイトルバー部、802…メニューバー部、803…ツールバー部、804…被検者情報部、805…データ情報部、805−1…解析データ部、805−2…参照波形部、806…操作領域部、808−1〜15…アイコン。
Claims (16)
- 被検者の心臓から発せられる磁場を複数の位置で計測する心臓磁場計測装置において、
健常者の集団から得た計測信号情報からなる基準データベースを記憶させた基準データベース記憶手段と、
新たに計測した信号情報を記憶する計測信号情報記憶手段と、
該計測信号情報記憶手段に記憶された信号情報から少なくともR波ピーク近傍の一定時間ごとの電流方向を含む複数の特徴パラメータを算出する特徴パラメータ演算手段と、
該演算手段により算出された特徴パラメータを、前記基準データベース記憶手段に記憶された基準データベースをもとに健常者の集団の特徴パラメータと比較する比較手段を備え、
更に該比較手段による比較結果に基づき被検者が健常であるかどうかを識別するための識別情報を算出する識別情報算出手段を備えたことを特徴とする心臓磁場計測装置。 - 請求項1記載の心臓磁場計測装置において、
前記特徴パラメータ演算手段で信号情報から算出される特徴パラメータは、更に
計測した磁場強度のうちの最も振幅の大きな時刻を代表時刻として、該時刻での等磁場線図に基づいて抽出されたパラメータを含むことを特徴とする心臓磁場計測装置。 - 請求項1または2に記載の心臓磁場計測装置において、
識別情報算出手段が、前記特徴パラメータ演算手段により算出された特徴パラメータに基づきマハラノビス空間を構築するマハラノビス空間構築手段と、
該マハラノビス空間構築手段により構築されたマハラノビス空間における前記被検者の、複数の特徴パラメータに基づいたマハラノビスの距離を計算するマハラノビスの距離計算手段と、
該マハラノビスの距離計算手段により得られたマハラノビスの距離の情報に基づいて被検者が健常であるかどうかを識別するための識別情報を算出する識別情報算出手段であることを特徴とする心臓磁場計測装置。 - 請求項1〜3のいずれかに記載の心臓磁場計測装置において、
前記識別情報算出手段により算出された識別情報と、被検者が健常であるかどうかを識別するため予め記憶された識別基準情報とを比較し、被検者が健常であるかどうかを識別する識別手段を備えたことを特徴とする心臓磁場計測装置。 - 請求項4記載の心臓磁場計測装置において、
前記識別基準情報を操作者が任意に設定可能であることを特徴とする心臓磁場計測装置。 - 請求項4記載の心臓磁場計測装置において、
前記識別手段により被検者が健常でないと識別されたときに、前記複数の特徴パラメータに基づき、疾患の候補を選択する疾患推定機能を備えたことを特徴とする心臓磁場計測装置。 - 請求項3〜6のいずれかに記載の心臓磁場計測装置において、
複数の被検者毎の前記特徴パラメータ,マハラノビスの距離を記憶する記憶手段を備え、
選択した被検者の少なくとも1つの特徴パラメータ,マハラノビスの距離を表示する機能を備えたことを特徴とする心臓磁場計測装置。 - 請求項3〜6のいずれかに記載の心臓磁場計測装置において、
複数の被検者毎の前記特徴パラメータ,マハラノビスの距離を記憶する記憶手段を備え、
選択した被検者の少なくとも1つの特徴パラメータ,マハラノビスの距離と、
前記構築されたマハラノビス空間におけるマハラノビスの距離及び特徴パラメータを同時に表示する機能を備えたことを特徴とする心臓磁場計測装置。 - 請求項7または8記載の心臓磁場計測装置において、
選択した被検者の特徴パラメータと、上記パラメータのマハラノビス空間を定義した健常者の集団における平均値と標準偏差とを表示する機能を備えたことを特徴とする心臓磁場計測装置。 - 請求項7または8記載の心臓磁場計測装置において、
選択した被検者の特徴パラメータと基準化した特徴パラメータを併記する機能を備えたことを特徴とする心臓磁場計測装置。 - 請求項7または8記載の心臓磁場計測装置において、
マハラノビスの距離が、被検者が健常であるかどうかを識別するため予め記憶された識別基準情報以上となるデータを有する被検者のデータを強調表示する機能を備えたことを特徴とする心臓磁場計測装置。 - 請求項7または8記載の心臓磁場計測装置において、
マハラノビスの距離が、被検者が健常であるかどうかを識別するため予め記憶された識別基準情報以上となるデータの時に、基準化した特徴パラメータが予め定められたしきい値以上の特徴パラメータを強調表示する機能を備えたことを特徴とする心臓磁場計測装置。 - 請求項7または8記載の心臓磁場計測装置において、
記憶された複数の被検者のデータを、マハラノビスの距離の大きさによって表示順を昇順あるいは降順に並び替える機能を備えたことを特徴とする心臓磁場計測装置。 - 請求項3記載の心臓磁場計測装置において、
複数のマハラノビス空間データを記憶する記憶手段を備え、被検者情報に基づいて使用する該空間データを選択できる機能を備えたことを特徴とする心臓磁場計測装置。 - 請求項14記載の心臓磁場計測装置において、
前記空間データごとに特徴パラメータを定義できる機能を備えたことを特徴とする心臓磁場計測装置。 - 請求項1〜3のいずれかに記載の心臓磁場計測装置において、
被検者の生体内から発せられる磁場が主に心臓から発せられる心臓磁場であり、かつ心臓磁場信号情報のうちの、P波・QRS波・T波近傍の等磁線図を1画面上に表示する機能を備えたことを特徴とする心臓磁場計測装置。
Priority Applications (4)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2001350920A JP3944383B2 (ja) | 2001-11-16 | 2001-11-16 | 心臓磁場計測装置 |
US10/255,037 US20030097056A1 (en) | 2001-11-16 | 2002-09-26 | Biomagnetic field measuring apparatus |
EP02022202A EP1312304B1 (en) | 2001-11-16 | 2002-10-01 | Biomagnetic field measuring apparatus |
DE60228670T DE60228670D1 (de) | 2001-11-16 | 2002-10-01 | Messgerät für ein biomagnetisches Feld |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2001350920A JP3944383B2 (ja) | 2001-11-16 | 2001-11-16 | 心臓磁場計測装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2003144406A JP2003144406A (ja) | 2003-05-20 |
JP3944383B2 true JP3944383B2 (ja) | 2007-07-11 |
Family
ID=19163316
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2001350920A Expired - Lifetime JP3944383B2 (ja) | 2001-11-16 | 2001-11-16 | 心臓磁場計測装置 |
Country Status (4)
Country | Link |
---|---|
US (1) | US20030097056A1 (ja) |
EP (1) | EP1312304B1 (ja) |
JP (1) | JP3944383B2 (ja) |
DE (1) | DE60228670D1 (ja) |
Families Citing this family (28)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6473518B1 (en) * | 1997-10-02 | 2002-10-29 | Hitachi, Ltd. | Method of measuring a biomagnetic field, method of analyzing a measured biomagnetic field, method of displaying biomagnetic field data, and apparatus therefor |
JP3518502B2 (ja) * | 2000-10-19 | 2004-04-12 | 株式会社日立製作所 | 生体磁場計測装置 |
WO2003073931A1 (fr) * | 2002-03-05 | 2003-09-12 | Dainippon Pharmaceutical Co., Ltd. | Appareil de carte d'electrocardiographie et procede associe |
JP3736640B2 (ja) * | 2004-01-19 | 2006-01-18 | 株式会社エイティング | 人身監視システム |
JP2005285001A (ja) * | 2004-03-30 | 2005-10-13 | Sunstar Inc | 健康支援方法及びそのシステム |
JP4922564B2 (ja) * | 2005-03-16 | 2012-04-25 | 株式会社リコー | 異常判定装置及び画像形成装置 |
JP4641492B2 (ja) * | 2005-11-16 | 2011-03-02 | 日本電信電話株式会社 | 特異映像検出装置、特異映像検出方法およびプログラム |
EP1818651A1 (en) * | 2006-02-10 | 2007-08-15 | Rijksuniversiteit Groningen | System and a method for determining one or more parameters of a source of a potential-energy field |
JP4972443B2 (ja) * | 2007-03-30 | 2012-07-11 | 株式会社日立ハイテクノロジーズ | 生体磁場計測装置 |
US20090118632A1 (en) * | 2007-11-05 | 2009-05-07 | Goepp Julius G | Effort-Independent, Portable, User-Operated Capnograph Devices And Related Methods |
KR100980331B1 (ko) * | 2008-05-28 | 2010-09-06 | 한국표준과학연구원 | 심자도의 파라미터 조합을 이용한 심근경색 진단법 |
JP5611546B2 (ja) * | 2009-06-30 | 2014-10-22 | 株式会社東芝 | 自動診断支援装置、超音波診断装置及び自動診断支援プログラム |
JP5746480B2 (ja) * | 2010-05-13 | 2015-07-08 | 横河電機株式会社 | 解析結果表示装置および解析結果表示方法 |
US9015093B1 (en) | 2010-10-26 | 2015-04-21 | Michael Lamport Commons | Intelligent control with hierarchical stacked neural networks |
US8775341B1 (en) | 2010-10-26 | 2014-07-08 | Michael Lamport Commons | Intelligent control with hierarchical stacked neural networks |
UA104073C2 (uk) * | 2012-07-13 | 2013-12-25 | Илья Анатольевич Чайковский | Спосіб оцінки ступеня ушкодження міокарда на основі аналізу змін у часі показників щільності струму |
US9500464B2 (en) * | 2013-03-12 | 2016-11-22 | Adidas Ag | Methods of determining performance information for individuals and sports objects |
US9684069B2 (en) * | 2014-01-10 | 2017-06-20 | Raytheon Company | High angular resolution low frequency radar with small antenna and high resolution low frequency ground penetrating radar |
US10772520B2 (en) | 2015-06-25 | 2020-09-15 | DePuy Synthes Products, Inc. | Intraoperative magnetometry monitoring system |
US11191506B2 (en) * | 2016-07-08 | 2021-12-07 | Ricoh Company, Ltd. | Diagnosis support system, diagnosis support apparatus, and recording medium |
JP6890448B2 (ja) * | 2017-03-30 | 2021-06-18 | 株式会社日立ハイテク | 磁場データ処理システム |
US11134877B2 (en) | 2017-08-09 | 2021-10-05 | Genetesis, Inc. | Biomagnetic detection |
US20200269128A1 (en) * | 2017-09-28 | 2020-08-27 | Neofect Co., Ltd. | Pegboard training method and program therefor |
WO2020138170A1 (ja) | 2018-12-26 | 2020-07-02 | 旭化成エレクトロニクス株式会社 | 磁場計測装置 |
US11585869B2 (en) * | 2019-02-08 | 2023-02-21 | Genetesis, Inc. | Biomagnetic field sensor systems and methods for diagnostic evaluation of cardiac conditions |
JP2021060378A (ja) * | 2019-10-02 | 2021-04-15 | 株式会社リコー | 磁場計測装置 |
WO2022178314A1 (en) * | 2021-02-22 | 2022-08-25 | Genetesis, Inc. | Biomagnetic field sensor systems and methods for diagnostic evaluation of cardiac conditions |
JP2022168642A (ja) | 2021-04-26 | 2022-11-08 | Tdk株式会社 | 生体信号処理システムおよび生体信号計測システム |
Family Cites Families (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4974598A (en) * | 1988-04-22 | 1990-12-04 | Heart Map, Inc. | EKG system and method using statistical analysis of heartbeats and topographic mapping of body surface potentials |
US4951674A (en) * | 1989-03-20 | 1990-08-28 | Zanakis Michael F | Biomagnetic analytical system using fiber-optic magnetic sensors |
US4977896A (en) * | 1989-05-26 | 1990-12-18 | Biomagnetic Technologies, Inc. | Analysis of biological signals using data from arrays of sensors |
JPH0767440B2 (ja) * | 1990-09-30 | 1995-07-26 | 工業技術院長 | 心磁図同期検出用の基準部分検出方法、心磁図同期加算方法、心磁図同期検出装置および心磁図同期加算装置 |
JPH0542119A (ja) * | 1991-08-12 | 1993-02-23 | Fujitsu Ltd | 生体磁気計測装置 |
JPH063427A (ja) * | 1992-06-23 | 1994-01-11 | Fujitsu Ltd | 磁気検出装置 |
JP4208275B2 (ja) * | 1997-10-30 | 2009-01-14 | 株式会社東芝 | 心臓内電気現象の診断装置およびその現象の表示方法 |
US6195576B1 (en) * | 1998-03-09 | 2001-02-27 | New York University | Quantitative magnetoencephalogram system and method |
US6217525B1 (en) * | 1998-04-30 | 2001-04-17 | Medtronic Physio-Control Manufacturing Corp. | Reduced lead set device and method for detecting acute cardiac ischemic conditions |
US6275719B1 (en) * | 1998-09-09 | 2001-08-14 | Hitachi, Ltd. | Biomagnetic field measurement apparatus |
WO2000037917A2 (en) * | 1998-12-23 | 2000-06-29 | Medispectra, Inc. | Systems and methods for optical examination of samples |
JP3712348B2 (ja) * | 2000-07-18 | 2005-11-02 | 独立行政法人科学技術振興機構 | 心房粗動および心房細動の心臓磁界診断装置およびその作動方法 |
US7106057B2 (en) * | 2004-03-17 | 2006-09-12 | University Of Maryland | High frequency scanning SQUID microscope and method of measuring high frequency magnetic fields |
-
2001
- 2001-11-16 JP JP2001350920A patent/JP3944383B2/ja not_active Expired - Lifetime
-
2002
- 2002-09-26 US US10/255,037 patent/US20030097056A1/en not_active Abandoned
- 2002-10-01 EP EP02022202A patent/EP1312304B1/en not_active Expired - Fee Related
- 2002-10-01 DE DE60228670T patent/DE60228670D1/de not_active Expired - Lifetime
Also Published As
Publication number | Publication date |
---|---|
EP1312304B1 (en) | 2008-09-03 |
JP2003144406A (ja) | 2003-05-20 |
US20030097056A1 (en) | 2003-05-22 |
DE60228670D1 (de) | 2008-10-16 |
EP1312304A3 (en) | 2004-01-14 |
EP1312304A2 (en) | 2003-05-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP3944383B2 (ja) | 心臓磁場計測装置 | |
JP5085155B2 (ja) | 生体磁場計測装置 | |
JP3997224B2 (ja) | 磁場計測装置の表示装置 | |
JP3563624B2 (ja) | 生体磁場計測装置 | |
JP3514655B2 (ja) | 生体磁場計測装置 | |
JP2004249124A (ja) | 生体磁場計測装置 | |
JP3595660B2 (ja) | 生体の磁場の計測装置 | |
JP3518495B2 (ja) | 生体磁場計測装置 | |
JP3722140B2 (ja) | 生体磁場計測装置 | |
JP3228702B2 (ja) | 生体の磁場の計測方法および磁場図表示装置 | |
JP3228703B2 (ja) | 生体が発する磁場の計測データの処理方法および装置 | |
JP4054011B2 (ja) | 磁場計測装置 | |
JP3720960B2 (ja) | 生体が発する磁場の計測データの処理方法 | |
JP3992705B2 (ja) | 磁場計測装置の表示装置 | |
JP3714346B2 (ja) | 生体磁場計測装置 | |
JP3518494B2 (ja) | 生体磁場計測装置 | |
JP2004089724A (ja) | 生体磁場計測装置 | |
JP3058619B2 (ja) | 生体の磁場図形成装置 | |
JP3768178B2 (ja) | 磁場図表示方法 | |
JP2004041788A (ja) | 生体磁場計測装置 | |
JP4044529B2 (ja) | 磁場計測装置の表示装置 | |
JP3140731B2 (ja) | 生体の磁場図表示装置 | |
JP3140735B2 (ja) | 生体の磁場図表示方法 | |
JPH11104101A (ja) | 生体の磁場図表示装置 | |
JP2004041786A (ja) | 生体磁場計測装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20040908 |
|
RD02 | Notification of acceptance of power of attorney |
Free format text: JAPANESE INTERMEDIATE CODE: A7422 Effective date: 20060512 |
|
RD04 | Notification of resignation of power of attorney |
Free format text: JAPANESE INTERMEDIATE CODE: A7424 Effective date: 20060512 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20060912 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20061109 |
|
TRDD | Decision of grant or rejection written | ||
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20070327 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20070409 |
|
R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 Ref document number: 3944383 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110413 Year of fee payment: 4 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120413 Year of fee payment: 5 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120413 Year of fee payment: 5 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130413 Year of fee payment: 6 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20140413 Year of fee payment: 7 |
|
S531 | Written request for registration of change of domicile |
Free format text: JAPANESE INTERMEDIATE CODE: R313531 |
|
S533 | Written request for registration of change of name |
Free format text: JAPANESE INTERMEDIATE CODE: R313533 |
|
R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
EXPY | Cancellation because of completion of term |