CN110786842B - 测量舒张期血流速度的方法、装置、系统及存储介质 - Google Patents
测量舒张期血流速度的方法、装置、系统及存储介质 Download PDFInfo
- Publication number
- CN110786842B CN110786842B CN201911066599.8A CN201911066599A CN110786842B CN 110786842 B CN110786842 B CN 110786842B CN 201911066599 A CN201911066599 A CN 201911066599A CN 110786842 B CN110786842 B CN 110786842B
- Authority
- CN
- China
- Prior art keywords
- frame
- flow velocity
- blood flow
- image
- blood vessel
- 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
Images
Classifications
-
- 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
- A61B5/026—Measuring blood flow
- A61B5/0275—Measuring blood flow using tracers, e.g. dye dilution
Abstract
本申请提供了一种测量舒张期血流速度的方法、装置、系统及存储介质。测量舒张期血流速度的方法包括:读取至少一个体位的冠状动脉二维造影图像组;从所述冠状动脉二维造影图像组中,提取感兴趣的血管段;提取所述血管段的中心线;对任意两帧冠状动脉二维造影图像中的造影剂流经所述血管段的时间做差,差值为Δt,且对所述分段中心线做差,差值为ΔL;根据所述ΔL与所述Δt的比值,求解血流速度;选取所述血流速度的最大值,即为舒张期的血流速度。本申请无需通过无创血压仪测量,如此测量的舒张期的血流速度不会受到情绪等外界因素的影响,提高了测量舒张期血流速度的准确度。
Description
技术领域
本发明涉及冠状动脉技术领域,特别是涉及一种测量舒张期血流速度的方法、装置、血管评定参数的计算方法、冠状动脉分析系统及计算机存储介质。
背景技术
世界卫生组织统计,心血管疾病已经成为人类健康的“头号杀手”。近些年,使用血流动力学分析心血管疾病的生理和病理行为也已经成为心血管疾病诊断的一个非常重要的手段。
血液流量和流速作为血流动力学的非常重要的参数,如何准确、便捷地测量血液流量和流速成为广大研究学者研究的重点。
血管评定参数包括:冠状动脉舒张期的血流速度IFR,以及冠状动脉舒张期的微循环阻力指数IFMR等;而IFR和IFMR均需要基于冠状动脉舒张期的血流速度,目前舒张期的血流速度均是通过无创血压仪获得的,存在舒张期的血流速度测量不准确的问题。
发明内容
本发明提供了一种测量舒张期血流速度的方法、装置、血管评定参数的计算方法、冠状动脉分析系统及计算机存储介质,以解决通过无创血压仪获得的舒张期的血流速度不准确的问题。
为实现上述目的,第一方面,本申请提供了一种测量舒张期血流速度的方法,包括:
读取至少一个体位的冠状动脉二维造影图像组;
从所述冠状动脉二维造影图像组中,提取感兴趣的血管段;
提取所述血管段的中心线;
对任意两帧冠状动脉二维造影图像中的造影剂流经所述血管段的时间做差,差值为Δt,且对所述分段中心线做差,差值为ΔL;
根据所述ΔL与所述Δt的比值,求解血流速度;
选取所述血流速度的最大值,即为舒张期的血流速度。
可选地,上述的测量舒张期血流速度的方法,所述读取至少一个体位的冠状动脉二维造影图像组的方法包括:
通过无线或者有线方式从造影图像拍摄装置或者医院平台上,直接读取至少一个体位的冠状动脉二维造影图像组;或
通过存储装置读取至少一个体位的冠状动脉二维造影图像组。
可选地,上述的测量舒张期血流速度的方法,所述从所述冠状动脉二维造影图像组中,提取感兴趣的血管段的方法,包括:
从所述冠状动脉二维造影图像组中选取N帧冠状动脉二维造影图像;
在所述冠状动脉二维造影图像上,拾取感兴趣的所述血管的首末点,获取感兴趣的所述血管段。
可选地,上述的测量舒张期血流速度的方法,所述提取所述血管段的中心线的方法包括:
从所述冠状动脉二维造影图像中提取血管骨架;
依据所述血管段的延伸方向,以及两点之间获取最短路径的原则;
沿着所述血管骨架,提取所述血管段的中心线。
可选地,上述的测量舒张期血流速度的方法,所述沿着所述血管骨架,提取所述血管段的中心线的方法还包括:
在感兴趣的所述血管段上添加至少一个种子点;
根据所述首末点、种子点,沿着所述血管骨架,重新生成所述血管中心线。
可选地,上述的测量舒张期血流速度的方法,所述对任意两帧冠状动脉二维造影图像中的造影剂流经所述血管段的时间做差,差值为Δt,且对所述分段中心线做差,差值为ΔL;根据所述ΔL与所述Δt的比值,求解血流速度的方法,包括:
以造影剂流至冠状动脉入口,即所述血管段的首点时的冠脉造影图像作为第一帧图像,以造影剂流至所述血管段的末点时的冠脉造影图像作为第N帧图像;
依次求解第N帧图像至第N-1帧,...,第N-b帧,...,第N-a帧,...,第1帧图像的时间差和中心线长度差,时间差分别为Δt1,...,Δtb,…,Δta,…,ΔtN-1;中心线长度差分别为ΔL1,…,ΔLb,...,ΔLa,...,ΔLN-1;
Δt=m×fps,由于每组冠状动脉二维造影图像组中含有连续播放的多帧冠状动脉二维造影图像,因此m表示每组冠状动脉二维造影图像组中,选取的两帧冠状动脉二维造影图像所处帧数的差值,fps表示相邻两帧图像之间切换的间隔时间,优选地,fps=1/15秒。
根据v=ΔL/Δt,其中,v表示血流速度,分别得到第N帧图像至第N-1帧,...,第N-b帧,...,第N-a帧,...,第1帧图像的血流速度,血流速度分别为v1,...,vb,...,va,...,vN-1。
可选地,上述的测量舒张期血流速度的方法,所述对任意两帧冠状动脉二维造影图像中的造影剂流经所述血管段的时间做差,差值为Δt,且对所述分段中心线做差,差值为ΔL;根据所述ΔL与所述Δt的比值,求解血流速度的方法,包括:
依次求解第N帧至第b帧,第N-1帧至第b-1帧,...,第N-b-a帧至第N-a帧,...,第N-b+1至第1帧的图像的时间差和中心线长度差;
根据v=ΔL/Δt,其中,v表示血流速度,分别得到第N帧至第b帧,第N-1帧至第b-1帧,...,第N-b-a帧至第N-a帧,...,第N-b+1至第1帧图像的血流速度。
可选地,上述的测量舒张期血流速度的方法,所述选取所述血流速度的最大值,即为舒张期的血流速度的方法,包括:
通过递归算法或者冒泡算法,选取所述血流速度的最大值,即为舒张期的血流速度;或
通过递归算法或者冒泡算法,选取所述血流速度的最大值,即为舒张期的血流速度;选取所述血流速度的最小值,即为收缩期的血流速度。
可选地,上述的测量舒张期血流速度的方法,在所述提取所述血管段的中心线的方法之后,在所述对任意两帧冠状动脉二维造影图像中的造影剂流经所述血管段的时间做差,差值为Δt,且对所述分段中心线做差,差值为ΔL的方法之前,还包括:
读取至少两个体位的冠状动脉二维造影图像组;
获取所述血管段的几何结构信息;
对感兴趣的所述血管段进行图形处理;
提取所述血管段的血管轮廓线;
根据所述血管段的几何结构信息,将至少两个体位的提取了血管的所述中心线、轮廓线的冠状动脉二维造影图像投影在三维平面上,合成三维血管模型。
可选地,上述的测量舒张期血流速度的方法,所述根据所述ΔL与所述Δt的比值,求解血流速度的方法包括:
对任意两帧冠状动脉二维造影图像中的造影剂流经血管段的时间做差,差值为Δt;根据三维血管模型,获取三维血管模型的中心线,对通过冠状动脉二维造影图像提取的中心线进行修正,且对修正后的分段中心线做差,差值为ΔL’;
根据ΔL’与Δt的比值,求解血流速度v。
第二方面,本申请提供了一种血管评定参数的计算方法,包括:上述任一项所述的测量舒张期血流速度的方法。
第三方面,本申请提供了一种测量舒张期血流速度的装置,用于上述任一项所述的测量舒张期血流速度的方法,包括:依次连接的图像读取单元、血管段提取单元、中心线提取单元,与所述图像读取单元连接的时间差单元,分别与所述时间差单元、所述中心线差单元连接的血流速度获取单元;所述中心线差单元与所述中心线提取单元连接;
所述图像读取单元,用于读取至少一个体位的冠状动脉二维造影图像组;
所述血管段提取单元,用于接收所述图像读取单元发送的冠状动脉二维造影图像,提取所述图像中感兴趣的血管段;
所述中心线提取单元,用于接收所述血管段提取单元发送的血管段,提取所述血管段的中心线;
所述时间差单元,用于接收所述图像读取单元发送的任意两帧冠状动脉二维造影图像,对所述两帧冠状动脉二维造影图像中的造影剂流经所述血管段的时间做差,差值为Δt;
所述中心线差单元,用于接收所述中心线提取单元发送的两帧冠状动脉二维造影图像中的造影剂流经血管段的分段中心线,对所述分段中心线做差,差值为ΔL;
所述血流速度获取单元,包括血流速度计算模块、舒张期血流速度计算模块,所述血流速度计算模块分别与所述时间差单元、所述中心线差单元连接,所述舒张期血流速度计算模块与所述血流速度计算模块连接;
血流速度计算模块,用于接收所述时间差单元和所述中心线差单元发送的所述ΔL与所述Δt,根据所述ΔL与所述Δt的比值,求解血流速度;
所述舒张期血流速度计算模块,用于接收所述血流速度计算模块发送的血流速度值,选取所述血流速度的最大值,即为舒张期的血流速度。
可选地,上述的测量舒张期血流速度的装置,还包括:均与所述图像读取单元连接的血管骨架提取单元、几何信息获取单元、三维血管重建单元,与所述血管骨架提取单元连接的轮廓线提取单元,均与所述中心线提取单元连接的几何信息获取单元、所述三维血管重建单元;所述三维血管重建单元与所述几何信息获取单元连接;
所述血管骨架提取单元,用于接收所述图像读取单元发送的冠状动脉二维造影图像,提取所述图像中的血管骨架;
所述轮廓线提取单元,用于接收所述血管骨架提取单元的血管骨架,根据所述血管骨架,提取感兴趣的所述血管段的轮廓线;
所述几何信息获取单元,用于接收所述图像读取单元的冠状动脉二维造影图像,接收所述中心线提取单元的中心线,接收所述轮廓线提取单元的轮廓线,获取所述血管段的几何结构信息;
所述三维血管重建单元,用于接收所述轮廓线提取单元、所述几何信息获取单元、所述中心线提取单元发送的轮廓线、几何结构信息和中心线,接收所述图像读取单元发送的冠状动脉二维造影图像,用于根据所述血管段的几何结构信息,将至少两个体位的提取了血管的所述中心线、轮廓线的冠状动脉二维造影图像投影在三维平面上,合成三维血管模型。
第四方面,本申请提供了一种冠状动脉分析系统,包括:上述的测量舒张期血流速度的装置。
第五方面,本申请提供了一种计算机存储介质,计算机程序被处理器执行时实现上述任一项所述的测量舒张期血流速度的方法。
本申请实施例提供的方案带来的有益效果至少包括:
本申请提供了一种测量舒张期血流速度的方法,通过对任意两帧冠状动脉二维造影图像中的造影剂流经所述血管段的时间做差,差值为Δt,且对所述分段中心线做差,差值为ΔL;根据所述ΔL与所述Δt的比值,求解血流速度;选取所述血流速度的最大值,即为舒张期的血流速度;无需通过无创血压仪测量,如此测量的舒张期的血流速度不会受到情绪等外界因素的影响,提高了测量舒张期血流速度的准确度。
附图说明
此处所说明的附图用来提供对本发明的进一步理解,构成本发明的一部分,本发明的示意性实施例及其说明用于解释本发明,并不构成对本发明的不当限定。在附图中:
图1为本申请的测量舒张期血流速度的方法的实施例1的流程图;
图2为本申请的S200的流程图;
图3为本申请的S300的流程图;
图4为本申请的S330的流程图;
图5为本申请的S400的一个实施例的流程图;
图6为本申请的S400的另一实施例的流程图;
图7为本申请的测量舒张期血流速度的方法的实施例3的流程图;
图8为本申请的S700的流程图;
图9为本申请的测量舒张期血流速度的装置的一个实施例的结构框图;
图10为申请的测量舒张期血流速度的装置的另一实施例的结构框图;
下面对附图标记进行说明:
图像读取单元1,血管段提取单元2,中心线提取单元3,时间差单元4,中心线差单元5,血流速度获取单元6,血流速度计算模块610,舒张期血流速度计算模块620,血管骨架提取单元7,几何信息获取单元8,三维血管重建单元9,轮廓线提取单元10。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面将结合本发明具体实施例及相应的附图对本发明技术方案进行清楚、完整地描述。显然,所描述的实施例仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
以下将以图式揭露本发明的多个实施方式,为明确说明起见,许多实务上的细节将在以下叙述中一并说明。然而,应了解到,这些实务上的细节不应用以限制本发明。也就是说,在本发明的部分实施方式中,这些实务上的细节是非必要的。此外,为简化图式起见,一些习知惯用的结构与组件在图式中将以简单的示意的方式绘示之。
血管评定参数包括:冠状动脉舒张期的血流速度IFR,以及冠状动脉舒张期的微循环阻力指数IFMR等;而IFR和IFMR均需要基于冠状动脉舒张期的血流速度,目前舒张期的血流速度均是通过无创血压仪获得的,存在舒张期的血流速度测量不准确的问题。
实施例1:
为了解决上述问题,如图1所示,本申请提供了一种测量舒张期血流速度的方法,包括:
S100,读取至少一个体位的冠状动脉二维造影图像组;
S200,从冠状动脉二维造影图像组中,提取感兴趣的血管段;
S300,提取血管段的中心线;
S400,对任意两帧冠状动脉二维造影图像中的造影剂流经血管段的时间做差,差值为Δt,且对分段中心线做差,差值为ΔL;
S500,根据ΔL与Δt的比值,求解血流速度;
S600,选取血流速度的最大值,即为舒张期的血流速度。
本申请提供了一种测量舒张期血流速度的方法,通过对任意两帧冠状动脉二维造影图像中的造影剂流经所述血管段的时间做差,差值为Δt,且对所述分段中心线做差,差值为ΔL;根据所述ΔL与所述Δt的比值,求解血流速度;选取所述血流速度的最大值,即为舒张期的血流速度;无需通过无创血压仪测量,如此测量的舒张期的血流速度不会受到情绪等外界因素的影响,提高了测量舒张期血流速度的准确度。
实施例2:
如图1所示,本申请提供了一种测量舒张期血流速度的方法,包括:
S100,读取至少一个体位的冠状动脉二维造影图像组,具体包括:
通过无线或者有线方式从造影图像拍摄装置或者医院平台上,直接读取至少一个体位的冠状动脉二维造影图像组;或通过存储装置读取至少一个体位的冠状动脉二维造影图像组;
S200,从冠状动脉二维造影图像组中,提取感兴趣的血管段,具体如图2所示,包括:
S210,从冠状动脉二维造影图像组中选取N帧冠状动脉二维造影图像;
S220,在冠状动脉二维造影图像上,拾取感兴趣的血管的首末点,获取感兴趣的血管段;
S300,提取血管段的中心线,具体如图3所示,包括:
S310,从冠状动脉二维造影图像中提取血管骨架;
S320,依据血管段的延伸方向,以及两点之间获取最短路径的原则;
S330,沿着血管骨架,提取血管段的中心线,具体如图4所示,包括:
S331,在感兴趣的血管段上添加至少一个种子点;
S332,根据首末点、种子点,沿着血管骨架,重新生成血管中心线;
S400,对任意两帧冠状动脉二维造影图像中的造影剂流经血管段的时间做差,差值为Δt,且对分段中心线做差,差值为ΔL;
S500,根据ΔL与Δt的比值,求解血流速度;
S600,选取血流速度的最大值,即为舒张期的血流速度;优选地,本申请通过递归算法或者冒泡算法,选取血流速度的最大值,即为舒张期的血流速度。
本申请的一个实施例中,S200还包括:将有导管出现的第一帧冠状动脉二维造影图像定义为参考图像,将有完整冠状动脉出现的第k帧冠状动脉二维造影图像定义为目标图像,k为大于1的正整数;将参考图像减去目标图像,提取导管的特征点O;优选地,除去部分静态噪声;进一步地,采用均值滤波,除去部分动态噪声;以及通过灰度直方图分析,利用阈值进一步去噪;将目标图像减去的参考图像,提取冠状动脉所处位置的区域图像;区域图像以导管的特征点作为种子点进行动态生长,获得感兴趣的血管段图像。
本申请的一个实施例中,S400包括两种获取方法,A方法如图5所示,包括:
S410A,以造影剂流至冠状动脉入口,即血管段的首点时的冠脉造影图像作为第一帧图像,以造影剂流至血管段的末点时的冠脉造影图像作为第N帧图像;
S420A,依次求解第N帧图像至第N-1帧,...,第N-b帧,...,第N-a帧,...,第1帧图像的时间差和中心线长度差,时间差分别为Δt1,...,Δtb,...,Δta,...,ΔtN-1;中心线长度差分别为ΔL1,...,ΔLb,...,ΔLa,...,ΔLN-1;
Δt=m×fps,由于每组冠状动脉二维造影图像组中含有连续播放的多帧冠状动脉二维造影图像,因此m表示每组冠状动脉二维造影图像组中,选取的两帧冠状动脉二维造影图像所处帧数的差值,fps表示相邻两帧图像之间切换的间隔时间,优选地,fps=1/15秒;
S430A,根据v=ΔL/Δt,其中,v表示血流速度,分别得到第N帧图像至第N-1帧,...,第N-b帧,...,第N-a帧,...,第1帧图像的血流速度,血流速度分别为v1,...,vb,...,va,...,vN-1。
本申请的一个实施例中,S400包括两种获取方法,B方法如图6所示,包括:
S410B,以造影剂流至冠状动脉入口,即血管段的首点时的冠脉造影图像作为第一帧图像,以造影剂流至血管段的末点时的冠脉造影图像作为第N帧图像;
S420B,依次求解第N帧至第b帧,第N-1帧至第b-1帧,...,第N-b-a帧至第N-a帧,...,第N-b+1至第1帧的图像的时间差和中心线长度差;
Δt=m×fps,由于每组冠状动脉二维造影图像组中含有连续播放的多帧冠状动脉二维造影图像,因此m表示每组冠状动脉二维造影图像组中,选取的两帧冠状动脉二维造影图像所处帧数的差值,fps表示相邻两帧图像之间切换的间隔时间,优选地,fps=1/15秒;
S430B,根据v=ΔL/Δt,其中,v表示血流速度,分别得到第N帧至第b帧,第N-1帧至第b-1帧,...,第N-b-a帧至第N-a帧,...,第N-b+1至第1帧图像的血流速度。
本申请也可以通过递归算法或者冒泡算法,选取血流速度的最小值,即为收缩期的血流速度。
实施例3:
如图7所示,本申请通过三维建模获取血管段中心线长度和舒张期的血流速度,包括:
S100,读取至少两个体位的冠状动脉二维造影图像组;
S200,从冠状动脉二维造影图像组中,提取感兴趣的血管段;
S300,提取血管段的中心线;
S700,将至少两个体位的冠状动脉二维造影图像投影在三维平面上,合成三维血管模型;
S400,对任意两帧冠状动脉二维造影图像中的造影剂流经血管段的时间做差,差值为Δt;根据三维血管模型,获取三维血管模型的中心线,对通过冠状动脉二维造影图像提取的中心线进行修正,且对修正后的分段中心线做差,差值为ΔL’;根据ΔL’与Δt的比值,求解血流速度v。
S500,根据ΔL’与Δt的比值,求解血流速度;
S600,选取血流速度的最大值,即为舒张期的血流速度。
如图8所示,本申请的一个实施例中,S700包括:
S710,获取血管段的几何结构信息;
S720,对感兴趣的血管段进行图形处理;
S730,提取血管段的血管轮廓线;
S740,根据血管段的几何结构信息,将至少两个体位的提取了血管的中心线、轮廓线的冠状动脉二维造影图像投影在三维平面上,合成三维血管模型。
实施例4:
本申请提供了一种血管评定参数的计算方法,包括:上述任一项的测量舒张期血流速度的方法。
血管评定参数包括:IMR、IFR等。
实施例5:
如图9所示,本申请提供了一种测量舒张期血流速度的装置,用于上述任一项的测量舒张期血流速度的方法,包括:依次连接的图像读取单元1、血管段提取单元2、中心线提取单元3,与图像读取单元1连接的时间差单元4,分别与时间差单元4、中心线差单元5连接的血流速度获取单元6,中心线差单元5与中心线提取单元3连接;图像读取单元1,用于读取至少一个体位的冠状动脉二维造影图像组;血管段提取单元2,用于接收图像读取单元1发送的冠状动脉二维造影图像,提取图像中感兴趣的血管段;中心线提取单元3,用于接收血管段提取单元2发送的血管段,提取血管段的中心线;时间差单元4,用于接收图像读取单元1发送的任意两帧冠状动脉二维造影图像,对两帧冠状动脉二维造影图像中的造影剂流经血管段的时间做差,差值为Δt;中心线差单元5,用于接收中心线提取单元3发送的两帧冠状动脉二维造影图像中的造影剂流经血管段的分段中心线,对分段中心线做差,差值为ΔL;血流速度获取单元6,包括血流速度计算模块610、舒张期血流速度计算模块620,血流速度计算模块610分别与时间差单元4、中心线差单元5连接,舒张期血流速度计算模块620与血流速度计算模块610连接;血流速度计算模块610,用于接收时间差单元4和中心线差单元5发送的ΔL与Δt,根据ΔL与Δt的比值,求解血流速度v;舒张期血流速度计算模块620,用于接收血流速度计算模块610发送的血流速度v值,选取血流速度的最大值vmax,即为舒张期的血流速度。
如图10所示,本申请的一个实施例中,还包括:均与图像读取单元1连接的血管骨架提取单元7、几何信息获取单元8、三维血管重建单元9,与血管骨架提取单元7连接的轮廓线提取单元10,均与中心线提取单元3连接的几何信息获取单元8、三维血管重建单元9;三维血管重建单元9与几何信息获取单元8连接;血管骨架提取单元7,用于接收图像读取单元1发送的冠状动脉二维造影图像,提取图像中的血管骨架;轮廓线提取单元10,用于接收血管骨架提取单元7的血管骨架,根据血管骨架,提取感兴趣的血管段的轮廓线;几何信息获取单元8,用于接收图像读取单元1的冠状动脉二维造影图像,接收中心线提取单元3的中心线,接收轮廓线提取单元10的轮廓线,获取血管段的几何结构信息;三维血管重建单元9,用于接收轮廓线提取单元10、几何信息获取单元8、中心线提取单元3发送的轮廓线、几何结构信息和中心线,接收图像读取单元1发送的冠状动脉二维造影图像,用于根据血管段的几何结构信息(包括体位拍摄角度、患者性别、年龄等),将至少两个体位的提取了血管的中心线以及中心线长度、轮廓线的冠状动脉二维造影图像投影在三维平面上,合成三维血管模型。
本申请提供了一种冠状动脉分析系统,包括:上述的测量舒张期血流速度的装置。
本申请提供了一种计算机存储介质,计算机程序被处理器执行时实现上述任一项的测量舒张期血流速度的方法。
所属技术领域的技术人员知道,本发明的各个方面可以实现为系统、方法或计算机程序产品。因此,本发明的各个方面可以具体实现为以下形式,即:完全的硬件实施方式、完全的软件实施方式(包括固件、驻留软件、微代码等),或硬件和软件方面结合的实施方式,这里可以统称为“电路”、“模块”或“系统”。此外,在一些实施例中,本发明的各个方面还可以实现为在一个或多个计算机可读介质中的计算机程序产品的形式,该计算机可读介质中包含计算机可读的程序代码。本发明的实施例的方法和/或系统的实施方式可以涉及到手动地、自动地或以其组合的方式执行或完成所选任务。
例如,可以将用于执行根据本发明的实施例的所选任务的硬件实现为芯片或电路。作为软件,可以将根据本发明的实施例的所选任务实现为由计算机使用任何适当操作系统执行的多个软件指令。在本发明的示例性实施例中,由数据处理器来执行如本文的根据方法和/或系统的示例性实施例的一个或多个任务,诸如用于执行多个指令的计算平台。可选地,该数据处理器包括用于存储指令和/或数据的易失性储存器和/或用于存储指令和/或数据的非易失性储存器,例如,磁硬盘和/或可移动介质。可选地,也提供了一种网络连接。可选地也提供显示器和/或用户输入设备,诸如键盘或鼠标。
可利用一个或多个计算机可读的任何组合。计算机可读介质可以是计算机可读信号介质或计算机可读存储介质。计算机可读存储介质例如可以是——但不限于——电、磁、光、电磁、红外线、或半导体的系统、装置或器件,或者任意以上的组合。计算机可读存储介质的更具体的例子(非穷举列表)将包括以下各项:
具有一个或多个导线的电连接、便携式计算机盘、硬盘、随机存取存储器(RAM)、只读存储器(ROM)、可擦式可编程只读存储器(EPROM或闪存)、光纤、便携式紧凑盘只读存储器(CD-ROM)、光存储器件、磁存储器件、或者上述的任意合适的组合。在本文件中,计算机可读存储介质可以是任何包含或存储程序的有形介质,该程序可以被指令执行系统、装置或者器件使用或者与其结合使用。
计算机可读的信号介质可以包括在基带中或者作为载波一部分传播的数据信号,其中承载了计算机可读的程序代码。这种传播的数据信号可以采用多种形式,包括但不限于电磁信号、光信号或上述的任意合适的组合。计算机可读的信号介质还可以是计算机可读存储介质以外的任何计算机可读介质,该计算机可读介质可以发送、传播或者传输用于由指令执行系统、装置或者器件使用或者与其结合使用的程序。
计算机可读介质上包含的程序代码可以用任何适当的介质传输,包括(但不限于)无线、有线、光缆、RF等等,或者上述的任意合适的组合。
例如,可用一个或多个编程语言的任何组合来编写用于执行用于本发明的各方面的操作的计算机程序代码,包括诸如Java、Smalltalk、C++等面向对象编程语言和常规过程编程语言,诸如"C"编程语言或类似编程语言。程序代码可以完全地在用户计算机上执行、部分地在用户计算机上执行、作为一个独立的软件包执行、部分在用户计算机上部分在远程计算机上执行、或者完全在远程计算机或服务器上执行。在涉及远程计算机的情形中,远程计算机可以通过任意种类的网络--包括局域网(LAN)或广域网(WAN)-连接到用户计算机,或者,可以连接到外部计算机(例如利用因特网服务提供商来通过因特网连接)。
应当理解,流程图和/或框图的每个方框以及流程图和/或框图中各方框的组合,都可以由计算机程序指令实现。这些计算机程序指令可以提供给通用计算机、专用计算机或其它可编程数据处理装置的处理器,从而生产出一种机器,使得这些计算机程序指令在通过计算机或其它可编程数据处理装置的处理器执行时,产生了实现流程图和/或框图中的一个或多个方框中规定的功能/动作的装置。
也可以把这些计算机程序指令存储在计算机可读介质中,这些指令使得计算机、其它可编程数据处理装置、或其它设备以特定方式工作,从而,存储在计算机可读介质中的指令就产生出包括实现流程图和/或框图中的一个或多个方框中规定的功能/动作的指令的制造品(article of manufacture)。
还可将计算机程序指令加载到计算机(例如,冠状动脉分析系统)或其它可编程数据处理设备上以促使在计算机、其它可编程数据处理设备或其它设备上执行一系列操作步骤以产生计算机实现过程,使得在计算机、其它可编程装置或其它设备上执行的指令提供用于实现在流程图和/或一个或多个框图方框中指定的功能/动作的过程。
本发明的以上的具体实例,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上仅为本发明的具体实施例而已,并不用于限制本发明,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (11)
1.一种测量舒张期血流速度的方法,其特征在于,包括:
读取至少一个体位的冠状动脉二维造影图像组;
从所述冠状动脉二维造影图像组中,提取感兴趣的血管段;
从所述冠状动脉二维造影图像中提取血管骨架;
依据所述血管段的延伸方向,以及两点之间获取最短路径的原则;
沿着所述血管骨架,提取所述血管段的中心线;
根据血管段的几何结构信息,将至少两个体位的提取了血管的所述中心线、轮廓线的冠状动脉二维造影图像投影在三维平面上,合成三维血管模型;
对任意两帧冠状动脉二维造影图像中的造影剂流经所述血管段的时间做差,差值为Δt,Δt=m×fps,由于每组冠状动脉二维造影图像组中含有连续播放的多帧冠状动脉二维造影图像,因此m表示每组冠状动脉二维造影图像组中,选取的两帧冠状动脉二维造影图像所处帧数的差值,fps表示相邻两帧图像之间切换的间隔时间,包括:以造影剂流至冠状动脉入口,即所述血管段的首点时的冠脉造影图像作为第一帧图像,以造影剂流至所述血管段的末点时的冠脉造影图像作为第N帧图像;依次求解第N帧图像至第N-1帧,...,第N-b帧,...,第N-a帧,...,第1帧图像的时间差,时间差分别为Δt1,...,Δtb,...,Δta,...,ΔtN-1;或依次求解第N帧至第b帧,第N-1帧至第b-1帧,...,第N-b-a帧至第N-a帧,...,第N-b+1至第1帧的图像的时间差;
对分段的所述中心线做差,差值为ΔL,包括:依次求解第N帧图像至第N-1帧,...,第N-b帧,...,第N-a帧,...,第1帧图像的中心线长度差,分别为ΔL1,...,ΔLb,...,ΔLa,...,ΔLN-1;或依次求解第N帧至第b帧,第N-1帧至第b-1帧,...,第N-b-a帧至第N-a帧,...,第N-b+1至第1帧的图像的中心线长度差;
根据三维血管模型,获取三维血管模型的中心线,对通过冠状动脉二维造影图像提取的中心线进行修正,且对修正后的分段中心线做差,差值为ΔL’;
根据ΔL’与Δt的比值,求解血流速度v;其中,v表示血流速度,分别得到第N帧图像至第N-1帧,...,第N-b帧,...,第N-a帧,...,第1帧图像的血流速度,血流速度分别为v1,...,vb,...,va,...,vN-1;或分别得到第N帧至第b帧,第N-1帧至第b-1帧,...,第N-b-a帧至第N-a帧,...,第N-b+1至第1帧图像的血流速度;
选取所述血流速度的最大值,即为舒张期的血流速度。
2.根据权利要求1所述的测量舒张期血流速度的方法,其特征在于,所述读取至少一个体位的冠状动脉二维造影图像组的方法包括:
通过无线或者有线方式从造影图像拍摄装置或者医院平台上,直接读取至少一个体位的冠状动脉二维造影图像组;或
通过存储装置读取至少一个体位的冠状动脉二维造影图像组。
3.根据权利要求1所述的测量舒张期血流速度的方法,其特征在于,所述从所述冠状动脉二维造影图像组中,提取感兴趣的血管段的方法,包括:
从所述冠状动脉二维造影图像组中选取N帧冠状动脉二维造影图像;
在所述冠状动脉二维造影图像上,拾取感兴趣的所述血管的首末点,获取感兴趣的所述血管段。
4.根据权利要求3所述的测量舒张期血流速度的方法,其特征在于,所述沿着所述血管骨架,提取所述血管段的中心线的方法还包括:
在感兴趣的所述血管段上添加至少一个种子点;
根据所述首末点、种子点,沿着所述血管骨架,重新生成血管中心线。
5.根据权利要求1所述的测量舒张期血流速度的方法,其特征在于,所述选取所述血流速度的最大值,即为舒张期的血流速度的方法,包括:
通过递归算法或者冒泡算法,从权利要求1中选取所述血流速度的最大值,即为舒张期的血流速度;或
通过递归算法或者冒泡算法,从权利要求1中选取所述血流速度的最大值,即为舒张期的血流速度;选取所述血流速度的最小值,即为收缩期的血流速度。
6.根据权利要求1所述的测量舒张期血流速度的方法,其特征在于,在所述提取所述血管段的中心线的方法之后,在所述对任意两帧冠状动脉二维造影图像中的造影剂流经所述血管段的时间做差,差值为Δt,且对所述分段中心线做差,差值为ΔL的方法之前,还包括:
读取至少两个体位的冠状动脉二维造影图像组;
获取所述血管段的几何结构信息;
对感兴趣的所述血管段进行图形处理;
提取所述血管段的血管轮廓线。
7.一种血管评定参数的计算方法,其特征在于,包括:权利要求1~6任一项所述的测量舒张期血流速度的方法。
8.一种测量舒张期血流速度的装置,用于权利要求1~6任一项所述的测量舒张期血流速度的方法,其特征在于,包括:依次连接的图像读取单元、血管段提取单元、中心线提取单元,与所述图像读取单元连接的时间差单元,分别与所述时间差单元、所述中心线差单元连接的血流速度获取单元;所述中心线差单元与所述中心线提取单元连接;
所述图像读取单元,用于读取至少一个体位的冠状动脉二维造影图像组;
所述血管段提取单元,用于接收所述图像读取单元发送的冠状动脉二维造影图像,提取所述图像中感兴趣的血管段;
所述中心线提取单元,用于接收所述血管段提取单元发送的血管段,提取所述血管段的中心线;
所述时间差单元,用于接收所述图像读取单元发送的任意两帧冠状动脉二维造影图像,对所述两帧冠状动脉二维造影图像中的造影剂流经所述血管段的时间做差,差值为Δt,包括:以造影剂流至冠状动脉入口,即所述血管段的首点时的冠脉造影图像作为第一帧图像,以造影剂流至所述血管段的末点时的冠脉造影图像作为第N帧图像;依次求解第N帧图像至第N-1帧,...,第N-b帧,...,第N-a帧,...,第1帧图像的时间差,分别为Δt1,...,Δtb,...,Δta,...,ΔtN-1;或依次求解第N帧至第b帧,第N-1帧至第b-1帧,...,第N-b-a帧至第N-a帧,...,第N-b+1至第1帧的图像的时间差;
所述中心线差单元,用于接收所述中心线提取单元发送的两帧冠状动脉二维造影图像中的造影剂流经血管段的分段中心线,对所述分段中心线做差,差值为ΔL,包括:依次求解第N帧图像至第N-1帧,...,第N-b帧,...,第N-a帧,...,第1帧图像的时间差和中心线长度差,分别为ΔL1,...,ΔLb,...,ΔLa,...,ΔLN-1,或依次求解第N帧至第b帧,第N-1帧至第b-1帧,...,第N-b-a帧至第N-a帧,...,第N-b+1至第1帧的图像的中心线长度差;
所述血流速度获取单元,包括血流速度计算模块、舒张期血流速度计算模块,所述血流速度计算模块分别与所述时间差单元、所述中心线差单元连接,所述舒张期血流速度计算模块与所述血流速度计算模块连接;
血流速度计算模块,用于接收所述时间差单元和所述中心线差单元发送的所述ΔL与所述Δt,根据所述ΔL与所述Δt的比值,求解血流速度,包括:根据v=ΔL/Δt,其中,v表示血流速度,分别得到第N帧图像至第N-1帧,...,第N-b帧,...,第N-a帧,...,第1帧图像的血流速度,血流速度分别为v1,...,vb,...,va,...,vN-1;或分别得到第N帧至第b帧、第N-1帧至第b-1帧,...,第N-b-a帧至第N-a帧,...,第N-b+1至第1帧图像的血流速度;
所述舒张期血流速度计算模块,用于接收所述血流速度计算模块发送的血流速度值,选取所述血流速度的最大值,即为舒张期的血流速度。
9.根据权利要求8所述的测量舒张期血流速度的装置,其特征在于,还包括:均与所述图像读取单元连接的血管骨架提取单元、几何信息获取单元、三维血管重建单元,与所述血管骨架提取单元连接的轮廓线提取单元,均与所述中心线提取单元连接的几何信息获取单元、所述三维血管重建单元;所述三维血管重建单元与所述几何信息获取单元连接;
所述血管骨架提取单元,用于接收所述图像读取单元发送的冠状动脉二维造影图像,提取所述图像中的血管骨架;
所述轮廓线提取单元,用于接收所述血管骨架提取单元的血管骨架,根据所述血管骨架,提取感兴趣的所述血管段的轮廓线;
所述几何信息获取单元,用于接收所述图像读取单元的冠状动脉二维造影图像,接收所述中心线提取单元的中心线,接收所述轮廓线提取单元的轮廓线,获取所述血管段的几何结构信息;
所述三维血管重建单元,用于接收所述轮廓线提取单元、所述几何信息获取单元、所述中心线提取单元发送的轮廓线、几何结构信息和中心线,接收所述图像读取单元发送的冠状动脉二维造影图像,用于根据所述血管段的几何结构信息,将至少两个体位的提取了血管的所述中心线、轮廓线的冠状动脉二维造影图像投影在三维平面上,合成三维血管模型。
10.一种冠状动脉分析系统,其特征在于,包括:权利要求8或9所述的测量舒张期血流速度的装置。
11.一种计算机存储介质,其特征在于,计算机程序被处理器执行时实现权利要求1~6任一项所述的测量舒张期血流速度的方法。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911066599.8A CN110786842B (zh) | 2019-11-04 | 2019-11-04 | 测量舒张期血流速度的方法、装置、系统及存储介质 |
PCT/CN2019/116637 WO2021087961A1 (zh) | 2019-11-04 | 2019-11-08 | 测量舒张期血流速度的方法、装置、系统及存储介质 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911066599.8A CN110786842B (zh) | 2019-11-04 | 2019-11-04 | 测量舒张期血流速度的方法、装置、系统及存储介质 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110786842A CN110786842A (zh) | 2020-02-14 |
CN110786842B true CN110786842B (zh) | 2022-10-04 |
Family
ID=69442600
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911066599.8A Active CN110786842B (zh) | 2019-11-04 | 2019-11-04 | 测量舒张期血流速度的方法、装置、系统及存储介质 |
Country Status (2)
Country | Link |
---|---|
CN (1) | CN110786842B (zh) |
WO (1) | WO2021087961A1 (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112785580B (zh) * | 2021-01-28 | 2024-02-02 | 深圳睿心智能医疗科技有限公司 | 确定血管流速的方法及装置 |
CN113269806B (zh) * | 2021-04-30 | 2023-04-14 | 北京阅影科技有限公司 | 测量血管内部血流流量的方法、装置与处理器 |
CN115337096A (zh) * | 2021-05-13 | 2022-11-15 | 苏州润迈德医疗科技有限公司 | 用于内脏动脉去交感神经的消融方法、系统及存储介质 |
CN116035549A (zh) * | 2022-11-23 | 2023-05-02 | 杭州脉流科技有限公司 | 获得冠脉血流储备分数的方法和计算机设备 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108186038A (zh) * | 2018-02-11 | 2018-06-22 | 杭州脉流科技有限公司 | 基于动脉造影影像计算冠脉血流储备分数的系统 |
Family Cites Families (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN2792415Y (zh) * | 2005-04-25 | 2006-07-05 | 中山医科大学生物医学工程开发中心 | 血管保健康复治疗仪 |
US20180344174A9 (en) * | 2009-09-23 | 2018-12-06 | Lightlab Imaging, Inc. | Lumen Morphology and Vascular Resistance Measurements Data Collection Systems, Apparatus and Methods |
CN103110420A (zh) * | 2013-01-31 | 2013-05-22 | 深圳先进技术研究院 | 磁共振血管成像方法和系统 |
CN106327487B (zh) * | 2016-08-18 | 2018-01-02 | 苏州润迈德医疗科技有限公司 | 基于x射线冠脉造影图像的冠状动脉血流储备分数计算方法 |
EP3403582B1 (en) * | 2017-05-15 | 2023-06-07 | Pie Medical Imaging BV | Method and apparatus for determining blood velocity in x-ray angiography images |
CN108245178A (zh) * | 2018-01-11 | 2018-07-06 | 苏州润迈德医疗科技有限公司 | 一种基于x射线冠脉造影图像的血液流动速度计算方法 |
CN108133478B (zh) * | 2018-01-11 | 2021-11-23 | 苏州润迈德医疗科技有限公司 | 一种提取冠状动脉血管中心线的方法 |
CN110226923B (zh) * | 2018-03-05 | 2021-12-14 | 苏州润迈德医疗科技有限公司 | 一种无需血管扩张剂测量血流储备分数的方法 |
CN108550189A (zh) * | 2018-05-03 | 2018-09-18 | 苏州润迈德医疗科技有限公司 | 基于造影图像和流体力学模型的微循环阻力指数计算方法 |
CN110384494A (zh) * | 2018-09-19 | 2019-10-29 | 苏州润迈德医疗科技有限公司 | 测量微循环阻力指数的方法 |
CN109770888A (zh) * | 2019-03-19 | 2019-05-21 | 苏州润迈德医疗科技有限公司 | 基于压力传感器和造影图像计算瞬时无波形比率的方法 |
-
2019
- 2019-11-04 CN CN201911066599.8A patent/CN110786842B/zh active Active
- 2019-11-08 WO PCT/CN2019/116637 patent/WO2021087961A1/zh active Application Filing
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108186038A (zh) * | 2018-02-11 | 2018-06-22 | 杭州脉流科技有限公司 | 基于动脉造影影像计算冠脉血流储备分数的系统 |
Non-Patent Citations (1)
Title |
---|
基于图像处理的冠脉造影图像二维分析研究;小卧;《道客巴巴》;20140222;摘要,第28-37页 * |
Also Published As
Publication number | Publication date |
---|---|
WO2021087961A1 (zh) | 2021-05-14 |
CN110786842A (zh) | 2020-02-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110786842B (zh) | 测量舒张期血流速度的方法、装置、系统及存储介质 | |
US10803994B2 (en) | Vascular flow assessment | |
US11779294B2 (en) | Method, device and system for calculating microcirculation indicator based on image and pressure sensor | |
CN108038848B (zh) | 基于医学影像序列斑块稳定性指标的快速计算方法及系统 | |
CN110889896B (zh) | 获取血管狭窄病变区间及三维合成方法、装置和系统 | |
JP6745875B2 (ja) | 圧力又は流量測定及び血管造影法からの流量、抵抗又は圧力の推定 | |
CN109805949B (zh) | 基于压力传感器和造影图像计算血流储备分数的方法 | |
CN110786841B (zh) | 基于微循环阻力指数调节最大充血状态流速的方法及装置 | |
CN112384137B (zh) | 基于造影图像获取静息态下血管评定参数的方法及装置 | |
CN108665449B (zh) | 预测血流矢量路径上的血流特征的深度学习模型及装置 | |
CN112164020A (zh) | 精确提取血管中心线的方法、装置、分析系统和存储介质 | |
CN110786840B (zh) | 基于生理参数获取血管评定参数的方法、装置及存储介质 | |
CN112155580B (zh) | 基于造影图像修正血流速度和微循环参数的方法及装置 | |
CN112384138B (zh) | 获取心表大动脉的血流量的方法、装置、系统及存储介质 | |
CN110929604B (zh) | 基于造影图像的流速的筛选方法、装置、系统和存储介质 |
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 |