(57)【要約】  【課題】  簡単なハードウエア構成で、信頼性のある経穴の判定を行う。  【解決手段】  電流発生部１から発生させた電流を、被検体の皮膚３０の測定点＃１〜＃ｎへ電流印加電極３1、３2〜３nを通じて印加し、通電した電流と、通電によって皮膚に生じた電圧を測定部６にて測定し、周波数解析部７にて測定された電流と電圧の周波数解析を行い、各測定点での皮膚インピーダンスZ(jf)を算出し、特徴量生成部８にて測定点＃１〜＃ｎでの皮膚インピーダンスZ(jf)の周波数特性から、皮膚インピーダンスZ(jf)の実数部RZ(f)とXZ(f)虚数部の比の周波数特性K(ｆ)＝XZ(f)／RZ(f)を算出し、判定部９にてその周波数特性K(f)の違いに基づき各測定点のなかから経穴位置を判定する。【選択図】  図１
2006296809.tif 000002
【特許請求の範囲】
【請求項１】
少なくとも１つの周波数成分から構成される電流を発生する電流発生部と、
前記電流発生部からの出力電流を生体の皮膚に配置した複数の電極により複数の異なる皮膚上の測定点にほぼ同時に通電する電極系と、
前記複数の測定点にそれぞれ通電された電流を検出する電流検出器と、
前記電流検出器により検出された電流と、前記電極系の通電によって前記複数の測定点の皮膚に生じた電圧を測定する測定部と、
前記測定部により測定された電流と電圧の両方の周波数解析を行い、各測定点での皮膚インピーダンスの周波数特性を得る周波数解析部と、
前記周波数解析部により解析された前記複数の測定点での皮膚インピーダンスの周波数特性から、前記各測定点での皮膚インピーダンスの特徴量を生成する特徴量生成部と、
前記特徴量生成部により生成された前記各測定点での特徴量から経穴位置を判定する判定部と、
前記特徴量生成部により生成された前記各測定点での特徴量と前記判定部により判定された経穴位置とを表示する表示手段とを、
備えた経穴位置評価装置。
【請求項２】
前記周波数解析部は、前記各測定点での皮膚インピーダンスの周波数特性において Cole-Cole の円弧則が成立するか否かを確認することを特徴とする請求項１に記載の経穴位置評価装置。
【請求項３】
前記電流発生部において発生させる通電電流は、１つの最低周波数成分と、その整数倍の周波数成分を含むことを特徴とする請求項１又は２に記載の経穴位置評価装置。
【請求項４】
前記電流発生部は、通電電流に含ませる周波数成分の大きさを、周波数ごと並びに測定位置ごとに変更可能とすることを特徴とする請求項１から３のいずれか１つに記載の経穴位置評価装置。
【請求項５】
前記特徴量生成部において生成される特徴量は、少なくとも１つの周波数における皮膚インピーダンスの実数部と虚数部との比であることを特徴とする請求項１から４のいずれか１つに記載の経穴位置評価装置。
【請求項６】
前記特徴量生成部は、少なくとも１つの周波数における皮膚インピーダンスの実数部と虚数部との比のなかで前記各測定点で差が顕著となる周波数での比を特徴量とすることを特徴とする請求項１から４のいずれか１つに記載の経穴位置評価装置。
【請求項７】
前記特徴量生成部において生成される特徴量は、通電電流に含まれる最低周波数での皮膚インピーダンスの実数部と虚数部との比であることを特徴とする請求項１から４のいずれか１つに記載の経穴位置評価装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、鍼治療などにおいて、鍼刺激の効果と経穴（ツボ）と呼ばれる刺激点の位置を非侵襲的かつ客観的に、しかも高精度に判定できる経穴位置評価装置に関する。
【背景技術】
【０００２】
古くから皮膚表面上に周辺より皮膚抵抗が低い部位が小さな領域で存在し、その分布は鍼灸治療における刺激点である経穴の分布と比較的よく一致することが知られている。このような皮膚抵抗が周辺より低い部位（皮膚抵抗減弱点と言う）は、良導絡自律神経調整法や低周波治療器などにおいて治療に利用されており、これらの皮膚抵抗減弱点を探索する装置も市販され、実際に臨床での診断や治療に用いられている。
【０００３】
従来、これらの皮膚抵抗減弱点である経穴を探索（判定）する装置としては、経穴周辺の皮膚の電気的特性として皮膚の直流抵抗を測定し、その違いから経穴を判定する方法が多用されてきたが、これら直流抵抗の測定による経穴の判定では、測定結果が測定電極の接触状態に依存したり、皮膚電圧の自発的な変動との区別が明確にできないなどの問題があり、測定結果の信頼性・再現性に乏しい。
【０００４】
そこで、近年においては、これら測定結果の信頼性・再現性を向上させる手法として、経穴周辺の皮膚の電気的特性として、経穴周辺の皮膚への交流電流の印加により生成する電圧に基づいて皮膚インピーダンスの周波数特性の複素平面上での分布（皮膚インピーダンス軌跡）を測定し、皮膚インピーダンス軌跡の違いに基づいて経穴の判定を行うものがある。例えば皮膚インピーダンス軌跡の形状を特徴付ける４つのパラメータ（後述するZ0、Z∞、β、τm）のうち、測定電極の接触状態に大きく依存せず、経穴とそうでない部位での電気的特性の差を明確に反映するパラメータとして Cole-Cole の円弧則における中心緩和時間τmに着目し、そのパラメータτmの違いに基づいて経穴の判定を行っている。このようにすることで、従来技術が有する問題点を回避できるようになっている（例えば下記の特許文献１参照）。
【特許文献１】特開２００４−３３７３４９号公報（要約書、請求項２）
【発明の開示】
【発明が解決しようとする課題】
【０００５】
しかしながら、上記従来技術では、経穴とそうでない部位での電気的特性の差を評価する際に、皮膚インピーダンス軌跡の形状を特徴付けるパラメータとして、Cole-Cole の円弧則における中心緩和時間τmを算出しているが、以下のような問題点を有していた。
皮膚インピーダンス軌跡の形状を特徴付ける中心緩和時間τmを算出する際、皮膚インピーダンスの周波数特性において Cole-Cole の円弧則が成立している必要があるが、その成立の可否を確認する手段を有していない。また、測定された複数の周波数における皮膚インピーダンスから中心緩和時間τmを非線形最小２乗法などによって推定する必要があり、推定するための手順が煩雑で、これを実現するハードウエアの構成が複雑になり、装置の小型化には不向きであった。さらに、中心緩和時間τmを精度良く推定するには、図４に示すように、皮膚インピーダンス軌跡上に測定された皮膚インピーダンスデータが偏らず、まんべんなく分布している必要がある。
【０００６】
しかし、実際の皮膚で測定される皮膚インピーダンスデータの分布は、図５に示すように、インピーダンス軌跡上の一部分に偏って存在することが多く、又、その分布の仕方は測定点ごとに大きく異なっている場合が多い。このような場合、皮膚インピーダンス軌跡の形状を特徴付ける中心緩和時間τmを精度良く推定するには、極めて低い周波数も含めて広範囲の周波数範囲の皮膚インピーダンスを測定する必要があった。そのため、時間分解能が低下するほか、測定対象周波数が低くなればなるほど温度ドリフトや体動などによる外乱の影響を受けやすくなり、測定結果の信頼性・再現性を保証するのが困難であった。
【０００７】
本発明は上記従来技術における問題点を解決するもので、簡単なハードウエア構成で、信頼性・再現性のある経穴の判定を行うことができる経穴位置評価装置を提供することを目的とする。
【課題を解決するための手段】
【０００８】
上記従来の問題点を解決するために、本発明に係る経穴位置評価装置は、
少なくとも１つの周波数成分から構成される電流を発生する電流発生部と、
前記電流発生部からの出力電流を生体の皮膚に配置した複数の電極により複数の異なる皮膚上の測定点にほぼ同時に通電する電極系と、
前記複数の測定点にそれぞれ通電された電流を検出する電流検出器と、
前記電流検出器により検出された電流と、前記電極系の通電によって前記複数の測定点の皮膚に生じた電圧を測定する測定部と、
前記測定部により測定された電流と電圧の両方の周波数解析を行い、各測定点での皮膚インピーダンスの周波数特性を得る周波数解析部と、
前記周波数解析部により解析された前記複数の測定点での皮膚インピーダンスの周波数特性から、前記各測定点での皮膚インピーダンスの特徴量を生成する特徴量生成部と、
前記特徴量生成部により生成された前記各測定点での特徴量から経穴位置を判定する判定部と、
前記特徴量生成部により生成された前記各測定点での特徴量と前記判定部により判定された経穴位置とを表示する表示手段とを備えた構成とした。
この構成により、評価に用いる特徴量の算出手順が簡単になり、測定結果の信頼性・再現性を従来と比べ高くすることができるほか、評価に用いる特徴量は従来着目していた中心緩和時間τmと相関があるので、測定電極の接触状態に大きく依存せず、経穴とそうでない部位での電気的特性の差を明確に反映する評価パラメータの違いに基づいて経穴の判定を行うことが可能になる。
【０００９】
また、本発明に係る経穴位置評価装置は、
前記周波数解析部が、前記各測定点での皮膚インピーダンスの周波数特性において Cole-Cole の円弧則が成立するか否かを確認することを特徴とする。
この構成により、測定された皮膚インピーダンスにおいて Cole-Cole の円弧則の成立可否を確認でき、測定結果の信頼性・再現性を従来と比べ高くすることができる。
【００１０】
また、本発明に係る経穴位置評価装置は、
前記電流発生部において発生させる通電電流には、１つの最低周波数成分と、その整数倍の周波数成分が含まれることを特徴とする。
この構成により、離散フーリエ変換をした場合の隣接する周波数へのスペクトルの洩れは無く、時間分解能の下限は最低周波数のみに依存して決まるので、高時間分解能で皮膚インピーダンス軌跡を測定することができる。
【００１１】
また、本発明に係る経穴位置評価装置は、
前記電流発生部が、通電電流に含ませる周波数成分の大きさを、周波数ごと並びに測定位置ごとに変更可能とすることを特徴とする。
この構成により、各周波数、各測定点ごとに最適な信号対雑音比を確保でき、高い精度で皮膚インピーダンス軌跡を測定することができる。
【００１２】
また、本発明に係る経穴位置評価装置は、
前記特徴量生成部において生成される特徴量が、少なくとも１つの周波数における皮膚インピーダンスの実数部と虚数部との比であることを特徴とする。
この構成により、パラメータの推定手順の煩雑さを回避することができる。
【００１３】
また、本発明に係る経穴位置評価装置は、
前記特徴量生成部において、少なくとも１つの周波数における皮膚インピーダンスの実数部と虚数部との比のなかで、前記各測定点で差が顕著となる周波数での比を特徴量とする。
この構成により、パラメータの推定手順の煩雑さを回避することができるほか、判定部での経穴位置の判定の信頼性がより高くなる。
【００１４】
また、本発明に係る経穴位置評価装置は、
前記特徴量生成部において生成される特徴量が、通電電流に含まれる最低周波数での皮膚インピーダンスの実数部と虚数部との比であることを特徴とする。
この構成により、パラメータの推定手順の煩雑さを回避することができるほか、経穴位置の判定を迅速に行える。
【発明の効果】
【００１５】
本発明によれば、測定された皮膚インピーダンスにおいて Cole-Cole の円弧則が成立するか否かを確認でき、評価に用いる特徴量の算出手順を簡単なハードウエア構成で実現できる。また、測定結果の信頼性・再現性を従来と比べ高くすることができるほか、評価に用いる特徴量は従来着目していた中心緩和時間τmと相関があるので、測定電極の接触状態に大きく依存せず、経穴とそうでない部位での電気的特性の差を明確に反映する評価パラメータの違いに基づいて経穴の判定を行うことが可能になる。したがって、簡単なハードウエア構成でより信頼性・再現性の高い経穴位置装置を提供することができる。
【発明を実施するための最良の形態】
【００１６】
以下、図面を参照して本発明の実施の形態について説明する。図１は本発明の経穴位置評価装置の第１の実施の形態の概略構成を示すブロック図である。電流発生部１から発生された交流電流Ｉ1〜Ｉnはそれぞれ電流検出器２1〜２n、電流印加電極３1〜３nを介して、被検体の皮膚３０の各測定点＃１〜＃ｎへ印加され、接地電極４へと流れていく。この通電によって生じた電流印加電極３1〜３nとマイナス電極５（−）との間の皮膚に生じた電圧Ｖ1〜Ｖnの降下を接地電極４の電位を基準にして測定部６により測定する。このような電極系により測定する手法は３電極法と呼ばれ、電流印加電極３1〜３nの直下、つまり、測定点＃１〜＃ｎの直下における皮膚インピーダンスを測定する際によく用いられる。
【００１７】
図２に電流発生部１並びに測定部６のより詳細な構成例を示すブロック図を示す。電流発生部１は少なくとも１つの電流源１1〜１nから構成されており、電流源１1〜１nはそれぞれ各測定点＃１〜＃ｎへ通電する交流電流Ｉ1〜Ｉnを発生する。ここで、各電流源１1〜１nから発生される電流Ｉ1〜Ｉnに含ませる周波数成分の数、周波数、大きさは、制御部２０からの制御信号２１０により設定可能になっている。測定部６は、通電によって各測定点＃１〜＃ｎ直下の皮膚に生じた電圧Ｖ1〜Ｖnを測定するための少なくとも１つの差動増幅器６１1〜６１nと、少なくとも１つのハイパスフィルタ（ＨＰＦ）６２1〜６２n（並びに６７1〜６７n）と、少なくとも１つのプログラマブルゲインアンプ（ＰＧＡ）６３1〜６３n（並びに６８1〜６８n）と、少なくとも１つのローパスフィルタ（ＬＰＦ）６４1〜６４n（並びに６９1〜６９n）と、少なくとも１つのＡ／Ｄ変換器６５1〜６５n（並びに７０1〜７０n）とから構成されている。
【００１８】
通電によって測定点＃１〜＃ｎの皮膚に生じた電圧Ｖ1〜Ｖnはそれぞれ差動増幅器６１1〜６１nにより測定されるが、測定された電圧Ｖ1〜Ｖnには通電に起因しない不要な信号成分が含まれている。これは周波数解析部７での正確な周波数解析の妨げになるため、ハイパスフィルタ６２1〜６２nによりそれぞれ差動増幅器６１1〜６１nの出力から不要な低周波成分を除去する。不要な低周波成分が除去された測定点＃１〜＃ｎの電圧Ｖ1〜Ｖnはそれぞれ、必要に応じてプログラマブルゲインアンプ６３1〜６３nにより増幅され、次いでローパスフィルタ６４1〜６４nにより不要な高周波成分が除去される。
【００１９】
また、測定点＃１〜＃ｎの皮膚に通電された電流Ｉ1〜Ｉnは、それぞれ電流検出器２1〜２nにより測定されるが、測定点＃１〜＃ｎでの電圧Ｖ1〜Ｖnと同じ信号処理を、測定点＃１〜＃ｎに通電した電流Ｉ1〜Ｉnに対して行なわなければ、図１に示す周波数解析部７での周波数解析の手順が繁雑になり、解析結果も不正確になってしまう。そのため、測定点＃１〜＃ｎの電圧Ｖ1〜Ｖnと同じ信号処理を行なうためにハイパスフィルタ６７1〜６７n、プログラマブルゲインアンプ６８1〜６８n、ローパスフィルタ６９1〜６９nが必要となる。ここで、プログラマブルゲインアンプ６３1〜６３n、６８1〜６８nの増幅率は制御部２０からの制御信号２１１、２１２によって制御可能にしておく。
【００２０】
なお、本発明は、測定された電流と電圧に対して行う信号処理の順序や手段を限定するものではなく、周波数解析部７において正確な周波数解析が行える処理であればどのような順序や手段を用いてもよい。このようにして測定された、測定点＃１〜＃ｎにそれぞれ印加した電流Ｉ1〜Ｉn、各測定点での電圧Ｖ1〜Ｖnは、Ａ／Ｄ変換器７０1〜７０n、６５1〜６５nによってディジタル信号へと変換され、周波数解析部７へと送出される。
【００２１】
図１に示す周波数解析部７では、測定点＃１〜＃ｎの皮膚にそれぞれ通電された電流Ｉ1〜Ｉnと、測定点＃１〜＃ｎの皮膚の電圧Ｖ1〜Ｖnとに対し周波数解析を行い、その結果を基に各測定点＃１〜＃ｎ直下の皮膚インピーダンスの周波数特性を算出する。以下ではその原理について述べる。今、皮膚インピーダンスが通電電流密度に依存せず、皮膚が線形システムとみなせると仮定する。被検体の皮膚に印加された電流が次式（１）のような電流i(t)である場合、皮膚が線形システムと仮定できる場合、通電によって皮膚に生じた電圧は同様に次式（２）のような電圧v(t)になる。
【００２２】
【数１】
2006296809.tif 000003
【００２３】
ここで、ω=２πf（fは測定周波数）、|In|、|Vn|はそれぞれ各周波数成分の電流と電圧の実効値であり、θin、θvnはそれぞれ各周波数成分fnの電流及び電圧の位相である。式（１）、（２）におけるi(t)、v(t)のフーリエ変換をI(jω)、V(jω)とすると、
【００２４】
【数２】
2006296809.tif 000004
【００２５】
であるから、皮膚インピーダンスZ(jω)は
【００２６】
【数３】
2006296809.tif 000005
【００２７】
と表すことができる。ただし、
【００２８】
【数４】
2006296809.tif 000006
【００２９】
とする。
【００３０】
つまり、必要な周波数fnをi(t)に複数含ませておけば、i(t)、ｖ(t)をフーリエ変換し、それぞれの実数部と虚数部から各周波数におけるインピーダンスが求まる。また、電流波形に含ませる周波数を最低周波数fminの整数倍にすることで、離散フーリエ変換をした場合の隣接する周波数へのスペクトルの洩れは無く、時間分解能の下限は最低周波数のみに依存して決まる。よって、高時間分解能で皮膚インピーダンス軌跡が測定できる。
【００３１】
ところで、皮膚インピーダンスZ(jω)=Rz(ω)-jXz（ω）は次式(9)で表される Cole-Cole の円弧則を満たし、周波数特性の複素平面上での軌跡は、図３に示すような円弧上に
分布することが知られている。
【００３２】
【数５】
2006296809.tif 000007
【００３３】
ただし、ω=2πf（fは測定周波数）、τmは中心緩和時間、βは円弧の中心角をあらわすパラメータであり、
【００３４】
【数６】
2006296809.tif 000008
【００３５】
とする。一般に皮膚インピーダンスの場合、高周波領域のインピーダンスは非常に小さいことからZ∞≒0とみなせる。式(9)を角周波数ωではなく周波数fで表現すると、
【００３６】
【数７】
2006296809.tif 000009
【００３７】
となり、fm = 1/2πτmを特性周波数と呼ぶ。式(9)と(12)は等価であるが、以下の説明では式(12)を用いる。
【００３８】
ここで、測定された皮膚インピーダンスの周波数特性の複素平面上での軌跡が円弧上に分布していても、Cole-Cole の円弧則を満たさない場合がある。その場合も含めた一般的な皮膚インピーダンスを表現すると、上式(12)にもう１つのパラメータαを加えた次式(13)(14)にて表現される。
【００３９】
【数８】
2006296809.tif 000010
【００４０】
ここで、皮膚インピーダンスが測定電流密度に依存し、皮膚が線形システムとみなせない場合が Cole-Cole の円弧則を満たさない例として挙げられる。従来例においては、測定電極の接触状態に大きく依存せず、経穴とそうでない部位での電気的特性の差を明確に反映するパラメータに着目している。具体的には、測定された皮膚インピーダンスが Cole-Cole の円弧則を満足すると仮定し、式(12)中の４つのパラメータ（Z0、Z∞、β、τm=2π/fm）のうち、Cole-Cole の円弧則における中心緩和時間τmに着目し、その違いに基づいて経穴の判定を行っている。しかし、前述したように以下の問題点があった。
・測定された皮膚インピーダンスが Cole-Cole の円弧則の成立を確認する手段を設けていない。
・中心緩和時間τmの推定手順が煩雑である。
・測定データが皮膚インピーダンス軌跡上に一様に分布しない。
【００４１】
そこで、周波数解析部７で算出された各測定点＃１〜＃ｎでの皮膚インピーダンスが測定電流密度に依存しないことをあらかじめ確認してから、判定部９での経穴位置の判定を行なうようにする。こうすることで、測定された皮膚インピーダンスにおいて Cole-Cole の円弧則が成立するか否かを間接的に確認できるので、特徴量生成部８で生成される特徴量の信頼性・再現性を従来例より高くすることができる。
【００４２】
また、特徴量生成部８においては、周波数解析部７から送出された皮膚インピーダンスZ(jf)の周波数特性から、皮膚インピーダンスZ(jf)の実数部RZ(f)と虚数部XZ(f)の比の周波数特性K(ｆ)を算出し、制御部２０より設定される固定周波数f0でのK(f0)を各測定点ごとに生成する。このK(f0)を、従来例におけるパラメータτmの代わりに評価パラメータとして採用する。このような構成にすることで、パラメータの推定手順の煩雑さが回避できるほか、前記K(ｆ)と従来例において着目しているパラメータτmの大小と正の相関があるので、従来例と同様、測定電極３、４、５の接触状態に大きく依存せず、経穴とそうでない部位での電気的特性の差を明確に反映する評価パラメータの違いに基づいて経穴の判定を行える。以下では、その理由について述べる。まず、式(14)を変形すると、次式(15)(16)が得られる。
【００４３】
【数９】
2006296809.tif 000011
【００４４】
ここで、皮膚インピーダンスが Cole-Cole の円弧則を満たし、高周波領域での皮膚インピーダンスが０とみなせる場合、つまり、α=βとZ∞=0が成り立つ場合、上式(15)(16)は次式(17)(18)のように表される。
【００４５】
【数１０】
2006296809.tif 000012
【００４６】
さらに、式(17)(18)の実数部＝RZ(f) と虚数部＝XZ(f) の比K(f)＝XZ(f)／RZ(f) は次式(19)のようになる。
【００４７】
【数１１】
2006296809.tif 000013
【００４８】
よって、βが一定であれば、同一測定周波数fでのK(f)の大小はτmの大小と正の相関を持ち、fmの大小と負の相関を持つ。
【００４９】
ここで、上記説明が成立するのは、βが一定で、かつ、Z∞=0でなければならない。一般に皮膚インピーダンスの場合、βは0.65〜0.9と極めて限定された範囲の値をとることが知られており、さらに、測定電極３、４、５の面積が同一で、測定点＃１〜＃ｎが比較的隣接していれば、ほぼ一定とみなすことができる。また、高周波領域における皮膚インピーダンスの大きさは低周波領域と比較すると小さく、特に、上記３電極法によって測定された皮膚インピーダンスの場合、極めて小さいため、Z∞=0とみなしてよい。
【００５０】
よって、特徴量生成部８において、例えば電流発生部１から発生された電流Ｉには１つの周波数f0の周波数成分を含ませておき、皮膚インピーダンスの実数部と虚数部の比K(f0)を各測定点＃１〜＃ｎごとに算出し、判定部９において、K(f0)が他の測定点よりも低い値になった測定点が経穴であると判定することで、従来例と比べ簡単なハードウエア構成で、かつ、測定電極３、４、５の接触状態に大きく依存せず、経穴とそうでない部位での電気的特性の差を明確に反映する評価パラメータの違いに基づいて経穴の判定を行えることになる。また、上記K(f0)はある固定周波数f0のみの値であり、従来例のように直接、パラメータτmを精度良く求める際のように、極めて低い周波数も含む広範囲の皮膚インピーダンスを測定する必要がない。
【００５１】
なお、前述の判定手法においては、電流Ｉに含ませる周波数はf0１つのみとして説明したが、本発明はそれを限定するものではなく、K(ｆ0)が経穴とそうでない部位での電気的特性の差を明確に反映する周波数であればその数を限定しない。例えばあらかじめ電流発生部１において複数の周波数成分を含む電流Ｉを発生させておいて、その中から各測定点＃１〜＃ｎでのK(f)の値の差が顕著に見られる周波数を特徴量生成部８にて選択し、その周波数でのK(f)を判定部９に送出してもよい。こうすることで、簡単なハードウエア構成で、判定部９での経穴位置の判定の信頼性をより高くすることができるほか、前述したβやZ∞に関する仮定の成立可否を確認できる。
【００５２】
なお、式(19)中のｆm=τm/2πは虚数部XZの絶対値が最大値をとる周波数である。よって特徴量生成部８において生成する特徴量を、周波数解析部７において算出される皮膚インピーダンスの周波数特性の虚数部XZの絶対値が最大値をとる周波数ｆmとしてもよい。その場合、判定部９では、ｆmが他の測定点よりも大きくなった点を経穴と判定するようにすればよい。こうすることで、従来例と比べ簡単なハードウエア構成で、かつ、測定電極３、４、５の接触状態に大きく依存せず、経穴とそうでない部位での電気的特性の差を明確に反映する評価パラメータの違いに基づいて経穴の判定を行えることになる。以上のようにして特徴量生成部８で生成された各測定点＃１〜＃ｎでの特徴量や、判定部９での判定結果は表示部１０へと送出され、モニタなどの表示手段により適宜表示される。
【産業上の利用可能性】
【００５３】
以上のように、本発明にかかる経穴位置評価装置は、測定された皮膚インピーダンスにおいて Cole-Cole の円弧則の成立可否を確認でき、評価に用いる特徴量の算出手順が簡単である。また測定結果の信頼性・再現性を従来と比べ高くすることができるほか、評価に用いる特徴量は、従来着目していたパラメータとしてCole-Cole の円弧則における中心緩和時間τmと相関があるので、測定電極の接触状態に大きく依存せず、経穴とそうでない部位での電気的特性の差を明確に反映する評価パラメータの違いに基づいて経穴の判定を行うことが可能になる。よって、従来例に比べ、より簡単なハードウエア構成で、より信頼性・再現性の高い経穴位置判定を行えることから、鍼治療などにおいて、鍼刺激の効果と経穴（ツボ）と呼ばれる刺激点の位置探索や良導点探索など、医療分野にて皮膚の電気的特性の差を非侵襲的かつ客観的に評価する際に有用である。
【図面の簡単な説明】
【００５４】
【図１】本発明の経穴位置評価装置の第１の実施の形態における概略構成を示すブロック図
【図２】図１の電極系と測定部のより詳細な構成例を示すブロック図
【図３】Cole-Cole の円弧則を説明する説明図
【図４】測定データがインピーダンス軌跡上にまんべんなく分布することを説明する図
【図５】測定データがインピーダンス軌跡上の一部分に偏って分布することを説明する図
【符号の説明】
【００５５】
１  電流発生部
１1、１2〜１n  電流源
２1、２2〜２n  電流検出器
３1、３2〜３n  電流印加電極
４  接地電極
５  マイナス電極
６  測定部
７  周波数解析部
８  特徴量生成部
９  判定部
１０  表示部
２０  制御部
６１1〜６１n  差動増幅器
６２1〜６２n、６７1〜６７n  ハイパスフィルタ（ＨＰＦ）
６３1〜６３n、６８1〜６８n  プログラマブルゲインアンプ（ＰＧＡ）
６４1〜６４n、６９1〜６９n  ローパスフィルタ（ＬＰＦ）
６５1〜６５n、７０1〜７０n  Ａ／Ｄ変換器
２１１、２１２  制御信号
【図１】
2006296809.tif 000014
【図２】
2006296809.tif 000015
【図３】
2006296809.tif 000016
【図４】
2006296809.tif 000017
【図５】
2006296809.tif 000018