JP5623759B2 - 筋電図信号に基づいた脳内の並列運動制御機能の同定及び評価法 - Google Patents
筋電図信号に基づいた脳内の並列運動制御機能の同定及び評価法 Download PDFInfo
- Publication number
- JP5623759B2 JP5623759B2 JP2010042301A JP2010042301A JP5623759B2 JP 5623759 B2 JP5623759 B2 JP 5623759B2 JP 2010042301 A JP2010042301 A JP 2010042301A JP 2010042301 A JP2010042301 A JP 2010042301A JP 5623759 B2 JP5623759 B2 JP 5623759B2
- Authority
- JP
- Japan
- Prior art keywords
- motion
- frequency
- electromyogram
- velocity
- joint
- 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
- 210000004556 brain Anatomy 0.000 title claims description 41
- 238000011156 evaluation Methods 0.000 title description 17
- 230000033001 locomotion Effects 0.000 claims description 219
- 210000003205 muscle Anatomy 0.000 claims description 91
- 230000000694 effects Effects 0.000 claims description 40
- 230000001133 acceleration Effects 0.000 claims description 37
- 238000000034 method Methods 0.000 claims description 33
- 230000001364 causal effect Effects 0.000 claims description 30
- 238000012545 processing Methods 0.000 claims description 6
- 230000000926 neurological effect Effects 0.000 claims description 3
- 210000003857 wrist joint Anatomy 0.000 description 42
- 238000004364 calculation method Methods 0.000 description 38
- 206010049816 Muscle tightness Diseases 0.000 description 37
- 230000006870 function Effects 0.000 description 25
- 238000000926 separation method Methods 0.000 description 13
- 208000012902 Nervous system disease Diseases 0.000 description 12
- 208000025966 Neurological disease Diseases 0.000 description 11
- 238000005259 measurement Methods 0.000 description 11
- 238000013500 data storage Methods 0.000 description 10
- 210000000707 wrist Anatomy 0.000 description 9
- 238000004458 analytical method Methods 0.000 description 8
- 210000000056 organ Anatomy 0.000 description 7
- 208000006011 Stroke Diseases 0.000 description 6
- 230000005540 biological transmission Effects 0.000 description 6
- 238000010586 diagram Methods 0.000 description 6
- 208000015879 Cerebellar disease Diseases 0.000 description 5
- 230000007659 motor function Effects 0.000 description 5
- 238000002474 experimental method Methods 0.000 description 4
- 210000001503 joint Anatomy 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 206010034010 Parkinsonism Diseases 0.000 description 3
- 230000002490 cerebral effect Effects 0.000 description 3
- 230000009017 pursuit movement Effects 0.000 description 3
- 206010003694 Atrophy Diseases 0.000 description 2
- 206010008025 Cerebellar ataxia Diseases 0.000 description 2
- 208000012661 Dyskinesia Diseases 0.000 description 2
- 208000018737 Parkinson disease Diseases 0.000 description 2
- 208000027089 Parkinsonian disease Diseases 0.000 description 2
- 230000009471 action Effects 0.000 description 2
- 230000036982 action potential Effects 0.000 description 2
- 230000037444 atrophy Effects 0.000 description 2
- 210000003169 central nervous system Anatomy 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000004590 computer program Methods 0.000 description 2
- 230000007850 degeneration Effects 0.000 description 2
- 238000003745 diagnosis Methods 0.000 description 2
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 230000003902 lesion Effects 0.000 description 2
- 210000000337 motor cortex Anatomy 0.000 description 2
- 210000002161 motor neuron Anatomy 0.000 description 2
- 238000010606 normalization Methods 0.000 description 2
- 230000001575 pathological effect Effects 0.000 description 2
- 238000004445 quantitative analysis Methods 0.000 description 2
- 208000011580 syndromic disease Diseases 0.000 description 2
- 208000024827 Alzheimer disease Diseases 0.000 description 1
- 208000009017 Athetosis Diseases 0.000 description 1
- 208000003174 Brain Neoplasms Diseases 0.000 description 1
- 206010008027 Cerebellar atrophy Diseases 0.000 description 1
- 206010012289 Dementia Diseases 0.000 description 1
- 206010067889 Dementia with Lewy bodies Diseases 0.000 description 1
- 208000014094 Dystonic disease Diseases 0.000 description 1
- 201000011240 Frontotemporal dementia Diseases 0.000 description 1
- 208000002972 Hepatolenticular Degeneration Diseases 0.000 description 1
- 208000023105 Huntington disease Diseases 0.000 description 1
- 208000009829 Lewy Body Disease Diseases 0.000 description 1
- 206010027146 Melanoderma Diseases 0.000 description 1
- 208000016285 Movement disease Diseases 0.000 description 1
- 208000001089 Multiple system atrophy Diseases 0.000 description 1
- 240000007817 Olea europaea Species 0.000 description 1
- 206010033799 Paralysis Diseases 0.000 description 1
- 208000000609 Pick Disease of the Brain Diseases 0.000 description 1
- 208000010112 Spinocerebellar Degenerations Diseases 0.000 description 1
- 208000018839 Wilson disease Diseases 0.000 description 1
- 210000000544 articulatio talocruralis Anatomy 0.000 description 1
- 210000004227 basal ganglia Anatomy 0.000 description 1
- 201000002922 basal ganglia calcification Diseases 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 208000016791 bilateral striopallidodentate calcinosis Diseases 0.000 description 1
- 210000001638 cerebellum Anatomy 0.000 description 1
- 230000000739 chaotic effect Effects 0.000 description 1
- GTKRFUAGOKINCA-UHFFFAOYSA-M chlorosilver;silver Chemical compound [Ag].[Ag]Cl GTKRFUAGOKINCA-UHFFFAOYSA-M 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000001054 cortical effect Effects 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 238000002059 diagnostic imaging Methods 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000035475 disorder Diseases 0.000 description 1
- 208000010118 dystonia Diseases 0.000 description 1
- 210000002310 elbow joint Anatomy 0.000 description 1
- 210000000245 forearm Anatomy 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 210000004394 hip joint Anatomy 0.000 description 1
- 230000001771 impaired effect Effects 0.000 description 1
- 238000002955 isolation Methods 0.000 description 1
- 210000000629 knee joint Anatomy 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 230000004807 localization Effects 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 201000006417 multiple sclerosis Diseases 0.000 description 1
- 230000004118 muscle contraction Effects 0.000 description 1
- 230000003387 muscular Effects 0.000 description 1
- 210000002346 musculoskeletal system Anatomy 0.000 description 1
- 210000001087 myotubule Anatomy 0.000 description 1
- 210000005036 nerve Anatomy 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000007170 pathology Effects 0.000 description 1
- 208000027232 peripheral nervous system disease Diseases 0.000 description 1
- 201000002212 progressive supranuclear palsy Diseases 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 210000000323 shoulder joint Anatomy 0.000 description 1
- 210000002435 tendon Anatomy 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/45—For evaluating or diagnosing the musculoskeletal system or teeth
- A61B5/4528—Joints
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/103—Detecting, measuring or recording devices for testing the shape, pattern, colour, size or movement of the body or parts thereof, for diagnostic purposes
- A61B5/11—Measuring movement of the entire body or parts thereof, e.g. head or hand tremor, mobility of a limb
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/103—Detecting, measuring or recording devices for testing the shape, pattern, colour, size or movement of the body or parts thereof, for diagnostic purposes
- A61B5/11—Measuring movement of the entire body or parts thereof, e.g. head or hand tremor, mobility of a limb
- A61B5/1121—Determining geometric values, e.g. centre of rotation or angular range of movement
-
- 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/316—Modalities, i.e. specific diagnostic methods
- A61B5/389—Electromyography [EMG]
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/40—Detecting, measuring or recording for evaluating the nervous system
- A61B5/4058—Detecting, measuring or recording for evaluating the nervous system for evaluating the central nervous system
- A61B5/4064—Evaluating the brain
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/68—Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient
- A61B5/6801—Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient specially adapted to be attached to or worn on the body surface
- A61B5/6813—Specially adapted to be attached to a specific body part
- A61B5/6824—Arm or wrist
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Medical Informatics (AREA)
- Surgery (AREA)
- Veterinary Medicine (AREA)
- Public Health (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Engineering & Computer Science (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- General Health & Medical Sciences (AREA)
- Molecular Biology (AREA)
- Animal Behavior & Ethology (AREA)
- Dentistry (AREA)
- Physiology (AREA)
- Oral & Maxillofacial Surgery (AREA)
- Neurology (AREA)
- Geometry (AREA)
- Psychology (AREA)
- Neurosurgery (AREA)
- Orthopedic Medicine & Surgery (AREA)
- Rheumatology (AREA)
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
- Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
Description
第一に、手関節運動に関わる二十数個の筋肉からわずか4個の主要筋を測定することで運動指令を解析することができた点である。
第二に、4個の筋活動を皮膚表面に貼り付けた電極で無痛かつ非侵襲的に記録することができた点である。
この結果、これまでの研究室レベルの運動指令の解析を簡便で非侵襲的に臨床現場で行えるようになった。
さらに第三に、以上の通り簡便で非侵襲的に記録した4個の筋活動と手首の動きの因果関係とを、関節トルクのレベルで同定する方程式を考案した点である(WO2009/028221)。
その結果、神経疾患における異常運動への個々の筋肉の関与を定量的に分析できるようになった。
(1)動く目標を追跡させる運動計測装置を用いて被験者に対し目標追跡運動を行なったときに測定された、関節の主動筋の筋電図情報並びに関節の位置、速度及び加速度情報から、脳内の運動制御機能を評価するシステムであって、
(a) 前記筋電図情報の周波数、並びに前記位置、速度及び加速度情報の周波数を複数の周波数成分に分離する手段、
(b) 前記筋電図情報並びに前記位置、速度及び加速度情報を下記運動方程式(1):
(c) 前記B/K比を指標として、脳内運動制御機能と目標追跡運動との因果関係を評価する手段
を含む前記システム。
(2)動く目標を追跡させる運動計測装置を用いて被験者に対し目標追跡運動を行なったときに測定された、関節の主動筋の筋電図情報並びに関節の位置、速度及び加速度情報から、脳内の運動制御機能を評価するためのプログラムであって、コンピュータに、
(a) 前記筋電図情報の周波数、並びに前記位置、速度及び加速度情報の周波数を複数の周波数成分に分離する手段、
(b) 前記筋電図情報並びに前記位置、速度及び加速度情報を下記運動方程式(1):
(c) 前記B/K比を指標として、脳内運動制御機能と目標追跡運動との因果関係を評価する手段
を実行させるための前記プログラム。
(3)上記(2)に記載のプログラムを記録したコンピュータ読み取り可能な記録媒体。
(4)動く目標を追跡させる運動計測装置を用いて被験者に対し目標追跡運動を行なったときに測定された、関節の主動筋の筋電図情報並びに関節の位置、速度及び加速度情報から、脳内の運動制御機能を評価する方法であって、
(a) 前記筋電図情報の周波数、並びに前記位置、速度及び加速度情報の周波数を複数の周波数成分に分離するステップ、
(b) 前記筋電図情報並びに前記位置、速度及び加速度情報を下記運動方程式(1):
(c) 前記B/K比を指標として、脳内運動制御機能と目標追跡運動との因果関係を評価するステップ
を含む前記方法。
(5)動く目標を追跡させる運動計測装置を用いて被験者に対し目標追跡運動を行なったときに測定された、関節の主動筋の筋電図情報並びに関節の位置、速度及び加速度情報を処理する方法であって、
(a) 前記筋電図情報の周波数、並びに前記位置、速度及び加速度情報の周波数を複数の周波数成分に分離するステップ、並びに
(b) 前記筋電図情報並びに前記位置、速度及び加速度情報を下記運動方程式(1):
を含む前記方法。
(6)本発明において、筋電図の大きさは、関節トルクの大きさに比例させて正規化し、当該正規化された筋電図を全波整流し、ローパスフィルターを通して求めることが好ましい。周波数成分としては、例えば低周波数成分及び高周波数成分の2つの周波数成分が挙げられる。また、2つの周波数成分に分けたときは、その境界は0.3〜0.8Hzの範囲内の周波数であることが好ましい。
本発明において、被験者としては例えば神経疾患患者である。
また、本発明において、測定の対象となる関節の主動筋は、例えば腕の筋である。腕の筋としては、例えば橈側手根伸筋(ECR)、尺側手根伸筋(ECU)、尺側手根屈筋(FCU)及び橈側手根屈筋(FCR))から選ばれる少なくとも1つを挙げることができる。
1.概要
本発明は、関節運動において関節を動かす所定個数の筋肉の筋活動と、その動きの3成分(位置、速度、加速度)との因果関係に基づき、運動中の筋活動から脳内の運動制御機能又は運動制御器の状態を推定する新規方法を提供するものである。
本発明の方法は、大きく2つのステップに分かれる。
(1)目標追跡運動中の筋活動と運動方程式との間の因果関係を関節トルクレベルで同定し、その際に求められる運動方程式のパラメータに、筋活動と動きの3成分の関係を集約して抽出する。
(2)そのパラメータに基づき、脳内の予測制御器とフィードバック制御器の機能的な評価を行う。
上記2つのステップを実現するために、本発明は、動く目標を追跡させる運動計測装置(マニピュランダム)を用いて被験者に対し目標追跡運動を行なったときに測定された、関節の主動筋の筋電図情報並びに関節の位置、速度及び加速度情報から、脳内の運動制御機能を評価するシステムを提供する。また、本発明は、このような手段をコンピュータに実行させるためのコンピュータプログラム、及びこのような手段を実行する機能評価方法を提供する。
本発明のシステム等においては、
(a) 関節の主動筋の筋電図情報の周波数、並びに関節の位置、速度及び加速度情報の周波数を複数の周波数成分に分離する手段又はステップ、
(b) 前記筋電図情報並びに前記位置、速度及び加速度情報を下記運動方程式(1):
(c) 前記B/K比を指標として、目標追跡運動と脳内運動制御器との因果関係を評価する手段又はステップが含まれる。
以前の出願(WO2009/28221)では、方程式の中の下記加速度項、速度項及び位置項の係数、即ちM、B及びKは他の研究者による先行技術(文献名:Gielen and Houk 1984; Grey 1997; de Serres and Milner 1991; Milner and Cloutier 1998)で得られた値をそのまま定数として使用していた。
この「考慮」のパターンは、理論的には少なくとも2つ想定される。
1つは、単純な関節の位置制御であり、他の1つは、関節の位置及び速度の両方を考慮するハイブリッド制御である。前者は工学的なフィードバック制御器に対応し、最も原始的な制御である。後者は工学的な予測制御器に対応し、より高級な制御である。
この理論的な着想を実際に実施させるために、本発明においては、第2の工夫として、「予測制御」及び「フィードバック制御」と呼ばれる制御機能の概念を考案し、脳から出力される運動指令が予測制御及びフィードバック制御の2つに分離されやすい運動課題をデザインした。そのために、目標が滑らかに等速で一定の軌道上を動き、その目標を関節の動きに連動するカーソルで追跡する運動課題を用いた。
例えば、図1に示すように、WO2009/28221に記載のシステムに従って運動機能を評価する場合を考える。WO2009/28221に記載のシステムにおいては、マニピュランダムに被験者の腕が固定され、手首だけが動かせるようにしてある。腕には、いくつかの筋肉の筋電図を測定するために電極が配置されている。PC画面上には、予め所定の文字や図形を描くように動作が設定されたターゲット画像(例えば「○」印で表示されるターゲット画像)と、被験者の手首関節の動きに対応して動くカーソルが表示される。そして、ターゲット画像(○印)の動きに合わせて、被験者がターゲット画像内にカーソルを収めるようにマニピュランダムで手首(手首関節)を動かす。例えば、図2に示すように、ターゲット画像が数字の「2」を描く動きに合わせて、被験者は手首関節を動かしてターゲットの○印内にカーソルを収めつつ、数字の「2」を描くように運動する(なお、図2では数字の「2」が破線で示されているが、実際の試行時には、破線は表示されない)。
このようにして所定の筋肉の筋電図を測定したとき、本発明者は、その測定された筋電図の波形から、所定の文字等を予測的に等速で描く動きは波形がゆるやかとなり、ターゲット画像内にカーソルを収めようと修正するする動きは波形が小刻みになる点に着目した。そして、前者は、(i) 目的の文字等を等速で描くときの運動に反映される周波数、後者は(ii) ターゲット画像内にカーソルを収めようとする動きに反映される周波数と考えて、それぞれの周波数に分離することを考えた。
被験者は、マニピュランダムによる試行の際に、数回の練習を行うことができるので、被験者は、実際の試行の際にカーソルの動きを予測することができる(例えば数字の「2」を描くように動かすことを予め知ることができる)。従って、上記(i)の周波数は、脳内において、動きの予測を制御する機能として反映され、上記(ii)の周波数は、正しい位置にカーソルを保持させて位置のずれに対する補正を行なおうとする機能として反映されるといえる。本発明では、この動きの予測を制御する機能を「予測制御」、正しい位置にカーソルを保持させようとする機能を「フィードバック制御」と定義し、それぞれの制御を担う脳内器官を「予測制御器」、「フィードバック制御器」という。
上記のことから、予測的運動指令は、目標の動きと同じ滑らかな低周波成分の運動指令として現れ、位置のずれに対する補正を行うフィードバック的運動指令は、細かな、より高周波の運動指令として現れると期待される。その結果、予測的運動指令とフィードバック的運動指令は周波数の違いで容易に分離でき、どちらの運動指令も、BとKとの比率により評価することができる。
実施例に示す通り、二つの運動制御器の出力は、予想通りの低周波の予測的運動指令と高周波のフィードバック的運動指令として分離され、本発明により、2つの運動制御器を分離して個別に評価する方法論の確立に至ったのである。
これに対し、本発明においては、その長年の仮説に対し、初めて実験的裏付けを与え、同時に、複雑な脳の運動制御系の状態を2つの仮想的制御器(「予測制御器」及び「フィードバック制御器」)に帰着させて評価する方法を確立したことになる。例えば、予測制御を担う脳内器官として、小脳などが優位に働き、フィードバック制御を担う脳内器官として、大脳皮質運動野などが優位に働く。従って、どちらの制御機能が優位又は劣位であるかによって、運動に携わる脳内器官の機能を解析することができる。
さらに、本発明のシステムでは、MRIやMEG等の高価な計測機器を使用せず、非侵襲的に神経疾患の患者にとっても負担が少なく、簡易にデータを得ることができ、ベッドサイドでの利用に適している。この点は実用上特筆すべき利点である。
前記マニピュランダムを用いた試行(実技)により得られたデータは、数字、文字、図形、記号等を追跡する際に手首が行なう動き成分と、その動きを行うための腕の筋肉の筋張力成分とに分けて解析することができる。
手首関節を動かす筋肉(主動筋)において、筋電図測定に使用する筋としては、特に限定されるものではないが、例えば橈側手根伸筋(ECR)、尺側手根伸筋(ECU)、尺側手根屈筋(FCU)及び橈側手根屈筋(FCR))であることが好ましく、これらの筋の1つ又は2つ以上を組合せて測定することができる。
目標追跡運動課題では、決まった軌道を等速で動くターゲットを追跡するが、健康成人であれば、ターゲットの動きを予測しながら追跡運動を行うことができる。従って、ターゲットの速度と一致する低周波領域の運動指令には、既知の軌道と速度で動いているターゲットを予測して追跡する予測制御の成分が主に入っていると解釈できる。一方、高周波の運動指令はターゲットの動きとの相関はなく、予測的な制御とは無関係である。
上記運動方程式からBとKの意味を分析してみると、粘性係数Bが高いことは、筋活動の線形和と速度成分が高い相関を持つことを意味し、筋活動の中に速度制御の運動指令が多く含まれることを示す。一方、弾性係数Kが大きいことは、筋活動の線形和と位置成分が高い相関を持つことを意味し、筋活動の中に位置制御の運動指令が多く入っていることを示す。そこで2つ成分のB/K比の違いを分析してみると、運動成分の低周波領域についてはBとKがほぼ同じ比率であり(図7B)、筋活動が位置成分と速度成分の両方と深い関係がある。他方、運動成分の高周波領域ではKが優位であり(図8B)、筋活動が手首の位置成分のみと深い関係があることを示す。
以上をまとめてみると、低周波領域の運動指令は、既知の軌道と速度で動くターゲットの位置及び速度の両方にマッチする関節の動きをリアルタイムで生成し、課題の要求に応えるための最も重要な運動指令と考えられる。一方、高周波領域の運動指令は、カーソルの位置がターゲットの中心から偏倚しかけた時に素早くオンラインで補正する、フィードバック的又は補足的な役割を担うといえる。
ここで指標追跡運動における2つの運動指令成分を、それを生成する運動制御器という観点から考えてみる。低周波領域の運動指令は位置と速度成分の両方を一定の割合で含み、位置と速度両方を明示的に指定している。これに対して、高周波領域の運動指令はほぼ位置成分のみであり、速度指定はない。以上のように、2つの運動指令は位置成分と速度成分のブレンド比が全く異なり、制御のシステムが根本的に違う。さらに、両者の運動指令の間には全く相関がないことから、それぞれを生成する運動制御器は独立して並列に働いていることが強く示される。
以上の考察に基づいて目標追跡運動の運動制御メカニズムを考えると、図9のようなモデルが考えられる。つまり、目標追跡運動には、既知の軌道と速度で動くターゲットの動きを予測し、それと同じ関節の動きを位置、速度の両方を指定して再現する予測的制御器に加えて、ターゲットとの関節(カーソル)位置の誤差のみを修正するフィードバック制御器という、2つの並列制御器が存在する。二つの制御器からの出力は運動ニューロンまでの中枢神経系のどこかで加算され、最終的に同じ筋肉を駆動することになる。
従って、各周波数成分についてB/K比を算出し、さらに正常パターンと比較することで、どのような疾患の患者でどちらの制御器がどのように障害されるかが分かり、医師が患者の病態を理解し、今後の治療方針やスケジュールを設定する上で重要な情報となる。
ここで、本発明のシステムに適用可能な被験者は、健常人のほか、神経系疾患患者が挙げられる。本発明のシステムは、例えば、神経疾患を治療する際に使用される。治療には、例えば、神経疾患患者における運動機能のリハビリテーションが含まれる。また前記神経疾患は、例えば、運動障害を伴う神経疾患が挙げられ、具体的には、パーキンソン病、パーキンソン症候群、ハンチントン病、アテトーゼ、ジストニー、小脳・脊髄萎縮症(小脳疾患、脊髄小脳変性症を含む)、多系統萎縮症、線状体黒質変性症、オリーブ橋小脳萎縮症、シャイ・ドレーガー症候群、皮質基底核変性症、進行性核上性麻痺、大脳基底核石灰化症、パーキンソニズム痴呆症候群、びまん性レヴィ小体症、アルツハイマー病、ピック病、ウィルソン病、多発性硬化症、末梢神経疾患、脳腫瘍及び脳卒中からなる群より選ばれる少なくとも1種が挙げられ、なかでも特に、パーキンソン病、パーキンソン症候群、小脳・脊髄萎縮症、脳卒中が好ましい。
図12は、本発明のシステムの構成図である。図12において、本発明のシステムは計算手段10とデータベース20から構成され、計算手段10は、(i) 測定データ蓄積手段11、(ii) 関節トルク及び筋張力計算手段12、(iii) 因果関係計算手段13、及び(iv)計算結果表示手段14を備える。
(i) 測定データ蓄積手段11
測定データ蓄積手段11は、マニピュランダムにより計測された筋電図並びに関節位置及び動きデータを蓄積するための手段であり、蓄積された情報は、グラフ又は表により確認することができる。
(ii) 周波数分離、並びに関節トルク及び筋張力計算手段12
周波数分離並びに関節トルク及び筋張力計算手段12は、測定データ蓄積手段11又はデータベース20から、周波数分離を実行するとともに、各分離された周波数成分について、運動方程式を用いて関節トルク及び筋張力を計算する手段である。この手段では、マニピュランダムからの情報に従い、動き成分である関節の位置、速度、加速度を計算することができる。
関節トルクの計算式(2)は、以下の通りである。
因果関係計算手段13は、分離された周波数成分、B/K比の計算及び試行回数から、筋活動と動き成分との因果関係、目標追跡運動における手首の速度成分と筋電図から得られた周波数成分との関係、あるいは予測制御とフィードバック制御の割合などを計算する手段である。
計算結果出力手段14は、測定データ、計算された前記関節トルク、動き成分及び筋張力、B/K比、並びに因果関係を出力する手段であり、上記(ii)により計算された関節トルク、筋張力をアニメーション又はグラフ表示する。また、これと同時に(iii)により計算された因果関係やB/K比をグラフ表示することができる。
入力された測定データと計算結果は関連付けられてデータ蓄積手段としてデータベース20に保存される。
保存された計算条件と計算結果は、再度データベース20から、あるいは蓄積手段11、計算手段12及び13、並びに計算結果表示手段14から読み込むことができる。
本発明のプログラムにおいて、コンピュータを実行させるための手段を示す構成例を図13に示す。
図13は、本発明のプログラムを実行させるためのシステム100の詳細構成図である。図13において、システム100は、図12に示す計算部10及びデータベース(以下「DB」という)20を備え、さらに、制御部101、送信/受信部102、入力部103、出力部104、ROM105、RAM106、ハードディスクドライブ(HDD)107、CD-ROMドライブ108により構成される。
図14は、本発明のプログラムの動作を説明するフローチャートである。
本発明のプログラムは、入力された測定データをもとに、
(a) 周波数分離、並びに関節トルク及び筋張力計算手段12による計算、並びに
(b) 因果関係計算手段13による計算を行う。
計算の起動は制御部101からの指令に基づいて行うことができる。
計算は、例えば以下の順に実施される。
(i)測定データ蓄積手段11又はデータベース20から手関節位置データを読み取る(S201)。
(ii)周波数分離、並びに関節トルク及び筋張力計算手段12で、手関節位置データを周波数分離する(S202)。
(iii)周波数分離、並びに関節トルク及び筋張力計算手段12で、位置データの各周波数成分について、動き成分である手関節の位置、速度、加速度を計算する(S203)。
(iv)周波数分離、並びに関節トルク及び筋張力計算手段12で、位置データの各周波数成分について、動き成分(手関節の位置、速度、加速度)から運動方程式を用いて関節トルクを計算する(S204)。
(v)他方、測定データ蓄積手段11又はデータベース20から筋電図信号データを読み取る(S205)。
(vi)周波数分離、並びに関節トルク及び筋張力計算手段12で、筋電図信号から筋張力を求める(S206)。(図4参照)
(vii)周波数分離、並びに関節トルク及び筋張力計算手段12で、筋張力を周波数分離する(S207)。
(viii)周波数分離、並びに関節トルク及び筋張力計算手段12で、各周波数成分について、計算された関節トルクと筋張力の線形和との間で最適近似を行う(S208)。(図5参照)
(ix)因果関係計算手段13で、各周波数成分について、様々なBとKにおいて手順(iii)〜手順(viii)による筋活動と動き成分の因果関係の同定を行い、最適なBとKの比率(B/K比)を計算する(S209)。
(x)計算結果出力手段14で、計算された関節トルク、動き成分及び筋張力、B/K比、並びに因果関係を出力する(S210)。表示結果の一例を図4、図6、図7及び図8に示す。
(xi)計算結果は、データ蓄積手段のデータベース20に逐次保存される(S210)。
計算結果は、データ蓄積手段のデータベースに逐次保存される。
本発明のプログラムは、例えばC言語、Java(登録商標)、Perl、Fortran、Pascal等で書くことができ、そしてクロスプラットフォームに対応できるように設計されている。従って、このソフトウエアはWindows(登録商標)95/98/2000/XP/Vista/7、Linux、UNIX(登録商標)、Macintosh等で作動させることが可能である。
本発明は、被験者の様々な関節運動における運動指令としての筋活動と、その結果である関節の動きの因果関係に基づいて解析する方法であり、手関節運動を利用した運動指令解析評価診断システム(WO2009/028221)を用いた。図1に、本発明で使った運動指令解析評価診断システムの概要図と実験の様子を示した。
被験者はカーソルとターゲットが表示されているPC画面の前に座り、前腕を支持台に載せ、手で手関節マニピュランダムを操作する。手関節の2自由度の動きは、2個の位置センサーで計測され、PC画面上のカーソル(直径2ミリの黒点)の動きに反映される。ターゲットは丸い円で表示され、その直径(1cm)は手関節の4.5度の動きに相当し、その位置は手関節の動きの指標となる。
実験には神経疾患の既往のない健康成人4名(44-63歳)が被験者として参加した。被験者には等速(平均速度:6.2deg/sec)で動くターゲットを追跡する手関節運動を実験課題として行わせた(図2)。運動の開始点はモニターの左上(X=-10°,Y=8°)である。最初に、モニターの左上に円形のターゲットが表示されると、被験者は手関節を動かして連動するカーソルを開始点に保持する。3秒後にターゲットが等速で数字2の軌道を描きながら動き始める。被験者は手関節を動かして動いているターゲットの中に可能な限りカーソルを保ち続ける。各被験者はこのタスクを十分に理解できるように2〜3回練習した後、5回ずつの本番の試行を行った。
運動課題を行っている際に、手関節の位置(X,Y)と手関節運動に関わる4個の主動筋(橈側手根伸筋(ECR)、尺側手根伸筋(ECU)、尺側手根屈筋(FCU)、橈側手根屈筋(FCR))からの表面筋電図信号を同時記録した。手関節の速度と加速度は位置信号をそれぞれ1階、2階微分して求め、手関節の関節トルクを推定するのに使用した。
筋電図の記録は、一対の銀塩化銀表面電極を用いて皮膚表面電位を双極誘導し、 2kHz,12bit でサンプリングした。電極の直径は5mmとし、電極間の距離が10mm となるように筋線維に沿って貼った。測定した筋肉と電極のおよその位置を図3Aに示す。表面電極による筋電図記録では、記録されている筋肉の同定が重要であり、実験精度を左右する。本実験系における同定の信頼性を評価するために、1名の被験者において4個の筋肉について表面筋電図(図3C左)とワイヤー電極による表面電極直下の筋肉からの直接記録(図3C右)の同時記録を行い比較した。
その結果、両者の波形は非常によく一致し、それに基づいて計算された筋活動のピークの方向(preferred direction, PD)もほぼ一致した(図3B)。このようにして、筋肉の同定と表面筋電図による記録の有効性を確認した。
中枢神経系からの運動指令は脊髄運動ニューロンから筋肉に活動電位を誘発し、筋肉を収縮させる。この筋収縮の活動電位を皮膚上で計測したものが表面筋電図(以下筋電図)であり、その大きさは筋張力と比例関係にある(Basmajian and De Luca 1985)。しかし、計測された筋電図は皮膚抵抗や筋肉に対する電極の相対的位置により記録される信号の絶対値が異なるため、そのままでは定量的な分析には適さない。そこで図4に示すように、筋電図の大きさを関節トルクの大きさに比例させて調整する正規化を行った。すなわち、一定の力を発生している筋活動の大きさを、異なる被験者間、又は同一被験者の異なる日時の記録の間で一定とするように調整した。次いで、正規化された筋電図を全波整流(記録された筋電図の絶対値計算)した後、最終的に2次のローパスフィルター(遮断周波数: 3 Hz)を通して筋張力を求めた。「ローパスフィルター」(law-pass filter)とは、全波整流された筋電図波形を筋張力に変換する畳込み積分を意味する。
本発明では図5に示したように手関節運動における筋活動と手関節の動き成分の間の因果関係の同定を関節トルクレベルで行い、その際に求まる運動方程式の粘性係数Bと弾性係数Kの大小から、原因の筋活動が動きのどの成分と関係が深いかを調べた。
具体的な計算としては、式(1)の中辺の運動方程式により関節トルク(図5中央の青線)を計算し、次いでそれを4個の筋活動(筋張力)の線形和(式(1)の右辺)(図5中央の赤線)で最適近似した。これは関節トルクを媒介にして筋活動と運動方程式の間の因果関係を同定したことになる。
手関節の運動中のBとKの値は実測例がないため、安静時のBとKの値(Gielen and Houk 1984; Grey 1997; de Serres and Milner 1991; Milner and Cloutier 1998)に基づき範囲を生理的な範囲に制限して[B: 0-0.5Nms/rad; K: 0-0.4Nm/rad]最適値を求めた。特に、手関節の場合、慣性モーメントM(定数)がB,Kに比べて非常に小さく、加速度項が事実上無視できる。そのため未定係数B,Kの絶対値を特定することは困難であるので、本発明ではBとKの値をその比率(B/K比)として評価することにした。その結果、BとKの絶対値を求める技術的に困難な問題を回避することが可能となった。この点も技術的に重要な工夫である。
目標追跡運動中の手関節(カーソル)の動きをみると、ターゲットを追随する滑らかな動きに加えて、小刻みに振動するような、より高周波の動きが見えた。そこで、この現象を目標追跡運動の速度(X成分, Y成分)と接線速度から分析した(図6A)。図6Aを見ると、手関節の速度変化(青線)はターゲットの緩やかな速度変化(緑線)を中心に、小刻みに振動している。この振動現象は全ての被験者に共通して見られ、全被験者に対して速度成分の周波数分析を行った(図6B)。
その結果、およそ0.5Hzを境界に、より低い周波数領域と高い周波数領域に分離する傾向が見られた。そこで速度成分を低周波(0〜0.5Hz)の「F1領域」と高周波(0.5〜3Hz)の「F2領域」に分離した(図6C)。
低周波のF1領域(図6C(a))では目標の速度(緑線)とぴったり一致する手関節の速度成分(青線)だけが残り、高周波のF2領域(図6C(b))では主に1.4Hz程度の周波数で振動する目標の動きとは無相関な速度成分が残った。
この結果は次の通り解釈することができる。すなわち、目標追跡運動課題では決まった軌道を等速で動くターゲットを追跡する。健康成人であれば2〜3回の練習でターゲットの動きを完全に理解でき、ターゲットの動きを予測しながらほぼ100%近い成功率で追跡運動が行える。従って、ターゲットの速度と一致するF1領域の運動指令には、既知の軌道と速度で動いているターゲットを予測して追跡する予測制御の成分が主に入っていると解釈できる。一方、F2領域の高周波の運動指令はターゲットの動きとの相関はなく、予測的な制御とは無関係と考えられる。
図6に示されたように、目標追跡運動における2つの成分は明らかに機能的に異なる役割を果たす。そこで、2つの成分のそれぞれについて、運動指令と手関節の動きの間に最適な近似を与えるB/K比を計算し、その結果から各成分の機能的な意義を検討した。
まず4個の筋張力とキネマティクスの各成分をF1領域(0〜0.5Hz)とF2領域(0.5〜3Hz)に分離し、運動方程式(式(1))を用いて先述の筋活動と手関節の動きの間の最適近似を行った。
F1領域とF2領域に対する最適近似の例を、それぞれ図7Aと図8Aに示す。
どちらの成分についても(図7Aと図8Aの最下段)、筋張力により手関節のトルクが高い相関で近似できていることがわかる(F1領域:R = 0.98, F2領域:R = 0.70)。F2領域の相関が若干低いが、健常被験者ではこの領域の筋活動がF1領域に比べて小さく、S/N比が劣るためと推定される。驚くべきことに、2つの運動指令成分におけるB/K比は好対照を示した。低周波のF1領域ではB/K比が高く(図7B,C)、対照的に高周波のF2領域のB/K比は非常に低い(図8B,C)。
運動方程式からBとKの意味を分析してみると、粘性係数Bが高いことは、筋活動の線形和と速度成分が高い相関を持つことを意味し、筋活動の中に速度制御の運動指令が多く含まれることを示す。一方、弾性係数Kが大きいことは、筋活動の線形和と位置成分が高い相関を持つことを意味し、筋活動の中に位置制御の運動指令が多く入っていることを示す。この解釈に基づき2つ成分のB/K比の違いを分析してみると、BとKがほぼ同じ比率であるF1領域では、筋活動が位置成分と速度成分の両方と深い関係があり、Kが優位なF2領域では、筋活動が手首の位置成分のみと深い関係があることを示す。
1.2つの並列運動制御器の臨床応用1 - 小脳疾患
本発明の方法は、神経疾患の混沌とした病態を、2つの運動制御器という機能的意味が明瞭なモジュールに帰着させることにより、直観的な解釈が容易な運動機能の評価ができる点に最大の利用価値がある。
例えば、健常者では、小脳の予測制御器を十分に使って運動を行えるため、フィードバック制御器は補足的に使うと予想される(図10A)が、小脳疾患では、小脳の予測制御器の障害を代償するために、運動野のフィードバック制御器を多用すると予想される(図10B)。
この予想を実際のデータで裏付けたのが図10Cである。
この図では、赤で示す予測制御の運動指令と青で示すフィードバック制御の運動指令がどのような割合で含まれるかを正常者と小脳疾患患者で比較したものである。
予想通りに、健常被験者は主に予測制御器に依存して運動を行っているが、小脳疾患では障害の重さに逆相関して予測制御の割合が減少し、代わりにフィードバック制御を多く使うことが分かった。
このように、本発明の方法を用いると、神経疾患が脳の運動制御系にどの様な影響(定性的)をどの程度(定量的)与えているかを、予測制御とフィードバック制御という、二つの運動制御器の観点から明快に分析できる。
本発明の方法を応用し、脳卒中患者の複雑な脳内運動制御系の状態を、予測制御器とフィードバック制御器という2つの仮想的運動制御器の指標として定量化し、リハビリ可能性を評価し、また、麻痺からの回復過程を記述した。
図11Aには、ある脳卒中患者が指標追跡運動を5回繰り返し行ったときの、1回目と5回目の手関節の軌跡を示す。1回目は千々に乱れた動きの連続で意図した動きが全くできていないが、5回目には明らかな改善が見られ、数字の2の形が現れた。この軌道の改善に対応して、回を重ねるにつれて予測制御の運動指令の割合が増加し(図11C)、同時に同領域の運動指令のB/K比が上昇し、滑らかな速度制御を行う患者の意図が運動指令に現れてきていることがわかる(図11B)。この結果は、少なくとも一部の片麻痺患者の運動制御器が柔軟な学習能力を保持していることを示すものである。
本発明の方法により、脳内の2つの運動制御器の状態を、脳の信号を記録することなく簡便・非侵襲的に定量評価できるようになった。このような機能的評価法は従来のリハビリテーションの評価や、神経学的診断・評価法の範疇には存在しなかった。
近年、MRI等による画像診断技術の進歩により、脳内の病変の位置と広がりに関しては、ミリ単位で、小型の神経核程度の精度で評価することは今や常識である。画像的な局在診断は一つの到達点にある。ところが、病変の機能的意義を評価する方法論は、誇張でなくこの100年間大きな進歩がなかった。近年開発が進んだ運動の定量的分析法は、運動の記述を定量化しただけであり、脳内の運動制御器に関する評価は全くできなかった。本発明の方法は、脳内の運動制御器を機能的かつ定量的に評価できる唯一のシステムである。
101:制御部、102:送信/受信部、103:入力部、104:出力部、
105:ROM、106:RAM、107:ハードディスクドライブ、108:CD-ROMドライブ
111:インターネット回線、120:CD-ROM
Claims (11)
- 動く目標を追跡させる運動計測装置を用いて被験者に対し目標追跡運動を行なったときに測定された、関節の主動筋の筋電図情報並びに関節の位置、速度及び加速度情報から、脳内の運動制御機能を評価するシステムであって、
(a) 前記筋電図情報の周波数、並びに前記位置、速度及び加速度情報の周波数を複数の周波数成分に分離する手段、
(b) 前記筋電図情報並びに前記位置、速度及び加速度情報を下記運動方程式(1):
(c) 前記B/K比を指標として、目標追跡運動と脳内運動制御機能との因果関係を評価する手段
を含む前記システム。 - 筋電図の大きさは、関節トルクの大きさに比例させて正規化し、当該正規化された筋電図を全波整流し、ローパスフィルターを通して求めるものである請求項1に記載のシステム。
- 周波数成分が、低周波数成分及び高周波数成分の2つの周波数成分である請求項1に記載のシステム。
- 2つの周波数成分の境界が0.3〜0.8Hzの範囲内の周波数である請求項3に記載のシステム。
- 被験者が神経疾患患者である請求項1に記載のシステム。
- 関節の主動筋が腕の筋である請求項1に記載のシステム。
- 腕の筋が、橈側手根伸筋(ECR)、尺側手根伸筋(ECU)、尺側手根屈筋(FCU)及び橈側手根屈筋(FCR))から選ばれる少なくとも1つである請求項6に記載のシステム。
- 動く目標を追跡させる運動計測装置を用いて被験者に対し目標追跡運動を行なったときに測定された、関節の主動筋の筋電図情報並びに関節の位置、速度及び加速度情報から、脳内の運動制御機能を評価するためのプログラムであって、コンピュータに、
(a) 前記筋電図情報の周波数、並びに前記位置、速度及び加速度情報の周波数を複数の周波数成分に分離する手段、
(b) 前記筋電図情報並びに前記位置、速度及び加速度情報を下記運動方程式(1):
(c) 前記B/K比を指標として、目標追跡運動と脳内運動制御機能との因果関係を評価する手段
を実行させるための前記プログラム。 - 請求項8に記載のプログラムを記録したコンピュータ読み取り可能な記録媒体。
- 動く目標を追跡させる運動計測装置を用いて被験者に対し目標追跡運動を行なったときに測定された、関節の主動筋の筋電図情報並びに関節の位置、速度及び加速度情報から、脳内の運動制御機能を評価する方法であって、
(a) 前記筋電図情報の周波数、並びに前記位置、速度及び加速度情報の周波数を複数の周波数成分に分離するステップ、
(b) 前記筋電図情報並びに前記位置、速度及び加速度情報を下記運動方程式(1):
(c) 目標追跡運動と脳内の運動制御機能との因果関係を評価する手段が、前記B/K比を指標として、粘性係数Bが大きい場合は筋活動の中に速度制御の運動指令が相対的に多く含まれると判定し、弾性係数Kが大きい場合は筋活動の中に位置制御の運動指令が相対的に多く含まれると判定するステップ
を含む前記方法。
- 動く目標を追跡させる運動計測装置を用いて被験者に対し目標追跡運動を行なったときに測定された、関節の主動筋の筋電図情報並びに関節の位置、速度及び加速度情報を処理する方法であって、
(a) 前記筋電図情報の周波数、並びに前記位置、速度及び加速度情報の周波数を複数の周波数成分に分離するステップ、並びに
(b) 前記筋電図情報並びに前記位置、速度及び加速度情報を下記運動方程式(1):
を含む前記方法。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2010042301A JP5623759B2 (ja) | 2010-02-26 | 2010-02-26 | 筋電図信号に基づいた脳内の並列運動制御機能の同定及び評価法 |
US12/807,861 US8660641B2 (en) | 2010-02-26 | 2010-09-15 | Method for identifying and evaluating parallel motor control function in the brain based on electromyogram signals |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2010042301A JP5623759B2 (ja) | 2010-02-26 | 2010-02-26 | 筋電図信号に基づいた脳内の並列運動制御機能の同定及び評価法 |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2011177228A JP2011177228A (ja) | 2011-09-15 |
JP5623759B2 true JP5623759B2 (ja) | 2014-11-12 |
Family
ID=44505657
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2010042301A Expired - Fee Related JP5623759B2 (ja) | 2010-02-26 | 2010-02-26 | 筋電図信号に基づいた脳内の並列運動制御機能の同定及び評価法 |
Country Status (2)
Country | Link |
---|---|
US (1) | US8660641B2 (ja) |
JP (1) | JP5623759B2 (ja) |
Families Citing this family (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9936891B2 (en) * | 2011-06-08 | 2018-04-10 | Precision Biometrics, Inc. | Systems and methods for providing biometric related to performance of a physical movement |
US10327675B2 (en) * | 2013-04-30 | 2019-06-25 | Tokyo Metropolitan Institute Of Medical Science | Motor function analysis system and operational method of system |
JP5898134B2 (ja) * | 2013-06-27 | 2016-04-06 | トヨタ自動車株式会社 | リハビリテーション装置、制御方法及び制御プログラム |
US9265685B1 (en) * | 2014-05-01 | 2016-02-23 | University Of South Florida | Compliant bimanual rehabilitation device and method of use thereof |
JP5836429B2 (ja) * | 2014-05-23 | 2015-12-24 | 日本電信電話株式会社 | 身体状態提示装置、方法及びプログラム |
CN104492033B (zh) * | 2014-12-17 | 2017-07-21 | 中国科学院自动化研究所 | 基于sEMG的单关节主动训练控制方法及相应的康复机器人 |
CA2880623A1 (en) * | 2015-02-02 | 2016-08-02 | Jian Chen | Emg diagnostic method for traumatic brain injury |
EP3315177A4 (en) * | 2016-06-30 | 2019-06-12 | Shanghai Fourier Intelligence Co., Ltd. | SUPERIOR MEMBER REHABILITATION DRIVE MACHINE |
KR20180031176A (ko) * | 2016-09-19 | 2018-03-28 | 삼성전자주식회사 | 근육의 특성에 대한 파라미터 식별 방법, 이를 이용한 보행 보조 방법, 및 이를 수행하는 장치들 |
JP2019162288A (ja) * | 2018-03-20 | 2019-09-26 | 公立大学法人公立諏訪東京理科大学 | 判定値算出装置、方法、及びプログラム、マイクロエラー発生判断装置、方法、及びプログラム |
CN110584790B (zh) * | 2019-08-13 | 2022-05-06 | 广东工业大学 | 基于手臂刚度的手术机器人遥操作比例控制方法 |
KR102322105B1 (ko) * | 2019-11-18 | 2021-11-05 | 주식회사 씨앰지 | 동작인식분석 기반의 뇌 인지기능 평가 방법 및 평가 시스템 |
Family Cites Families (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
AU4744796A (en) | 1994-12-29 | 1996-07-24 | Burgan, Pinny | System and method for detection and quantification of parkinson's disease |
JP2000279463A (ja) | 1999-03-31 | 2000-10-10 | Sanyo Electric Co Ltd | 上肢機能回復訓練装置 |
JP4524054B2 (ja) | 2001-06-14 | 2010-08-11 | 株式会社ユニメック | 脳・神経系疾患の病状診断、投薬処方及び機能回復訓練のための支援システム |
JP4045487B2 (ja) * | 2002-06-06 | 2008-02-13 | 日本電信電話株式会社 | 筋骨格系機械インピーダンス計測方法及び装置 |
JP3660330B2 (ja) * | 2002-08-13 | 2005-06-15 | 独立行政法人科学技術振興機構 | インピーダンス測定装置及び運動学習支援装置 |
JP4078419B2 (ja) * | 2003-06-18 | 2008-04-23 | 独立行政法人産業技術総合研究所 | ポインティングデバイス |
JP2005185557A (ja) | 2003-12-25 | 2005-07-14 | Yunimekku:Kk | 運動能力評価システム |
JP2007202612A (ja) * | 2006-01-31 | 2007-08-16 | Hiroshi Kato | 身体動作能力の質を定量化するシステム装置 |
JP4702860B2 (ja) * | 2007-07-18 | 2011-06-15 | キリンホールディングス株式会社 | 嚥下感覚の簡易測定方法 |
US8792977B2 (en) * | 2007-08-31 | 2014-07-29 | Tokyo Metropolitan Institute Of Medical Science | Quantitative motor function evaluation system |
-
2010
- 2010-02-26 JP JP2010042301A patent/JP5623759B2/ja not_active Expired - Fee Related
- 2010-09-15 US US12/807,861 patent/US8660641B2/en not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
US20110213267A1 (en) | 2011-09-01 |
JP2011177228A (ja) | 2011-09-15 |
US8660641B2 (en) | 2014-02-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP5623759B2 (ja) | 筋電図信号に基づいた脳内の並列運動制御機能の同定及び評価法 | |
JP5154558B2 (ja) | 定量的運動機能評価システム及び運動機能評価用プログラム | |
Gohel et al. | Review on electromyography signal acquisition and processing | |
Jerbi et al. | Inferring hand movement kinematics from MEG, EEG and intracranial EEG: From brain-machine interfaces to motor rehabilitation | |
JP5784709B2 (ja) | 神経障害を診断および/または治療する際に使う装置 | |
JP2020034849A (ja) | 作業支援装置、作業支援方法及び作業支援プログラム | |
KR20190077305A (ko) | 적응적 응답 기한 절차들에서의 신호 검출 메트릭들을 구현하기 위한 플랫폼들 | |
CN107961135A (zh) | 康复训练*** | |
Memar et al. | Objective assessment of human workload in physical human-robot cooperation using brain monitoring | |
CN112955066A (zh) | 治疗空间评估 | |
Lu et al. | Evaluating rehabilitation progress using motion features identified by machine learning | |
Lambercy et al. | Robots for measurement/clinical assessment | |
Garro et al. | Neuromechanical biomarkers for robotic neurorehabilitation | |
JP7070253B2 (ja) | パフォーマンス計測装置、パフォーマンス計測方法及びパフォーマンス計測プログラム | |
CN110931104A (zh) | 基于机器学习的上肢康复机器人智能训练***及方法 | |
Scano et al. | Robotic assistance for upper limbs may induce slight changes in motor modules compared with free movements in stroke survivors: a cluster-based muscle synergy analysis | |
CN115054250B (zh) | 镜像溢出运动检测分析方法和*** | |
Pogrzeba et al. | Analysis and quantification of repetitive motion in long-term rehabilitation | |
Morishita et al. | Brain-machine interface to control a prosthetic arm with monkey ECoGs during periodic movements | |
Chen et al. | Association between finger-to-nose kinematics and upper extremity motor function in subacute stroke: a principal component analysis | |
CN113614751A (zh) | 脑电图信号的识别和提取 | |
Hajian et al. | An investigation of dimensionality reduction techniques for EMG-based force estimation | |
Hajian et al. | Multimodal estimation of endpoint force during quasi-dynamic and dynamic muscle contractions using deep learning | |
Gunawardane et al. | A fusion algorithm for saccade eye movement enhancement with EOG and lumped-element models | |
Hajian et al. | PCA-based channel selection in high-density EMG for improving force estimation |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A711 | Notification of change in applicant |
Free format text: JAPANESE INTERMEDIATE CODE: A711 Effective date: 20121015 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A821 Effective date: 20121015 |
|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20130215 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20131118 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20131203 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20131227 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20140218 |
|
RD03 | Notification of appointment of power of attorney |
Free format text: JAPANESE INTERMEDIATE CODE: A7423 Effective date: 20140220 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20140310 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20140228 |
|
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: 20140902 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20140925 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 5623759 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
LAPS | Cancellation because of no payment of annual fees |