CN111696687B - 一种采用手机摄像头实现微循环及脉搏同步监测的方法 - Google Patents

一种采用手机摄像头实现微循环及脉搏同步监测的方法 Download PDF

Info

Publication number
CN111696687B
CN111696687B CN202010722872.4A CN202010722872A CN111696687B CN 111696687 B CN111696687 B CN 111696687B CN 202010722872 A CN202010722872 A CN 202010722872A CN 111696687 B CN111696687 B CN 111696687B
Authority
CN
China
Prior art keywords
formula
microcirculation
channel
image
speckle
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
Application number
CN202010722872.4A
Other languages
English (en)
Other versions
CN111696687A (zh
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.)
University of Shanghai for Science and Technology
Original Assignee
University of Shanghai for 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 University of Shanghai for Science and Technology filed Critical University of Shanghai for Science and Technology
Priority to CN202010722872.4A priority Critical patent/CN111696687B/zh
Publication of CN111696687A publication Critical patent/CN111696687A/zh
Application granted granted Critical
Publication of CN111696687B publication Critical patent/CN111696687B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16HHEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H80/00ICT specially adapted for facilitating communication between medical practitioners or patients, e.g. for collaborative diagnosis, therapy or health monitoring
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/0002Inspection of images, e.g. flaw detection
    • G06T7/0012Biomedical image inspection
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30004Biomedical image processing
    • G06T2207/30101Blood vessel; Artery; Vein; Vascular

Landscapes

  • Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Medical Informatics (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Epidemiology (AREA)
  • Pathology (AREA)
  • Primary Health Care (AREA)
  • Biomedical Technology (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • Quality & Reliability (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)

Abstract

本发明提供一种采用手机摄像头实现微循环及脉搏同步监测的方法,包括如下步骤:将两束波长不相同的激光通过扩束镜发射到待测部位,利用波长较短的激光照射在皮表组织所产生的散斑图像进行分析得到由心脏搏动所产生的脉搏波信号,利用波长较长的激光采集因血红细胞位移所产生的散斑图像,进而分析得到血流灌注量和节律周期,采用手机相机采集连续时间的多帧的散斑图像,然后对每帧图像红绿通道进行处理,并基于散斑信号的空间分布特性进行上采样插值来增强散斑信号,最后对红色通道的图像采用散斑衬比分析,对绿色通道的图像采用灰度共生矩阵分析,从而得到微循环血流灌注量和节律周期及脉搏的变化规律。

Description

一种采用手机摄像头实现微循环及脉搏同步监测的方法
技术领域
本发明属于生物医学成像领域,具体涉及一种采用手机摄像头实现微循环及脉搏同步监测的方法。
背景技术
监测心脏搏动及微循环系统对人体的生理状态的研究有着重要的意义。人体是一个复杂的机体,各组织细胞紧密相连,对人体生理状态的监测应将两者联系起来同步监测。
现有的监测心率的方法主要有:血氧法,该方法通过氧合血红蛋白和还原血红蛋白对不同光源的吸收率不同从而得到血氧信息和心率,但是需要穿透组织,故只能在人体指尖和耳垂这种较薄皮肤组织处使用。超声波测量技术,利用超声波在传播途中遇到阻碍而反弹,从而在超声波传感器出接受到带有心率信息的信号,该方法由于超声波在传播途中有很大的衰减和受环境影响较大的原因,准确性较低。电信号法,通过监测窦房结的电信号来监测心率,该方法是接触式检测且需要人体两个监测点。本专利监测心率的方法是,通过使用灰度共生矩阵对皮肤震动产生的连续时间散斑图像分析得到脉搏波信号,进而监测心率变化。
监测微循环的方法主要有:正交偏振光谱技术,它用波长550nm左右的线形偏振光,该波段为血红蛋白的等消光点,借此显像流经该血管的红血球,该方法的不足是使用半定量的分析。测流暗视野成像技术,该方法是正交偏振光谱技术的一种衍生技术,对微循环中的白细胞和红细胞有更高的分辨率,但数据分析部分仍然比较耗时。激光多普勒技术,通过对组织发射激光发生散射后,击中血细胞的激光波长发生了多普勒频移,从而监测血流变化,但是该方法是单点测量且只能测量1mm直径的血管。
若使用多个机器,分别采集心率信息和微循环信息,将获取到的信息进行分析和整合,这会导致操作复杂、价格昂贵、便携性较差,而且由于测量方法和部位的不同在比较两种生理信息的时候可能出现误差。因此,需要设计一种能够解决上述问题的方法。
发明内容
本发明是为了解决上述问题而进行的,目的在于提供一种采用手机摄像头实现微循环及脉搏同步监测的方法。
本发明提供了一种采用手机摄像头实现微循环及脉搏同步监测的方法,具有这样的特征,包括如下步骤:步骤1,同时采用一束长波长的第一激光和一束短波长的第二激光分别通过第一扩束镜和第二扩束镜对待测部位进行照射,而后用手机相机对激光照射的部位进行拍照,采集连续N帧的图像,生成连续时间轴的带有散斑信息的图像,并记录采集时间t和帧数N;步骤2,将采集到的图像输入电脑,而后采用matlab软件进行分析,对采集到的图像分别提取红色通道即r通道和绿色通道即g通道,得到N帧r通道灰度图像和N帧g通道灰度图像,并对r通道灰度图像和g通道灰度图像使用公式(1)进行上采样插值,得到插值后的r通道灰度图像和插值后的g通道灰度图像;步骤3,选定一个感兴趣区域ROI进行分析,对通过步骤(2)处理后的N帧r通道灰度图像的感兴趣区域ROI,使用统计学方法进行计算,通过公式(2)计算出血流相对速度V,并以时间t为横坐标,速度V为纵坐标绘制曲线,该曲线为微循环的血流速度变化曲线,记为Vn;步骤4,对经过步骤2处理后的N帧g通道灰度图像的感兴趣区域ROI,计算当步长为1,方向为0°,45°,90°,135°时的灰度共生矩阵,分别记为A1,A2,A3,A4;步骤5,采用公式(3)计算灰度共生矩阵A1,A2,A3,A4的各元素值的平方和,分别记为Asm1,Asm2,Asm3,Asm4;步骤6,通过公式(4)计算灰度共生矩阵A1,A2,A3,A4的相关性特征值,分别记为Cor1,Cor2,Cor3,Cor4;步骤7,将Cor1,Cor2,Cor3,Cor4进行归一化,将其结果映射到(0,1)之间,分别记为Co1,Co2,Co3,Co4;步骤8,采用公式(5)计算脉搏波信号强度Ch,根据N帧g通道灰度图像的感兴趣区域ROI内Ch的变化,以时间t为横坐标,Ch为纵坐标绘制曲线,该曲线为脉搏波信号曲线,记为Mn;步骤9,将脉搏波信号曲线Mn中峰值对应的时间点间隔作为脉搏波周期即心动周期,记为TMn,将血流速度变化曲线Vn中峰值对应的时间点间隔作为微循环节律周期,记为TVn;步骤10,根据心动周期TMn以及微循环节律周期TVn计算人体循环系统的指标Zn,而后根据该指标对微循环及脉搏进行同步监测。
在本发明提供的采用手机摄像头实现微循环及脉搏同步监测的方法中,还可以具有这样的特征:其中,步骤1中的第一激光的波长为632nm,第二激光的波长为530nm。
在本发明提供的采用手机摄像头实现微循环及脉搏同步监测的方法中,还可以具有这样的特征:其中,步骤2中的公式(1)为:
Figure BDA0002600650050000041
式(1)中,L1为使用双线性插值法计算出的插值点,N1为使用最近邻插值法计算出的插值点,e为自然常数,σ为插值点周围像素的灰度值的标准差,将σ作为两种插值法的权值进行计算,G1为计算出的新插值点。
在本发明提供的采用手机摄像头实现微循环及脉搏同步监测的方法中,还可以具有这样的特征:其中,步骤3中的所述公式(2)为:
Figure BDA0002600650050000042
式(2)中,
Figure BDA0002600650050000043
为散斑空间强度均值,σs为散斑空间强度的标准偏差。
在本发明提供的采用手机摄像头实现微循环及脉搏同步监测的方法中,还可以具有这样的特征:其中,步骤5中的公式(3)为:
Figure BDA0002600650050000044
式(3)中,P(i,j)为灰度共生矩阵中的元素值的平方和,i为该元素所在的行数,j为该元素所在的列数,k为灰度共生矩阵的总行数和总列数。
在本发明提供的采用手机摄像头实现微循环及脉搏同步监测的方法中,还可以具有这样的特征:其中,步骤6中的公式(4)为:
Figure BDA0002600650050000051
式(4)中,
Figure BDA0002600650050000052
Figure BDA0002600650050000053
且k为灰度共生矩阵的总行数和总列数,μx
Figure BDA0002600650050000054
的均值即Px(1),Px(2),……,Px(k)的均值,μy
Figure BDA0002600650050000055
的均值即Py(1),Py(2),……,Py(k)的均值,σx为Px(i)的方差,σy为Py(j)的方差,其中,x为横向计算的标记,y为纵向计算的标记。
在本发明提供的采用手机摄像头实现微循环及脉搏同步监测的方法中,还可以具有这样的特征:其中,步骤8中的公式(5)为:
Ch=(Co1Asm1+Co2Asm2+Co3Asm3+Co4Asm4)*0.0001 (5)。
在本发明提供的采用手机摄像头实现微循环及脉搏同步监测的方法中,还可以具有这样的特征:其中,步骤10中的人体循环系统的指标Zn的计算公式为:
Figure BDA0002600650050000056
发明的作用与效果
根据本发明所涉及的采用手机摄像头实现微循环及脉搏同步监测的方法,实现了微循环和脉搏波监测的时间同步和空间同步,使用便携的手机设备作为采集设备,并根据散斑信号空间分布特性进行上采样插值来增强散斑信号,进而得到微循环的血流速度及节律周期,然后再利用灰度共生矩阵对带有脉搏波信息的散斑图像进行分析得到脉搏波信号,还利用散斑血流速度变化得到微循环节律周期进而和脉搏波周期进行比较。此外,在人体某些皮肤组织较厚的位置也可以使用本实施例中监测脉搏波的方法。
进一步地,本发明所涉及的采用手机摄像头实现微循环及脉搏同步监测的方法,利用激光散斑的方法,仅需一个手机相机即可对心率和微循环做同步的监测,以较低的成本获取多种生理信息,并且针对手机摄像头像素尺寸较小的问题,采用了上采样插值来提高分辨率,而后再进行散斑信号的分析。
因此,采用本发明的采用手机摄像头实现微循环及脉搏同步监测的方法对微循环及脉搏进行同步监测时,操作简便、便携性高、误差小、准确性较高。
附图说明
图1是本发明的实施例中采用手机摄像头实现微循环及脉搏同步监测的方法的流程图;
图2是本发明的实施例中采用手机摄像头实现微循环及脉搏同步监测的方法的装置示意图。
具体实施方式
为了使本发明实现的技术手段与功效易于明白了解,以下结合实施例及附图对本发明作具体阐述。
实施例:
如图2所示,本实施例中采用的装置包括用于发出第一激光的第一激光器1、用于第一激光通过的第一扩束镜2、用于发出第二激光的第二激光器、用于第二激光通过的第二扩束镜、用于进行拍照的手机相机5以及对图像进行分析处理的计算机6。
如图1所示,本实施例的一种采用手机摄像头实现微循环及脉搏同步监测的方法,包括如下步骤:
步骤1,采用第一激光器射出632nm,15mW的第一激光,第二激光器射出530nm,15mW的第二激光分别通过第一扩束镜和第二扩束镜后,同时打在指端上,手指静止置于平面上,而后用手机相机对激光照射的部位进行拍照,采集连续N帧的图像,生成连续时间轴的带有散斑信息的图像,并记录采集时间t和帧数N。
本实施例中,波长较短的第二激光照射在皮表组织所产生的散斑图像进行分析得到由心脏搏动所产生的脉搏波信号,利用第一激光采集因血红细胞位移所产生的散斑图像,进而分析得到血流灌注量和节律周期。
步骤2,将采集到的图像输入电脑,而后采用matlab软件进行分析,并对采集到的图像分别提取红色通道即r通道和绿色通道即g通道,得到N帧r通道灰度图像和N帧g通道灰度图像,并对r通道灰度图像和g通道灰度图像使用公式(1)进行上采样插值,得到插值后的r通道灰度图像和插值后的g通道灰度图像。
本实施例中,公式(1)为:
Figure BDA0002600650050000081
式中,L1为使用双线性插值法计算出的插值点,N1为使用最近邻插值法计算出的插值点,e为自然常数,约为2.718,σ为插值点周围像素的灰度值的标准差,将σ作为两种插值法的权值进行计算,G1为计算出的新插值点。
步骤3,选定一个感兴趣区域ROI进行分析,对通过步骤(2)处理后的N帧r通道灰度图像的感兴趣区域ROI,使用统计学方法进行计算,通过公式(2)计算出血流相对速度V,并以时间t为横坐标,速度V为纵坐标绘制曲线,该曲线为微循环的血流速度变化曲线,记为Vn
本实施例中,公式(2)为:
Figure BDA0002600650050000082
式中,
Figure BDA0002600650050000083
为散斑空间强度均值,σs为散斑空间强度的标准偏差。
步骤4,对经过步骤2处理后的N帧g通道灰度图像的感兴趣区域ROI,计算当步长为1,方向为0°,45°,90°,135°时的灰度共生矩阵,分别记为A1,A2,A3,A4
步骤5,采用公式(3)计算灰度共生矩阵A1,A2,A3,A4的各元素值的平方和,分别记为Asm1,Asm2,Asm3,Asm4
本实施例中,公式(3)为:
Figure BDA0002600650050000084
式中,P(i,j)为灰度共生矩阵中的元素值的平方和,i为该元素所在的行数,j为该元素所在的列数,k为灰度共生矩阵的总行数和总列数。
步骤6,采用公式(4)计算灰度共生矩阵A1,A2,A3,A4的相关性特征值,分别记为Cor1,Cor2,Cor3,Cor4
本实施例中,公式(4)为:
Figure BDA0002600650050000091
式中,
Figure BDA0002600650050000092
Figure BDA0002600650050000093
且k为灰度共生矩阵的总行数和总列数,μx
Figure BDA0002600650050000094
的均值即Px(1),Px(2),……,Px(k)的均值,μy
Figure BDA0002600650050000095
的均值即Py(1),Py(2),……,Py(k)的均值,σx为Px(i)的方差,σy为Py(j)的方差,其中,x为横向计算的标记,y为纵向计算的标记。
步骤7,将Cor1,Cor2,Cor3,Cor4进行归一化,将其结果映射到(0,1)之间,分别记为Co1,Co2,Co3,Co4
步骤8,采用公式(5)计算···Ch,根据N帧g通道灰度图像的感兴趣区域ROI内Ch的变化,以时间t为横坐标,Ch为纵坐标绘制曲线,该曲线为脉搏波信号曲线,记为Mn
步骤9,将脉搏波信号曲线Mn中峰值对应的时间点间隔作为脉搏波周期即心动周期,记为TMn,将血流速度变化曲线Vn中峰值对应的时间点间隔作为微循环节律周期,记为TVn
本实施例中,脉搏波信号曲线Mn体现了心动周期的过程,血流速度变化曲线Vn体现了微循环的血流速度变化,且公式(5)为:
Ch=(Co1Asm1+Co2Asm2+Co3Asm3+Co4Asm4)*0.0001。
步骤10,根据心动周期TMn以及微循环节律周期TVn计算人体循环系统的指标Zn,而后根据该指标对微循环及脉搏进行同步监测。
本实施例中,人体循环系统的指标Zn的计算公式为:
Figure BDA0002600650050000101
本实施例中,通过分析Zn的变化区间以及Mn曲线和Vn曲线变化的区间,进而监测人体的生理状态,具体地,当Zn变小时,心脏对末梢循环提供的动力相对降低或微循环对末梢循环提供的动力相对增加,当Zn变大时,心脏对末梢循环提供的动力相对增加或微循环对末梢循环提供的动力相对降低。
实施例的作用与效果
根据本实施例所涉及的采用手机摄像头实现微循环及脉搏同步监测的方法,实现了微循环和脉搏波监测的时间同步和空间同步,使用便携的手机设备作为采集设备,并根据散斑信号空间分布特性进行上采样插值来增强散斑信号,进而得到微循环的血流速度及节律周期,然后再利用灰度共生矩阵对带有脉搏波信息的散斑图像进行分析得到脉搏波信号,还利用散斑血流速度变化得到微循环节律周期进而和脉搏波周期进行比较。此外,在人体某些皮肤组织较厚的位置也可以使用本实施例中监测脉搏波的方法。
进一步地,本实施例所涉及的采用手机摄像头实现微循环及脉搏同步监测的方法,利用激光散斑的方法,仅需一个手机相机即可对心率和微循环做同步的监测,以较低的成本获取多种生理信息,并且针对手机摄像头像素尺寸较小的问题,采用了上采样插值来提高分辨率,而后再进行散斑信号的分析。
因此,采用本实施例的采用手机摄像头实现微循环及脉搏同步监测的方法对微循环及脉搏进行同步监测时,操作简便、便携性高、误差小、准确性较高。
上述实施方式为本发明的优选案例,并不用来限制本发明的保护范围。

Claims (2)

1.一种采用手机摄像头实现微循环及脉搏同步监测的方法,其特征在于,包括如下步骤:
步骤1,同时采用一束长波长的第一激光和一束短波长的第二激光分别通过第一扩束镜和第二扩束镜对待测部位进行照射,而后用手机相机对激光照射的部位进行拍照,采集连续N帧的图像,生成连续时间轴的带有散斑信息的图像,并记录采集时间t和帧数N;
步骤2,将采集到的所述图像输入电脑,而后采用matlab软件进行分析,并对采集到的所述图像分别提取红色通道即r通道和绿色通道即g通道,得到N帧r通道灰度图像和N帧g通道灰度图像,并对所述r通道灰度图像和所述g通道灰度图像使用公式(1)进行上采样插值,得到插值后的r通道灰度图像和插值后的g通道灰度图像;
步骤3,选定一个感兴趣区域ROI进行分析,对通过步骤2处理后的N帧所述r通道灰度图像的所述感兴趣区域ROI,使用统计学方法进行计算,通过公式(2)计算出血流相对速度V,并以时间t为横坐标,速度V为纵坐标绘制曲线,该曲线为微循环的血流速度变化曲线,记为Vn
步骤4,对经过步骤2处理后的N帧所述g通道灰度图像的所述感兴趣区域ROI,计算当步长为1,方向为0°,45°,90°,135°时的灰度共生矩阵,分别记为A1,A2,A3,A4
步骤5,采用公式(3)计算所述灰度共生矩阵A1,A2,A3,A4的各元素值的平方和,分别记为Asm1,Asm2,Asm3,Asm4
步骤6,通过公式(4)计算所述灰度共生矩阵A1,A2,A3,A4的相关性特征值,分别记为Cor1,Cor2,Cor3,Cor4
步骤7,将Cor1,Cor2,Cor3,Cor4进行归一化,将其结果映射到(0,1)之间,分别记为Co1,Co2,Co3,Co4
步骤8,采用公式(5)计算脉搏波信号强度Ch,根据N帧所述g通道灰度图像的所述感兴趣区域ROI内Ch的变化,以时间t为横坐标,Ch为纵坐标绘制曲线,该曲线为脉搏波信号曲线,记为Mn
步骤9,将所述脉搏波信号曲线Mn中峰值对应的时间点间隔作为脉搏波周期即心动周期,记为TMn,将所述血流速度变化曲线Vn中峰值对应的时间点间隔作为微循环节律周期,记为TVn
步骤10,根据所述心动周期TMn以及所述微循环节律周期TVn计算人体循环系统的指标Zn,而后根据该指标对微循环及脉搏进行同步监测,
其中,所述步骤2中的所述公式(1)为:
Figure FDA0003501269540000021
式(1)中,L1为使用双线性插值法计算出的插值点,N1为使用最近邻插值法计算出的插值点,e为自然常数,σ为插值点周围像素的灰度值的标准差,将σ作为两种插值法的权值进行计算,G1为计算出的新插值点,
所述步骤3中的所述公式(2)为:
Figure FDA0003501269540000022
式(2)中,
Figure FDA0003501269540000023
为散斑空间强度均值,σs为散斑空间强度的标准偏差,
所述步骤5中的公式(3)为:
Figure FDA0003501269540000031
式(3)中,P(i,j)为灰度共生矩阵中的元素值的平方和,i为该元素所在的行数,j为该元素所在的列数,k为灰度共生矩阵的总行数和总列数,
所述步骤6中的公式(4)为:
Figure FDA0003501269540000032
式(4)中,
Figure FDA0003501269540000033
Figure FDA0003501269540000034
且k为灰度共生矩阵的总行数和总列数,μx
Figure FDA0003501269540000035
的均值即Px(1),Px(2)……Px(k)的均值,μy
Figure FDA0003501269540000036
的均值即Py(1),Py(2)……Py(k)的均值,σx为Px(i)的方差,σy为Py(j)的方差,其中,x为横向计算的标记,y为纵向计算的标记,
所述步骤8中的公式(5)为:
Ch=(Co1Asm1+Co2Asm2+Co3Asm3+Co4Asm4)*0.0001 (5),
所述步骤10中的所述人体循环系统的所述指标Zn的计算公式为:
Figure FDA0003501269540000037
2.根据权利要求1所述的采用手机摄像头实现微循环及脉搏同步监测的方法,其特征在于:
其中,所述步骤1中的所述第一激光的波长为632nm,第二激光的波长为530nm。
CN202010722872.4A 2020-07-24 2020-07-24 一种采用手机摄像头实现微循环及脉搏同步监测的方法 Active CN111696687B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010722872.4A CN111696687B (zh) 2020-07-24 2020-07-24 一种采用手机摄像头实现微循环及脉搏同步监测的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010722872.4A CN111696687B (zh) 2020-07-24 2020-07-24 一种采用手机摄像头实现微循环及脉搏同步监测的方法

Publications (2)

Publication Number Publication Date
CN111696687A CN111696687A (zh) 2020-09-22
CN111696687B true CN111696687B (zh) 2022-03-25

Family

ID=72486802

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010722872.4A Active CN111696687B (zh) 2020-07-24 2020-07-24 一种采用手机摄像头实现微循环及脉搏同步监测的方法

Country Status (1)

Country Link
CN (1) CN111696687B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113782228B (zh) * 2021-09-28 2024-02-27 西安航空学院 一种远程中医诊疗系统

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105653876A (zh) * 2016-01-22 2016-06-08 中国科学院微电子研究所 基于三维成像的脉搏波图像采集与存储方法
TWM542444U (zh) * 2016-11-08 2017-06-01 黃勝榮 穿戴式多功能脈搏波診斷與分析裝置
WO2020006454A1 (en) * 2018-06-28 2020-01-02 Children's National Medical Center Methods and system for dye-free visualization of blood flow and tissue perfusion in laparoscopy
CN110706826A (zh) * 2019-10-29 2020-01-17 江苏大学 一种基于视频图像的非接触式的实时多人心率和血压测量方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2012096878A2 (en) * 2011-01-10 2012-07-19 East Carolina University Methods, systems and computer program products for non-invasive determination of blood flow distribution using speckle imaging techniques and hemodynamic modeling

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105653876A (zh) * 2016-01-22 2016-06-08 中国科学院微电子研究所 基于三维成像的脉搏波图像采集与存储方法
TWM542444U (zh) * 2016-11-08 2017-06-01 黃勝榮 穿戴式多功能脈搏波診斷與分析裝置
WO2020006454A1 (en) * 2018-06-28 2020-01-02 Children's National Medical Center Methods and system for dye-free visualization of blood flow and tissue perfusion in laparoscopy
CN110706826A (zh) * 2019-10-29 2020-01-17 江苏大学 一种基于视频图像的非接触式的实时多人心率和血压测量方法

Also Published As

Publication number Publication date
CN111696687A (zh) 2020-09-22

Similar Documents

Publication Publication Date Title
US9741111B2 (en) Subject information acquiring device and subject information acquiring method
US20190133451A1 (en) Apparatus for acquiring biofunctional information, method for acquiring biofunctional information, and program therefor
KR101273585B1 (ko) 초음파 영상 장치 및 초음파 영상 표시방법
US11471130B2 (en) Method and ultrasound system for shear wave elasticity imaging
CN103327904B (zh) 超声波摄像装置和超声波摄像方法
CA2319458A1 (en) Transabdominal examination, monitoring and imaging of tissue
CN1596832A (zh) 超声波诊断装置和图象处理装置
US20220087651A1 (en) Systems and Methods for Microvessel Ultrasound Imaging
CN109752377B (zh) 一种分光式双模态投影层析组织血管成像装置及方法
US20220104794A1 (en) Method and ultrasound system for shear wave elasticity imaging
US20230172565A1 (en) Systems, devices, and methods for developing a model for use when performing oximetry and/or pulse oximetry and systems, devices, and methods for using a fetal oximetry model to determine a fetal oximetry value
CN111696687B (zh) 一种采用手机摄像头实现微循环及脉搏同步监测的方法
JP6415650B2 (ja) 機能情報取得装置、機能情報取得方法、及びプログラム
Zhang et al. Noncontact heart rate measurement using a webcam, based on joint blind source separation and a skin reflection model: For a wide range of imaging conditions
Noda et al. Ultrasound imaging with a flexible probe based on element array geometry estimation using deep neural network
AU2021350834A9 (en) Systems, devices, and methods for developing a fetal oximetry model for use to determine a fetal oximetry value
CN112168210B (zh) 医学图像处理终端、超声诊断设备和胎儿图像的处理方法
KR20230097482A (ko) Rf신호로부터 혈류의 속력을 측정하는 방법
CN105050505B (zh) 用于超声微钙化检测的波束形成技术
Wu et al. Non-invasive assessment of arterial distension waveforms using gradient-based Hough transform and power Doppler ultrasound imaging
Hafthorsdottir et al. Automatic extraction of muscle fascicle pennation angle from raw ultrasound data
WO2023047601A1 (ja) 画像生成方法、画像生成プログラムおよび画像生成装置
WO2018165992A1 (zh) 一种多通道脑功能成像装置及方法
WO2019163172A1 (ja) 超音波信号処理装置、超音波診断装置および超音波信号演算処理方法
Sathyanarayana et al. Blood flow recovery from subsampled data in photoacoustic microscopy

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