CN112837306A - 基于深度学习和中智理论的冠状动脉病变功能学定量方法 - Google Patents
基于深度学习和中智理论的冠状动脉病变功能学定量方法 Download PDFInfo
- Publication number
- CN112837306A CN112837306A CN202110192837.0A CN202110192837A CN112837306A CN 112837306 A CN112837306 A CN 112837306A CN 202110192837 A CN202110192837 A CN 202110192837A CN 112837306 A CN112837306 A CN 112837306A
- Authority
- CN
- China
- Prior art keywords
- coronary artery
- coronary
- region
- image
- blood flow
- 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
Images
Classifications
-
- 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
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
- G06F18/241—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/004—Artificial life, i.e. computing arrangements simulating life
- G06N3/006—Artificial life, i.e. computing arrangements simulating life based on simulated virtual individual or collective life forms, e.g. social simulations or particle swarm optimisation [PSO]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/08—Learning methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/11—Region-based segmentation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/187—Segmentation; Edge detection involving region growing; involving region merging; involving connected component labelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/194—Segmentation; Edge detection involving foreground-background segmentation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/60—Analysis of geometric attributes
- G06T7/62—Analysis of geometric attributes of area, perimeter, diameter or volume
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/40—Extraction of image or video features
- G06V10/44—Local feature extraction by analysis of parts of the pattern, e.g. by detecting edges, contours, loops, corners, strokes or intersections; Connectivity analysis, e.g. of connected components
-
- 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/10—Image acquisition modality
- G06T2207/10072—Tomographic images
- G06T2207/10081—Computed x-ray tomography [CT]
-
- 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/20—Special algorithmic details
- G06T2207/20081—Training; Learning
-
- 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/20—Special algorithmic details
- G06T2207/20212—Image combination
- G06T2207/20221—Image fusion; Image merging
-
- 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
- G06T2207/30104—Vascular flow; Blood flow; Perfusion
-
- 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/30172—Centreline of tubular or elongated structure
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Data Mining & Analysis (AREA)
- Life Sciences & Earth Sciences (AREA)
- Artificial Intelligence (AREA)
- General Health & Medical Sciences (AREA)
- Health & Medical Sciences (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Physics (AREA)
- Molecular Biology (AREA)
- Software Systems (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- Computational Linguistics (AREA)
- Computing Systems (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Multimedia (AREA)
- Evolutionary Biology (AREA)
- Geometry (AREA)
- Medical Informatics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Quality & Reliability (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
本发明提供了一种基于深度学习和中智理论的冠状动脉病变功能学定量方法,属于生物医学领域,结合CTA和CAG图像的各自特点,首先利用CTA图像来获得冠脉各个心动周期的三维图像数据,然后将其配准并投影到CAG图像的冠脉区域上,根据对比剂的信息确定冠脉狭窄病变区域的位置。该方法克服了单凭CAG二维图像重建血管精度上的不足,能较准确地自动识别冠脉区域,进而根据对比剂的跟踪路线和CAG视频序列的时间,计算出血液在不同冠脉区域的流速,进而计算出冠脉各点的血流量,从而得到病变远端与近端血流量的比值,获得计算机辅助诊断的FFRCAD,实现对冠脉病变功能学的全方位评价,解决冠状动脉病变功能学无创定量测量难题。
Description
技术领域
本发明属于生物医学领域,具体涉及一种基于深度学习和中智理论的冠状动脉病变功能学定量方法。
背景技术
从1990年开始,针对心血管的冠状动脉造影(coronary artery angiography,CAG)图像及其疾病诊断研究逐渐被国内外学者广泛关注。然而,由于心脏特殊运动模式和图像快门速度等因素的限制,使得采集到的图像相邻帧之间区别很大,而且CAG图像灰度对比度较低,冠脉分支互相重叠和交叉,使冠脉的精确检测和测量成为一直困扰国内外研究者的难题。目前针对冠脉检测主要归纳为七大类:
(1)传统模式识别法:是利用传统模式识别手段对血管或血管特征的自动检测和分类,较常用的有血管边缘特征提取法、区域增长法等,但这些方法对CAG图像要求较高,对于复杂冠脉识别较差。
(2)跟踪法:是从初始点开始探测血管中心线或通过分析正交于跟踪方向的像素来检测边缘,是目前应用最广泛、最直观的方法。但由于心脏跳动模式和血管交叉重叠,使跟踪方法的准确性降低。
(3)模型法:是用能够明确描述血管特征的模型或模板来提取血管的方法,该方法缺点是需要较多人工参与,并且距离比较近的血管边缘或中心线会出现干扰,运算量大,耗时长。
(4)专家系统法:将专家系统等人工智能的方法融入到血管的提取方法中。该方法适用于血管灰度与背景灰度差异较大的造影图像,而且在血管交叉和血管阻塞部位得不到满意效果。
(5)神经网络法:Nekovei等提出采用BP神经网络来提取造影图像中的心血管。网络训练采用多幅已经做过血管和背景标记的造影图像。该方法只是在血管和背景像素的区分上作了尝试,不能有效提取血管信息。
(6)波导法:该方法由Caselles和Malladi提出,模仿水波传播特征来提取血管的中心线。提取的中心线连续性比较好,在一定程度上,能反映出血管的拓扑结构,但与血管的整体形状特征和血管走向出入较大,无法满足三维重建要求。
(7)时间域分割法:Uchida等利用DSA图像序列血管的运动信息,对图像序列中冠脉运动小于边缘宽度的连续三幅DSA边缘图像进行相减累加操作得到血管边缘图,进而分割出血管。此种方法只能针对图像背景少,图像目标明显的DSA进行血管分割。这种分割方法过程复杂,对图像序列做多次操作,耗时长,效率低。
近年基于影像学的冠脉狭窄病变无创功能学评价方法成为研究热点,包括基于冠脉计算机断层扫描血管造影(computed tomography angiography,CTA)的FFR检测(coronary CTA–derived FFR,FFRCT)和基于冠脉造影的FFR检测(angiography-derivedFFR)等冠脉功能学评估方法。
FFRCT是基于静息状态心室舒张期冠脉CTA图像的无创冠脉狭窄病变功能学评估方法,不需要特殊扫描技术和额外用药。在心室舒张期冠脉微循环阻力最小,且压力与血流量成正比,通过流体力学计算方法,模拟冠脉内血流和压力,根据冠脉狭窄病变近段与远段血流量和压力阶差,推算管腔横截面,再经过复杂的图像处理和运算过程,包括图像分割和冠脉树提取、左心室质量估算、微循环阻力估算、冠脉最大血流量估算及流体力学计算等,形成管腔横截面和血管轮廓,重建冠脉树的三维解剖模型,从而获取冠脉树任意点的FFRCT值。多项临床研究结果证实:CTA结合无创FFRCT对判断冠脉狭窄病变是否引起心肌缺血的敏感性、特异性及准确性均优于单独应用CTA。
基于冠脉造影的FFR评估(angiography-derived FFR)是应用冠脉造影的快速血流动力学分析方法,结合冠状动脉造影三维重建定量分析(3-dimensional quantitativecoronary angiography,3D QCA)和TIMI计帧法,重建冠状动脉三维结构模型。FFR QCA≤0.80诊断冠脉狭窄病变引起心肌缺血的准确率达88%,整个分析所需时间<10min。FAVORII China Study显示:基于冠脉造影的定量血流分数(quantitative flow ratio,QFR)识别有血流动力学意义冠脉狭窄病变的敏感性达94.6%,特异性达91.7%,较FFR QCA更有实用价值。
冠脉CTA成像获得的是数百层横断面重建的间接冠脉图像,尚不能完全模拟真实的冠脉最大充血状态,判断严重钙化病变和迂曲病变功能学意义的准确性较低,而且模型重建和计算时间长,不能在导管室实时在线使用。另外,CTA只能观察冠脉静态特征,不能提供实时动态冠脉信息。
虽然FFR QCA和QFR提高了冠脉造影的诊断性能,但仍然有其局限性:①获取造影图像时仍需使用腺苷诱导冠脉达最大充血状态;②需要对所有边支血管进行重建,分支血管重叠和交叉对结果影响较大,使用相对复杂;③对分支血管病变功能学判断准确性较低。
针对临床缺乏无创实时评估冠脉狭窄病变功能的影像学方法,本发明的目的是基于深度学习和中智理论,研制开发一套CTA和CAG影像组学的冠脉病变功能学定量分析系统与平台,来解决冠脉病变自动检测和功能学定量评估等关键问题。
发明内容
为了克服上述现有技术存在的不足,本发明提供了一种基于深度学习和中智理论的冠状动脉病变功能学定量方法。
为了实现上述目的,本发明提供如下技术方案:
基于深度学习和中智理论的冠状动脉病变功能学定量方法,包括以下步骤:
步骤1、自动提取CTA中冠脉区域并进行区域处理
步骤1.1、提取冠状动脉CAT图像;
步骤1.2、利用多尺度非局部积分滤波算法对CAT图像进行预处理;
步骤1.3、根据临床经验总结的冠脉知识特征构建多通道信息;
步骤1.4、利用变通道多尺度区域候选分类网络提取CAT图像中的三维冠脉区域;
步骤2、提取CTG图像并将CTA冠脉区域与CTG图像进行投射与融合
步骤2.1、提取冠状动脉CAG图像;
步骤2.2、将CTA中三维冠脉区域通过仿射变换与CAG图像中的对应冠脉进行配准,并投影到二维CAG上,确定冠脉对比剂导航区域;
步骤2.3、基于中智粒子滤波追踪算法追踪对比剂在导航区域内的流动情况;
步骤2.4、自动检测并定位提取CAG冠脉对比剂区域;
步骤2.5、将CTA冠脉区域与CTG图像进行投射与融合;
步骤3、定量评估狭窄病变
步骤3.1、提取融合影像中冠脉的中心线,定位冠脉区域;
步骤3.2、计算冠脉血流速度和载面积,获得局部血流量;
步骤3.3、计算病变远端及近端血流量比值,获得FFRCAD;
步骤3.4、根据FFRCAD值定量评估狭窄病变。
优选地,将所述步骤1.2中多尺度非局部积分滤波算法的函数定义为:
σs(x,y)=f(I(x,y))=aI(x,y)+b (1)
其中,σs是在尺度S下高斯函数核的标准方差,应用其对原始CTA图像进行滤波运算;然后构建冠状动脉知识通道图像如下;
优选地,所述步骤2.2包括以下步骤:
步骤2.2.1、冠脉造影对比剂导航区域确定
利用CAG影像采集时的角度和分辨率,确定仿射变换函数的参数,将CTA中冠脉分割结果进行仿射变换,并将变换后的三维冠脉图像投影到二维CAG图像上,作为冠脉的可能候选区域,也就是对比剂要经过的导航区域;
步骤2.2.2、冠脉造影对比剂区域自动检测和提取:
根据造影对比剂区域在CAG视频的灰度和纹理信息,定义其中智不确定隶属度,从而作为粒子的属性,并定义当前区域及其邻域的中智相似度,来作为粒子滤波的代价函数,利用粒子滤波算法来追踪对比剂在导航区域内的流动情况,从而提取CAG图像中的冠脉区域。
优选地,所述步骤2.3包括以下步骤:
步骤2.3.1、提取冠脉目标特征:通过CTA图像中获得冠脉作为造影的导航区域,人机交互选定待跟踪的初始目标,计算跟踪目标的特征,作为描述目标的特征;
步骤2.3.2、特征搜索阶段:初始化提取了目标的特征后,在目标的附近以高斯分布来撒粒子,并根据导航区域来确定所撒的粒子的范围;每个粒子按照初始化时候得到的目标特征计算出它所处的位置的灰度特征,并计算该直方图向量与目标直方图向量的相似度;然后对每个粒子与目标计算出的相似度做归一化,使得所有粒子的相似度加起来和等于1;
步骤2.3.3、决策阶段:根据每个粒子与目标的相似程度做加权平均;根据最大的相似度候选区域的位置,则为目标的最可能位置;
步骤2.3.4、粒子重采样:根据粒子权重对粒子进行筛选,筛选过程中,既要大量保留权重大的粒子,又要有一小部分权重小的粒子,进而重复特征搜索阶段和决策阶段。
优选地,所述步骤3.1包括以下步骤:
步骤3.1.1、先在冠脉区域上做距离变换,计算各个点到边界的距离;
步骤3.1.2、根据各个点的距离,定义脊点作为冠脉中心线的候选点,并定义各个点之间的距离代价函数;
步骤3.13、利用最短路径Dijkstra算法来通过脊点寻找最短路径,作为冠脉中心线;
步骤3.1.4、在得到冠脉的中心线后,建立冠脉形状模型,根据先验知识构建正常冠脉的形状和直径;
步骤3.1.5、与实际图像计算出的冠脉进行比对,从而检测出冠脉异常区域,定义为冠脉病变区域。
优选地,所述步骤3.2包括以下步骤:
步骤3.2.1、冠脉病变区域血流速度计算
冠脉病变区域的血流速度指血液在单位时间内通过的距离,根据对比剂的跟踪路线和CAG视频序列的时间,可以计算出血液在不同区域的速度(v),从速度上定量分析血液在流经病变区域的变化情况,其计算过程如下:
其中pp和pd是分别在狭窄区域近端和远端的取样点;Disk-1,Disk和Disk+1分别是取样点在第k-1,k和k+1帧上的距离;
步骤3.2.2、冠脉病变区域血流量计算
冠脉病变区域的血流量(Q)是指单位时间内流经病变区截面积的血液体积;利用极坐标系,考虑血管的某一截面积上任一微面积ds上血流速度为v(r,θ,t),则通过该微截面的血流量应为:
(5)式为血流量计算的基本公式,将(4)式变形后得出:
步骤3.2.3、根据冠脉对比剂在冠脉中的横截宽度和CTA的冠脉区域分割结果,利用样条函数来精确拟合出冠脉横截面;
步骤3.2.4、根据获得的最终横截面,可以获得冠脉病变狭窄远端的血流量Qd和病变近端无狭窄处血流量Qp,其计算过程如下:
其中p是取样点,S(p)为p点横截面,V(p)为取样点的血流速度,t为时间间隔;Qp和Qd分别为近端和远端取样点的血流量。
优选地,所述步骤3.3通过计算机辅助获得FFRCAD,其公式为:
FFRCAD=Qd/Qp。 (11)
本发明提供的基于深度学习和中智理论的冠状动脉病变功能学定量方法具有以下有益效果:
本发明的优点是结合CTA和CAG图像的各自特点,首先利用CTA图像来获得冠脉各个心动周期的三维图像数据,然后将其配准并投影到CAG图像的冠脉区域上,根据对比剂的信息确定冠脉狭窄病变区域的位置。该方法克服了单凭CAG二维图像重建血管精度上的不足,能较准确地自动识别冠脉区域,进而根据对比剂的跟踪路线和CAG视频序列的时间,计算出血液在不同冠脉区域的流速,进而计算出冠脉各点的血流量,从而得到病变远端与近端血流量的比值,获得计算机辅助诊断的FFRCAD,实现对冠脉病变功能学的全方位评价,解决冠状动脉病变功能学无创定量测量难题。
附图说明
为了更清楚地说明本发明实施例及其设计方案,下面将对本实施例所需的附图作简单地介绍。下面描述中的附图仅仅是本发明的部分实施例,对于本领域普通技术人员来说,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例1的基于深度学习和中智理论的冠状动脉病变功能学定量方法的流程图;
图2为变通道多尺度区域候选分类网络结构图;
图3为CTA中冠脉RCA分支分割结果三维示意图;
图4为CTA冠脉三维RCA分支和CAG二维图像仿射变换结果示意图;
图5为CAG分割结果;
图6为冠脉分割的三维建模结果。
具体实施方式
为了使本领域技术人员更好的理解本发明的技术方案并能予以实施,下面结合附图和具体实施例对本发明进行详细说明。以下实施例仅用于更加清楚地说明本发明的技术方案,而不能以此来限制本发明的保护范围。
实施例1
本发明提供了一种基于深度学习和中智理论的冠状动脉病变功能学定量方法,具体如图1所示,包括以下步骤:
步骤1、自动提取CTA中冠脉区域并进行区域处理
步骤1.1、提取冠状动脉CAT DICOM图像;
步骤1.2、利用多尺度非局部积分滤波算法对CAT图像进行预处理;
与传统滤波算法相比,多尺度非局部积分滤波算法不仅考虑到当前点和其邻域点特征信息,而且还引入非局部积分的概念,根据邻域点位置特征来确定滤波算子的权值,能更好的保持滤波点和周围邻域点的关系特征,进一步精确冠状动脉特征信息。
多尺度非局部积分滤波算子的核函数可应用高斯函数定义为:
其中,σs是在尺度S下高斯函数核的标准方差;通过σs来确定滤波算子的权值比例。该算法中,将σs定义为与图像不确定性特征的线性函数。如果当前点的不确定性比较高,则σs比较大,需要考虑更广的邻域;如果当前点的不确定性不高,则σs较小,仅需要考虑较小的邻域特征;
σs(x,y)=f(I(x,y))=aI(x,y)+b (1)
应用新定义的多尺度非局部积分滤波算子,对原始CTA图像进行滤波运算。其结果图像中,冠状动脉的线状特征信息得到增强,可用作构建冠状动脉知识通道图像;
步骤1.3、根据临床经验总结的冠脉知识特征构建多通道信息;
具体为,根据多尺度非局部积分滤波后图像,计算出造影区域的灰度和梯度特征,定义各个区域与理想造影区域的相似度,作为下一步对比剂候选区域确定算法的输入;
步骤1.4、利用变通道多尺度区域候选分类网络提取CAT图像中的三维冠脉区域;
起源于神经网络的深度学习网络,由于其在分类问题上优秀的学习和泛化能力,已成功应用于计算机视觉目标检测。在CTA三维图像中,冠状动脉也可以作为一种检测目标,通过训练不同的深度网络将冠状动脉从背景中分割出来。然而,仅根据灰度和纹理信息分割冠状动脉会包含很多假阳性区域。为了更精确地分割出冠状动脉图像,本实施例提出一种新颖深度学习网络——变通道多尺度区域候选分类网络,并结合根据临床经验总结的冠状动脉知识特征构建多通道信息,从而进行CTA冠状动脉区域的检测和分类,如图2所示。
在变通道多尺度区域候选分类网络,将原始CTA图像和知识特征通道送入深度学习网络中的区域候选模块;利用知识特征通道在表述冠状动脉区域的优势,确定候选冠状动脉区域;然后结合原始图像在候选冠状动脉区域的特征,利用网络中的特征提取和分类模块,进行分类和假阳性去除,分类网络输出最终的冠状动脉检测结果。
步骤2、如图3和图4所示,提取CTG DICOM图像并将CTA冠脉区域与CTG图像进行投射与融合
步骤2.1、提取冠状动脉CAG DICOM图像;
步骤2.2、将CTA中三维冠脉区域通过仿射变换与CAG图像中的对应冠脉进行配准,并投影到二维CAG上,确定冠脉对比剂导航区域;
具体的,本实施例中,步骤2.2包括以下步骤:
步骤2.2.1、冠脉造影对比剂导航区域确定
利用CAG影像采集时的角度和分辨率,确定仿射变换函数的参数,将CTA中冠脉分割结果进行仿射变换,并将变换后的三维冠脉图像投影到二维CAG图像上,作为冠脉的可能候选区域,也就是对比剂要经过的导航区域;
步骤2.2.2、冠脉造影对比剂区域自动检测和提取:
根据造影对比剂区域在CAG视频的灰度和纹理等信息,定义其中智不确定隶属度,从而作为粒子的属性,并定义当前区域及其邻域的中智相似度,来作为粒子滤波的代价函数,利用粒子滤波算法来追踪对比剂在导航区域内的流动情况,从而提取CAG图像中的冠脉区域。
步骤2.3、基于中智粒子滤波追踪算法追踪对比剂在导航区域内的流动情况;
具体的,本实施例中,步骤2.3包括以下步骤:
步骤2.3.1、提取冠脉目标特征:通过CTA图像中获得冠脉作为造影的导航区域,人机交互选定待跟踪的初始目标,计算跟踪目标的特征,作为描述目标的特征;
步骤2.3.2、特征搜索阶段:初始化提取了目标的特征后,在目标的附近以高斯分布来撒粒子,并根据导航区域来确定所撒的粒子的范围;每个粒子按照初始化时候得到的目标特征计算出它所处的位置的灰度特征,并计算该直方图向量与目标直方图向量的相似度;然后对每个粒子与目标计算出的相似度做归一化,使得所有粒子的相似度加起来和等于1;
步骤2.3.3、决策阶段:根据每个粒子与目标的相似程度做加权平均;根据最大的相似度候选区域的位置,则为目标的最可能位置;
步骤2.3.4、粒子重采样:根据粒子权重对粒子进行筛选,筛选过程中,既要大量保留权重大的粒子,又要有一小部分权重小的粒子,进而重复特征搜索阶段和决策阶段。
步骤2.4、自动检测并定位提取CAG冠脉对比剂区域;
步骤2.5、将CTA冠脉区域与CTG图像进行投射与融合;
步骤3、定量评估狭窄病变
步骤3.1、提取融合影像中冠脉的中心线,定位冠脉区域;
具体的,本实施例中,步骤3.1包括以下步骤:
步骤3.1.1、先在冠脉区域上做距离变换,计算各个点到边界的距离;
步骤3.1.2、根据各个点的距离,定义脊点作为冠脉中心线的候选点,并定义各个点之间的距离代价函数;
步骤3.13、利用最短路径Dijkstra算法来通过脊点寻找最短路径,作为冠脉中心线;
步骤3.1.4、在得到冠脉的中心线后,建立冠脉形状模型,根据先验知识构建正常冠脉的形状和直径;
步骤3.1.5、与实际图像计算出的冠脉进行比对,从而检测出冠脉异常区域,定义为冠脉病变区域。
步骤3.2、计算冠脉血流速度和载面积,获得局部血流量;
具体的,本实施例中,步骤3.2包括以下步骤:
步骤3.2.1、冠脉病变区域血流速度计算
冠脉病变区域的血流速度指血液在单位时间内通过的距离,根据对比剂的跟踪路线和CAG视频序列的时间,可以计算出血液在不同区域的速度(v),从速度上定量分析血液在流经病变区域的变化情况,其计算过程如下:
其中pp和pd是分别在狭窄区域近端和远端的取样点;Disk-1,Disk和Disk+1分别是取样点在第k-1,k和k+1帧上的距离;
步骤3.2.2、冠脉病变区域血流量计算
冠脉病变区域的血流量(Q)是指单位时间内流经病变区截面积的血液体积;利用极坐标系,考虑血管的某一截面积上任一微面积ds上血流速度为v(r,θ,t),则通过该微截面的血流量应为:
(6)式为血流量计算的基本公式,将(4)式变形后得出:
步骤3.2.3、根据冠脉对比剂在冠脉中的横截宽度和CTA的冠脉区域分割结果,利用样条函数来精确拟合出冠脉横截面;
由于获得的截面积个数有限(2-3个CAG角度),仅用CAG来重建冠脉的三维体积,精度不高,影响计算结果。而仅用CTA来计算冠脉体积,会受到冠脉钙化病变的影响,降低冠脉体积计算的准确性。如图3所示,本实施例结合CAG和CTA各自的优点,根据冠脉对比剂在冠脉中的横截宽度和CTA的冠脉区域分割结果,利用样条函数来精确拟合出冠脉横截面。
其中图5中的短虚线和长虚线分别为根据两次不同角度生成的对比剂的横断面和CTA冠脉分割结果横截面,图6为最后拟合结果;
步骤3.2.4、根据获得的最终横截面,可以获得冠脉病变狭窄远端的血流量Qd和病变近端无狭窄处血流量Qp,其计算过程如下:
其中p是取样点,S(p)为p点横截面,V(p)为取样点的血流速度,t为时间间隔;Qp和Qd分别为近端和远端取样点的血流量。
步骤3.3、计算病变远端及近端血流量比值,获得FFRCAD;
步骤3.4、根据FFRCAD值定量评估狭窄病变。
具体的,本实施例中,步骤3.3通过计算机辅助获得FFRCAD,理论上应该更能代表冠脉真实的血流储备分数,其公式为:
FFRCAD=Qd/Qp。 (11)
经过对本实施例提供轨道方法模拟机实验,该方法在实际中完全可行。与同类型的算法相比,本发明结果优异,结合CTA和CAG图像的各自特点,能实现对冠脉病变功能学的定量精确测量。
以上所述实施例仅为本发明较佳的具体实施方式,本发明的保护范围不限于此,任何熟悉本领域的技术人员在本发明披露的技术范围内,可显而易见地得到的技术方案的简单变化或等效替换,均属于本发明的保护范围。
Claims (7)
1.一种基于深度学习和中智理论的冠状动脉病变功能学定量方法,其特征在于,包括以下步骤:
步骤1、自动提取CTA中冠脉区域并进行区域处理
步骤1.1、提取冠状动脉CAT图像;
步骤1.2、利用多尺度非局部积分滤波算法对CAT图像进行预处理;
步骤1.3、根据临床经验总结的冠脉知识特征构建多通道信息;
步骤1.4、利用变通道多尺度区域候选分类网络提取CAT图像中的三维冠脉区域;
步骤2、提取CTG图像并将CTA冠脉区域与CTG图像进行投射与融合
步骤2.1、提取冠状动脉CAG图像;
步骤2.2、将CTA中三维冠脉区域通过仿射变换与CAG图像中的对应冠脉进行配准,并投影到二维CAG上,确定冠脉对比剂导航区域;
步骤2.3、基于中智粒子滤波追踪算法追踪对比剂在导航区域内的流动情况;
步骤2.4、自动检测并定位提取CAG冠脉对比剂区域;
步骤2.5、将CTA冠脉区域与CTG图像进行投射与融合;
步骤3、定量评估狭窄病变
步骤3.1、提取融合影像中冠脉的中心线,定位冠脉区域;
步骤3.2、计算冠脉血流速度和载面积,获得局部血流量;
步骤3.3、计算病变远端及近端血流量比值,获得FFRCAD;
步骤3.4、根据FFRCAD值定量评估狭窄病变。
3.根据权利要求1所述的基于深度学习和中智理论的冠状动脉病变功能学定量方法,其特征在于,所述步骤2.2包括以下步骤:
步骤2.2.1、冠脉造影对比剂导航区域确定
利用CAG影像采集时的角度和分辨率,确定仿射变换函数的参数,将CTA中冠脉分割结果进行仿射变换,并将变换后的三维冠脉图像投影到二维CAG图像上,作为冠脉的可能候选区域,也就是对比剂要经过的导航区域;
步骤2.2.2、冠脉造影对比剂区域自动检测和提取:
根据造影对比剂区域在CAG视频的灰度和纹理信息,定义其中智不确定隶属度,从而作为粒子的属性,并定义当前区域及其邻域的中智相似度,来作为粒子滤波的代价函数,利用粒子滤波算法来追踪对比剂在导航区域内的流动情况,从而提取CAG图像中的冠脉区域。
4.根据权利要求3所述的基于深度学习和中智理论的冠状动脉病变功能学定量方法,其特征在于,所述步骤2.3包括以下步骤:
步骤2.3.1、提取冠脉目标特征:通过CTA图像中获得冠脉作为造影的导航区域,人机交互选定待跟踪的初始目标,计算跟踪目标的特征,作为描述目标的特征;
步骤2.3.2、特征搜索阶段:初始化提取了目标的特征后,在目标的附近以高斯分布来撒粒子,并根据导航区域来确定所撒的粒子的范围;每个粒子按照初始化时候得到的目标特征计算出它所处的位置的灰度特征,并计算该直方图向量与目标直方图向量的相似度;然后对每个粒子与目标计算出的相似度做归一化,使得所有粒子的相似度加起来和等于1;
步骤2.3.3、决策阶段:根据每个粒子与目标的相似程度做加权平均;根据最大的相似度候选区域的位置,则为目标的最可能位置;
步骤2.3.4、粒子重采样:根据粒子权重对粒子进行筛选,筛选过程中,既要大量保留权重大的粒子,又要有一小部分权重小的粒子,进而重复特征搜索阶段和决策阶段。
5.根据权利要求1所述的基于深度学习和中智理论的冠状动脉病变功能学定量方法,其特征在于,所述步骤3.1包括以下步骤:
步骤3.1.1、先在冠脉区域上做距离变换,计算各个点到边界的距离;
步骤3.1.2、根据各个点的距离,定义脊点作为冠脉中心线的候选点,并定义各个点之间的距离代价函数;
步骤3.13、利用最短路径Dijkstra算法来通过脊点寻找最短路径,作为冠脉中心线;
步骤3.1.4、在得到冠脉的中心线后,建立冠脉形状模型,根据先验知识构建正常冠脉的形状和直径;
步骤3.1.5、与实际图像计算出的冠脉进行比对,从而检测出冠脉异常区域,定义为冠脉病变区域。
6.根据权利要求1所述的基于深度学习和中智理论的冠状动脉病变功能学定量方法,其特征在于,所述步骤3.2包括以下步骤:
步骤3.2.1、冠脉病变区域血流速度计算
冠脉病变区域的血流速度指血液在单位时间内通过的距离,根据对比剂的跟踪路线和CAG视频序列的时间,可以计算出血液在不同区域的速度v,从速度上定量分析血液在流经病变区域的变化情况,其计算过程如下:
其中pp和pd是分别在狭窄区域近端和远端的取样点;Disk-1,Disk和Disk+1分别是取样点在第k-1,k和k+1帧上的距离;
步骤3.2.2、冠脉病变区域血流量计算
冠脉病变区域的血流量Q是指单位时间内流经病变区截面积的血液体积;利用极坐标系,考虑血管的某一截面积上任一微面积ds上血流速度为v(r,θ,t),则通过该微截面的血流量应为:
(5)式为血流量计算的基本公式,将(4)式变形后得出:
步骤3.2.3、根据冠脉对比剂在冠脉中的横截宽度和CTA的冠脉区域分割结果,利用样条函数来精确拟合出冠脉横截面;
步骤3.2.4、根据获得的最终横截面,可以获得冠脉病变狭窄远端的血流量Qd和病变近端无狭窄处血流量Qp,其计算过程如下:
其中p是取样点,S(p)为p点横截面,V(p)为取样点的血流速度,t为时间间隔;Qp和Qd分别为近端和远端取样点的血流量。
7.根据权利要求1所述的基于深度学习和中智理论的冠状动脉病变功能学定量方法,其特征在于,所述步骤3.3通过计算机辅助获得FFRCAD,其公式为:
FFRCAD=Qd/Qp。 (11)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110192837.0A CN112837306B (zh) | 2021-02-20 | 2021-02-20 | 基于深度学习和中智理论的冠状动脉病变功能学定量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110192837.0A CN112837306B (zh) | 2021-02-20 | 2021-02-20 | 基于深度学习和中智理论的冠状动脉病变功能学定量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112837306A true CN112837306A (zh) | 2021-05-25 |
CN112837306B CN112837306B (zh) | 2022-11-22 |
Family
ID=75934058
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110192837.0A Active CN112837306B (zh) | 2021-02-20 | 2021-02-20 | 基于深度学习和中智理论的冠状动脉病变功能学定量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112837306B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113408152A (zh) * | 2021-07-23 | 2021-09-17 | 上海友脉科技有限责任公司 | 冠脉旁路移植仿真系统、方法、介质及电子设备 |
CN114862850A (zh) * | 2022-07-06 | 2022-08-05 | 深圳科亚医疗科技有限公司 | 一种血管医学图像的目标检测方法、装置及介质 |
TWI790179B (zh) * | 2022-07-27 | 2023-01-11 | 台灣整合心臟醫學協會 | 心導管影像辨識及評估方法 |
CN116649925A (zh) * | 2023-07-28 | 2023-08-29 | 杭州脉流科技有限公司 | 颅内动脉狭窄功能学评估的方法和装置 |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CA2361849A1 (en) * | 1999-03-08 | 2000-09-14 | Genentech, Inc. | Promotion or inhibition of angiogenesis and cardiovascularization |
CN202146329U (zh) * | 2011-07-19 | 2012-02-22 | 李悦 | 自动程控压力定量控制止血装置 |
CN102411795A (zh) * | 2011-08-16 | 2012-04-11 | 中国人民解放军第三军医大学第一附属医院 | 一种关于冠状动脉的多模图像的联合显示方法 |
AU2013202912A1 (en) * | 2007-06-08 | 2013-05-02 | Firmenich Incorporated | Modulation of chemosensory receptors and ligands associated therewith |
CN105228508A (zh) * | 2013-03-08 | 2016-01-06 | 新加坡健康服务有限公司 | 一种测定用于分类的危险评分的系统和方法 |
CN108010041A (zh) * | 2017-12-22 | 2018-05-08 | 数坤(北京)网络科技有限公司 | 基于深度学习神经网络级联模型的人体心脏冠脉提取方法 |
US20180146941A1 (en) * | 2016-11-28 | 2018-05-31 | Toshiba Medical Systems Corporation | Medical image-processing apparatus, x-ray ct apparatus, and medical image-processing method |
CN111227822A (zh) * | 2018-11-28 | 2020-06-05 | 苏州润心医疗器械有限公司 | 基于心肌血流量和ct图像的冠状动脉血流储备分数计算方法 |
CN111657883A (zh) * | 2020-06-03 | 2020-09-15 | 北京理工大学 | 基于序列造影的冠脉syntax分数自动计算方法与系统 |
US20210189342A1 (en) * | 2018-06-29 | 2021-06-24 | Verseau Therapeutics, Inc. | Compositions and methods for modulating monocyte and macrophage inflammatory phenotypes and immunotherapy uses thereof |
-
2021
- 2021-02-20 CN CN202110192837.0A patent/CN112837306B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CA2361849A1 (en) * | 1999-03-08 | 2000-09-14 | Genentech, Inc. | Promotion or inhibition of angiogenesis and cardiovascularization |
AU2013202912A1 (en) * | 2007-06-08 | 2013-05-02 | Firmenich Incorporated | Modulation of chemosensory receptors and ligands associated therewith |
CN202146329U (zh) * | 2011-07-19 | 2012-02-22 | 李悦 | 自动程控压力定量控制止血装置 |
CN102411795A (zh) * | 2011-08-16 | 2012-04-11 | 中国人民解放军第三军医大学第一附属医院 | 一种关于冠状动脉的多模图像的联合显示方法 |
CN105228508A (zh) * | 2013-03-08 | 2016-01-06 | 新加坡健康服务有限公司 | 一种测定用于分类的危险评分的系统和方法 |
US20180146941A1 (en) * | 2016-11-28 | 2018-05-31 | Toshiba Medical Systems Corporation | Medical image-processing apparatus, x-ray ct apparatus, and medical image-processing method |
CN108010041A (zh) * | 2017-12-22 | 2018-05-08 | 数坤(北京)网络科技有限公司 | 基于深度学习神经网络级联模型的人体心脏冠脉提取方法 |
US20210189342A1 (en) * | 2018-06-29 | 2021-06-24 | Verseau Therapeutics, Inc. | Compositions and methods for modulating monocyte and macrophage inflammatory phenotypes and immunotherapy uses thereof |
CN111227822A (zh) * | 2018-11-28 | 2020-06-05 | 苏州润心医疗器械有限公司 | 基于心肌血流量和ct图像的冠状动脉血流储备分数计算方法 |
CN111657883A (zh) * | 2020-06-03 | 2020-09-15 | 北京理工大学 | 基于序列造影的冠脉syntax分数自动计算方法与系统 |
Non-Patent Citations (8)
Title |
---|
CHRISTIANTESCHEMD 等: "Influence of Coronary Calcium on Diagnostic Performance of Machine Learning CT-FFR: Results From MACHINE Registry", 《JACC: CARDIOVASCULAR IMAGING》, vol. 13, no. 3, 31 March 2020 (2020-03-31), pages 760 - 770, XP086073618, DOI: 10.1016/j.jcmg.2019.06.027 * |
DAMINI DEY: "Machine-Learning CT-FFR and Extensive Coronary Calcium: Overcoming the Achilles Heel of Coronary Computed Tomography Angiography", 《JACC: CARDIOVASCULAR IMAGING》, vol. 13, no. 3, 31 March 2020 (2020-03-31), pages 771 - 773, XP086073620, DOI: 10.1016/j.jcmg.2019.08.011 * |
GUOHONG 等: "e0154 Expression of tumour necrosis factor-α converting enzyme and tumour necrosis factor-α in rats with alcoholic cardiomyopathy", 《HEART》, vol. 96, no. 3, 31 December 2010 (2010-12-31), pages 49 - 50 * |
NAHOKO KATO 等: "Relative atherosclerotic plaque volume by CT coronary angiography trumps conventional stenosis assessment for identifying flow-limiting lesions", 《THE INTERNATIONAL JOURNAL OF CARDIOVASCULAR IMAGING》, vol. 33, no. 11, 7 June 2017 (2017-06-07), pages 1847 - 1855, XP036356806, DOI: 10.1007/s10554-017-1186-z * |
安梦楠 等: "冠状动脉狭窄病变功能学评估新进展", 《中国介入心脏病学杂志》, vol. 26, no. 11, 30 November 2018 (2018-11-30), pages 638 - 640 * |
李浚利: "人工智能评估CCTA冠状动脉狭窄准确性的临床研究", 《中国优秀博硕士学位论文全文数据库(硕士)医药卫生科技辑》, no. 02, 15 February 2021 (2021-02-15), pages 062 - 583 * |
翁婷雯 等: "计算机断层成像血流储备分数对冠状动脉中度病变的稳定型冠状动脉性心脏病患者预后的预测作用", 《上海医学》, vol. 43, no. 1, 15 January 2020 (2020-01-15), pages 14 - 17 * |
薛竟宜: "联合应用EPO和G-CSF促进梗死心肌血管再生及其机制的实验研究", 《中国优秀博硕士学位论文全文数据库(博士)医药卫生科技辑》, no. 01, 15 January 2018 (2018-01-15), pages 062 - 17 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113408152A (zh) * | 2021-07-23 | 2021-09-17 | 上海友脉科技有限责任公司 | 冠脉旁路移植仿真系统、方法、介质及电子设备 |
CN113408152B (zh) * | 2021-07-23 | 2023-07-25 | 上海友脉科技有限责任公司 | 冠脉旁路移植仿真系统、方法、介质及电子设备 |
CN114862850A (zh) * | 2022-07-06 | 2022-08-05 | 深圳科亚医疗科技有限公司 | 一种血管医学图像的目标检测方法、装置及介质 |
CN114862850B (zh) * | 2022-07-06 | 2022-09-20 | 深圳科亚医疗科技有限公司 | 一种血管医学图像的目标检测方法、装置及介质 |
TWI790179B (zh) * | 2022-07-27 | 2023-01-11 | 台灣整合心臟醫學協會 | 心導管影像辨識及評估方法 |
CN116649925A (zh) * | 2023-07-28 | 2023-08-29 | 杭州脉流科技有限公司 | 颅内动脉狭窄功能学评估的方法和装置 |
CN116649925B (zh) * | 2023-07-28 | 2023-10-31 | 杭州脉流科技有限公司 | 颅内动脉狭窄功能学评估的方法和装置 |
Also Published As
Publication number | Publication date |
---|---|
CN112837306B (zh) | 2022-11-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112837306B (zh) | 基于深度学习和中智理论的冠状动脉病变功能学定量方法 | |
Balocco et al. | Standardized evaluation methodology and reference database for evaluating IVUS image segmentation | |
Haris et al. | Model-based morphological segmentation and labeling of coronary angiograms | |
Friman et al. | Multiple hypothesis template tracking of small 3D vessel structures | |
Li et al. | Vessels as 4-D curves: Global minimal 4-D paths to extract 3-D tubular surfaces and centerlines | |
Cheng et al. | Accurate vessel segmentation with constrained B-snake | |
CN110448319B (zh) | 基于造影影像及冠状动脉的血流速度计算方法 | |
EP3723038A1 (en) | Fast calculation method and system employing plaque stability index of medical image sequence | |
CN105741299B (zh) | 一种冠状动脉ct血管造影图像分割方法 | |
Cheng et al. | Automated delineation of calcified vessels in mammography by tracking with uncertainty and graphical linking techniques | |
CN112597982B (zh) | 基于人工智能的图像分类方法、装置、设备和介质 | |
Meiburger et al. | Carotid Ultrasound Boundary Study (CUBS): Technical considerations on an open multi-center analysis of computerized measurement systems for intima-media thickness measurement on common carotid artery longitudinal B-mode ultrasound scans | |
Krishnan et al. | Analysis of time-dependent flow-sensitive PC-MRI data | |
Blondel et al. | Automatic trinocular 3D reconstruction of coronary artery centerlines from rotational X-ray angiography | |
Xiao et al. | Automatic vasculature identification in coronary angiograms by adaptive geometrical tracking | |
CN111445473A (zh) | 基于血管内超声图像序列多角度重建的血管膜精确分割方法及系统 | |
Samet et al. | A new methodology for blood vessel segmentation on lung CT images | |
Chakshu et al. | Automating fractional flow reserve (FFR) calculation from CT scans: A rapid workflow using unsupervised learning and computational fluid dynamics | |
Wang et al. | Segmentation of intravascular ultrasound images based on convex–concave adjustment in extreme regions | |
Sebbe | Computer-aided Diagnosis of Pulmonary Embolism in Opacified CT Images | |
Saif et al. | Computer Vision-based Efficient Segmentation Method for Left Ventricular Epicardium and Endocardium using Deep Learning | |
Sakellarios et al. | IVUS image processing methodologies | |
Baral et al. | A review on machine learning for arterial extraction and quantitative assessment on invasive coronary angiograms | |
Popov et al. | A Review of Modern Approaches for Coronary Angiography Imaging Analysis | |
Bas et al. | Automated extraction of blood vessel networks from 3d microscopy image stacks via multi-scale principal curve tracing |
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 |