(57)【特許請求の範囲】
【請求項１】
生体表面に配置される電極に生じる電圧によって、生体内の信号源の位置を検出する生体内信号源位置検出装置であって、
生体表面に配置する少なくとも３個の電極と、
前記各電極とグランド電位との間に、第１の外部抵抗及び第２の外部抵抗をそれぞれ並列接続する接続手段と、
前記各電極を生体表面に配置した状態で、前記接続手段により、前記各電極に前記第１の外部抵抗を並列接続したときに各電極に生じる第１の出力電圧Ｖｉ（ｉ＝１，２，３）、及び各電極に前記第２の外部抵抗を並列接続したときに各電極に生じる第２の出力電圧Ｖ’ｉ（ｉ＝１，２，３）を測定する測定手段と、
前記第１の出力電圧Ｖｉ及び前記第２の出力電圧Ｖ’ｉの測定を１サイクルとして、繰り返し行い、前記第１の出力電圧Ｖｉ及び前記第２の出力電圧Ｖ’ｉの時系列な測定データから、前記第１の出力電圧Ｖｉ及び前記第２の出力電圧Ｖ’ｉの時間変化を近似する近似波形を取得する近似手段と、
各サイクルにおいて、測定時刻に時間差を有する前記第１の出力電圧Ｖｉ及び前記第２の出力電圧Ｖ’ｉに対して、前記近似手段で取得した近似波形を用いて、同一の測定時刻における前記第１の出力電圧Ｖｉ及び前記第２の出力電圧Ｖ’ｉの補正値Ｖｃｉ（ｉ＝１，２，３）及びＶｃ’ｉ（ｉ＝１，２，３）を算出する算出手段と、
サイクル毎に、前記算出された補正値Ｖｃｉ及びＶｃ’ｉの比Ｖｃｉ／Ｖｃ’ｉ（ｉ＝１，２，３）を求め、これら３つの比Ｖｃｉ／Ｖｃ’ｉ（ｉ＝１，２，３）に基づいて、生体内の信号源の位置を検出する検出手段と
を備えた、生体内信号源位置検出装置。
【請求項２】
前記接続手段は、前記各電極とグランド電位との間に、前記第１の外部抵抗及び第２の外部抵抗を、所定の切り替え時間で切り替えて並列接続する切替手段からなり、
前記算出手段は、各サイクルにおいて、測定時刻に時間差を有する前記第１の出力電圧Ｖｉ及び前記第２の出力電圧Ｖ’ｉに対して、前記近似手段で取得した近似波形を用いて、前記切り替え時間に応じた同一の測定時刻における前記第１の出力電圧Ｖｉ及び前記第２の出力電圧Ｖ’ｉの補正値Ｖｃｉ及びＶｃ’ｉを算出する、請求項１に記載の生体内信号源位置検出装置。
【請求項３】
前記算出手段は、各サイクルにおいて、測定時刻に時間差を有する前記第１の出力電圧Ｖｉ及び前記第２の出力電圧Ｖ’ｉに対して、前記近似手段で取得した各近似波形の特徴ポイントを合わせることによって、同一の測定時刻における前記第１の出力電圧Ｖｉ及び前記第２の出力電圧Ｖ’ｉの補正値Ｖｃｉ及びＶｃ’ｉを算出する、請求項１に記載の生体内信号源位置検出装置。
【請求項４】
生体表面に配置される電極に生じる電圧によって、生体内の信号源の位置を検出する生体内信号源位置検出方法であって、
生体表面に３個の電極を配置する配置工程と、
前記各電極とグランド電位との間に、第１の外部抵抗及び第２の外部抵抗を、それぞれ並列接続する接続工程と、
前記各電極を生体表面に配置した状態で、前記各電極に前記第１の外部抵抗を並列接続したときに各電極に生じる第１の出力電圧Ｖｉ（ｉ＝１，２，３）、及び各電極に前記第２の外部抵抗を並列接続したときに各電極に生じる第２の出力電圧Ｖ’ｉ（ｉ＝１，２，３）を測定する測定工程と、
前記第１の出力電圧Ｖｉ及び前記第２の出力電圧Ｖ’ｉの測定を１サイクルとして、繰り返し行い、前記第１の出力電圧Ｖｉ及び前記第２の出力電圧Ｖ’ｉの時系列な測定データから、前記第１の出力電圧Ｖｉ及び前記第２の出力電圧Ｖ’ｉの時間変化を近似する近似波形を取得する取得工程と、
各サイクルにおいて、測定時刻に時間差を有する前記第１の出力電圧Ｖｉ及び前記第２の出力電圧Ｖ’ｉに対して、前記近似波形を用いて、同一の測定時刻における前記第１の出力電圧Ｖｉ及び前記第２の出力電圧Ｖ’ｉの補正値Ｖｃｉ（ｉ＝１，２，３）及びＶｃ’ｉ（ｉ＝１，２，３）を算出する算出工程と、
サイクル毎に、前記算出された補正値Ｖｃｉ及びＶｃ’ｉの比Ｖｃｉ／Ｖｃ’ｉ（ｉ＝１，２，３）を求め、これら３つの比Ｖｃｉ／Ｖｃ’ｉ（ｉ＝１，２，３）に基づいて、生体内の信号源の位置を検出する検出工程と
を含む、生体内信号源位置検出方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、生体内の信号源の位置を検出する生体内信号源位置検出装置、及び生体内信号源位置検出方法に関する。
【背景技術】
【０００２】
生体表面に電極を取り付けて、その電極に生じた電圧を測定することによって、心電図等の生体内の電気活動を計測することが一般に行われている。
【０００３】
例えば、特許文献１には、生体と所定平面との交線（閉曲線）上の各点において、表面電位を測定することによって、その平面を通る生体の断面内の電位分布を求める方法が開示されている。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開平１１−１１３８６７号公報
【特許文献２】特許第５１７８９０９号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
しかしながら、特許文献１に開示された方法では、生体に多数の電極を隙間なく配置して測定する必要があるため、生体への負担が大きい。また、電極数が少なければ、生体への負担は軽減されるが、解像度の低い電位分布しか得ることができない。
【０００６】
本願出願人は、このような課題を解決するために、少ない数の電極を用いて、生体内の信号源の位置を検出することができる方法を、特許文献２及びＰＣＴ／ＪＰ２０１４／００５７３２の出願明細書に開示している。ここで開示した方法は、生体表面に配置した電極とグランド電位との間に複数の抵抗を切り替えて接続し、各抵抗に接続したときに電極に生じる電圧を測定することによって、生体内の信号源の位置を検出するものである。
【０００７】
しかしながら、上記の検出方法では、各抵抗に接続したときに電極に生じる電圧を測定する際、抵抗を切り替える時間を要するため、電圧の測定時刻に時間差が生じてしまう。この場合、生体内の信号源からの電圧変化が大きいと、生体内の信号源の位置検出に誤差が生じてしまうという問題がある。
【０００８】
本発明は、上記課題に鑑みなされたもので、その主な目的は、少ない数の電極を用いて、生体内の信号源の位置を精度よく検出することができる生体内信号源位置検出装置、及び生体内信号源位置検出方法を提供することにある。
【課題を解決するための手段】
【０００９】
本発明に係る生体内信号源位置検出装置は、生体表面に配置される電極に生じる電圧によって、生体内の信号源の位置を検出する生体内信号源位置検出装置であって、
生体表面に配置するｍ個の電極と、
各電極とグランド電位との間に、ｎ個の外部抵抗をそれぞれ並列接続する接続手段と、
各電極を生体表面に配置した状態で、接続手段により、各電極に各外部抵抗を並列接続したときに各電極に生じるｍ×ｎ個の出力電圧Ｖｉ（ｉは１〜ｍ×ｎの整数）を測定する測定手段と、
ｍ×ｎ個の出力電圧Ｖｉの測定を１サイクルとして、繰り返し行い、各出力電圧Ｖｉの時系列な測定データから、各出力電圧Ｖｉの時間変化を近似する近似波形を取得する近似手段と、
各サイクルにおいて、測定時刻に時間差を有する各出力電圧Ｖｉに対して、近似手段で取得した近似波形を用いて、同一の測定時刻における各出力電圧Ｖｉの補正値Ｖｃｉ（ｉは１〜ｍ×ｎの整数）を算出する算出手段と、
サイクル毎に、算出された各補正値Ｖｃｉに基づいて、生体内の信号源の位置を検出する検出手段とを備えたことを特徴とする。
【００１０】
本発明に係る生体内信号源位置検出方法は、生体表面に配置される電極に生じる電圧によって、生体内の信号源の位置を検出する生体内信号源位置検出方法であって、
生体表面にｍ個の電極を配置する配置工程と、
各電極とグランド電位との間に、ｎ個の外部抵抗を、それぞれ並列接続する接続工程と、
各電極を生体表面に配置した状態で、各電極に各外部抵抗を並列接続したときに各電極に生じるｍ×ｎ個の出力電圧Ｖｉ（ｉは１〜ｍ×ｎの整数）を測定する測定工程と、
ｍ×ｎ個の出力電圧Ｖｉの測定を１サイクルとして、繰り返し行い、各出力電圧Ｖｉの時系列な測定データから、各出力電圧Ｖｉの時間変化を近似する近似波形を取得する取得工程と、
各サイクルにおいて、測定時刻に時間差を有する各出力電圧Ｖｉに対して、近似波形を用いて、同一の測定時刻における各出力電圧Ｖｉの補正値Ｖｃｉ（ｉは１〜ｍ×ｎの整数）を算出する算出工程と、
サイクル毎に、算出された各補正値Ｖｃｉに基づいて、生体内の信号源の位置を検出する検出工程とを含むことを特徴とする。
【発明の効果】
【００１１】
本発明によれば、少ない数の電極を用いて、生体内の信号源の位置を精度よく検出することができる生体内信号源位置検出装置、及び生体内信号源位置検出方法を提供することができる。
【図面の簡単な説明】
【００１２】
【図１】生体内の信号源の位置を検出する方法の概要を示したブロック図である。
【図２】生体内の信号源の位置を検出する方法を例示した図である。
【図３】各電極と外部抵抗との接続状態の切り替えステップを示した表である。
【図４】信号源の位置座標を求める方法を示した図である。
【図５】（ａ）、（ｂ）は、生体内信号源位置検出方法を用いて、心臓の電位を測定した結果のイメージ図である。
【図６】本発明の一実施形態における生体内信号源位置検出装置の構成を示した図である。
【図７】各電極と外部抵抗との接続状態の切り替えステップを示した表である。
【図８】各電極の出力電圧を、同一の測定時刻におけるデータに補正する方法を説明した図である。
【図９】各電極の出力電圧を、同一の測定時刻におけるデータに補正する方法を説明した図である。
【図１０】各電極の出力電圧を、同一の測定時刻におけるデータに補正する他の方法を説明した図である。
【図１１】各電極の出力電圧を、同一の測定時刻におけるデータに補正する他の方法を説明した図である。
【図１２】本発明の他の実施形態における生体内信号源位置検出装置の構成を示した図である。
【発明を実施するための形態】
【００１３】
本発明を説明する前に、本願出願人が、特許文献２に開示した、生体内の信号源の位置を検出する方法の概要を、図１を参照しながら説明する。
【００１４】
図１に示すように、生体１０の表面に、測定電極２１と接地電極２０とが配置されている。測定電極２１に生じる電圧は、アンプ４０で増幅されて、出力電圧Ｖoutとして測定される。測定電極２１と接地電極２０との間には、インピーダンス切替手段３０が配置されている。ここでは、複数の抵抗素子Ｚ１〜Ｚｎが並列に配置され、スイッチＳＷによって、抵抗素子Ｚ１〜Ｚｎのうち１つが選択される。そして、測定電極２１に選択された抵抗素子が並列接続されたときに、測定電極２１に生じる電圧Ｖoutが測定される。
【００１５】
このような方法により、スイッチＳＷで、ｎ個の抵抗素子Ｚ１〜Ｚｎを切り替えることによって、時刻ｔにおけるｎ個の測定データが得られる。なお、測定電極２１は、１個に限らず、複数であってもよい。その場合は、ｍ個の測定電極２１と、ｎ個の抵抗素子との組み合わせにより、ｍ×ｎ個の測定データが得られる。これにより、比較的少ない測定電極を用いて、生体１０内の信号源の位置検出に必要な数の測定データを得ることができる。
【００１６】
生体１０内の信号源の位置は、抵抗素子Ｚ１〜Ｚｎを切り替えることによって得られた時刻ｔにおける測定データを用いて、例えば、以下の方法により検出することができる。
【００１７】
測定電極２１を配置した位置を含む生体１０の断層画像を取得し、この断層画像を格子状に分割して、各格子点の周りに、該当する組織に対するアドミッタンスを配置し、インピーダンス切替手段３０で切り替えたインピーダンスを含む回路網を生成する。そして、この回路網の何れかの格子点と接地との間に信号源を配置した条件で、生体１０内の信号源の位置を求めることができる。さらに、時系列に測定した測定データを用いることによって、生体１０内の信号源の時間変化を検出することができる。
【００１８】
図２〜図５は、本願出願人が、ＰＣＴ／ＪＰ２０１４／００５７３２の出願明細書に開示した、生体内の信号源の位置を検出する方法を説明した図である。
【００１９】
図２に示すように、生体１０の表面に、３つの電極２１、２２、２３を配置する。また、各電極２１、２２、２３とグランド電極２０との間に、相互に切り替え可能な第１の外部抵抗及び第２の外部抵抗を並列接続する。ここでは、第１の外部抵抗の抵抗値を無限大とし、第２の外部抵抗の抵抗値をＲgとしている。これにより、各電極２１、２２、２３とグランド電極２０との間は、切り替えスイッチＳＷによって、外部抵抗が接続されていない場合と、外部抵抗Ｒgが接続されている場合とに切り替えられる。ここでは、生体１０の表面に、グランド電極２０を配置し、これをグランド電位にしているが、必ずしも、生体１０の表面に、グランド電極２０を配置する必要はない。
【００２０】
生体１０の表面に配置した各電極２１、２２、２３には、生体１０内の信号源Ｖsからの電圧が発生し、その電圧をアンプ４０で増幅して、出力電圧Ｖoutが出力される。各電極２１、２２、２３とアンプ４０との間には、スイッチＳ１、Ｓ２、Ｓ３がそれぞれ配置され、各スイッチＳ１、Ｓ２、Ｓ３を順次、導通させることにより、各電極２１、２２、２３に生じた電圧が出力電圧Ｖoutとして測定される。
【００２１】
図３に示すように、ステップ１〜ステップ６で、切り替えスイッチＳＷ、及びスイッチＳ１、Ｓ２、Ｓ３を、それぞれ切り替えることにより、各電極２１、２２、２３とグランド電位との間に、外部抵抗が接続されていないときに各電極２１、２２、２３に生じる第１の電圧Ｖ１、Ｖ２、Ｖ３、及び各電極２１、２２、２３とグランド電位との間に、外部抵抗Ｒgが接続されているときに各電極２１、２２、２３に生じる第２の電圧Ｖ’１、Ｖ’２、Ｖ’３が測定される。なお、図３では、各電極２１、２２、２３を、それぞれ、チャネルｃｈ１、ｃｈ２、ｃｈ３と表示している。
【００２２】
詳細な説明は省略するが、生体１０内の信号源Ｖsと、各電極２１、２２、２３との距離Ｌ１、Ｌ２、Ｌ３は、それぞれ、第１の電圧と第２の電圧との比Ｖ１／Ｖ’１、Ｖ２／Ｖ’２、Ｖ３／Ｖ’３を用いて、式（１）、（２）、（３）で表される。
【００２３】
【数１】
0006518084.tif 000002
【００２４】
【数２】
0006518084.tif 000003
【００２５】
【数３】
0006518084.tif 000004
【００２６】
ここで、βは、生体１０の導電率等で定まる定数である。また、Ｒｂ１、Ｒｂ２、Ｒｂ３は、それぞれ、生体１０内の信号源Ｖsと、電極２１、２２、２３との間の内部抵抗の抵抗値を表す。
【００２７】
信号源Ｖsは、図４（ａ）、（ｂ）に示すように、各電極２１、２２、２３を中心とする半径Ｌ１、Ｌ２、Ｌ３の球体Ｑ１、Ｑ２、Ｑ３の交点に存在すると考えられる。従って、信号源Ｖsの３次元的な位置座標（ｘ、ｙ、ｚ）は、３つの球体Ｑ１、Ｑ２、Ｑ３の式（４）、（５）、（６）を解くことによって求めることができる。ここで、各電極２１、２２、２３の位置座標を、（ａ１、ｂ１、ｃ１）、（ａ２、ｂ２、ｃ２）、（ａ３、ｂ３、ｃ３）としている。
【００２８】
【数４】
0006518084.tif 000005
【００２９】
【数５】
0006518084.tif 000006
【００３０】
【数６】
0006518084.tif 000007
【００３１】
このように、生体１０の表面に配置した３つの電極２１、２２、２３と、グランド電位との間に外部抵抗を並列接続し、その接続状態を切り替えて、各電極２１、２２、２３に生じる電圧の比Ｖｉ／Ｖ’ｉ（ｉ＝１，２，３）を測定することによって、生体１０内の信号源Ｖsの３次元位置を検出することができる。これにより、少ない数の電極を用いて、生体内の信号源Ｖsの３次元位置を精度よく検出することができる。
【００３２】
また、電圧比Ｖｉ／Ｖ’ｉ（ｉ＝１，２，３）の測定を１サイクルとして、繰り返し行うことにより、各サイクルにおける電圧比の時系列な測定データから、生体内の信号源Ｖsの３次元位置の移動軌跡をリアルタイムに検出することができる。
【００３３】
図５（ａ）、（ｂ）は、上記の方法を用いて、心臓の電位を測定した結果のイメージ図である。図５（ａ）は、各電極で測定された信号源の電圧波形（心電図）を示す。また、図５（ｂ）は、各電極での電圧比Ｖｉ／Ｖ’ｉ（ｉ＝１，２，３）に基づいて算出された信号源の３次元位置の移動軌跡を破線で示したものである。図５（ｂ）の破線で示した移動軌跡において、ポイントＰ１、Ｐ２、Ｐ３、Ｒ、Ｓ、Ｔは、それぞれ、図５（ａ）に示した電圧波形のＰ１、Ｐ２、Ｐ３、Ｒ、Ｓ、Ｔに対応する信号源の３次元位置を示している。図５（ｂ）に示すように、心臓の電気活動における信号源が、心房から心室へと移動している様子をリアルタイムに検出することができる。
【００３４】
ところで、上記の方法では、切り替えスイッチＳＷ、及びスイッチＳ１、Ｓ２、Ｓ３を切り替えて、各電極２１、２２、２３における第１の電圧Ｖ１、Ｖ２、Ｖ３、及び第２の電圧Ｖ’１、Ｖ’２、Ｖ’３を順次測定している。従って、これらの切り替え時間内に、信号源Ｖsの電位が変化すると、信号源Ｖsの位置検出に誤差が生じるおそれがある。例えば、図５（ａ）に示した心電図において、ＱＲＳ波のように、信号源Ｖsの電位が急激に変化する場合、ＱＲＳ波が発生する信号源Ｖsの位置を精度よく検出することが難しくなる。
【００３５】
本発明は、各抵抗に接続したときに電極に生じる出力電圧の測定時刻に時間差が生じた場合でも、生体内の信号源の位置を精度よく検出することができる生体内信号源位置検出装置、及び生体内信号源位置検出方法を提供するものである。
【００３６】
以下、本発明の実施形態を図面に基づいて詳細に説明する。なお、本発明は、以下の実施形態に限定されるものではない。また、本発明の効果を奏する範囲を逸脱しない範囲で、適宜変更は可能である。なお、以下の説明において、特に断らない限り、「電極」は、生体表面に取り付ける部材をいい、「電位」は、電気的レベルをいい、「電圧」は、測定された電気的レベルをいう。
【００３７】
図６は、本発明の一実施形態における生体内信号源位置検出装置の構成を模式的に示した図である。
【００３８】
図６に示すように、本実施形態における生体内信号源位置検出装置は、生体１０の表面に配置する３つの電極２１、２２、２３と、各電極２１、２２、２３とグランド電極２０との間に、相互に切り替え可能な第１の外部抵抗及び第２の外部抵抗を並列接続する切替手段（接続手段）３０とを有している。ここでは、第１の外部抵抗の抵抗値を無限大とし、第２の外部抵抗の抵抗値をＲgとしている。これにより、各電極２１、２２、２３とグランド電極２０との間は、切り替えスイッチＳＷによって、外部抵抗が接続されていない場合と、外部抵抗Ｒgが接続されている場合とに切り替えられる。ここでは、生体１０の表面に、グランド電極２０を配置し、これをグランド電位にしているが、必ずしも、生体１０の表面に、グランド電極２０を配置する必要はない。
【００３９】
生体１０の表面に配置した各電極２１、２２、２３には、生体１０内の信号源Ｖsからの電圧が発生し、その電圧をアンプ（測定手段）４０で増幅して、出力電圧Ｖoutが出力される。各電極２１、２２、２３とアンプ４０との間には、スイッチＳ１、Ｓ２、Ｓ３がそれぞれ配置され、各スイッチＳ１、Ｓ２、Ｓ３を順次、導通させることにより、各電極２１、２２、２３に生じた電圧が出力電圧Ｖoutとして測定される。
【００４０】
図７は、各電極２１、２２、２３と外部抵抗との接続状態の切り替えステップを示した表である。
【００４１】
図７に示すように、ステップ１〜ステップ６で、切り替えスイッチＳＷ、及びスイッチＳ１、Ｓ２、Ｓ３を、それぞれ切り替えることにより、各電極２１、２２、２３とグランド電位との間に、外部抵抗が接続されていないときに各電極２１、２２、２３に生じる第１の出力電圧Ｖ１、Ｖ２、Ｖ３、及び各電極２１、２２、２３とグランド電位との間に、外部抵抗Ｒgが接続されているときに各電極２１、２２、２３に生じる第２の出力電圧Ｖ’１、Ｖ’２、Ｖ’３が測定される。なお、図７では、各電極２１、２２、２３を、それぞれ、チャネルｃｈ１、ｃｈ２、ｃｈ３と表示している。
【００４２】
図７には、各ステップにおける切り替え時間を示している。切り替え時間は、切り替え手段４０の切り替え特性によって決まるが、ここでは、０．５ｍｓを例示している。この場合、各電極２１、２２、２３に外部抵抗を接続したときに得られる電圧Ｖ１、Ｖ２、Ｖ３、及びＶ’１、Ｖ’２、Ｖ’３の測定時刻は、それぞれ、０．５ｍｓの時間差がある。すなわち、ステップ１〜６で測定データを得る１サイクルの時間は、３ｍｓ（３３３Ｈｚ）となる。
【００４３】
このように、各抵抗に接続したときに各電極２１、２２、２３に生じる出力電圧の測定時刻に時間差が生じた場合でも、信号源Ｖsの位置検出に誤差が生じないようにするためには、各出力電圧を、同一の測定時刻におけるデータに補正する必要がある。
【００４４】
図８及び図９は、各抵抗に接続したときに各電極２１、２２、２３に生じる出力電圧を、同一の測定時刻におけるデータに補正する方法を説明した図である。
【００４５】
図８は、外部抵抗が接続されていないときに各電極２１、２２、２３（チャネルｃｈ１、ｃｈ２、ｃｈ３）に生じた第１の出力電圧Ｖ１、Ｖ２、Ｖ３の測定データを示している。各チャネルｃｈ１、ｃｈ２、ｃｈ３において、代表する３つの黒丸で示した点が、３サイクルの測定で得られた時系列な測定データを示す。なお、チャネルｃｈ１、ｃｈ２、ｃｈ３に生じた第２の出力電圧Ｖ’１、Ｖ’２、Ｖ’３の測定データは省略している。
【００４６】
図８に示すように、１サイクルの測定で、チャネルｃｈ１、ｃｈ２、ｃｈ３で得られた測定データは、Ｖ１（ｔ）、Ｖ２（ｔ＋Δｔ２）、Ｖ３（ｔ＋Δｔ３）と表される。ここで、Δｔ２は、チャネルｃｈ１とｃｈ２との測定時間差、Δｔ３は、チャネルｃｈ１とｃｈ３との測定時間差を表す。
【００４７】
図９は、チャネルｃｈ３において、Ｖ３の時系列な測定データから、Ｖ３の時間変化を近似するフィッティングカーブ（近似波形）Ｃを描いた図である。図９に示すように、このフィッティングカーブＣを用いることによって、切り替え時間（Δｔ３）に応じたチャネルｃｈ１の測定データＶ１（ｔ）と同一の測定時刻ｔにおけるチャネルｃｈ３の測定データＶ３（ｔ＋Δｔ３）の補正値Ｖｃ３（ｔ）（図中の白丸の点）を得ることができる。
【００４８】
同様の方法で、チャネルｃｈ２の測定データＶ２（ｔ＋Δｔ２）についても、Ｖ２の時系列な測定データから求めたフィッティングカーブを用いて、切り替え時間（Δｔ２）に応じた補正値Ｖｃ２（ｔ）を得ることができる。さらに、チャネルｃｈ１、ｃｈ２、ｃｈ３に生じた第２の出力電圧Ｖ’１、Ｖ’２、Ｖ’３の測定データについても、同様の方法で、チャネルｃｈ１の測定データＶ１（ｔ）と同一の測定時刻ｔにおける補正値Ｖｃ’１（ｔ）、Ｖｃ’２（ｔ）、Ｖｃ’３（ｔ）を得ることができる。
【００４９】
このようにして、サイクル毎に、Ｖｉ（ｉ＝１，２，３）及びＶ’ｉ（ｉ＝１，２，３）の時系列な測定データから得られたフィッティングカーブ（近似波形）を用いて、同一の測定時刻における補正値Ｖｃｉ及びＶｃ’ｉを得ることができる。
【００５０】
このようにして得られた補正値Ｖｃｉ及びＶｃ’ｉから、比Ｖｃｉ／Ｖｃ’ｉ（ｉ＝１，２，３）を求めることによって、式（１）〜（３）、及び式（４）〜（５）を用いて、生体内の信号源Ｖsの位置を精度よく検出することができる。
【００５１】
なお、フィッティングカーブ（近似波形）を求める方法は、特に限定されず、例えば、多次元補完だけでなく、線形補間でもよい。
【００５２】
本実施形態における生体内信号源位置検出装置は、図６に示すように、図２に示した構成に対して、近似手段５０、算出手段６０、及び検出手段をさらに備えたものである。
【００５３】
アンプ（測定手段）４０に接続された近似手段５０では、第１の出力電圧Ｖｉ（ｉ＝１，２，３）及び第２の出力電圧Ｖ’ｉ（ｉ＝１，２，３）の測定を１サイクルとして、繰り返し行い、第１の出力電圧Ｖｉ及び第２の出力電圧Ｖ’ｉの時系列な測定データから、第１の出力電圧Ｖｉ及び第２の出力電圧Ｖ’ｉの時間変化を近似するフィッティングカーブ（近似波形）を取得する。
【００５４】
また、近似手段５０に接続された算出手段６０では、各サイクルにおいて、測定時刻に時間差を有する第１の出力電圧Ｖｉ及び第２の出力電圧Ｖ’ｉに対して、近似手段５０で取得したフィッティングカーブを用いて、同一の測定時刻における第１の出力電圧Ｖｉ及び第２の出力電圧Ｖ’ｉの補正値Ｖｃｉ及びＶｃ’ｉを算出する。
【００５５】
また、算出手段６０に接続された検出手段７０では、サイクル毎に算出された補正値Ｖｃｉ及びＶｃ’ｉの比Ｖｃｉ／Ｖｃ’ｉ（ｉ＝１，２，３）を求め、これら３つの比Ｖｃｉ／Ｖｃ’ｉに基づいて、生体内の信号源Ｖsの位置を検出する。
【００５６】
本実施形態によれば、近似手段５０で取得したフィッティングカーブ（近似波形）を用いて、同一の測定時刻における出力電圧Ｖｉ及びＶ’ｉ（ｉ＝１，２，３）の補正値Ｖｃｉ及びＶｃ’ｉ（ｉ＝１，２，３）を算出し、算出した補正値Ｖｃｉ及びＶｃ’ｉに基づいて、生体内の信号源Ｖsの位置を検出することによって、各抵抗に接続したときに各電極２１、２２、２３に生じる出力電圧Ｖｉ及びＶ’ｉの測定時刻に時間差が生じた場合でも、信号源Ｖsの位置を精度よく検出することができる。
【００５７】
本実施形態において、各抵抗に接続したときに各電極２１、２２、２３に生じる出力電圧Ｖｉ及びＶ’ｉの測定時刻の時間差は、切替手段３０で、電極２１、２２、２３及び外部抵抗を切り替えるときの切り替え時間によって生じる。この場合、切り替え時間（例えば０．５ｍｓ）が、測定時間差になるため、図９に示したように、フィッティングカーブ（近似波形）を用いて、測定時間差の分だけ時刻をずらすことによって、出力電圧Ｖｉ及びＶ’ｉの補正値Ｖｃｉ及びＶｃ’ｉを、容易に算出することができる。
【００５８】
ところで、信号源Ｖｓと、各電極２１、２２、２３から構成されるそれぞれの回路は生体のコンデンサ成分の影響を受ける。そのため、信号源Ｖｓで発生した電圧が、各電極２１、２２、２３に伝わる時間に差が生じることもある。
【００５９】
この時間差が無視できない場合には、各電極２１、２２、２３に生じる出力電圧Ｖｉ及びＶ’ｉ（ｉ＝１，２，３）の測定時刻に時間差は、外部抵抗を切り替えるときの切り替え時間の他に、生体のコンデンサ成分の影響による時間差が加わることになる。
【００６０】
この場合、外部抵抗を切り替えるときの切り替え時間は知ることができるが、信号源Ｖsで発生した電圧から生体のコンデンサ成分の影響を知ることは難しい。従って、図９に示したような方法では、出力電圧Ｖｉ及びＶ’ｉ（ｉ＝１，２，３）の補正値Ｖｃｉ及びＶｃ’ｉ（ｉ＝１，２，３）を求めることは難しい。
【００６１】
図１０及び図１１は、各電極２１、２２、２３に生じる出力電圧Ｖｉ及びＶ’ｉ（ｉ＝１，２，３）の測定時刻に時間差を知ることが難しい場合に、出力電圧Ｖｉ及びＶ’ｉの補正値Ｖｃｉ及びＶｃ’ｉ（ｉ＝１，２，３）を求める方法を説明した図である。
【００６２】
図１０は、外部抵抗が接続されていないときに各電極２１、２３（チャネルｃｈ１、ｃｈ３）に生じた第１の電圧Ｖ１、Ｖ３の測定データを示している。各チャネルｃｈ１、ｃｈ３において、代表する２つの黒丸で示した点が、２サイクルの測定で得られた時系列な測定データを示す。
【００６３】
図１０に示すように、１サイクルの測定で、チャネルｃｈ１、ｃｈ３で得られた測定データは、Ｖ１（ｔ）、Ｖ３（ｔ＋Δｔ３）と表される。ここで、Δｔ３は、チャネルｃｈ１とｃｈ３との測定時間差を表す。
【００６４】
図１１は、チャネルｃｈ１、ｃｈ３において、Ｖ１及びＶ３の時系列な測定データから、Ｖ１及びＶ３の時間変化を近似するフィッティングカーブ（近似波形）Ｃ１及びＣ３を描いた図である。ここで、フィッティングカーブＣ３は、図１１に示すように、フィッティングカーブＣ１、Ｃ３の特徴ポイント、例えば、フィッティングカーブＣ１、Ｃ３のピークＫ１、Ｋ３が合うように時間軸をずらしてある。
【００６５】
ここで、同一の測定時刻ｔにおいて、チャネルｃｈ１、ｃｈ３で測定される出力電圧の変化は、信号源Ｖsでの電圧変化と近似すると考えられる。従って、フィッティングカーブＣ３において、Ｖ１（ｔ）と同じ時刻ｔに対する点（図中の白丸の点）が、同一の測定時刻ｔにおけるチャネルｃｈ３の測定データＶ３（ｔ＋Δｔ３）の補正値Ｖｃ３（ｔ）とすることができる。
【００６６】
同様の方法で、チャネルｃｈ２の測定データＶ２（ｔ＋Δｔ２）についても、Ｖ２の時系列な測定データから得られたフィッティングカーブの特徴ポイントを、Ｖ１のフィッティングカーブの特徴ポイントと合わせることによって、測定データＶ１（ｔ）と同一の測定時刻ｔにおける補正値Ｖｃ２（ｔ）を得ることができる。さらに、チャネルｃｈ１、ｃｈ２、ｃｈ３に生じた第２の出力電圧Ｖ’１、Ｖ’２、Ｖ’３の測定データについても、同様の方法で、チャネルｃｈ１の測定データＶ１（ｔ）と同一の測定時刻ｔにおける補正値Ｖｃ’１（ｔ）、Ｖｃ’２（ｔ）、Ｖｃ’３（ｔ）を得ることができる。
【００６７】
このようにして、サイクル毎に、Ｖｉ（ｉ＝１，２，３）及びＶ’ｉ（ｉ＝１，２，３）の時系列な測定データから得られたフィッティングカーブ（近似波形）の特徴ポイントを合わせることによって、同一の測定時刻における補正値Ｖｃｉ及びＶｃ’ｉを得ることができる。
【００６８】
以上、本発明を好適な実施形態により説明してきたが、こうした記述は限定事項ではなく、もちろん、種々の改変が可能である。
【００６９】
例えば、上記実施形態では、生体１０表面に配置する電極が３個で、外部抵抗が２個の場合を例に説明したが、勿論、これに限定されず、電極及び外部抵抗は、任意の個数であってもよい。
【００７０】
すなわち、生体表面に配置する電極の個数をｍ個、各電極とグランド電位との間に並列接続する外部抵抗の個数をｎ個としたとき、本発明における生体内信号源位置検出装置は、各電極に各外部抵抗を並列接続したときに各電極に生じるｍ×ｎ個の出力電圧Ｖｉ（ｉは１〜ｍ×ｎの整数）を測定する測定手段と、ｍ×ｎ個の出力電圧Ｖｉの測定を１サイクルとして、繰り返し行い、各出力電圧Ｖｉの時系列な測定データから、各出力電圧Ｖｉの時間変化を近似する近似波形を取得する近似手段と、各サイクルにおいて、測定時刻に時間差を有する各出力電圧Ｖｉに対して、近似手段で取得した近似波形を用いて、同一の測定時刻における各出力電圧Ｖｉの補正値Ｖｃｉ（ｉは１〜ｍ×ｎの整数）を算出する算出手段と、サイクル毎に、算出された各補正値Ｖｃｉに基づいて、生体内の信号源の位置を検出する検出手段とを備える。
【００７１】
また、上記実施形態では、各電極２１、２２、２３とグランド電極２０との間に、第１の外部抵抗及び第２の外部抵抗を、切替手段３０により、相互に切り替えて並列接続したが、必ずしも、切替手段３０を用いなくてもよい。
【００７２】
すなわち、各電極２１、２２、２３は、隣接した電極２１ａ、２１ｂと、電極２２ａ、２２ｂと、電極２３ａ、２３ｂとで構成されており、一方側の電極２１ａ、２２ａ、２３ａに、第１の外部抵抗Ｒg1を接続することによって、１つの回路を形成し、他方側の電極２１ｂ、２２ｂ、２３ｂに、第２の外部抵抗Ｒg2を接続することによって、別の回路を形成する。
【００７３】
これにより、各電極２１，２２，２３に生じる第１及び第２の電圧を測定し、生体内の信号源Ｖsの位置を検出することができる。
【００７４】
図１２は、切り替え手段３０を用いない場合の生体内信号源位置検出装置の構成を例示した図である。なお、図中において、信号源Ｖs、並びに信号源Ｖsと各電極２１（２１ａ、２１ｂ）、２２（２２ａ、２２ｂ）、２３（２３ａ、２３ｂ）及びグランド電極２０との間の内部抵抗Ｒｂ1、Ｒｂ２、Ｒｂ３、Ｒｂ０は、省略している。
【００７５】
図１２に示すように、生体１０の表面に配置された３つの電極２１（２１ａ、２１ｂ）、２２（２２ａ、２２ｂ）、２３（２３ａ、２３ｂ）と、グランド電極２０との間に、第１の外部抵抗Ｒg1及び第２の外部抵抗Ｒg2が、それぞれ並列接続されている。そして、電極２１ａとグランド電極２０との間に、第１の外部抵抗Ｒg1を並列接続したときに生じる第１の出力電圧Ｖ１、及び電極２１ｂとグランド電極２０との間に、第２の外部抵抗Ｒg2を並列接続したときに生じる第２の出力電圧Ｖ’１が、それぞれアンプ４０Ａ１、４０Ｂ１で増幅されて測定される。同様に、電極２２ａとグランド電極２０との間に、第１の外部抵抗Ｒg1を並列接続したときに生じる第１の出力電圧Ｖ２、及び電極２２ｂとグランド電極２０との間に、第２の外部抵抗Ｒg2を並列接続したときに生じる第２の出力電圧Ｖ’２が、それぞれアンプ４０Ａ２、４０Ｂ２で増幅されて測定され、電極２３ａとグランド電極２０との間に、第１の外部抵抗Ｒg1を並列接続したときに生じる第１の出力電圧Ｖ３、及び電極２３ｂとグランド電極２０との間に、第２の外部抵抗Ｒg2を並列接続したときに生じる第２の出力電圧Ｖ’３が、それぞれアンプ４０Ａ３、４０Ｂ３で増幅されて測定される。なお、この場合、各電極２１（２１ａ、２１ｂ）、２２（２２ａ、２２ｂ）、２３（２３ａ、２３ｂ）とグランド電極２０との間に、第１の外部抵抗Ｒg1及び第２の外部抵抗Ｒg2を、それぞれ並列接続する接続手段としては、配線等により行うことができる。
【００７６】
なお、この場合、各抵抗に接続したときに各電極２１（２１ａ、２１ｂ）、２２（２２ａ、２２ｂ）、２３（２３ａ、２３ｂ）に生じる出力電圧Ｖｉ及びＶ’ｉの測定時刻の時間差は、外部抵抗を切り替えるときの切り替え時間によるものはなく、信号源Ｖsから各電極２１（２１ａ、２１ｂ）、２２（２２ａ、２２ｂ）、２３（２３ａ、２３ｂ）までの距離の違いによる時間差となる。従って、この場合は、図１０、１１に示した方法で、出力電圧Ｖｉ及びＶ’ｉの補正値Ｖｃｉ及びＶｃ’ｉ（ｉ＝１，２，３）を求めることができる。
【符号の説明】
【００７７】
１０      生体 
２０      接地電極 
２１、２２、２３      電極 
３０      切替手段 
４０      アンプ（測定手段） 
５０      近似手段 
６０      算出手段 
７０      検出手段 
【図１】
0006518084.tif 000008
【図２】
0006518084.tif 000009
【図３】
0006518084.tif 000010
【図４】
0006518084.tif 000011
【図５】
0006518084.tif 000012
【図６】
0006518084.tif 000013
【図７】
0006518084.tif 000014
【図８】
0006518084.tif 000015
【図９】
0006518084.tif 000016
【図１０】
0006518084.tif 000017
【図１１】
0006518084.tif 000018
【図１２】
0006518084.tif 000019