JP3726825B2 - Heart rate signal correction method - Google Patents

Heart rate signal correction method Download PDF

Info

Publication number
JP3726825B2
JP3726825B2 JP2003178559A JP2003178559A JP3726825B2 JP 3726825 B2 JP3726825 B2 JP 3726825B2 JP 2003178559 A JP2003178559 A JP 2003178559A JP 2003178559 A JP2003178559 A JP 2003178559A JP 3726825 B2 JP3726825 B2 JP 3726825B2
Authority
JP
Japan
Prior art keywords
signal
reference value
heartbeat signal
heartbeat
interval
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
Application number
JP2003178559A
Other languages
Japanese (ja)
Other versions
JP2004000660A (en
Inventor
九州男 岩永
浩行 井邊
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Panasonic Electric Works Co Ltd
Original Assignee
Matsushita Electric Works Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Matsushita Electric Works Ltd filed Critical Matsushita Electric Works Ltd
Priority to JP2003178559A priority Critical patent/JP3726825B2/en
Publication of JP2004000660A publication Critical patent/JP2004000660A/en
Application granted granted Critical
Publication of JP3726825B2 publication Critical patent/JP3726825B2/en
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Images

Landscapes

  • Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)

Description

【0001】
【発明の属する技術分野】
本発明は測定された心拍信号の基準値を用いて行う心拍信号の補正方法に関するものである。
【0002】
【従来の技術】
心拍信号を測定してこの心拍情報を利用するにあたり、電気的ノイズや被験者の動きなどによって不正確な情報しか得られないケースが多々ある。このために測定された心拍信号にノイズがのっていないかどうかを検査しなくてはならないのであるが、従来は人間の心拍間隔として妥当な上限値と下限値とを用いて、測定された心拍信号の妥当性をチェックしているにとどまっている。
【0003】
また、測定された心拍信号にたとえばノイズがのっていたり欠損していることが明らかであるにもかかわらず、連続した心拍信号が必要である場合、測定された心拍信号からノイズを除去する補正を行ったり、心拍信号の欠損部分を補完する補正を行うことが必要となる。この時の補正のための基準値としては、従来、心拍間隔の平均値を用いている。
【0004】
【発明が解決しようとする課題】
しかし、上記上限値と下限値とを用いた妥当性のチェックだけではノイズの多いシステムにおいては正確な検査が行われているとは言い難い。また心拍信号の補正についても、人間の心拍間隔の条件を満たすタイミングでノイズが偶然混入している時には、これに対処することができない。加えるに、基準値の作成に際して不正確な情報が混入している可能性も存在しており、このような基準値に基づいた補正では、補正された後の心拍信号も不正確なものとなってしまう。
【0005】
本発明はこのような点に鑑みなされたものであり、その目的とするところは心拍信号検査用基準値として適切な値で導出した基準値を用いて心拍信号の補正を適切に行うことができる心拍信号補正方法を提供するにある。
【0006】
【課題を解決するための手段】
しかして本発明は、間隔のばらつきが所定%以内である心拍信号が所定回数以上連続して測定された時の心拍信号間隔の平均値を心拍信号検査用の基準値とすることに特徴を有している。連続した心拍間隔のばらつきに注目して、ばらつきが所定の範囲に収まっていなければ人間の心拍間隔として妥当な範囲であってもキャンセルして、基準値に反映されないようにするのである。
【0007】
そして本発明に係る心拍信号補正方法は、測定された心拍信号のうち、この基準値の所定%以下の間隔の信号はキャンセルすることに第1の特徴を有し、測定された心拍信号のうち、上記導出方法で得られた基準値の所定%以上の間隔の信号についてはその間に補完のための信号を挿入することに第2の特徴を有している。
【0008】
上記両補正方法を共に行うことが好ましいのはもちろんであり、この場合、基準値の所定%以下の間隔の信号のキャンセルを行い、次いで基準値の所定%以上の間隔の信号についてはその間に補完のための信号を挿入するものとする。
【0009】
信号の挿入に際しては、信号を挿入することによってできるそれぞれの間隔が基準値に最も近くなるような本数及び位置に配置することが好ましい。
【0010】
【発明の実施の形態】
以下本発明について詳述すると、図2は基準値の導出方法及び心拍信号補正方法を実装するためのシステム構成の一例を示している。被験者1からピックアップ2によって検出された心拍による物理的な変化は検出回路3によって電気信号に変換されて二値化回路4によってデジタル信号に変換された後、マイクロコンピュータ(CPU)5に入力される。上記導出方法及び補正方法はCPU5に適した方法で実装され、具体的な作業はCPU5内で行われる。以下このシステムを想定して導出方法及び補正方法について説明する。なお上記システムでは信号の立ち上がりエッジを心拍と認識するようなシステム及び制御プログラムの構成になっているとする。
【0011】
まず基準値の導出についてであるが、これは図1に示すように、ステップ202で検出されたデータをステップ203にて間隔データとする。次いでステップ204において間隔データの妥当性について検査する。これは人間の心拍間隔として正常であるかということについて調べるものであり、従来例でも述べたように、人間の心拍間隔の上限値(H-limit msec)と下限値(L-limit msec)との間に収まっているかどうかをチェックする。なお、この両値については、システムの用途やシステムの対象者の状態などによって多少変化すると考えられる。
【0012】
ステップ204において異常と判断されれば検出の累積数をキャンセルして次の間隔データを待つ。正常であればステップ206において間隔のばらつきが±R%以内であるかどうかを検査する。同一の人物が同一の状態である場合、心拍間隔のばらつきの範囲は、ステップ204で検査した正常範囲よりもはるかに狭くなくてはならず、ばらつきの大きい連続した信号についてはノイズであると判断することができる。ここでの閾値Rについても対象者やシステムに依存するため一概に決めることはできないが、20〜30%の値を好適に使用することができる。
【0013】
そしてステップ206において異常と判断された時には、上記ステップ204の場合と同様に処理され、正常であればデータの検出数に1を加える(ステップ207)。以上の処理を繰り返すことによって、連続するX個の正常な間隔データが検出されたならば、それらの平均値を求めて基準値とする。このときのXについては、少なすぎると値の信頼性が低くなり、多すぎると基準値を得るのに時間がかかりすぎてしまうことから、3〜10ぐらい、好ましくは5程度が適当である。
【0014】
次に測定された心拍信号を上記基準値を用いて補正する補正方法について説明する。まず本来の心拍信号の間にノイズが入った場合のような誤報については、図3に示すように、測定した心拍信号の間隔データが基準値のA%より小さい場合には異常と判定して、今回得られた立ち上がりエッジをキャンセルして検出されなかったものとする。また前回の立ち上がりエッジに関しても間隔データとして正常であると判定はされているが、微妙な前後へのシフトの可能性があるのでキャンセルし前々回の立ち上がりエッジを信頼できる最後のエッジと設定し次の処理を行う。
【0015】
次に本来の心拍信号が欠損している場合のような失報については、図4に示すように、測定した心拍信号の間隔データが基準値のB%より大きい場合、今回と前回の立ち上がりのエッジの間に検出されるべき信号が検出されなかったとして、次に述べる信号の補完を行う。
【0016】
図5に補完ルーチンを示す。ステップ502において、その間に複数の検出されてない信号を含むと考えられる間隔データを基準値で除す。得られた答についてはステップ503及びステップ504において四捨五入などの適切な方法で整数化したのち1を引く。こうして得られた値を補完すべき信号の数とし、その本数の補完エッジ(心拍信号)を間隔データ内に配置するのであるが、この配置にあたってはこれらが等間隔となるようにする。
【0017】
上記補正にあたっての閾値A,Bの値としては、A=80%〜60%、B=125%〜150%の値を好適に用いることができるが、これらの値にしても、この補正方法が実装されるシステムに依存するので一概に決めることはできない。例えばだんだん心拍間隔が広くなることが予測されるシステム(例えば眠りを誘うようなシステム)においてはBの値を大きめに、Aの値を小さめに設定するべきであり、またこのようにシステムの特性を加味することによってより精度の高い補正を行うことができる。
【0018】
上記誤報と失報とに対する補正は、組み合わせて実行することが好ましいのはもちろんであり、この場合のフローを図6に示す。誤報に対する補正の後に、失報に対する補正を行う。基準値と比較して短すぎる間隔データをキャンセルして間隔の広すぎる間隔データとして扱って補完を行うということができて、アルゴリズムの単純化及び高速化を期待することができる。
【0019】
以上のシステムは説明の一例としてリアルタイム処理系を示したが、それに限定するものではなくバッチ処理系にも応用することができる。
【0020】
【発明の効果】
以上のように心拍信号検査用基準値の導出は、ばらつき範囲を設定しているために、短時間でより正確な基準値を導出することができる。そして本発明の心拍信号補正方法においては、基準値の所定%以下の間隔の信号はキャンセルし、基準値の所定%以上の間隔の信号についてはその間に基準値から導いた補完のための信号を挿入するにあたり、正確な基準値を用いるために、より正確な補正を行うことができるものである。特に基準値と比較して短すぎる間隔データは修正するのではなくキャンセルして間隔の広すぎる間隔データとして扱って補完を行うことにより、アルゴリズムの単純化及び高速化を期待することができる。
【0021】
また、信号の挿入による補完に際し、信号を挿入することによってできるそれぞれの間隔が基準値に最も近くなるような本数及び位置に配置することで、実際の心拍信号により近い信号の補完を行うことができる。
【図面の簡単な説明】
【図1】本発明における基準値導出に関するフローチャートである。
【図2】同上のシステム構成図である。
【図3】誤報に対する補正に関するフローチャートである。
【図4】失報に対する補正に関するフローチャートである。
【図5】補完補正についてのフローチャートである。
【図6】補正に関するフローチャートである。
【符号の説明】
1 被験者
2 ピックアップ
3 検出回路
4 二値化回路
5 CPU
[0001]
BACKGROUND OF THE INVENTION
The present invention relates to a heartbeat signal correction method performed using a measured reference value of a heartbeat signal.
[0002]
[Prior art]
In measuring the heartbeat signal and using this heartbeat information, there are many cases where only inaccurate information can be obtained due to electrical noise, subject movement, and the like. For this reason, it is necessary to check whether the measured heartbeat signal is free of noise, but in the past, it was measured using the upper and lower limits that are reasonable for human heartbeat intervals. It only checks the validity of the heart rate signal.
[0003]
Also, a correction that removes noise from the measured heartbeat signal when a continuous heartbeat signal is required, even though it is clear that the measured heartbeat signal has, for example, noise. It is necessary to perform correction to compensate for the missing portion of the heartbeat signal. Conventionally, an average value of heartbeat intervals is used as a reference value for correction at this time.
[0004]
[Problems to be solved by the invention]
However, it is difficult to say that an accurate inspection is performed in a noisy system only by checking validity using the upper limit value and the lower limit value. Further, the correction of the heartbeat signal cannot be dealt with when noise is accidentally mixed in at the timing satisfying the condition of the human heartbeat interval. In addition, there is a possibility that inaccurate information is mixed in the creation of the reference value, and in the correction based on such a reference value, the corrected heart rate signal is also inaccurate. End up.
[0005]
The present invention has been made in view of the above points, and an object of the present invention is to appropriately correct a heartbeat signal using a reference value derived as an appropriate value as a reference value for heartbeat signal inspection. It is to provide a heartbeat signal correction method.
[0006]
[Means for Solving the Problems]
Therefore, the present invention is characterized in that an average value of heartbeat signal intervals when a heartbeat signal having a variation in interval within a predetermined percentage is continuously measured a predetermined number of times or more is used as a reference value for heartbeat signal inspection. are doing. Focusing on the variation in the continuous heartbeat interval, if the variation is not within a predetermined range, even if it is within a reasonable range for the human heartbeat interval, it is canceled and not reflected in the reference value.
[0007]
The heartbeat signal correction method according to the present invention has a first feature in canceling a signal having an interval equal to or less than a predetermined percentage of the reference value among the measured heartbeat signals. A second feature resides in that a signal for complementation is inserted between signals having intervals of a predetermined percentage or more of the reference value obtained by the derivation method.
[0008]
Of course, it is preferable to perform both of the above correction methods. In this case, signals having an interval equal to or less than a predetermined percentage of the reference value are canceled, and then signals having an interval equal to or greater than the predetermined value of the reference value are complemented in the meantime. It is assumed that a signal for is inserted.
[0009]
When inserting the signals, it is preferable to arrange them in such numbers and positions that the intervals generated by inserting the signals are closest to the reference value.
[0010]
DETAILED DESCRIPTION OF THE INVENTION
Hereinafter, the present invention will be described in detail. FIG. 2 shows an example of a system configuration for implementing a reference value deriving method and a heartbeat signal correcting method. A physical change caused by a heartbeat detected from the subject 1 by the pickup 2 is converted into an electrical signal by the detection circuit 3, converted into a digital signal by the binarization circuit 4, and then input to the microcomputer (CPU) 5. . The derivation method and the correction method are implemented by a method suitable for the CPU 5, and specific work is performed in the CPU 5. A derivation method and a correction method will be described below assuming this system. It is assumed that the system has a system and control program configuration that recognizes the rising edge of a signal as a heartbeat.
[0011]
First, regarding the derivation of the reference value, as shown in FIG. 1, the data detected in step 202 is used as interval data in step 203. Step 204 then checks the validity of the interval data. This is to check whether the human heart rate interval is normal. As described in the conventional example, the upper limit value (H-limit msec) and lower limit value (L-limit msec) Check whether it fits in between. It should be noted that these two values are considered to change somewhat depending on the usage of the system and the state of the target person of the system.
[0012]
If it is determined as abnormal in step 204, the accumulated number of detections is canceled and the next interval data is awaited. If normal, it is checked in step 206 whether the variation in interval is within ± R%. If the same person is in the same state, the variation range of the heart rate interval must be much narrower than the normal range examined in step 204, and it is determined that the continuous signal having a large variation is noise. can do. Since the threshold value R here depends on the subject and the system, it cannot be determined unconditionally, but a value of 20 to 30% can be suitably used.
[0013]
If it is determined in step 206 that there is an abnormality, the same processing as in step 204 is performed, and if normal, 1 is added to the number of detected data (step 207). By repeating the above processing, if X consecutive normal interval data are detected, the average value of them is obtained as a reference value. As for X at this time, if the amount is too small, the reliability of the value is lowered, and if it is too large, it takes too much time to obtain the reference value, so about 3 to 10, preferably about 5 is appropriate.
[0014]
Next, a correction method for correcting the measured heartbeat signal using the reference value will be described. First, as shown in FIG. 3, regarding false alarms such as when noise enters between the original heartbeat signals, it is determined that there is an abnormality when the measured heartbeat signal interval data is smaller than A% of the reference value. It is assumed that the rising edge obtained this time is canceled and not detected. Also, although it is determined that the previous rising edge is normal as the interval data, there is a possibility of a slight forward / backward shift, so cancel the previous rising edge and set the last rising edge as the reliable last edge. Process.
[0015]
Next, as shown in FIG. 4, when the original heartbeat signal is missing, if the measured heartbeat signal interval data is greater than B% of the reference value, as shown in FIG. Assuming that a signal to be detected during the edge is not detected, the following signal complement is performed.
[0016]
FIG. 5 shows a complement routine. In step 502, the interval data that is considered to contain a plurality of undetected signals in between is divided by the reference value. The obtained answer is converted to an integer by an appropriate method such as rounding off in steps 503 and 504, and then 1 is subtracted. The value obtained in this way is used as the number of signals to be complemented, and the number of complement edges (heartbeat signals) are arranged in the interval data. In this arrangement, these are set to be equally spaced.
[0017]
As the values of the thresholds A and B in the correction, values of A = 80% to 60% and B = 125% to 150% can be preferably used. Since it depends on the system to be implemented, it cannot be decided unconditionally. For example, in a system in which the heartbeat interval is predicted to gradually increase (for example, a system that induces sleep), the value of B should be set larger and the value of A should be set smaller. By taking into account, correction with higher accuracy can be performed.
[0018]
Needless to say, the correction for the false alarm and the misreport is preferably performed in combination, and the flow in this case is shown in FIG. After correction for false alarms, correction for false alarms is performed. It is possible to cancel the interval data that is too short compared to the reference value, treat it as interval data that is too wide, and perform complementation, so that simplification and speeding up of the algorithm can be expected.
[0019]
The above system shows a real-time processing system as an example of the description, but the present invention is not limited to this and can be applied to a batch processing system.
[0020]
【The invention's effect】
As described above, since the derivation of the reference value for heartbeat signal inspection sets a variation range, a more accurate reference value can be derived in a short time. In the heartbeat signal correction method of the present invention, signals having an interval equal to or less than a predetermined percentage of the reference value are canceled, and signals having an interval equal to or greater than the predetermined percentage of the reference value are interpolated with signals for complementation derived from the reference value. In insertion, more accurate correction can be performed in order to use an accurate reference value. In particular, interval data that is too short compared to the reference value is not corrected, but is canceled and treated as interval data that is too wide.
[0021]
In addition, when complementing by inserting signals, signals that are closer to the actual heartbeat signal can be complemented by arranging them at the number and position where each interval that can be obtained by inserting the signal is closest to the reference value. it can.
[Brief description of the drawings]
FIG. 1 is a flowchart regarding reference value derivation in the present invention.
FIG. 2 is a system configuration diagram of the above.
FIG. 3 is a flowchart regarding correction for false alarms.
FIG. 4 is a flowchart relating to correction for missing information.
FIG. 5 is a flowchart for complementary correction.
FIG. 6 is a flowchart regarding correction.
[Explanation of symbols]
1 Test subject 2 Pickup 3 Detection circuit 4 Binarization circuit 5 CPU

