JP2009095552A - 心拍変動解析方法および解析装置 - Google Patents
心拍変動解析方法および解析装置 Download PDFInfo
- Publication number
- JP2009095552A JP2009095552A JP2007271434A JP2007271434A JP2009095552A JP 2009095552 A JP2009095552 A JP 2009095552A JP 2007271434 A JP2007271434 A JP 2007271434A JP 2007271434 A JP2007271434 A JP 2007271434A JP 2009095552 A JP2009095552 A JP 2009095552A
- Authority
- JP
- Japan
- Prior art keywords
- heart rate
- action potential
- rate variability
- data
- fluctuation
- 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.)
- Granted
Links
- 238000004458 analytical method Methods 0.000 title claims abstract description 32
- 230000001133 acceleration Effects 0.000 claims abstract description 31
- 230000002107 myocardial effect Effects 0.000 claims description 50
- 230000033764 rhythmic process Effects 0.000 claims description 29
- 210000000038 chest Anatomy 0.000 claims description 6
- 238000000034 method Methods 0.000 claims description 6
- 230000002123 temporal effect Effects 0.000 claims 1
- 238000005259 measurement Methods 0.000 abstract description 11
- 230000000241 respiratory effect Effects 0.000 abstract description 9
- 230000000694 effects Effects 0.000 abstract description 5
- 210000000115 thoracic cavity Anatomy 0.000 abstract description 4
- 230000036982 action potential Effects 0.000 abstract 3
- 230000000747 cardiac effect Effects 0.000 abstract 3
- 238000012544 monitoring process Methods 0.000 abstract 1
- 230000007704 transition Effects 0.000 description 6
- 230000006870 function Effects 0.000 description 5
- 238000004891 communication Methods 0.000 description 4
- 238000012790 confirmation Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 201000010099 disease Diseases 0.000 description 4
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 4
- 230000008035 nerve activity Effects 0.000 description 4
- 238000005070 sampling Methods 0.000 description 4
- 206010003119 arrhythmia Diseases 0.000 description 3
- 230000006793 arrhythmia Effects 0.000 description 3
- 230000002159 abnormal effect Effects 0.000 description 2
- 230000002567 autonomic effect Effects 0.000 description 2
- 239000000919 ceramic Substances 0.000 description 2
- 230000004217 heart function Effects 0.000 description 2
- 210000004072 lung Anatomy 0.000 description 2
- 210000005037 parasympathetic nerve Anatomy 0.000 description 2
- 201000002859 sleep apnea Diseases 0.000 description 2
- 238000001228 spectrum Methods 0.000 description 2
- 230000002889 sympathetic effect Effects 0.000 description 2
- 208000024172 Cardiovascular disease Diseases 0.000 description 1
- 208000031226 Hyperlipidaemia Diseases 0.000 description 1
- 206010020772 Hypertension Diseases 0.000 description 1
- 208000008589 Obesity Diseases 0.000 description 1
- 206010042434 Sudden death Diseases 0.000 description 1
- 208000033774 Ventricular Remodeling Diseases 0.000 description 1
- 210000001015 abdomen Anatomy 0.000 description 1
- 230000005856 abnormality Effects 0.000 description 1
- 210000003403 autonomic nervous system Anatomy 0.000 description 1
- 238000003287 bathing Methods 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000002612 cardiopulmonary effect Effects 0.000 description 1
- 230000009084 cardiovascular function Effects 0.000 description 1
- 206010008118 cerebral infarction Diseases 0.000 description 1
- 208000026106 cerebrovascular disease Diseases 0.000 description 1
- 230000000739 chaotic effect Effects 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 206010012601 diabetes mellitus Diseases 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 235000020824 obesity Nutrition 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 210000002345 respiratory system Anatomy 0.000 description 1
- 230000004622 sleep time Effects 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
- 208000024891 symptom Diseases 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
- 238000010998 test method Methods 0.000 description 1
- 210000003813 thumb Anatomy 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/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
- A61B5/024—Detecting, measuring or recording pulse rate or heart rate
- A61B5/02405—Determining heart rate variability
-
- 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/318—Heart-related electrical modalities, e.g. electrocardiography [ECG]
- A61B5/333—Recording apparatus specially adapted therefor
-
- 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/318—Heart-related electrical modalities, e.g. electrocardiography [ECG]
- A61B5/346—Analysis of electrocardiograms
- A61B5/349—Detecting specific parameters of the electrocardiograph cycle
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7203—Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal
- A61B5/7207—Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal of noise induced by motion artifacts
- A61B5/7214—Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal of noise induced by motion artifacts using signal cancellation, e.g. based on input of two identical physiological sensors spaced apart, or based on two signals derived from the same sensor, for different optical wavelengths
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B2562/00—Details of sensors; Constructional details of sensor housings or probes; Accessories for sensors
- A61B2562/02—Details of sensors specially adapted for in-vivo measurements
- A61B2562/0219—Inertial sensors, e.g. accelerometers, gyroscopes, tilt switches
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Cardiology (AREA)
- Engineering & Computer Science (AREA)
- Surgery (AREA)
- Medical Informatics (AREA)
- Veterinary Medicine (AREA)
- Physics & Mathematics (AREA)
- Public Health (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- General Health & Medical Sciences (AREA)
- Molecular Biology (AREA)
- Animal Behavior & Ethology (AREA)
- Physiology (AREA)
- Signal Processing (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Psychiatry (AREA)
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
- Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
Abstract
【課題】 ２４時間心電計の機能を有すると共に、高精度のＲ−Ｒ間隔心拍計の機能を有し、心拍変動に対する呼吸運動による影響を排除することができる心拍変動解析方法および解析装置を提供する。また、詳細な心拍変動の解析が可能な心拍変動解析装置を提供する。
【解決手段】 心拍変動解析装置は測定部１と解析表示部２からなる。測定部１は胸郭部位と横隔膜部位にそれぞれに取り付けられる第１、第２の心筋活動電位測定手段１１、１２と、３軸加速度測定手段１３、１４、測定された心筋活動電位信号と加速度信号を解析表示部２に送信する送信機１５、１６からなる。
解析表示部２は、送信機１５、１６から送信された信号を受信する受信機２１と、受信機２１で受信された信号をコンピュータに接続するインターフェース２２と、入力された心筋活動電位信号と加速度信号から心拍変動を解析し表示および／または印刷することが可能なコンピュータ２３からなる。
【選択図】 図１
【解決手段】 心拍変動解析装置は測定部１と解析表示部２からなる。測定部１は胸郭部位と横隔膜部位にそれぞれに取り付けられる第１、第２の心筋活動電位測定手段１１、１２と、３軸加速度測定手段１３、１４、測定された心筋活動電位信号と加速度信号を解析表示部２に送信する送信機１５、１６からなる。
解析表示部２は、送信機１５、１６から送信された信号を受信する受信機２１と、受信機２１で受信された信号をコンピュータに接続するインターフェース２２と、入力された心筋活動電位信号と加速度信号から心拍変動を解析し表示および／または印刷することが可能なコンピュータ２３からなる。
【選択図】 図１
Description
本発明は、高性能の時間分解能で測定された胸郭部位および横隔膜部位からの心筋活動電位信号と、加速度計で測定された加速度信号とから、心拍変動を解析し表示および／または印刷することが可能な心拍変動解析方法および解析装置に関する。
ホルター心電図は、通常の日常生活にあるがままの状況下で、１拍毎の心電図を長時間連続して記録することにより心機能を評価する検査法であり、現在は、生活習慣病に映る心血管機能の異常を精確かつ忠実に抽出する手法として発展しつつある。
最近、ホルター心電計の時間分解能の向上により、測定される心拍にゆらぎが存在することが見いだされており、この心拍変動を測定することにより、患者のレベルで言えば、例えば、突然死や、心機能低下、心室リモデリング、脳梗塞などの循環器疾患、糖尿病や肥満症、高脂血症、高血圧などの生活習慣病、の重症化予測と重症化予防が期待され、また、健康人レベルでは未病の治癒が期待されている。
また、心拍変動は特に自律神経系によってその変動が強く支配されており、心拍の変動を解析することによって、自律神経機能を充分な精度をもって類推することが可能になると期待されている。
最近、ホルター心電計の時間分解能の向上により、測定される心拍にゆらぎが存在することが見いだされており、この心拍変動を測定することにより、患者のレベルで言えば、例えば、突然死や、心機能低下、心室リモデリング、脳梗塞などの循環器疾患、糖尿病や肥満症、高脂血症、高血圧などの生活習慣病、の重症化予測と重症化予防が期待され、また、健康人レベルでは未病の治癒が期待されている。
また、心拍変動は特に自律神経系によってその変動が強く支配されており、心拍の変動を解析することによって、自律神経機能を充分な精度をもって類推することが可能になると期待されている。
心拍変動を解析する心拍変動測定装置としては、統計手法に基づいて異常間隔の不整脈やアーチファクトの影響を回避できるようにしたものが提案されている（特許文献１）。
このものは、心拍信号検出手段と、心拍間隔計時手段と、心拍間隔データを大きさの順に並べるソート手段と、中間領域を例えば四分位範囲法により設定して第１四分位及び第３四分位間の差に相当する四分位範囲を心拍変動の評価用指標として演算する心拍データ解析手段と、指標を表示もしくはプリントアウトする出力手段とを備えている。
このものは、心拍信号検出手段と、心拍間隔計時手段と、心拍間隔データを大きさの順に並べるソート手段と、中間領域を例えば四分位範囲法により設定して第１四分位及び第３四分位間の差に相当する四分位範囲を心拍変動の評価用指標として演算する心拍データ解析手段と、指標を表示もしくはプリントアウトする出力手段とを備えている。
然しながら、その後、ホルター心電計の時間分解能が改良されている（８msec）とはいえ、まだ十分ではなく、Ｒ波の頂点の位置および形状を精確に確認できるものではなく、従って、Ｒ−Ｒ間隔の精確な測定ができなかった。
また、心拍変動には呼吸系のものもあるが、肺（呼吸運動による）の影響を排除する心電計は従来無く、肺の影響を排除しようとすれば、従来は、入院して面倒な検査をする必要があった。
また、心拍変動には呼吸系のものもあるが、肺（呼吸運動による）の影響を排除する心電計は従来無く、肺の影響を排除しようとすれば、従来は、入院して面倒な検査をする必要があった。
本発明は、如上の事情に鑑みてなされたもので、２４時間心電計の機能を有すると共に、高精度のＲ−Ｒ間隔心拍計の機能を有し、心拍変動に対する呼吸運動による影響を排除することができる心拍変動解析方法および解析装置を提供することを目的とする。
また、詳細な心拍変動の解析が可能な心拍変動解析装置を提供することを目的とする。
また、詳細な心拍変動の解析が可能な心拍変動解析装置を提供することを目的とする。
本発明の心拍変動解析方法は、胸郭部位から３msec以下の時間分解能で測定された心筋活動電位信号と、横隔膜部位から３msec以下の時間分解能で測定された心筋活動電位信号と、身体の加速度信号とから心拍変動を解析することを特徴とする。ここで、解析データからの病気（未病を含む）の診断を容易にするために、心筋活動電位信号の解析により、瞬時心拍数時系列データと心拍数変動リズム比率データ、心拍変動揺らぎ周波数ＨＦ成分時系列データ、心拍変動揺らぎ周波数ＨＦ成分リズム比率データ、心拍変動揺らぎ周波数ＬＦ／ＨＦ成分時系列データ、心拍変動揺らぎ周波数ＬＦ／ＨＦ成分リズム比率データ、の６種類のデータが演算され表示されるのが好ましい。また、使いやすさとコストの面から、心筋活動電位信号は時間分解能０．１〜３msecのホルター心電計で測定されるのが好ましい。
また、本発明の心拍変動解析装置は、胸郭部位からの心筋活動電位信号を測定するための３msec以下の時間分解能を有する第１の心筋活動電位測定手段と、横隔膜部位からの心筋活動電位信号を測定するための３msec以下の時間分解能を有する第２の心筋活動電位測定手段と、身体の加速度信号を測定するための３軸加速度測定手段とを含んでなる測定部と、該測定部から送信される心筋活動電位信号と加速度信号を受信して心拍変動を解析し表示する解析表示部、を含んでなることを特徴とする。ここで、心筋活動電位信号の解析により、瞬時心拍数時系列データと心拍数変動リズム比率データ、心拍変動揺らぎ周波数ＨＦ成分時系列データ、心拍変動揺らぎ周波数ＨＦ成分リズム比率データ、心拍変動揺らぎ周波数ＬＦ／ＨＦ成分時系列データ、心拍変動揺らぎ周波数ＬＦ／ＨＦ成分リズム比率データ、の６種類のデータが求められ表示される様にするのが好ましい。また、心筋活動電位信号は、時間分解能０．１〜３msecのホルター心電計で測定される様にするのが好ましい。また、必要ならば、ホルター心電計に傾斜計を内蔵させてもよい。
以上、一般的に本発明を記述したが、より一層の理解は、いくつかの特定の実施例を参照することによって得ることが出来る。これらの実施例は本明細書に例示の目的のためにのみ提供されるものであり、他の旨が特定されない限り、限定的なものではない。
以上、一般的に本発明を記述したが、より一層の理解は、いくつかの特定の実施例を参照することによって得ることが出来る。これらの実施例は本明細書に例示の目的のためにのみ提供されるものであり、他の旨が特定されない限り、限定的なものではない。
本発明によれば、以下のような効果が期待できる。すなわち、心拍変動には心臟系のものと呼吸系のものとがあるが、本発明の心拍変動解析装置は、呼吸運動を精確に捉えるために、胸郭と横隔膜（上腹部）の２カ所から加速度信号を測定できるように構成しているので、肺の影響（呼吸運動による影響）を排除することができる。また、加速度計を使用することにより、睡眠時無呼吸を検知して呼吸運動を物理的に測定することが可能になった。
また、３msec以下の高い時間分解能で心筋活動電位を測定しているので、詳細な心拍変動解析が可能である。また、必要なデータをグラフ等のレポートとして出力できるように構成しているので、これまで経験則と所見の組み合わせによってのみ判断されてきた自律神経機能の状態を、客観的に判断するための手がかりとして利用可能である。
また、３msec以下の高い時間分解能で心筋活動電位を測定しているので、詳細な心拍変動解析が可能である。また、必要なデータをグラフ等のレポートとして出力できるように構成しているので、これまで経験則と所見の組み合わせによってのみ判断されてきた自律神経機能の状態を、客観的に判断するための手がかりとして利用可能である。
胸郭部位からの心筋活動電位信号を測定するための時間分解能１msecの第１の心筋活動電位測定手段と加速度信号を測定するための３軸加速度計と送信機からなる第１の心筋活動電位測定計と、横隔膜部位からの心筋活動電位信号を測定するための時間分解能１msecの第２の心筋活動電位測定手段と加速度信号を測定するための３軸加速度計と送信機からなる第２の心筋活動電位測定計、とからなるホルター心電計と、このホルター心電計から送信される心筋活動電位信号と加速度信号を受信して心拍変動を解析し表示する解析表示部、を含む。解析によって求められたデータは、瞬時心拍数時系列データと心拍数変動リズム比率データ、心拍変動揺らぎ周波数ＨＦ成分時系列データ、心拍変動揺らぎ周波数ＨＦ成分リズム比率データ、心拍変動揺らぎ周波数ＬＦ／ＨＦ成分時系列データ、心拍変動揺らぎ周波数ＬＦ／ＨＦ成分リズム比率データ、の６種類のデータに分けられて表示される様になっている。
先ず、実施例１について図１〜図３を用いて説明する。
図１は本発明の心拍変動装置のブロック図であり、図２は本発明のホルター心電計の説明図、図３は図１の受信機とインターフェースを内蔵する筐体を示す説明図である。
実施例1の心拍変動解析装置は、図１に示すように、測定部１と解析表示部２からなる。
測定部１は胸郭部位に取り付けられる第１の心筋活動電位測定手段１１と横隔膜部位に取り付けられる第２の心筋活動電位測定手段１２、３軸加速度測定手段１３、１４、心筋活動電位測定手段１１で測定された心筋活動電位信号と加速度測定手段１３で測定された加速度信号、心筋活動電位測定手段１２で測定された心筋活動電位信号と加速度測定手段１４で測定された加速度信号をそれぞれ解析表示部２に送信する送信機１５、１６からなる。
また、解析表示部２は、送信機１５、１６から送信された心筋活動電位信号と加速度信号を受信する受信機２１と、受信機２１で受信された信号をコンピュータに接続するインターフェース２２と、このインターフェース２２を介して入力された心筋活動電位信号と加速度信号から心拍変動を解析し表示および／または印刷することが可能なコンピュータ２３からなる。
図１は本発明の心拍変動装置のブロック図であり、図２は本発明のホルター心電計の説明図、図３は図１の受信機とインターフェースを内蔵する筐体を示す説明図である。
実施例1の心拍変動解析装置は、図１に示すように、測定部１と解析表示部２からなる。
測定部１は胸郭部位に取り付けられる第１の心筋活動電位測定手段１１と横隔膜部位に取り付けられる第２の心筋活動電位測定手段１２、３軸加速度測定手段１３、１４、心筋活動電位測定手段１１で測定された心筋活動電位信号と加速度測定手段１３で測定された加速度信号、心筋活動電位測定手段１２で測定された心筋活動電位信号と加速度測定手段１４で測定された加速度信号をそれぞれ解析表示部２に送信する送信機１５、１６からなる。
また、解析表示部２は、送信機１５、１６から送信された心筋活動電位信号と加速度信号を受信する受信機２１と、受信機２１で受信された信号をコンピュータに接続するインターフェース２２と、このインターフェース２２を介して入力された心筋活動電位信号と加速度信号から心拍変動を解析し表示および／または印刷することが可能なコンピュータ２３からなる。
第１、第２の心筋活動電位測定手段１１、１２は３msec以下の時間分解能を有しており、測定された心筋活動電位信号は、Ｒ波認識アルゴリズムを用いてＲ波の頂点が自動認識され、Ｒ−Ｒ間隔時系列がサンプリング間隔３msec以下の時間分解能でＡ／Ｄ変換される。解析表示部２へはこのＡ／Ｄ変換された時系列データ(生データ)が送信される。
心筋活動電位信号を測定するために採用されるホルター心電計Ｈ１、Ｈ２は、図２に示すように、筐体１１１（１２１）と、これに内蔵された心電回路１１２（１２２）、電極１１３（１２３）、３軸加速度計（３軸加速度測定手段）１３（１４）、および送信機１５（１６）からなる。ホルター心電計Ｈ１、Ｈ２で測定された心筋活動電位信号と加速度信号は、送信機１５（１６）によって解析表示部２に送信されるようになっている。図中、１１４（１２４）は動作確認用インジケータ、１１３（１２３）は電極であり、この電極１１３（１２３）と心電回路１１２（１２２）で心筋活動電位測定手段１１（１２）を構成している。
心筋活動電位信号を測定するために採用されるホルター心電計Ｈ１、Ｈ２は、図２に示すように、筐体１１１（１２１）と、これに内蔵された心電回路１１２（１２２）、電極１１３（１２３）、３軸加速度計（３軸加速度測定手段）１３（１４）、および送信機１５（１６）からなる。ホルター心電計Ｈ１、Ｈ２で測定された心筋活動電位信号と加速度信号は、送信機１５（１６）によって解析表示部２に送信されるようになっている。図中、１１４（１２４）は動作確認用インジケータ、１１３（１２３）は電極であり、この電極１１３（１２３）と心電回路１１２（１２２）で心筋活動電位測定手段１１（１２）を構成している。
ホルター心電計Ｈ１、Ｈ２としては時間分解能０．１〜３msecのものが使用されており、これを小型化することにより、電極シール（図示していない）などを用いて皮膚に貼着できるようになっている。また、電極シールとして防水のものを用いると、入浴が可能になるので、心筋活動電位の２４時間測定が可能になる。尚、心電計は、時間分解能が０．１〜３msecのものであればどの様なものでも採用可能であるが、皮膚に添付可能なサイズと形状のものが好ましい。時間分解能０．１msec以下では装置が大きくなりすぎるので実用的ではなく、時間分解能３msec以上では心拍の間隔（Ｒ−Ｒ間隔として測定）を測定するためのＲ波の頂点の位置が不明確なので心拍変動を精確に測定することができない。３軸加速度計１３、１４としては、圧電式加速度センサーが採用される。加速度計は必ずしも心電計１１、１２に内蔵させなくともよい。送信機１５としては、セラミックモジュールの２．４ＧＨｚ帯高度化小電力データ通信システムを採用しており、通信方式はＴＤＭＡを採用している。
解析表示部２は、測定部１から送信された心筋活動電位信号と加速度信号を受信する受信機２１と、受信機２１で受信した心筋活動電位信号と加速度信号をコンピュータ２３に接続するインターフェース２と、コンピュータ２３からなる。
受信機２１は、図３に示すように、インターフェース２２と一緒に筐体２１１に収容されており、スイッチ２１２をオンにすると駆動するようになっている。受信した信号はインターフェース２２の記憶媒体に記憶され、随時コンピュータ２３に接続できるようになっている。尚、２１３は動作確認用インジケータである。
受信機２１としては、セラミックモジュールの２．４ＧＨｚ帯高度化小電力データ通信システムを採用しており、通信方式はＴＤＭＡを採用している。コンピュータ２３には、胸郭部位および横隔膜部位からの心筋活動電位信号と加速度信号とから心拍変動を解析するためのプログラムが入力されており、解析された心拍変動は、瞬時心拍数時系列データと心拍数変動リズム比率データ、心拍変動揺らぎ周波数ＨＦ成分時系列データ、心拍変動揺らぎ周波数ＨＦ成分リズム比率データ、心拍変動揺らぎ周波数ＬＦ／ＨＦ成分時系列データ、心拍変動揺らぎ周波数ＬＦ／ＨＦ成分リズム比率データ、の６種類のデータとして表示し印刷することが可能になっている。
受信機２１は、図３に示すように、インターフェース２２と一緒に筐体２１１に収容されており、スイッチ２１２をオンにすると駆動するようになっている。受信した信号はインターフェース２２の記憶媒体に記憶され、随時コンピュータ２３に接続できるようになっている。尚、２１３は動作確認用インジケータである。
受信機２１としては、セラミックモジュールの２．４ＧＨｚ帯高度化小電力データ通信システムを採用しており、通信方式はＴＤＭＡを採用している。コンピュータ２３には、胸郭部位および横隔膜部位からの心筋活動電位信号と加速度信号とから心拍変動を解析するためのプログラムが入力されており、解析された心拍変動は、瞬時心拍数時系列データと心拍数変動リズム比率データ、心拍変動揺らぎ周波数ＨＦ成分時系列データ、心拍変動揺らぎ周波数ＨＦ成分リズム比率データ、心拍変動揺らぎ周波数ＬＦ／ＨＦ成分時系列データ、心拍変動揺らぎ周波数ＬＦ／ＨＦ成分リズム比率データ、の６種類のデータとして表示し印刷することが可能になっている。
図4−１、図４−２、図４−３は、心拍変動の解析プログラムを示すフローチャートである。
コンピュータ２３にＲ−Ｒ間隔生データ(２４時間)が入力されると（Ｓ１）、先ず不整脈などの異常値が除去された後線形補間が行われ（Ｓ２）、Ｒ−Ｒ間隔フィルター後データ(２４時間)が出力される（Ｓ３）。このフィルター後データは、５分間毎の平均心拍数計算により（Ｓ４）心拍数の推移時系列（Ａ）として出力される（Ｓ５）一方、線形補間を用いて０．５秒の等間隔に再サンプリングされて（Ｓ６）、再サンプリング後のＲ−Ｒ間隔時系列データ（２４時間）として出力される（Ｓ７）。この再サンプリング後データは５分間毎に分割された（Ｓ８）後、周波数解析されて（Ｓ９）、５分間毎のパワースペクトル密度（ＰＳＤ）として出力され（Ｓ１０）、更に０．１５〜０．４Ｈｚのパワー：ＨＦと０．０４〜０．１５Ｈｚのパワー：ＬＦに分割され（Ｓ１１）、ＨＦは√ＨＦの推移時系列（Ｂ）として出力され（Ｓ１２）、ＬＦはＨＦで割られて（Ｓ１３）ＬＦ／ＨＦの推移時系列（Ｃ）として出力される（Ｓ１４）。また、Ｓ７で出力された再サンプリング後データは、5分間ずつずらして１００００秒（２．７時間）のデータに分割された（Ｓ１５）後、周波数解析されて（Ｓ１６）、１００００秒間のパワースペクトル密度（ＰＳＤ）として出力され（Ｓ１７）、このパワースペクトル密度はさらに傾き計算され（Ｓ１８）、ＰＳＤの傾きの推移時系列（Ｄ）として出力される（Ｓ１９）。
コンピュータ２３にＲ−Ｒ間隔生データ(２４時間)が入力されると（Ｓ１）、先ず不整脈などの異常値が除去された後線形補間が行われ（Ｓ２）、Ｒ−Ｒ間隔フィルター後データ(２４時間)が出力される（Ｓ３）。このフィルター後データは、５分間毎の平均心拍数計算により（Ｓ４）心拍数の推移時系列（Ａ）として出力される（Ｓ５）一方、線形補間を用いて０．５秒の等間隔に再サンプリングされて（Ｓ６）、再サンプリング後のＲ−Ｒ間隔時系列データ（２４時間）として出力される（Ｓ７）。この再サンプリング後データは５分間毎に分割された（Ｓ８）後、周波数解析されて（Ｓ９）、５分間毎のパワースペクトル密度（ＰＳＤ）として出力され（Ｓ１０）、更に０．１５〜０．４Ｈｚのパワー：ＨＦと０．０４〜０．１５Ｈｚのパワー：ＬＦに分割され（Ｓ１１）、ＨＦは√ＨＦの推移時系列（Ｂ）として出力され（Ｓ１２）、ＬＦはＨＦで割られて（Ｓ１３）ＬＦ／ＨＦの推移時系列（Ｃ）として出力される（Ｓ１４）。また、Ｓ７で出力された再サンプリング後データは、5分間ずつずらして１００００秒（２．７時間）のデータに分割された（Ｓ１５）後、周波数解析されて（Ｓ１６）、１００００秒間のパワースペクトル密度（ＰＳＤ）として出力され（Ｓ１７）、このパワースペクトル密度はさらに傾き計算され（Ｓ１８）、ＰＳＤの傾きの推移時系列（Ｄ）として出力される（Ｓ１９）。
次に、データＡ、Ｂ、Ｃ、Ｄ（Ｓ２０）は全評価領域、覚醒時間帯、睡眠時間帯別の平均値が算出され（Ｓ２１）、表形式で表示される（Ｓ２２）一方、データＡ、Ｂ、Ｃ、Ｄ（Ｓ２０）は周期５時間（5.55e-5 Hz）のローパスフィルタにかけられ（Ｓ２３）、時系列の波形とローパスフィルタ後の波形が重ねて表示される（Ｓ２４）。更にデータＡ、Ｂ、Ｃ、Ｄ（Ｓ２０）からは、各周波数成分の基底リズムが検出され、リズム周期解析されて（Ｓ２５）、その結果が４〜５種類の周期に分類され（Ｓ２６）、各々の周期成分毎のエネルギー量が算出され（Ｓ２７）、Ａ、Ｂ、Ｃ、Ｄ別に各成分のパーセントが計算されて、円グラフ表示される（Ｓ２８）。更にまた、データＡ、Ｂ、Ｃ、Ｄ（Ｓ２０）は、カオス解析（ターケンスの埋め込み、アトラクタ構成、リアプノフ指数計算）され（Ｓ２９）、リアプノフ指数がＡ、Ｂ、Ｃ、Ｄ別にカオティックスコアとして表示される（Ｓ３０）。
尚、Ｓ２４で波形表示されたデータＡ、Ｂ、Ｃはそれぞれ、瞬時心拍数時系列データと心拍変動揺らぎ周波数ＨＦ成分時系列データ、心拍変動揺らぎ周波数ＬＦ／ＨＦ成分時系列データであり、図５〜７にそれぞれの表示例を示す。また、Ｓ２８で円グラフ表示されたデータＡ、Ｂ、Ｃはそれぞれ心拍数変動リズム比率データと心拍変動揺らぎ周波数ＨＦ成分リズム比率データ、心拍変動揺らぎ周波数ＬＦ／ＨＦ成分リズム比率データであり、図８〜１０にそれぞれの表示例を示す。
これらのデータのうち、瞬時心拍数時系列データからは不整脈及び心悸亢進症状の有無の推測が行われ、心拍変動揺らぎ周波数ＨＦ成分時系列データからは副交感神経活動の時系列推移の推測が、心拍変動揺らぎ周波数ＬＦ／ＨＦ成分時系列データからは交感神経活動の時系列推移の推測が行われる。また、心拍数変動リズム比率データからは心肺機能の活動安定性の推測が行われ、また、心拍変動揺らぎ周波数ＨＦ成分リズム比率データからは副交感神経活動の安定性の推測が、揺らぎ周波数ＬＦ／ＨＦ成分リズム比率データからは交感神経活動の安定性の推測が行われる。
これらのデータのうち、瞬時心拍数時系列データからは不整脈及び心悸亢進症状の有無の推測が行われ、心拍変動揺らぎ周波数ＨＦ成分時系列データからは副交感神経活動の時系列推移の推測が、心拍変動揺らぎ周波数ＬＦ／ＨＦ成分時系列データからは交感神経活動の時系列推移の推測が行われる。また、心拍数変動リズム比率データからは心肺機能の活動安定性の推測が行われ、また、心拍変動揺らぎ周波数ＨＦ成分リズム比率データからは副交感神経活動の安定性の推測が、揺らぎ周波数ＬＦ／ＨＦ成分リズム比率データからは交感神経活動の安定性の推測が行われる。
一方、加速度信号については、加速度信号（３軸加速度計データ）がコンピュータ２３に入力されると（Ｓ３１）、２つの加速度計データが時系列に同期され（Ｓ３２）、呼吸運動解析が行われて（Ｓ３３）心拍変動への呼吸運動の関与が推測される（Ｓ３４）一方、就寝時間解析が行われて（Ｓ３５）、睡眠時無呼吸発作の推測が行われる（Ｓ３６）。
１ 測定部
１１ 第１の心筋活動電位測定手段
１２ 第２の心筋活動電位測定手段
１３、１４ ３軸加速度測定手段（３軸加速度計）
１５、１６ 送信機
Ｈ１、Ｈ２ ホルター心電図
１１１、１２１ 筐体
１１２、１２２ 心電回路
１１３、１２３ 電極
１１４、１２４ 動作確認用インジケータ
２ 解析表示部
２１ 受信機
２１１ 筐体
２１２ スイッチ
２１３ 動作確認用インジケータ
２２ インターフェース
２３ コンピュータ
１１ 第１の心筋活動電位測定手段
１２ 第２の心筋活動電位測定手段
１３、１４ ３軸加速度測定手段（３軸加速度計）
１５、１６ 送信機
Ｈ１、Ｈ２ ホルター心電図
１１１、１２１ 筐体
１１２、１２２ 心電回路
１１３、１２３ 電極
１１４、１２４ 動作確認用インジケータ
２ 解析表示部
２１ 受信機
２１１ 筐体
２１２ スイッチ
２１３ 動作確認用インジケータ
２２ インターフェース
２３ コンピュータ
Claims (6)
- 胸郭部位から３msec以下の時間分解能で測定された心筋活動電位信号と、横隔膜部位から３msec以下の時間分解能で測定された心筋活動電位信号と、身体の加速度信号とから心拍変動を解析する心拍変動解析方法。
- 心筋活動電位信号の解析により、瞬時心拍数時系列データと心拍数変動リズム比率データ、心拍変動揺らぎ周波数ＨＦ成分時系列データ、心拍変動揺らぎ周波数ＨＦ成分リズム比率データ、心拍変動揺らぎ周波数ＬＦ／ＨＦ成分時系列データ、心拍変動揺らぎ周波数ＬＦ／ＨＦ成分リズム比率データ、の６種類のデータが演算され表示される請求項１に記載の方法。
- 心筋活動電位信号が時間分解能０．１〜３msecのホルター心電計で測定される請求項1または２に記載の方法。
- 胸郭部位からの心筋活動電位信号を測定するための３msec以下の時間分解能を有する第１の心筋活動電位測定手段と、横隔膜部位からの心筋活動電位信号を測定するための３msec以下の時間分解能を有する第２の心筋活動電位測定手段と、加速度信号を測定するための３軸加速度測定手段とを含んでなる測定部と、該測定部から送信される心筋活動電位信号と加速度信号を受信して心拍変動を解析し表示する解析表示部、を含んでなる心拍変動解析装置。
- 心筋活動電位信号の解析により、瞬時心拍数時系列データと心拍数変動リズム比率データ、心拍変動揺らぎ周波数ＨＦ成分時系列データ、心拍変動揺らぎ周波数ＨＦ成分リズム比率データ、心拍変動揺らぎ周波数ＬＦ／ＨＦ成分時系列データ、心拍変動揺らぎ周波数ＬＦ／ＨＦ成分リズム比率データ、の６種類のデータが求められ表示される請求項４に記載の装置。
- 心筋活動電位信号が時間分解能０．１〜３msecのホルター心電計で測定される請求項４または５に記載の装置。
Priority Applications (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2007271434A JP5221926B2 (ja) | 2007-10-18 | 2007-10-18 | 心拍変動解析方法および解析装置 |
| EP08018066A EP2050387A1 (en) | 2007-10-18 | 2008-10-15 | Heart-rate variability analysis method and analysis device |
| US12/254,468 US20090105601A1 (en) | 2007-10-18 | 2008-10-20 | Heart-rate variability analysis method and analysis device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2007271434A JP5221926B2 (ja) | 2007-10-18 | 2007-10-18 | 心拍変動解析方法および解析装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2009095552A true JP2009095552A (ja) | 2009-05-07 |
| JP5221926B2 JP5221926B2 (ja) | 2013-06-26 |
Family
ID=40350205
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2007271434A Active JP5221926B2 (ja) | 2007-10-18 | 2007-10-18 | 心拍変動解析方法および解析装置 |
Country Status (3)
| Country | Link |
|---|---|
| US (1) | US20090105601A1 (ja) |
| EP (1) | EP2050387A1 (ja) |
| JP (1) | JP5221926B2 (ja) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2014510545A (ja) * | 2010-10-27 | 2014-05-01 | コーニンクレッカ フィリップス エヌ ヴェ | 対象者のサーカディアンリズムの特徴を決定するための方法 |
| JP2019103614A (ja) * | 2017-12-12 | 2019-06-27 | 株式会社テクノアクセルネットワークス | 心拍データ出力装置、心拍データ出力方法、心拍データ出力プログラム及びペースメーカ |
| JP2021079035A (ja) * | 2019-11-22 | 2021-05-27 | 株式会社Ｂｄａｓｈ | 自律神経機能評価装置及び自律神経機能評価プログラム |
Families Citing this family (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US9042969B2 (en) * | 2009-08-12 | 2015-05-26 | Angel Medical Systems, Inc. | Heart rate correction system and methods for the detection of cardiac events |
| NZ630602A (en) * | 2012-09-19 | 2016-10-28 | Resmed Sensor Technologies Ltd | System and method for determining sleep stage |
| US9737258B2 (en) * | 2013-01-31 | 2017-08-22 | Massachusetts Institute Of Technology | Systems and methods for detecting a physiological abnormality in a patient by using cardiac or other chaos in combination with non-increasing parasympathetic modulation |
| US10022068B2 (en) | 2013-10-28 | 2018-07-17 | Covidien Lp | Systems and methods for detecting held breath events |
| CN110946573A (zh) * | 2019-11-01 | 2020-04-03 | 东软集团股份有限公司 | 心脏停搏的检测装置、检测模型的训练装置、方法及设备 |
| CN114129145A (zh) * | 2021-12-16 | 2022-03-04 | 成都怡康科技有限公司 | 一种跟踪学生体育课堂跳绳心率变化的方法 |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH1156801A (ja) * | 1997-08-22 | 1999-03-02 | Mitsuaki Yamamoto | 携帯型プログラマブル生体情報長期計測・記憶システム |
| JPH11178804A (ja) * | 1997-12-25 | 1999-07-06 | Matsushita Electric Ind Co Ltd | 体調チェック装置 |
| WO2003015630A1 (fr) * | 2001-08-10 | 2003-02-27 | Dainippon Pharmaceutical Co., Ltd. | Systeme de communication de donnees de signal biologique par discrimination et procede de communication associe |
| JP2004121864A (ja) * | 1995-09-13 | 2004-04-22 | Seiko Epson Corp | 健康状態管理装置 |
| JP2007044442A (ja) * | 2005-08-12 | 2007-02-22 | Nippon Sekijiyuujishiya | 周波数解析装置および自律神経活動評価装置 |
| JP3132458U (ja) * | 2007-03-30 | 2007-06-07 | 静岡県 | 携帯型微弱信号収録装置 |
| JP2007521852A (ja) * | 2003-11-18 | 2007-08-09 | ヴィーヴォメトリックス インコーポレイテッド | 移動式生理学的モニタリングからのデータを処理する方法及びシステム |
Family Cites Families (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5205295A (en) * | 1991-06-03 | 1993-04-27 | Del Mar Avionics | Method and apparatus for holter recorder with high resolution signal averaging capability for late potential analysis |
| JPH05317277A (ja) | 1992-05-22 | 1993-12-03 | Yoshihiko Suzuki | 心拍変動解析方法及び測定装置 |
| US5840038A (en) * | 1997-05-29 | 1998-11-24 | Marquette Medical Systems, Inc. | Method and apparatus for signal averaging and analyzing high resolution P wave signals from an electrocardiogram |
| US6370423B1 (en) * | 1998-10-05 | 2002-04-09 | Juan R. Guerrero | Method for analysis of biological voltage signals |
| CA2522006C (en) * | 2003-04-10 | 2013-05-28 | Vivometrics, Inc. | Systems and methods for respiratory event detection |
| EP1609415A1 (en) * | 2004-06-23 | 2005-12-28 | Instrumentarium Corporation | Method and a device for recording signals |
| US7668588B2 (en) * | 2006-03-03 | 2010-02-23 | PhysioWave, Inc. | Dual-mode physiologic monitoring systems and methods |
| EP2020913A4 (en) * | 2006-05-08 | 2010-02-24 | A M P S Llc | METHOD AND DEVICE FOR EXTRACTING OPTICAL ECG HOLTER READING |
| US20090024176A1 (en) * | 2007-07-17 | 2009-01-22 | Joonkyoo Anthony Yun | Methods and devices for producing respiratory sinus arrhythmia |
-
2007
- 2007-10-18 JP JP2007271434A patent/JP5221926B2/ja active Active
-
2008
- 2008-10-15 EP EP08018066A patent/EP2050387A1/en not_active Withdrawn
- 2008-10-20 US US12/254,468 patent/US20090105601A1/en not_active Abandoned
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2004121864A (ja) * | 1995-09-13 | 2004-04-22 | Seiko Epson Corp | 健康状態管理装置 |
| JPH1156801A (ja) * | 1997-08-22 | 1999-03-02 | Mitsuaki Yamamoto | 携帯型プログラマブル生体情報長期計測・記憶システム |
| JPH11178804A (ja) * | 1997-12-25 | 1999-07-06 | Matsushita Electric Ind Co Ltd | 体調チェック装置 |
| WO2003015630A1 (fr) * | 2001-08-10 | 2003-02-27 | Dainippon Pharmaceutical Co., Ltd. | Systeme de communication de donnees de signal biologique par discrimination et procede de communication associe |
| JP2007521852A (ja) * | 2003-11-18 | 2007-08-09 | ヴィーヴォメトリックス インコーポレイテッド | 移動式生理学的モニタリングからのデータを処理する方法及びシステム |
| JP2007044442A (ja) * | 2005-08-12 | 2007-02-22 | Nippon Sekijiyuujishiya | 周波数解析装置および自律神経活動評価装置 |
| JP3132458U (ja) * | 2007-03-30 | 2007-06-07 | 静岡県 | 携帯型微弱信号収録装置 |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2014510545A (ja) * | 2010-10-27 | 2014-05-01 | コーニンクレッカ フィリップス エヌ ヴェ | 対象者のサーカディアンリズムの特徴を決定するための方法 |
| JP2019103614A (ja) * | 2017-12-12 | 2019-06-27 | 株式会社テクノアクセルネットワークス | 心拍データ出力装置、心拍データ出力方法、心拍データ出力プログラム及びペースメーカ |
| JP2021079035A (ja) * | 2019-11-22 | 2021-05-27 | 株式会社Ｂｄａｓｈ | 自律神経機能評価装置及び自律神経機能評価プログラム |
Also Published As
| Publication number | Publication date |
|---|---|
| EP2050387A1 (en) | 2009-04-22 |
| US20090105601A1 (en) | 2009-04-23 |
| JP5221926B2 (ja) | 2013-06-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP5221926B2 (ja) | 心拍変動解析方法および解析装置 | |
| FI126600B (en) | Detection of sleep phenomena using ballistocardiography | |
| US7664606B2 (en) | Apparatus and method for monitoring biological information, and computer program product | |
| US9655532B2 (en) | Wearable physiological monitoring and notification system based on real-time heart rate variability analysis | |
| CN104055499B (zh) | 连续监控人体生理体征的可穿戴式智能手环及方法 | |
| US9265430B2 (en) | Method, system and software product for the measurement of heart rate variability | |
| JP5336803B2 (ja) | 脈波計測装置 | |
| US20060200011A1 (en) | Health management apparatus, health management system, health management method and computer program product | |
| US20100125217A1 (en) | Method and Apparatus for Presenting Heart Rate Variability by Sound and/or Light | |
| JP2007117591A (ja) | 脈波解析装置 | |
| US20110137189A1 (en) | Physiological signal sensing system without time and place contraint and its method | |
| WO2009150744A1 (ja) | 睡眠状態モニタリング装置、モニタリングシステムおよびコンピュータプログラム | |
| JP2007289540A (ja) | ストレスセンサシステム | |
| JP2011120618A (ja) | 自律神経機能評価装置およびプログラム | |
| CN104027109A (zh) | 心房颤动解析装置以及程序 | |
| US9980664B2 (en) | Biological signal measurement system, apparatus, method, and computer program thereof | |
| JP2003079591A (ja) | 生体情報表示方法及び生体情報表示装置 | |
| WO2019138327A1 (en) | Wearable ecg and auscultation monitoring system with sos and remote monitoring | |
| GB2469547A (en) | Measurement of heart rate variability | |
| TWI469765B (zh) | 無線量測睡眠深度的裝置與方法 | |
| US11291406B2 (en) | System for determining a set of at least one cardio-respiratory descriptor of an individual during sleep | |
| US9486154B2 (en) | Device and method for recording physiological signal | |
| JPH0670898A (ja) | 睡眠状態モニター装置 | |
| CN1228290A (zh) | 新生儿生理参数监测仪 | |
| CN114040705A (zh) | 规则性评价信息的分析方法、监护设备及监护系统 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20101013 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20120629 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20120925 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20121121 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20121211 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20130206 |
|
| 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: 20130226 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20130308 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20160315 Year of fee payment: 3 |
|
| R150 | Certificate of patent or registration of utility model |
Ref document number: 5221926 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |