JP2017144236A - 生体信号処理方法及びその装置並びにコンピュータ読出可能記録媒体 - Google Patents
生体信号処理方法及びその装置並びにコンピュータ読出可能記録媒体 Download PDFInfo
- Publication number
- JP2017144236A JP2017144236A JP2017010317A JP2017010317A JP2017144236A JP 2017144236 A JP2017144236 A JP 2017144236A JP 2017010317 A JP2017010317 A JP 2017010317A JP 2017010317 A JP2017010317 A JP 2017010317A JP 2017144236 A JP2017144236 A JP 2017144236A
- Authority
- JP
- Japan
- Prior art keywords
- reference point
- biological signal
- determining
- pulse direction
- point
- 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
- 238000003672 processing method Methods 0.000 title claims abstract description 33
- 238000000034 method Methods 0.000 claims abstract description 38
- 239000000284 extract Substances 0.000 claims description 15
- 230000008859 change Effects 0.000 claims description 14
- 238000001914 filtration Methods 0.000 claims description 2
- 238000012545 processing Methods 0.000 description 76
- 230000008569 process Effects 0.000 description 18
- 210000000707 wrist Anatomy 0.000 description 14
- 230000007935 neutral effect Effects 0.000 description 6
- 238000007781 pre-processing Methods 0.000 description 6
- 230000036772 blood pressure Effects 0.000 description 5
- 238000001514 detection method Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 5
- 230000007423 decrease Effects 0.000 description 4
- 230000036541 health Effects 0.000 description 4
- 238000002565 electrocardiography Methods 0.000 description 3
- 230000029058 respiratory gaseous exchange Effects 0.000 description 3
- 230000017531 blood circulation Effects 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 230000003862 health status Effects 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 230000002093 peripheral effect Effects 0.000 description 2
- 210000002321 radial artery Anatomy 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000001154 acute effect Effects 0.000 description 1
- 206010003119 arrhythmia Diseases 0.000 description 1
- 230000006793 arrhythmia Effects 0.000 description 1
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 description 1
- 210000000748 cardiovascular system Anatomy 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 230000002526 effect on cardiovascular system Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 210000003811 finger Anatomy 0.000 description 1
- 230000004927 fusion Effects 0.000 description 1
- 208000019622 heart disease Diseases 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 229910052760 oxygen Inorganic materials 0.000 description 1
- 239000001301 oxygen Substances 0.000 description 1
- 230000000241 respiratory effect Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 210000003813 thumb Anatomy 0.000 description 1
- 230000002792 vascular Effects 0.000 description 1
Images
Classifications
-
- 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/7271—Specific aspects of physiological measurement analysis
- A61B5/7282—Event detection, e.g. detecting unique waveforms indicative of a medical condition
-
- 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/02007—Evaluating blood vessel condition, e.g. elasticity, compliance
-
- 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/021—Measuring pressure in heart or blood vessels
- A61B5/02108—Measuring pressure in heart or blood vessels from analysis of pulse wave characteristics
-
- 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/02416—Detecting, measuring or recording pulse rate or heart rate using photoplethysmograph signals, e.g. generated by infrared radiation
-
- 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
-
- 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/6887—Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient mounted on external non-worn devices, e.g. non-medical devices
- A61B5/6898—Portable consumer electronic devices, e.g. music players, telephones, tablet computers
-
- 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
-
- 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/7235—Details of waveform analysis
-
- 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/7235—Details of waveform analysis
- A61B5/725—Details of waveform analysis using specific filters therefor, e.g. Kalman or adaptive filters
-
- 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/02438—Detecting, measuring or recording pulse rate or heart rate with portable devices, e.g. worn by the patient
-
- 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/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/6802—Sensor mounted on worn items
- A61B5/681—Wristwatch-type devices
-
- 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/7225—Details of analog processing, e.g. isolation amplifier, gain or sensitivity adjustment, filtering, baseline or drift compensation
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Molecular Biology (AREA)
- Animal Behavior & Ethology (AREA)
- Veterinary Medicine (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Public Health (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- General Health & Medical Sciences (AREA)
- Surgery (AREA)
- Physics & Mathematics (AREA)
- Physiology (AREA)
- Cardiology (AREA)
- Signal Processing (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Psychiatry (AREA)
- Vascular Medicine (AREA)
- Multimedia (AREA)
- Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
- Power Engineering (AREA)
Abstract
Description
特に、事務室及び家庭などのように日常生活内でユーザの健康状態をリアルタイムにモニタリングするモバイルヘルスケア分野に対する研究が盛んに行われている。
モバイルヘルスケアでは、時間及び空間の制約なしでユーザの生体信号を測定し、測定された生体信号を分析してユーザの健康状態を推定する。
例えば、モバイルヘルスケアでは、心臓疾患のユーザの心電図(Electrocardiography:ECG)をリアルタイムに測定及び分析して、不整脈が発生した時間又は不整脈の種類を推定することができ、これに関する情報をユーザに提供する。
また、本発明の他の目的は、モバイル環境ではリソースを少なく使用しながら生体信号の特徴点を早く抽出することのできる生体信号処理方法及びその装置を提供することにある。
前記候補基準点は、極大点に対応する複数の第1候補基準点及び極小点に対応する複数の第2候補基準点を含み、前記基準点を決定するステップは、互いに隣接する第1候補基準点で信号値が最も大きい第1候補基準点を第1基準点として決定するステップと、互いに隣接する第2候補基準点で信号値が最も小さい第2候補基準点を第2基準点として決定するステップと、を含むことが好ましい。
前記候補基準点を抽出するステップは、前記生体信号の波形の所定の時間区間内で決定された最も大きい信号値又は最も小さい信号値に基づいて前記候補基準点を決定するステップを含むことが好ましい。
前記パルス方向を決定するステップは、前記第1基準点及び第2基準点のそれぞれに対応する尖度を決定するステップと、前記尖度に基づいて前記パルス方向を決定するステップと、を含むことが好ましい。
前記尖度を決定するステップは、前記生体信号の波形の中の前記第1基準点から第1基準点に隣接する第1候補基準点の間の距離に基づいて前記第1基準点に対応する尖度を決定し、前記生体信号の波形の中の前記第2基準点から第2基準点に隣接する第2候補基準点の間の距離に基づいて前記第2基準点に対応する尖度を決定することが好ましい。
前記第1基準点に対応する尖度を決定するステップは、前記生体信号の波形から前記第1基準点と第1基準点に隣接する第1候補基準点とを連結する第1ラインと、前記第1基準点と第1基準点に隣接する他の第1候補基準点とを連結する第2ラインとが形成する角度に基づいて前記第1基準点に対応する尖度が決定されることが好ましい。
前記パルス方向を決定するステップは、前記第1基準点の尖度と前記第2基準点の尖度に基づいたパルス方向の強度に基づいて前記パルス方向を決定することが好ましい。
前記パルス方向を決定するステップは、現在の基準点が第1基準点であるか、あるいは第2基準点であるか否か、及び前記現在の基準点に対応する尖度に基づいて現在のパルス方向の強度を決定するステップと、次の基準点が第1基準点であるか、あるいは第2基準点であるか否か、及び前記次の基準点に対応する尖度と前記現在のパルス方向の強度とに基づいて次のパルス方向の強度を決定するステップと、を含むことが好ましい。
前記生体信号の特徴点を決定するステップは、前記生体信号のパルス方向が前記第1パルス方向であると決定された場合、前記第1パルス方向に対応する第1特徴点の決定方式により前記生体信号の特徴点を決定し、前記生体信号のパルス方向が前記第2パルス方向であると決定された場合、前記第2パルス方向に対応する第2特徴点の決定方式により前記生体信号の特徴点を決定することが好ましい。
前記パルス方向を決定するステップは、複数の第1基準点の尖度の和と複数の第2基準点の尖度の和との比較結果に基づいて前記パルス方向を決定することが好ましい。
前記生体信号から直流成分を除去するステップをさらに含み、前記基準点を抽出するステップは、前記直流成分が除去された生体信号の波形から基準点を抽出することが好ましい。
前記直流成分が除去された生体信号をフィルタリングして前記生体信号から高周波ノイズを除去するステップをさらに含み、前記基準点を抽出するステップは、前記高周波ノイズが除去された生体信号の波形から基準点を抽出すること。
前記生体信号の特徴点を決定するステップは、前記パルス方向に応じて第1基準点又は第2基準点のうちいずれか1つに基づいて前記生体信号の特徴点を決定することが好ましい。
前記生体信号の特徴点を決定するステップは、前記生体信号の波形に低域通過フィルタを適用するステップと、前記パルス方向に基づいて前記低域通過フィルタが適用された生体信号の波形に対して時間による傾きの和を決定するステップと、前記傾きの和と閾値に基づいて前記生体信号の特徴点を決定するステップと、を含むことが好ましい。
前記閾値は、前記生体信号の変化により適応的に変化することが好ましい。
前記候補基準点は、極大点に対応する第1候補基準点及び極小点に対応する第2候補基準点を含むことが好ましい。
また、モバイル環境ではリソースを少なく使用しながら生体信号の特徴点を早く抽出することができるという効果がある。
本実施形態で用いる用語は、単に特定の実施形態を説明するために用いるものであって、実施形態を限定しようとする意図はない。
単数の表現は、文脈上、明白に異なる意味をもたない限り複数の表現を含む。
本明細書において、「含む」又は「有する」等の用語は明細書上に記載した特徴、数字、ステップ、動作、構成要素、部品又はこれらを組み合わせたものが存在することを示すものであって、1つ又はそれ以上の他の特徴や数字、ステップ、動作、構成要素、部品、又はこれを組み合わせたものなどの存在又は付加の可能性を予め排除しないものとして理解しなければならない。
一般的に用いられる予め定義された用語は、関連技術の文脈上で有する意味と一致する意味を有するものと解釈すべきであって、本明細書で明白に定義しない限り、理想的又は過度に形式的な意味として解釈されることはない。
図面を参照して説明する際に、図面符号に関係なく同一の構成要素には同一の参照符号を付与し、それに対する重複説明を省略する。本実施形態の説明において関連する公知技術に対する具体的な説明が実施形態の要旨を不要に曖昧にすると判断される場合、その詳細な説明を省略する。
本明細書では、説明の便宜のために光電式容積脈波(Photoplethysmogram:以下、PPGと記す)から信号特徴を抽出する実施形態を中心に説明するが、実施形態の範囲がこれに限定されることはない。
例えば、実施形態は、心電図(Electrocardiography:ECG)、酸素飽和度、心弾(動)図(Ballistocardiogram:BCG)、脳電図(Electroencephalogram:EEG)、筋電図(Electromyogram:EMG)などの生体信号だけではなく、他のタイプの生体信号から信号特徴を抽出する過程にも適用され得る。
バンドタイプ又は腕時計タイプのウェアラブル装置100は、ユーザの手首に着用されてユーザの体から生体信号を測定して分析することができる。
ウェアラブル装置100は、非侵襲的方式(noninvasive manner)で生体信号を継続的にモニタリングする。
生体信号は、ユーザ身体に対する生体的情報を含む信号として、PPG、ECG、又はEMGなどの形態に測定される。
測定された生体信号に基づいてユーザの健康に関する様々な生体情報を推定することができる。
センサは、ウェアラブル装置100のストラップに位置して手首の橈骨動脈(radial artery)内の血流量の変化をPPG信号として測定する。
PPG信号は、心拍による血流量の変化情報を含む。
ウェアラブル装置100は、PPG信号の変化を分析してユーザの心拍数を算出し、算出された心拍数に基づいて心血管系情報を推定する。
又は、ウェアラブル装置100は、PPG信号に基づいて決定された心拍数、血圧などに基づいてユーザに運動コーチング情報を提供することもできる。
また、モバイル環境では少ないリソースを使用し、生体信号の特徴点を早く抽出することが要求される。
下記で説明する実施形態は、上記の要求事項を全て満たし得るソリューションを提供する。
PPG信号の測定は、モーションアーチファクトのようなノイズに敏感である。
センサを用いて測定されたPPG信号には、目的とする生体情報に関する信号成分だけではなく、呼吸活動及び身体の動きによるノイズ成分が含まれることがある。
図1に示すように、手首にウェアラブル装置100を着用したユーザが手首を動かす場合(例えば、手首を回転する場合)にPPG信号の逆転現象が生じ得る。
橈骨動脈上の皮膚表面とPPG信号を検出するためのセンサとの間の圧力が変化した時、又はセンサの検出位置が皮膚表面上で変更した時にPPG信号の逆転現象が発生する可能性がある。
PPG信号のパルス方向が上側に向かっている。
図2Bは、ユーザが手首を中立状態で回転させた時に測定されたPPG信号の波形を示す。
図2Bに示されたPPG信号の波形は、図2Aとは相違して、PPG信号のパルス方向が下側に向かっている。
このように、PPG信号の波形は手首の姿勢などにより逆転し、PPG信号から信号特徴を正確に抽出するためにはPPG信号のパルスがどちらの方に向かっているかを決定することが重要である。
以下で説明する実施形態によれば、生体信号から信号特徴の検出においてパルス方向を考慮することによって信号特徴をより正確に抽出するようにする。
本発明の一実施形態に係る生体信号処理方法は、図14に示す生体信号処理装置1400によって実行される。
図3及び図4に示した動作は、図に示すように順次に実行されたり又は説明された実施形態の範囲及び技術的な思想から離脱することなく、一部の動作が実行されないか、動作の順序が変更され得る。
図3及び図4に示した動作は、並列的又は同時に実行されてもよい。
また、上記で説明した図1及び図2に対する内容は、下記の図3及び図4に対する説明にも適用される。
生体信号処理装置は、生体信号から直流成分(DC component)を除去し、直流成分が除去された生体信号をフィルタリングして生体信号から高周波ノイズを除去する。
基準点は、測定された生体信号波形で生体情報を含むものとして考慮される。
一実施形態によれば、生体信号処理装置は、生体信号波形から候補基準点を一次的に抽出し、候補基準点のうち一定の基準を満たす基準点を決定する。
ピークは信号値が増加してから減少する点であり、トラフは信号値が減少してから増加する点である。
例えば、生体信号処理装置は、生体信号波形の所定の時間区間内で決定された最も大きい信号値に基づいて第1候補基準点を決定し、所定の時間区間内で決定された最も小さい信号値に基づいて第2候補基準点を決定する。
生体信号処理装置は、互いに隣接する第1候補基準点のうち信号値が最も大きい第1候補基準点を第1基準点として決定し、互いに隣接する第2候補基準点のうち信号値が最も小さい第2候補基準点を第2基準点として決定する。
例えば、生体信号処理装置は、生体信号のパルスが第1パルス方向(例えば、上側)及び第2パルス方向(例えば、下側)のいずれか一方に向かっているかを決定する。
パルス方向を決定するために、生体信号処理装置は、各基準点に対応する尖度(pulse sharpness)を用いることができる。
尖度は、基準点を中心に生体信号のパルスがどれ程急激に変化するか、又はどれ程尖っているかを示す。
生体信号処理装置は、第1基準点及び第2基準点のそれぞれに対応する尖度を決定し、決定された尖度に基づいてパルス方向を決定する。
他の実施形態によれば、生体信号処理装置は、第1基準点の尖度及び第2基準点の尖度に基づいてパルス方向の強度を算出し、パルス方向の強度を用いて生体信号のパルス方向を決定する。
これについては、図4を参照して以下で詳しく説明する。
生体信号処理装置は、生体信号波形で第1基準点から該当の第1基準点に隣接する第1候補基準点の間の距離に基づいて該当の第1基準点に対応する尖度を決定する。
ここで、尖度は、第1基準点と該当の第1基準点に隣接する第1候補基準点を連結する第1ラインと、該当の第1基準点と該当の第1基準点に隣接する他の第1候補基準点を連結する第2ラインとによって形成される角度によって決定される。
これと同様に、生体信号処理装置は、生体信号波形で第2基準点から該当の第2基準点に隣接する第2候補基準点の間の距離に基づいて該当の第2基準点に対応する尖度を決定する。
ここで、尖度は、第1基準点と該当の第1基準点に隣接する第2候補基準点を連結する第1ラインと、該当の第1基準点と該当の第1基準点に隣接する他の第2候補基準点を連結する第2ラインとによって形成される角度によって決定される。
これと同様に、生体信号処理装置は、生体信号波形で第2基準点から該当の第2基準点に隣接する第1候補基準点の間の距離に基づいて該当の第2基準点に対応する尖度を決定する。
一実施形態によれば、生体信号処理装置は、現在の基準点が第1基準点であるか又は第2基準点であるか否か、及び現在の基準点に対応する尖度に基づいて現在のパルス方向の強度を決定する。
その後、生体信号処理装置は、次の基準点が第1基準点であるか又は第2基準点であるか否か、及び次の基準点に対応する尖度及び現在のパルス方向の強度に基づいて次のパルス方向の強度を決定する。
このような工程により生体信号処理装置は、生体信号波形で第1基準点及び第2基準点が示されるごとにパルス方向の強度を継続的にアップデートする。
例えば、パルス方向の強度は、数式(4)に基づいて算出され、これについては後述にて説明する。
一実施形態によれば、生体信号処理装置は、パルス方向の強度が閾値(例えば、0)よりも大きい場合には生体信号波形が第1パルス方向であると決定し、パルス方向の強度が閾値以下である場合には生体信号波形が第1パルス方向と異なる第2パルス方向であると決定する。
例えば、生体信号処理装置は、生体信号のパルス方向が第1パルス方向であると決定された場合には、第1パルス方向に対応する第1特徴点の決定方式により生体信号の特徴点を決定し、生体信号のパルス方向が第2パルス方向であると決定された場合には、第2パルス方向に対応する第2特徴点の決定方式により生体信号の特徴点を決定する。
例えば、第1特徴点の決定方式と第2特徴点の決定方式は、全体的な特徴点の決定過程に類似するが、適用されるパラメータ値に差が存在したり、又は決定方式が全く異なることがある。
傾きの和は、数式(5)に基づいて算出され、これについては後述にて説明することにする。
生体信号処理装置は、傾きの和と閾値とを比較し、比較結果の傾きの和が閾値よりも大きい時点を識別する。
ここで、閾値は、生体信号の変化に応じてその大きさが適応的に変化し得る。
例えば、閾値の大きさは生体信号の大きさに比例して変化する。
そのため、生体信号の変化が大きいとしても該当の変化に適応的に閾値も変化するため、生体信号波形から特徴点をより正確に検出することができる。
生体信号処理装置は、識別された時点の周辺区間で生体信号の信号値が最も大きいか又は最も小さい信号値を探索し、探索された信号値を生体信号の特徴点として決定する。
例えば、パルス方向が上側方向にあると決定された場合、生体信号処理装置は、以前に決定された第1基準点のうち尖度が閾値よりも大きい第1基準点を特徴点として決定する。
反対に、パルス方向が下側であると決定された場合、生体信号処理装置は、以前に決定された第2基準点のうち尖度が閾値よりも大きい第2基準点を特徴点として決定する。
上記のような工程により、複雑な算出なしに生体信号の特徴点を迅速に決定でき、低周波のノイズを完全に除去しなくても生体信号の特徴点を正確に決定きる。
また、生体信号の変動及び逆転に強い特徴点を決定できる。
図5は、本発明の一実施形態に係る様々な状態で測定されたPPG信号波形の一例を示すグラフである。
PPG信号波形は、ユーザの状態に応じて多様に示され得る。
図5を参照すると、(a)は手首を中立位置にした時に測定されたPPG信号、(b)は手首を中立状態で回転した時に測定されたPPG信号、(c)は息を平常時より深くして手首を回転した位置から再び中立状態に戻す時に測定されたPPG信号、(d)は息を平常時より深くして手首を中立状態で回転させた時に測定されたPPG信号波形を示す。
(c)及び(d)は、PPG信号波形が(a)及び(b)に比べて相対的に激しく示される。
生体信号処理装置は、PPG信号が測定されれば、測定されたPPG信号から直流成分を除去する。
例えば、生体信号処理装置は、センサを用いて継続的にサンプリングされるPPG信号値の平均値を元(raw)のPPG信号から除去することでPPG信号から直流成分を除去することができる。
図6において、波形610は、各(a)、(b)、(c)、(d)の区間でセンサを用いて測定された元のPPG信号波形を示し、波形620は元のPPG信号から直流成分を除去した結果波形を示す。
直流成分を除去した後、生体信号処理装置は、低域通過フィルタを直流成分の除去されたPPG信号に適用して高周波ノイズを除去し得る。
例えば、生体信号処理装置は、10Hzのカットオフ周波数を有するハミングウィンドウ(Hamming window)を直流成分の除去されたPPG信号に適用する。
ハミングウィンドウによってPPG信号から10Hz以上の高周波成分が除去される。
生体信号処理装置は、測定されたPPG信号波形710から候補基準点として局部的なピーク及びトラフを抽出する。
ここで、ピークは第1候補基準点720に対応し、トラフは第2候補基準点730に対応する。
x(t)が第1候補基準点720又は第2候補基準点730である場合、p(t)はx(t)の値を有し、その他の場合には、p(t)が「0」の値を有する。
ただし、iの範囲は上記の数式(1)に定義された範囲に限定されることなく可変され得る。
このような工程により抽出された候補基準点のうちの一部は、PPG信号のパルス成分に対応しないことがある。
したがって、抽出された候補基準点で目的とするPPG信号のパルス成分に対応する基準点を追加的に決定する工程が求められ、これについては図8を参照して説明することにする。
生体信号処理装置は、PPG信号波形810から抽出された各候補基準点と隣接する他の候補基準点とを比較することによって不要な候補基準点を除外し、パルス方向を決定するための基準となる基準点を決定する。
実際のPPG信号のパルス成分に該当する第1候補基準点の信号値は、隣接する他の第1候補基準点の信号値よりも大きい。
また、生体信号処理装置は、隣接する第2候補基準点(840、845、850)のうち最も信号値の小さい第2候補基準点845を第2基準点として決定する。
実際のPPG信号のパルス成分に該当する第2候補基準点の信号値は、隣接する他の第2候補基準点の信号値よりも小さい。
図8には、上記のような工程により決定された第1基準点(830、855、860、865、870)と第2基準点(845、875、880、885、890)が示されている。
例えば、第1基準点830と第1基準点830に隣接する各第1候補基準点(825、835)を連結するラインが形成され、2つのラインが形成している角度に基づいて第1基準点830の尖度が算出される。
これと同様に、第2基準点845と第2基準点845に隣接する第2候補基準点(840、850)を連結するラインが形成され、2つのラインが形成している角度に基づいて第2基準点845の尖度が算出され得る。
尖度を算出する工程については、以下の図9を参照してより詳しく説明する。
一実施形態によれば、各基準点の尖度はコサイン法則を用いて算出され得る。
例えば、PPG信号波形910で決定された第1基準点920と、第1基準点920に隣接する第1候補基準点930を連結したラインb、第1基準点920と第1基準点920に隣接する他の第1候補基準点940を連結したラインc、及び第1候補基準点930と他の第1候補基準点940を連結したラインaによって形成された三角形を仮定する。
三角形を用いて第1基準点920の尖度が量的に算出され得る。
「−100」から「+100」の間の範囲を有するよう、角度θに100を乗算してθ値を整数化し、θ値を整数化することによって消費電力をより低減させる。
整数化されたθ値が「−100」に近いときにθは鈍角を示し、反対に、整数化されたθ値が「+100」に近いときはθが鋭角を示す。
整数化されたθ値に基づいて第1基準点920に対する尖度が決定される。
整数化されたθ値が大きいほど尖度は大きくなり、ここで、尖度が大きいことはパルスがより尖っていることを示す。
反対に、整数化されたθ値が小さいほど尖度が小さくなり、ここで、尖度が小さいことはパルスがより緩やかであることを示す。
そのため、算出量がより節減される。そして、求められたθ値に基づいて第1基準点920に対する尖度を決定することができる。
図10では、便宜上、各基準点の尖度に100を乗算して尖度が「−100」〜「+100」の間の値を有するようにした。
基準点の尖度が「+100」に近いほど基準点を含んでいるパルスの成分はより尖る。
図10に示したPPG信号波形1010は、PPG信号のパルス方向が上側方向にある場合である。
パルス方向が上側方向にある場合、第1基準点(1020、1025、1030、1035、1040)の尖度は、主に、第2基準点(1045、1050、1055、1060、1065)の尖度よりも大きい値を有する。
反対に、パルス方向が下側の場合、第2基準点(1045、1050、1055、1060、1065)の尖度は、主に、第1基準点(1020、1025、1030、1035、1040)の尖度よりも小さい値を有する。
第1尖度和が第2尖度和よりも大きい場合、生体信号処理装置はPPG信号波形1010のパルス方向が上側に向かうものと決定する。
反対に、第2尖度和が第1尖度和よりも大きい場合、生体信号処理装置はPPG信号波形1010のパルス方向が下側に向かうものと決定する。
パルス方向の強度は、PPG信号波形1010でパルスが上側又は下側のいずれかの方向にどれ程強く形成されたかを示す。
βが大きいほどdnewはdoldに影響をより少なく受ける。
数式(4)において、分母である(β+1)はdnew値を正規化するために用いられる。
sは現在の基準点が第1基準点であれば「+1」の値を有し、第2基準点であれば「−1」の値を有する。
生体信号処理装置は、dnewの値が「0」よりも大きい場合にはPPG信号波形1010のパルス方向が上側に向かうものと決定し、反対にdnew値が「0」よりも小さい場合にはPPG信号波形1010のパルス方向が下側に向かうものと決定する。
図11を参照すると、(a)区間で、PPG信号波形1110のパルス方向は上側に向かっているため、パルス方向の強度1120が「0」よりも大きい値を示す。
(b)区間で、PPG信号波形1110のパルス方向が下側に向うものとして逆転し、(a)区間におけるパルス方向の強度1120が少しずつ減少し、結局は「0」よりも小さくなることを示している。
(c)、(d)区間で、PPG信号波形1110は深い呼吸などにより激しく揺れていることを示す。
しかし、生体信号処理装置は、PPG信号波形1110の急激な揺動に関係なく、(c)、(d)区間でパルス方向の強度1120を正確に算出し、これによってPPG信号波形1110のパルス方向も正確に決定することができる。
PPG信号のパルス方向が決定された後に、パルス方向を考慮して生体信号(又は、図6を参照して説明した前処理工程が実行された生体信号)に対して特徴点が決定される。
一実施形態によれば、生体信号処理装置は、PPG信号波形に低域通過フィルタを適用し、生体信号のパルス方向に基づいてPPG信号波形に対する傾きの和を決定する。
Δykは「yk−yk−1」を示す。
傾きの和を用いることによって、別途の高域通過フィルタを使用しなくてもPPG信号波形から特徴点を正確に決定することができる。
ここで、閾値は、PPG信号の変化により適応的に変化し得る。
例えば、閾値が現在のPPG信号の信号値に比例して変化する。
Vprevは以前のパルス信号値を示し、stdPPGはPPG信号の標準偏差を示す。
Fsはサンプリング周波数を示す。
数式(6)によれば、PPG信号が閾値よりも大きくなるまで閾値は小さくなる。
そのため、呼吸や身体的な動きによってPPG信号の信号値が大きく変動する場合でも正確に特徴点を検出することができる。
PPG信号波形1210の変化により、数式(6)に基づいた閾値1230が適応的に決定され、閾値1230に基づいてPPG信号波形1210から特徴点1240を決定する。
例えば、傾きの和1220から閾値1230よりも大きい位置が識別され、PPG信号波形1210から該当の識別された位置を中心にする一定時間区間で最も大きい信号値を有する信号の位置が特徴点として決定される。
その後、図12と同様に、PPG信号波形1310の変化により数式(6)に基づいた閾値1330が適応的に決定され、閾値1330に基づいてPPG信号波形1310から特徴点1340が決定される。
例えば、傾きの和1320から閾値1330よりも大きい位置が識別され、PPG信号波形1310から該当の識別された位置を中心にする一定時間区間で最も小さい信号値を有する信号の位置が特徴点として決定される。
上記のような工程により、PPG信号波形のパルス方向がいずれの方向に向かってもPPG信号波形から特徴点を正確かつ迅速に検出することができる。
図14を参照すると、生体信号処理装置1400は、センサ1450からセンサ1450によって検知された生体信号を受信して、受信した生体信号から特徴点を検出する。
検出された特徴点に関する情報は、生体情報推定装置1440に伝えられ、生体情報推定装置1440は、特徴点に関する情報に基づいて血圧、心拍数などのような生体情報を推定する。
例えば、血圧及び心拍数のような生体情報はディスプレイ1430上に表示される。
プロセッサ1410は、図1〜図13を参照して前述した1つ以上の動作を行う。
例えば、プロセッサ1410は生体信号の波形から基準点を抽出し、抽出された基準点に基づいて生体信号のパルス方向を決定した後、パルス方向に対応する特徴点の決定方式により生体信号の特徴点を決定する。
プロセッサ1410は、基準点のそれぞれに対応する尖度を決定し、尖度に基づいて生体信号の特徴点を決定する。
このようなプロセッサ1410は、複数の論理ゲートのアレイに実現され、他の形態のハードウェアでも具現できることは本実施形態が属する技術分野で通常の知識を有する者であれば理解できるのであろう。
一部の実施形態において、メモリ1420は、非一時的なコンピュータで読み取り可能な記録媒体、例えば、高速ランダムアクセスメモリ及び/又は不揮発性コンピュータで読み取り可能な記録媒体(例えば、1つ以上のディスク記録デバイス、フラッシュメモリデバイス、又は、その他の非揮発性固体メモリデバイス)を含み得る。
ディスプレイ1430は、ユーザ入力を受信したりユーザインタフェースを実現する1つ以上のハードウェアコンポーネントを含み得る。
ディスプレイ1430は、ディスプレイ領域、ジェスチャキャプチャー領域、タッチ検出ディスプレイ、及び構成可能な領域の内のいずれかの組み合わせを含み得る。
ディスプレイ1430は、生体信号処理装置1400に内蔵されたり、又は生体信号処理装置1400から取り外され外部周辺機器であり得る。
例えば、本実施形態で説明した装置及び構成要素は、例えば、プロセッサ、コントローラ、ALU(arithmetic logic unit)、デジタル信号プロセッサ(digital signal processor)、マイクロコンピュータ、FPA(field programmable array)、PLU(programmable logic unit)、マイクロプロセッサ、又は命令(instruction)を実行して応答する異なる装置のように、1つ以上の汎用コンピュータ又は特殊目的コンピュータを用いて具現される。
また、処理装置は、ソフトウェアの実行に応答してデータをアクセス、格納、操作、処理、及び生成する。
理解の便宜のために、処理装置は1つが使用されるものとして説明する場合もあるが、当該技術分野で通常の知識を有する者は、処理装置が複数の処理要素(processing element)及び/又は複数類型の処理要素を含むことが分かる。
例えば、処理装置は、複数のプロセッサ又は1つのプロセッサ及び1つのコントローラを含む。
また、並列プロセッサ(parallel processor)のような、他の処理構成も可能である。
ソフトウェア及び/又はデータは、処理装置によって解釈され、処理装置に命令又はデータを提供するためのあらゆる類型の機械、構成要素、物理的装置、仮想装置、コンピュータ格納媒体又は装置、或いは送信される信号波を介して永久的又は一時的に具現化される。
ソフトウェアは、ネットワークに接続されたコンピュータシステム上に分散され、分散された方法で格納されるか又は実行され得る。
ソフトウェア及びデータは、1つ以上のコンピュータ読み取り可能な記録媒体に格納される。
記録媒体は、プログラム命令、データファイル、データ構造などを単独又は組み合わせて含む。
記録媒体及びプログラム命令は、本発明の目的のために特別に設計して構成されたものでもよく、コンピュータソフトウェア分野の技術を有する当業者にとって公知のものであり使用可能なものであってもよい。
ハードウェア装置は、本発明の動作を実行するために1つ以上のソフトウェアモジュールとして作動するように構成してもよく、その逆も同様である。
1400 生体信号処理装置
1410 プロセッサ
1420 メモリ
1430 ディスプレイ
1440 生体情報推定装置
1450 センサ
Claims (21)
- 生体信号の波形から基準点を抽出するステップと、
前記抽出された基準点に基づいて前記生体信号のパルス方向を決定するステップと、
前記パルス方向に対応する特徴点の決定方式により前記生体信号の特徴点を決定するステップと、を有することを特徴とする生体信号処理方法。 - 前記基準点を抽出するステップは、前記生体信号の波形から極大点及び極小点のうち少なくとも1つに対応する複数の候補基準点を抽出するステップと、
前記複数の候補基準点の各々の信号値に基づいて前記候補基準点の内から前記基準点を決定するステップと、を含むことを特徴とする請求項1に記載の生体信号処理方法。 - 前記候補基準点は、極大点に対応する複数の第1候補基準点及び極小点に対応する複数の第2候補基準点を含み、
前記基準点を決定するステップは、互いに隣接する第1候補基準点で信号値が最も大きい第1候補基準点を第1基準点として決定するステップと、
互いに隣接する第2候補基準点で信号値が最も小さい第2候補基準点を第2基準点として決定するステップと、を含むことを特徴とする請求項2に記載の生体信号処理方法。 - 前記候補基準点を抽出するステップは、前記生体信号の波形の所定の時間区間内で決定された最も大きい信号値又は最も小さい信号値に基づいて前記候補基準点を決定するステップを含むことを特徴とする請求項2又は3に記載の生体信号処理方法。
- 前記パルス方向を決定するステップは、前記第1基準点及び第2基準点のそれぞれに対応する尖度を決定するステップと、
前記尖度に基づいて前記パルス方向を決定するステップと、を含むことを特徴とする請求項3又は4に記載の生体信号処理方法。 - 前記尖度を決定するステップは、前記生体信号の波形の中の前記第1基準点から第1基準点に隣接する第1候補基準点の間の距離に基づいて前記第1基準点に対応する尖度を決定し、
前記生体信号の波形の中の前記第2基準点から第2基準点に隣接する第2候補基準点の間の距離に基づいて前記第2基準点に対応する尖度を決定することを特徴とする請求項5に記載の生体信号処理方法。 - 前記第1基準点に対応する尖度を決定するステップは、前記生体信号の波形から前記第1基準点と第1基準点に隣接する第1候補基準点とを連結する第1ラインと、前記第1基準点と第1基準点に隣接する他の第1候補基準点とを連結する第2ラインとが形成する角度に基づいて前記第1基準点に対応する尖度が決定されることを特徴とする請求項6に記載の生体信号処理方法。
- 前記パルス方向を決定するステップは、前記第1基準点の尖度と前記第2基準点の尖度に基づいたパルス方向の強度に基づいて前記パルス方向を決定することを特徴とする請求項5乃至7のいずれか一項に記載の生体信号処理方法。
- 前記パルス方向を決定するステップは、現在の基準点が第1基準点であるか、あるいは第2基準点であるか否か、及び前記現在の基準点に対応する尖度に基づいて現在のパルス方向の強度を決定するステップと、
次の基準点が第1基準点であるか、あるいは第2基準点であるか否か、及び前記次の基準点に対応する尖度と前記現在のパルス方向の強度とに基づいて次のパルス方向の強度を決定するステップと、を含むことを特徴とする請求項8に記載の生体信号処理方法。 - 前記パルス方向を決定するステップは、前記パルス方向の強度が閾値よりも大きい場合、前記生体信号の波形が第1パルス方向であると決定し、
前記パルス方向の強度が前記閾値以下である場合、前記生体信号の波形が第2パルス方向であると決定することを特徴とする請求項8又は9に記載の生体信号処理方法。 - 前記生体信号の特徴点を決定するステップは、前記生体信号のパルス方向が前記第1パルス方向であると決定された場合、前記第1パルス方向に対応する第1特徴点の決定方式により前記生体信号の特徴点を決定し、
前記生体信号のパルス方向が前記第2パルス方向であると決定された場合、前記第2パルス方向に対応する第2特徴点の決定方式により前記生体信号の特徴点を決定することを特徴とする請求項10に記載の生体信号処理方法。 - 前記パルス方向を決定するステップは、複数の第1基準点の尖度の和と複数の第2基準点の尖度の和との比較結果に基づいて前記パルス方向を決定することを特徴とする請求項5乃至11のいずれか一項に記載の生体信号処理方法。
- 前記生体信号から直流成分を除去するステップをさらに含み、
前記基準点を抽出するステップは、前記直流成分が除去された生体信号の波形から基準点を抽出することを特徴とする請求項1乃至12のいずれか一項に記載の生体信号処理方法。 - 前記直流成分が除去された生体信号をフィルタリングして前記生体信号から高周波ノイズを除去するステップをさらに含み、
前記基準点を抽出するステップは、前記高周波ノイズが除去された生体信号の波形から基準点を抽出することを特徴とする請求項13に記載の生体信号処理方法。 - 前記生体信号の特徴点を決定するステップは、前記パルス方向に応じて第1基準点又は第2基準点のうちいずれか1つに基づいて前記生体信号の特徴点を決定することを特徴とする請求項3乃至14のいずれか一項に記載の生体信号処理方法。
- 前記生体信号の特徴点を決定するステップは、前記生体信号の波形に低域通過フィルタを適用するステップと、
前記パルス方向に基づいて前記低域通過フィルタが適用された生体信号の波形に対して時間による傾きの和を決定するステップと、
前記傾きの和と閾値に基づいて前記生体信号の特徴点を決定するステップと、を含むことを特徴とする請求項1乃至15のいずれか一項に記載の生体信号処理方法。 - 前記閾値は、前記生体信号の変化により適応的に変化することを特徴とする請求項16に記載の生体信号処理方法。
- 前記候補基準点は、極大点に対応する第1候補基準点及び極小点に対応する第2候補基準点を含むことを特徴とする請求項2に記載の生体信号処理方法。
- 請求項1乃至請求項18のいずれか一項に記載の生体信号処理方法を実行させるためのプログラムが記録されたことを特徴とするコンピュータ読出可能記録媒体。
- 少なくとも1つのプロセッサを含み、
前記プロセッサは、生体信号の波形から基準点を抽出する動作と、
前記抽出された基準点に基づいて前記生体信号のパルス方向を決定する動作と、
前記パルス方向に対応する特徴点の決定方式により前記生体信号の特徴点を決定する動作と、を実行することを特徴とする生体信号処理装置。 - 少なくとも1つのプロセッサを含み、
前記プロセッサは、
生体信号の波形から基準点を抽出する動作と、
前記基準点のそれぞれに対応する尖度を決定する動作と、
前記尖度に基づいて前記生体信号の特徴点を決定する動作と、を実行し、
前記尖度は、前記基準点を中心に前記生体信号の波形がその程度急激に変化するかを示すものであることを特徴とする生体信号処理装置。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
KR10-2016-0016901 | 2016-02-15 | ||
KR1020160016901A KR102551184B1 (ko) | 2016-02-15 | 2016-02-15 | 생체신호 처리 방법 및 생체신호 처리 장치 |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2017144236A true JP2017144236A (ja) | 2017-08-24 |
JP6842306B2 JP6842306B2 (ja) | 2021-03-17 |
Family
ID=57794047
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2017010317A Active JP6842306B2 (ja) | 2016-02-15 | 2017-01-24 | 生体信号処理方法及びその装置並びにコンピュータ読出可能記録媒体 |
Country Status (5)
Country | Link |
---|---|
US (2) | US10335093B2 (ja) |
EP (1) | EP3205262B1 (ja) |
JP (1) | JP6842306B2 (ja) |
KR (1) | KR102551184B1 (ja) |
CN (1) | CN107080526B (ja) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20210028414A (ko) * | 2019-09-04 | 2021-03-12 | 아주대학교산학협력단 | 데이터 획득 및 데이터 변환 방법 및 이를 위한 장치 |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2020029031A1 (zh) * | 2018-08-06 | 2020-02-13 | 高维度(深圳)生物信息智能应用有限公司 | 一种信号处理方法、系统及计算机存储介质 |
US11445927B2 (en) | 2019-02-13 | 2022-09-20 | Viavi Solutions Inc. | Baseline correction and extraction of heartbeat profiles |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS6257538A (ja) * | 1985-05-16 | 1987-03-13 | テルモ株式会社 | 血圧測定装置 |
JPH05184548A (ja) * | 1992-01-08 | 1993-07-27 | Nippon Colin Co Ltd | 脈拍数測定装置 |
US5243992A (en) * | 1990-03-30 | 1993-09-14 | Colin Electronics Co., Ltd. | Pulse rate sensor system |
JP2002010987A (ja) * | 2000-06-29 | 2002-01-15 | Takeshi Matsumura | 生体情報信号処理方法と装置およびそのプログラムを記憶した記憶媒体 |
JP2008061824A (ja) * | 2006-09-07 | 2008-03-21 | Omron Healthcare Co Ltd | 医療用測定器、生体信号波形抽出方法、および生体信号波形抽出プログラム |
JP2008212391A (ja) * | 2007-03-05 | 2008-09-18 | Isuzu Motors Ltd | ストレス検出装置 |
JP2009189562A (ja) * | 2008-02-14 | 2009-08-27 | Citizen Holdings Co Ltd | 脈波測定装置 |
US9042973B2 (en) * | 2011-03-24 | 2015-05-26 | Draeger Medical Systems, Inc. | Apparatus and method for measuring physiological signal quality |
US20160022152A1 (en) * | 2014-07-28 | 2016-01-28 | Murata Manufacturing Co., Ltd. | Method and system for monitoring stress |
Family Cites Families (24)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3965444A (en) * | 1975-01-03 | 1976-06-22 | Raytheon Company | Temperature compensated surface acoustic wave devices |
JP3487829B2 (ja) * | 2001-02-22 | 2004-01-19 | 晴子 高田 | 血管老化の評価装置 |
JP2005080712A (ja) | 2003-09-04 | 2005-03-31 | Medical Bridge Kk | 心臓健康指数の算出方法及び心電図特定波分類方法 |
JP4696574B2 (ja) * | 2005-02-02 | 2011-06-08 | 株式会社ニコン | 半導体検査装置 |
JP4469746B2 (ja) | 2005-03-29 | 2010-05-26 | 株式会社東芝 | 心拍計測装置及び心拍計測装置の作動方法 |
JP2006280809A (ja) | 2005-04-04 | 2006-10-19 | Shibuya Kogyo Co Ltd | 血管狭窄検査装置 |
JP5234078B2 (ja) | 2010-09-29 | 2013-07-10 | 株式会社デンソー | 脈波解析装置および血圧推定装置 |
JP5516428B2 (ja) | 2010-10-14 | 2014-06-11 | 株式会社村田製作所 | 拍動周期算出装置およびこれを備えた生体センサ |
JP2012139433A (ja) | 2011-01-04 | 2012-07-26 | Sony Corp | 検出装置、検出方法およびプログラム |
JP5605253B2 (ja) | 2011-02-09 | 2014-10-15 | セイコーエプソン株式会社 | 脈波計測装置、およびプログラム |
KR20120113530A (ko) | 2011-04-05 | 2012-10-15 | 삼성전자주식회사 | 생체 신호의 피크 검출 방법 및 장치 |
JP2012242338A (ja) * | 2011-05-23 | 2012-12-10 | Toshiba Corp | 診断用パルス信号を備える制御システム、及びその制御装置 |
KR20130072728A (ko) | 2011-12-22 | 2013-07-02 | 삼성전자주식회사 | R-피크 검출 장치 및 그 제어 방법 |
KR101308609B1 (ko) | 2012-01-20 | 2013-09-23 | 경북대학교 산학협력단 | 적응적 국부 임계화를 이용한 심전도 신호의 r 피크 검출방법 |
US9936889B2 (en) | 2012-03-08 | 2018-04-10 | Samsung Electronics Co., Ltd. | Apparatus and method of controlling threshold for detecting peaks of physiological signals |
CN104507384A (zh) | 2012-07-30 | 2015-04-08 | 三菱化学控股株式会社 | 检体信息检测单元、检体信息处理装置、电动牙刷装置、电动剃须刀装置、检体信息检测装置、老龄化度评价方法及老龄化度评价装置 |
JP5526298B1 (ja) | 2012-09-10 | 2014-06-18 | パナソニック株式会社 | 電子機器、情報処理装置、情報処理方法およびプログラム |
JP6251997B2 (ja) * | 2012-09-18 | 2017-12-27 | カシオ計算機株式会社 | 脈拍データ検出装置、脈拍データ検出方法、および脈拍データ検出プログラム |
US8954140B2 (en) | 2012-09-27 | 2015-02-10 | Samsung Electronics Co., Ltd. | Method and system for determining QRS complexes in electrocardiogram signals |
KR101409800B1 (ko) | 2012-11-08 | 2014-06-24 | 상지대학교산학협력단 | 맥파로 얻어진 맥파형의 주파수 분석을 통한 호흡성분 추출장치 및 방법 |
KR102023991B1 (ko) * | 2012-11-12 | 2019-09-23 | 삼성전자주식회사 | 생체 신호 전송기, 생체 신호 수신기, 및 생체 신호 송수신 방법 |
FR3007553B1 (fr) | 2013-06-25 | 2015-07-17 | Morpho | Procede de detection d'un vrai visage |
KR102168657B1 (ko) | 2013-07-25 | 2020-10-22 | 삼성전자주식회사 | R 피크 검출을 위한 심전도 장치 및 방법 |
US20150359447A1 (en) * | 2014-06-17 | 2015-12-17 | Kabushiki Kaisha Toshiba | Electronic device, method, and computer program product |
-
2016
- 2016-02-15 KR KR1020160016901A patent/KR102551184B1/ko active IP Right Grant
- 2016-09-21 US US15/271,255 patent/US10335093B2/en active Active
- 2016-11-14 CN CN201611019661.4A patent/CN107080526B/zh active Active
- 2016-12-19 EP EP16205005.8A patent/EP3205262B1/en active Active
-
2017
- 2017-01-24 JP JP2017010317A patent/JP6842306B2/ja active Active
-
2019
- 2019-05-22 US US16/419,297 patent/US11911189B2/en active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS6257538A (ja) * | 1985-05-16 | 1987-03-13 | テルモ株式会社 | 血圧測定装置 |
US5243992A (en) * | 1990-03-30 | 1993-09-14 | Colin Electronics Co., Ltd. | Pulse rate sensor system |
JPH05184548A (ja) * | 1992-01-08 | 1993-07-27 | Nippon Colin Co Ltd | 脈拍数測定装置 |
JP2002010987A (ja) * | 2000-06-29 | 2002-01-15 | Takeshi Matsumura | 生体情報信号処理方法と装置およびそのプログラムを記憶した記憶媒体 |
JP2008061824A (ja) * | 2006-09-07 | 2008-03-21 | Omron Healthcare Co Ltd | 医療用測定器、生体信号波形抽出方法、および生体信号波形抽出プログラム |
JP2008212391A (ja) * | 2007-03-05 | 2008-09-18 | Isuzu Motors Ltd | ストレス検出装置 |
JP2009189562A (ja) * | 2008-02-14 | 2009-08-27 | Citizen Holdings Co Ltd | 脈波測定装置 |
US9042973B2 (en) * | 2011-03-24 | 2015-05-26 | Draeger Medical Systems, Inc. | Apparatus and method for measuring physiological signal quality |
US20160022152A1 (en) * | 2014-07-28 | 2016-01-28 | Murata Manufacturing Co., Ltd. | Method and system for monitoring stress |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20210028414A (ko) * | 2019-09-04 | 2021-03-12 | 아주대학교산학협력단 | 데이터 획득 및 데이터 변환 방법 및 이를 위한 장치 |
KR102353390B1 (ko) * | 2019-09-04 | 2022-01-19 | 아주대학교 산학협력단 | 데이터 획득 및 데이터 변환 방법 및 이를 위한 장치 |
Also Published As
Publication number | Publication date |
---|---|
US20170231579A1 (en) | 2017-08-17 |
US11911189B2 (en) | 2024-02-27 |
US20190269373A1 (en) | 2019-09-05 |
CN107080526B (zh) | 2021-04-02 |
US10335093B2 (en) | 2019-07-02 |
EP3205262B1 (en) | 2019-08-28 |
KR102551184B1 (ko) | 2023-07-04 |
EP3205262A1 (en) | 2017-08-16 |
JP6842306B2 (ja) | 2021-03-17 |
KR20170095483A (ko) | 2017-08-23 |
CN107080526A (zh) | 2017-08-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
TWI720215B (zh) | 提供即時訊號分段和基準點對準架構的系統與方法 | |
Nathan et al. | Particle filtering and sensor fusion for robust heart rate monitoring using wearable sensors | |
US9681814B2 (en) | Devices and methods for real-time denoising of electrocardiograms | |
EP3379997B1 (en) | Method to quantify photoplethysmogram (ppg) signal quality | |
Castells-Rufas et al. | Simple real-time QRS detector with the MaMeMi filter | |
EP3110313B1 (en) | Method and apparatus for generating assessments using physical activity and biometric parameters | |
JP6190466B2 (ja) | 生体信号測定器及び接触状態推定方法 | |
KR20170113252A (ko) | 저전력 센서를 이용한 심박수 및 호흡수 추정을 위한 방법 및 장치 | |
TWI702937B (zh) | 從具雜訊之心電圖資料決定心跳率的設備、系統和方法 | |
WO2017118127A1 (zh) | 一种心跳信号处理方法、装置和系统 | |
WO2012158984A1 (en) | Methods and apparatus for assessment of atypical brain activity | |
CN112165897A (zh) | 光体积描记图数据分析与呈现 | |
JP2016531629A (ja) | 生理学的信号を処理するデバイス、方法及びシステム | |
JP6869691B2 (ja) | 信号特徴抽出方法及びその装置並びにコンピュータ読取可能記録媒体 | |
JP6842306B2 (ja) | 生体信号処理方法及びその装置並びにコンピュータ読出可能記録媒体 | |
Qiao et al. | Revise: Remote vital signs measurement using smartphone camera | |
KR20140089897A (ko) | 저전력 생체신호 모니터링 방법 및 장치 | |
KR102445561B1 (ko) | 심박수 추정을 위한 데이터 수집 및 심박 추정 알고리즘 및 이를 위한 장치 | |
Anwer et al. | Evaluation of data processing and artifact removal approaches used for physiological signals captured using wearable sensing devices during construction tasks | |
WO2016053444A1 (en) | System and method for motion artifact reduction using surface electromyography | |
KR20230081154A (ko) | 심박주기 감지를 위한 데이터 처리 장치 및 그 동작 방법 | |
WO2020012807A1 (ja) | 生体情報処理装置及び生体情報処理方法 | |
Oweis et al. | Heart sounds segmentation utilizing teager energy operator | |
JP7327417B2 (ja) | 状態推定装置、状態推定方法、及びプログラム | |
Kher et al. | RPCA-based detection and quantification of motion artifacts in ECG signals |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20191002 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20200630 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20200707 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20201006 |
|
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: 20210216 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20210219 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 6842306 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |