CN109419497A - 基于热成像的关脉识别方法 - Google Patents
基于热成像的关脉识别方法 Download PDFInfo
- Publication number
- CN109419497A CN109419497A CN201710872547.4A CN201710872547A CN109419497A CN 109419497 A CN109419497 A CN 109419497A CN 201710872547 A CN201710872547 A CN 201710872547A CN 109419497 A CN109419497 A CN 109419497A
- Authority
- CN
- China
- Prior art keywords
- image
- wrist
- pixel
- radial artery
- edge
- 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.)
- Granted
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/01—Measuring temperature of body parts ; Diagnostic temperature sensing, e.g. for malignant or inflamed tissue
- A61B5/015—By temperature mapping of body part
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/01—Measuring temperature of body parts ; Diagnostic temperature sensing, e.g. for malignant or inflamed tissue
-
- 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
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/48—Other medical applications
- A61B5/4854—Diagnosis based on concepts of traditional oriental medicine
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/48—Other medical applications
- A61B5/4887—Locating particular structures in or on the body
- A61B5/489—Blood vessels
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61H—PHYSICAL THERAPY APPARATUS, e.g. DEVICES FOR LOCATING OR STIMULATING REFLEX POINTS IN THE BODY; ARTIFICIAL RESPIRATION; MASSAGE; BATHING DEVICES FOR SPECIAL THERAPEUTIC OR HYGIENIC PURPOSES OR SPECIFIC PARTS OF THE BODY
- A61H39/00—Devices for locating or stimulating specific reflex points of the body for physical therapy, e.g. acupuncture
- A61H39/02—Devices for locating such points
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Medical Informatics (AREA)
- Veterinary Medicine (AREA)
- Public Health (AREA)
- Engineering & Computer Science (AREA)
- General Health & Medical Sciences (AREA)
- Animal Behavior & Ethology (AREA)
- Biomedical Technology (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Heart & Thoracic Surgery (AREA)
- Physics & Mathematics (AREA)
- Pathology (AREA)
- Biophysics (AREA)
- Rehabilitation Therapy (AREA)
- Physiology (AREA)
- Cardiology (AREA)
- Alternative & Traditional Medicine (AREA)
- Epidemiology (AREA)
- Pain & Pain Management (AREA)
- Physical Education & Sports Medicine (AREA)
- Vascular Medicine (AREA)
- Image Analysis (AREA)
- Image Processing (AREA)
Abstract
一种基于热成像的关脉识别方法,包括:获取被测者小臂和腕部的热成像影像,并通过边缘检测算法提取所述热成像影像中小臂和手腕的轮廓线;对得到的轮廓线进行预处理,提取桡骨突茎特征点的x坐标;对桡动脉进行图像分割,将桡动脉热成像图拟合成一条直线函数,该直线函数中桡骨突茎特征点所在x坐标和对应的y坐标就是关脉最终位置。本发明的方法首次提出了能够通过热成像图像处理自动识别关脉位置的算法,将识别关脉客观化、数字化,数据具有较好重复性和再现性,比超声阵列和其他类型传感器阵列接触式判别关脉方法实现上更有可行性。
Description
技术领域
本发明涉及脉络检测技术领域,尤其是涉及一种基于热成像的关脉识 别方法。
背景技术
在中医脉诊诊断过程中,关脉的定位极为重要。能否找准关脉是能否 准确判断病人疾病的前提。中医脉诊中经典的判断关脉的生物特征点是桡 骨突茎(中医叫掌后高骨)和腕横纹,经典古籍中认为桡骨突茎上的桡动 脉即为关脉。另一部分古籍描述了以腕横纹为原点的情况下,关脉与腕横 纹之间的距离关系。尽管有着两个明显的特征点,几千年来,中医医生的 诊断仍然只是依靠手的触感来判断关脉位置,由于不同人的感觉差异,造 成关脉判别的主观性过强,关脉判别的再现性较差,相同人对关脉定位重 复性也较差,无法保证判别精度。现存的中医自动化诊断仪器,关脉定位 步骤也都是由人判断完成后,再将脉搏波探头放入寸关尺对应位置,关脉 定位的客观性无法保证。
发明内容
有鉴于此,本发明的主要目的在于提供一种关脉识别算法,以解决上 述技术问题中的至少之一。
为了实现上述目的,本发明提供了一种关脉识别方法,包括以下步骤:
步骤S1,获取被测者小臂和腕部的轮廓线;
步骤S2,对步骤S1得到的轮廓线进行预处理,提取桡骨突茎特征点 的x坐标;
步骤S3,获取被测者小臂和腕部的热成像图,将热成像图中的桡动脉 影像拟合成一条直线函数,所述桡骨突茎特征点的x坐标和该x坐标在所 述直线函数上对应的y坐标就是关脉最终位置。
其中,步骤S1具体包括:
获取被测者小臂和腕部的热成像影像,并通过边缘检测算法提取所 述热成像影像中小臂和手腕的轮廓线。
其中,步骤S2中对步骤S1得到的轮廓线进行预处理的步骤具体包括:
步骤S21,对步骤S1得到的轮廓线进行连通域识别,判断所述轮 廓线中的最大连通域是否贯穿所述红外影像左右两侧的边界,即连通域是 否存在断点,如果不存在断点,则所述最大连通域为所要提取的所述小臂 和腕部的轮廓线,跳转到步骤S23;
步骤S22,以最大连通域左侧断点为原点对该点上、左上、左、左 下、下5个方向2个像素范围内寻找边缘片段,若在搜寻范围内存在其他 连通域,则将两个连通域相连,中间断点像素通过插值在两个片断之间补 像素点,最终形成一个新的连通域,并以新的连通域左侧断点为原点进一 步寻找其他边缘片段,直至到达图像左侧边界;
以最大连通域右侧断点为原点对该点上、右上、右、右 下、下5个方向2个像素范围内寻找边缘片段,若在搜寻范围内存在其他 连通域,则将两个连通域相连,中间断点像素通过插值在两个片断之间补 像素点,最终形成一个新的连通域,并以新的连通域右侧断点为原点进一 步寻找其他边缘片段,直至到达图像右侧边界;
其中,上述向左、向右的搜寻和接续断点的步骤不分先 后;
步骤S23,将基于上述最大连通域得到的二维图像的边缘轮廓线转 为一维曲线,消除转换过程中产生的阶跃点,使转换的一维曲线更平滑, 突出小臂和腕部的边缘特征。
其中,步骤S2中所述提取桡骨突茎特征点的步骤具体包括:
步骤S24,对转换的小臂和腕部边缘轮廓线的一维曲线进行特征提 取,寻找手部与桡骨突茎之间的最低凹陷处;
步骤S25,以所述最低凹陷处为基准,在该基准左侧0~3cm范围内, 在轮廓线对应的曲率波形中寻找第一个波峰,如果存在,则该波峰即识别 为关脉的x坐标,如果不存在,则所述最低凹陷处即识别为关脉的x坐标。
其中,在步骤S24的寻找手部与桡骨突茎之间的最低凹陷处的步骤中, 是基于桡骨突茎在手腕凹陷到小臂的边缘轮廓线上存在一个曲率最大点。
其中,步骤S3中将热成像图中的桡动脉影像拟合成一条直线函数的 步骤具体包括:
步骤S31,对所述热成像图构建桡动脉区域,设置阈值,用于为二 值化所述桡动脉区域提供阈值参考;
步骤S32,二值化所述桡动脉区域,将桡动脉的图像与热成像图中 的其它部分分离;
步骤S33,将步骤S32中得到的桡动脉图像拟合成直线函数。
其中,步骤S31中对热成像图构建桡动脉区域,设置阈值的步骤具体 包括:
用步骤S2中预处理后生成的轮廓线上的每一个像素为原点构建一 个区域;
根据桡动脉边界位置的区域像素均值和方差的统计规律,设定均值 和方差的阈值。
其中,步骤S32中的二值化所述桡动脉区域的步骤具体包括:
计算每个边缘像素区域中像素的均值和方差;
将生成的每个边缘像素区域中像素的均值和方差逐次与阈值作比 较,二值化符合阈值条件的区域。
其中,步骤S33中的将步骤S32中得到的桡动脉图像拟合成直线函数 的步骤具体包括:
对二值化的桡动脉图像像素纵坐标求平均,获得描述桡动脉图像的 曲线;
对曲线进行直线拟合,得到包含桡动脉走势的直线函数。
基于上述技术方案可知,本发明的关脉识别方法相对于现有技术具有 如下有益效果:
1、依据的是中医认可的两个生理特征,为中医脉诊提供了新的关脉 识别方法;
2、首次提出了能够通过热成像图像处理自动识别关脉位置的算法;
3、首次将识别关脉客观化、数字化,数据具有较好重复性和再现性;
4、在关脉判断上,消除了人的主观干扰因素;
5、非接触式的关脉判别方法,比超声阵列和其他类型传感器阵列接 触式判别关脉方法实现上更有可行性;
6、其中的边缘检测算法可以实现任意角度边缘检测;
7、首次提出了桡骨突茎判别方法。
附图说明
图1是携带桡动脉信息的手臂腕部的红外热成像图;
图2是手臂边缘断点连接示意图;
图3是手臂和腕部边缘图像;
图4是转化为一维曲线的边缘和经过滤波或高阶多项式拟合过的手臂 腕部曲线;
图5是手臂腕部边缘与对应的曲率曲线图;
图6是带有桡动脉信息的手臂腕部边缘图像;
图7是分割出的桡动脉图像;
图8是桡动脉像素纵坐标平均化和直线拟合曲线;
图9是桡动脉的坐标显示图;
图10是本发明的红外热成像的非接触式手腕成像平台的结构示意图;
图11是腕横纹与关脉x坐标的关系示意图;
图12是图像相邻像素紧凑连接的示意图;
图13是图像相邻像素宽松连接的示意图;
图14是图像边缘检测任意角度组成形式的示意图;
图15是超出图像边界的部分对图像进行补0扩增的示意图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白,以下结合具体实 施例,并参照附图,对本发明作进一步的详细说明。
本发明的目的是提供一种能够去除人主观性差异的关脉自动识别方 法,该方法通过图像边缘识别技术结合桡动脉热成像特点,对人体手腕关 脉进行非接触式识别。在对手臂腕部热成像时,红外热成像光谱范围远于 可见光谱范围,可见光谱范围内的图像均不显示,生成的图像有效的简化 图像背景的复杂度。因为桡动脉温度高于其他皮肤温度,所以热成像光谱 可以凸显桡动脉,为桡动脉图像分割和直线拟合提供了图像基础,如图1 所示。
为实现上述目的,本发明公开了一种关脉识别方法,包括以下步骤:
1、边缘检测算法提取手腕轮廓;
2、对手腕边缘进行预处理,提取桡骨突茎特征点x坐标;
3、对桡动脉进行图像分割,将桡动脉热成像图拟合成一条直线函数, 该直线函数中桡骨突茎特征点所在x坐标和对应的y坐标就是关脉最终位 置。
更具体地,本发明的关脉识别方法包括以下步骤:
1、对待检测的手臂和腕部的边缘进行识别,生成手臂和腕部的边缘 线条。识别手臂和腕部边缘的算法可以为后文所述的本申请的任意角度边 缘检测算法,或者也可以是现有的其它类型的边缘检测算法。其中,为了 描述的方便,拍摄时以左手为例进行描述,将手腕置于照片的右侧,连接 肘部的手臂位于照片的左侧,以此建立上下左右的平面方位坐标。但需要 说明的是,这仅是为了描述的方便,并不用于限制本发明。
2、对手臂和腕部的边缘预处理,将手臂和腕部的边缘进一步优化, 为后续腕部关脉的识别提供保障。该步骤具体包括识别手臂边缘最大连通 域、手臂边缘断点连接、手臂腕部曲线拟合,如下所示:
(1)识别手臂边缘最大连通域:
对生成的边缘图像进行连通域识别,找出图像右侧边界最大连通域。 若最大连通域贯穿图像左右两侧边界,即连通域不存在断点,该最大连通 域即可认为是手臂腕部边缘。
(2)手臂边缘断点连接:
将手臂边缘片段连接起来,形成一个贯穿图像左右边界的手臂腕部 整体边缘。在边缘存在断点情况下,最大连通域只是手臂腕部边缘的一部 分,因此需要将其他手臂腕部边缘片段连接起来。以最大连通域左侧断点 为原点对该点上、左上、左、左下、下5个方向2个像素范围内寻找边缘 片段。若在搜寻范围内存在其他连通域,则将两个连通域相连,中间断点 像素通过插值或者其它拟合方式在两个片断之间补像素点,最终形成一个 新的连通域,并以新的连通域左侧断点为原点进一步寻找其他边缘片段, 直至到达图像左侧边界;此外,也包括以最大连通域右侧断点为原点对该 点上、右上、右、右下、下5个方向2个像素范围内寻找边缘片段,若在 搜寻范围内存在其他连通域,则将两个连通域相连,中间断点像素通过插 值在两个片断之间补像素点,最终形成一个新的连通域,并以新的连通域右侧断点为原点进一步寻找其他边缘片段,直至到达图像右侧边界;上述 向左、向右的搜寻和接续断点的步骤不分先后。
(3)手臂腕部曲线拟合:
消除二维图像边缘转为一维曲线过程中产生的阶跃点,使转换的一 维手臂边缘曲线更平滑,突出手臂腕部边缘特征。
3、识别桡骨突茎算法用于识别桡骨突茎特征点:
首先对提取的手臂腕部边缘进行特征提取,识别手部与桡骨突茎之间 的凹陷处,寻找凹陷处的最低点。桡骨突茎在表皮顶部曲率变化特点是手 腕凹陷到其左侧两个手指的距离范围(大概是0~3cm或0~4cm)内存在 几个曲率局部峰值点,即边界变化弯曲幅度较大的点。其次,寻找该距离 里离凹陷处最远的那个曲率局部峰值点(也就是几个曲率局部峰值点里的 第一个峰值点)。最后,在该曲率局部峰值点识别为关脉x坐标;如果没 有峰,则该凹陷处识别为关脉x坐标。
4、桡动脉图像分割和关脉识别用于分割桡动脉图像并拟合成能够反 映桡动脉走势的直线函数,具体步骤包括:
(1)区域构建和阈值设置,用于为二值化桡动脉提供阈值参考;
(2)二值化桡动脉区域,用于将桡动脉图像和其他图像分离;
(3)桡动脉直线拟合用于获得反映桡动脉走势的直线函数和最终的 关脉坐标。
在一个具体实施方式中,本发明的关脉识别方法包括以下步骤:
首先对整个图像进行边缘识别,生成手臂腕部边缘。
之后对手臂腕部边缘进行预处理,进一步优化手臂腕部边缘,为后续 腕部关脉识别提供保障。预处理过程包括识别手臂边缘最大连通域、手臂 边缘断点连接、手臂腕部曲线拟合。
(1)识别手臂边缘最大连通域:对生成的边缘图像进行连通域识别, 找出图像右侧边界最大连通域。若最大连通域贯穿图像左右两侧边界,即 连通域不存在断点,该最大连通域即可认为是手臂腕部边缘。
(2)如图2所示,手臂边缘断点连接包括以下步骤:将手臂边缘片 段连接起来,形成一个贯穿图像左右边界的手臂腕部整体边缘。在边缘存 在断点情况下,最大连通域只是手臂腕部边缘的一部分,因此需要将其他 手臂腕部边缘片段连接起来。以最大连通域左侧断点为原点对该点上、左 上、左、左下、下5个方向2个像素范围内寻找边缘片段。若在搜寻范围 内存在其他连通域,则将两个连通域相连,中间断点像素通过插值或者其 它拟合方式在两个片断之间补像素点,最终形成一个新的连通域,并以新 的连通域左侧断点为原点进一步寻找其他边缘片段,直至到达图像左侧边 界。
(3)如图3和图4所示,手臂腕部曲线拟合包括以下步骤:用低通 滤波器或者多项式曲线拟合消除二维图像边缘转为一维曲线过程中产生 的阶跃点,使转换的一维手臂边缘曲线更平滑,突出手臂腕部边缘特征。
识别桡骨突茎算法用于识别桡骨突茎特征点。如图5所示,首先对提 取的手臂腕部边缘进行特征提取,识别手部与桡骨突茎之间的凹陷处,寻 找凹陷处的最低点。桡骨突茎在表皮顶部曲率变化特点是手腕凹陷到其左 侧两个手指该距离范围内存在几个曲率局部峰值点,即边界变化弯曲幅度 较大的点。其次,寻找该距离里离凹陷处最远的那个曲率局部峰值点(也 就是几个曲率局部峰值点里的第一个峰值点)。最后,在该曲率局部峰值点识别为关脉x坐标。
接着进行桡动脉图像分割和关脉识别。用之前生成的边缘图像(图6) 中每一个像素为原点构建一个区域。根据桡动脉边界位置的区域像素均值 和方差的统计规律,设定均值和方差的阈值。计算每个边缘像素区域中像 素的均值和方差。将生成的每个边缘像素区域中像素的均值和方差逐次与 阈值作比较,二值化符合阈值条件的区域(图7)。对二值化的桡动脉图像 像素纵坐标求平均,获得描述桡动脉图像的曲线。对曲线进行二次多项式直线拟合,得到包含桡动脉走势的直线函数(图8),将关脉x坐标代入该 直线函数得到关脉的纵坐标。关脉在图像中位置即可确定,如图9所示。
本发明还公开了一种基于红外热成像的非接触式关脉自动识别系统, 包括基于红外热成像的非接触式手腕成像平台和关脉提取算法单元;其中, 该基于红外热成像的非接触式手腕成像平台的固定架结构能够保证测量 人每次手腕以固定姿势进入成像区域,平台的一字激光标志提供了手腕腕 横纹对准参照位置;该基于红外热成像的非接触式手腕成像平台,将携带 桡动脉图像信息的手腕热成像图像传送给关脉提取单元,关脉提取单元采 用上述的关脉识别方法步骤提取关脉。
具体地,该基于红外热成像的非接触式手腕成像平台包括三维移动台、 红外热成像仪、点状激光、一字线激光、手与腕部固定架、仪器支架和显 示器。其中仪器支架与三维移动台、红外热成像仪、点状激光、一字线激 光、手与腕部固定架一侧相连,用于支持整个系统。
三维移动台另一侧与点状激光一侧相连用于搭载点状激光指示手腕 关脉位置,其中三维移动台控制部分可以与计算机通信获得关脉坐标,将 点状激光移动到关脉坐标处。三维移动台还可以连接脉搏波传感器,用于 将脉搏波传感器送到关脉位置采集脉搏波。红外热成像仪一侧与仪器支架 相连用于拍摄携带桡动脉信息的手腕图像。一字线激光一侧与仪器支架相 连,用于为腕横纹提供对准位置。手与腕部固定架一侧与仪器支架连接用 于保证被测者每次测试手型以及手臂位置基本不变。计算机用于计算和显 示最终关脉位置并传送坐标信号给三维移动平台。
在一个优选实施例中,本发明的基于红外热成像的非接触式关脉自动 识别系统如图10所示,被试者手腕放在手与腕部固定架上并将腕横纹与 一字线激光对齐。手与腕部固定架包含了一个把手状或指套装置,被试者 在每次测试中需要握住把手状装置,该装置可使被试者手部和腕部在每次 测试中保持相同姿势,并在图像中保持相同的位置。一字线激光投射位置 为热成像图像右侧边界,这个热成像图像坐标是最左侧为原点,腕横纹对应的是图像右侧最大坐标,被测者手腕腕横纹与一字线激光对齐,手腕腕 横纹即与图像右侧边界基本对应,如图11所示。启动关脉识别程序,热 成像仪采集腕部图像和桡动脉图像。由于红外热成像光谱范围远于可见光 谱范围,可见光谱范围内的图像均不显示,有效的简化图像背景的复杂度。 另外因为桡动脉温度高于其他皮肤温度,所以热成像光谱可以对桡动脉进 行成像,为桡动脉图像分割和直线拟合提供了图像基础,如图10所示。 之后热成像仪将携带桡动脉热成像信息的腕部图像传送给计算机。计算机 的关脉提取单元对图像中关脉位置进行识别,并将生成的关脉坐标传送给 三维移动台的控制部分,控制部分再控制三维移动平台携带点激光移动到 图像关脉坐标对应的实际位置上,点状激光所指示的位置就是计算出的关 脉实际位置。采用该方法计算出在图像中关脉的x坐标与图像右侧边界坐 标的差值就是以腕横纹为原点距离关脉的距离了。另外,该平台有效保证 了在该系统的测试结果的重复性。
对于识别手臂和腕部边缘的算法,本发明还公开了一种优选的任意角 度的边缘检测方法,其包括以下步骤:
(1)构建边缘检测角度区间边界
将图像的相邻像素点关系分为紧凑连接和宽松连接,以每行两个像素 为例,紧凑连接如图12所示:从图像最左侧列的像素为起点,相邻行的 像素点首尾像素垂直相连,每两行组成一个紧凑连接单元。根据这种方式 若干紧凑连接单元连接成一条线直到图像边界,这条直线与它的y轴方向 投影的夹角就是边缘检测角度。其矩阵Qθ2L表示方式为:
其边缘检测角度
宽松连接如图13所示:从图像左上角顶点起,以最左侧列和最上侧 行的像素为起点,相邻行的像素点首尾像素对角相连,每两行组成一个宽 松连接单元。根据这种方式若干宽松连接单元连接成一条线直到图像边界, 这条直线与它的y轴方向投影的夹角就是边缘检测方向。其矩阵Qθ2R表示 方式为:
其边缘检测角度
因此两像素的紧凑连接单元组成的边缘检测角度为该 段角度区间的左边界。两像素的宽松连接单元组成的边缘检测角度 为该段角度区间的右边界。所以该角度区间为(θ2L,θ2R)。
当像素个数为k个像素时,其矩阵QθkL表示方式为:
k个像素的紧凑连接单元组成的边缘检测角度区间的左边界
其矩阵QθkR表示方式为:
k个像素的紧凑连接单元组成的边缘检测角度区间的右边界
因此,检测角度区间边界角度为θ2L,θ2R,θ3L,θ3R,…θnL,θnR,该角度边界形 成的角度区间为若干子区间的并集(θ2L,θ2R)∪(θ3L,θ3R)∪…∪(θnL,θnR)。该并集 的范围为(45°,90°)。
(2)构建边缘检测角度区间中的任意角度
以两像素相连的单元为例,区间中的任意角度组成形式如下:
如图14所示,明确需要进行边缘检测的角度,构建一个单元。该单 元由i个紧凑连接和j个宽松连接组成,其中i和j的个数与检测的角度对 应。该单元重复出现r次直至超过图像边界,其行数m和列数n与i、j 和r关系为:
r(i+j)+1=m (1)
ri(k-1)+krj+k=n (2)
因此每个角度区间里的边缘检测角度
另外,每个边界条件也符合上式。
所以通过上述方法能够将图像中的像素根据所需角度进行组合,对于 算法实现超出图像边界的部分对图像进行补0扩增,如图15所示。
(3)将按上述方法构建的若干像素直线分别与高斯函数的一阶导数 fσ(t)作卷积运算,并对卷积运算结果取绝对值,并对绝对值取局部极大值; 以紧凑连接生成的边缘检测上界角度为例解释若干像素直线,以左侧边界 为起点生成若干像素直线为X1、X2…Xm,以上侧边界为起点生成的若干像 素直线为Y1…Ym-1。其中m为行,k为连接像素个数。
以宽松连接生成的边缘检测下界角度为例解释若干像素直线,以左侧 边界为起点生成若干像素直线为X′1、X′2…X′m,以上侧边界为起点生成的若 干像素直线为Y′1…Y′m-1其中m为行,k为连接像素个数。
每一个像素直线都分别与高斯函数的一阶导数fσ(t)作卷积运算,并对 卷积运算结果取绝对值得到:|fσ(t)*X1|,|fσ(t)*X2|,…|fσ(t)*Xm| |fσ(t)*X′1|,|fσ(t)*X′2|,…|fσ(t)*X′m|和|fσ(t)*Y1|,…|fσ(t)*Ym-1||fσ(t)*Y′1|,…|fσ(t)*Y′m-1|。 通过对构建的若干像素直线作卷积和取绝对值的运算,使边缘检测角度从 [0°,360°]缩减到[0°,180°]。因此只需对边缘检测角度[0°,180°]的区间里对图像 进行处理。
(4)对得到的|fσ(t)*X1|,|fσ(t)*X2|,…|fσ(t)*Xm|和|fσ(t)*Y1|,…|fσ(t)*Ym-1|进 行局部极大值运算赋灰度值,其他非局部极大值像素灰度设为0,其灰度 值为(255/边缘检测角度个数)。根据像素下标将具有灰度值的图像像素替 换到原图像中相同像素下标位置上;
(5)将不同边缘检测角度方向得到的若干图像进行灰度叠加,根据 实际所需边缘图像要求对多次叠加后图像的灰度设二值化阈值,根据该二 值化阈值对图像进行二值化处理。(这里的阈值也是没有固定要求的)最 终得到所需边缘。
上述边缘检测角度范围为(45°,90°),45°边缘检测角度就是一个像素 依次连接组成的像素直线,即k=1时候。90°方向就是垂直分割图像,每 列像素分别组成像素直线。因此,该检测角度范围[45°,90°]可以实现。
通过将图像矩阵转置和翻转可以将角度范围为[45°,90°]映射到[0°,45°],[90°,135°]和[135°,180°]。具体方法如下:
将图像矩阵水平翻转,边缘检测角度区间为从[45°,90°]映射为[90°,135°]。 将图像矩阵转置后,边缘检测角度区间为从[45°,90°]映射为[135°,180°]。将 图像矩阵水平翻转和转置后边缘检测角度区间为从[45°,90°]映射为[0°,45°]。 基于以上方法,实现[0°,360°]角度区间的边缘检测只需应用[45°,90°]边缘检 测角度即可实现。
以上所述的具体实施例,对本发明的目的、技术方案和有益效果进行 了进一步详细说明,应理解的是,以上所述仅为本发明的具体实施例而已, 并不用于限制本发明,凡在本发明的精神和原则之内,所做的任何修改、 等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (9)
1.一种基于热成像的关脉识别方法,包括以下步骤:
步骤S1,获取被测者小臂和腕部的轮廓线;
步骤S2,对步骤S1得到的轮廓线进行预处理,提取桡骨突茎特征点的x坐标;
步骤S3,获取被测者小臂和腕部的热成像图,将热成像图中的桡动脉影像拟合成一条直线函数,所述桡骨突茎特征点的x坐标和该x坐标在所述直线函数上对应的y坐标就是关脉最终位置。
2.根据权利要求1所述的方法,其特征在于,步骤S1具体包括:
获取被测者小臂和腕部的热成像影像,并通过边缘检测算法提取所述热成像影像中小臂和手腕的轮廓线。
3.根据权利要求1所述的方法,其特征在于,步骤S2中对步骤S1得到的轮廓线进行预处理的步骤具体包括:
步骤S21,对步骤S1得到的轮廓线进行连通域识别,判断所述轮廓线中的最大连通域是否贯穿所述红外影像左右两侧的边界,即连通域是否存在断点,如果不存在断点,则所述最大连通域为所要提取的所述小臂和腕部的轮廓线,跳转到步骤S23;
步骤S22,以最大连通域左侧断点为原点对该点上、左上、左、左下、下5个方向2个像素范围内寻找边缘片段,若在搜寻范围内存在其他连通域,则将两个连通域相连,中间断点像素通过插值在两个片断之间补像素点,最终形成一个新的连通域,并以新的连通域左侧断点为原点进一步寻找其他边缘片段,直至到达图像左侧边界;
以最大连通域右侧断点为原点对该点上、右上、右、右下、下5个方向2个像素范围内寻找边缘片段,若在搜寻范围内存在其他连通域,则将两个连通域相连,中间断点像素通过插值在两个片断之间补像素点,最终形成一个新的连通域,并以新的连通域右侧断点为原点进一步寻找其他边缘片段,直至到达图像右侧边界;
其中,上述向左、向右的搜寻和接续断点的步骤不分先后;
步骤S23,将基于上述最大连通域得到的二维图像的边缘轮廓线转为一维曲线,消除转换过程中产生的阶跃点,使转换的一维曲线更平滑,突出小臂和腕部的边缘特征。
4.根据权利要求3所述的方法,其特征在于,步骤S2中所述提取桡骨突茎特征点的步骤具体包括:
步骤S24,对转换的小臂和腕部边缘轮廓线的一维曲线进行特征提取,寻找手部与桡骨突茎之间的最低凹陷处;
步骤S25,以所述最低凹陷处为基准,在该基准左侧0~4cm范围内,在轮廓线对应的曲率波形中寻找第一个波峰,如果存在,则该波峰即识别为关脉的x坐标,如果不存在,则所述最低凹陷处即识别为关脉的x坐标。
5.根据权利要求4所述的方法,其特征在于,在步骤S24的寻找手部与桡骨突茎之间的最低凹陷处的步骤中,是基于桡骨突茎在手腕凹陷到小臂的边缘轮廓线上存在一个边界变化弯曲幅度最大点来实现的。
6.根据权利要求1所述的方法,其特征在于,步骤S3中将热成像图中的桡动脉影像拟合成一条直线函数的步骤具体包括:
步骤S31,对所述热成像图构建桡动脉区域,设置阈值,用于为二值化所述桡动脉区域提供阈值参考;
步骤S32,二值化所述桡动脉区域,将桡动脉的图像与热成像图中的其它部分分离;
步骤S33,将步骤S32中得到的桡动脉图像拟合成直线函数。
7.根据权利要求6所述的方法,其特征在于,步骤S31中对热成像图构建桡动脉区域,设置阈值的步骤具体包括:
用步骤S2中预处理后生成的轮廓线上的每一个像素为原点构建一个区域;
根据桡动脉边界位置的区域像素均值和方差的统计规律,设定均值和方差的阈值。
8.根据权利要求7所述的方法,其特征在于,步骤S32中的二值化所述桡动脉区域的步骤具体包括:
计算每个边缘像素区域中像素的均值和方差;
将生成的每个边缘像素区域中像素的均值和方差逐次与阈值作比较,二值化符合阈值条件的区域。
9.根据权利要求7所述的方法,其特征在于,步骤S33中的将步骤S32中得到的桡动脉图像拟合成直线函数的步骤具体包括:
对二值化的桡动脉图像像素纵坐标求平均,获得描述桡动脉图像的曲线;
对曲线进行直线拟合,得到包含桡动脉走势的直线函数。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710778521 | 2017-08-31 | ||
CN2017107785213 | 2017-08-31 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109419497A true CN109419497A (zh) | 2019-03-05 |
CN109419497B CN109419497B (zh) | 2021-12-17 |
Family
ID=65513656
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710872547.4A Active CN109419497B (zh) | 2017-08-31 | 2017-09-25 | 基于热成像的关脉识别方法 |
Country Status (2)
Country | Link |
---|---|
CN (1) | CN109419497B (zh) |
WO (1) | WO2019041412A1 (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111297344A (zh) * | 2020-03-24 | 2020-06-19 | 陈恬慧 | 一种脉诊视觉阵列采集方法 |
CN111803071A (zh) * | 2020-06-23 | 2020-10-23 | 山东大学 | 耳朵外观轮廓三维扫描系统及方法 |
CN113155288A (zh) * | 2020-11-30 | 2021-07-23 | 齐鲁工业大学 | 一种光伏电池热斑的图像识别方法 |
CN113171062A (zh) * | 2021-04-29 | 2021-07-27 | 中国科学院微电子研究所 | 识别寸、关、尺的方法、终端、系统、介质及计算机设备 |
CN113303771A (zh) * | 2021-07-30 | 2021-08-27 | 天津慧医谷科技有限公司 | 脉搏采集点的确定方法、装置和电子设备 |
US11346938B2 (en) | 2019-03-15 | 2022-05-31 | Msa Technology, Llc | Safety device for providing output to an individual associated with a hazardous environment |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110781791B (zh) * | 2019-10-18 | 2023-03-24 | 河南大学 | 一种基于重心对齐的步态能量图优化合成方法 |
GB202005213D0 (en) | 2020-04-08 | 2020-05-20 | Bevan Heba | Remote temperature measurement system |
CN111951152B (zh) * | 2020-08-10 | 2023-07-25 | 河南大学 | 一种兼顾原边连续性和断边碴口匹配度的甲骨缀合方法 |
CN112101365A (zh) * | 2020-09-10 | 2020-12-18 | 国网辽宁省电力有限公司电力科学研究院 | 基于红外热像处理的电力设备关键特征提取方法及系统 |
CN113469185B (zh) * | 2021-06-04 | 2023-03-24 | 四川轻化工大学 | 一种变压器油枕的油位高度检测方法 |
CN113850800B (zh) * | 2021-10-15 | 2024-04-30 | 郑州磨料磨具磨削研究所有限公司 | 一种硬脆材料划切缝崩边检测方法 |
CN115376694B (zh) * | 2022-10-24 | 2023-02-28 | 佛山科学技术学院 | 一种肢体功能障碍的多模态自适应评估系统及方法 |
CN115797356B (zh) * | 2023-02-09 | 2023-05-23 | 山东第一医科大学附属省立医院(山东省立医院) | 一种核磁共振肿瘤区域提取方法 |
Citations (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20080045847A1 (en) * | 2006-06-30 | 2008-02-21 | University Of Louisville Research Foundation, Inc. | Non-contact and passive measurement of arterial pulse through thermal IR imaging, and analysis of thermal IR imagery |
CN101551851A (zh) * | 2008-03-31 | 2009-10-07 | 中国科学院沈阳自动化研究所 | 一种红外图像目标识别方法 |
CN102446349A (zh) * | 2010-09-29 | 2012-05-09 | 奥林巴斯株式会社 | 图像处理装置和图像处理方法 |
CN103006195A (zh) * | 2013-01-10 | 2013-04-03 | 浙江大学 | 一种基于图像处理的非接触式生命体征数据监测系统和监测方法 |
CN104027094A (zh) * | 2014-05-30 | 2014-09-10 | 中国中医科学院医学实验中心 | 一种用于中医脉诊信息采集的寸、关、尺定位方法 |
CN104268853A (zh) * | 2014-03-06 | 2015-01-07 | 上海大学 | 一种红外图像与可见光图像配准方法 |
CN104484667A (zh) * | 2014-12-30 | 2015-04-01 | 华中科技大学 | 一种基于亮度特征和轮廓完整性的轮廓提取方法 |
US20150133791A1 (en) * | 2012-05-29 | 2015-05-14 | National University Corporation Kochi University | Artery visualization device and artery imaging device |
CN104799815A (zh) * | 2015-03-26 | 2015-07-29 | 山东大学 | 一种基于图像引导的血气采集装置及方法 |
CN105147261A (zh) * | 2015-08-03 | 2015-12-16 | 刘垚 | 一种中医脉诊仪及其定位寸、关、尺脉搏点的方法 |
WO2016018463A2 (en) * | 2013-12-23 | 2016-02-04 | Jason Case | Methods and systems using led sources to create contrast agents for mid-infrared imaging of blood vessels |
CN105469414A (zh) * | 2015-12-15 | 2016-04-06 | 凌云光技术集团有限责任公司 | 轮廓连接方法及装置 |
CN106456056A (zh) * | 2015-05-13 | 2017-02-22 | 普勒斯麦公司 | 动脉可视化装置 |
CN106859956A (zh) * | 2017-01-13 | 2017-06-20 | 北京奇虎科技有限公司 | 一种人体穴位识别按摩方法、装置及ar设备 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE69837526T9 (de) * | 1997-11-19 | 2008-04-10 | Seiko Epson Corp. | Verfahren und vorrichtung zum nachweis von pulswellen sowie verfahren zur anzeige der lage von arterien |
CN101773383B (zh) * | 2010-03-09 | 2011-05-25 | 哈尔滨工业大学 | 桡动脉位置检测方法 |
CN103479510B (zh) * | 2013-09-26 | 2015-10-21 | 深圳先进技术研究院 | 一种穴位定位方法和系统 |
CN105433915B (zh) * | 2014-08-12 | 2018-01-26 | 上海华博信息服务有限公司 | 一种中医脉诊仪及其脉象检测方法 |
-
2017
- 2017-09-25 CN CN201710872547.4A patent/CN109419497B/zh active Active
- 2017-09-26 WO PCT/CN2017/103401 patent/WO2019041412A1/zh active Application Filing
Patent Citations (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20080045847A1 (en) * | 2006-06-30 | 2008-02-21 | University Of Louisville Research Foundation, Inc. | Non-contact and passive measurement of arterial pulse through thermal IR imaging, and analysis of thermal IR imagery |
CN101551851A (zh) * | 2008-03-31 | 2009-10-07 | 中国科学院沈阳自动化研究所 | 一种红外图像目标识别方法 |
CN102446349A (zh) * | 2010-09-29 | 2012-05-09 | 奥林巴斯株式会社 | 图像处理装置和图像处理方法 |
US20150133791A1 (en) * | 2012-05-29 | 2015-05-14 | National University Corporation Kochi University | Artery visualization device and artery imaging device |
CN103006195A (zh) * | 2013-01-10 | 2013-04-03 | 浙江大学 | 一种基于图像处理的非接触式生命体征数据监测系统和监测方法 |
WO2016018463A2 (en) * | 2013-12-23 | 2016-02-04 | Jason Case | Methods and systems using led sources to create contrast agents for mid-infrared imaging of blood vessels |
CN104268853A (zh) * | 2014-03-06 | 2015-01-07 | 上海大学 | 一种红外图像与可见光图像配准方法 |
CN104027094A (zh) * | 2014-05-30 | 2014-09-10 | 中国中医科学院医学实验中心 | 一种用于中医脉诊信息采集的寸、关、尺定位方法 |
CN104484667A (zh) * | 2014-12-30 | 2015-04-01 | 华中科技大学 | 一种基于亮度特征和轮廓完整性的轮廓提取方法 |
CN104799815A (zh) * | 2015-03-26 | 2015-07-29 | 山东大学 | 一种基于图像引导的血气采集装置及方法 |
CN106456056A (zh) * | 2015-05-13 | 2017-02-22 | 普勒斯麦公司 | 动脉可视化装置 |
CN105147261A (zh) * | 2015-08-03 | 2015-12-16 | 刘垚 | 一种中医脉诊仪及其定位寸、关、尺脉搏点的方法 |
CN105469414A (zh) * | 2015-12-15 | 2016-04-06 | 凌云光技术集团有限责任公司 | 轮廓连接方法及装置 |
CN106859956A (zh) * | 2017-01-13 | 2017-06-20 | 北京奇虎科技有限公司 | 一种人体穴位识别按摩方法、装置及ar设备 |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US11346938B2 (en) | 2019-03-15 | 2022-05-31 | Msa Technology, Llc | Safety device for providing output to an individual associated with a hazardous environment |
CN111297344A (zh) * | 2020-03-24 | 2020-06-19 | 陈恬慧 | 一种脉诊视觉阵列采集方法 |
CN111803071A (zh) * | 2020-06-23 | 2020-10-23 | 山东大学 | 耳朵外观轮廓三维扫描系统及方法 |
CN111803071B (zh) * | 2020-06-23 | 2021-10-15 | 山东大学 | 耳朵外观轮廓三维扫描系统及方法 |
CN113155288A (zh) * | 2020-11-30 | 2021-07-23 | 齐鲁工业大学 | 一种光伏电池热斑的图像识别方法 |
CN113155288B (zh) * | 2020-11-30 | 2022-09-06 | 齐鲁工业大学 | 一种光伏电池热斑的图像识别方法 |
CN113171062A (zh) * | 2021-04-29 | 2021-07-27 | 中国科学院微电子研究所 | 识别寸、关、尺的方法、终端、系统、介质及计算机设备 |
CN113171062B (zh) * | 2021-04-29 | 2024-03-26 | 中国科学院微电子研究所 | 识别寸、关、尺的方法、终端、系统、介质及计算机设备 |
CN113303771A (zh) * | 2021-07-30 | 2021-08-27 | 天津慧医谷科技有限公司 | 脉搏采集点的确定方法、装置和电子设备 |
CN113303771B (zh) * | 2021-07-30 | 2021-11-05 | 天津慧医谷科技有限公司 | 脉搏采集点的确定方法、装置和电子设备 |
Also Published As
Publication number | Publication date |
---|---|
WO2019041412A1 (zh) | 2019-03-07 |
CN109419497B (zh) | 2021-12-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109419497A (zh) | 基于热成像的关脉识别方法 | |
CN109427065A (zh) | 基于热成像的关脉识别系统 | |
CN105188522B (zh) | 脉搏波传播速度的测定方法、测定系统及摄像装置 | |
Sikdar et al. | Computer-vision-guided human pulse rate estimation: a review | |
Clarke et al. | Using high frame rate CMOS sensors for three-dimensional eye tracking | |
CN112465772B (zh) | 眼底彩照图像血管评估方法、装置、计算机设备和介质 | |
Hosp et al. | RemoteEye: An open-source high-speed remote eye tracker: Implementation insights of a pupil-and glint-detection algorithm for high-speed remote eye tracking | |
CN113171062B (zh) | 识别寸、关、尺的方法、终端、系统、介质及计算机设备 | |
CN102670168A (zh) | 眼科设备及其控制方法 | |
TW202038133A (zh) | 利用深度學習快速定位虹膜之裝置與方法 | |
CA2753818A1 (en) | A method for reconstructing iris scans through novel inpainting techniques and mosaicing of partial collections | |
CN102542242B (zh) | 非接触式采集图像的生物特征区域定位方法和装置 | |
CN103559489A (zh) | 一种非接触成像方式下手掌特征提取方法 | |
CN108292433A (zh) | 血管内图像中的阴影的检测与验证 | |
CN110276271A (zh) | 融合ippg和深度信息抗噪声干扰的非接触心率估计法 | |
TW200947262A (en) | Non-contact type cursor control method using human eye, pupil tracking system and storage media | |
CN112164043A (zh) | 一种用于多眼底图像拼接的方法及系统 | |
JPWO2010041731A1 (ja) | 照合装置、照合方法、及びプログラム | |
JP3720839B2 (ja) | 超音波診断装置 | |
Fukuda et al. | Model-based eye-tracking method for low-resolution eye-images | |
WO2018181363A1 (ja) | 診断画像処理装置、評価支援方法及びプログラム | |
JP6452236B2 (ja) | 眼球識別装置及び眼球識別方法 | |
Chen et al. | iVoiding: A thermal-image based artificial intelligence dynamic voiding detection system | |
Rafiqi et al. | Work-in-progress, PupilWare-M: Cognitive load estimation using unmodified smartphone cameras | |
Cao et al. | Gaze tracking on any surface with your phone |
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 |