CN107095655A - 基于三维数字图像相关的脉诊多维信息数字化系统及方法 - Google Patents
基于三维数字图像相关的脉诊多维信息数字化系统及方法 Download PDFInfo
- Publication number
- CN107095655A CN107095655A CN201710344008.3A CN201710344008A CN107095655A CN 107095655 A CN107095655 A CN 107095655A CN 201710344008 A CN201710344008 A CN 201710344008A CN 107095655 A CN107095655 A CN 107095655A
- Authority
- CN
- China
- Prior art keywords
- pulse
- measured
- value
- image
- acupuncture point
- 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/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/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/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/103—Detecting, measuring or recording devices for testing the shape, pattern, colour, size or movement of the body or parts thereof, for diagnostic purposes
- A61B5/11—Measuring movement of the entire body or parts thereof, e.g. head or hand tremor, mobility of a limb
- A61B5/1126—Measuring movement of the entire body or parts thereof, e.g. head or hand tremor, mobility of a limb using a particular sensing technique
- A61B5/1128—Measuring movement of the entire body or parts thereof, e.g. head or hand tremor, mobility of a limb using a particular sensing technique using image 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0012—Biomedical image inspection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30101—Blood vessel; Artery; Vein; Vascular
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- General Health & Medical Sciences (AREA)
- Physics & Mathematics (AREA)
- Molecular Biology (AREA)
- Pathology (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Biophysics (AREA)
- Physiology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Artificial Intelligence (AREA)
- Signal Processing (AREA)
- Psychiatry (AREA)
- Cardiology (AREA)
- Dentistry (AREA)
- Oral & Maxillofacial Surgery (AREA)
- Alternative & Traditional Medicine (AREA)
- Quality & Reliability (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Finger-Pressure Massage (AREA)
- Measuring And Recording Apparatus For Diagnosis (AREA)
Abstract
本发明提供一种基于三维数字图像相关的脉诊多维信息数字化系统及方法,该系统中支撑架与通过气泵和水银差分气压计充入目标气压值气体的透明气袋对待测体上的预制有散斑的待测穴位进行挤压,双目摄像机采集待测穴位的图像组并发送至上位机,上位机接收图像组,并基于三维数字图像相关法计算待测穴位区域内的全场三维物点在目标气压值下接收时间处因脉搏搏动引起的振动位移值,并且对设定时段内计算得到的全部全场三维物点的振动位移值按照相应目标气压值以及相应接收时间进行分析,得到待测穴位的脉诊多维信息。这就实现了脉诊的多维信息数字化,为中医脉诊数字化提供便利客观的诊疗手段,使医学大数据采集生成可能。
Description
技术领域
本发明涉及图像处理及中医领域,更具体地说,涉及一种基于三维数字图像相关的脉诊多维信息数字化系统及方法。
背景技术
脉诊是中医四诊的重要组成部分,是我过传统医院最早创用,并且至今一直沿用的最具中医特色的诊断方法。中医脉诊通常的做法就是中医师用手指触压患者寸口桡动脉处寸、关、尺三个部位,通过手指施加浮、中、沉等不同的压力来感受患者脉搏的变化。而对于浮、中、沉等压力产生的一系列包含了脉搏的宽度、长度以及振幅等信息的脉搏波,我们称之为脉象。从脉象信息中,中医师能够了解病人的疾病信息。但是这种传统的中医脉诊有很大的主观性,准确掌握和运用有一定难度,因此不能保证脉诊结果真实准确。
为了解决传统中医脉诊的缺点,脉象传感器应运而生,而大部分的脉象传感器是普通压力传感器,由于普通压力传感器通常只能测量一个固定压力下的脉搏信号,这就使得脉象传感器不能针对浮、中、沉等压力采集脉搏信号,从而导致所获取的脉象信息的精度很低。
发明内容
有鉴于此,本发明提供一种基于三维数字图像相关的脉诊多维信息数字化系统及方法,以解决脉象传感器不能针对浮、中、沉等压力采集脉搏信号,从而导致所获取的脉象信息的精度很低的问题。技术方案如下:
一种基于三维数字图像相关的脉诊多维信息数字化系统,包括:上位机、与所述上位机相连的双目摄像机、气泵、支撑架、与所述气泵一端相连的水银差分气压计和与所述气泵另一端相连的透明气袋;
其中,所述支撑架用于支撑被测体、为充气过程中的所述透明气袋提供位移约束,并与充入目标气压值气体的所述透明气袋对所述被测体上的待测穴位进行挤压,所述待测穴位区域内皮肤表面上预制有散斑;
所述双目摄像机,用于当接收到所述上位机发送的图像采集通知时,控制左相机和右相机同时采集所述被测穴位的图像,并实时将包含有所述左相机采集的左图像以及所述右相机采集的右图像的图像组发送至所述上位机;
所述上位机,用于当接收到所述图像组时,记录接收时间,并根据三维数字图像相关法和该图像组计算所述待测穴位区域内的全场三维物点在所述目标气压值下所述接收时间处因脉搏搏动引起的振动位移值,以及对设定时段内计算得到的全部所述全场三维物点的所述振动位移值按照相应目标气压值以及相应接收时间进行分析,得到所述待测穴位在所述设定时段内的脉诊多维信息。
优选的,所述透明气袋由三个平行排列的透明柱体构成,并且各个所述透明柱体均与所述气泵相连。
优选的,用于根据三维数字图像相关法和该图像组计算所述待测穴位区域内的全场三维物点在所述目标气压值下所述接收时间处因脉搏搏动引起的振动位移值的所述上位机,具体用于:
从该图像组中包含的所述左图像和所述右图像中选取一个作为基准图像,另一个则作为配准图像,以及对于所述待测穴位区域内的各个全场三维物点,从所述基准图像上选取相应像素点,并从所述配准图像上查找与所述像素点相匹配的匹配点,以及计算所述像素点与所述匹配点在相应图像坐标下的像素位移值,并将所述像素位移值确定为所述像素点的视差值,以及根据所述视差值计算所述像素点在三维空间内的深度值,以及依据所述深度值计算所述全场三维物点在所述目标气压值下所述接收时间处的振动位移值。
优选的,用于对设定时段内计算得到的全部所述全场三维物点的所述振动位移值按照相应目标气压值以及相应接收时间进行分析,得到所述待测穴位在所述设定时段内的脉诊多维信息的所述上位机,具体用于:
根据设定时段内计算得到的各个所述全场三维物点的所述振动位移值,生成所述待测穴位在相应目标气压值下相应接收时间处的全场位移分布图,以及对于目标气压值相同的全部所述全场位移分布图,从中选取振动位移值最大的所述全场位移分布图,并根据选取的所述全场位移分布图计算所述待测穴位在相应目标气压值下的脉搏宽度和脉搏长度,以及对于目标气压值相同的全部所述全场位移分布图,计算所述待测穴位在相应目标气压值下的振动位移均值,并根据相应接收时间生成所述待测穴位在相应目标气压值下的用于表征时间与振动位移关系的第一脉搏振幅曲线,以及当生成至少两组不同目标气压值下的所述第一脉搏振幅曲线时,从各个所述第一脉搏振幅曲线上获取相应最大振幅值,并根据相应目标气压值生成所述待测穴位的用于表征气压与最大振幅关系的第二脉搏振幅曲线,以及当生成至少两组相同目标气压值下的所述第一脉搏振幅曲线时,从全部所述第一脉搏振幅曲线上获取各个接收时间的最大振幅值,并生成所述待测穴位的用于表征时间与最大振幅关系的第三脉搏振幅曲线和/或波峰位移分布图。
一种基于三维数字图像相关的脉诊多维信息数字化方法,应用于上述技术方案任意一项所述的系统,该系统包括上位机、双目摄像机、气泵、支撑架、水银差分气压计和透明气袋,所述方法包括:
当所述双目摄像机接收到所述上位机发送的图像采集通知时,控制左相机和右相机同时采集所述被测穴位的图像,并实时将包含有所述左相机采集的左图像以及所述右相机采集的右图像的图像组发送至所述上位机;
当所述上位机接收到所述图像组时,记录接收时间,并根据三维数字图像相关法和该图像组计算所述待测穴位区域内的全场三维物点在所述目标气压值下所述接收时间处因脉搏搏动引起的振动位移值,以及对设定时段内计算得到的全部所述全场三维物点的所述振动位移值按照相应目标气压值以及相应接收时间进行分析,得到所述待测穴位在所述设定时段内的脉诊多维信息。
优选的,所述上位机根据三维数字图像相关法和该图像组计算所述待测穴位区域内的全场三维物点在所述目标气压值下所述接收时间处因脉搏搏动引起的振动位移值,包括:
从该图像组中包含的所述左图像和所述右图像中选取一个作为基准图像,另一个则作为配准图像;
对于所述待测穴位区域内的各个全场三维物点,从所述基准图像上选取相应像素点,并从所述配准图像上查找与所述像素点相匹配的匹配点;
计算所述像素点与所述匹配点在相应图像坐标下的像素位移值,并将所述像素位移值确定为所述像素点的视差值;
根据所述视差值计算所述像素点在三维空间内的深度值;
依据所述深度值计算所述全场三维物点在所述目标气压值下所述接收时间处的振动位移值。
优选的,所述上位机对设定时段内计算得到的全部所述全场三维物点的所述振动位移值按照相应目标气压值以及相应接收时间进行分析,得到所述待测穴位在所述设定时段内的脉诊多维信息,包括:
根据设定时段内计算得到的各个所述全场三维物点的所述振动位移值,生成所述待测穴位在相应目标气压值下相应接收时间处的全场位移分布图;
对于目标气压值相同的全部所述全场位移分布图,从中选取振动位移值最大的所述全场位移分布图,并根据选取的所述全场位移分布图计算所述待测穴位在相应目标气压值下的脉搏宽度和脉搏长度;
对于目标气压值相同的全部所述全场位移分布图,计算所述待测穴位在相应目标气压值下的振动位移均值,并根据相应接收时间生成所述待测穴位在相应目标气压值下的用于表征时间与振动位移关系的第一脉搏振幅曲线;
当生成至少两组不同目标气压值下的所述第一脉搏振幅曲线时,从各个所述第一脉搏振幅曲线上获取相应最大振幅值,并根据相应目标气压值生成所述待测穴位的用于表征气压与最大振幅关系的第二脉搏振幅曲线;
当生成至少两组相同目标气压值下的所述第一脉搏振幅曲线时,从全部所述第一脉搏振幅曲线上获取各个接收时间的最大振幅值,并生成所述待测穴位的用于表征时间与最大振幅关系的第三脉搏振幅曲线和/或波峰位移分布图。
相较于现有技术,本发明实现的有益效果为:
本发明公开一种基于三维数字图像相关的脉诊多维信息数字化系统及方法,该系统中支撑架与通过气泵和水银差分气压计充入目标气压值气体的透明气袋对待测体上的预制有散斑的待测穴位进行挤压,双目摄像机采集待测穴位的图像组并发送至上位机,上位机接收图像组,并基于三维数字图像相关法计算待测穴位区域内的全场三维物点在目标气压值下接收时间处因脉搏搏动引起的振动位移值,并且对设定时段内计算得到的全部全场三维物点的振动位移值按照相应目标气压值以及相应接收时间进行分析,得到待测穴位的脉诊多维信息。这就实现了脉诊的多维信息数字化,为中医脉诊数字化提供便利客观的诊疗手段,使医学大数据采集生成可能。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据提供的附图获得其他的附图。
图1为本发明实施例提供的基于三维数字图像相关的脉诊多维信息数字化系统的结构示意图;
图2为本发明实施例提供的全场位移分布图;
图3为本发明实施例提供的又一全场位移分布图;
图4为本发明实施例提供的第一脉搏振幅曲线;
图5为本发明实施例提供的第二脉搏振幅曲线;
图6为本发明实施例提供的第三脉搏振幅曲线以及波峰位移分布图;
图7为本发明实施例提供的基于三维数字图像相关的脉诊多维信息数字化方法的方法流程图;
图8为本发明实施例提供的基于三维数字图像相关的脉诊多维信息数字化方法的部分方法流程图;
图9为本发明实施例提供的基于三维数字图像相关的脉诊多维信息数字化方法的又一部分方法流程图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明实施例公开一种基于三维数字图像相关的脉诊多维信息数字化系统,包括:上位机101、与上位机101相连的双目摄像机102、气泵103、支撑架104、与气泵103一端相连的水银差分气压计105和与气泵103另一端相连的透明气袋106;
其中,支撑架104用于支撑被测体、为充气过程中的透明气袋106提供位移约束,并与充入目标气压值气体的透明气袋106对被测体上的待测穴位进行挤压,待测穴位区域内皮肤表面上预制有散斑;
双目摄像机102,用于当接收到上位机101发送的图像采集通知时,控制左相机和右相机同时采集被测穴位的图像,并实时将包含有左相机采集的左图像以及右相机采集的右图像的图像组发送至上位机101;
上位机101,用于当接收到图像组时,记录接收时间,并根据三维数字图像相关法和该图像组计算待测穴位区域内的全场三维物点在目标气压值下接收时间处因脉搏搏动引起的振动位移值,以及对设定时段内计算得到的全部全场三维物点的振动位移值按照相应目标气压值以及相应接收时间进行分析,得到待测穴位在设定时段内的脉诊多维信息。
在开启上述脉诊多维信息数字化系统之前,为使待测穴位表面具有足够的、不重复的并且各项同性的纹理,可在待测穴位区域内的表面皮肤上制备散斑,本实施例中可使用水转印散斑贴于皮肤表面,根据待测穴位在支撑架中的位置选取双目摄像机的左相机和右相机,并对双目摄像机进行标定,也就是根据摄像机参数调整左相机以及右相机的位置,以将待测穴位成像于左相机与右相机的画幅之中,进一步的,为保证上位机生成多维脉象信息的高效性,可将上位机的脉象信息生成方式设置为高帧率的后处理生成方式,也可设置为低帧率的实时生成方式,优选的,高帧率可为472帧率,低帧率可为10帧率。
可选的,上位机101包括但不局限于手机、电脑、平板或者笔记本,用户可根据实际需要具体设置,本实施例不做任何限定。
可选的,气泵103包括气囊,用户可根据实际需要具体设置,本实施例不做任何限定。
可选的,支撑架104可为透明盒子,用户可根据实际需要具体设置,本实施例不做任何限定。
可选的,透明气袋106可为透明的3微米厚的5cm*7cm是气袋,气袋的规格参数本实施例不做任何限定;
另外,为真实模拟三根手指对待测穴位的按压,优选的,透明气袋106可由三个平行排列的透明柱体构成,并且各个透明柱体均与气泵103相连,具体的,该透明柱体可为半径为0.5cm、长度为5cm的圆柱体。
图1为被测体为患者寸口桡动脉、上位机101为电脑、气泵103为气囊,并且支撑架104可为透明盒子时,采用本实施例公开的脉诊多维信息数字化系统对患者寸口桡动脉处进行脉象采集的示意图。
在脉象采集过程中,医生可通过观看水银差分气压计105的读数,使用气囊103将透明气袋106的气压值调整目标气压值,并与透明盒子104完成对待测穴位的按压,进一步,为保证患者的舒适性,如图1所示,还可在透明盒子104底部设置一个脉枕。
可选的,用于根据三维数字图像相关法和该图像组计算待测穴位区域内的全场三维物点在目标气压值下接收时间处因脉搏搏动引起的振动位移值的上位机101,具体用于:
从该图像组中包含的左图像和右图像中选取一个作为基准图像,另一个则作为配准图像,以及对于待测穴位区域内的各个全场三维物点,从基准图像上选取相应像素点,并从配准图像上查找与像素点相匹配的匹配点,以及计算像素点与匹配点在相应图像坐标下的像素位移值,并将像素位移值确定为像素点的视差值,以及根据视差值计算像素点在三维空间内的深度值,以及依据深度值计算全场三维物点在目标气压值下接收时间处的振动位移值。
本实施例中,对于待测穴位区域内的各个全场三维物点,在左相机采集的左图像以及右相机采集的右图像上均存在该全场三维物点的像素点,这种由同一全场三维物点在不同成像平面上所成的像素点可构成一对匹配点;
在上位机接收到第一帧图像组时,可将该帧图像组设置为参考帧,按照预设规则从图像组中选取基准图像以及配准图像,选取方式可根据时候需要具体选送,例如,选取左相机采集的左图像作为基准图像,选取右相机采集的右图像作为配准图像,对于每一个全场三维物点,可从基准图像上选取相应像素点,然后从配准图像上查找与该像素点匹配的匹配点,并根据像素点及其匹配点在各自图像坐标下的像素坐标计算像素坐标位移值,并将该像素坐标位移值确定为视差值,进而根据视差值以及双目摄像机的标定参数计算像素点在三维空间内的深度值;最后,在接收到下一帧图像组时,以第一帧图像组的深度值作为参考值,计算下一帧图像组的深度差值,并将计算得到的深度差值作为此全场三维物点在目标气压值下相应接收时间处的振动位移值。
可选的,用于对设定时段内计算得到的全部全场三维物点的振动位移值按照相应目标气压值以及相应接收时间进行分析,得到待测穴位在设定时段内的脉诊多维信息的上位机101,具体用于:
根据设定时段内计算得到的各个全场三维物点的振动位移值,生成待测穴位在相应目标气压值下相应接收时间处的全场位移分布图,以及对于目标气压值相同的全部全场位移分布图,从中选取振动位移值最大的全场位移分布图,并根据选取的全场位移分布图计算待测穴位在相应目标气压值下的脉搏宽度和脉搏长度,以及对于目标气压值相同的全部全场位移分布图,计算待测穴位在相应目标气压值下的振动位移均值,并根据相应接收时间生成待测穴位在相应目标气压值下的用于表征时间与振动位移关系的第一脉搏振幅曲线,以及当生成至少两组不同目标气压值下的第一脉搏振幅曲线时,从各个第一脉搏振幅曲线上获取相应最大振幅值,并根据相应目标气压值生成待测穴位的用于表征气压与最大振幅关系的第二脉搏振幅曲线,以及当生成至少两组相同目标气压值下的第一脉搏振幅曲线时,从全部第一脉搏振幅曲线上获取各个接收时间的最大振幅值,并生成待测穴位的用于表征时间与最大振幅关系的第三脉搏振幅曲线和/或波峰位移分布图。
本实施例中,全场位移分布图可采用不同颜色区别各个全场三维物点的振动位移值,例如,振动位移值越大,则颜色越深;
由于脉搏搏动引起各全场三维物点到达最大振动位移值的时间是一致的,因此,可从目标气压值相同的全部全场位移分布图中选取各全场三维物点振动位移值最大的全场位移分布图,选取的位移值最大的全场位移分布图可为一个或多个,对于选取的每一个全场位移分布图,可生成图中包含全部全场三维物点的最小矩形,并计算该最小矩形的长度和宽度,因此,如果选取的全场位移分布图为一个时,计算得到的宽度就为脉搏宽度,计算得到的长度就为脉搏长度,而如果选取的全场位移分布图为至少两个时,可计算最小矩形的长度平均值以及宽度平均值,并将长度平均值作为脉搏长度,将长度宽度值作为脉搏宽度。
需要说明的是,本实施公开的第一脉搏振幅曲线、第二脉搏振幅曲线、第三脉搏振幅曲线以及波峰位移分布图仅仅是根据全场位移分布图分析得到的几种多维脉象信息,用于还可根据生成的全场位移分布图控制其中一个或几个变量进行分析,得到其他的脉象多维信息,本实施例不做任何限定。
本发明实施例公开基于三维数字图像相关的脉诊多维信息数字化系统中,支撑架与通过气泵和水银差分气压计充入目标气压值气体的透明气袋对待测体上的预制有散斑的待测穴位进行挤压,双目摄像机采集待测穴位的图像组并发送至上位机,上位机接收图像组,并基于三维数字图像相关法计算待测穴位区域内的全场三维物点在目标气压值下接收时间处因脉搏搏动引起的振动位移值,并且对设定时段内计算得到的全部全场三维物点的振动位移值按照相应目标气压值以及相应接收时间进行分析,得到待测穴位的脉诊多维信息。这就实现了脉诊的多维信息数字化,为中医脉诊数字化提供便利客观的诊疗手段,使医学大数据采集生成可能。
利用上述图1提供的基于三维数字图像相关的脉诊多维信息数字化系统对患者寸口桡动脉处寸、关、尺三个穴位进行挤压;
如图2所示,为脉诊多维信息数字化系统以60mmHg目标气压值对寸、关、尺三个穴位进行挤压,生成的Tc时刻下的全场位移分布图;
如图3所示,为脉诊多维信息数字化系统以60mmHg目标气压值对寸、关、尺三个穴位进行挤压,选取的位移值最大的Td时刻下全场位移分布图,其中,Td时刻下全场位移分布图中最小矩形的长度值为6cm、宽度值为1.2cm;
如图4所示,为脉诊多维信息数字化系统以60mmHg目标气压值对寸、关、尺三个穴位进行挤压,寸、关、尺三个穴位在60mmHg目标气压值下的第一脉搏振幅曲线,其中,寸穴位用cun表示、关穴位用guan表示、尺穴位为用chi表示;
如图5所示,寸、关、尺三个穴位在目标气压值下的第二脉搏振幅曲线,其中,寸穴位用cun表示、关穴位用guan表示、尺穴位为用chi表示;
如图6所示,为脉诊多维信息数字化系统以30mmHg目标气压值对寸、关、尺三个穴位进行挤压,寸、关、尺三个穴位在30mmHg目标气压值下的第三脉搏振幅曲线以及波峰位移分布图,其中,寸穴位用cun表示、关穴位用guan表示、尺穴位为用chi表示。
基于上述实施例提供的基于三维数字图像相关的脉诊多维信息数字化系统,本发明又一实施例公开一种基于三维数字图像相关的脉诊多维信息数字化方法,方法流程图如图7所示,包括如下步骤:
S10,当双目摄像机接收到上位机发送的图像采集通知时,控制左相机和右相机同时采集被测穴位的图像,并实时将包含有左相机采集的左图像以及右相机采集的右图像的图像组发送至上位机;
S20,当上位机接收到图像组时,记录接收时间,并根据三维数字图像相关法和该图像组计算待测穴位区域内的全场三维物点在目标气压值下接收时间处因脉搏搏动引起的振动位移值,以及对设定时段内计算得到的全部全场三维物点的振动位移值按照相应目标气压值以及相应接收时间进行分析,得到待测穴位在设定时段内的脉诊多维信息。
在具体实现过程中,步骤S20中“上位机根据三维数字图像相关法和该图像组计算待测穴位区域内的全场三维物点在目标气压值下接收时间处因脉搏搏动引起的振动位移值”可以具体采用以下步骤,方法流程图如图8所示:
S201,从该图像组中包含的左图像和右图像中选取一个作为基准图像,另一个则作为配准图像;
S202,对于待测穴位区域内的各个全场三维物点,从基准图像上选取相应像素点,并从配准图像上查找与像素点相匹配的匹配点;
S203,计算像素点与匹配点在相应图像坐标下的像素位移值,并将像素位移值确定为像素点的视差值;
S204,根据视差值计算像素点在三维空间内的深度值;
S205,依据深度值计算全场三维物点在目标气压值下接收时间处的振动位移值。
在具体实现过程中,步骤S20中“上位机对设定时段内计算得到的全部全场三维物点的振动位移值按照相应目标气压值以及相应接收时间进行分析,得到待测穴位在设定时段内的脉诊多维信息”可以具体采用以下步骤,方法流程图如图9所示:
S206,根据设定时段内计算得到的各个全场三维物点的振动位移值,生成待测穴位在相应目标气压值下相应接收时间处的全场位移分布图;
S207,对于目标气压值相同的全部全场位移分布图,从中选取振动位移值最大的全场位移分布图,并根据选取的全场位移分布图计算待测穴位在相应目标气压值下的脉搏宽度和脉搏长度;
S208,对于目标气压值相同的全部全场位移分布图,计算待测穴位在相应目标气压值下的振动位移均值,并根据相应接收时间生成待测穴位在相应目标气压值下的用于表征时间与振动位移关系的第一脉搏振幅曲线;
S209,当生成至少两组不同目标气压值下的第一脉搏振幅曲线时,从各个第一脉搏振幅曲线上获取相应最大振幅值,并根据相应目标气压值生成待测穴位的用于表征气压与最大振幅关系的第二脉搏振幅曲线;
S210,当生成至少两组相同目标气压值下的第一脉搏振幅曲线时,从全部第一脉搏振幅曲线上获取各个接收时间的最大振幅值,并生成待测穴位的用于表征时间与最大振幅关系的第三脉搏振幅曲线和/或波峰位移分布图。
以上步骤S201~步骤S205仅仅是本申请实施例公开的步骤S20中“上位机根据三维数字图像相关法和该图像组计算待测穴位区域内的全场三维物点在目标气压值下接收时间处因脉搏搏动引起的振动位移值”过程的一种优选的实现方式,有关此过程的具体实现方式可根据自己的需求任意设置,在此不做限定。
以上步骤S206~步骤S210仅仅是本申请实施例公开的步骤S20中“上位机对设定时段内计算得到的全部全场三维物点的振动位移值按照相应目标气压值以及相应接收时间进行分析,得到待测穴位在设定时段内的脉诊多维信息”过程的一种优选的实现方式,有关此过程的具体实现方式可根据自己的需求任意设置,在此不做限定。
本发明实施例公开的基于三维数字图像相关的脉诊多维信息数字化方法,包括:双目摄像机采集待测穴位的图像组并发送至上位机,上位机接收图像组,并基于三维数字图像相关法计算待测穴位区域内的全场三维物点在目标气压值下接收时间处因脉搏搏动引起的振动位移值,并且对设定时段内计算得到的全部全场三维物点的振动位移值按照相应目标气压值以及相应接收时间进行分析,得到待测穴位的脉诊多维信息。这就实现了脉诊的多维信息数字化,为中医脉诊数字化提供便利客观的诊疗手段,使医学大数据采集生成可能。
以上对本发明所提供的一种基于三维数字图像相关的脉诊多维信息数字化系统及方法进行了详细介绍,本文中应用了具体个例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处,综上所述,本说明书内容不应理解为对本发明的限制。
需要说明的是,本说明书中的各个实施例均采用递进的方式描述,每个实施例重点说明的都是与其他实施例的不同之处,各个实施例之间相同相似的部分互相参见即可。对于实施例公开的系统而言,由于其与实施例公开的方法相对应,所以描述的比较简单,相关之处参见方法部分说明即可。
还需要说明的是,在本文中,诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来,而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、物品或者设备所固有的要素,或者是还包括为这些过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下,由语句“包括一个……”限定的要素,并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。
对所公开的实施例的上述说明,使本领域专业技术人员能够实现或使用本发明。对这些实施例的多种修改对本领域的专业技术人员来说将是显而易见的,本文中所定义的一般原理可以在不脱离本发明的精神或范围的情况下,在其它实施例中实现。因此,本发明将不会被限制于本文所示的这些实施例,而是要符合与本文所公开的原理和新颖特点相一致的最宽的范围。
Claims (7)
1.一种基于三维数字图像相关的脉诊多维信息数字化系统,其特征在于,包括:上位机、与所述上位机相连的双目摄像机、气泵、支撑架、与所述气泵一端相连的水银差分气压计和与所述气泵另一端相连的透明气袋;
其中,所述支撑架用于支撑被测体、为充气过程中的所述透明气袋提供位移约束,并与充入目标气压值气体的所述透明气袋对所述被测体上的待测穴位进行挤压,所述待测穴位区域内皮肤表面上预制有散斑;
所述双目摄像机,用于当接收到所述上位机发送的图像采集通知时,控制左相机和右相机同时采集所述被测穴位的图像,并实时将包含有所述左相机采集的左图像以及所述右相机采集的右图像的图像组发送至所述上位机;
所述上位机,用于当接收到所述图像组时,记录接收时间,并根据三维数字图像相关法和该图像组计算所述待测穴位区域内的全场三维物点在所述目标气压值下所述接收时间处因脉搏搏动引起的振动位移值,以及对设定时段内计算得到的全部所述全场三维物点的所述振动位移值按照相应目标气压值以及相应接收时间进行分析,得到所述待测穴位在所述设定时段内的脉诊多维信息。
2.根据权利要求1所述的系统,其特征在于,所述透明气袋由三个平行排列的透明柱体构成,并且各个所述透明柱体均与所述气泵相连。
3.根据权利要求1所述的系统,其特征在于,用于根据三维数字图像相关法和该图像组计算所述待测穴位区域内的全场三维物点在所述目标气压值下所述接收时间处因脉搏搏动引起的振动位移值的所述上位机,具体用于:
从该图像组中包含的所述左图像和所述右图像中选取一个作为基准图像,另一个则作为配准图像,以及对于所述待测穴位区域内的各个全场三维物点,从所述基准图像上选取相应像素点,并从所述配准图像上查找与所述像素点相匹配的匹配点,以及计算所述像素点与所述匹配点在相应图像坐标下的像素位移值,并将所述像素位移值确定为所述像素点的视差值,以及根据所述视差值计算所述像素点在三维空间内的深度值,以及依据所述深度值计算所述全场三维物点在所述目标气压值下所述接收时间处的振动位移值。
4.根据权利要求1所述的系统,其特征在于,用于对设定时段内计算得到的全部所述全场三维物点的所述振动位移值按照相应目标气压值以及相应接收时间进行分析,得到所述待测穴位在所述设定时段内的脉诊多维信息的所述上位机,具体用于:
根据设定时段内计算得到的各个所述全场三维物点的所述振动位移值,生成所述待测穴位在相应目标气压值下相应接收时间处的全场位移分布图,以及对于目标气压值相同的全部所述全场位移分布图,从中选取振动位移值最大的所述全场位移分布图,并根据选取的所述全场位移分布图计算所述待测穴位在相应目标气压值下的脉搏宽度和脉搏长度,以及对于目标气压值相同的全部所述全场位移分布图,计算所述待测穴位在相应目标气压值下的振动位移均值,并根据相应接收时间生成所述待测穴位在相应目标气压值下的用于表征时间与振动位移关系的第一脉搏振幅曲线,以及当生成至少两组不同目标气压值下的所述第一脉搏振幅曲线时,从各个所述第一脉搏振幅曲线上获取相应最大振幅值,并根据相应目标气压值生成所述待测穴位的用于表征气压与最大振幅关系的第二脉搏振幅曲线,以及当生成至少两组相同目标气压值下的所述第一脉搏振幅曲线时,从全部所述第一脉搏振幅曲线上获取各个接收时间的最大振幅值,并生成所述待测穴位的用于表征时间与最大振幅关系的第三脉搏振幅曲线和/或波峰位移分布图。
5.一种基于三维数字图像相关的脉诊多维信息数字化方法,其特征在于,应用于权利要求1~4任意一项所述的系统,该系统包括上位机、双目摄像机、气泵、支撑架、水银差分气压计和透明气袋,所述方法包括:
当所述双目摄像机接收到所述上位机发送的图像采集通知时,控制左相机和右相机同时采集所述被测穴位的图像,并实时将包含有所述左相机采集的左图像以及所述右相机采集的右图像的图像组发送至所述上位机;
当所述上位机接收到所述图像组时,记录接收时间,并根据三维数字图像相关法和该图像组计算所述待测穴位区域内的全场三维物点在所述目标气压值下所述接收时间处因脉搏搏动引起的振动位移值,以及对设定时段内计算得到的全部所述全场三维物点的所述振动位移值按照相应目标气压值以及相应接收时间进行分析,得到所述待测穴位在所述设定时段内的脉诊多维信息。
6.根据权利要求5所述的方法,其特征在于,所述上位机根据三维数字图像相关法和该图像组计算所述待测穴位区域内的全场三维物点在所述目标气压值下所述接收时间处因脉搏搏动引起的振动位移值,包括:
从该图像组中包含的所述左图像和所述右图像中选取一个作为基准图像,另一个则作为配准图像;
对于所述待测穴位区域内的各个全场三维物点,从所述基准图像上选取相应像素点,并从所述配准图像上查找与所述像素点相匹配的匹配点;
计算所述像素点与所述匹配点在相应图像坐标下的像素位移值,并将所述像素位移值确定为所述像素点的视差值;
根据所述视差值计算所述像素点在三维空间内的深度值;
依据所述深度值计算所述全场三维物点在所述目标气压值下所述接收时间处的振动位移值。
7.根据权利要求5所述的方法,其特征在于,所述上位机对设定时段内计算得到的全部所述全场三维物点的所述振动位移值按照相应目标气压值以及相应接收时间进行分析,得到所述待测穴位在所述设定时段内的脉诊多维信息,包括:
根据设定时段内计算得到的各个所述全场三维物点的所述振动位移值,生成所述待测穴位在相应目标气压值下相应接收时间处的全场位移分布图;
对于目标气压值相同的全部所述全场位移分布图,从中选取振动位移值最大的所述全场位移分布图,并根据选取的所述全场位移分布图计算所述待测穴位在相应目标气压值下的脉搏宽度和脉搏长度;
对于目标气压值相同的全部所述全场位移分布图,计算所述待测穴位在相应目标气压值下的振动位移均值,并根据相应接收时间生成所述待测穴位在相应目标气压值下的用于表征时间与振动位移关系的第一脉搏振幅曲线;
当生成至少两组不同目标气压值下的所述第一脉搏振幅曲线时,从各个所述第一脉搏振幅曲线上获取相应最大振幅值,并根据相应目标气压值生成所述待测穴位的用于表征气压与最大振幅关系的第二脉搏振幅曲线;
当生成至少两组相同目标气压值下的所述第一脉搏振幅曲线时,从全部所述第一脉搏振幅曲线上获取各个接收时间的最大振幅值,并生成所述待测穴位的用于表征时间与最大振幅关系的第三脉搏振幅曲线和/或波峰位移分布图。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710344008.3A CN107095655B (zh) | 2017-05-16 | 2017-05-16 | 基于三维数字图像相关的脉诊多维信息数字化系统及方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710344008.3A CN107095655B (zh) | 2017-05-16 | 2017-05-16 | 基于三维数字图像相关的脉诊多维信息数字化系统及方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107095655A true CN107095655A (zh) | 2017-08-29 |
CN107095655B CN107095655B (zh) | 2020-05-12 |
Family
ID=59669750
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710344008.3A Active CN107095655B (zh) | 2017-05-16 | 2017-05-16 | 基于三维数字图像相关的脉诊多维信息数字化系统及方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107095655B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107898626A (zh) * | 2017-12-20 | 2018-04-13 | 大连交通大学 | 一种双目视觉穴位定位仪及其工作方法 |
CN108784654A (zh) * | 2018-03-28 | 2018-11-13 | 广州中康科技发展有限公司 | 基于激光三维扫描的脉象五维信息自动获取的设备和方法 |
CN109124595A (zh) * | 2018-07-20 | 2019-01-04 | 南开大学 | 基于动态三维脉搏波图像的智能脉搏诊断方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5365930A (en) * | 1991-03-18 | 1994-11-22 | Sony Corporation | Apparatus for measuring and evaluating blood circulation |
CN1491611A (zh) * | 2002-10-24 | 2004-04-28 | 成都三九数码医药设备有限公司 | 脉象信息的采集和处理方法 |
CN102151126A (zh) * | 2010-05-28 | 2011-08-17 | 兰州理工大学 | 脉搏动态图像信息的检测装置 |
CN106595528A (zh) * | 2016-11-10 | 2017-04-26 | 华中科技大学 | 一种基于数字散斑的远心显微双目立体视觉测量方法 |
-
2017
- 2017-05-16 CN CN201710344008.3A patent/CN107095655B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5365930A (en) * | 1991-03-18 | 1994-11-22 | Sony Corporation | Apparatus for measuring and evaluating blood circulation |
CN1491611A (zh) * | 2002-10-24 | 2004-04-28 | 成都三九数码医药设备有限公司 | 脉象信息的采集和处理方法 |
CN102151126A (zh) * | 2010-05-28 | 2011-08-17 | 兰州理工大学 | 脉搏动态图像信息的检测装置 |
CN106595528A (zh) * | 2016-11-10 | 2017-04-26 | 华中科技大学 | 一种基于数字散斑的远心显微双目立体视觉测量方法 |
Non-Patent Citations (1)
Title |
---|
SHAO X, DAI X, CHEN Z, ET AL: "Real-time 3D digital image correlation method and its application in human pulse monitoring", 《APPLIED OPTICS》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107898626A (zh) * | 2017-12-20 | 2018-04-13 | 大连交通大学 | 一种双目视觉穴位定位仪及其工作方法 |
CN108784654A (zh) * | 2018-03-28 | 2018-11-13 | 广州中康科技发展有限公司 | 基于激光三维扫描的脉象五维信息自动获取的设备和方法 |
CN108784654B (zh) * | 2018-03-28 | 2023-05-30 | 广州中康资讯股份有限公司 | 基于激光三维扫描的脉象五维信息自动获取的设备和方法 |
CN109124595A (zh) * | 2018-07-20 | 2019-01-04 | 南开大学 | 基于动态三维脉搏波图像的智能脉搏诊断方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107095655B (zh) | 2020-05-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Ghosh et al. | Three-dimensional reconstruction of echocardiographic images using the rotation method | |
CN106600596B (zh) | 一种心脏腔室动态显示方法及系统 | |
CN107427236A (zh) | 用于心血管监测的心血管传感器同步 | |
CN103126655B (zh) | 一种非约束目标的非接触式脉搏波采集系统及采集方法 | |
CN107095655A (zh) | 基于三维数字图像相关的脉诊多维信息数字化系统及方法 | |
CN110325141A (zh) | 影像整合装置及影像整合方法 | |
JP6442695B2 (ja) | 生体情報測定装置および生体情報測定装置もしくは生体情報測定システムの作動方法 | |
CN107092771A (zh) | 采用定期数据整合的个性化模型 | |
CN101129261A (zh) | 一种获取识别脉动信息与舌诊信息的装置及其方法 | |
US20150011894A1 (en) | System for and method of quantifying on-body palpitation for improved medical diagnosis | |
US20230330507A1 (en) | Muscle training method and system for providing visual feedback by using ultrasonic imaging | |
CN110279404A (zh) | 一种脉搏波采集装置及方法 | |
US20230218178A1 (en) | Construction method and application of digital human cardiovascular system based on hemodynamics | |
CN105718730B (zh) | 对受试者的疼痛量化评估方法以及实施该方法的系统 | |
Catherwood et al. | Validation of a passive stereophotogrammetry system for imaging of the breast: a geometric analysis | |
EP3188120B1 (en) | Registering first image data of a first stream with second image data of a second stream | |
CN102579017A (zh) | 无创血流动力学参数分析仪 | |
Campo et al. | Digital image correlation for full-field time-resolved assessment of arterial stiffness | |
Azhari et al. | An Analytical Descrptor of Three-Dimensional Geometry: Application to the Analysis of the Left Ventricle Shape and Contraction | |
CN114159091B (zh) | 一种基于可穿戴传感器阵列的心音传播关系检测系统 | |
CN102217952B (zh) | 基于心肌运动轨迹的向量环图生成方法和装置 | |
Stork | Simulation of ECG, blood pressure and ballistocardiographic signals | |
CN107038716A (zh) | 一种基于椭球模型的非刚性医学图像抗旋转配准方法 | |
CN110720896A (zh) | 脉象信号处理方法及装置 | |
Mahr | Three-dimensional, in-silico breast phantom for multimodality image simulations |
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 |