CN113409447B - 一种基于多切片组合的冠状动脉分割方法和装置 - Google Patents
一种基于多切片组合的冠状动脉分割方法和装置 Download PDFInfo
- Publication number
- CN113409447B CN113409447B CN202110644581.2A CN202110644581A CN113409447B CN 113409447 B CN113409447 B CN 113409447B CN 202110644581 A CN202110644581 A CN 202110644581A CN 113409447 B CN113409447 B CN 113409447B
- Authority
- CN
- China
- Prior art keywords
- dimensional
- slice
- segmentation
- coronary artery
- slices
- 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
- 230000011218 segmentation Effects 0.000 title claims abstract description 141
- 210000004351 coronary vessel Anatomy 0.000 title claims abstract description 69
- 238000000034 method Methods 0.000 title claims abstract description 43
- 230000002526 effect on cardiovascular system Effects 0.000 claims abstract description 59
- 238000012549 training Methods 0.000 claims abstract description 16
- 238000012545 processing Methods 0.000 claims description 13
- 230000006870 function Effects 0.000 claims description 10
- 230000010354 integration Effects 0.000 claims description 9
- 238000012307 MRI technique Methods 0.000 claims description 5
- 239000003086 colorant Substances 0.000 claims description 5
- 150000001875 compounds Chemical class 0.000 claims description 3
- 230000007547 defect Effects 0.000 abstract description 3
- 210000004204 blood vessel Anatomy 0.000 description 5
- 238000003384 imaging method Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 208000024172 Cardiovascular disease Diseases 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 230000006835 compression Effects 0.000 description 2
- 238000007906 compression Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000018109 developmental process Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 210000001835 viscera Anatomy 0.000 description 2
- 238000013528 artificial neural network Methods 0.000 description 1
- 239000008280 blood Substances 0.000 description 1
- 210000004369 blood Anatomy 0.000 description 1
- 210000005242 cardiac chamber Anatomy 0.000 description 1
- 238000002591 computed tomography Methods 0.000 description 1
- 238000013170 computed tomography imaging Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 208000029078 coronary artery disease Diseases 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- PCHJSUWPFVWCPO-UHFFFAOYSA-N gold Chemical compound [Au] PCHJSUWPFVWCPO-UHFFFAOYSA-N 0.000 description 1
- 238000003709 image segmentation Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000004445 quantitative analysis Methods 0.000 description 1
- 238000002601 radiography Methods 0.000 description 1
- 238000011524 similarity measure Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
-
- 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
- G06T5/00—Image enhancement or restoration
- G06T5/50—Image enhancement or restoration using two or more images, e.g. averaging or subtraction
-
- 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/90—Determination of colour characteristics
-
- 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/10—Image acquisition modality
- G06T2207/10072—Tomographic images
- G06T2207/10088—Magnetic resonance imaging [MRI]
-
- 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/20021—Dividing image into blocks, subimages or windows
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Software Systems (AREA)
- Computing Systems (AREA)
- Data Mining & Analysis (AREA)
- Mathematical Physics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Artificial Intelligence (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- Computational Linguistics (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- General Health & Medical Sciences (AREA)
- Molecular Biology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Geometry (AREA)
- Computer Graphics (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Image Analysis (AREA)
Abstract
本发明公开了一种基于多切片组合的冠状动脉分割方法和装置,所述方法包括:获取三维心血管图像;对三维心血管图像进行切片处理,获得二维切片;将二维切片以m种组合方式组合后,输入分割模型中进行训练,获得m个训练好的分割模型;将m个训练好的分割模型进行集成,获得集成分割模型;以相同的m种组合方式,将待分割的三维心血管图像的二维切片组合后输入集成分割模型中进行分割,获得二维冠状动脉预测图片;对二维冠状动脉预测图片进行三维重构,获得冠状动脉三维图像。本发明充分考虑了二维切片间的邻域信息,弥补了单一切片方式输入时邻域信息的缺失,跨切片进行一致性预测,提升了冠状动脉分割结果的准确性。
Description
技术领域
本发明涉及冠状动脉分割的技术领域,更具体地,涉及一种基于多切片组合 的冠状动脉分割方法和装置。
背景技术
冠状动脉疾病是最常见的心血管疾病之一,冠状动脉的分割对心血管疾病的 准确定位和定量分析具有重要意义。近年来,非侵入性成像技术的发展使诊断的 准确性发生了革命性的变化,现代设备可以对内部器官的亚毫米细节进行成像。 随着计算机技术的发展,三维成像得以实现,较多用于心脏部位的计算机断层扫 描造影可以根据强度很好的签别内脏器官。在心脏部位的计算机断层扫描造影中, 充血的血管相比于周围的组织更加明亮一些,这使得手动跟踪冠脉结构具有合理 的准确性。目前对于心血管三维图像的传统分割方法的黄金标准是由临床专家手 工完成的,但是手工分割不仅耗时,也存在一定程度主观上的偏差和错误,造成 分割结果不准确;利用计算机技术对三维心血管图像分割通常是在2D切片上逐 个完成,再将这些2D切片在多个截面上分组成几何上一致的3D心血管图像,构建分层的区域图是计算密集型的,不仅需要精确分割的2D切片,还需要跨切 片进行一致的预测,但该方法提取的区域特征没有包含足够的邻域信息,分割准 确性不高。
2017年12月29日公开的中国专利CN107527341A提供了一种血管造影图 像的处理方法,包括:在血管造影图像中,确定包含冠状动脉的目标区域;对所 述包含冠状动脉的目标区域进行线性增强;确定冠状动脉的预设灰度阈值范围; 保存灰度值在预设灰度阈值范围的体素点;处理所述体素点形成冠状动脉图像, 从而将冠状动脉图像与心脏腔室图像剥离。该方法基于血管分割算法,将血管分 割问题看作像素分类问题,将每个像素判断为血管或非血管,该算法的设计复杂 性与应用局限性较高,并且特征提取算法与分类器结合的多样性使分割准确性不 高。
发明内容
本发明为克服上述现有技术对三维心血管图像分割结果不准确的缺陷,提供 一种基于多切片组合的冠状动脉分割方法和装置,充分考虑了邻域信息,提高了 分割结果的准确性。
为解决上述技术问题,本发明的技术方案如下:
本发明提供一种基于多切片组合的冠状动脉分割方法,所述方法包括以下步 骤:
S1:获取三维心血管图像;
S2:对三维心血管图像进行切片处理,获得二维切片;
S3:将二维切片以m种组合方式组合后,输入分割模型中进行训练,获得 m个训练好的分割模型;
S4:将m个训练好的分割模型进行集成,获得集成分割模型;
S5:将待分割的三维心血管图像进行切片处理,得到待分割的三维心血管图 像的二维切片,将其以与S3中相同的m种组合方式进行组合后,输入集成分割 模型中进行分割,获得二维冠状动脉预测图片;
S6:对二维冠状动脉预测图片进行三维重构,获得冠状动脉三维图像。
优选地,所述三维心血管图像由CT或MRI技术获得。CT或MRI技术均属 于非侵入性成像技术,安全性高。
优选地,所述对三维心血管图像进行切片处理的具体方法为:
对三维心血管图像构建最大外切长方体,大小为a*b*c,a、b、c分别为三 维心血管图像最大外切长方体长、宽、高数值;将三维心血管图像的最大外切长 方体置于三维坐标系中,在Z轴方向进行切片处理,获得平行于XoY平面的二 维切片,共c张,大小为a*b,则Z轴方向上的第i张二维切片记为Ii,i= 1,…,p,p+1,p+2,…,c-q-1,c-q,…,c。
优选地,对二维切片组合前还需舍去Z轴方向上的正数前p张和倒数后q 张二维切片。心脏周围的三维心血管图像的二维切片,正数前p张和倒数后q张 通常不包含冠状动脉,舍去p+q张图片减少了训练时的数据量,且不会影响分割 结果的准确性。
优选地,所述S3中,所述S3中,二维切片以m种组合方式中的m=3,具 体组合方式为:
将第i张二维切片记为Ii,相应的切片真实标签记为Li,舍去Z轴方向上的 正数前p张和倒数后q张二维切片后,i=p+1,p+2,…,c-q-1,c-q;以单 切片的方式输入分割模型时,输入Ii;以三切片组合的方式输入分割模型时,输 入[Ii-1,Ii,Ii+1];以五切片组合的方式输入分割模型时,输入[Ii-2,Ii-1,Ii,Ii+1,Ii+2]。
优选地,所述分割模型为U-Net模型,获得m个训练好的分割模型的具体 方法为:
将二维切片以m种组合方式组合后分别输入m个U-Net模型中,每个U-Net 模型均设置二分类交叉熵损失函数,损失函数收敛时,获得m个训练好的分割 模型;二分类交叉熵损失函数为:
优选地,所述步骤S4中,利用投票集成法将m个训练好的分割模型进行集 成。投票集成法属于集成学习,可以有效提高模型的泛化能力,减少模型的错误 率。
优选地,所述S5中,获得二维冠状动脉预测图片的具体方法为:
S5.1:将待分割的三维心血管图像的二维切片以单切片、三切片组合、五切片 组合的方式输入集成分割模型中,获得二维切片上像素点的第一预测概率、第二 预测概率和第三预测概率;
S5.2:比较像素点的第一预测概率、第二预测概率和第三预测概率,选择数值 最大的预测概率作为该像素点的最终预测概率;
S5.3:若像素点的最终预测概率大于0.5,则记为1,该像素点在二维冠状动 脉预测图片显示为白色;反之,则记为0,该像素点在二维冠状动脉预测图片显 示为黑色;
S5.4:对二维切片上所有像素点执行S5.1-S5.3,所有确定颜色的像素点构成 二维冠状动脉预测图片。
优选地,所述获得二维切片上像素点的第一预测概率、第二预测概率和第三 预测概率的公式为:
本发明还提供了一种基于多切片组合的冠状动脉分割装置,所述装置包括:
图像获取模块,用于获取三维心血管图像;
切片模块,用于对三维心血管图像进行切片处理,获得二维切片;
组合训练模块,用于将二维切片以m种组合方式组合后,输入分割模型中 进行训练,获得m个训练好的分割模型;
模型集成模块,用于将m个训练好的分割模型进行集成,获得集成分割模 型;
待分割图像输入模块,用于将待分割的三维心血管图像进行切片处理,得到 待分割的三维心血管图像的二维切片,将其以与S3中相同的m种组合方式进行 组合后,输入集成分割模型中进行分割,获得二维冠状动脉预测图片;
三维重构模块,用于对二维冠状动脉预测图片进行三维重构,获得冠状动脉 三维图像。
与现有技术相比,本发明技术方案的有益效果是:
本发明采用了不同组合方式的切片对现有分割模型进行训练,训练好的分割 模型具有不同的参数,集成后获得集成分割模型,提高了模型的分割准确性;将 待分割的三维心血管图像的二维切片以不同组合方式组合后输入集成分割模型 中进行分割,组合后二维切片包括了丰富的邻域信息,弥补了单一切片方式输入 时切片间邻域信息的缺失的缺陷,跨切片进行一致性预测,提升了冠状动脉分割 结果的准确性。
附图说明
图1为实施例1所述的一种基于多切片组合的冠状动脉分割方法的流程图;
图2为实施例1所述的二维切片不同组合方式示意图;
图3为实施例1所述的不同组合方式二维切片获得的预测图片示意图;
图4为实施例2所述的一种基于多切片组合的冠状动脉分割转置的示意图。
具体实施方式
附图仅用于示例性说明,不能理解为对本专利的限制;
为了更好说明本实施例,附图某些部件会有省略、放大或缩小,并不代表实 际产品的尺寸;
对于本领域技术人员来说,附图中某些公知结构及其说明可能省略是可以理 解的。
下面结合附图和实施例对本发明的技术方案做进一步的说明。
实施例1
本实施例提供一种基于多切片组合的冠状动脉分割方法,如图1所示,所述 方法包括以下步骤:
S1:获取三维心血管图像;
所述三维心血管图像由CT或MRI技术获得。CT或MRI技术均属于非侵入 性成像技术,安全性高。
S2:对三维心血管图像进行切片处理,获得二维切片;
对三维心血管图像构建最大外切长方体,大小为a*b*c,a、b、c分别为三 维心血管图像最大外切长方体长、宽、高数值;将三维心血管图像的最大外切长 方体置于三维坐标系中,在Z轴方向进行切片处理,获得平行于XoY平面的二 维切片,共c张,大小为a*b,则Z轴方向上的第i张二维切片记为Ii,i= 1,…,p,p+1,p+2,…,c-q-1,c-q,…,c。
舍去Z轴方向上的正数前p张和倒数后q张二维切片;心脏周围的三维心血 管图像的二维切片,正数前p张以及倒数后q张通常不包含冠状动脉,舍去这 p+q张图片减少了训练时的数据量,且不会影响分割结果的准确性。
S3:将二维切片以m种组合方式组合后,输入分割模型中进行训练,获得 m个训练好的分割模型;
本实施例中,m=3,将二维切片以3种组合方式组合;将第i张二维切片记 为Ii,相应的切片真实标签记为Li,舍去Z轴方向上的正数前p张和倒数后q张 二维切片后,i=p+1,p+2,…,c-q-1,c-q;以单切片的方式输入分割模型 时,输入Ii;以三切片组合的方式输入分割模型时,输入[Ii-1,Ii,Ii+1];以五切片 组合的方式输入分割模型时,输入[Ii-2,Ii-1,Ii,Ii+1,Ii+2]。
本实施例中,分割模型为U-Net模型。
获得m个训练好的分割模型的具体方法为:
将二维切片以m种组合方式组合后分别输入m个U-Net模型中,每个U-Net 模型均设置二分类交叉熵损失函数,损失函数收敛时,获得m个训练好的分割 模型;二分类交叉熵损失函数为:
S4:采用投票集成法将m个训练好的分割模型进行集成,获得集成分割模 型;
S5:将待分割的三维心血管图像进行切片处理,得到待分割的三维心血管图 像的二维切片,将其以与S3中相同的m种组合方式进行组合后,输入集成分割 模型中进行分割,获得二维冠状动脉预测图片;具体方法为:
S5.1:将待分割的三维心血管图像的二维切片以单切片、三切片组合、五切片 组合的方式输入集成分割模型中,获得二维切片上像素点的第一预测概率、第二 预测概率和第三预测概率;公式为:
S5.2:比较像素点的第一预测概率、第二预测概率和第三预测概率,选择数值 最大的预测概率作为该像素点的最终预测概率;
S5.3:若像素点的最终预测概率大于0.5,则记为1,该像素点在二维冠状动 脉预测图片显示为白色;反之,则记为0,该像素点在二维冠状动脉预测图片显 示为黑色;
S5.4:对二维切片上所有像素点执行S5.1-S5.3,确定颜色的所有像素点构成 二维冠状动脉预测图片;
S6:对二维冠状动脉预测图片进行三维重构,获得冠状动脉三维图像。
在具体实施过程中,
利用CT技术获得三维心血管图像,对三维心血管图像构建最大外切长方体, 大小为512×512×275,在Z轴方向进行切片处理,获得平行于XoY平面的二 维切片,获得275张大小为512×512的二维切片,则Z轴方向上的第i张二维 切片记为Ii,i=1,2,…,275,相应的切片真实标签记为Li;本实施例证,因为心 脏周围的三维心血管图像的二维切片,正数前两张以及倒数两张通常不包含冠状 动脉,p=q=2,舍去I1,I2,I274,I275这4张二维切片,此时i=3,4,…,273;
如图2所示,将二维切片以三种组合方式进行组合后,输入U-Net模型中进 行训练,获得训练好的第一分割模型、第二分割模型和第三分割模型;具体为: 以单切片Ii的方式输入U-Net模型时,大小为512×512×1;以三切片组合 [Ii-1,Ii,Ii+1]的方式输入U-Net模型时,大小为512×512×3,;以五切片组合 [Ii-2,Ii-1,Ii,Ii+1,Ii+2]的方式输入U-Net模型时,大小为512×512×5,;采用切 片多种方式组合输入到分割模型中进行训练,有利于模型获得邻域信息,提高分 割结果的准确性;
选取分割模型为U-Net模型,该模型是基于全卷积神经网络(FCN)进行改 进得到的,主要由捕获上下文信息的压缩路径和实现精确定位的与压缩路径对称 的扩展路径组成;这种网络充分利用少量的图像进行端到端的训练,可以对每个 像素点进行分割,获得更高的分割准确率,同时这种网络结构能够更好的学习图 像深层特征。采用投票集成法将训练好的第一分割模型、第二分割模型和第三分 割模型进行集成,获得集成分割模型;
将待分割的三维心血管图像的二维切片以单切片方式输入集成分割模型中, 对单切片输入的二维切片上的每个像素点都计算获得相应的第一预测概率;将待 分割的三维心血管图像的二维切片以三切片组合的方式输入集成分割模型中,对 三切片组合输入的二维切片上的每个像素点都计算获得相应的获得第二预测概 率;将待分割的三维心血管图像的二维切片以五切片组合的方式输入集成分割模 型中,对五切片组合输入二维切片上的每个像素点都计算获得相应的获得第三预 测概率;比较像素点的第一预测概率、第二预测概率和第三预测概率,选择数值 最大的预测概率作为该像素点的最终预测概率;若像素点的最终预测概率大于 0.5,则记为1,该像素点在二维冠状动脉预测图片显示为白色;反之,则记为0, 该像素点在二维冠状动脉预测图片显示为黑色,所有确定颜色的像素点构成二维 冠状动脉预测图片。如图3所示,单切片预测图片由单切片输入后获得,三切片预测图片由三切片组合输入后获得,五切片预测图片由五切片组合输入后获得; 将所有二维切片的二维冠状动脉预测图片进行三维重构,获得冠状动脉三维图像。
本实施例中所述的单切片、三切片组合、五切片组合的三种组合方式仅仅是 为清楚地说明本发明所作的举例,也可以是七切片组合、九切片组合,也可以是 更多种组合方式组合后对分割模型进行训练。
为了定量分析切片不同组合方式以及集成分割模型的分割效果,本实施例使 用了Dice系数作为指标对分割效果进行评估。Dice系数是一种集合相似度度量 函数,通常用于计算两个样本的相似度(值范围为[0,1]),其计算公式为:
式中,s表示Dice系数,|X|表示预测图片的面积,|Y|表示二维切片真实图 片的面积,|X∩Y|表示预测图片和真实图片的交集,Dice系数数值越大,则预测 图与真实图片越相似;使用Dice系数对单切片、三切片组合、五切片组合和集 成分割模型的分割结果与二维切片对应的切片真实标签Li之间的相似度进行评 估,评估结果如下表所示:
分割方法 | Dice系数 |
单切片 | 0.5612 |
三切片组合 | 0.5942 |
五切片组合 | 0.5973 |
集成分割模型 | 0.6118 |
本实施例训练出的集成分割模型的Dice系数高于单一切片方式输入的Dice 系数,说明本实施例提供的方法具有以下优点:
(1)相较于传统的手工分割方法,本实施例不需要临床专家介入分割过程, 便可获得较好的分割效果,减少了经历和时间的投入;
(2)弥补了单一切片方式输入对邻域信息的缺失,跨切片进行一致性预测, 提升了冠状动脉分割结果的准确性;
(3)采用了不同组合的切片对现有分割模型U-Net进行训练,并对训练得 到的、具有不同的参数的U-Net模型进行了集成,提高了模型的分割准确性。
实施例2
本实施例提供一种基于多切片组合的冠状动脉分割装置,如图4所示,所述 装置包括:
图像获取模块,用于获取三维心血管图像;
切片模块,用于对三维心血管图像进行切片处理,获得二维切片;
组合训练模块,用于将二维切片以m种组合方式组合后,输入分割模型中 进行训练,获得m个训练好的分割模型;
模型集成模块,用于将m个训练好的分割模型进行集成,获得集成分割模 型;
待分割图像输入模块,用于将待分割的三维心血管图像进行切片处理,得到 待分割的三维心血管图像的二维切片,将其以与S3中相同的m种组合方式进行 组合后,输入集成分割模型中进行分割,获得二维冠状动脉预测图片;;
三维重构模块,用于对二维冠状动脉预测图片进行三维重构,获得冠状动脉 三维图像。
显然,本发明的上述实施例仅仅是为清楚地说明本发明所作的举例,而并非 是对本发明的实施方式的限定。对于所属领域的普通技术人员来说,在上述说明 的基础上还可以做出其它不同形式的变化或变动。这里无需也无法对所有的实施 方式予以穷举。凡在本发明的精神和原则之内所作的任何修改、等同替换和改进 等,均应包含在本发明权利要求的保护范围之内。
Claims (8)
1.一种基于多切片组合的冠状动脉分割方法,其特征在于,所述方法包括以下步骤:
S1:获取三维心血管图像;
S2:对三维心血管图像进行切片处理,获得二维切片;
S3:将二维切片以m种组合方式组合后,输入分割模型中进行训练,获得m个训练好的分割模型;
二维切片以m种组合方式中的m=3,具体组合方式为:
将第i张二维切片记为Ii,相应的切片真实标签记为Li,舍去Z轴方向上的正数前p张和倒数后q张二维切片后,i=p+1,p+2,…,c-q-1,c-q;以单切片的方式输入分割模型时,输入Ii;以三切片组合的方式输入分割模型时,输入[Ii-1,Ii,Ii+1];以五切片组合的方式输入分割模型时,输入[Ii-2,Ii-1,Ii,Ii+1,Ii+2];
S4:将m个训练好的分割模型进行集成,获得集成分割模型;
S5:将待分割的三维心血管图像进行切片处理,得到待分割的三维心血管图像的二维切片,将其以与S3中相同的m种组合方式进行组合后,输入集成分割模型中进行分割,获得二维冠状动脉预测图片;具体方法为:
S5.1:将待分割的三维心血管图像的二维切片以单切片、三切片组合、五切片组合的方式输入集成分割模型中,获得二维切片上像素点的第一预测概率、第二预测概率和第三预测概率;
S5.2:比较像素点的第一预测概率、第二预测概率和第三预测概率,选择数值最大的预测概率作为该像素点的最终预测概率;
S5.3:若像素点的最终预测概率大于0.5,则记为1,该像素点在二维冠状动脉预测图片显示为白色;反之,则记为0,该像素点在二维冠状动脉预测图片显示为黑色;
S5.4:对二维切片上所有像素点执行S5.1-S5.3,确定颜色的所有像素点构成二维冠状动脉预测图片;
S6:对二维冠状动脉预测图片进行三维重构,获得冠状动脉三维图像。
2.根据权利要求1所述的基于多切片组合的冠状动脉分割方法,其特征在于,所述三维心血管图像由CT或MRI技术获得。
3.根据权利要求2所述的基于多切片组合的冠状动脉分割方法,其特征在于,所述对三维心血管图像进行切片处理的具体方法为:
对三维心血管图像构建最大外切长方体,大小为a*b*c,a、b、c分别为三维心血管图像最大外切长方体长、宽、高数值;将三维心血管图像的最大外切长方体置于三维坐标系中,在Z轴方向进行切片处理,获得平行于XoY平面的二维切片,共c张,大小为a*b,则Z轴方向上的第i张二维切片记为Ii,i=1,…,p,p+1,p+2,…,c-q-1,c-q,…,c。
4.根据权利要求3所述的基于多切片组合的冠状动脉分割方法,其特征在于,对二维切片组合前还需舍去Z轴方向上的正数前p张和倒数后q张二维切片。
6.根据权利要求5所述的基于多切片组合的冠状动脉分割方法,其特征在于,所述步骤S4中,利用投票集成法将m个训练好的分割模型进行集成。
8.一种基于多切片组合的冠状动脉分割装置,其特征在于,所述装置包括:
图像获取模块,用于获取三维心血管图像;
切片模块,用于对三维心血管图像进行切片处理,获得二维切片;
组合训练模块,用于将二维切片以m种组合方式组合后,输入分割模型中进行训练,获得m个训练好的分割模型;
二维切片以m种组合方式中的m=3,具体组合方式为:
将第i张二维切片记为Ii,相应的切片真实标签记为Li,舍去Z轴方向上的正数前p张和倒数后q张二维切片后,i=p+1,p+2,…,c-q-1,c-q;以单切片的方式输入分割模型时,输入Ii;以三切片组合的方式输入分割模型时,输入[Ii-1,Ii,Ii+1];以五切片组合的方式输入分割模型时,输入[Ii-2,Ii-1,Ii,Ii+1,Ii+2];
模型集成模块,用于将m个训练好的分割模型进行集成,获得集成分割模型;
待分割图像输入模块,用于将待分割的三维心血管图像进行切片处理,得到待分割的三维心血管图像的二维切片,将其以与组合训练模块中相同的m种组合方式进行组合后,输入集成分割模型中进行分割,获得二维冠状动脉预测图片;具体方法为:
将待分割的三维心血管图像的二维切片以单切片、三切片组合、五切片组合的方式输入集成分割模型中,获得二维切片上像素点的第一预测概率、第二预测概率和第三预测概率;
比较像素点的第一预测概率、第二预测概率和第三预测概率,选择数值最大的预测概率作为该像素点的最终预测概率;
若像素点的最终预测概率大于0.5,则记为1,该像素点在二维冠状动脉预测图片显示为白色;反之,则记为0,该像素点在二维冠状动脉预测图片显示为黑色;
对二维切片上所有像素点执行上述方法,确定颜色的所有像素点构成二维冠状动脉预测图片
三维重构模块,用于对二维冠状动脉预测图片进行三维重构,获得冠状动脉三维图像。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110644581.2A CN113409447B (zh) | 2021-06-09 | 2021-06-09 | 一种基于多切片组合的冠状动脉分割方法和装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110644581.2A CN113409447B (zh) | 2021-06-09 | 2021-06-09 | 一种基于多切片组合的冠状动脉分割方法和装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113409447A CN113409447A (zh) | 2021-09-17 |
CN113409447B true CN113409447B (zh) | 2023-01-24 |
Family
ID=77683287
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110644581.2A Active CN113409447B (zh) | 2021-06-09 | 2021-06-09 | 一种基于多切片组合的冠状动脉分割方法和装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113409447B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115984251B (zh) * | 2023-02-14 | 2023-05-09 | 成都泰莱生物科技有限公司 | 基于肺部ct与多基因甲基化的肺结节分类方法及产品 |
Family Cites Families (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9514530B2 (en) * | 2014-04-16 | 2016-12-06 | Heartflow, Inc. | Systems and methods for image-based object modeling using multiple image acquisitions or reconstructions |
CN108447062B (zh) * | 2018-02-01 | 2021-04-20 | 浙江大学 | 一种基于多尺度混合分割模型的病理切片非常规细胞的分割方法 |
CN109671076A (zh) * | 2018-12-20 | 2019-04-23 | 上海联影智能医疗科技有限公司 | 血管分割方法、装置、电子设备及存储介质 |
CN111476796B (zh) * | 2020-03-10 | 2023-04-18 | 西北大学 | 一种结合多种网络的半监督冠状动脉分割系统及分割方法 |
CN111784653B (zh) * | 2020-06-28 | 2023-08-01 | 西安电子科技大学 | 基于形状约束的多尺度网络mri胰腺轮廓定位方法 |
CN112508827B (zh) * | 2020-11-06 | 2022-04-22 | 中南大学湘雅医院 | 一种基于深度学习的融合多场景的危及器官的分割方法 |
CN112330682B (zh) * | 2020-11-09 | 2022-08-26 | 重庆邮电大学 | 一种基于深度卷积神经网络的工业ct图像分割方法 |
CN112561860B (zh) * | 2020-11-23 | 2022-05-03 | 重庆邮电大学 | 一种基于先验形状约束的BCA-UNet肝脏分割方法 |
-
2021
- 2021-06-09 CN CN202110644581.2A patent/CN113409447B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN113409447A (zh) | 2021-09-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US11562491B2 (en) | Automatic pancreas CT segmentation method based on a saliency-aware densely connected dilated convolutional neural network | |
CN111798462B (zh) | 一种基于ct图像的鼻咽癌放疗靶区自动勾画方法 | |
CN111161290B (zh) | 一种图像分割模型的构建方法、图像分割方法及系统 | |
CN111932559B (zh) | 基于深度学习的新冠肺炎肺部病灶区域分割系统 | |
CN103249358B (zh) | 医用图像处理装置 | |
CN111445481A (zh) | 一种基于尺度融合的腹部ct多器官的分割方法 | |
US20150029184A1 (en) | Three-dimensional model data generation device, method and program | |
CN113554669B (zh) | 一种改进注意力模块的Unet网络脑肿瘤MRI图像分割方法 | |
CN108921851A (zh) | 一种基于3d对抗网络的医学ct图像分割方法 | |
CN111340827A (zh) | 一种肺部ct影像数据处理分析方法及系统 | |
CN110874860B (zh) | 基于混合损失函数的对称监督模型的目标提取方法 | |
CN110265119A (zh) | 骨龄评估与身高预测模型、其系统及其预测方法 | |
CN107767362A (zh) | 一种基于深度学习的肺癌早筛装置 | |
CN113034507A (zh) | 基于ccta影像的冠状动脉三维分割方法 | |
CN111275706A (zh) | 一种基于剪切波弹性成像超声组学深度分析方法及系统 | |
CN111402231B (zh) | 一种用于肺部ct影像质量的自动评估系统及方法 | |
CN111311626A (zh) | 基于ct图像的颅骨骨折自动检测方法及电子介质 | |
CN115471512A (zh) | 一种基于自监督对比学习的医学影像分割方法 | |
CN113409447B (zh) | 一种基于多切片组合的冠状动脉分割方法和装置 | |
Somasundaram et al. | Brain extraction method for T1-weighted magnetic resonance scans | |
CN117523204A (zh) | 面向医疗场景的肝脏肿瘤图像分割方法、装置和可读存储介质 | |
CN111145183A (zh) | 一种用于透明隔腔超声图像的分割系统及方法 | |
CN116664590A (zh) | 基于动态对比增强磁共振图像的自动分割方法及装置 | |
CN116503548A (zh) | 牙齿和牙槽骨图像分割重构方法及系统 | |
CN114066901A (zh) | 腹部ct图像的肝脏分割方法及其ct成像方法 |
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 |