CN107260173B - 一种基于摄像头与球状标记物的呼吸测量方法 - Google Patents
一种基于摄像头与球状标记物的呼吸测量方法 Download PDFInfo
- Publication number
- CN107260173B CN107260173B CN201710519852.5A CN201710519852A CN107260173B CN 107260173 B CN107260173 B CN 107260173B CN 201710519852 A CN201710519852 A CN 201710519852A CN 107260173 B CN107260173 B CN 107260173B
- Authority
- CN
- China
- Prior art keywords
- image
- circular areas
- border circular
- camera
- marker
- 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
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/08—Detecting, measuring or recording devices for evaluating the respiratory organs
- A61B5/0816—Measuring devices for examining respiratory frequency
-
- 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/0062—Arrangements for scanning
- A61B5/0064—Body surface scanning
-
- 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/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/7253—Details of waveform analysis characterised by using transforms
- A61B5/7257—Details of waveform analysis characterised by using transforms using Fourier transforms
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/40—Extraction of image or video features
- G06V10/44—Local feature extraction by analysis of parts of the pattern, e.g. by detecting edges, contours, loops, corners, strokes or intersections; Connectivity analysis, e.g. of connected components
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/40—Extraction of image or video features
- G06V10/46—Descriptors for shape, contour or point-related descriptors, e.g. scale invariant feature transform [SIFT] or bags of words [BoW]; Salient regional features
- G06V10/478—Contour-based spectral representations or scale-space representations, e.g. by Fourier analysis, wavelet analysis or curvature scale-space [CSS]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/08—Feature extraction
- G06F2218/10—Feature extraction by analysing the shape of a waveform, e.g. extracting parameters relating to peaks
Abstract
本发明提供一种基于摄像头与球状标记物的呼吸测量方法,球状标记物附着在人体胸腹部,利用摄像头采集标记物的图像,标记物与摄像机之间的距离变化会引起标记物在摄像机内的成像大小发生变化,通过提取标记物的成像区域轮廓半径并记录连续图像内区域轮廓半径的大小变化即可反应人体的呼吸运动,跟踪人体的呼吸运动模式和计算呼吸频率。本发明通过简单的球状标记物和摄像头实现了呼吸模式的跟踪与呼吸频率的计算,避免昂贵的传感器连接及其给人体带来的不适感,此外也避免其他图像测量呼吸算法中复杂的跟踪与计算算法,有助于日常生活中进行呼吸频率的测量与呼吸模式的跟踪。
Description
技术领域
本发明涉及生理信号监测与分析技术领域,涉及到人体呼吸信号的模式检测与呼吸率测量,特别涉及到利用球状标记物与摄像头进行呼吸模式追踪与呼吸率测量的方法。
背景技术
针对呼吸信号的监测与测量,主要利用相应的阻抗传感器,温度传感器与人体直接连接,测量人体呼吸的气流温度、流速及身体的运动,往往需要复杂的连接与使用,并会给人体带来一定不适感。针对这些问题,产生了一种基于彩色摄像头的呼吸测量方法,其中Stefan Wiesner通过测量附着在人体腹部的多个块状标记物在彩色图像中色彩饱和度变化的方法测量呼吸,然而此方法只适用于卧姿,并且需要多个标记物辅助测量。Zhang和Shao分别利用不同的方法计算彩色图像中人体的细微运动变化从而测量呼吸,然而这两种方法需要复杂的计算与区域跟踪算法和昂贵的CCD相机,此外对于人体的细微运动十分敏感,容易产生较大的误差。
发明内容
本发明的主要目的是为了解决现有技术问题,提供一种利用摄像头和球状标记物进行呼吸测量的方法。
本发明的技术方案为:
一种基于摄像头与球状标记物的呼吸测量方法,包括以下步骤:
第一步,采集球状标记物图像
1.1)由于呼吸运动时人体的胸腹部运动变化最为明显,将球状标记物用双面胶附着在人体胸腹部。
1.2)调整摄像头位置,保证球状标记物位于成像区域内。
1.3)打开摄像头,开始采集图像,采集时间为1分钟;
第二步,求解该图像中球状标记物的圆形区域的半径
2.1)利用球状标记物的颜色特征,通过图像中R,G,B颜色通道的差值进行标记物的区域提取,提取得到图像中球状标记物的圆形区域;圆形区域的提取公式如下:
R-G>threshold&R-B>threshold
其中,R,G,B为彩色图像的三种颜色通道;threshold为根据相机分辨率,实验光照条件所定义的阈值。
随后提取该圆形区域的轮廓获得圆形区域轮廓点集的坐标(Xi,Yi),i∈(1,2,3...N),其中,X为点在图像内的横坐标,Y为点在图像内的纵坐标,N为圆形区域轮廓点集内点的数量。
2.2)利用最小二乘法对提取到的圆形区域轮廓进行圆形拟合,计算该圆形区域的半径。采用最小二乘法拟合的圆形曲线为:
R2=(x-A)2+(y-B)2
其中,(x,y)代表圆形上的点坐标;R代表圆形半径;(A,B)代表圆心坐标。
令a=-2A,b=-2B,c=A2+B2-R2则圆形曲线方程的另一个形式为:
x2+y2+ax+by+c=0
将A,B,R的求解问题转换为求解:参数a,b,c:
通过步骤2.1)获取到的圆形区域轮廓点集(Xi,Yi),则样本集
(Xi,Yi),i∈(1,2,3...N)中点到圆心的距离为di:
其与圆半径平方R2的差为:
δi=di 2-R2=(Xi-A)2+(Yi-B)2-R2=Xi 2+Yi 2+aXi+bYi+c
令最小作为判定条件,采用最小二乘法即可求解出a,b,c,从而求解得到圆形区域的半径R。
第三步,通过连续多帧图像处理获得呼吸波形,利用傅里叶变换计算呼吸频率
3.1)对连续多帧图像进行第二步所述的处理,获得随时间连续变化的圆形区域半径波形;人体呼吸运动能够改变球状标记物与相机之间的距离,根据相机成像原理,由于物体与相机间距离的改变,物体在相机中的成像大小也会随之改变,距离越近,成像区域越大,距离越远,成像区域越小,因此人体呼吸运动变化反映彩色图像中球状标记物的区域大小,两者具有相同的周期性,因而随时间变化的圆形区域半径波形即可代表人体的呼吸波形。
3.2)对步骤3.1)中获得的呼吸波形进行傅里叶变换,得到相应的呼吸频谱,其中,频谱图中的最高峰即代表受试者在该段时间内的呼吸频率。
本发明的有益效果为:本发明通过提取球状标记物在彩色图像中的圆形区域轮廓,并利用拟合算法求出该轮廓的半径,通过记录连续图像中轮廓半径的变化即可反应人体的呼吸模式并进行呼吸频率的计算。本发明可以实现远程,低成本的人体呼吸测量,通过简单的区域轮廓提取避免了复杂的传感器连接与人体运动跟踪算法即可显示人体呼吸波形,从而有助于日常生活中的方便使用。
附图说明
图1是本发明的呼吸测量与检测的系统框图;
图2是图像中提取到的球状标记物区域轮廓;
图3(a)是提取到的呼吸模式波形;
图3(b)是提取到的受试者正常呼吸波形;
图4是呼吸波形的频谱图。
具体实施方式
为使本发明实施的目的、技术方案和优点更加清晰,以下结合本发明的技术方案和附图作进一步详细描述:
利用球状标记物及摄像头测量呼吸的方法,其系统总体框图如图1所示。该方法可分为三个环节,分别为采集球状标记物图像、提取图像中球状标记物圆形区域轮廓半径与建立呼吸测量模型。具体步骤如下:
步骤A.采集球状标记物的图像主要包含如下步骤:
A1.由于人体呼吸运动时胸腹部的运动最为明显,测量时将球状标记物用双面胶附着在人体的胸部或腹部,受试者可以呈现坐姿和卧姿。
A2.利用电脑端程序打开摄像头,根据电脑中的图像调整摄像头的位置,保证其位于球状标记物前方,并能采集到完整的球状标记物区域。
A3.利用电脑端编写的程序打开摄像头并采集相应的图像存储到电脑中,图像采集的频率为25Hz,采集到的图像分辨率为640×480,采集时间为1分钟。
步骤B.提取球状标记物的圆形区域轮廓,利用最小二乘法拟合进行圆形区域半径计算。所述步骤B具体包括如下步骤:
B1.利用球状标记物的颜色特征,通过图像R,G,B颜色通道的差值进行区域提取,基本公式如下:
R-G>threshold&R-B>threshold
其中,R,G,B分别为彩色图像的三种颜色通道,threshold为根据相机分辨率,实验光照条件所定义的阈值,本次实验设置为50,随后对提取到的区域进行相应的膨胀腐蚀操作,从而保证区域选择的正确性,最终提取到的区域轮廓如图2所示,随后提取圆形区域轮廓点集的坐标(Xi,Yi),i∈(1,2,3...N),其中,X为点在图像内的横坐标,Y为点在图像内的纵坐标,N为区域圆形区域轮廓点集的数量。
B2.利用最小二乘法进行圆形拟合,并求解相应的圆形区域半径。最小二乘法拟合的圆形曲线为:
R2=(x-A)2+(y-B)2
其中(x,y)代表圆形上的点坐标,R代表圆形半径,(A,B)代表圆心坐标,令a=-2A,b=-2B,c=A2+B2-R2则可以得到圆形曲线方程的另一个形式:
x2+y2+ax+by+c=0
只要求出参数a,b,c就可以求出相应的A,B,R:
通过步骤B1获取到的圆形区域轮廓点集(Xi,Yi),则样本集
(Xi,Yi),i∈(1,2,3...N)中点到圆心的距离为di:
其与圆半径平方R2的差为:
δi=di 2-R2=(Xi-A)2+(Yi-B)2-R2=Xi 2+Yi 2+aXi+bYi+c
令最小作为判定条件,采用最小二乘法即可求解出a,b,c,从而求解得到圆形区域的半径R。
步骤C.通过多帧图像处理获得呼吸波形,利用傅里叶变换计算呼吸频率。
步骤C具体包括:
C1.对多帧图像进行步骤B的处理得到连续变化的轮廓区域半径如图3(a),图3(b)所示,其中横轴代表时间,纵轴代表区域半径大小,根据摄像头成像原理可知该区域半径的变化波形即可反应人体的呼吸波形,图3(a)为人体进行不同呼吸状态时的呼吸波形,可以看出很好的反应出人体的呼吸模式,包括正常呼吸,快速呼吸,深度呼吸,屏住呼吸,图3(b)显示人体正常呼吸的波形。
C2.对图3(b)提取到的呼吸波形作傅里叶变换,得到呼吸波形的频谱如图4所示,频谱的最大峰值为0.21Hz,代表人体的呼吸频率为0.35×60=21次每分钟。
以上所述,仅为本发明较佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,根据本发明的技术方案及其发明构思加以等同替换或改变,都应涵盖在本发明的保护范围之内。
Claims (2)
1.一种基于摄像头与球状标记物的呼吸测量方法,其特征在于以下步骤:
第一步,采集球状标记物图像
1.1)将球状标记物贴附在人体胸腹部;
1.2)调整摄像头位置,保证球状标记物位于成像区域内;
1.3)打开摄像头,开始采集图像,采集时间为1分钟;
第二步,求解该图像中球状标记物的圆形区域的半径
2.1)利用球状标记物的颜色特征,通过图像中R,G,B颜色通道的差值进行球状标记物的区域提取,得到图像中球状标记物的圆形区域;圆形区域的提取公式如下:
R-G>threshold&R-B>threshold
其中,R,G,B为彩色图像的三种颜色通道;threshold为根据相机分辨率,光照条件所定义的阈值;
提取该圆形区域的轮廓,获得圆形区域轮廓点集坐标(Xi,Yi),i∈(1,2,3...N),其中,X为点在图像内的横坐标,Y为点在图像内的纵坐标,N为圆形区域轮廓点集内点的数量;
2.2)利用最小二乘法对提取到的圆形区域轮廓进行圆形拟合,计算得到该圆形区域的半径;
第三步,通过连续多帧图像处理获得呼吸波形,利用傅里叶变换计算呼吸频率
3.1)对连续多帧图像进行第二步所述的处理,获得随时间连续变化的圆形区域半径波形,且该随时间变化的圆形区域半径波形代表人体呼吸波形;
3.2)对步骤3.1)中获得的呼吸波形进行傅里叶变换,得到相应的呼吸频谱,其中,频谱图中的最高峰即代表受试者在该段时间内的呼吸频率。
2.根据权利要求1所述的一种基于摄像头与球状标记物的呼吸测量方法,其特征在于,所述的步骤2.2)得到圆形区域的半径的步骤为:
采用最小二乘法拟合的圆形曲线为:
R2=(x-A)2+(y-B)2
其中,(x,y)代表圆形上的点坐标;R代表圆形半径;(A,B)代表圆心坐标;
令a=-2A,b=-2B,c=A2+B2-R2则圆形曲线方程的另一个形式为:
x2+y2+ax+by+c=0
将A,B,R的求解问题转换为求解:参数a,b,c:
A=-a/2,B=-b/2,
通过步骤2.1)获取到的圆形区域轮廓点集(Xi,Yi),i∈(1,2,3...N)中点到圆心的距离为di:
其与圆半径平方R2的差为:
δi=di 2-R2=(Xi-A)2+(Yi-B)2-R2=Xi 2+Yi 2+aXi+bYi+c
令最小作为判定条件,采用最小二乘法求解得到a,b,c,进而得到圆形区域的半径R。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710519852.5A CN107260173B (zh) | 2017-06-30 | 2017-06-30 | 一种基于摄像头与球状标记物的呼吸测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710519852.5A CN107260173B (zh) | 2017-06-30 | 2017-06-30 | 一种基于摄像头与球状标记物的呼吸测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107260173A CN107260173A (zh) | 2017-10-20 |
CN107260173B true CN107260173B (zh) | 2019-05-10 |
Family
ID=60070406
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710519852.5A Active CN107260173B (zh) | 2017-06-30 | 2017-06-30 | 一种基于摄像头与球状标记物的呼吸测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107260173B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110236547B (zh) * | 2018-03-09 | 2022-10-28 | 浙江清华柔性电子技术研究院 | 呼吸频率的检测方法 |
CN110051356A (zh) * | 2019-04-04 | 2019-07-26 | 北京贝麦克斯科技有限公司 | 人体呼吸状态信息的获取方法和装置 |
CN110269624B (zh) * | 2019-07-16 | 2024-02-06 | 浙江伽奈维医疗科技有限公司 | 一种基于rgbd相机的呼吸监测装置及其呼吸监测方法 |
CN114332147B (zh) * | 2021-12-30 | 2023-01-24 | 合肥的卢深视科技有限公司 | 散斑图预处理的方法、装置、电子设备和存储介质 |
CN114596279B (zh) * | 2022-03-08 | 2023-09-22 | 江苏省人民医院(南京医科大学第一附属医院) | 一种基于计算机视觉的非接触式呼吸检测方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103824057A (zh) * | 2014-02-18 | 2014-05-28 | 江苏大学 | 一种基于面积特征算子的猪的呼吸频率检测方法 |
JP2015229101A (ja) * | 2014-06-03 | 2015-12-21 | 株式会社イデアクエスト | 呼吸運動測定装置 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8792969B2 (en) * | 2012-11-19 | 2014-07-29 | Xerox Corporation | Respiratory function estimation from a 2D monocular video |
FR3023699B1 (fr) * | 2014-07-21 | 2016-09-02 | Withings | Procede et dispositif de surveillance d'un bebe et d'interaction |
-
2017
- 2017-06-30 CN CN201710519852.5A patent/CN107260173B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103824057A (zh) * | 2014-02-18 | 2014-05-28 | 江苏大学 | 一种基于面积特征算子的猪的呼吸频率检测方法 |
JP2015229101A (ja) * | 2014-06-03 | 2015-12-21 | 株式会社イデアクエスト | 呼吸運動測定装置 |
CN105283128A (zh) * | 2014-06-03 | 2016-01-27 | 创探科技股份公司 | 呼吸运动测定装置 |
Also Published As
Publication number | Publication date |
---|---|
CN107260173A (zh) | 2017-10-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107260173B (zh) | 一种基于摄像头与球状标记物的呼吸测量方法 | |
US10506952B2 (en) | Motion monitor | |
US11589776B2 (en) | Non-contact breathing activity monitoring and analyzing through thermal and CO2 imaging | |
US11712176B2 (en) | Systems and methods for video-based non-contact tidal volume monitoring | |
Hu et al. | Synergetic use of thermal and visible imaging techniques for contactless and unobtrusive breathing measurement | |
JP6371837B2 (ja) | 被験者のバイタルサインを取得するデバイスと方法 | |
US11776146B2 (en) | Edge handling methods for associated depth sensing camera devices, systems, and methods | |
WO2018170009A1 (en) | Imaging-based spirometry systems and methods | |
CN105869144B (zh) | 一种基于深度图像数据的非接触式呼吸监测方法 | |
CN105007806A (zh) | 用于确定来自对象的生命体征的装置和方法 | |
Chatterjee et al. | Real-time respiration rate measurement from thoracoabdominal movement with a consumer grade camera | |
Lukáč et al. | Contactless recognition of respiration phases using web camera | |
Gwak et al. | Motion-based respiratory rate estimation with motion artifact removal using video of face and upper body | |
Tan et al. | Lightweight video-based respiration rate detection algorithm: An application case on intensive care | |
CN105046209B (zh) | 一种基于典型相关分析的非接触式心率测量方法 | |
CN113456061A (zh) | 一种基于无线信号睡眠姿势监测方法及系统 | |
Chatterjee et al. | Real-time visual respiration rate estimation with dynamic scene adaptation | |
Chiang et al. | The non-contact respiratory monitoring system using thermal image processing | |
US20230329590A1 (en) | Non-Contact Monitoring System and Method | |
Ta-Chi et al. | The Non-Contact Respiratory Monitoring System Using Thermal Image Processing | |
Khan et al. | Joint use of a low thermal resolution thermal camera and an RGB camera for respiration measurement | |
WO2023199193A1 (en) | Non-contact monitoring system and method | |
Chatterjee et al. | Real-time Visual Respiration Rate Estimation with Automatic Scene Adaptation | |
Fatema-Tuz-Zohra Khanam et al. | Integrating RGB-Thermal image sensors for non-contact automatic respiration |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |