CN105769162B - 一种刮痧疗效判断方法和装置 - Google Patents

一种刮痧疗效判断方法和装置 Download PDF

Info

Publication number
CN105769162B
CN105769162B CN201410802495.XA CN201410802495A CN105769162B CN 105769162 B CN105769162 B CN 105769162B CN 201410802495 A CN201410802495 A CN 201410802495A CN 105769162 B CN105769162 B CN 105769162B
Authority
CN
China
Prior art keywords
blood flow
scraping
region
target
flow area
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
CN201410802495.XA
Other languages
English (en)
Other versions
CN105769162A (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.)
Shenzhen Institute of Advanced Technology of CAS
Original Assignee
Shenzhen Institute of Advanced Technology of CAS
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 Shenzhen Institute of Advanced Technology of CAS filed Critical Shenzhen Institute of Advanced Technology of CAS
Priority to CN201410802495.XA priority Critical patent/CN105769162B/zh
Publication of CN105769162A publication Critical patent/CN105769162A/zh
Application granted granted Critical
Publication of CN105769162B publication Critical patent/CN105769162B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)

Abstract

本发明实施例公开了一种刮痧疗效判断方法,包括:多次获取针对目标刮痧区域刮痧前和刮痧后的血流变化量;分析所述目标刮痧区域的所述血流变化量的变化趋势;根据所述变化趋势判断所述目标刮痧区域的疗效。相应地,本发明实施例还公开了一种刮痧疗效判断装置。采用本发明实施例,可以通过客观的依据判断刮痧疗效,提高诊断的可靠性。

Description

一种刮痧疗效判断方法和装置
技术领域
本发明涉及医疗诊断技术领域,尤其涉及一种刮痧疗效判断方法和装置。
背景技术
“刮痧”是中医外治疗法之一。其机理是通过器械作用,刺激穴位、皮肤经络,将皮下乃至深层组织、内脏之邪气呈现于表、通达于外,从而达到祛除邪气、疏通经络、行气活血、增强脏腑功能、调动卫气的作用。从现代医学观点来看,其机理可能在于它使微细血管扩大或破裂,造成充血或瘀血,化为一种刺激源,形成一种反射,以活跃生理机能,改善血液循环,促进细胞代谢,促进毒素排泄,增强机体免疫力,以达到抗御疾病,促进身体康复的作用。
激光散斑血流成像技术是基于激光散斑对比分析技术(英文:Laser SpeckleContrast Analysis,简称:LASCA),可对大面积组织进行实时的血流动态成像监测。现有的刮痧疗效判断多采用用户的体验和医生对刮痧处的肉眼观察,存在的缺点主要是缺乏客观的依据,容易发生误判。
发明内容
本发明实施例提供一种刮痧疗效判断方法和装置,可以通过客观的依据判断刮痧疗效,提高诊断的可靠性。
本发明实施例提供了一种刮痧疗效判断方法,包括:
多次获取针对目标刮痧区域刮痧前和刮痧后的血流变化量;
分析所述目标刮痧区域的所述血流变化量的变化趋势;
根据所述变化趋势判断所述目标刮痧区域的疗效。
其中,所述获取针对目标刮痧区域刮痧前和刮痧后的血流变化量包括:
分别获取所述目标刮痧区域刮痧前和刮痧后的激光散斑血流量图;
根据获取的所述目标刮痧区域刮痧前和刮痧后的激光散斑血流量图,计算所述目标刮痧区域刮痧前和刮痧后的血流变化量。
进一步地,所述分别获取所述目标刮痧区域刮痧前和刮痧后的激光散斑血流量图包括:
通过激光发射装置发射激光,分别照射刮痧前和刮痧后的所述目标刮痧区域;
通过图像采集装置,分别获取由刮痧前和刮痧后的所述目标刮痧区域反射回的图像信号;
通过激光散斑成像算法,分别将刮痧前和刮痧后的所述目标刮痧区域的原始图转换为激光散斑血流图,所述目标刮痧区域的原始图根据所述目标刮痧区域的图像信号确定。
再进一步地,所述目标刮痧区域的原始图包括多个像素区域;
所述通过激光散斑成像算法,分别将刮痧前和刮痧后的所述目标刮痧区域的原始图转换为激光散斑血流图包括:
根据公式分别计算刮痧前和刮痧后的所述原始图中的每个像素区域的激光散斑衬比值K,其中,所述σs为所述像素区域的光强的平均标准偏差,所述<I>为所述像素区域的光强的平均强度;
根据公式和Pn=(2T/τc)*Μ,分别计算刮痧前和刮痧后的所述原始图的每个像素区域的灰度值,所述灰度值用于指示所述像素区域的血流速度,T为所述目标刮痧区域的曝光时间,τc为散斑强度波动的自相关时间,Pn为所述像素区域的灰度值,Μ为预设常量;
根据刮痧前和刮痧后的所述原始图中的每个像素区域的灰度值,分别将刮痧前和刮痧后的所述原始图转换为激光散斑血流图。
其中,所述激光散斑血流量图包括多个像素区域,每个所述像素区域对应一个灰度值;
所述根据获取的所述目标刮痧区域刮痧前和刮痧后的激光散斑血流量图,计算所述目标刮痧区域刮痧前和刮痧后的血流变化量包括:
在所述目标刮痧区域刮痧后的所述激光散斑血流量图上,确定高血流区域和低血流区域,所述高血流区域和所述低血流区域包括多个像素区域;
根据所述高血流区域和所述低血流区域的每个所述像素区域的灰度值,分别计算所述高血流区域和所述低血流区域的平均血流灌注量;
在所述目标刮痧区域刮痧前的所述激光散斑血流量图上,分别确定与所述目标刮痧区域刮痧后的所述高血流区域和所述低血流区域位置相同的第一血流区域和第二血流区域;
根据所述第一血流区域和所述第二血流区域的每个所述像素区域的灰度值,分别计算所述第一血流区域和所述第二血流区域的平均血流灌注量;
根据公式计算所述目标刮痧区域刮痧前和刮痧后的血流变化量,其中,所述Ε为所述血流变化量,所述PR2T2为刮痧后的所述高血流区域的平均血流灌注量,所述PR1T2为刮痧后的所述低血流区域的平均血流灌注量,所述PR2T1为刮痧前的所述第一血流区域的平均血流灌注量,所述PR1T1为刮痧前的所述第二血流区域的平均血流灌注量。
进一步地,所述在所述目标刮痧区域刮痧后的所述激光散斑血流量图上,确定高血流区域包括:
显示所述目标刮痧区域刮痧后的所述激光散斑血流量图;
获取针对显示的所述激光散斑血流量图输入的高血流区域选择指令;
根据所述高血流区域选择指令,生成兴趣区域,所述兴趣区域包括核心区域和扩散区域,所述扩散区域在所述核心区域的四周;
分别计算所述核心区域和所述扩散区域的平均血流灌注量;
若所述扩散区域的平均血流灌注量小于或等于所述核心区域的平均血流灌注量的预设比值,则显示第一预设颜色;
若所述扩散区域的平均血流灌注量大于所述核心区域的平均血流灌注量的所述预设比值,则显示第二预设颜色;
接收针对所述兴趣区域输入的确认指令,并将所述兴趣区域的核心区域在所述目标刮痧区域刮痧后的所述激光散斑血流量图上映射的区域确定为所述高血流区域。
再进一步地,所述分析所述目标刮痧区域的所述血流变化量的变化趋势包括:
以针对所述目标刮痧区域的刮痧次数为横坐标,以所述目标刮痧区域的血流变化量为纵坐标绘制曲线;
所述根据所述变化趋势判断所述目标刮痧区域的疗效包括:
若所述曲线单调下降,则判定所述目标刮痧区域处于恢复状态;
若所述曲线保持水平发展,则判定所述目标刮痧区域趋于稳定。
相应地,本发明实施例还提供了一种刮痧疗效判断装置,包括:
血流变化量获取模块,用于获取针对目标刮痧区域刮痧前和刮痧后的血流变化量;
变化趋势分析模块,用于分析所述血流变化量获取模块多次获取的目标刮痧区域的所述血流变化量的变化趋势;
疗效诊断模块,用于根据所述变化趋势判断所述目标刮痧区域的疗效。
其中,所述血流变化量获取模块包括:
激光散斑血流量图获取单元,用于获取所述目标刮痧区域刮痧前和刮痧后的激光散斑血流量图;
血流变化量获取单元,用于根据获取的所述目标刮痧区域刮痧前和刮痧后的激光散斑血流量图,计算所述目标刮痧区域刮痧前和刮痧后的血流变化量。
进一步地,所述激光散斑血流量图获取单元包括:
激光发射子单元,用于通过激光发射装置发射激光,分别照射刮痧前和刮痧后的所述目标刮痧区域;
图像采集子单元,用于通过图像采集装置,分别获取由刮痧前和刮痧后的所述目标刮痧区域反射回的图像信号;
图像处理子单元,用于通过激光散斑成像算法,分别将刮痧前和刮痧后的所述目标刮痧区域的原始图转换为激光散斑血流图,所述目标刮痧区域的原始图根据所述目标刮痧区域的图像信号确定。
再进一步地,所述目标刮痧区域的原始图包括多个像素区域;
所述图像处理子单元具体用于:
根据公式分别计算刮痧前和刮痧后的所述原始图中的每个像素区域的激光散斑衬比值K,其中,所述σs为所述像素区域的光强的平均标准偏差,所述<I>为所述像素区域的光强的平均强度;
根据公式和Pn=(2T/τc)*Μ,分别计算刮痧前和刮痧后的所述原始图的每个像素区域的灰度值,所述灰度值用于指示所述像素区域的血流速度,T为所述目标刮痧区域的曝光时间,τc为散斑强度波动的自相关时间,Pn为所述像素区域的灰度值,Μ为预设常量;
根据刮痧前和刮痧后的所述原始图中的每个像素区域的灰度值,分别将刮痧前和刮痧后的所述原始图转换为激光散斑血流图。
其中,所述激光散斑血流量图包括多个像素区域,每个所述像素区域对应一个灰度值;
所述血流变化量获取单元包括:
区域确定子单元,用于在所述目标刮痧区域刮痧后的所述激光散斑血流量图上,确定高血流区域和低血流区域,所述高血流区域和所述低血流区域包括多个像素区域;
平均血流灌注量计算子单元,用于根据所述高血流区域和所述低血流区域的每个所述像素区域的灰度值,分别计算所述高血流区域和所述低血流区域的平均血流灌注量;
所述区域确定子单元还用于:
在所述目标刮痧区域刮痧前的所述激光散斑血流量图上,分别确定与所述目标刮痧区域刮痧后的所述高血流区域和所述低血流区域位置相同的第一血流区域和第二血流区域;
所述平均血流灌注量计算子单元还用于:
根据所述第一血流区域和所述第二血流区域的每个所述像素区域的灰度值,分别计算所述第一血流区域和所述第二血流区域的平均血流灌注量;
血流变化量计算子单元,用于根据公式计算所述目标刮痧区域刮痧前和刮痧后的血流变化量,其中,所述Ε为所述血流变化量,所述PR2T2为刮痧后的所述高血流区域的平均血流灌注量,所述PR1T2为刮痧后的所述低血流区域的平均血流灌注量,所述PR2T1为刮痧前的所述第一血流区域的平均血流灌注量,所述PR1T1为刮痧前的所述第二血流区域的平均血流灌注量。
进一步地,所述区域确定子单元具体用于:
显示所述目标刮痧区域刮痧后的所述激光散斑血流量图;
获取针对显示的所述激光散斑血流量图输入的高血流区域选择指令;
根据所述高血流区域选择指令,生成兴趣区域,所述兴趣区域包括核心区域和扩散区域,所述扩散区域在所述核心区域的四周;
分别计算所述核心区域和所述扩散区域的平均血流灌注量;
若所述扩散区域的平均血流灌注量小于或等于所述核心区域的平均血流灌注量的预设比值,则显示第一预设颜色;
若所述扩散区域的平均血流灌注量大于所述核心区域的平均血流灌注量的所述预设比值,则显示第二预设颜色;
接收针对所述兴趣区域输入的确认指令,并将所述兴趣区域的核心区域在所述目标刮痧区域刮痧后的所述激光散斑血流量图上映射的区域确定为所述高血流区域。
再进一步地,所述变化趋势分析模块用于:
以针对所述目标刮痧区域的刮痧次数为横坐标,以所述目标刮痧区域的血流变化量为纵坐标绘制曲线;
所述疗效诊断模块具体用于:
若所述曲线单调下降,则判定所述目标刮痧区域处于恢复状态;
若所述曲线保持水平发展,则判定所述目标刮痧区域趋于稳定。
由上可见,本发明实施例通过多次获取针对目标刮痧区域刮痧前和刮痧后的血流变化量,分析所述目标刮痧区域的所述血流变化量的变化趋势,根据所述变化趋势判断所述目标刮痧区域的疗效,实现了通过客观的依据判断刮痧疗效,提高了诊断的可靠性。
附图说明
为了更清楚地说明本发明实施例中的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1是本发明实施例提供的一种刮痧疗效判断方法的流程图;
图2是本发明实施例提供的一种激光散斑血流量图获取方法的流程图;
图3是本发明实施例提供的一种刮痧疗效判断装置的结构示意图;
图4是本发明实施例提供的一种血流变化量获取模块的结构示意图;
图5是本发明实施例提供的一种激光散斑血流量图获取单元的结构示意图;
图6是本发明实施例提供的一种血流变化量获取单元的结构示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
请参阅图1,图1是本发明实施例提供的一种刮痧疗效判断方法的流程图。如图1所示本实施例中的刮痧疗效判断方法的流程可以包括:
S101,多次获取针对目标刮痧区域刮痧前和刮痧后的血流变化量。
其中,每次获取针对目标刮痧区域刮痧前和刮痧后的血流变化量具体可以包括以下步骤:
分别获取所述目标刮痧区域刮痧前和刮痧后的激光散斑血流量图;所述激光散斑血流量图可以反映所述目标刮痧区域内,浅表组织血流灌注量的分布情况,通过所述激光散斑血流量图可以直观地观测所述目标刮痧区域不同部位的血流量高低;
根据获取的所述目标刮痧区域刮痧前和刮痧后的激光散斑血流量图,计算所述目标刮痧区域刮痧前和刮痧后的血流变化量。
其中,获取所述目标刮痧区域刮痧后的激光散斑血流量图具体可以包括以下步骤:
通过激光发射装置发射激光,照射刮痧后的所述目标刮痧区域。具体地,所述激光发射装置可以采用低功率的近红外激光发射装置,作为成像的照射光源,波长可以为785nm左右,使入射光束均匀扩束,在距离所述目标刮痧区域约15cm左右的位置,入射光近似均匀的照射在直径约20cm左右的圆形范围内;
通过图像采集装置,获取由刮痧后的所述目标刮痧区域反射回的图像信号。具体地,所述图像采集装置主要可以由CCD激光探测相机组成,可以采用高分辨率和高采样频率的CCD相机,接收由所述目标刮痧区域反射回来的图像信号,保证了成像的高时空响应特性;
通过激光散斑成像算法,将刮痧后的所述目标刮痧区域的原始图转换为激光散斑血流图,所述目标刮痧区域的原始图根据所述目标刮痧区域的图像信号确定;需要指出的是,如何根据图像信号获取原始图是本领域技术人员可理解的,在此不再赘述;
进一步地,所述目标刮痧区域的原始图包括多个像素区域,所述通过激光散斑成像算法,将刮痧后的所述目标刮痧区域的原始图转换为激光散斑血流图具体可以包括以下步骤:
根据公式计算刮痧后的所述原始图中的每个像素区域的激光散斑衬比值K,其中,所述σs为所述像素区域的光强的平均标准偏差,所述<I>为所述像素区域的光强的平均强度;
具体地,激光照射在血细胞上,反射回来的光信号由于相互干涉,会呈现明暗相间的模糊样式,称为激光散斑血流图;
血细胞静止时,所述激光散斑血流图不会变化;而血细胞流动起来后,所述激光散斑血流图会随时间变化,在曝光时间固定的情况下,所述激光散斑血流图模糊程度越低,血流速度越快;所述激光散斑血流图模糊程度越高,血流速度越慢。用一个空间统计量来描述所述激光散斑血流图的模糊程度,即光强的平均标准偏差与平均强度的比值,即激光散斑衬比值K;
根据公式和Pn=(2T/τc)*Μ,计算刮痧后的所述原始图的每个像素区域的灰度值,所述灰度值用于指示所述像素区域的血流速度,T为所述目标刮痧区域的曝光时间,τc为散斑强度波动的自相关时间,Pn为所述像素区域的灰度值,Μ为预设常量;
根据刮痧后的所述原始图中的每个像素区域的灰度值,分别将刮痧前和刮痧后的所述原始图转换为激光散斑血流图。
可理解的是,获取所述目标刮痧区域刮痧前和刮痧后的激光散斑血流量图的方法流程一样,本发明实施例只具体描述获取目标刮痧区域刮痧后的激光散斑血流量图,本领域技术人员很容易根据获取目标刮痧区域刮痧后的激光散斑血流量图的方法,获取目标刮痧区域刮痧前的激光散斑血流量图。
其中,所述激光散斑血流量图包括多个像素区域,每个所述像素区域对应一个灰度值,所述根据获取的所述目标刮痧区域刮痧前和刮痧后的激光散斑血流量图,计算所述目标刮痧区域刮痧前和刮痧后的血流变化量具体可以包括以下步骤:
在所述目标刮痧区域刮痧后的所述激光散斑血流量图上,确定高血流区域和低血流区域,所述高血流区域和所述低血流区域包括多个像素区域;
根据所述高血流区域和所述低血流区域的每个所述像素区域的灰度值,分别计算所述高血流区域和所述低血流区域的平均血流灌注量;具体地,首先可以根据公式Pnt=Pn*T分别计算所述高血流区域和所述低血流区域中每个像素区域的血流灌注量,再根据公式PRT=ΣPnt/n分别计算所述高血流区域和所述低血流区域的平均血流灌注量,其中,所述Pn为像素区域的灰度值,T为目标刮痧区域的曝光时间,Pnt为像素区域的血流灌注量,PRT为单位面积内的平均血流灌注量,n为所述单位面积内的像素区域的个数,所述单位面积可以是所述高血流区域或低血流区域;
在所述目标刮痧区域刮痧前的所述激光散斑血流量图上,分别确定与所述目标刮痧区域刮痧后的所述高血流区域和所述低血流区域位置相同的第一血流区域和第二血流区域;具体地,当在所述目标刮痧区域刮痧后的所述激光散斑血流量图上,确定了高血流区域和低血流区域之后,再通过映射在所述目标刮痧区域刮痧前的所述激光散斑血流量图上,分别确定与所述目标刮痧区域刮痧后的所述高血流区域和所述低血流区域位置相同的第一血流区域和第二血流区域;
根据所述第一血流区域和所述第二血流区域的每个所述像素区域的灰度值,分别计算所述第一血流区域和所述第二血流区域的平均血流灌注量;需要指出的是,计算所述第一血流区域和所述第二血流区域的平均血流灌注量与计算所述高血流区域和所述低血流区域的平均血流灌注量的方法相同,在此不再赘述;
根据公式计算所述目标刮痧区域刮痧前和刮痧后的血流变化量,其中,所述Ε为所述血流变化量,所述PR2T2为刮痧后的所述高血流区域的平均血流灌注量,所述PR1T2为刮痧后的所述低血流区域的平均血流灌注量,所述PR2T1为刮痧前的所述第一血流区域的平均血流灌注量,所述PR1T1为刮痧前的所述第二血流区域的平均血流灌注量。
其中,所述在所述目标刮痧区域刮痧后的所述激光散斑血流量图上,确定高血流区域具体可以包括以下步骤:
显示所述目标刮痧区域刮痧后的所述激光散斑血流量图;
获取针对显示的所述激光散斑血流量图输入的高血流区域选择指令;
根据所述高血流区域选择指令,生成兴趣区域,所述兴趣区域包括核心区域和扩散区域,所述扩散区域在所述核心区域的四周;具体地,当获取到针对显示的所述激光散斑血流量图输入的高血流区域选择指令时,可以在当前显示的所述激光散斑血流量图上产生一个兴趣区域,所述兴趣区域包括核心区域和扩散区域,所述扩散区域在所述核心区域的四周;其中,所述兴趣区域在所述激光散斑血流量图上的位置可以随着光标的移动而移动,可选地,所述兴趣区域可以是以所述光标为中心,以5mm为半径的圆形区域;
分别计算所述核心区域和所述扩散区域的平均血流灌注量;需要指出的是,计算所述核心区域和所述扩散区域的平均血流灌注量与与计算所述高血流区域和所述低血流区域的平均血流灌注量的方法相同,在此不再赘述;
若所述扩散区域的平均血流灌注量小于或等于所述核心区域的平均血流灌注量的预设比值,则显示第一预设颜色;具体地,若所述扩散区域的平均血流灌注量小于或等于所述核心区域的平均血流灌注量的预设比值,则说明所述核心区域的平均血流灌注量较高,也就相当于提示医生所述核心区域可以作为高血流区域;
若所述扩散区域的平均血流灌注量大于所述核心区域的平均血流灌注量的所述预设比值,则显示第二预设颜色,具体地,若所述扩散区域的平均血流灌注量大于所述核心区域的平均血流灌注量的预设比值,则说明所述核心区域的平均血流灌注量较低,也就相当于提示医生所述核心区域可以作为低血流区域;
所述第一预设颜色和所述第二预设颜色可以根据实际应用或用户的设置而确定,本发明不做限定;
接收针对所述兴趣区域输入的确认指令,并将所述兴趣区域的核心区域在所述目标刮痧区域刮痧后的所述激光散斑血流量图上映射的区域确定为所述高血流区域。
S102,分析所述目标刮痧区域的所述血流变化量的变化趋势。
一种可选的实施方式,可以以针对所述目标刮痧区域的刮痧次数为横坐标,以所述目标刮痧区域的血流变化量为纵坐标绘制曲线,根据绘制的所述曲线的变化趋势判断疗效。
S103,根据所述变化趋势判断所述目标刮痧区域的疗效。
若所述曲线单调下降,则可以判定所述目标刮痧区域处于恢复状态,可以向医生输出继续治疗以及有效建议;
若所述曲线保持水平发展,则判定所述目标刮痧区域趋于稳定,可以向医生输出暂缓治疗以及康复建议。
在图1所示的实施例中,可以通过多次获取针对目标刮痧区域刮痧前和刮痧后的血流变化量,分析所述目标刮痧区域的所述血流变化量的变化趋势,根据所述变化趋势判断所述目标刮痧区域的疗效,实现了通过客观的依据判断刮痧疗效,提高了诊断的可靠性。
请参阅图2,图2是本发明实施例提供的一种激光散斑血流量图获取方法的流程图。本发明实施例可以用于获取目标刮痧区域刮痧前或刮痧后的激光散斑血流量图,为了便与描述,本发明实施例以获取目标刮痧区域刮痧后的激光散斑血流量图为例。如图2所示本实施例中的激光散斑血流量图获取方法的流程可以包括:
S201,通过激光发射装置发射激光,照射刮痧后的目标刮痧区域。
具体地,所述激光发射装置可以采用低功率的近红外激光发射装置,作为成像的照射光源,波长可以为785nm左右,使入射光束均匀扩束,在距离所述目标刮痧区域约15cm左右的位置,入射光近似均匀的照射在直径约20cm左右的圆形范围内。
S202,通过图像采集装置,获取由刮痧后的所述目标刮痧区域反射回的图像信号。
具体地,所述图像采集装置主要可以由CCD激光探测相机组成,可以采用高分辨率和高采样频率的CCD相机,接收由所述目标刮痧区域反射回来的图像信号,保证了成像的高时空响应特性。
S203,通过激光散斑成像算法,将刮痧后的所述目标刮痧区域的原始图转换为激光散斑血流图,所述目标刮痧区域的原始图根据所述目标刮痧区域的图像信号确定。
所述目标刮痧区域的原始图包括多个像素区域;
首先,根据公式计算刮痧后的所述原始图中的每个像素区域的激光散斑衬比值K,其中,所述σs为所述像素区域的光强的平均标准偏差,所述<I>为所述像素区域的光强的平均强度;
激光照射在血细胞上,反射回来的光信号由于相互干涉,会呈现明暗相间的模糊样式,称为激光散斑血流图;
血细胞静止时,所述激光散斑血流图不会变化;而血细胞流动起来后,所述激光散斑血流图会随时间变化,在曝光时间固定的情况下,所述激光散斑血流图模糊程度越低,血流速度越快;所述激光散斑血流图模糊程度越高,血流速度越慢。用一个空间统计量来描述所述激光散斑血流图的模糊程度,即光强的平均标准偏差与平均强度的比值,即激光散斑衬比值K;
其次,根据公式和Pn=(2T/τc)*Μ,计算刮痧后的所述原始图的每个像素区域的灰度值,所述灰度值用于指示所述像素区域的血流速度,T为所述目标刮痧区域的曝光时间,τc为散斑强度波动的自相关时间,Pn为所述像素区域的灰度值,Μ为预设常量;
根据刮痧后的所述原始图中的每个像素区域的灰度值,分别将刮痧前和刮痧后的所述原始图转换为激光散斑血流图。
可理解的是,获取所述目标刮痧区域刮痧前和刮痧后的激光散斑血流量图的方法流程一样,本发明实施例只具体描述获取目标刮痧区域刮痧后的激光散斑血流量图,本领域技术人员很容易根据获取目标刮痧区域刮痧后的激光散斑血流量图的方法,获取目标刮痧区域刮痧前的激光散斑血流量图。
在图2所示的实施例中,可以通过激光散斑成像算法,将刮痧后的目标刮痧区域的原始图转换为激光散斑血流图,实现了非接触的监测其微循环血流信号。
请参阅图3,图3是本发明实施例提供的一种刮痧疗效判断装置的结构示意图。如图3所示本实施例中的刮痧疗效判断装置3至少可以包括:血流变化量获取模块31、变化趋势分析模块32以及疗效诊断模块33,其中:
血流变化量获取模块31,用于获取针对目标刮痧区域刮痧前和刮痧后的血流变化量。
所述血流变化量获取模块31如图4所示进一步可以包括:激光散斑血流量图获取单元311以及血流变化量获取单元312,其中:
激光散斑血流量图获取单元311,用于获取所述目标刮痧区域刮痧前和刮痧后的激光散斑血流量图。
血流变化量获取单元312,用于根据获取的所述目标刮痧区域刮痧前和刮痧后的激光散斑血流量图,计算所述目标刮痧区域刮痧前和刮痧后的血流变化量。
进一步地,所述激光散斑血流量图获取单元311如图5所示进一步可以包括:激光发射子单元3111、图像采集子单元3112以及图像处理子单元3113,其中:
激光发射子单元3111,用于通过激光发射装置发射激光,分别照射刮痧前和刮痧后的所述目标刮痧区域。
图像采集子单元3112,用于通过图像采集装置,分别获取由刮痧前和刮痧后的所述目标刮痧区域反射回的图像信号。
图像处理子单元3113,用于通过激光散斑成像算法,分别将刮痧前和刮痧后的所述目标刮痧区域的原始图转换为激光散斑血流图,所述目标刮痧区域的原始图根据所述目标刮痧区域的图像信号确定。
其中,所述目标刮痧区域的原始图包括多个像素区域;
所述图像处理子单元3113具体可以用于:
根据公式分别计算刮痧前和刮痧后的所述原始图中的每个像素区域的激光散斑衬比值K,其中,所述σs为所述像素区域的光强的平均标准偏差,所述<I>为所述像素区域的光强的平均强度;
根据公式和Pn=(2T/τc)*Μ,分别计算刮痧前和刮痧后的所述原始图的每个像素区域的灰度值,所述灰度值用于指示所述像素区域的血流速度,T为所述目标刮痧区域的曝光时间,τc为散斑强度波动的自相关时间,Pn为所述像素区域的灰度值,Μ为预设常量;
根据刮痧前和刮痧后的所述原始图中的每个像素区域的灰度值,分别将刮痧前和刮痧后的所述原始图转换为激光散斑血流图。
进一步地,所述激光散斑血流量图包括多个像素区域,每个所述像素区域对应一个灰度值,所述血流变化量获取单元312如图6所示进一步可以包括:区域确定子单元3121、平均血流灌注量计算子单元3122以及血流变化量计算子单元3123,其中:
区域确定子单元3121,用于在所述目标刮痧区域刮痧后的所述激光散斑血流量图上,确定高血流区域和低血流区域,所述高血流区域和所述低血流区域包括多个像素区域。
具体地,区域确定子单元3121可以显示所述目标刮痧区域刮痧后的所述激光散斑血流量图,获取针对显示的所述激光散斑血流量图输入的高血流区域选择指令,根据所述高血流区域选择指令,生成兴趣区域,所述兴趣区域包括核心区域和扩散区域,所述扩散区域在所述核心区域的四周,分别计算所述核心区域和所述扩散区域的平均血流灌注量,若所述扩散区域的平均血流灌注量小于或等于所述核心区域的平均血流灌注量的预设比值,则显示第一预设颜色,若所述扩散区域的平均血流灌注量大于所述核心区域的平均血流灌注量的所述预设比值,则显示第二预设颜色,接收针对所述兴趣区域输入的确认指令,并将所述兴趣区域的核心区域在所述目标刮痧区域刮痧后的所述激光散斑血流量图上映射的区域确定为所述高血流区域。
区域确定子单元3121,同样可以通过以上方式在所述目标刮痧区域刮痧后的所述激光散斑血流量图上确定低血流区域,在此不再赘述。
平均血流灌注量计算子单元3122,用于根据所述高血流区域和所述低血流区域的每个所述像素区域的灰度值,分别计算所述高血流区域和所述低血流区域的平均血流灌注量;
所述区域确定子单元3121还用于:
在所述目标刮痧区域刮痧前的所述激光散斑血流量图上,分别确定与所述目标刮痧区域刮痧后的所述高血流区域和所述低血流区域位置相同的第一血流区域和第二血流区域;
所述平均血流灌注量计算子单元3122还用于:
根据所述第一血流区域和所述第二血流区域的每个所述像素区域的灰度值,分别计算所述第一血流区域和所述第二血流区域的平均血流灌注量;
血流变化量计算子单元3123,用于根据公式计算所述目标刮痧区域刮痧前和刮痧后的血流变化量,其中,所述Ε为所述血流变化量,所述PR2T2为刮痧后的所述高血流区域的平均血流灌注量,所述PR1T2为刮痧后的所述低血流区域的平均血流灌注量,所述PR2T1为刮痧前的所述第一血流区域的平均血流灌注量,所述PR1T1为刮痧前的所述第二血流区域的平均血流灌注量。
变化趋势分析模块32,用于分析所述血流变化量获取模块多次获取的目标刮痧区域的所述血流变化量的变化趋势。
一种可选的实施方式,所述变化趋势分析模块32可以以针对所述目标刮痧区域的刮痧次数为横坐标,以所述目标刮痧区域的血流变化量为纵坐标绘制曲线,根据绘制的所述曲线分析所述目标刮痧区域的血流变化量的变化趋势。
疗效诊断模块33,用于根据所述变化趋势判断所述目标刮痧区域的疗效。
若所述曲线单调下降,疗效诊断模块33可以判定所述目标刮痧区域处于恢复状态;
若所述曲线保持水平发展,疗效诊断模块33判定所述目标刮痧区域趋于稳定。
可理解的是,本实施例的刮痧疗效判断装置3的各功能模块的功能可根据上述方法实施例中的方法具体实现,可以具体对应参考图1至图2方法实施例的相关描述,此处不再赘述。
本发明实施例方法中的步骤可以根据实际需要进行顺序调整、合并和删减。
本发明实施例装置中的模块可以根据实际需要进行合并、划分和删减。
在本说明书的描述中,参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中,对上述术语的示意性表述不是必须针对相同的实施例或示例。而且,描述的具体特征、结构、材料或者特点可以在任一个或多个实施例或示例中以合适的方式结合。此外,在不相互矛盾的情况下,本领域的技术人员可以将本说明书中描述的不同实施例或示例以及不同实施例或示例的特征进行结合和组合。
此外,术语“第一”、“第二”仅用于描述目的,而不能理解为指示或暗示相对重要性或者隐含指明所指示的技术特征的数量。由此,限定有“第一”、“第二”的特征可以明示或者隐含地包括至少一个该特征。在本发明的描述中,“多个”的含义是至少两个,例如两个,三个等,除非另有明确具体的限定。
流程图中或在此以其他方式描述的任何过程或方法描述可以被理解为,表示包括一个或多个用于实现特定逻辑功能或过程的步骤的可执行指令的代码的模块、片段或部分,并且本发明的优选实施方式的范围包括另外的实现,其中可以不按所示出或讨论的顺序,包括根据所涉及的功能按基本同时的方式或按相反的顺序,来执行功能,这应被本发明的实施例所属技术领域的技术人员所理解。
在流程图中表示或在此以其他方式描述的逻辑和/或步骤,例如,可以被认为是用于实现逻辑功能的可执行指令的程序列表,可以具体实现在任何计算机可读介质中,以供指令执行系统、装置或设备(如基于计算机的系统、包括处理器的系统或其他可以从指令执行系统、装置或设备取指令并执行指令的系统)使用,或结合这些指令执行系统、装置或设备而使用。就本说明书而言,"计算机可读介质"可以是任何可以包含、存储、通信、传播或传输程序以供指令执行系统、装置或设备或结合这些指令执行系统、装置或设备而使用的装置。计算机可读介质的更具体的示例(非穷尽性列表)包括以下:具有一个或多个布线的电连接部(电子装置),便携式计算机盘盒(磁装置),随机存取存储器(RAM),只读存储器(ROM),可擦除可编辑只读存储器(EPROM或闪速存储器),光纤装置,以及便携式光盘只读存储器(CDROM)。另外,计算机可读介质甚至可以是可在其上打印所述程序的纸或其他合适的介质,因为可以例如通过对纸或其他介质进行光学扫描,接着进行编辑、解译或必要时以其他合适方式进行处理来以电子方式获得所述程序,然后将其存储在计算机存储器中。
本技术领域的普通技术人员可以理解实现上述实施例方法携带的全部或部分步骤是可以通过程序来指令相关的硬件完成,所述的程序可以存储于一种计算机可读存储介质中,该程序在执行时,包括方法实施例的步骤之一或其组合。
此外,在本发明各个实施例中的各功能单元可以集成在一个处理模块中,也可以是各个单元单独物理存在,也可以两个或两个以上单元集成在一个模块中。上述集成的模块既可以采用硬件的形式实现,也可以采用软件功能模块的形式实现。所述集成的模块如果以软件功能模块的形式实现并作为独立的产品销售或使用时,也可以存储在一个计算机可读取存储介质中。
上述提到的存储介质可以是只读存储器,磁盘或光盘等。尽管上面已经示出和描述了本发明的实施例,可以理解的是,上述实施例是示例性的,不能理解为对本发明的限制,本领域的普通技术人员在本发明的范围内可以对上述实施例进行变化、修改、替换和变型。

Claims (5)

1.一种刮痧疗效判断装置,其特征在于,所述装置包括:
血流变化量获取模块,用于获取针对目标刮痧区域刮痧前和刮痧后的血流变化量;
变化趋势分析模块,用于分析所述血流变化量获取模块多次获取的目标刮痧区域的所述血流变化量的变化趋势;
疗效诊断模块,用于根据所述变化趋势判断所述目标刮痧区域的疗效;
其中,所述血流变化量获取模块具体用于在目标刮痧区域刮痧前和目标刮痧区域刮痧后的激光散斑血流量图上,确定高血流区域和低血流区域,获取所述高血流区域的平均血流灌注量和所述低血流区域的平均血流灌注量;根据所述高血流区域的平均血流灌注量和所述低血流区域的平均血流灌注量计算所述目标刮痧区域刮痧前和刮痧后的血流变化量;
其中,所述血流变化量获取模块包括:
激光散斑血流量图获取单元,用于获取所述目标刮痧区域刮痧前和刮痧后的激光散斑血流量图;
血流变化量获取单元,用于根据获取的所述目标刮痧区域刮痧前和刮痧后的激光散斑血流量图,计算所述目标刮痧区域刮痧前和刮痧后的血流变化量;
其中,所述激光散斑血流量图包括多个像素区域,每个所述像素区域对应一个灰度值;
所述血流变化量获取单元包括:
区域确定子单元,用于在所述目标刮痧区域刮痧后的所述激光散斑血流量图上,确定高血流区域和低血流区域,所述高血流区域和所述低血流区域包括多个像素区域;
平均血流灌注量计算子单元,用于根据所述高血流区域和所述低血流区域的每个所述像素区域的灰度值,分别计算所述高血流区域和所述低血流区域的平均血流灌注量;
所述区域确定子单元还用于:
在所述目标刮痧区域刮痧前的所述激光散斑血流量图上,分别确定与所述目标刮痧区域刮痧后的所述高血流区域和所述低血流区域位置相同的第一血流区域和第二血流区域;
所述平均血流灌注量计算子单元还用于:
根据所述第一血流区域和所述第二血流区域的每个所述像素区域的灰度值,分别计算所述第一血流区域和所述第二血流区域的平均血流灌注量;
血流变化量计算子单元,用于根据公式计算所述目标刮痧区域刮痧前和刮痧后的血流变化量,其中,所述E为所述血流变化量,所述PR2T2为刮痧后的所述高血流区域的平均血流灌注量,所述PR1T2为刮痧后的所述低血流区域的平均血流灌注量,所述PR2T1为刮痧前的所述第一血流区域的平均血流灌注量,所述PR1T1为刮痧前的所述第二血流区域的平均血流灌注量。
2.如权利要求1所述的装置,其特征在于,所述激光散斑血流量图获取单元包括:
激光发射子单元,用于通过激光发射装置发射激光,分别照射刮痧前和刮痧后的所述目标刮痧区域;
图像采集子单元,用于通过图像采集装置,分别获取由刮痧前和刮痧后的所述目标刮痧区域反射回的图像信号;
图像处理子单元,用于通过激光散斑成像算法,分别将刮痧前和刮痧后的所述目标刮痧区域的原始图转换为激光散斑血流图,所述目标刮痧区域的原始图根据所述目标刮痧区域的图像信号确定。
3.如权利要求2所述的装置,其特征在于,所述目标刮痧区域的原始图包括多个像素区域;
所述图像处理子单元具体用于:
根据公式分别计算刮痧前和刮痧后的所述原始图中的每个像素区域的激光散斑衬比值K,其中,所述σs为所述像素区域的光强的平均标准偏差,所述<I>为所述像素区域的光强的平均强度;
根据公式计算所述目标刮痧区域的曝光时间T和所述目标刮痧区域的自相关时间τc的比值,根据所述比值和Pn=(2/Tτc)*Μ,分别计算刮痧前和刮痧后的所述原始图的每个像素区域的灰度值,所述灰度值用于指示所述像素区域的血流速度,T为所述目标刮痧区域的曝光时间,τc为散斑强度波动的自相关时间,Pn为所述像素区域的灰度值,Μ为预设常量;
根据刮痧前和刮痧后的所述原始图中的每个像素区域的灰度值,分别将刮痧前和刮痧后的所述原始图转换为激光散斑血流图。
4.如权利要求1所述的装置,其特征在于,所述区域确定子单元具体用于:
显示所述目标刮痧区域刮痧后的所述激光散斑血流量图;
获取针对显示的所述激光散斑血流量图输入的高血流区域选择指令;
根据所述高血流区域选择指令,生成兴趣区域,所述兴趣区域包括核心区域和扩散区域,所述扩散区域在所述核心区域的四周;
分别计算所述核心区域和所述扩散区域的平均血流灌注量;
若所述扩散区域的平均血流灌注量小于或等于所述核心区域的平均血流灌注量的预设比值,则显示第一预设颜色;
若所述扩散区域的平均血流灌注量大于所述核心区域的平均血流灌注量的所述预设比值,则显示第二预设颜色;
接收针对所述兴趣区域输入的确认指令,并将所述兴趣区域的核心区域在所述目标刮痧区域刮痧后的所述激光散斑血流量图上映射的区域确定为所述高血流区域。
5.如权利要求1-4任一项所述的装置,其特征在于,
所述变化趋势分析模块用于:
以针对所述目标刮痧区域的刮痧次数为横坐标,以所述目标刮痧区域的血流变化量为纵坐标绘制曲线;
所述疗效诊断模块具体用于:
若所述曲线单调下降,则判定所述目标刮痧区域处于恢复状态;
若所述曲线保持水平发展,则判定所述目标刮痧区域趋于稳定。
CN201410802495.XA 2014-12-18 2014-12-18 一种刮痧疗效判断方法和装置 Active CN105769162B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410802495.XA CN105769162B (zh) 2014-12-18 2014-12-18 一种刮痧疗效判断方法和装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410802495.XA CN105769162B (zh) 2014-12-18 2014-12-18 一种刮痧疗效判断方法和装置

