一种非接触式桡动脉波形提取方法技术

本发明专利技术公开了一种非接触式桡动脉波形提取的方法，主要解决真实场景下非接触式腕部桡动脉脉搏波提取准确性和健壮性的问题，其实现过程是：（1）通过计算机控制摄像头采集30s正对着人体手腕的视频；（2）对采集的视频进行拉普拉斯金字塔，以实现空间分解，然后对金字塔分解的子带进行IIR滤波，以完成时域滤波，该过程实现视频中桡动脉搏动信号放大；（3）对放大后的视频分块处理，分析所有分块的功率谱，选取信噪比最高的区域作为感兴趣的区域；（4）选取感兴趣的区域作为输入信号，将该区域分为3个信号源，进行独立成分分析，最终得到桡动脉波形。本发明专利技术利用桡动脉的搏动特性，通过视频运动放大提取腕部桡动脉搏动信号，而后进行感兴趣区域精确选取，降低了环境噪声的干扰，同时使用独立成分分析，进一步提高了脉搏波波形的准确性。

A Non-contact Waveform Extraction Method for Radial Artery

The invention discloses a method of non-contact radial artery waveform extraction, which mainly solves the problem of accuracy and robustness of non-contact radial artery pulse wave extraction in real scenes. The realization process is: (1) collecting 30 s video facing human wrist through a computer-controlled camera; (2) Laplacian pyramid is used to realize spatial decomposition of the collected video. Secondly, IIR filtering is applied to the subbands decomposed by pyramid to complete the time domain filtering, which enlarges the radial artery pulsation signal in the video; (3) After the enlarged video is partitioned, the power spectrum of all the subbands is analyzed, and the region with the highest signal-to-noise ratio is selected as the region of interest; (4) The region of interest is selected as the input signal, and the region is divided into three signal sources. Independent component analysis was carried out to obtain the radial artery waveform. The invention utilizes the pulsation characteristics of the radial artery, extracts the wrist radial artery pulsation signal by video motion amplification, then precisely selects the region of interest, reduces the interference of environmental noise, and further improves the accuracy of the pulse waveform by using independent component analysis.

【技术实现步骤摘要】
一种非接触式桡动脉波形提取方法
本专利技术属于计算机视觉领域，特别涉及一种非接触式桡动脉波形提取方法。
技术介绍
近年来，随着计算机视觉技术的深入探索，越来越多的计算机视觉理论知识应用到医学领域，一定程度上推动了医学的发展。计算机视觉借助计算机和摄像头等工具作为媒介，通过相关技术，对医学领域特定问题进行建模和研究，并从计算机视觉的角度探索解决问题方案，为医学的发展提供了一个独辟蹊径的思路。脉搏波是伴随着心脏周期性收缩的动脉血管波动的轨迹，它综合了心脏泵血活动和压力波沿动脉树传递途中携带的各种信息。通过记录脉搏波信号，并获取信息(如周期、幅值和速度等)来评判动脉性能，是心血管疾病防治领域前沿性课题之一。脉搏波信号的准确获取，是获取人体病例信息相关的特征参数，并为临床诊断提供客观的依据的前提。因此，寻求一种简单、可靠的脉搏波检测手段十分关键。目前，非临床手段进行脉搏波检测的方式很多，按照位置分类，可以分为颈动脉、股动脉、指动脉和桡动脉。其中，桡动脉最容易测量也最容易获取稳定的波形。桡动脉是人体前臂主要动脉之一，也是中医切脉的主要依据。桡动脉脉搏波常用于广义传递函数转换为主动脉脉搏波，通过计算力学指标来评估动脉的相关性能。按照检测原理分类，脉搏波检测又可以分为平面张力法和光电容积描记法。平面张力法，即将压力传感器固定于表面动脉，通过传感器记录信号以此记录脉搏波波形，目前国际上常用的传感器有美国的Millar公司的SPT-301，国内的合肥华科电子有限公司的HK-2000B应用较广。光电容积描记法常用于指尖式脉搏波信号的采集与分析，如基于指夹的脉搏波检测仪等。上述各种方法均能够获取桡动脉的脉搏波波形，在不同的应用场景下具有各自的优势。但是接触式的脉搏波提取方法需要在一定专业人员的指导下进行操作，使用流程较为繁琐。其中基于压力传感器的脉搏波获取方式要求与被测人员的皮肤表层亲密接触，对于某些皮肤疾病患者而言，该方式极为不便。对于光电容积描记法而言，需要辅助光管照射人体并将发射器和接收器夹紧人体手指指尖，这在一定程度上会给受测者造成不适。因此，这些方法都无法很好的解决现实场景下桡动脉波形提取的问题。
技术实现思路
本专利技术针对以上问题，提出了一种非接触式桡动脉波形提取方法，该方法无需压力传感器接触人体，也无需辅助光源对人体进行照射，只需要借助普通摄像头捕获桡动脉视频，通过视频放大技术对桡动脉的运动信息进行获取和放大，然后通过进一步降噪处理，从而获取准确的桡动脉波形。本专利技术的技术思路是：首先使用摄像头获取一段30s的桡动脉视频，然后通过欧拉视频放大(EulerianVideoMagnification，EVM)技术对视频进行运动放大，使视频中桡动脉的搏动信号得以增强，并降低运动噪声，再通过分块筛选，选出信噪比最高的区域作为桡动脉感兴趣的区域(Regionofinterest，ROI)，最后将该区域作为待处理的信号，使用独立成分分析(IndependentComponentAnalysis，ICA)做进一步的降噪处理，得到最终的桡动脉波形。本专利技术的采用的技术方案如下：一种非接触式桡动脉波形提取方法，其特征在于，包括如下步骤：(1)视频采集：将一台普通的摄像机设备放置于正对于人的腕部0.3m位置，将摄像机与计算机连接，操作者利用MATLAB软件进行控制摄像头拍摄视频，确保手的腕部桡动脉在所录制的视频帧画面以内；所有的视频以640x480的分辨率、30帧/秒的帧率和RGB颜色空间进行采集，视频长度为30秒，视频编码为MJPEG，视频格式为MP4；(2)视频预处理：读取采集的腕部视频，将视频帧图像转换成矩阵，原始视频像素点在RGB颜色空间内采用8位整型存储，取值范围为[0,255]，为了保持较高的运算精度，进行归一化处理，像素点的取值范围映射到[0,1]区间，使用64位双精度类型存储，每一帧图像Fi到矩阵Mi的映射公式为：(3)视频放大：对预处理后的视频帧序列采用欧拉视频放大方法进行处理，使视频中包含的腕部桡动脉搏动信号可视化；a)空间分解使用拉普拉斯金字塔进行视频帧图像的空间分解，拉普拉斯金字塔的理念是基于高斯金字塔建立一个多分辨率与尺度变化的图像组，高分辨率的图像包含了图像的细节内容，低分辨率图像表象图像的结构信息，因此拉普拉斯金字塔可以良好地表现图像的多种信息；从高斯金字塔的底层开始为图像的原始尺寸图像G0，逐级向上为第i层图像Gi，每一级所在的图像都是在下一级图像的下采样与高斯滤波得到，即第i级图像；那么第i层拉普拉斯金字塔的定义为，第i层的Gi减去第i+1层做过上采样的Gi+1，两者的残差即为第i层的拉普拉斯金字塔图像Li:其中，UP表示上采样，表示卷积运算，g5x5表示5x5的高斯滤波核，基于运算速度和滤波效果的综合考虑，选用5x5大小的滤波核；b)时域滤波利用无限脉冲(InfiniteImpulseResponse，IIR)滤波器对拉普拉斯金字塔图像分解后的图像组，进行时域滤波，以获取桡动脉搏动信号；其中，IIR滤波器设计实现如下，IIR滤波器差分方程为：Y(z)＝(1-r)Y(n-1)z-1+rX(z)(3)IIR滤波器系统函数为：其中0≤r≤1，可以看出(1-r)n逼近x轴方向，是个单边信号；从频域上看，该IIR滤波器的设计实现是一个低通滤波器，它的传输函数，在时域上是一个单边函数，频域上具有低筒滤波的特性；两个所述低通滤波器相减就构成带通滤波器，其对应的带通就是所要获取的桡动脉运动信号；c)放大与重构对时域滤波后的信号进行放大，将R通道的幅值保持不变，G和B通道分别放大100倍，即对其通道内的每个像素点值乘以100；每一帧放大之后都扩展到原始视频帧大小，并与原始视频帧图像相加，得到放大后的视频帧，使运动信号的变化情况能被肉眼所观察到；(4)感兴趣区域分块选取为了进一步选取精准的感兴趣区域(RegionOfInterest，ROI)，对放大后的视频进行分块，分块区域大小为20x20分辨率，640x480的视频帧区域被分成768个区域，每个区域的记为Bi，计算并分析每个分块在时间序列上的频率谱和功率谱，并求信噪比(SignaltonoiseRatio，SNR)；求出每个区域的信噪比之后，选取信噪比最高的区域作为下一步待处理的ROI；运算过程描述如下：ROI＝max{0≤i≤768|SNG(Bi)}(5)其中，SNG表示求功率谱的函数；(5)独立成分分析获取桡动脉运动信号在获取ROI之后，对ROI区域内的信号进一步降噪，使用独立成分分析(IndependentComponentAlalysis,ICA)分离信号和噪声；ICA是基于信号高阶统计特性的分析方法，经ICA分解出的各信号分量之间是相互独立的；设x(t)＝(x1(t),x2(t),...,xm(t))T,t＝1,2,...n,为ROI区域内的m个时间长度为n的观测信号，另设s(t)＝(s1(t),s2(t),...,sm(t))T,t＝1,2,...,n,为m个未知的相互独立的源信号矢量，则瞬时线性混叠盲分离混合方程如下：x＝As(6)其中，A＝(a1,a2,...,an)T为未知时不变的n×m混合矩阵，使用ICA算法从观测信号矢量x(t)中通过矩阵本文档来自技高网...

