(57)【要約】【課題】生体信号と指関節角度の関係の時間変化に対応可能であって、複数の指関節角度を連続的にリアルタイムに推定することができる指関節角度推定装置を提供する。【解決手段】手指関節角度推定装置１０は、各指関節が動いているときの複数の表面筋電信号を入力して、積分処理部１４により各表面筋電信号の特徴量を抽出し、正規化処理部１５により正規化した該特徴量を一次元に並べて特徴量ベクトルとする。ニューラルネットワーク１６は、過去から現在までの時系列の特徴量ベクトルを一次元に並べたフィードフォワード成分と、演算周期毎に現在まで推定していた時系列の出力ベクトルからなるフィードバック成分からなる入力ベクトルを入力して、指関節角度を推定する。【選択図】図１
2010125287.tif 000002
【特許請求の範囲】
【請求項１】
各指関節が動いているときの単数又は複数種類の生体信号を入力して、各生体信号の特徴量を抽出し、抽出した該特徴量を一次元に並べて特徴量ベクトルとする特徴量抽出手段と、
演算周期毎に入力ベクトルに基づいて、各指関節の角度（以下、指関節角度という）の推定を行い、推定した前記各指関節角度を出力ベクトルとして出力するニューラルネットワーク手段とを含み、
前記入力ベクトルは、前記特徴量抽出手段が抽出した過去から現在までの時系列の特徴量ベクトルを一次元に並べたフィードフォワード成分と、出力のフィードバック成分で構成され、
前記フィードバック成分は、前記ニューラルネットワーク手段が、前記演算周期毎に現在まで推定していた時系列の前記出力ベクトルで構成されていることを特徴とする指関節角度推定装置。
【請求項２】
前記生体信号は、表面筋電信号であることを特徴とする請求項１に記載の指関節角度推定装置。
【請求項３】
前記特徴量抽出手段は、前記表面筋電信号に対して積分区間を演算周期毎にずらして積分し、積分して得られた積分値を特徴量として抽出するものである請求項２に記載の指関節角度推定装置。
【請求項４】
前記特徴量抽出手段は、前記表面筋電信号のパワースペクトルの所定のバンド幅毎の区間和を特徴量として抽出するものである請求項２に記載の指関節角度推定装置。
【請求項５】
前記特徴量抽出手段は、前記表面筋電信号のパワースペクトルの平均周波数又は中央周波数を特徴量として抽出するものである請求項２に記載の指関節角度推定装置。
【請求項６】
前記生体信号は、脳波であり、
前記特徴量抽出手段は、前記脳波のα波及びβ波の少なくともいずれか１つの波の増減変化による事象関連同期化、及び事象関連脱同期化を特徴量として抽出するものである請求項１に記載の指関節角度推定装置。
【請求項７】
前記生体信号は、脳部位のＭＲＩ信号であり、
前記特徴量抽出手段は、運動を司る部位における血流中の還元ヘモグロビンの変化から活性化している場所を特徴量として抽出するものである請求項１に記載の指関節角度推定装置。
【請求項８】
前記生体信号は、頭部又は指関節に関係する筋肉を透過させた近赤外光を検出した信号であり、
前記特徴量抽出手段は、頭部又は指関節に関係する筋肉組織を流れている血液中のヘモグロビン酸素化状態の変化を特徴量として抽出するものである請求項１に記載の指関節角度推定装置。
【請求項９】
前記生体信号は、脳磁図であり、
前記特徴量抽出手段は、前記脳磁図の磁場の変化を特徴量として抽出するものである請求項１に記載の指関節角度推定装置。
【請求項１０】
前記生体信号は、指関節の動きによって腱が移動するのを捉えた信号であり、
前記特徴量抽出手段は、前記腱の移動の変化を特徴量として抽出するものである請求項１に記載の指関節角度推定装置。
【請求項１１】
前記生体信号は、指関節の動きによって筋肉の膨らみが変化するのを捉えた信号であり、
前記特徴量抽出手段は、前記筋肉の膨らみの変化を特徴量として抽出するものである請求項１に記載の指関節角度推定装置。
【請求項１２】
前記生体信号は、指関節の動きによって脳から、当該指の筋肉細胞までに伝わる神経信号を検出する信号であり、
前記特徴量抽出手段は、前記指関節の動きに応じて変化する神経信号の変化を特徴量として抽出するものである請求項１に記載の指関節角度推定装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、指関節角度推定装置に関する。
【背景技術】
【０００２】
従来、ヒューマン・ロボットシステムにおいて、表面筋電信号を利用した手指形状を識別する方法が提案されている（非特許文献１）。又、非特許文献２では、表面筋電信号−手指関節角度の関係を、モデル式によりモデル化して、関数近似を行い、学習データから最小二乗法によりパラメータ推定を行う方法が提案されている。又、特許文献１では、表面筋電信号と指関節角度を一対一の対応付けで、複数の指の角度推定を行う方法が提案されている。
【特許文献１】特開平５−１１１８８１号公報
【非特許文献１】福田  修, 卜  楠, 辻  敏夫，「生ＥＭＧ信号による電動義手の制御」，計測自動制御学会論文集, Vol.40, No.11, pp. 1124-1131, 2004.
【非特許文献２】北村 徹, 辻内 伸好, 小泉 孝之, 「筋電信号による動作推定に基づくロボットハンドのマニピュレーション」, 日本機械学会論文集(C編), 73巻735号, pp.152-158, 2007
【発明の開示】
【発明が解決しようとする課題】
【０００３】
前記非特許文献１では、関節角度について、離散的なパターン識別を行うものであり、手指関節角度の推定を行うものではない。非特許文献２では、特定の指関節において単一の角度推定を行うものであり、多数の指関節への応用の可能性については具体的に提案されていない。特許文献１では、表面筋電信号と指関節角度の一対一の対応付けのみが行われており、表面筋電信号と手指関節角度の関係において、時間変換への対応の可能性については具体的に提案されていない。
【０００４】
従来は、表面筋電信号等の生体信号と、手指関節角度の関係の時間変化に対応可能であって、手指関節角度の推定を、複数の指関節、例えば５指関節について行われた指関節角度推定装置について提案されているものがない。
【０００５】
本発明の目的は、生体信号と指関節角度の関係の時間変化に対応可能であって、複数の指関節角度を連続的にリアルタイムに推定することができる指関節角度推定装置を提供することにある。
【課題を解決するための手段】
【０００６】
上記問題点を解決するために請求項１の発明は、各指関節が動いているときの単数又は複数種類の生体信号を入力して、各生体信号の特徴量を抽出し、抽出した該特徴量を一次元に並べて特徴量ベクトルとする特徴量抽出手段と、演算周期毎に入力ベクトルに基づいて、各指関節の角度（以下、指関節角度という）の推定を行い、推定した前記各指関節角度を出力ベクトルとして出力するニューラルネットワーク手段とを含み、前記入力ベクトルは、前記特徴量抽出手段が抽出した過去から現在までの時系列の特徴量ベクトルを一次元に並べたフィードフォワード成分と、出力のフィードバックとが一次元に並べられて構成され、前記フィードバック成分は、前記ニューラルネットワーク手段が、前記演算周期毎に現在まで推定していた時系列の前記出力ベクトルで構成されていることを特徴とする指関節角度推定装置を要旨とするものである。
【０００７】
なお、本明細書で生体信号とは、身体活動する際に起こる生理的な情報の信号であり、センサなどの検出手段により感知できる信号のことをいう。
請求項２の発明は、請求項１において、前記生体信号は、表面筋電信号であることを特徴とする。
【０００８】
請求項３の発明は、請求項２において、前記特徴量抽出手段は、前記表面筋電信号に対して積分区間を演算周期毎にずらして積分し、積分して得られた積分値を特徴量として抽出するものである。
【０００９】
請求項４の発明は、請求項２において、前記特徴量抽出手段は、前記表面筋電信号のパワースペクトルの所定のバンド幅毎の区間和を特徴量として抽出するものである。
請求項５の発明は、請求項２において、前記特徴量抽出手段は、前記表面筋電信号のパワースペクトルの平均周波数又は中央周波数を特徴量として抽出するものである。
【００１０】
請求項６の発明は、請求項１において、前記生体信号は、脳波であり、前記特徴量抽出手段は、前記脳波のα波及びβ波の少なくともいずれか１つの波の増減変化による事象関連同期化、及び事象関連脱同期化を特徴量として抽出するものである。
【００１１】
請求項７の発明は、請求項１において、前記生体信号は、脳部位のＭＲＩ信号であり、前記特徴量抽出手段は、運動を司る部位における血流中の還元ヘモグロビンの変化から活性化している場所を特徴量として抽出するものである。
【００１２】
請求項８の発明は、請求項１１において、前記生体信号は、頭部又は指関節に関係する筋肉を透過させた近赤外光を検出した信号であり、前記特徴量抽出手段は、頭部又は指関節に関係する筋肉組織を流れている血液中のヘモグロビン酸素化状態の変化を特徴量として抽出するものである。
【００１３】
請求項９の発明は、請求項１において、前記生体信号は、脳磁図であり、前記特徴量抽出手段は、前記脳磁図の磁場の変化を特徴量として抽出するものである。
請求項１０の発明は、請求項１において、前記生体信号は、指関節の動きによって腱が移動するのを捉えた信号であり、前記特徴量抽出手段は、前記腱の移動の変化を特徴量として抽出するものである。
【００１４】
請求項１１の発明は、請求項１において、前記生体信号は、指関節の動きによって筋肉の膨らみが変化するのを捉えた信号であり、前記特徴量抽出手段は、前記筋肉の膨らみの変化を特徴量として抽出するものである。
【００１５】
請求項１２の発明は、請求項１において、前記生体信号は、指関節の動きによって脳から、当該指の筋肉細胞までに伝わる神経信号を検出する信号であり、前記特徴量抽出手段は、前記指関節の動きに応じて変化する神経信号の変化を特徴量として抽出するものである。
【発明の効果】
【００１６】
請求項１乃至請求項１２の発明によれば、生体信号と指関節角度の関係の時間変化に対応可能であって、複数の指関節角度を連続的にリアルタイムに推定することができる指関節角度推定装置を提供することができる。
【発明を実施するための最良の形態】
【００１７】
以下、本発明の指関節角度推定装置を手指関節角度推定装置に具体化した一実施形態を図１〜図１１を参照して説明する。図１は手指関節角度推定装置の概略ブロック図である。
【００１８】
手指関節角度推定装置１０は、コンピュータからなり、Ａ／Ｄ１１、整流処理部１２、平滑処理部１３、積分処理部１４、正規化処理部１５、ニューラルネットワーク１６、及び変換処理部１９を備える。図１中、前記整流処理部１２、平滑処理部１３、積分処理部１４、正規化処理部１５、ニューラルネットワーク１６、及び変換処理部１７は、コンピュータが備える中央処理装置の機能ブロックを示しており、ハード構成を示すものではない。
【００１９】
手指関節角度推定装置１０は、人の前腕部Ｚに取り付けされた複数の電極２０及びＡ／Ｄ１１を介して表面筋電信号を入力する。なお、各電極２０から入力した信号をチャンネル信号という。Ａ／Ｄ１１は、各チャンネル毎に、共通のサンプリング周波数で取得した前記表面筋電信号（アナログ信号）をデジタル信号に変換する。
【００２０】
整流処理部１２は、Ａ／Ｄ１１からのデジタル信号を各チャンネル毎に全波整流する。平滑処理部１３はローパスフィルタからなり、整流処理部１２が全波整流した信号を、各チャンネル毎に所定周波数（例えば、50[Hz]）以上をカットオフして平滑する。
【００２１】
積分処理部１４は、平滑処理部１３から入力した各チャンネル信号に対して、積分区間を演算周期毎にずらしながら、該積分区間の振幅を積分して積分値を取得する（図９（ｂ）参照）。この積分値が特徴量に相当する。
【００２２】
正規化処理部１５は、積分処理部１４により積分処理された各チャンネル信号に対して、正規化を行う。本実施形態では、正規化処理は、各チャンネル信号の振幅を０〜１の範囲にするものであるが、この範囲に限定されるものではない。ここで、整流処理部１２、平滑処理部１３、積分処理部１４、及び正規化処理部１５は、特徴量抽出手段に相当する。
【００２３】
ニューラルネットワーク１６は、素子伝達関数としてシグモイド関数を用いた３層パーセプトロンニューラルネットワークから構成されている。３層パーセプトロンニューラルネットワークは図２に示すように入力層１６ａ、中間層１６ｂ、出力層１６ｃを備えている。ニューラルネットワーク１６におけるネットワークの学習は、５本指の各指関節角度が推定可能に一般的なバックプロパゲーション法により行われている。後述するフィードバック成分Ｂは、前記学習時には、過去の教師信号が用いられ、リアルタイムの手指関節角度の推定時には、過去にネットワークが実際に出力した値が用いられる。
【００２４】
ニューラルネットワーク１６の入力層１６ａを構成するユニットは、フィードフォワード成分Ａと、フィードバック成分Ｂからなる入力ベクトルＩｋ（図５参照）を入力するための個数分用意されている。又、出力層１６ｃを構成するユニットは、推定された５本指の各指の各指関節角度を１次元に並べた出力ベクトルＯｋとして出力する。なお、出力ベクトルＯｋに含まれる、各指の各指関節角度は範囲（０，１）に正規化されている。出力層１６ｃのユニットは、５本指の場合、本実施形態では、例えば、１本の指が３つの関節を有しているため、この３つの関節の屈曲角度（指関節角度）と、内外転の関節角度（指関節角度）を検出するため、各指毎の各関節角度の推定値が得られるように２０個分用意されているが、この数値に限定されるものではない。
【００２５】
図２、図３に示すように、フィードフォワード成分Ａは、Ａ０，……，Ａｎからなる。Ａ０は正規化処理部１５から入力された、現在（すなわち、最新）の各チャンネル信号が正規化された積分値（特徴量）Ｆｋを１次元に並べた特徴量ベクトルからなる。
【００２６】
Ａｎのｎは、現在からｎ（なお、ｎは正の整数）ステップ前の演算周期を意味する。従って、Ａｎは、ｎステップ前に正規化処理部１５から入力され、各チャンネル信号が正規化された積分値（特徴量）を１次元に並べた特徴量ベクトルである。具体的に説明すると、図３（ａ）に示すように、現在時刻ｔのフィードフォワード成分Ａは、次の演算周期の時刻ｔ＋１（図３（ｂ）参照）では、正規化処理部１５から特徴量ベクトルが入力されて、Ａ０が更新されるとともに、現在時刻ｔ時のＡ０の内容は、Ａ１の内容として更新される。
【００２７】
以下、同様にして、Ａ２〜Ａｎが更新される。なお、演算周期の時刻ｔのＡｎの内容は、演算周期の時刻ｔ＋１時では、破棄される。以下、図３（ｃ）に示すように、次の時刻ｔ＋２以降においても、同様に、Ａ０〜Ａｎの内容は更新される。上記のようにして、フィードフォワード成分Ａは、正規化された過去から現在までの時系列の特徴量ベクトルを一次元に並べたものとなる。
【００２８】
図２、図４に示すように、フィードバック成分Ｂは、Ｂ１，……，Ｂｈからなる。Ｂ１は出力層１６ｃからフィードバック入力された、現在の演算周期の１ステップ前の、各指の各関節角度を１次元に並べた出力ベクトルＯｋからなる。Ｂｈのｈは、現在からｈ（なお、ｈは正の整数）ステップ前の演算周期を意味する。従って、Ｂｈは、出力層１６ｃからフィードバックされた、現在の演算周期のｈステップ前の、各指の各関節角度を１次元に並べた出力ベクトルＯｋからなる。
【００２９】
具体的に説明すると、図４（ａ）に示すように、現在時刻ｔのフィードバック成分Ｂは、次の演算周期の時刻ｔ＋１（図４（ｂ）参照）では、出力層１６ｃから出力ベクトルＯｋがフィードバックされて、Ｂ１が更新されるとともに、現在時刻ｔ時のＢ１の内容は、Ｂ２の内容として更新される。以下、同様にして、Ｂ３〜Ｂｈが更新される。なお、演算周期の時刻ｔのＢｈの内容は、演算周期の時刻ｔ＋１時では、破棄される。以下、図４（ｃ）に示すように、次の時刻ｔ＋２以降においても、同様に、Ｂ１〜Ｂｈの内容は更新される。上記のようにして、フィードバック成分Ｂは、現在までに推定された時系列の出力ベクトルとなる。
【００３０】
変換処理部１７は、ニューラルネットワーク１６が推定した出力ベクトルＯｋを入力し、逆正規化処理して、各指の各指関節角度Ｎに変換する。式（１）は、逆正規化処理する場合の式である。
【００３１】
【数１】
2010125287.tif 000003
【００３２】
ここで、ＭＩＮは指関節角度の最小値であり、ＭＡＸは指関節角度の最大値である。
式（１）により、逆正規化処理して角度に変換した後、手指関節角度推定装置１０は、出力装置３０に出力する。出力装置３０は、例えば、表示ディプレイ、或いはプリンタ等にて構成されている。
【００３３】
（実施形態の作用）
さて、上記のように構成された手指関節角度推定装置１０の作用を図６〜図１１を参照して説明する。
【００３４】
手指関節角度推定装置１０に、電極２０を介して生体信号として入力された表面筋電信号は、所定のサンプリング周波数で取得される。このサンプリング周波数は、例えば１［kHz］であるが、この数値に限定されるものではない。図６は、Ａ／Ｄ変換されたチャンネル信号のチャートの例である。なお、図６の例では、サンプリング周波数は１［kHz］である。この信号は、脱力状態から握り、そのまま数秒経過した後にまた脱力状態に戻した時の、拇指屈曲に関与する筋肉の付近で計測した信号である。
【００３５】
次に、整流処理部１２において、入力された各チャンネル信号を全波整流化する。図７は、図６のチャンネル信号が全波整流処理されたチャートの例である。次に、平滑処理部１３において、各チャンネル信号を平滑処理する。図８は、図７のチャンネル信号が平滑処理されたチャートの例である。
【００３６】
次に、積分処理部１４において、平滑処理部１３から入力した各チャンネル信号に対して、積分区間を演算周期毎にずらしながら、該積分区間の振幅を積分して積分値を取得する。図９（ａ）は、図８の平滑処理されたチャンネル信号と同じであり、図９（ｂ）は、積分処理されたチャンネル信号の例のチャートである。正規化処理部１５は、積分処理部１４から出力された各チャンネル信号を０〜１の範囲となるように正規化する。図１０は、図９に示す信号が正規化処理されたチャンネル信号のチャートである。
【００３７】
なお、図６〜図９において、縦軸は、振幅、横軸は時間軸である。図１０において、縦軸は、正規化時の振幅、横軸は時間軸である。
ニューラルネットワーク１６は、フィードフォワード成分Ａ、フィードバック成分Ｂにて構成される入力ベクトルＩｋを演算周期毎に入力層１６ａに入力することにより、正規化されている５本指の各指関節角度を推定し、出力層１６ｃから前記出力ベクトルＯｋを出力する。この出力ベクトルＯｋを、変換処理部１７にて逆正規化して角度変換し、出力装置３０に出力する。
【００３８】
以下には、手に、指関節角度を計測するデータグローブ（図示しない）を取着して、５本の指の各指関節の角度を個々に検出可能とした状態で、かつ、前腕部Ｚに８個の電極２０を付けて、表面筋電信号を手指関節角度推定装置１０に入力して、指関節角度を推定した結果を図１１（ａ）〜（ｄ）に示す。
【００３９】
各図において、縦軸は角度、横軸は時間［ms］であり、点線は前記データグローブにより検出測定した実測値のチャートであり、実線は手指関節角度推定装置１０により推定した推定値のチャートである。
【００４０】
図１１（ａ）は、人指し指の内外転の角度のチャート、図１１（ｂ）は、人指し指の第１関節の屈曲角度のチャート、図１１（ｃ）は、人指し指の第２関節の屈曲角度のチャート、図１１（ｄ）は、人指し指の第３関節の屈曲角度のチャートを表わしている。各図に示すように、いずれも、手指関節角度推定装置１０にて推定した指関節角度は、実測値とほぼ一致し、又、リアルタイムに各関節角度がほぼ正確に得られていることが分かる。なお、図示はしないが、人指し指以外の他の指の各指関節角度も同様にほぼ正確に得られた。
【００４１】
比較例のために、ニューラルネットワーク１６の出力をフィードバック成分とせず、さらに現時刻前の生体信号を用いない場合の従来のニューラルネットワークの概念図を図１３に示す。
【００４２】
この比較例の場合、手にデータグローブ（図示しない）を取着して、５本の指の各指関節の角度を個々に検出可能とした状態で、かつ、前腕部Ｚに８個の電極２０を付けて、表面筋電信号を手指関節角度推定装置１０に入力して、指関節角度を推定した結果を図１４（ａ），（ｂ）に示す。図１４（ａ）、（ｂ）において、縦軸は角度、横軸は時間［ms］であり、点線は前記データグローブにより検出測定した実測値のチャートであり、実線は手指関節角度推定装置１０により推定した推定値のチャートである。
【００４３】
図１４（ａ）は、人指し指の内外転の角度のチャート、図１４（ｂ）は、人指し指の第１関節の屈曲角度のチャートである。同図に示すように推定した指関節角度は、実測値とは大きくかけ離れた値となっている。
【００４４】
以上詳述した本実施形態によれば、以下に記載する各効果を得ることができる。
（１） 本実施形態の手指関節角度推定装置１０は、各指関節が動いているときの複数の表面筋電信号（生体信号）を入力して、各表面筋電信号の特徴量を抽出し、抽出した該特徴量を一次元に並べて特徴量ベクトルとする整流処理部１２、平滑処理部１３、積分処理部１４、及び正規化処理部１５（特徴量抽出手段）を備える。又、手指関節角度推定装置１０は、演算周期毎に入力ベクトルＩｋに基づいて、各指関節角度の推定を行い、推定した各指関節角度を出力ベクトルとして出力するニューラルネットワーク１６（ニューラルネットワーク手段）とを備える。そして、入力ベクトルＩｋは、整流処理部１２、平滑処理部１３、積分処理部１４、及び正規化処理部１５が抽出した過去から現在までの時系列の特徴量ベクトルを一次元に並べたフィードフォワード成分Ａと、フィードバック成分Ｂとが一次元に並べられて構成され、フィードバック成分Ｂは、ニューラルネットワーク１６が、演算周期毎に現在まで推定していた時系列の出力ベクトルＯｋで構成されている。この結果、本実施形態の手指関節角度推定装置１０によれば、表面筋電信号と指関節角度の関係の時間変化に対応可能であって、５本の指の各指関節角度を連続的にリアルタイムに正確に推定することができる。
【００４５】
（２）  本実施形態の手指関節角度推定装置１０において、整流処理部１２、平滑処理部１３、積分処理部１４、及び正規化処理部では、表面筋電信号に対して積分区間を演算周期毎にずらし積分し、積分して得られた積分値を特徴量として抽出するようにした。この結果、演算周期毎に表面筋電信号の特徴量である積分値を取得することができる。
【００４６】
なお、本発明の実施形態は以下のように変更してもよい。
○  前記実施形態では、手指関節角度を推定する際に、生体信号として、前腕部において計測した表面筋電信号（ｓＥＭＧ）から得られる積分値を特徴量として用いたが、表面筋電信号の特徴量は、積分値に限定されるものではない。
【００４７】
例えば、表面筋電信号のパワースペクトルを使用してもよい。
この場合、前腕部において測定部位が複数ある場合、測定部位のチャンネル番号をｄとする。そして、各チャンネル毎に、かつサンプリング周期毎に得られたパワースペクトルを所定のバンド幅毎の区間和を求めて、ベクトル化すると、ある時に、あるチャンネルでサンプリングされたパワースペクトル（ベクトル）は式（２）のベクトルＦkdで表わすことができる。なお、所定のバンド幅は、限定されるものではなく、適宜に決定すればよい。前記パワースペクトルの区間和が特徴量となる。
【００４８】
【数２】
2010125287.tif 000004
【００４９】
式（２）中、Ｔは転置行列、ｋは時系列における順位、ｄはチャンネル番号、ｐは周波数の所定のバンド幅で区間分けした区分数である。
なお、パワースペクトルは、全波整流された表面筋電信号をフーリエ変換することにより求めることができる。なお、パワースペクトルの算出方法は、限定されるものではないが、フーリエ変換による算出法が一般的である。
【００５０】
図１２は、表面筋電信号をフーリエ変換した後のパワースペクトルを表わし、横軸が周波数、縦軸がパワースペクトルの強度である。パワースペクトルには、例えば、筋疲労の影響を反映しうる等、前記積分値を特徴量とする場合よりも得られ難い情報を含んでいる。
【００５１】
○  パワースペクトルは、波形であるため、定量的な評価が難しい場合がある。この場合、表面筋電信号の特徴量として、平均周波数（Mean  power  frequency：以下、ＭＮＦという）、或いは中央周波数（Median  power  frequency：以下、ＭＤＦという）を使用してもよい。
【００５２】
ＭＮＦの方が数学的に素直な量であるが、ＭＤＦには高周波成分に含まれるノイズの影響を受けにくい利点がある。
ＭＮＦ及びＭＤＦの定義式を次式（３）、式（４）で示す。
【００５３】
【数３】
2010125287.tif 000005
【００５４】
【数４】
2010125287.tif 000006
【００５５】
ここで、Ｐ（ｘ）はパワースペクトルである。ｆは周波数である。
○  又、生体信号は、表面筋電信号に限定されるものではない。表面筋電信号（ｓＥＭＧ）は、脳から発せられる神経信号に基づいて発生しているため、その基となる脳情報も，手指関節角度の推定に有用に使用することが可能である。なお、表面筋電信号（ｓＥＭＧ）は、筋繊維の活動電位の総和であり信号の振幅に大きな意味合いがあるが、脳情報に関しては各信号によって推奨されるべき信号処理法が異なる．
以下、脳情報を含むものとして、脳波（ＥＥＧ）、機能的磁気共鳴画像を形成するためのＭＲＩ信号、近赤外分光法（ＮＩＲＳ）、脳磁図（Magneto-Encephalo-Graphy :ＭＥＧ）を挙げることができる。
【００５６】
（１）脳波（Electroencephalogram：ＥＥＧ）
ＥＥＧは、古くからある非侵襲な計測手法で計測できる信号である。脳内の神経細胞の活動電位を、頭皮表面で計測した信号である。ここではα波及びβ波の少なくともいずれか１つの波の増減変化による事象関連同期化、及び事象関連脱同期化を特徴量として使用することができる。例えば、運動野が活動するとα波が弱まる事象関連脱同期化という現象があり、これは実際に運動をせず、想像する際にも生じる。又、α波が大きくなった時の事象関連同期化がある。又、ＥＥＧは後述するｆＭＲＩよりも時間分解能が高い利点がある。ＥＥＧで重要であるのが、前記事象関連脱同期化という現象である。この現象は特定の部位の脳波の、特定の周波数成分を変化させる。そのため、脳波から運動に関する要素を抽出するためには、高速フーリエ変換（FFT）やウェーブレット変換などの周波数解析手法によってパワースペクトルや周波数成分の分布を得られるような信号処理を施す。そして、本実施形態では、これらの信号処理をコンピュータにて行わせ、その信号処理結果をフィードフォワード成分とする。そして、前記実施形態のニューラルネットワーク１６に該フィードフォワード成分と、フィードバック成分を入力ベクトルとして入力するようにする。
【００５７】
（２）機能的磁気共鳴画像(ｆＭＲＩ)
神経活動がより大きく生じている脳部位を、血流中の還元ヘモグロビンの変化から推定する方法であり、ＭＲＩで測定し、脳部位から得られたＭＲＩ信号を生体信号として使用する。ここでは、脳部位において、運動を司る部位における前記血流中の還元ヘモグロビンの変化から活性化している場所を特徴量とするものである。
【００５８】
時間分解能は前述したＥＥＧに劣るが、空間分解能が格段に高い利点がある。計測機器は大がかりであり高価であるが、神経情報を詳細に読み取り機械制御信号に変換する方法として、非侵襲式の中では最も高性能である。
【００５９】
そして、本実施形態では、前記特徴量を有するＭＲＩ信号をフィードフォワード成分とする。そして、前記実施形態のニューラルネットワーク１６に該フィードフォワード成分と、フィードバック成分を入力ベクトルとして入力するようにする。
【００６０】
（３）近赤外分光法(ＮＩＲＳ)
近赤外分光法(ＮＩＲＳ)では、頭部や筋肉などの生体組織に対して透過性が高い近赤外光 (波長７００ｎｍ〜１０００ｎｍ) を外部から照射し、組織を透過してきた光を分析することにより、組織を流れている血液中のヘモグロビン酸素化状態を外部から調べる装置を使用する。生体信号は、該装置が検出した信号となる。又、前記組織を流れている血液中のヘモグロビン酸素化状態の変化を特徴量とする。
【００６１】
本変形例では、この脳が活動したときに一過性に生じる局所的な血流量の変化を光センサで検出するものである。ＥＥＧとｆＭＲＩの中間程度の特性を持つ。計測システムはＥＥＧと同等の規模であるが、非接触であるため、ＥＥＧと異なり、生体信号を身体に取付けする端子が必要でなく、該端子の取り付けの煩わしさがない利点がある。そして、本実施形態では、前記装置が検出した特徴量を有する信号をフィードフォワード成分とする。そして、前記実施形態のニューラルネットワーク１６に該フィードフォワード成分と、フィードバック成分を入力ベクトルとして入力するようにする。
【００６２】
（４）脳磁図（Magnetoencephalography：ＭＥＧ）
脳磁図（ＭＥＧ）は、前記ＥＥＧと似ているが、脳内の活動部位によって生じる磁場の変化を特徴量として検出するものである。電気信号はノイズが乗りやすいが、それに対して磁場は乗りにくい利点がある。
【００６３】
そして、本実施形態では、被検体頭部の周りに該磁場を検出するセンサ（すなわち、ピックアップコイル）を配置して、該磁場の変化を特徴量を有する信号を生成して、該信号に含まれる磁場の変化をフィードフォワード成分とする。そして、前記実施形態のニューラルネットワーク１６に該フィードフォワード成分と、フィードバック成分を入力ベクトルとして入力するようにする。
【００６４】
前記（２）〜（４）での検出は、基本的には脳内の活動部位の分布を検出するものである。ｆＭＲＩやＮＩＲＳは脳内の血流量の変化を、ＭＥＧは脳内磁場の分布を見ることができ、どの部位の神経細胞が活性化しているかを見ることができる。そのため、医学的な知見やあらかじめの計測において、手指の運動と、その際の脳内の活動部位との対応付けを行っておくものとする。
【００６５】
○  他の変形例として、下記のようにしてもよい。
関節は、筋肉の収縮により腱が移動し張力が発生することで駆動する。そのため、筋肉の収縮に伴う力こぶ等も、手指関節角度の推定に有効に利用できる。
【００６６】
（１） 腱の移動
手首の内側を見ながら指を動かすと、動かす関節によって異なる腱が移動している様子が伺える。腱の移動と関節の屈曲は直接的な関係があるため、腱の移動量から関節の回転角度を推定することが可能である。
【００６７】
このように、腱の移動を検出するために、すなわち、生体信号を得るために、２次元平面の圧力分布を検出するセンサ、例えば、分布型圧力センサを用いることができる。或いは、ビデオカメラを用いて、腱の移動を撮像し、該腱の移動に伴う影などの変化情報を含む画像信号を生体信号として取得する。
【００６８】
ここで、前記分布型圧力センサの検出信号により腱の移動の変化を特徴量としてもよい。この場合、前記特徴量を有する検出信号をフィードフォワード成分とする。そして、前記実施形態のニューラルネットワーク１６に該フィードフォワード成分と、フィードバック成分を入力ベクトルとして入力するようにする。
【００６９】
或いはビデオカメラが撮像した画像信号において、腱が移動したことを腱の移動を撮像した画像フレームの差分をとることにより、その変化を特徴量として得るようにしてもよい。この場合、前記変化がある部分を特徴量を有した画像信号が生体信号となる。そして、同様に、この場合、前記特徴量を有する生体信号をフィードフォワード成分とする。そして、前記実施形態のニューラルネットワーク１６に該フィードフォワード成分と、フィードバック成分を入力ベクトルとして入力するようにする。
【００７０】
（２）筋肉のふくらみ
腱の移動と筋肉の膨らみも密接な関係があるため、筋肉の膨らみもビデオカメラを使用して、手指関節の変化に応じた膨らみの変化を含む画像信号を生体信号として取得し、手指関節角度の推定に使用する。この場合、ビデオカメラが撮像した画像信号において、筋肉の膨らみの変化を撮像した画像フレームの差分をとることにより、その変化を特徴量として得るようにしてもよい。この場合、前記変化がある部分を特徴量を有した画像信号が生体信号となる。そして、同様に、この場合、前記特徴量を有する生体信号をフィードフォワード成分とする。そして、前記実施形態のニューラルネットワーク１６に該フィードフォワード成分と、フィードバック成分を入力ベクトルとして入力するようにする。
【００７１】
（３）神経信号
神経は微弱な電気信号を伝える。ｓＥＭＧを発生させる信号も電気信号であるため、指関節を変化させたときの脳から筋肉細胞まで伝わる信号を検出すれば、生体信号として有用な信号となる。具体的には、運動神経に端子を装着することにより、神経信号を検出することができる。神経を伝わる信号は、連続的なアナログ信号よりは、離散的なデジタル信号に近い。前腕部の筋群を支配する複数の神経を伝わる信号を計測し、どこの神経が、どれだけの時間信号を伝えたかを計測することにより、運動を反映しうる生体信号になる。そして、指関節の動きによって脳から、当該指の筋肉細胞までに伝わる神経信号を検出した信号において、その神経信号の変化が特徴量となる。
【００７２】
この場合、前記特徴量を有する生体信号をフィードフォワード成分とする。そして、前記実施形態のニューラルネットワーク１６に該フィードフォワード成分と、フィードバック成分を入力ベクトルとして入力するようにする。
【００７３】
○  前記実施形態では、手指関節角度を推定するようにしたが、手指関節角度に限定されるものではなく、足指関節角度を推定することも勿論可能である。
○  前記実施形態では、ニューラルネットワーク１６として、３層パーセプトロンニューラルネットワークにより構成したが、３層以上のパーセプトロンニューラルネットワークにより構成してもよい。又、ニューラルネットワーク１６はパーセプトロンに限定されるものではなく、他のニューラルネットワークにより構成してもよい。
【００７４】
○  前記実施形態では、１本の指が３つの関節を有しているため、これらの関節の屈曲角度（指関節角度）と、さらに、内外転の角度（指関節角度）を検出するために、すなわち、各関節角度を推定するためにニューラルネットワーク１６の出力層１６ｃのユニットを２０個としたが、推定に必要な関節角度に応じて、出力層１６ｃのユニットの数を減らしてもよい。
【００７５】
又、５本指の各関節角度の他に、手首関節等の他の関節角度を推定するために、逆に出力層１６ｃのユニットの数を増加させてもよい。
○  前記実施形態では、複数の生体信号（表面筋電信号）を取得するようにしたが、単数の表面筋電信号としてもよい。
【００７６】
○  前記実施形態では、単数種類の生体信号、すなわち、表面筋電信号を使用するようにしたが、上記各種の生体信号を複数種類組み合わせて使用することも可能である。
【図面の簡単な説明】
【００７７】
【図１】本発明を具体化した一実施形態の手指関節角度推定装置の概略ブロック図。
【図２】推定部のニューラルネットワークの構造の説明図。
【図３】（ａ）〜（ｃ）は、フィードフォワード成分Ａの説明図。
【図４】（ａ）〜（ｃ）は、フィードバック成分Ｂの説明図。
【図５】入力ベクトルの説明図。
【図６】Ａ／Ｄ変換されたチャンネル信号のチャート。
【図７】全波整流処理されたチャンネル信号のチャート。
【図８】平滑処理されたチャンネル信号のチャート。
【図９】（ａ）は図８と同じチャンネル信号のチャート、（ｂ）は積分処理されたチャンネル信号のチャート。
【図１０】正規化処理されたチャンネル信号のチャート。
【図１１】（ａ）は人指し指の内外転をデータグローブで検出測定した指関節角度のチャートと手指関節角度推定装置により推定した指関節角度のチャート、（ｂ）は人指し指の第１関節の角度をデータグローブで検出測定した指関節角度のチャートと手指関節角度推定装置により推定した指関節角度のチャート、（ｃ）は人指し指の第２関節の角度をデータグローブで検出測定した指関節角度のチャートと手指関節角度推定装置により推定した指関節角度のチャート、（ｄ）は人指し指の第３関節の角度をデータグローブで検出測定した指関節角度のチャートと手指関節角度推定装置により推定した指関節角度のチャート。
【図１２】表面筋電信号をフーリエ変換した後のパワースペクトルのチャート。
【図１３】比較例のニューラルネットワーク１６の出力をフィードバック成分としない場合の概念図。
【図１４】（ａ）は比較例における人指し指の内外転をデータグローブで検出測定した指関節角度のチャートと手指関節角度推定装置により推定した指関節角度のチャート、（ｂ）は比較例における人指し指の第１関節の角度をデータグローブで検出測定した指関節角度のチャートと手指関節角度推定装置により推定した指関節角度のチャート。
【符号の説明】
【００７８】
１０…手指関節角度推定装置（指関節角度推定装置）、
１２…整流処理部、１３…平滑処理部、１４…積分処理部、
１５…正規化処理部（平滑処理部１３、積分処理部１４とともに特徴量抽出手段を構成する）、
１６…ニューラルネットワーク（ニューラルネットワーク手段）、
１７…変換処理部、２０…電極、３０…出力装置。
【図１】
2010125287.tif 000007
【図２】
2010125287.tif 000008
【図３】
2010125287.tif 000009
【図４】
2010125287.tif 000010
【図５】
2010125287.tif 000011
【図６】
2010125287.tif 000012
【図７】
2010125287.tif 000013
【図８】
2010125287.tif 000014
【図９】
2010125287.tif 000015
【図１０】
2010125287.tif 000016
【図１１】
2010125287.tif 000017
【図１２】
2010125287.tif 000018
【図１３】
2010125287.tif 000019
【図１４】
2010125287.tif 000020