Publications (2)

Publication Number Publication Date
CN105769162A CN105769162A (zh) 2016-07-20
CN105769162B true CN105769162B (zh) 2019-05-17

Family

ID=56385931

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410802495.XA Active CN105769162B (zh) 2014-12-18 2014-12-18 一种刮痧疗效判断方法和装置

Country Status (1)

Country Link
CN (1) CN105769162B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108523872A (zh) * 2017-03-06 2018-09-14 上海健康医学院 一种理疗功效评估方法及系统
CN109077712A (zh) * 2018-06-25 2018-12-25 深圳市德力凯医疗设备股份有限公司 一种脑血流的血流变化幅度的显示方法及系统
CN109222952A (zh) * 2018-07-17 2019-01-18 上海健康医学院 一种激光散斑血流灌注成像方法
CN109717861A (zh) * 2018-12-28 2019-05-07 华中科技大学鄂州工业技术研究院 一种血流速度的激光散斑检测方法及检测装置

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102429650A (zh) * 2011-11-10 2012-05-02 华中科技大学 一种激光散斑血流成像衬比分析方法
US20130204112A1 (en) * 2012-02-07 2013-08-08 Laser Associated Sciences, LLC Perfusion assessment using transmission laser speckle imaging
CN103300841A (zh) * 2013-06-13 2013-09-18 上海理工大学 一种快速激光散斑血流成像系统和方法
US20140357990A1 (en) * 2013-05-29 2014-12-04 University Of Washington Through Its Center For Commercialization Methods for Laser Speckle Contrast Imaging of Blood Perfusion

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102429650A (zh) * 2011-11-10 2012-05-02 华中科技大学 一种激光散斑血流成像衬比分析方法
US20130204112A1 (en) * 2012-02-07 2013-08-08 Laser Associated Sciences, LLC Perfusion assessment using transmission laser speckle imaging
US20140357990A1 (en) * 2013-05-29 2014-12-04 University Of Washington Through Its Center For Commercialization Methods for Laser Speckle Contrast Imaging of Blood Perfusion
CN103300841A (zh) * 2013-06-13 2013-09-18 上海理工大学 一种快速激光散斑血流成像系统和方法