Claims (4)

間隔のばらつきが所定%以内である心拍信号が所定回数以上連続して測定された時の心拍信号間隔の平均値を心拍信号検査用の基準値とし、測定された心拍信号のうち、この基準値の所定%以下の間隔の信号はキャンセルすることを特徴とする心拍信号補正方法。The average value of the heartbeat signal intervals when the heartbeat signal whose interval variation is within a predetermined percentage is continuously measured a predetermined number of times or more is used as a reference value for heartbeat signal inspection, and this reference value among the measured heartbeat signals A method for correcting a heartbeat signal, wherein a signal having an interval equal to or less than a predetermined percentage is canceled. 間隔のばらつきが所定%以内である心拍信号が所定回数以上連続して測定された時の心拍信号間隔の平均値を心拍信号検査用の基準値とし、測定された心拍信号のうち、基準値の所定%以上の間隔の信号についてはその間に補完のための信号を挿入することを特徴とする心拍信号補正方法。The average value of the heartbeat signal interval when the heartbeat signal whose interval variation is within a predetermined percentage is continuously measured a predetermined number of times or more is used as a reference value for heartbeat signal inspection. A heartbeat signal correction method, wherein a signal for complementation is inserted between signals at intervals of a predetermined percentage or more. 間隔のばらつきが所定%以内である心拍信号が所定回数以上連続して測定された時の心拍信号間隔の平均値を心拍信号検査用の基準値とし、測定された心拍信号のうち、基準値の所定%以下の間隔の信号はキャンセルし、次いで基準値の所定%以上の間隔の信号についてはその間に補完のための信号を挿入することを特徴とする心拍信号補正方法。The average value of the heartbeat signal interval when the heartbeat signal whose interval variation is within a predetermined percentage is continuously measured a predetermined number of times or more is used as a reference value for heartbeat signal inspection. A heart rate signal correction method comprising: canceling a signal having an interval of a predetermined percentage or less; 信号の挿入にあたり、信号を挿入することによってできるそれぞれの間隔が基準値に最も近くなるような本数及び位置に配置することを特徴とする請求項2または3記載の心拍信号補正方法。4. The heartbeat signal correcting method according to claim 2, wherein the number and position of each signal are such that each interval generated by inserting the signal is closest to the reference value.
JP2003178559A 2003-06-23 2003-06-23 Heart rate signal correction method Expired - Fee Related JP3726825B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2003178559A JP3726825B2 (en) 2003-06-23 2003-06-23 Heart rate signal correction method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2003178559A JP3726825B2 (en) 2003-06-23 2003-06-23 Heart rate signal correction method

