JP6098257B2 - 信号処理装置、信号処理方法及び信号処理プログラム - Google Patents
信号処理装置、信号処理方法及び信号処理プログラム Download PDFInfo
- Publication number
- JP6098257B2 JP6098257B2 JP2013054165A JP2013054165A JP6098257B2 JP 6098257 B2 JP6098257 B2 JP 6098257B2 JP 2013054165 A JP2013054165 A JP 2013054165A JP 2013054165 A JP2013054165 A JP 2013054165A JP 6098257 B2 JP6098257 B2 JP 6098257B2
- Authority
- JP
- Japan
- Prior art keywords
- window
- signal
- waveform
- pulse wave
- correlation coefficient
- 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.)
- Active
Links
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/7235—Details of waveform analysis
- A61B5/7246—Details of waveform analysis using correlation, e.g. template matching or determination of similarity
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/0059—Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence
- A61B5/0077—Devices for viewing the surface of the body, e.g. camera, magnifying lens
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/0059—Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence
- A61B5/0082—Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence adapted for particular medical purposes
-
- 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
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Veterinary Medicine (AREA)
- Public Health (AREA)
- Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Physiology (AREA)
- Cardiology (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Signal Processing (AREA)
- Psychiatry (AREA)
- Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
Description
まず、本実施例に係る脈波検出装置の機能的構成について説明する。図1は、実施例1に係る脈波検出装置の機能的構成を示すブロック図である。図1に示す脈波検出装置10は、太陽光や室内光などの一般の環境光の下で生体に計測器具を接触させずに、被験者が撮影された画像を用いて被験者の脈波、すなわち心臓の拍動に伴う血液の体積の変動を検出する脈波検出処理を実行するものである。
Ra=Rs+Rn・・・(2)
k=Gn/Rn・・・(3)
S=Gs+Gn−k(Rs+Rn)・・・(5)
S=Gs−(Gn/Rn)Rs・・・(6)
v(i-1)>=0 and v(i)<0・・・(8)
Rmin= 60/240・・・(13)
Rmax= 60/42・・・(14)
Rmax= min(60/42,(T(J)-T(J-2*m))/m/γ)・・・(16)
Rmin≦w≦Rmax・・・(17)
Rmax= min(60/42,U/n/γ)・・・(19)
Rmin≦w≦Rmax・・・(20)
x(2)=v(I-X-2)
…
x(X)=v(I)
y(z,2)=v(I-X-2-z)
…
y(z,X)=v(I-z)
次に、本実施例に係る脈波検出装置の処理の流れについて説明する。なお、ここでは、脈波検出装置10によって実行される(1)波形検出処理について説明した後に、(2)相関係数算出処理について説明することとする。
図12は、実施例1に係る波形検出処理の手順を示すフローチャートである。この波形検出処理は、カメラ11から画像が取得する度に処理を起動し、画像が取得されなくなるまで繰り返し実行される処理である。なお、図示しない入力デバイス等を介して中断操作を受け付けた場合には、波形検出処理を中止することもできる。
図13は、実施例1に係る相関係数算出処理の手順を示すフローチャートである。この相関係数算出処理は、波形検出部15から脈波波形が入力される限り、繰り返し実行される処理である。
上述してきたように、本実施例に係る脈波検出装置10は、波形の類似を相関係数を用いて判断する場合に、基準窓内の振幅値の分散と比較窓内の振幅値の分散との差が大きくなるにしたがって値が小さくなる相関係数を算出する信号処理を実行する。このため、本実施例に係る脈波検出装置10では、互いの窓内の波形が相似の関係にある場合でも相関係数が低く算出される。したがって、本実施例に係る脈波検出装置10によれば、波形の類似性の判定精度を向上させることができる。
上記の実施例1では、波形の類似性、すなわち脈波の周期性を判定するのに相関係数を用いる場合を例示したが、類似性を判断するための指標は相関係数に限定されない。例えば、脈波検出装置10は、下記の式(24)を用いて、ユークリッド距離を算出したり、下記の式(25)を用いて、正規化されたユークリッド距離を算出することもできる。この場合には、全てのずらし範囲の中で最短距離を持つものが最も周期性が高いずらし範囲Zであると判断できる。このように、ユークリッド距離または正規化されたユークリッド距離を用いる場合にも、最短距離が所定の閾値THよりも大きい場合には、測定条件が不安定であると判断し、当該基準窓内の波形は脈波の検出には使用しない。
上記の実施例1及び実施例2では、入力信号としてR信号およびG信号の二種類を用いる場合を例示したが、異なる複数の光波長成分を持つ信号であれば任意の種類の信号および任意の数の信号を入力信号とすることができる。例えば、R、G、B、IRおよびNIRなどの光波長成分が異なる信号のうち任意の組合せの信号を2つ用いることもできるし、また、3つ以上用いることもできる。
上記の実施例1では、脈波検出装置10が上記の相関係数算出処理をスタンドアローンで実行する場合を例示したが、クライアントサーバシステムとして実装することもできる。例えば、脈波検出装置10は、相関係数算出処理を実行するWebサーバとして実装することとしてもよいし、アウトソーシングによって相関係数算出サービスを始めとするサービスを提供するクラウドとして実装することとしてもかまわない。このように、脈波検出装置10がサーバ装置として動作する場合には、スマートフォンや携帯電話機等の携帯端末装置やパーソナルコンピュータ等の情報処理装置をクライアント端末として収容することができる。これらクライアント端末からネットワークを介して被験者の顔が映った画像が取得された場合に波形検出処理を実行し、その脈波波形を用いて相関係数算出処理を実行した上でその検出結果やその検出結果を用いてなされた診断結果をクライアント端末へ応答することによって脈波検出サービス及び診断サービスを提供できる。
また、図示した各装置の各構成要素は、必ずしも物理的に図示の如く構成されていることを要しない。すなわち、各装置の分散・統合の具体的形態は図示のものに限られず、その全部または一部を、各種の負荷や使用状況などに応じて、任意の単位で機能的または物理的に分散・統合して構成することができる。例えば、転換判定部16a、極値判定部16b、窓幅設定部17a、ずらし範囲設定部17bまたは相関係数算出部18を脈波検出装置10の外部装置としてネットワーク経由で接続するようにしてもよい。また、転換判定部16a、極値判定部16b、窓幅設定部17a、ずらし範囲設定部17bまたは相関係数算出部18を別の装置がそれぞれ有し、ネットワーク接続されて協働することで、上記の脈波検出装置10の機能を実現するようにしてもよい。
また、上記の実施例で説明した各種の処理は、予め用意されたプログラムをパーソナルコンピュータやワークステーションなどのコンピュータで実行することによって実現することができる。そこで、以下では、図14を用いて、上記の実施例と同様の機能を有する信号処理プログラムを実行するコンピュータの一例について説明する。
11 カメラ
13 取得部
14 抽出部
15 波形検出部
16a 転換判定部
16b 極値判定部
17a 窓幅設定部
17b ずらし範囲設定部
18 相関係数算出部
19 脈波検出部
Claims (8)
- 心拍または脈波に関する信号を検出する信号検出部と、
前記信号検出部によって検出された信号を用いて、所定の時間幅を持つ第1の窓に含まれる信号の波形と、前記第1の窓内の波形と比較する第2の窓に含まれる信号の波形との間で、前記第1の窓内の振幅値の分散と前記第2の窓内の振幅値の分散との差が大きくなるにしたがって値が小さくなる相関係数を算出する相関係数算出部と、
を有することを特徴とする信号処理装置。 - 前記信号検出部は、
生体が撮影された画像から脈波に対応する周波数帯の信号を検出し、
前記相関係数算出部は、
前記信号検出部によって検出された信号を用いて、前記第1の窓内の波形と前記第2の窓内の波形との間で前記相関係数を算出することを特徴とする請求項1に記載の信号処理装置。 - 前記第1の窓の時間幅を設定する窓幅設定部と、
前記第2の窓をずらす範囲を設定する範囲設定部とをさらに有し、
前記相関係数算出部は、
前記窓幅設定部によって時間幅が設定された第1の窓内の波形と、前記範囲設定部によって設定された範囲内で複数回にわたってずらした第2の窓内の波形の各々との間で前記相関係数を算出することを特徴とする請求項2に記載の信号処理装置。 - 前記信号の振幅値の正負が転換するか否かを判定する転換判定部をさらに有し、
前記窓幅設定部は、
前記転換判定部によって振幅値の正負が転換すると判定された転換点から前記脈波の周期の整数倍と対応する個数になるまで転換点を遡った区間の長さを前記第1の窓の時間幅として設定することを特徴とする請求項3に記載の信号処理装置。 - 前記信号の振幅値が極値であるか否かを判定する極値判定部をさらに有し、
前記範囲設定部は、
前記極値判定部によって極大であると判定された極大点のうち前記窓幅設定部によって時間幅が設定された第1の窓に含まれる極大点の数が前記脈波の周期と対応する個数と等しい場合に、前記範囲の設定を実行することを特徴とする請求項3または4に記載の信号処理装置。 - 前記範囲設定部は、前記窓幅設定部によって第1の窓に設定された時間幅を用いて、前記範囲を設定することを特徴とする請求項3〜5のいずれか1つに記載の信号処理装置。
- コンピュータが、
心拍または脈波に関する信号を検出し、
検出された信号を用いて、所定の時間幅を持つ第1の窓に含まれる信号の波形と、前記第1の窓内の波形と比較する第2の窓に含まれる信号の波形との間で、前記第1の窓内の振幅値の分散と前記第2の窓内の振幅値の分散との差が大きくなるにしたがって値が小さくなる相関係数を算出する
処理を実行することを特徴とする信号処理方法。 - コンピュータに、
心拍または脈波に関する信号を検出し、
検出された信号を用いて、所定の時間幅を持つ第1の窓に含まれる信号の波形と、前記第1の窓内の波形と比較する第2の窓に含まれる信号の波形との間で、前記第1の窓内の振幅値の分散と前記第2の窓内の振幅値の分散との差が大きくなるにしたがって値が小さくなる相関係数を算出する
処理を実行させることを特徴とする信号処理プログラム。
Priority Applications (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2013054165A JP6098257B2 (ja) | 2013-03-15 | 2013-03-15 | 信号処理装置、信号処理方法及び信号処理プログラム |
US14/099,457 US9962126B2 (en) | 2013-03-15 | 2013-12-06 | Signal processor, signal processing method, and recording medium |
EP13197941.1A EP2777485B1 (en) | 2013-03-15 | 2013-12-18 | Signal processor, signal processing method, and signal processing program |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2013054165A JP6098257B2 (ja) | 2013-03-15 | 2013-03-15 | 信号処理装置、信号処理方法及び信号処理プログラム |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2014176584A JP2014176584A (ja) | 2014-09-25 |
JP6098257B2 true JP6098257B2 (ja) | 2017-03-22 |
Family
ID=49918405
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2013054165A Active JP6098257B2 (ja) | 2013-03-15 | 2013-03-15 | 信号処理装置、信号処理方法及び信号処理プログラム |
Country Status (3)
Country | Link |
---|---|
US (1) | US9962126B2 (ja) |
EP (1) | EP2777485B1 (ja) |
JP (1) | JP6098257B2 (ja) |
Families Citing this family (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP6191517B2 (ja) * | 2014-03-14 | 2017-09-06 | 富士通株式会社 | 検出装置、検出プログラム、及び検出方法 |
US10397082B2 (en) * | 2014-08-07 | 2019-08-27 | Citrix Systems, Inc. | Internet infrastructure measurement method and system adapted to session volume |
US9582865B2 (en) * | 2014-10-20 | 2017-02-28 | Microsoft Technology Licensing, Llc | Visualization for blood flow in skin image data |
JP6488722B2 (ja) * | 2015-01-26 | 2019-03-27 | 富士通株式会社 | 脈波検出装置、脈波検出方法及び脈波検出プログラム |
JP6406171B2 (ja) * | 2015-08-25 | 2018-10-17 | トヨタ自動車株式会社 | 瞬目検知装置 |
JP6607259B2 (ja) * | 2015-11-20 | 2019-11-20 | 富士通株式会社 | 脈波分析装置、脈波分析方法、および脈波分析プログラム |
CN105520724A (zh) * | 2016-02-26 | 2016-04-27 | 严定远 | 一种测量人体心跳速率和呼吸频率的方法 |
JP2017158675A (ja) * | 2016-03-08 | 2017-09-14 | パナソニックIpマネジメント株式会社 | 脈拍推定装置、脈拍推定システムおよび脈拍推定方法 |
KR102588906B1 (ko) * | 2017-12-01 | 2023-10-13 | 삼성전자주식회사 | 생체 신호 품질 평가 장치 및 방법 |
US11457872B2 (en) | 2017-12-01 | 2022-10-04 | Samsung Electronics Co., Ltd. | Bio-signal quality assessment apparatus and bio-signal quality assessment method |
WO2019146025A1 (ja) * | 2018-01-24 | 2019-08-01 | 富士通株式会社 | 脈波算出装置、脈波算出方法及び脈波算出プログラム |
WO2020070808A1 (ja) * | 2018-10-02 | 2020-04-09 | 富士通株式会社 | 脈波算出装置、脈波算出方法及び脈波算出プログラム |
CN112836546A (zh) * | 2019-11-22 | 2021-05-25 | 深圳市理邦精密仪器股份有限公司 | 检测生理信号质量的方法、装置及电子设备 |
Family Cites Families (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3029912B2 (ja) * | 1992-01-08 | 2000-04-10 | 日本コーリン株式会社 | 脈拍数測定装置 |
JP3519185B2 (ja) | 1995-09-22 | 2004-04-12 | 花王株式会社 | 血流解析装置 |
FR2794961B1 (fr) * | 1999-06-16 | 2001-09-21 | Global Link Finance | Procede de determination du decalage temporel entre les instants de passage d'une meme onde de pouls en deux points de mesure distincts d'un reseau arteriel d'un etre vivant et d'estimation de sa pression aortique |
US9820658B2 (en) | 2006-06-30 | 2017-11-21 | Bao Q. Tran | Systems and methods for providing interoperability among healthcare devices |
EP1877774A4 (en) | 2005-03-25 | 2011-01-12 | Cnoga Holdings Ltd | OPTICAL SENSOR DEVICE AND IMAGE PROCESSING UNIT FOR MEASURING CHEMICAL CONCENTRATIONS, CHEMICAL SATURATIONS AND BIOPHYSICAL PARAMETERS |
US7725187B1 (en) * | 2006-11-06 | 2010-05-25 | Pacesetter, Inc. | Motion detection for sensors sensitive to motion induced noise |
JP4992043B2 (ja) * | 2007-08-13 | 2012-08-08 | 株式会社国際電気通信基礎技術研究所 | 行動識別装置、行動識別システムおよび行動識別方法 |
US20100152600A1 (en) | 2008-04-03 | 2010-06-17 | Kai Sensors, Inc. | Non-contact physiologic motion sensors and methods for use |
CA2720871A1 (en) * | 2008-04-03 | 2009-10-08 | Kai Medical, Inc. | Non-contact physiologic motion sensors and methods for use |
JP4609539B2 (ja) * | 2008-07-04 | 2011-01-12 | トヨタ自動車株式会社 | 眠気検出装置 |
JP5582478B2 (ja) * | 2009-07-16 | 2014-09-03 | 株式会社デルタツーリング | 生体状態推定装置及びコンピュータプログラム |
JP5234078B2 (ja) * | 2010-09-29 | 2013-07-10 | 株式会社デンソー | 脈波解析装置および血圧推定装置 |
US10856752B2 (en) * | 2010-12-28 | 2020-12-08 | Sotera Wireless, Inc. | Body-worn system for continuous, noninvasive measurement of cardiac output, stroke volume, cardiac power, and blood pressure |
JP5718126B2 (ja) * | 2011-03-31 | 2015-05-13 | 沖電気工業株式会社 | 微細振動特徴量算出装置、微細振動特徴量算出方法及びプログラム |
JP5673341B2 (ja) | 2011-05-13 | 2015-02-18 | 富士通株式会社 | 心拍計測装置、心拍計測方法及び心拍計測プログラム |
US20140073869A1 (en) * | 2012-09-11 | 2014-03-13 | Nellcor Puritan Bennett Llc | Methods and systems for determining physiological information using status indicators |
-
2013
- 2013-03-15 JP JP2013054165A patent/JP6098257B2/ja active Active
- 2013-12-06 US US14/099,457 patent/US9962126B2/en active Active
- 2013-12-18 EP EP13197941.1A patent/EP2777485B1/en active Active
Also Published As
Publication number | Publication date |
---|---|
EP2777485B1 (en) | 2017-10-25 |
EP2777485A1 (en) | 2014-09-17 |
JP2014176584A (ja) | 2014-09-25 |
US20140276114A1 (en) | 2014-09-18 |
US9962126B2 (en) | 2018-05-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP6098257B2 (ja) | 信号処理装置、信号処理方法及び信号処理プログラム | |
US10292602B2 (en) | Blood flow index calculating method, blood flow index calculating apparatus, and recording medium | |
JP5915757B2 (ja) | 脈波検出方法、脈波検出装置及び脈波検出プログラム | |
JP6256488B2 (ja) | 信号処理装置、信号処理方法及び信号処理プログラム | |
JP6167614B2 (ja) | 血流指標算出プログラム、血流指標算出装置および血流指標算出方法 | |
WO2016006027A1 (ja) | 脈波検出方法、脈波検出プログラム及び脈波検出装置 | |
US20150366456A1 (en) | Pulse wave velocity measurement method | |
JP6393984B2 (ja) | 脈拍計測装置、脈拍計測方法及び脈拍計測プログラム | |
JP6115263B2 (ja) | 脈波検出装置、脈波検出方法及び脈波検出プログラム | |
JP5920465B2 (ja) | バイタルサイン検出方法、バイタルサイン検出装置及びバイタルサイン検出プログラム | |
Feng et al. | Motion artifacts suppression for remote imaging photoplethysmography | |
US9986923B2 (en) | Selecting a region of interest for extracting physiological parameters from a video of a subject | |
JP6135255B2 (ja) | 心拍測定プログラム、心拍測定方法及び心拍測定装置 | |
Lagido et al. | Using the smartphone camera to monitor heart rate and rhythm in heart failure patients | |
Chwyl et al. | SAPPHIRE: Stochastically acquired photoplethysmogram for heart rate inference in realistic environments | |
JPWO2020054122A1 (ja) | 情報処理装置、プログラム及び情報処理方法 | |
JP2021146061A (ja) | 生体情報取得装置、生体情報取得方法及びプログラム | |
JP6020015B2 (ja) | 脈波検出装置、脈波検出プログラム及び脈波検出方法 | |
JP2016202603A (ja) | 生体情報処理システム、プログラム、生体情報処理システムの制御方法 | |
JP6488722B2 (ja) | 脈波検出装置、脈波検出方法及び脈波検出プログラム | |
JP6167615B2 (ja) | 血流指標算出プログラム、端末装置および血流指標算出方法 | |
JP7106821B2 (ja) | 生体信号解析装置、生体信号解析方法及びプログラム | |
JP6167849B2 (ja) | 脈波検出装置、脈波検出方法及び脈波検出プログラム | |
Desai et al. | Heart rate monitor using mobile camera | |
Suryasari et al. | Illuminance Color Independent in Remote Photoplethysmography for Pulse Rate Variability and Respiration Rate Measurement |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20151007 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20160921 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20160927 |
|
A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20161128 |
|
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: 20170124 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20170206 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 6098257 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |