CN109512418A - 一种容积血流脉搏成像降低运动噪声方法 - Google Patents

一种容积血流脉搏成像降低运动噪声方法 Download PDF

Info

Publication number
CN109512418A
CN109512418A CN201811228837.6A CN201811228837A CN109512418A CN 109512418 A CN109512418 A CN 109512418A CN 201811228837 A CN201811228837 A CN 201811228837A CN 109512418 A CN109512418 A CN 109512418A
Authority
CN
China
Prior art keywords
signal
motion artifacts
blood flow
flow pulse
camera
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.)
Pending
Application number
CN201811228837.6A
Other languages
English (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.)
Harbin University of Science and Technology
Original Assignee
Harbin University of Science and Technology
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 Harbin University of Science and Technology filed Critical Harbin University of Science and Technology
Priority to CN201811228837.6A priority Critical patent/CN109512418A/zh
Publication of CN109512418A publication Critical patent/CN109512418A/zh
Pending legal-status Critical Current

Links

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/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/0059—Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence
    • 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/026—Measuring blood flow
    • A61B5/0295—Measuring blood flow using plethysmography, i.e. measuring the variations in the volume of a body part as modified by the circulation of blood therethrough, e.g. impedance plethysmography

Abstract

本发明涉及容积脉搏波临床的检测技术领域，具体为一种容积血流脉搏成像降低运动噪声方法，该方法包括如下步骤：步骤a：提取人脸区域的确定；步骤b：信号源和噪声源的提取；步骤c：检测者运动信号的确定；步骤d：噪音信号的处理；本发明通过摄像头提取容积血流脉搏信号的过程中，很好的降低运动噪声的功能，对带有运动噪声和容积血流脉搏波的输入信号，通过运动噪声参考信号进行滤波降噪，很好的降低了运动噪声的影响，从而为提取心率和呼吸率信息提高了准确度。

Description

一种容积血流脉搏成像降低运动噪声方法
技术领域
本发明涉及容积脉搏波临床的检测技术领域，尤其涉及一种容积血流脉搏成像降低运动噪声方法。
背景技术
容积脉搏波临床的检测方法是PPG，它是一种光学测量方法检测组织微血管中血容量变化。PPG信号的检测装置主要由一个照亮组织的光源和一个测量与组织内血容量相关光强度变化信息的光电检测器组成；指夹夹住手指端并屏蔽外界环境光，指夹内的光源发射一定波长的光照射到手指端表面，由于皮肤、血液等对光的吸收或反射能力不同，且心脏的舒张和收缩运动使动脉血容量呈周期性脉动变化，从而使得光电传感器接收到光的强度随血容积的改变也呈周期性的变化，可以得到容积脉搏波。经过半个多世纪的研究，PPG有着广泛的临床应用；PPG可以提取许多的生理信息，包括心率、呼吸率、SpO2和血压等，另外还可以评估血管功能参数，如动脉硬化强度、血管壁弹性、左心室射血时间和血压变异性等；PPG被当前临床医学界公认的检测容积脉搏波最可靠、最准确、最好的方法，并被当作容积脉搏波金标准的检测方法。
在环境光条件下使用摄像头可以提取容积脉搏波；它是借助容积脉搏成像(Photoplethysmographic imaging,PPGi)在活体组织中检测血液容积变化的方法。在活体组织中，皮肤、肌肉和组织等对光的吸收在整个血液循环中是保持恒定不变的，而皮肤内的血液容积在心脏收缩舒张作用下呈脉动性变化；当心脏收缩时外周血管血容量最多，光吸收量也最大，皮肤表面的光强也最小；而心脏舒张时外周血管血容量最少，光吸收量也最小，皮肤表面的光强也最大；摄像头对人体表面皮肤成像，成像亮度随着光强呈脉动性变化，通过对图像像素平均使二维信号转化为一维信号就可以获得容积脉搏波。
人的体表区域如人脸放到一个普通摄像头或摄像机镜头前1米左右的距离，利用图像分割技术提取体表区域的图像，并平均图像阵列像素使二维信号转化为一维信号提取容积脉搏波，再通过现代信号处理技术可以获得心率、心率变异性、脉率、脉率变异性、呼吸率等生命体征信息。
通过普通摄像头距离人1米左右采集人脸区域图像并把此区域图像的像素值平均，使得人脸区域的视频变为一个信号，此信号可以提取容积脉搏波，进而计算出人体心率，呼吸率等生理信息；然而，由于受到运动伪迹的限制，摄像头在检测容积脉搏波上受到极大干扰。
发明内容
本发明克服了上述现有技术的不足，提供了一种容积血流脉搏成像降低运动噪声方法。本发明通过摄像头提取容积血流脉搏信号的过程中，很好的降低运动噪声的功能，对带有运动噪声和容积血流脉搏波的输入信号，通过运动噪声参考信号进行滤波降噪，很好的降低了运动噪声的影响；从而为提取心率和呼吸率信息提高了准确度。
本发明的技术方案：
一种容积血流脉搏成像降低运动噪声方法，包括如下步骤：
步骤a：提取区域的确定：使提取摄像头与检测者的间距保持在0.5～1m；
步骤b：信号源和噪声源的提取：使检测者的脸部处于摄像头检测区域内部，确定人脸区域，并把此区域图像的像素值平均，使得人脸区域的视频变为一个信号，此信号为摄像头提取检测者容积脉搏波信号和检测者头部运动产生的噪声信号，记为d(n)，括号中的n表示第n个采样；
步骤c：检测者运动信号的确定：将检测者头部的运动按照x、y、z三个方向描述，当人脸检测区域面积大小变化时，认为人脸在z轴变化，当人脸检测区域在横向和纵向移动时，认为人脸在x和y轴变化；摄像头采集检测者头部在人脸检测区域变化在x、y、z三个方向上的物理变化，测量到运动产生的加速度，从而表征运动噪声，记为表示为：
其中，H为自适应滤波器系数数量；
步骤d：噪音信号的处理：通过检测计算机内的自适应滤波器将噪声信号进行消噪处理；对带有运动噪声和容积血流脉搏波的输入信号，通过运动噪声参考信号进行滤波降噪，很好的降低了运动噪声的影响；从而为提取心率和呼吸率信息提高了准确度。
上述的一种容积血流脉搏成像降低运动噪声方法，所述步骤d中自适应滤波器滤波的滤波方法具体为：自适应滤波的目的是从摄像头获取的信号中滤除运动噪声信号提取与容积脉搏血流相关信号，输出x'(n)，具体过程如下：
其中，WT(n)为自适应滤波系数向量，表示为：WT(n)＝[w1(n),w2(n),...,wH(n)]；
摄像头获得的人脸区域信号与运动噪声滤波后信号的差值如下所示：
自适应数字滤波器的滤波系数WT(n)的变化与信号e(n)有关，根据e(n)的值大小而自动调整WT(n)，使之适合下一时刻的输入u(n+1)，以便使输出x'(n+1)接近于所期望的容积血流脉搏相关信号。
上述的一种容积血流脉搏成像降低运动噪声方法，所述调整WT(n)的方法为最小均方根方法、最小二乘法、变换域方法或共扼梯度方法中的一种。
本发明相对于现有技术具有以下有益效果：
本发明基于摄像头提取容积血流脉搏信号的过程中，很好的降低运动噪声的功能，对带有运动噪声和容积血流脉搏波的输入信号，通过运动噪声参考信号进行滤波降噪，很好的降低了运动噪声的影响；从而为提取心率和呼吸率信息提高了准确度。
附图说明
图1为本发明的目的是从输出信号中提取真实准确的期望信号的工作方法示意图。
实施例
以下将结合附图对本发明进行详细说明。
实施例一：
本实施例的容积血流脉搏成像降低运动噪声方法，包括如下步骤：
步骤a：提取区域的确定：使提取摄像头与检测者的间距保持在0.5～1m；
步骤b：信号源和噪声源的提取：使检测者的脸部处于摄像头检测区域内部，确定人脸区域，并把此区域图像的像素值平均，使得人脸区域的视频变为一个信号，此信号为摄像头提取检测者容积脉搏波信号和检测者头部运动产生的噪声信号，记为d(n)，括号中的n表示第n个采样；
步骤c：检测者运动信号的确定：将检测者头部的运动按照x、y、z三个方向描述，当人脸检测区域面积大小变化时，认为人脸在z轴变化，当人脸检测区域在横向和纵向移动时，认为人脸在x和y轴变化；摄像头采集检测者头部在人脸检测区域变化在x、y、z三个方向上的物理变化，测量到运动产生的加速度，从而表征运动噪声，记为表示为：
其中，H为自适应滤波器系数数量；
步骤d：噪音信号的处理：通过检测计算机内的自适应滤波器将噪声信号进行消噪处理；对带有运动噪声和容积血流脉搏波的输入信号，通过运动噪声参考信号进行滤波降噪，很好的降低了运动噪声的影响；从而为提取心率和呼吸率信息提高了准确度。
实施例二：
本实施例的容积血流脉搏成像降低运动噪声方法，在实施例一的基础上，进一步限定所述步骤d中自适应滤波器滤波的滤波方法具体为：自适应滤波的目的是从摄像头获取的信号中滤除运动噪声信号提取与容积脉搏血流相关信号，输出x'(n)，具体过程如下：
其中，WT(n)为自适应滤波系数向量，表示为：WT(n)＝[w1(n),w2(n),...,wH(n)]；
摄像头获得的人脸区域信号与运动噪声滤波后信号的差值如下所示：
自适应数字滤波器的滤波系数WT(n)的变化与信号e(n)有关，根据e(n)的值大小而自动调整WT(n)，使之适合下一时刻的输入u(n+1)，以便使输出x'(n+1)接近于所期望的容积血流脉搏相关信号。
实施例三：
本实施例的容积血流脉搏成像降低运动噪声方法，在实施例二的基础上，进一步限定所述调整WT(n)的方法为最小均方根方法、最小二乘法、变换域方法或共扼梯度方法中的一种。
实施例四：
本实施例的容积血流脉搏成像降低运动噪声方法，在实施例三的基础上，进一步限定本实施例为通过最小均方根法调整滤波器系数WT(n)，最小均方根自适应滤波算法原理是通过调整滤波器系数WT(n)，使滤波信号与期望信号之间误差的均方值为最小；滤波器系数WT(n)的计算通过误差值在梯度方向上调整滤波器系数，并快速收敛使得误差的均方值为最小；其算法特点是复杂度低，不需要大量的计算，并且收敛速度块。最小均方根法定义目标函数为WT(n)E[e(n)2]达到最小时，即实现了最优滤波。所以公式可以变换如下：
为了获得最小误差的均方值J(n)，在滤波器系数WT(n)在梯度方向上进行优化，通过更新WT(n)，达到快速收敛，WT(n)优化公式如下所示：
以上实施例只是对本专利的示例性说明，并不限定它的保护范围，本领域技术人员还可以对其局部进行改变，只要没有超出本专利的精神实质，都在本专利的保护范围内。

Claims (3)

1.一种容积血流脉搏成像降低运动噪声方法，其特征在于，包括如下步骤：
步骤a：提取区域的确定：使提取摄像头与检测者的间距保持在0.5～1m；
步骤b：信号源和噪声源的提取：使检测者的脸部处于摄像头检测区域内部，确定人脸区域，并把此区域图像的像素值平均，使得人脸区域的视频变为一个信号，此信号为摄像头提取检测者容积脉搏波信号和检测者头部运动产生的噪声信号，记为d(n)，括号中的n表示第n个采样；
步骤c：检测者运动信号的确定：将检测者头部的运动按照x、y、z三个方向描述，当人脸检测区域面积大小变化时，认为人脸在z轴变化，当人脸检测区域在横向和纵向移动时，认为人脸在x和y轴变化；摄像头采集检测者头部在人脸检测区域变化在x、y、z三个方向上的物理变化，测量到运动产生的加速度，从而表征运动噪声，记为表示为：
其中，H为自适应滤波器系数数量；
步骤d：噪音信号的处理：通过检测计算机内的自适应滤波器将噪声信号进行消噪处理；对带有运动噪声和容积血流脉搏波的输入信号，通过运动噪声参考信号进行滤波降噪，很好的降低了运动噪声的影响；从而为提取心率和呼吸率信息提高了准确度。
2.如权利要求1所述的一种容积血流脉搏成像降低运动噪声方法，其特征在于，所述步骤d中自适应滤波器滤波的滤波方法具体为：自适应滤波的目的是从摄像头获取的信号中滤除运动噪声信号提取与容积脉搏血流相关信号，输出x'(n)，具体过程如下：
其中，WT(n)为自适应滤波系数向量，表示为：WT(n)＝[w1(n),w2(n),...,wH(n)]；
摄像头获得的人脸区域信号与运动噪声滤波后信号的差值如下所示：
自适应数字滤波器的滤波系数WT(n)的变化与信号e(n)有关，根据e(n)的值大小而自动调整WT(n)，使之适合下一时刻的输入u(n+1)，以便使输出x'(n+1)接近于所期望的容积血流脉搏相关信号。
3.如权利要求2所述的一种容积血流脉搏成像降低运动噪声方法，其特征在于，所述调整WT(n)的方法为最小均方根方法、最小二乘法、变换域方法或共扼梯度方法中的一种。
CN201811228837.6A 2018-10-22 2018-10-22 一种容积血流脉搏成像降低运动噪声方法 Pending CN109512418A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811228837.6A CN109512418A (zh) 2018-10-22 2018-10-22 一种容积血流脉搏成像降低运动噪声方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811228837.6A CN109512418A (zh) 2018-10-22 2018-10-22 一种容积血流脉搏成像降低运动噪声方法

Publications (1)

Publication Number Publication Date
CN109512418A true CN109512418A (zh) 2019-03-26

Family

ID=65772737

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811228837.6A Pending CN109512418A (zh) 2018-10-22 2018-10-22 一种容积血流脉搏成像降低运动噪声方法

Country Status (1)

Country Link
CN (1) CN109512418A (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110353648A (zh) * 2019-08-28 2019-10-22 华东师范大学 一种基于双摄像头的实时心率检测方法
CN111657936A (zh) * 2020-06-16 2020-09-15 河南大学 基于小波变换和全变差正则化的信号去噪方法
CN113397497A (zh) * 2021-06-10 2021-09-17 维沃移动通信有限公司 信号处理方法、装置和电子设备
CN114305370A (zh) * 2021-12-17 2022-04-12 南京理工大学 针对rPPG信号中运动干扰的自适应滤波方法及系统

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2012099538A1 (en) * 2011-01-20 2012-07-26 Nitto Denko Corporation A device and method for removal of ambient noise signal from a photoplethysmograph
CN103271743A (zh) * 2012-12-10 2013-09-04 中国人民解放军第一五二中心医院 一种基于成像设备的非接触式血氧饱和度测量装置
CN103271734A (zh) * 2012-12-10 2013-09-04 中国人民解放军第一五二中心医院 一种基于低端成像设备的心率测量方法
US20170209055A1 (en) * 2016-01-22 2017-07-27 Fitbit, Inc. Photoplethysmography-based pulse wave analysis using a wearable device
CN107334469A (zh) * 2017-07-24 2017-11-10 北京理工大学 基于支持向量机的非接触多人心率测量方法及装置
CN108478206A (zh) * 2018-02-02 2018-09-04 北京邮电大学 运动状态下基于脉搏波的心率监测方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2012099538A1 (en) * 2011-01-20 2012-07-26 Nitto Denko Corporation A device and method for removal of ambient noise signal from a photoplethysmograph
CN103271743A (zh) * 2012-12-10 2013-09-04 中国人民解放军第一五二中心医院 一种基于成像设备的非接触式血氧饱和度测量装置
CN103271734A (zh) * 2012-12-10 2013-09-04 中国人民解放军第一五二中心医院 一种基于低端成像设备的心率测量方法
US20170209055A1 (en) * 2016-01-22 2017-07-27 Fitbit, Inc. Photoplethysmography-based pulse wave analysis using a wearable device
CN107334469A (zh) * 2017-07-24 2017-11-10 北京理工大学 基于支持向量机的非接触多人心率测量方法及装置
CN108478206A (zh) * 2018-02-02 2018-09-04 北京邮电大学 运动状态下基于脉搏波的心率监测方法

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
SANG HYUN KIM ET AL.: "Adaptive Noise Cancellation Using Accelerometers for the PPG Signal from Forehead", 《2007 29TH ANNUAL INTERNATIONAL CONFERENCE OF THE IEEE ENGINEERING IN MEDICINE AND BIOLOGY SOCIETY》 *
TIM SCHACK ET AL.: "A new method for heart rate monitoring during physical exercise using photoplethysmographic signals", 《2015 23RD EUROPEAN SIGNAL PROCESSING CONFERENCE (EUSIPCO)》 *
曾周末，万柏坤: "《动态数据建模与处理》", 31 December 2005, 天津大学出版社 *
杨磊: "《网络视频监控技术》", 30 September 2017, 中国传媒大学出版社 *
雷玉堂: "《光电检测技术》", 30 June 2009, 中国计量出版社 *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110353648A (zh) * 2019-08-28 2019-10-22 华东师范大学 一种基于双摄像头的实时心率检测方法
CN111657936A (zh) * 2020-06-16 2020-09-15 河南大学 基于小波变换和全变差正则化的信号去噪方法
CN111657936B (zh) * 2020-06-16 2022-04-12 河南大学 基于小波变换和全变差正则化的信号去噪方法
CN113397497A (zh) * 2021-06-10 2021-09-17 维沃移动通信有限公司 信号处理方法、装置和电子设备
CN114305370A (zh) * 2021-12-17 2022-04-12 南京理工大学 针对rPPG信号中运动干扰的自适应滤波方法及系统

Similar Documents

Publication Publication Date Title
CN109512418A (zh) 一种容积血流脉搏成像降低运动噪声方法
CN105877723B (zh) 无创连续血压测量装置
CN110662480B (zh) 用于测量和处理对象的生理信号的设备、系统和方法
US20140275832A1 (en) Device and method for obtaining vital sign information of a subject
CN110647815A (zh) 一种基于人脸视频图像的非接触式心率测量方法及系统
US20170119304A1 (en) Device, system and method for detecting apnoea of a subject
WO2018210985A1 (en) Reliable acquisition of photoplethysmographic data
US20120157791A1 (en) Adaptive time domain filtering for improved blood pressure estimation
CN112233813A (zh) 一种基于ppg的非接触式无创心率呼吸测量方法及系统
CN110236508A (zh) 一种无创血压连续监测方法
JP6620999B2 (ja) 生体情報計測装置、生体情報計測プログラム、及び生体情報計測方法
CN113499059B (zh) 基于光纤传感非接触式的bcg信号处理系统及方法
EP3391340B1 (en) A method, information processing apparatus and server for determining a physiological parameter of an individual
Junior et al. Estimation of blood pressure and pulse transit time using your smartphone
CN112545472B (zh) 一种ppg信号质量评估的方法、装置、设备及储存介质
CN106618504A (zh) 一种具备自适应校准的高准确度连续无创血压测量装置
WO2023185873A1 (zh) 一种基于多阶多模态的非袖带式血压测量装置及方法
CN109498022A (zh) 一种基于光电容积脉搏波的呼吸频率提取方法
Volynsky et al. Blood peripheral circulation assessment method based on combined use of the video-capillaroscopy, imaging photoplethysmography, and electrocardiography
Wuerich et al. A feature-based approach on contact-less blood pressure estimation from video data
TW201524465A (zh) 動脈波分析方法及其系統
Chatterjee et al. Non-invasive cardiovascular monitoring
Chatterjee et al. Algorithm to Calculate Heart Rate and Comparison of Butterworth IIR and Savitzky-Golay FIR Filter
CN111588385B (zh) 一种血氧饱和度测量方法、装置、存储介质及电子设备
CN114331998A (zh) 非接触式的心肺耦合评价方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20190326