Non-Patent Citations (8)

* Cited by examiner, † Cited by third party
Title
Blood perfusion effect of acupuncture on acute facial paralysis observed by laser speckle technique;Wei-zheng Zhong等;《Journal of Acupuncture and Tuina Science》;20140630;第12卷(第3期);160-164 *
Laser speckle contrast imaging for measuring blood flow;J.David Briers;《Optica Applicata》;20071231;139-152 *
Qin-Yan Xu 等.The effects of Scraping therapy on local temperature and blood perfusion volume in healthy subjects.《Evidence-Based Complementary and Alternative Medicine》.2012,第2012卷 *
The effects of Scraping therapy on local temperature and blood perfusion volume in healthy subjects;Qin-Yan Xu 等;《Evidence-Based Complementary and Alternative Medicine》;20121231;第2012卷;1-7 *
刮痧治疗对正常人背部体表温度与血流灌注量的影响;徐青燕等;《2011中国针灸学会年会论文集(摘要)》;20111231;19-24 *
委中穴区刮痧对本经同侧经脉线上皮肤微循环血流灌注量的影响;徐青燕等;《Acupuncture Research》;20130228;第38卷(第1期);52-56 *
徐青燕等.刮痧治疗对正常人背部体表温度与血流灌注量的影响.《2011中国针灸学会年会论文集(摘要)》.2011, *
徐青燕等.委中穴区刮痧对本经同侧经脉线上皮肤微循环血流灌注量的影响.《Acupuncture Research》.2013,第38卷(第1期), *

