JP3735969B2 - Conducting action judging method and conducting action judging device - Google Patents
Conducting action judging method and conducting action judging device Download PDFInfo
- Publication number
- JP3735969B2 JP3735969B2 JP28646696A JP28646696A JP3735969B2 JP 3735969 B2 JP3735969 B2 JP 3735969B2 JP 28646696 A JP28646696 A JP 28646696A JP 28646696 A JP28646696 A JP 28646696A JP 3735969 B2 JP3735969 B2 JP 3735969B2
- Authority
- JP
- Japan
- Prior art keywords
- beat
- label
- label sequence
- hmm
- probability
- 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 - Fee Related
Links
- 238000000034 method Methods 0.000 claims description 22
- 238000004364 calculation method Methods 0.000 claims description 17
- 238000001514 detection method Methods 0.000 claims description 15
- 230000007704 transition Effects 0.000 description 23
- 238000010586 diagram Methods 0.000 description 15
- 238000004458 analytical method Methods 0.000 description 10
- 238000013528 artificial neural network Methods 0.000 description 4
- 238000005452 bending Methods 0.000 description 3
- 239000004020 conductor Substances 0.000 description 3
- 238000002372 labelling Methods 0.000 description 3
- 230000001133 acceleration Effects 0.000 description 2
- 238000007476 Maximum Likelihood Methods 0.000 description 1
- 238000010009 beating Methods 0.000 description 1
- 238000007621 cluster analysis Methods 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 235000015114 espresso Nutrition 0.000 description 1
- 230000002650 habitual effect Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012732 spatial analysis Methods 0.000 description 1
Images
Landscapes
- Electrophonic Musical Instruments (AREA)
- Image Analysis (AREA)
- Image Processing (AREA)
Description
【0001】
【発明の属する技術分野】
この発明は、動作に応じたセンサ出力をベクトル量子化し、これによって生成されたラベル系列を隠れマルコフモデル(HMM)が出力する確率を計算することにより、その動作の内容を判定する指揮動作判定方法および指揮動作判定装置に関する。
【0002】
【従来の技術】
人の動作を判定する装置や方法は従来より種々提案されている。その手法は、一般に、単純なフィルタや大小比較等の信号処理、二次元動作信号の角度や角度差などであった。
【0003】
【発明が解決しようとする課題】
しかし、人間の動作は極めて曖昧、且つ不安定なため、従来の単純な信号処理だけでは、検出や判別の精度、信頼性が低く、検出ミスや判別誤りが頻繁に発生した。
【0004】
このため、例えばこの動作判定を用いて音楽のテンポやダイナミクスを制御しようとした場合、
▲1▼指揮動作の認識率が極めて低いため、機械が認識し易い動作を操作者が習熟することが必要となり、ユーザがシステムに慣れるのに時間を要する。
▲2▼動作意図と異なる反応、すなわち誤認識が頻発し、安定した音楽演奏が困難である。
等の不都合が生じていた。
【0005】
この発明は、HMMを用いることにより指揮動作などの動作検出の精度や信頼性を向上した指揮動作判定方法および指揮動作判定装置を提供することを目的とする。
【0006】
【課題を解決するための手段】
この出願の請求項1の発明は、動作を検出するセンサから出力された検出値を一定時間毎にラベル化してラベル系列に編集し、特定の動作に対応するラベル系列を学習させた隠れマルコフモデルが該ラベル系列を出力する確率を計算し、この計算結果が所定の値を超えたとき前記センサが検出した動作の内容を判定することを特徴とする。
【0007】
この出願の請求項2の発明は、判定された動作の内容に基づいてラベルを生成し、これを後続の前記ラベル系列に挿入することを特徴とする。
【0008】
この出願の請求項3の発明は、動作を検出して検出値を出力するセンサと、該センサから入力される検出値を一定時間毎にラベル化する状態判定部と、前記状態判定部から一定時間毎に出力されるラベルをラベル系列に編集するラベル系列生成部と、特定の動作に対応するラベル系列を学習させた隠れマルコフモデルを記憶するHMM記憶部と、前記ラベル系列生成部からラベル系列を入力し、前記HMM記憶部に記憶されている隠れマルコフモデルが該ラベル系列を出力する確率を計算する確率計算部と、該確率計算部の計算結果が所定の値を超えたとき前記センサが検出した動作の内容を判定する動作判定部と、からなることを特徴とする。
【0009】
この出願の請求項4の発明は、前記ラベル系列生成部は、前記動作判定部が判定した動作の内容に基づくラベルを後続の前記ラベル系列に挿入する手段を含むことを特徴とする。
【0010】
【発明の実施の形態】
図1を参照してこの発明に用いられる隠れマルコフモデル(HMM)について説明する。HMMは、種々のラベル系列をそれぞれ所定の確率で出力するものであり、N個の状態(S1 ,S2 ,…,SN )をもち、一定周期ごとに状態を次々に遷移するとともに、その遷移の際にラベルを1つずつ出力するものである。つぎにどの状態に遷移するか、またその際にどのラベルを出力するか否かはそれぞれ「遷移確率」,「出力確率」によって確率的に決められている。また、HMMの出力側からみた特定のラベル系列が出力される確率を「ゆう度」といい、本願特許請求の範囲の「確率」はこの「ゆう度」を意味している。
【0011】
図1に示すHMMは、3つの状態S1 ,S2 ,S3 で構成され、2種類のラベルaとbのみからなるラベル系列を出力する。初期状態S1 からは0.3の確率でS1 自体に遷移する(その際に、ラベルaかbをそれぞれ0.8,0.2の確率で出力する)か、0.5の確率でS2 に遷移する(その際に常にaを出力する)か、0.2の確率で最終状態S3 に遷移する(その際に常にbを出力する)。状態S2 からは、0.4の確率でS2 自体に遷移する(その際に、aかbをそれぞれ0.3,0.7の確率で出力する)か、0.6の確率で最終状態S3 に遷移する(その際にaかbを等確率で出力する)。
【0012】
ここで、このHMMがラベル系列aabを出力する確率を考えてみる。このHMMが出力することのできる状態遷移の系列は、自己遷移を繰り返すことができるため、状態遷移数に制限を加えなければ無数にあるが、aabを出力する可能性のあるものはS1 S1 S2 S3 とS1 S2 S2 S3 とS1 S1 S1 S3 の3つだけであり、それぞれの状態遷移の系列をとる確率は、
0.3 ×0.8 ×0.5 ×1.0 ×0.6 ×0.5 =0.036
0.5 ×1.0 ×0.4 ×0.3 ×0.6 ×0.5 =0.018
0.3 ×0.8 ×0.3 ×0.8 ×0.2 ×1.0 =0.01152
であり、このHMMからラベル系列aabが出力される確率の合計、すなわち、ラベル系列aabのゆう度は、
0.036 +0.018 +0.01152 =0.06552
である。なお、HMMからaabが出力されても上記3種類の状態遷移の系列のいずれによるのもかは出力側から観測されないことから、この形態のマルコフモデルを隠れ(Hidden)マルコフモデルという。HMMは、従来より単語音声認識など音声認識の分野で用いられている。
【0013】
音声認識は、入力音声を10ms程度のフレーム時間毎にラベル化し、このラベル系列の出力確率をそれぞれ異なる単語の発音を出力するように学習された複数のHMMについて計算し、最もそのゆう度の高いHMMの単語をその入力音声の単語であると認識するというものである。この技術に関しては、電子情報通信学会誌1987年4月号352頁〜358頁「マルコフモデルによる音声認識(大河内正明)」に詳しい。
【0014】
次に、この実施形態においてHMMを用いて検出・判定する指揮法、すなわち、指揮棒の揺動動作について説明する。図2は3拍子を指揮するときの指揮棒の軌跡を示す図である。図3は4拍子を指揮する場合の指揮棒の軌跡を示す図である。図4は2拍子を指揮する場合の指揮棒の軌跡を示す図である。
【0015】
図2において同図(A)は通常(ノンレガート)、同図(B)はレガート、同図(C)は軽いスタッカート、同図(D)は強いスタッカートの場合の指揮棒の軌跡を示している。これらの図から明らかなように第1拍の指示動作は、上から下への振り下ろし動作であり、その動作の下端が第1拍の拍打点である。そして軽いスタッカート(同図(C))以外は振り下ろしの反動の短い振り上げ動作を伴う。また、第2拍の指示動作は、右方向への揺動動作であり、奏法によりその動作中に拍打点がある場合(ノンレガート(同図(A)),レガート(同図(B)))と、右向きの動作の右端に拍打点がある場合(スタッカート(同図(C),(D)))とがある。また、第3拍の指示動作は、右下から左上への振り上げ動作であり、軽いスタッカート(同図(C))がその終点(第1拍の振り下ろしの開始点)に拍打点があるのを除いては、その動作途中に拍打点がある。なお、図中拍打点をその拍番号の数字で示しているが□で囲んでいるものは指揮棒がその点で停止する拍打点であり、○で囲んでいるものは指揮棒がある速度でその点を通過または折り返しする拍打点である。
【0016】
図3において同図(A)は通常(ノンレガート)、同図(B)はレガート、同図(C)は軽いスタッカート、同図(D)は強いスタッカートの場合の指揮棒の軌跡を示している。4拍子の指揮法は3拍子の指揮法と類似しており、4拍子の第1拍指示動作と3拍子の第1拍指示動作、4拍子の第3拍指示動作と3拍子の第2拍指示動作、4拍子の第4拍指示動作と3拍子の第3拍指示動作がほぼ対応する。そして4拍子の第2拍指示動作は、第1拍指示動作の終了点から左方向への揺動動作であり、奏法によりその動作中に拍打点がある場合(ノンレガート(同図(A)),レガート(同図(B)))と、左向きの動作の左端に拍打点がある場合(スタッカート(同図(C),(D)))とがある。
【0017】
また、図4において、2拍子の拍指示はほぼ上下動で行われる。ノンレガート(同図(A)),レガート(同図(B))および強いスタッカート(同図(D))において、第1拍の指示動作は、上から下への振り下ろし動作と、その反動の短い振り上げ動作からなり,振り下ろし動作の下端が第1拍の拍打点である。また、第2拍の指示動作は、短い予備動作である振り下ろし動作と、下から上(第1拍指示動作の開始点)への振り上げ動作からなり,振り下ろし動作の下端が第2拍の拍打点である。
【0018】
HMMを用いて、指揮棒の揺動から上記指揮法を分析する指揮動作分析装置を図5に示す。また、図6は同指揮動作分析装置の状態判定部の状態判定動作を説明する図、図7は同指揮動作分析装置のHMM記憶部5に記憶されるHMMを例示する図である。
【0019】
この指揮動作分析装置は、センサ1,状態判定部2,レジスタ3,確率計算部4,HMM記憶部5,拍判定部6からなっており、拍判定部6が判定した内容は自動演奏装置7に入力される。センサ1は、操作者が手に持って上記指揮法にしたがって揺動させるコントローラに内蔵された角速度や加速度を検出するセンサである。指揮法に従って揺動させるコントローラとしては指揮棒状のものが一般的であるが、この他にハンドグリップ状のものや手に直接貼りつけるものなどが考えられる。このセンサ1は、図6(A)に示すような揺動角度の領域を検出するために用いられるため、たとえば、揺動の垂直方向成分を検出する垂直方向センサおよび揺動の水平方向成分を検出する水平方向センサからなっている。このような垂直・水平方向センサを有する指揮棒については、未公開であるが、先願の特願平7−110549号に詳細に説明されている。このセンサ1の検出値は状態判定部2に入力される。
【0020】
状態判定部2は、以下の動作を行う。▲1▼センサ1の出力を10msずつのフレームに分割する。▲2▼揺動速度(角速度)の領域を判断し、図6(A)の区分で各フレームをラベル付け(動作ラベル:l1〜l5)する。▲3▼ラベルをレジスタ3に入力する。この動作は10msのフレームクロック毎に繰り返し実行する。なお、図6(B)のラベルリストのうちl6〜l14 は拍判定部6が出力する拍ラベルである。
【0021】
なお、図6(A)の揺動動作のラベル付けは、一例であって本発明がこれに限定されるものではない。使用するセンシングシステム(コントローラやセンサの種類)、操作者、センサ部の把持方法等によって、入力動作に対応するセンサ出力が異なるため、これらに対応しラベル付けの精度を向上するためには、複数の操作者、複数のセンサ把持方法など様々な態様の拍指示動作のデータを大量に採取し、類似する動作のデータの代表点を求めてラベル付けをするようにすればよい。
【0022】
レジスタ3は図5(B)のような構成になっている。後述の拍判定部6が判定した拍判定情報(拍ラベル)を記憶する拍ラベルレジスタ30、前記状態判定部2から入力されるラベル(動作ラベル)を50個記憶する50段のシフトレジスタ31、および、前記拍ラベルと動作ラベルを連結して確率演算部4に入力する合成部32からなっている。シフトレジスタ31は、10msのフレームクロックで記憶内容をシフトしてゆき、状態判定部2から入力される動作ラベルを最新のものから50個(500ms分)記憶する。
【0023】
上記レジスタは、拍ラベルと動作ラベルを別に記憶し、拍ラベルが必ずラベル系列の先頭になるように拍ラベルと動作ラベルを連結するようにしているが、これは以下の理由による。
シフトレジスタ31の長さが1拍以上になると今回の拍打ち動作のみならず、前拍の拍打ち動作のラベルも含まれることになり分析が複雑になる。そこで、シフトレジスタ31の長さを500ms分に制限したが、そうすると、拍ラベルを動作ラベルと同様にシフトレジスタ31に時系列に入力したのでは、次の拍タイミングにはシフトレジスタ31から出てしまっていることになる。そこで、拍ラベルのみ別に記憶するようにした。
【0024】
ただし、同図(C)のようにシフトレジスタを十分1拍以上の長さに構成して拍ラベルを動作ラベルと同様に時系列に入力し、必ず前拍の拍ラベルと前拍打ちの動作ラベルがラベル系列中に時系列的に含まれるようにすることもできる。この場合には、前拍の拍打ち動作の動作も含めて分析することにより、分析は複雑になるがより正確に拍種類を判定することができる。
【0025】
ただし、シフトレジスタの段数やフレームクロックの周波数は上記実施形態に限定されるものではない。
【0026】
確率計算部4は、1個の拍ラベルと50個の動作ラベルからなる51個のラベルで構成されたラベル系列を出力する確率、すなわち、該ラベル系列のゆう度をHMM記憶部5に記憶されている全てのHMMについて計算する。HMM記憶部5には図7(A)に示すように各種拍打動作のラベル系列を出力する複数のHMMが記憶されている。これらのHMMは、認識対象である指揮動作が時間に変動する状態を、有限個の状態遷移確率でモデル化したものであり、自己遷移経路を持つ3〜4の状態数で構成され、学習によって状態遷移確率および各ラベルの出力確率を決定すればよい。確率計算部4が計算したゆう度は拍判定部6に出力される。
【0027】
ここで、図7(B),(C)は、4拍子の1拍目のラベル系列を学習したHMM、すなわち、M41の例を示す図である。同図(B)は、図5(B)の構成のレジスタから、必ず拍ラベルが先頭にあるラベル系列が入力される場合のHMMの構成を示し、同図(C)は図5(C)の構成のレジスタから、前拍の拍打ちの動作ラベル−前拍の拍ラベル−今回の拍打ちの動作ラベルからなるラベル系列が入力される場合のHMMの構成を示している。
図7(B)では拍ラベルは必ず先頭で1個のみであるため、HMMは、状態S1 から必ず(遷移確率=1で)状態S2 に遷移し、そのときいずれかの拍ラベルl6〜l14 を出力する。状態S2 および状態S3 では動作ラベルl1〜l5のみ出力する。
一方、同図(C)では前拍の拍打ちの動作ラベル,前拍の拍ラベルおよび今回の拍打ちの動作ラベルを分析するため状態数が4となり、全ての遷移(自己遷移を含む)時に、全てのラベルl1〜l14 を出力する可能性がある。
なお、HMMの構成は上記例に限定されるものではない。
【0028】
拍判定部6は、確率計算部4から入力された各HMMのゆう度を比較して最も高いものを選出し、その最も高いものが一定の閾値を越えている場合に、拍タイミングであり、そのときの拍種類(拍子および拍番号)は最高のゆう度を出力したHMMの拍種類であると判定する。入力されたゆう度のうち最も値の大きいものでも前記一定の閾値を越えていない場合には、拍タイミングでないとして何も出力しない。
【0029】
すなわち、レジスタ3は拍タイミングであるか否かに拘らず各フレームタイミング毎に51個のラベル系列を確率計算部4に出力し、確率計算部4はこの入力されたラベル系列に基づいて各HMMのゆう度をフレームタイミング毎に出力する。このように拍判定部6には拍タイミングであるか否かに拘らず全HMMのゆう度が入力されることになるが、拍タイミングに入力されたラベル系列に対するゆう度とそれ以外のタイミングに入力されたラベル系列に対するゆう度とでは絶対的な値が異なる。したがって、適当な閾値を定め、ゆう度がそれ以下の場合には、拍タイミングでないと判定することによって、拍判定部6は、最大ゆう度を出力したHMMに基づいて拍種類を判定するのに同期して、拍タイミングを検出することができる。
【0030】
拍判定部6が拍タイミングとその拍の拍種類を判定すると、その拍種類を自動演奏装置7に出力する。自動演奏装置7は、入力された拍タイミングとその拍種類情報に基づいて演奏の拍タイミング,拍種類が入力された拍タイミング,拍種類に一致するように演奏のテンポを制御する。さらに、拍判定部6はこの拍種類に応じた拍ラベル (l6〜l14)を生成してレジスタ3に入力する。レジスタ3は入力された拍ラベルを拍ラベルレジスタ30に記憶する。
【0031】
以上の構成により、指揮動作検出装置は、操作者の指揮による拍指示動作を検出して自動演奏装置7を制御することができる。この装置においては、拍判定部6の拍判定結果を拍ラベルに変換してレジスタ3に入力し、動作ラベルとは別のレジスタに記憶するようにしているが、動作ラベルと同様にシフトレジスタにその発生タイミングに応じた順位で記憶するようにしてもよい。
【0032】
なお、HMM記憶部5に記憶されるHMMに対してより高度な学習を施して認識に反映させることができる。なお、ここで学習とは「ある動作に対して付与されたラベル系列Lに対して当該動作を表現するHMMのモデルMのゆう度Pr(L:M)を最大にするように、モデルMのパラメータ(遷移確率と出力確率)を調整すること」である。この学習としては、
▲1▼ 特定個人向けのカスタマイズ:その個人のデータだけで代表点を計算しなおす方式
▲2▼ 一般化:より多くの人のデータを採取し、代表点を計算しなおす方式
▲3▼ 演奏中の微小チューニング:演奏者のデータが、相当するラベルに予め設定してある代表値から常にズレている場合は、代表値を定期的に微調整する。
【0033】
などがある。学習は、パラメータに適当な初期値を与えておき、データに基づく反復計算で収束させる方法で行われる。
【0034】
なお、HMMを用いた指揮法のモデル化は、ラベル,扱うパラメータの種類,HMMの構造,その他種々の要素を如何に決定するかにより種々の手法が考えられ、上記実施形態は、これら種々のモデル化手法の一例に過ぎない。
【0035】
また、扱うパラメータや状態ラベルを増やして、認識する動作の種類や音楽情報の種類や増やしたり、認識率を向上することも可能である。例えば、動作の大きさや速度からダイナミクスを認識したり、二次元平面内での動作軌跡曲率などを参照してレガート,スタッカート等を認識することもできる。すなわち、拍打時点での軌跡の曲率が小さい(滑らかな)動作であればレガート(スラー/エスプレッシーヴォ)であると判断することができ、曲率が大きい(クッキリした)動作であればスタッカートであると判断することができる。
【0036】
また、ラベル化に用いるパラメータとしては、揺動の方向と速度(角速度)を用いたが、指揮(拍指示)の軌跡の形状から主な動作方向の成分を求め、第1成分の軸を縦方向,第2成分の軸を横方向という変換を施して演奏者の把持方法の違いや、クセの成分を減らしたりしてもよい。
【0037】
また、揺動の方向,速度(角速度)以外に以下のようなパラメータを用いることもできる。
【0038】
二次元面内または三次元空間内の基準点からの角度,位置,速度,加速度等
時間域波形のピーク,谷,絶対値など、
前の拍の種類
前の拍打点(拍タイミング)からの差(角度,速度,位置など)
前の拍タイミングからの時間
前サンプルとの差
動作の中心から見た象限
これらのうちのいずれかを使ってもよいし、複数組み合わせて使ってもよい。また、これら複数パラメータの空間的偏りをクラスタ分析し、代表ベクトルを求めてラベルにしてもよい。
【0039】
また、上記指揮動作分析装置は、拍タイミングと拍種類を認識するための装置であるが、これの一段上位の階層の認識、すなわち、拍種類の連鎖による拍子の分析にHMMを適用することができる。例えば「2拍目→3拍目→1拍目と遷移したら今は3拍子を演奏中」などの認識を行わせる。この場合において、ラベルを出力せずに状態遷移するナル遷移を導入することで、全ての拍を打たなくても拍子を認識させることができる。例えば、3拍子の認識用のHMMにおいて、1拍目から3拍目へのナル遷移を許すと、必ずしも2拍目を打たなくても3拍子として認識可能となる。
【0040】
このように、この実施形態では、自動演奏のテンポなどを制御する指揮動作判定装置に本発明を適用しているが、指揮動作は時系列に一定のルールに基づいて繰り返される連続動作であるため、HMMの構造決定や学習が容易であり、高い判定精度を期待することができる。
【0041】
また、図5の装置は指揮動作の判定に用いるのみでなく、人や物の動作全般の判定に用いることが可能である。例えば、マルチメディアインタフェースへの応用として、バーチャルリアリティ(VR)における動作のインタフェースに適用することができる。VRで用いられるセンサとしては、三次元内での位置,角度なとを検出する三次元位置・角度センサ、指間接の曲げ角度を検出するグローブ型間接曲げ角度検出センサ(データグラブ)、スーツ型間接曲げ角度検出センサなどがある。また、カメラで撮影した動作画像を認識することもできる。図8は上記センサが検出した動作の状態をラベル化し(同図(A))、ゲーム内での動作を上記ラベルの時間遷移を学習させたHMM(同図(B))で判定する。例えば、パンチに対応するHMMは、手を握り(l8),肘を伸ばしつつあり(l2),肘を伸ばしきった(l4)という動作状態のラベル系列に対して高い確率を出力するように学習されている。
【0042】
また、図5の装置を手話認識に用いることも可能である。この場合においてもカメラやデータグラブなどで指の形状や手の位置を検出し、これを図9(A)のようにラベル化する。このラベル系列に基づいて同図(B)に示すようなHMMにより表現される単語を認識する。また、指の形状や手の位置に限定されず、体全体の動作に基づいて動作を検出するようにしてもよい。
【0043】
なお、動作認識の手法としてはHMM以外に、ファジィ制御やニューラルネットワークが知られているが、ファジィ制御は、検出/判別ルールを全て記述する必要があるのに対してHMMはルールの記述が不要で、学習によってルールを内蔵することができるためシステム構築が容易である。また、ニューラルネットワークでは、ネットワークに学習させるために非常に複雑な演算が必要であるが、HMMは単純な演算で学習を行わせることができ、学習が容易である。したがって、これらファジィ制御やニューラルネットワークよりも動作認識に有効である。
【0044】
また、HMMは、実際の動作の揺らぎをファジィ制御やニューラルネットワークよりも正確に反映することができる。これは、出力確率がその状態が発生する確率的な揺らぎに対応し、遷移確率が時間軸方向の揺らぎに対応しているからである。また構造が簡略であるため、統計理論や情報理論等による理論的展開がし易く、確率の概念を用いて上位階層レベルの認識に統合することが容易である。
【0045】
【発明の効果】
以上のようにこの発明によれば、HMMに特定の動作を学習させ、動作に基づくセンサ出力をラベル化して該HMMのゆう度を算出するようにしたことにより、動作の速度や動きのスペクトルまたはそれらの揺らぎを吸収して正確な動作判定が可能になる。
【0046】
また、動作判定結果に基づいてラベルを生成し、これをHMMに与えるラベル系列に含めることにより、過去の動作内容が次の動作内容の判定に影響を与えることになり、連続した動作を判定する場合に高い精度の動作判定が可能になる。
【図面の簡単な説明】
【図1】HMMの簡単な例を示す図
【図2】3拍子を指揮するときの指揮棒の軌跡を示す図
【図3】4拍子を指揮する場合の指揮棒の軌跡を示す図
【図4】2拍子を指揮する場合の指揮棒の軌跡を示す図
【図5】この発明の実施形態である指揮動作分析装置の構成を示す図
【図6】同指揮動作分析装置の状態判定部の状態判定動作を説明する図
【図7】同指揮動作分析装置のHMM記憶部に記憶されるHMMを例示する図
【図8】本発明の動作分析装置をゲームに適用した場合のラベルとHMMを示す図
【図9】本発明の動作分析装置を手話認識に適用した場合のラベルとHMMを示す図
【符号の説明】
1…センサ、2…状態判定部、3…レジスタ、4…確率計算部、
5…HMM記憶部、6…拍判定部、7…自動演奏装置[0001]
BACKGROUND OF THE INVENTION
This invention is vector-quantized sensor output corresponding to the operation, by this hidden Markov model generated label series by (HMM) to compute the probability of output, commanding operation determination method for determining the content of the operation And a command action determination device.
[0002]
[Prior art]
Various devices and methods for determining human movement have been proposed. The method is generally signal processing such as simple filter and size comparison, angle and angle difference of two-dimensional motion signal, and the like.
[0003]
[Problems to be solved by the invention]
However, since human operations are extremely ambiguous and unstable, the accuracy and reliability of detection and discrimination are low only with conventional simple signal processing, and detection errors and discrimination errors frequently occur.
[0004]
For this reason, for example, when trying to control the tempo and dynamics of music using this motion determination,
{Circle around (1)} Since the recognition rate of the commanding action is extremely low, it is necessary for the operator to become familiar with the action that the machine can easily recognize, and it takes time for the user to get used to the system.
{Circle around (2)} A reaction different from the intended operation, that is, misrecognition occurs frequently, and stable music performance is difficult.
Inconvenience such as occurred.
[0005]
An object of the present invention is to provide a command operation determination method and a command operation determination device that improve the accuracy and reliability of operation detection such as a command operation by using an HMM.
[0006]
[Means for Solving the Problems]
The invention of claim 1 of this application is a hidden Markov model in which detection values output from a sensor for detecting motion are labeled at regular intervals, edited into a label sequence, and a label sequence corresponding to a specific motion is learned. Calculates the probability of outputting the label sequence, and determines the content of the operation detected by the sensor when the calculation result exceeds a predetermined value .
[0007]
The invention of
[0008]
The invention of
[0009]
The invention of
[0010]
DETAILED DESCRIPTION OF THE INVENTION
A hidden Markov model (HMM) used in the present invention will be described with reference to FIG. The HMM outputs various label sequences with predetermined probabilities, has N states (S 1 , S 2 ,..., S N ), changes states one after another at regular intervals, At the time of the transition, labels are output one by one. Next, which state is to be transited and which label is to be output at that time are determined probabilistically by “transition probability” and “output probability”, respectively. Further, the probability that a specific label sequence as viewed from the output side of the HMM is referred to as “likelihood”, and “probability” in the claims of the present application means this “likelihood”.
[0011]
The HMM shown in FIG. 1 is composed of three states S 1 , S 2 , S 3 and outputs a label sequence consisting of only two types of labels a and b. Transition from the initial state S 1 to S 1 itself with a probability of 0.3 (in this case, labels a and b are output with a probability of 0.8 and 0.2, respectively) or with a probability of 0.5 Transition to S 2 (always outputs a) or transition to final state S 3 with a probability of 0.2 (always outputs b). From state S 2 , transition to S 2 itself with a probability of 0.4 (in which case a or b is output with a probability of 0.3 or 0.7, respectively), or a final probability of 0.6 Transition to the state S 3 (in which case a or b is output with equal probability).
[0012]
Now consider the probability that this HMM will output the label sequence aab. Since the sequence of state transitions that can be output by the HMM can repeat self-transitions, there are an infinite number without limiting the number of state transitions, but there is a possibility of outputting aab S 1 S There are only 1 S 2 S 3 and S 1 S 2 S 2 S 3 and S 1 S 1 S 1 S 3 , and the probability of taking each state transition sequence is as follows:
0.3 × 0.8 × 0.5 × 1.0 × 0.6 × 0.5 = 0.036
0.5 × 1.0 × 0.4 × 0.3 × 0.6 × 0.5 = 0.018
0.3 × 0.8 × 0.3 × 0.8 × 0.2 × 1.0 = 0.01152
The total probability that the label sequence aab is output from the HMM, that is, the likelihood of the label sequence aab is
0.036 +0.018 +0.01152 = 0.06552
It is. Note that even if aab is output from the HMM, it is not observed from the output side which of the above three types of state transition series, so this form of Markov model is referred to as a Hidden Markov model. The HMM is conventionally used in the field of speech recognition such as word speech recognition.
[0013]
In speech recognition, input speech is labeled every frame time of about 10 ms, and the output probability of this label sequence is calculated for a plurality of HMMs learned to output pronunciations of different words, and the likelihood is highest. The HMM word is recognized as the word of the input speech. This technology is described in detail in the Journal of the Institute of Electronics, Information and Communication Engineers, April 1987, pages 352 to 358 “Speech Recognition by Markov Model (Masaaki Okouchi)”.
[0014]
Next, a command method for detecting and determining using the HMM in this embodiment, that is, a swinging motion of the command rod will be described. FIG. 2 is a diagram showing the trajectory of the conductor when conducting triple time. FIG. 3 is a diagram showing the trajectory of a conductor when conducting a 4-beat. FIG. 4 is a diagram showing the trajectory of the conductor when conducting a double time signature.
[0015]
2A is a normal (non-legato), FIG. 2B is a legato, FIG. 2C is a light staccato, and FIG. 2D is a trajectory of a baton for a strong staccato. Yes. As is clear from these figures, the instruction operation for the first beat is a swing-down operation from top to bottom, and the lower end of the operation is the beat point of the first beat. Except for a light staccato ((C) in the same figure), a swing-up operation with a short recoil is accompanied. In addition, the instruction operation for the second beat is a swinging operation in the right direction, and when there is a beat point during the operation due to the playing style (non-legato (Fig. (A)), legato (Fig. (B)) ) And when there is a beat point at the right end of the rightward movement (staccato ((C), (D) in the figure)). In addition, the instruction operation for the third beat is a swing-up operation from the lower right to the upper left, and the light staccato ((C) in the figure) has a beat point at its end point (starting point of the first beat swinging down). Except for, there is a beat point in the middle of the operation. In the figure, the beat point is indicated by the number of the beat number, but the one surrounded by □ is the beat point where the baton stops at that point, and the one surrounded by ○ is the speed at which the baton is located. It is a beat point that passes or turns around the point.
[0016]
In FIG. 3, (A) is a normal (non-legato), (B) is a legato, (C) is a light staccato, and (D) is a trajectory of a baton for a strong staccato. Yes. The 4-beat command method is similar to the 3-beat command method, and the 4-beat first beat instruction operation, the 3-beat first beat instruction operation, the 4-beat third beat instruction operation, and the 3-beat second beat. The instruction operation, the fourth beat instruction operation of 4 beats, and the third beat instruction operation of 3 beats substantially correspond. Then, the second beat instructing action of 4 beats is a swinging action in the left direction from the end point of the first beat instructing action, and when there is a beat point during the action due to the rendition style (non-legato (FIG. (A)) ), Legato (FIG. (B))) and a beat point at the left end of the leftward movement (staccato (FIG. (C), (D))).
[0017]
Further, in FIG. 4, a two-beat instruction is almost vertically moved. In non-legato (Fig. (A)), legato (Fig. (B)), and strong staccato (Fig. (D)), the first beat indicating operation is a swing-down operation from top to bottom and its reaction The lower end of the swing-down motion is the beat point of the first beat. In addition, the second beat instruction operation includes a swing-down operation, which is a short preliminary operation, and a swing-up operation from the bottom to the top (starting point of the first beat instruction operation). The lower end of the swing-down operation is the second beat. This is the beat point.
[0018]
FIG. 5 shows a command action analysis apparatus that analyzes the command method from the swing of the command stick using the HMM. FIG. 6 is a diagram for explaining the state determination operation of the state determination unit of the command action analyzer, and FIG. 7 is a diagram illustrating the HMM stored in the HMM storage unit 5 of the command action analyzer.
[0019]
This commanding action analysis apparatus comprises a sensor 1, a
[0020]
The
[0021]
Note that the labeling of the swing motion in FIG. 6A is an example, and the present invention is not limited to this. The sensor output corresponding to the input operation differs depending on the sensing system (type of controller or sensor) used, the operator, the gripping method of the sensor unit, and so on. It is only necessary to collect a large amount of beat instruction motion data in various modes such as the operator and a plurality of sensor gripping methods, and to obtain a representative point of similar motion data for labeling.
[0022]
The
[0023]
The register stores the beat label and the action label separately, and connects the beat label and the action label so that the beat label is always at the head of the label series, for the following reason.
When the length of the
[0024]
However, as shown in (C) in the figure, the shift register is configured with a length of one or more beats, and beat labels are input in time series in the same way as operation labels. The labels can be included in the label series in time series. In this case, by analyzing including the operation of the beat operation of the previous beat, the analysis becomes complicated, but the beat type can be determined more accurately.
[0025]
However, the number of stages of the shift register and the frequency of the frame clock are not limited to the above embodiment.
[0026]
The
[0027]
Here, FIG. 7 (B), (C) is, HMM learned the quadruple time of the first beat of the label sequence, i.e., a diagram illustrating an example of M 41. FIG. 5B shows the configuration of the HMM when a label series having a beat label at the head is always input from the register having the configuration shown in FIG. 5B. FIG. 5C shows the configuration of FIG. The configuration of the HMM is shown in the case where a label series consisting of an operation label of the previous beat, a beat label of the previous beat, and an operation label of the current beat is input from the register having the configuration of FIG.
In FIG. 7B, since there is always only one beat label at the head, the HMM always transitions from state S 1 to state S 2 (with transition probability = 1), at which time any beat label l 6 ~ L 14 is output. Only the state S 2 and state S operation label l 1 to l 5 In 3 outputs.
On the other hand, in FIG. 5C, the number of states is 4 in order to analyze the previous beat beat action label, the previous beat beat label, and the current beat action label, and at all transitions (including self-transition). There is a possibility that all the labels l 1 to l 14 are output.
The configuration of the HMM is not limited to the above example.
[0028]
The
[0029]
That is, the
[0030]
When the
[0031]
With the above configuration, the conducting action detection device can control the
[0032]
Note that more advanced learning can be applied to the HMM stored in the HMM storage unit 5 and reflected in the recognition. Note that learning here means that “the likelihood M Pr (L: M) of the model M of the HMM that expresses the action for the label series L given to a certain action is maximized. Adjusting parameters (transition probabilities and output probabilities). As this learning,
(1) Customization for a specific individual: A method for recalculating representative points with only the individual data (2) Generalization: A method for collecting more people's data and recalculating the representative points (3) Performing Fine tuning: When the performer's data is always shifted from the representative value preset in the corresponding label, the representative value is finely adjusted periodically.
[0033]
and so on. Learning is performed by a method in which an appropriate initial value is given to a parameter and convergence is made by iterative calculation based on data.
[0034]
It should be noted that various methods can be considered for modeling the command method using the HMM depending on how the label, the type of parameter to be handled, the structure of the HMM, and other various elements are determined. It is just an example of a modeling technique.
[0035]
It is also possible to increase the number of parameters and status labels to be handled to increase the types of actions to be recognized, the types of music information, and the recognition rate. For example, the dynamics can be recognized from the magnitude and speed of the motion, or legato, staccato, etc. can be recognized by referring to the motion trajectory curvature in the two-dimensional plane. In other words, if the curvature of the trajectory at the time of beat is small (smooth), it can be determined that it is legato (slur / espresso), and if the curvature is large (clear), it is staccato. It can be judged.
[0036]
As the parameters used for labeling, the direction of swing and velocity (angular velocity) were used, but the main motion direction component was obtained from the shape of the command (beat instruction) trajectory, and the axis of the first component was The direction and the axis of the second component may be converted to the horizontal direction to reduce the difference in the player's gripping method and the habitual component.
[0037]
In addition to the swing direction and velocity (angular velocity), the following parameters can be used.
[0038]
Peak, valley, absolute value, etc. of time domain waveforms such as angle, position, velocity, acceleration from reference point in 2D plane or 3D space,
Difference from previous beat type (beat timing) (angle, speed, position, etc.)
Any of these quadrants as seen from the center of the difference operation from the previous beat timing from the previous beat timing may be used, or a plurality of them may be used in combination. Alternatively, a spatial analysis of these multiple parameters may be subjected to cluster analysis to obtain a representative vector as a label.
[0039]
In addition, the command action analysis device is a device for recognizing beat timing and beat type. However, HMM can be applied to recognition of a layer one level higher, that is, to analyze time signatures by chaining beat types. it can. For example, recognition is performed such as “when the transition from the second beat → the third beat → the first beat is performed, the three beats are now being played”. In this case, the time signature can be recognized without hitting all the beats by introducing a null transition that causes a state transition without outputting a label. For example, in the HMM for 3-beat recognition, if the null transition from the 1st beat to the 3rd beat is allowed, the 3-beat can be recognized without necessarily hitting the 2nd beat.
[0040]
As described above, in this embodiment, the present invention is applied to the commanding action determination device that controls the tempo of automatic performance, but the commanding action is a continuous action that is repeated based on a certain rule in time series. The structure determination and learning of the HMM are easy, and high determination accuracy can be expected.
[0041]
Further, the apparatus of FIG. 5 can be used not only for determining a commanding action, but also for determining a general action of a person or an object. For example, the present invention can be applied to an interface of operation in virtual reality (VR) as an application to a multimedia interface. Sensors used in VR include a three-dimensional position / angle sensor that detects the position and angle in three dimensions, a glove-type indirect bending angle detection sensor (data grab) that detects an indirect finger bending angle, and a suit type. There are indirect bending angle detection sensors. It is also possible to recognize an operation image taken by a camera. FIG. 8 labels the motion state detected by the sensor (FIG. 8A), and determines the motion in the game by the HMM (FIG. 5B) that learned the time transition of the label. For example, the HMM corresponding to the punch outputs a high probability for the label sequence in the operating state of grasping the hand (l 8 ), stretching the elbow (l 2 ), and fully extending the elbow (l 4 ). To be learned.
[0042]
It is also possible to use the apparatus shown in FIG. 5 for sign language recognition. Also in this case, the shape of the finger and the position of the hand are detected by a camera, a data grab, etc., and are labeled as shown in FIG. Based on this label sequence, a word expressed by an HMM as shown in FIG. Further, the movement is not limited to the shape of the finger or the position of the hand, and the movement may be detected based on the movement of the whole body.
[0043]
In addition to HMM, fuzzy control and neural network are known as motion recognition methods. However, fuzzy control needs to describe all detection / discrimination rules, whereas HMM does not require rule description. Therefore, it is easy to build a system because rules can be built in by learning. In addition, in the neural network, a very complicated operation is required for the network to learn, but in the HMM, the learning can be performed by a simple operation and the learning is easy. Therefore, it is more effective for motion recognition than these fuzzy controls and neural networks.
[0044]
Further, the HMM can reflect the fluctuation of the actual operation more accurately than the fuzzy control or the neural network. This is because the output probability corresponds to a probabilistic fluctuation in which the state occurs, and the transition probability corresponds to a fluctuation in the time axis direction. Moreover, since the structure is simple, it is easy to carry out theoretical development by statistical theory, information theory, etc., and it is easy to integrate into the recognition of the upper hierarchy level using the concept of probability.
[0045]
【The invention's effect】
As described above, according to the present invention, the HMM learns a specific operation, labels the sensor output based on the operation, and calculates the likelihood of the HMM. By absorbing these fluctuations, it is possible to accurately determine the operation.
[0046]
In addition, by generating a label based on the operation determination result and including it in the label sequence given to the HMM, the past operation content affects the determination of the next operation content, and the continuous operation is determined. In this case, it is possible to determine the operation with high accuracy.
[Brief description of the drawings]
[Fig. 1] A diagram showing a simple example of an HMM. [Fig. 2] A diagram showing a trajectory of a command bar when conducting a three time signature. [Fig. 3] A diagram showing a trajectory of a command rod when conducting a four time signature. [Fig. 4 is a diagram showing the trajectory of a command stick when commanding two beats. FIG. 5 is a diagram showing a configuration of a command action analysis apparatus according to an embodiment of the present invention. FIG. 6 is a diagram of a state determination unit of the command action analysis apparatus. FIG. 7 is a diagram illustrating the state determination operation. FIG. 7 is a diagram illustrating an HMM stored in the HMM storage unit of the command motion analysis device. FIG. 8 is a label and an HMM when the motion analysis device of the present invention is applied to a game. FIG. 9 is a diagram showing labels and HMMs when the motion analysis apparatus of the present invention is applied to sign language recognition.
DESCRIPTION OF SYMBOLS 1 ... Sensor, 2 ... State determination part, 3 ... Register, 4 ... Probability calculation part,
5 ... HMM storage unit, 6 ... beat determination unit, 7 ... automatic performance device
Claims (4)
該センサから入力される検出値を一定時間毎にラベル化する状態判定部と、
前記状態判定部から一定時間毎に出力されるラベルをラベル系列に編集するラベル系列生成部と、
特定の動作に対応するラベル系列を学習させた隠れマルコフモデルを記憶するHMM記憶部と、
前記ラベル系列生成部からラベル系列を入力し、前記HMM記憶部に記憶されている隠れマルコフモデルが該ラベル系列を出力する確率を計算する確率計算部と、
該確率計算部の計算結果が所定の値を超えたとき前記センサが検出した動作の内容を判定する動作判定部と、
からなる指揮動作判定装置。A sensor that detects movement and outputs a detection value;
A state determination unit that labels detection values input from the sensor at regular intervals;
A label sequence generation unit for editing labels output at regular intervals from the state determination unit into label sequences;
An HMM storage unit that stores a hidden Markov model in which a label sequence corresponding to a specific action is learned;
A probability calculation unit that inputs a label sequence from the label sequence generation unit and calculates a probability that a hidden Markov model stored in the HMM storage unit outputs the label sequence;
An operation determination unit that determines the content of the operation detected by the sensor when the calculation result of the probability calculation unit exceeds a predetermined value ;
A commanding action determination device consisting of
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP28646696A JP3735969B2 (en) | 1995-11-02 | 1996-10-29 | Conducting action judging method and conducting action judging device |
Applications Claiming Priority (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP7-285774 | 1995-11-02 | ||
| JP28577495 | 1995-11-02 | ||
| JP28646696A JP3735969B2 (en) | 1995-11-02 | 1996-10-29 | Conducting action judging method and conducting action judging device |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH09185716A JPH09185716A (en) | 1997-07-15 |
| JP3735969B2 true JP3735969B2 (en) | 2006-01-18 |
Family
ID=26556024
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP28646696A Expired - Fee Related JP3735969B2 (en) | 1995-11-02 | 1996-10-29 | Conducting action judging method and conducting action judging device |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP3735969B2 (en) |
Families Citing this family (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2004164563A (en) | 2002-09-26 | 2004-06-10 | Toshiba Corp | Image analysis method, image analysis device, image analysis program |
| JP5161643B2 (en) * | 2008-04-23 | 2013-03-13 | 富士重工業株式会社 | Safe driving support system |
| JP5717553B2 (en) * | 2011-06-19 | 2015-05-13 | 将文 萩原 | Suspicious person detection method and apparatus, abnormal event detection method and apparatus |
| GB2492748B (en) * | 2011-07-04 | 2014-05-07 | Jaguar Land Rover Ltd | Vehicle control system and method for controlling a vehicle |
| JP2015130070A (en) * | 2014-01-07 | 2015-07-16 | 富士通株式会社 | Detection program, detection method, and detection apparatus |
| JP7383943B2 (en) * | 2019-09-06 | 2023-11-21 | ヤマハ株式会社 | Control system, control method, and program |
| JP6258442B1 (en) * | 2016-10-28 | 2018-01-10 | 三菱電機インフォメーションシステムズ株式会社 | Action specifying device, action specifying method, and action specifying program |
-
1996
- 1996-10-29 JP JP28646696A patent/JP3735969B2/en not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| JPH09185716A (en) | 1997-07-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US5808219A (en) | Motion discrimination method and device using a hidden markov model | |
| JP4430368B2 (en) | Method and apparatus for analyzing gestures made in free space | |
| Morita et al. | A computer music system that follows a human conductor | |
| Deng et al. | A statistical approach to automatic speech recognition using the atomic speech units constructed from overlapping articulatory features | |
| Xia et al. | Autonomous robot dancing driven by beats and emotions of music | |
| US20090066641A1 (en) | Methods and Systems for Interpretation and Processing of Data Streams | |
| CN110245718A (en) | A kind of Human bodys' response method based on joint time-domain and frequency-domain feature | |
| CN110216677B (en) | Method for generating and optimizing mechanical dexterous hand playing piano action sequence | |
| US10846519B2 (en) | Control system and control method | |
| JP7383943B2 (en) | Control system, control method, and program | |
| JP3735969B2 (en) | Conducting action judging method and conducting action judging device | |
| Fang et al. | Signer-independent continuous sign language recognition based on SRN/HMM | |
| Wang et al. | Data-driven simulation framework for expressive piano playing by anthropomorphic hand with variable passive properties | |
| KR20230054522A (en) | Augmented reality rehabilitation training system applied with hand gesture recognition improvement technology | |
| CN121300636B (en) | A method and system for motion-sensing interaction recognition based on skeletal coordinate points | |
| Fang et al. | Signer-independent sign language recognition based on SOFM/HMM | |
| Gillian | Gesture recognition for musician computer interaction | |
| CN120183047A (en) | A sign language translation training system based on gesture tracking and a sign language tracking device | |
| CN114387616A (en) | A music performance system and method based on gesture action recognition | |
| Shi et al. | Optimized fingering planning for automatic piano playing using dual-arm robot system | |
| WO2024149988A1 (en) | Apparatus and technique for the rich representation of impacts and dynamic interaction on surfaces through sensors and real-time machine learning | |
| KR100463819B1 (en) | Synchronization Method of Rhythmic Motion with Music based on Beat Analysis | |
| JPH1153104A (en) | Pen-type input device and pattern recognition method for pen-type input device | |
| JP3627321B2 (en) | Performance control device | |
| CN114742121A (en) | Neural network training method, 3D piano playing generation method and system |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20050411 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20050426 |
|
| A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20050627 |
|
| 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: 20051004 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20051017 |
|
| R150 | Certificate of patent (=grant) or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| S531 | Written request for registration of change of domicile |
Free format text: JAPANESE INTERMEDIATE CODE: R313532 |
|
| R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
| FPAY | Renewal fee payment (prs date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20081104 Year of fee payment: 3 |
|
| FPAY | Renewal fee payment (prs date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20091104 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (prs date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20101104 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (prs date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20101104 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (prs date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20111104 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (prs date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20111104 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (prs date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20121104 Year of fee payment: 7 |
|
| FPAY | Renewal fee payment (prs date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20121104 Year of fee payment: 7 |
|
| FPAY | Renewal fee payment (prs date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20131104 Year of fee payment: 8 |
|
| LAPS | Cancellation because of no payment of annual fees |