JP6052005B2 - 脈波検出装置、脈波検出方法及び脈波検出プログラム - Google Patents
脈波検出装置、脈波検出方法及び脈波検出プログラム Download PDFInfo
- Publication number
- JP6052005B2 JP6052005B2 JP2013067667A JP2013067667A JP6052005B2 JP 6052005 B2 JP6052005 B2 JP 6052005B2 JP 2013067667 A JP2013067667 A JP 2013067667A JP 2013067667 A JP2013067667 A JP 2013067667A JP 6052005 B2 JP6052005 B2 JP 6052005B2
- Authority
- JP
- Japan
- Prior art keywords
- pulse wave
- region
- image
- wave detection
- regions
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
- 238000001514 detection method Methods 0.000 title claims description 101
- 230000008859 change Effects 0.000 claims description 93
- 238000012545 processing Methods 0.000 claims description 74
- 238000004364 calculation method Methods 0.000 claims description 48
- 238000000034 method Methods 0.000 claims description 38
- 230000036772 blood pressure Effects 0.000 claims description 27
- 230000008569 process Effects 0.000 claims description 24
- 238000005259 measurement Methods 0.000 claims description 17
- 238000003384 imaging method Methods 0.000 claims description 10
- 210000001367 artery Anatomy 0.000 claims description 2
- 238000005070 sampling Methods 0.000 description 84
- 238000010586 diagram Methods 0.000 description 21
- 230000000630 rising effect Effects 0.000 description 13
- 230000006870 function Effects 0.000 description 8
- 230000035488 systolic blood pressure Effects 0.000 description 8
- 230000000644 propagated effect Effects 0.000 description 6
- 238000004891 communication Methods 0.000 description 5
- 210000004204 blood vessel Anatomy 0.000 description 4
- 239000000284 extract Substances 0.000 description 4
- 210000003128 head Anatomy 0.000 description 4
- 239000008280 blood Substances 0.000 description 3
- 210000004369 blood Anatomy 0.000 description 3
- 238000003745 diagnosis Methods 0.000 description 3
- 230000001815 facial effect Effects 0.000 description 3
- 230000003287 optical effect Effects 0.000 description 3
- 239000004065 semiconductor Substances 0.000 description 3
- 210000000467 autonomic pathway Anatomy 0.000 description 2
- 230000004397 blinking Effects 0.000 description 2
- 238000012880 independent component analysis Methods 0.000 description 2
- 230000010354 integration Effects 0.000 description 2
- 239000010454 slate Substances 0.000 description 2
- 206010003210 Arteriosclerosis Diseases 0.000 description 1
- 238000012935 Averaging Methods 0.000 description 1
- 208000011775 arteriosclerosis disease Diseases 0.000 description 1
- 238000009530 blood pressure measurement Methods 0.000 description 1
- 230000000295 complement effect Effects 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 230000008094 contradictory effect Effects 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 230000003111 delayed effect Effects 0.000 description 1
- 210000005069 ears Anatomy 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 208000019622 heart disease Diseases 0.000 description 1
- 238000005286 illumination Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000010365 information processing Effects 0.000 description 1
- 229910044991 metal oxide Inorganic materials 0.000 description 1
- 150000004706 metal oxides Chemical class 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000012946 outsourcing Methods 0.000 description 1
- 210000005037 parasympathetic nerve Anatomy 0.000 description 1
- 230000001902 propagating effect Effects 0.000 description 1
- 230000010349 pulsation Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 230000002889 sympathetic effect Effects 0.000 description 1
Images
Landscapes
- Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
- Measuring And Recording Apparatus For Diagnosis (AREA)
Description
まず、本実施例に係る脈波検出装置の機能的構成について説明する。図1は、実施例1に係る脈波検出装置の機能的構成を示すブロック図である。図1に示す脈波検出装置10は、太陽光や室内光などの一般の環境光の下で生体に計測器具を接触させずに、被験者が撮影された画像を用いて被験者の脈波、すなわち心臓の拍動に伴う血液の体積の変動を検出する脈波検出処理を実行するものである。
続いて、本実施例に係る脈波検出装置の処理の流れについて説明する。図11は、実施例1に係る遅延量算出処理の手順を示すフローチャートである。この遅延量算出処理は、カメラ11によって画像が画像メモリ11aに保存される度に処理を起動し、画像メモリ11aに画像が保存されなくなるまで繰り返し実行される処理である。なお、図示しない入力デバイス等を介して中断操作を受け付けた場合には、遅延量算出処理を中止することもできる。
上述してきたように、本実施例に係る脈波検出装置10は、被験者の生体が撮影された画像を複数の領域に分割し、各領域で輝度が変化するタイミングの時間差から脈波の遅延量を算出する。このため、本実施例に係る脈波検出装置10では、脈波が伝搬するタイミングを測定する時間分解能を動画のフレーム周期以下に向上できる。したがって、本実施例に係る脈波検出装置10によれば、脈波の遅延量を始め、脈波伝搬速度、脈波到達時刻や脈波周期のゆらぎなどの検出精度を向上させることができる。
上記の実施例1では、各領域P1〜Pnの間で脈波の遅延量が一定であるとの想定の下、脈波の遅延量を算出する算出方法について説明したが、他の算出方法によって脈波の遅延量を算出することもできる。
また、脈波検出装置10は、収縮期血圧の測定値に応じて、各領域の遅延量を算出することもできる。一般に、脈波伝搬速度V1は、収縮期血圧P1に対して下記の1次式(3)で近似できることが知られている。ただし、下記の1次式(3)の傾きa1及び切片a0は、個人差によって値が異なる。このため、事前の測定によってキャリブレーションしておくことが好ましい。
なお、領域の分割数Nは、値を大きくするほど時間分解能を向上させることができるが、値を過度に大きくするとSN比が低下する場合もある。このため、例えば、フレーム間で輝度変化の立ち上がりにずれが発生する領域が少なくとも1つ発生する程度の数に分割数を設定することによって実施例1及び実施例2で説明する遅延量の算出を実現できる。
上記の実施例1及び実施例2では、入力信号としてR信号およびG信号の二種類を用いる場合を例示したが、異なる複数の光波長成分を持つ信号であれば任意の種類の信号および任意の数の信号を入力信号とすることができる。例えば、R、G、B、IRおよびNIRなどの光波長成分が異なる信号のうち任意の組合せの信号を2つ用いることもできるし、また、3つ以上用いることもできる。
上記の実施例1では、脈波検出装置10が上記の脈波検出処理をスタンドアローンで実行する場合を例示したが、クライアントサーバシステムとして実装することもできる。例えば、脈波検出装置10は、脈波検出サービスを提供するWebサーバとして実装することとしてもよいし、アウトソーシングによって脈波検出サービスを提供するクラウドとして実装することとしてもかまわない。このように、脈波検出装置10がサーバ装置として動作する場合には、スマートフォンや携帯電話機等の携帯端末装置やパーソナルコンピュータ等の情報処理装置をクライアント端末として収容することができる。これらクライアント端末からネットワークを介して被験者の顔が映った画像が取得された場合に脈波検出処理を実行し、脈波の検出結果やその検出結果を用いてなされた診断結果をクライアント端末へ応答することによって脈波検出サービスを提供できる。
また、図示した各装置の各構成要素は、必ずしも物理的に図示の如く構成されていることを要しない。すなわち、各装置の分散・統合の具体的形態は図示のものに限られず、その全部または一部を、各種の負荷や使用状況などに応じて、任意の単位で機能的または物理的に分散・統合して構成することができる。例えば、脈波検出装置10が有する各機能部を脈波検出装置10の外部装置としてネットワーク経由で接続するようにしてもよい。
また、各機能部を別の装置がそれぞれ有し、ネットワーク接続されて協働することで、上記の脈波検出装置10の機能を実現するようにしてもよい。
また、上記の実施例で説明した各種の処理は、予め用意されたプログラムをパーソナルコンピュータやワークステーションなどのコンピュータで実行することによって実現することができる。そこで、以下では、図14を用いて、上記の実施例と同様の機能を有する脈波検出プログラムを実行するコンピュータの一例について説明する。
11 カメラ
11a 画像メモリ
12 取得部
13g 第1の分割部
13r 第2の分割部
14g1,14g2,14gn,14r1,14r2,14rn 統計処理部
15g1,15g2,15gn,15r1,15r2,15rn BPF
16−1,16−2,16−n 差演算部
17−1,17−2,17−n 二値化部
18 測定部
19 算出部
Claims (10)
- 撮像装置によって被験者の生体が撮影された画像を取得する取得部と、
前記画像を複数の領域に分割する分割部と、
各領域に含まれる画素が持つ画素値に所定の統計処理を波長成分別に実行することによって各領域の代表値を波長成分別に算出する統計処理部と、
各波長成分間で合成がなされた領域別の代表値が前記画像のフレーム間で変化するタイミングを測定する測定部と、
各領域で代表値が変化するタイミングの時間差から脈波の遅延量を算出する算出部と
を有することを特徴とする脈波検出装置。 - 前記代表値が前記画像のフレーム間で変化するか否かによって代表値を二値化する二値化部をさらに有し、
前記算出部は、
各領域で代表値が変化するタイミングの時間差として、前記画像のフレーム間で二値が反転する領域の数を算出することを特徴とする請求項1に記載の脈波検出装置。 - 前記算出部は、
前記画像のフレーム周期を前記画像のフレーム間で二値が反転する領域の数で分配することによって各領域間の遅延量を算出することを特徴とする請求項1に記載の脈波検出装置。 - 前記領域ごとに同一のフレームで二値の反転が検出された領域の中で脈波の起点となる心臓から直近に位置する変化点の領域となった回数を計数する計数部をさらに有し、
前記算出部は、
前記変化点の領域となった回数の大小によって前記画像のフレーム周期を各領域へ分配する量を変えて各領域間の遅延量を算出することを特徴とする請求項3に記載の脈波検出装置。 - 前記算出部は、
各領域間の遅延量から、脈波伝搬速度、脈波伝搬時刻、脈拍周期、脈拍周期のゆらぎ、もしくは脈波伝搬速度、脈波伝搬時刻、脈拍周期及び脈拍周期のゆらぎの組合せをさらに算出することを特徴とする請求項1〜4のいずれか一つに記載の脈波検出装置。 - 第1の血圧を記憶する記憶部をさらに有し、
前記取得部は、
第1の部位及び前記第1の部位よりも脈波の起点となる心臓から遠い第2の部位を含む画像を取得し、
前記測定部は、
前記第1の部位及び前記第2の部位ごとに、各波長成分間で合成がなされた領域別の代表値が前記画像のフレーム間で変化するタイミングを測定し、
前記算出部は、
前記第1の血圧から導出される第1の脈波伝搬速度を用いて、前記第1の部位が分割された各領域間の遅延量を算出するとともに、前記第2の部位が分割された各領域間の遅延量を算出し、
前記第1の部位の各領域間の遅延量と、前記第1の部位の各領域で代表値が変化するタイミングの時間差とに基づいて当該第1の部位への脈波伝搬時刻を算出するとともに、前記第2の部位の各領域間の遅延量と、前記第2の部位の各領域で代表値が変化するタイミングの時間差とに基づいて当該第2の部位への脈波伝搬時刻を算出し、
前記第1の部位への脈波伝搬時刻及び前記第2の部位への脈波伝搬時刻の差から、前記第1の部位から前記第2の部位までの遅延量を算出し、
前記第1の部位から前記第2の部位までの遅延量と、前記第1の部位及び前記第2の部位の間の距離とから第2の脈波伝搬速度を算出し、
前記第2の脈波伝搬速度から導出される第2の血圧と、前記第1の血圧との差が所定の範囲内になるまで、前記第2の血圧を導出する処理を繰り返し実行することを特徴とする請求項5に記載の脈波検出装置。 - 前記分割部は、
前記被験者の動脈に沿って画像を分割することを特徴とする請求項1〜6のいずれか一つに記載の脈波検出装置。 - 前記分割部は、
前記画像から前記被験者の顔に含まれる第1の顔パーツ及び第2の顔パーツを抽出し、前記第1の顔パーツ及び前記第2の顔パーツによって定まる範囲の画像を分割することを特徴とする請求項7に記載の脈波検出装置。 - コンピュータが、
撮像装置によって被験者の生体が撮影された画像を取得し、
前記画像を複数の領域に分割し、
各領域に含まれる画素が持つ画素値に所定の統計処理を波長成分別に実行することによって各領域の代表値を波長成分別に算出し、
各波長成分間で合成がなされた領域別の代表値が前記画像のフレーム間で変化するタイミングを測定し、
各領域で代表値が変化するタイミングの時間差から脈波の遅延量を算出する
処理を実行することを特徴とする脈波検出方法。 - コンピュータに、
撮像装置によって被験者の生体が撮影された画像を取得し、
前記画像を複数の領域に分割し、
各領域に含まれる画素が持つ画素値に所定の統計処理を波長成分別に実行することによって各領域の代表値を波長成分別に算出し、
各波長成分間で合成がなされた領域別の代表値が前記画像のフレーム間で変化するタイミングを測定し、
各領域で代表値が変化するタイミングの時間差から脈波の遅延量を算出する
処理を実行させることを特徴とする脈波検出プログラム。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2013067667A JP6052005B2 (ja) | 2013-03-27 | 2013-03-27 | 脈波検出装置、脈波検出方法及び脈波検出プログラム |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2013067667A JP6052005B2 (ja) | 2013-03-27 | 2013-03-27 | 脈波検出装置、脈波検出方法及び脈波検出プログラム |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2014188237A JP2014188237A (ja) | 2014-10-06 |
JP6052005B2 true JP6052005B2 (ja) | 2016-12-27 |
Family
ID=51835176
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2013067667A Expired - Fee Related JP6052005B2 (ja) | 2013-03-27 | 2013-03-27 | 脈波検出装置、脈波検出方法及び脈波検出プログラム |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP6052005B2 (ja) |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2015129843A1 (ja) * | 2014-02-27 | 2015-09-03 | 京セラ株式会社 | センサ,センサ装置およびセンサ装置の駆動方法 |
JP6379899B2 (ja) * | 2014-09-12 | 2018-08-29 | 富士通株式会社 | 情報処理装置、脈波計測プログラムおよび脈波計測方法 |
US10736517B2 (en) | 2014-10-09 | 2020-08-11 | Panasonic Intellectual Property Management Co., Ltd. | Non-contact blood-pressure measuring device and non-contact blood-pressure measuring method |
JP6683367B2 (ja) * | 2015-03-30 | 2020-04-22 | 国立大学法人東北大学 | 生体情報計測装置、生体情報計測方法及び生体情報計測プログラム |
WO2016158624A1 (ja) * | 2015-03-30 | 2016-10-06 | 国立大学法人東北大学 | 生体情報計測装置、生体情報計測方法、生体情報表示装置及び生体情報表示方法 |
WO2016163019A1 (ja) * | 2015-04-10 | 2016-10-13 | 株式会社日立製作所 | 生体情報分析システム |
KR102403252B1 (ko) | 2015-07-16 | 2022-05-26 | 삼성전자주식회사 | 그립형 혈압 측정 장치 및 그의 동작 방법 |
EP3422931B1 (en) * | 2016-03-01 | 2021-12-29 | Koninklijke Philips N.V. | Device, system and method for determining a vital sign of a subject |
CN106798547A (zh) * | 2016-12-29 | 2017-06-06 | 乐普(北京)医疗器械股份有限公司 | 一种基于统计的臂踝脉搏波延时的计算方法 |
JP7068339B2 (ja) * | 2017-12-15 | 2022-05-16 | シャープ株式会社 | 血圧測定装置、および、血圧測定方法 |
CN111387966B (zh) * | 2020-03-20 | 2022-08-26 | 中国科学院深圳先进技术研究院 | 信号波重构方法及心率变异性信息的检测设备 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2001358984A (ja) * | 2000-06-15 | 2001-12-26 | Mitsubishi Heavy Ind Ltd | 動画像処理カメラ |
JP2004073398A (ja) * | 2002-08-14 | 2004-03-11 | Rikogaku Shinkokai | 画像処理による運動中の生理指標の推定方法 |
JP3742410B2 (ja) * | 2003-07-24 | 2006-02-01 | 富士通株式会社 | 移動物体の交通流監視システム |
JP5067024B2 (ja) * | 2007-06-06 | 2012-11-07 | ソニー株式会社 | 生体情報取得装置および生体情報取得方法 |
JP5332406B2 (ja) * | 2008-08-28 | 2013-11-06 | 富士通株式会社 | 脈拍計測装置、脈拍計測方法および脈拍計測プログラム |
JP2012169725A (ja) * | 2011-02-10 | 2012-09-06 | Sony Corp | 画像処理装置、画像処理方法、プログラムおよび画像表示装置 |
-
2013
- 2013-03-27 JP JP2013067667A patent/JP6052005B2/ja not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
JP2014188237A (ja) | 2014-10-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP6052005B2 (ja) | 脈波検出装置、脈波検出方法及び脈波検出プログラム | |
EP2979631B1 (en) | Blood flow index calculation method, blood flow index calculation program and blood flow index calculation device | |
WO2016006027A1 (ja) | 脈波検出方法、脈波検出プログラム及び脈波検出装置 | |
JP6125648B2 (ja) | 生体情報取得装置および生体情報取得方法 | |
JP6167614B2 (ja) | 血流指標算出プログラム、血流指標算出装置および血流指標算出方法 | |
JP6098257B2 (ja) | 信号処理装置、信号処理方法及び信号処理プログラム | |
JP6142664B2 (ja) | 脈波検出装置、脈波検出プログラム、脈波検出方法及びコンテンツ評価システム | |
JP6098304B2 (ja) | 脈波検出装置、脈波検出方法及び脈波検出プログラム | |
KR20150059631A (ko) | 심박수 측정 방법 및 이를 위한 장치 | |
JP6248780B2 (ja) | 脈波検出装置、脈波検出方法及び脈波検出プログラム | |
JP2014200390A (ja) | 脈波検出装置、脈波検出方法及び脈波検出プログラム | |
JP6135255B2 (ja) | 心拍測定プログラム、心拍測定方法及び心拍測定装置 | |
JP2014198200A (ja) | 脈波検出装置、脈波検出プログラムおよび脈波検出方法 | |
US20150243027A1 (en) | Image processing device, image processing method, and program | |
US20180242898A1 (en) | Viewing state detection device, viewing state detection system and viewing state detection method | |
JPWO2014002276A1 (ja) | バイタルサイン検出方法、バイタルサイン検出装置及びバイタルサイン検出プログラム | |
JP2015116368A (ja) | 脈拍計測装置、脈拍計測方法及び脈拍計測プログラム | |
US11989884B2 (en) | Method, apparatus and program | |
JP6020015B2 (ja) | 脈波検出装置、脈波検出プログラム及び脈波検出方法 | |
WO2020003910A1 (ja) | 心拍検出装置、心拍検出方法及びプログラム | |
JP6167849B2 (ja) | 脈波検出装置、脈波検出方法及び脈波検出プログラム | |
JP6488722B2 (ja) | 脈波検出装置、脈波検出方法及び脈波検出プログラム | |
JP6167615B2 (ja) | 血流指標算出プログラム、端末装置および血流指標算出方法 | |
JP2014104120A (ja) | 眼球運動解析システム、眼球運動診断装置および眼球運動解析方法 | |
US9978144B2 (en) | Biological information measurement apparatus, biological information measurement method, and computer-readable recording medium |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20151204 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20161021 |
|
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: 20161101 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20161114 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 6052005 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
LAPS | Cancellation because of no payment of annual fees |