Related Parent Applications (1)

Application Number Title Priority Date Filing Date
JP33724695A Division JP3460419B2 (en) 1995-12-25 1995-12-25 Derivation of reference value for heart rate signal test

Publications (2)

Publication Number Publication Date
JP2004000660A JP2004000660A (en) 2004-01-08
JP3726825B2 true JP3726825B2 (en) 2005-12-14

Family

ID=30438370

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2003178559A Expired - Fee Related JP3726825B2 (en) 2003-06-23 2003-06-23 Heart rate signal correction method

Country Status (1)

Country Link
JP (1) JP3726825B2 (en)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5315948B2 (en) * 2008-11-19 2013-10-16 富士通株式会社 Pulse rate calculation program, pulse rate calculation device, pulse rate calculation method, and portable terminal device
CN102449560B (en) 2009-04-26 2016-12-21 耐克创新有限合伙公司 Sports watch
US9141087B2 (en) 2009-04-26 2015-09-22 Nike, Inc. Athletic watch
KR101817587B1 (en) * 2012-01-04 2018-01-11 나이키 이노베이트 씨.브이. Athletic watch

Also Published As

Publication number Publication date
JP2004000660A (en) 2004-01-08

Similar Documents

Publication Publication Date Title
US20090112111A1 (en) Heart rate meter and method for removing noise of heart beat waveform
EP3576093A1 (en) Method and system for clinical effectiveness evaluation of artificial intelligence based medical device
RU2012101554A (en) ECG CONTROL WITH REDUCED FALSE ALARMS IN CONNECTION WITH ASYSTOLIA
CN114491383B (en) Abnormal data processing method and system for bridge monitoring
CN108324271A (en) Electrocardiosignal recognition methods, system and cardioelectric monitor equipment
JP3729143B2 (en) Pulse wave measuring device
KR102461702B1 (en) Processing method for electrocardiogram signal
JP3726825B2 (en) Heart rate signal correction method
JP4040096B2 (en) Device for detecting electromagnetic interference
JP2009125531A (en) Blood pressure measuring apparatus
JP3460419B2 (en) Derivation of reference value for heart rate signal test
JPH0326233A (en) Method and device for detecting qrs
US8005882B2 (en) System and method for determining application of adaptive filter
JP2002214185A (en) Detecting method and detecting device for sensor abnormality
CN110960203A (en) Cardiovascular characteristic parameter detection method and device
CN117426755A (en) Sleep monitoring method and device and intelligent equipment
CN111724386B (en) Data preprocessing method and system for multi-frequency imaging and multi-frequency imaging system
EP2938247B1 (en) Method and apparatus for reducing motion artifacts in ecg signals
CN110403593B (en) Heart rate detection method
JPS61266964A (en) Processor for period data on pulse signal
JPH06154342A (en) Pace maker pulse detecting circuit
CN111407275B (en) Method and device for monitoring skin condition in cupping device
WO2023032120A1 (en) Heartbeat detection method and heartbeat detection device
JP7160460B2 (en) Pulse rate calculation method, pulse rate calculation device, recording medium
CN118759316A (en) Active enhancement detection method and system for line fault characteristics of distribution box

Legal Events

Date Code Title Description
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: 20050906

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20050919

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20081007

Year of fee payment: 3

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20091007

Year of fee payment: 4

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20091007

Year of fee payment: 4

S533 Written request for registration of change of name

Free format text: JAPANESE INTERMEDIATE CODE: R313533

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20091007

Year of fee payment: 4

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20101007

Year of fee payment: 5

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20101007

Year of fee payment: 5

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20111007

Year of fee payment: 6

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20111007

Year of fee payment: 6

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20121007

Year of fee payment: 7

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20131007

Year of fee payment: 8

LAPS Cancellation because of no payment of annual fees