【技术保护点】
1.一种非接触式桡动脉波形提取方法，其特征在于，包括如下步骤：(1)视频采集：将摄像头固定于人的手腕上方约0.3的位置，使摄像头正对腕部区域，通过计算机控制摄像头获取一段长度约30s的mp4格式的视频；(2)视频预处理：读取采集的腕部视频，将视频帧图像转换成矩阵，原始视频像素点在RGB颜色空间内采用8位整型存储，取值范围为[0,255]，进行归一化处理，像素点的取值范围映射到[0,1]区间，使用64位双精度类型存储，每一帧图像Fi到矩阵Mi的映射公式为：

【技术特征摘要】
1.一种非接触式桡动脉波形提取方法，其特征在于，包括如下步骤：(1)视频采集：将摄像头固定于人的手腕上方约0.3的位置，使摄像头正对腕部区域，通过计算机控制摄像头获取一段长度约30s的mp4格式的视频；(2)视频预处理：读取采集的腕部视频，将视频帧图像转换成矩阵，原始视频像素点在RGB颜色空间内采用8位整型存储，取值范围为[0,255]，进行归一化处理，像素点的取值范围映射到[0,1]区间，使用64位双精度类型存储，每一帧图像Fi到矩阵Mi的映射公式为：(3)视频放大：对预处理后的视频帧序列采用欧拉视频放大方法进行处理，使视频中包含的腕部桡动脉搏动信号可视化；a)空间分解使用拉普拉斯金字塔进行视频帧图像的空间分解，从高斯金字塔的底层开始为图像的原始尺寸图像G0，逐级向上为第i层图像Gi，每一级所在的图像都是在下一级图像的下采样与高斯滤波得到，即第i级图像；那么第i层拉普拉斯金字塔的定义为，第i层的Gi减去第i+1层做过上采样的Gi+1，两者的残差即为第i层的拉普拉斯金字塔图像Li:其中，UP表示上采样，表示卷积运算，g5x5表示5x5的高斯滤波核，基于运算速度和滤波效果的综合考虑，选用5x5大小的滤波核；b)时域滤波利用IIR滤波器对拉普拉斯金字塔图像分解后的图像组，进行时域滤波，以获取桡动脉搏动信号；其中，IIR滤波器设计实现如下，IIR滤波器差分方程为：Y(z)＝(1-r)Y(n-1)z-1+rX(z)(3)IIR滤波器系统函数为：其中0≤r≤1，可以看出(1-r)n逼近x轴方向，是个单边信号；两个所述IIR滤波器相减就构成带通滤波器，其对应的带通就是所要获取的桡动脉运动信号；c)放大与重构对时域滤波后的信号进行放大，将R通道的幅值保持不变，G和B通道分别放大1...

【专利技术属性】
技术研发人员：杨学志，李江山，吴秀，方帅，金兢，
申请(专利权)人：合肥工业大学，
类型：发明
国别省市：安徽,34