Also Published As

Publication number Publication date
CN105769162A (zh) 2016-07-20

Similar Documents

Publication Publication Date Title
Kruger et al. Photoacoustic angiography of the breast
Nagae et al. Real-time 3D photoacoustic visualization system with a wide field of view for imaging human limbs
Kazmi et al. Expanding applications, accuracy, and interpretation of laser speckle contrast imaging of cerebral blood flow
Gratton et al. Shedding light on brain function: the event-related optical signal
Deán-Ben et al. Functional optoacoustic human angiography with handheld video rate three dimensional scanner
Deán-Ben et al. Functional optoacoustic imaging of moving objects using microsecond-delay acquisition of multispectral three-dimensional tomographic data
US10603134B2 (en) Touchless advanced image processing and visualization
WO2018102446A9 (en) Single-impulse panoramic photoacoustic computed tomography (sip-pact)
CN105769162B (zh) 一种刮痧疗效判断方法和装置
Volkov et al. Evaluation of blood microcirculation parameters by combined use of laser Doppler flowmetry and videocapillaroscopy methods
JP4495154B2 (ja) 哺乳動物での動的な熱力学的過程を同定し分類してそれらを識別するシステム
Reinhardt et al. VertiGo–a pilot project in nystagmus detection via webcam
He et al. Diffuse optical tomography to monitor the photocoagulation front during interstitial photothermal therapy: Numerical simulations and measurements in tissue-simulating phantoms
Daoudi et al. Two-dimensional spatiotemporal monitoring of temperature in photothermal therapy using hybrid photoacoustic–ultrasound transmission tomography
US20060264713A1 (en) Disease and therapy dissemination representation
Gruber et al. High-frequency ultrasound-guided fluorescence tomography of protoporphyrin IX in subcutaneous tumors
US20200205664A1 (en) Method and system for multi-modal imaging of tissue
Murakoshi et al. Quantitative analysis for peripheral vascularity assessment based on clinical photoacoustic and ultrasound images
Wang et al. Towards real-time diffuse optical tomography for imaging brain functions cooperated with Kalman estimator
Tremoleda et al. Heart-rate sensitive optical coherence angiography for measuring vascular changes due to posttraumatic brain injury in mice
Wang et al. An adaptive extended Kalman filter for fluorescence diffuse optical tomography of tumor pharmacokinetics
Sun et al. Multi-channel photon migration study in visible Chinese human muscle for optical detection of deep vein thrombosis
Abdalla et al. Challenges of medical image processing
Jauer et al. GPU-based real-time generation of large ultrasound volumes from freehand 3D sweeps
Vidal Rosas Advanced tomographic reconstruction methods for diffuse optical imaging of cerebral haemodynamic response.

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant