CN115423832B - 肺动脉分割模型构建方法、肺动脉分割方法和装置 - Google Patents

肺动脉分割模型构建方法、肺动脉分割方法和装置 Download PDF

Info

Publication number
CN115423832B
CN115423832B CN202211372949.5A CN202211372949A CN115423832B CN 115423832 B CN115423832 B CN 115423832B CN 202211372949 A CN202211372949 A CN 202211372949A CN 115423832 B CN115423832 B CN 115423832B
Authority
CN
China
Prior art keywords
pulmonary artery
edge
shape
image
map
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
CN202211372949.5A
Other languages
English (en)
Other versions
CN115423832A (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.)
Zhuhai Hengqin Shengao Yunzhi Technology Co ltd
Original Assignee
Zhuhai Hengqin Shengao Yunzhi Technology Co ltd
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 Zhuhai Hengqin Shengao Yunzhi Technology Co ltd filed Critical Zhuhai Hengqin Shengao Yunzhi Technology Co ltd
Priority to CN202211372949.5A priority Critical patent/CN115423832B/zh
Publication of CN115423832A publication Critical patent/CN115423832A/zh
Application granted granted Critical
Publication of CN115423832B publication Critical patent/CN115423832B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/12Edge-based segmentation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/194Segmentation; Edge detection involving foreground-background segmentation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00Arrangements for image or video recognition or understanding
    • G06V10/70Arrangements for image or video recognition or understanding using pattern recognition or machine learning
    • G06V10/74Image or video pattern matching; Proximity measures in feature spaces
    • G06V10/75Organisation of the matching processes, e.g. simultaneous or sequential comparisons of image or video features; Coarse-fine approaches, e.g. multi-scale approaches; using context analysis; Selection of dictionaries
    • G06V10/752Contour matching
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10072Tomographic images
    • G06T2207/10081Computed x-ray tomography [CT]
    • 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)
  • Computer Vision & Pattern Recognition (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Databases & Information Systems (AREA)
  • Artificial Intelligence (AREA)
  • Computing Systems (AREA)
  • Health & Medical Sciences (AREA)
  • Evolutionary Computation (AREA)
  • General Health & Medical Sciences (AREA)
  • Medical Informatics (AREA)
  • Software Systems (AREA)
  • Multimedia (AREA)
  • Image Processing (AREA)
  • Apparatus For Radiation Diagnosis (AREA)

Abstract

本发明提供一种肺动脉分割模型构建方法、肺动脉分割方法和装置,其中构建方法包括:基于肺动脉分割模型中共享编码器的掩膜预测分支、边缘预测分支和形状预测分支分别对样本肺部影像进行预测,得到样本肺部影像的肺动脉掩膜图、肺动脉边缘概率图和肺动脉形状图;基于肺动脉形状图确定样本肺部影像的融合距离图,基于肺动脉边缘概率图确定肺动脉边缘图,并基于肺动脉边缘图与融合距离图之间的匹配程度,确定边缘形状匹配损失;基于肺动脉掩膜图与标注掩膜图间的差异、肺动脉边缘概率图与标注边缘间的差异、肺动脉形状图与标注形状间的差异以及边缘形状匹配损失,更新肺动脉分割模型的参数。本发明提升了肺动脉分割的准确性。

Description

肺动脉分割模型构建方法、肺动脉分割方法和装置
技术领域
本发明涉及图像分割技术领域,尤其涉及一种肺动脉分割模型构建方法、肺动脉分割方法和装置。
背景技术
肺动脉分割是肺部疾病相关计算机辅助诊断系统(CADs)的一项重要任务,精准的动脉分割可以提供肺部疾病的具体位置,为决定治疗方案提供帮助,因此肺动脉分割在生物医学图像处理领域得到了广泛的研究。
然而,由于拍摄CT影像时,由于血管造影剂的剂量和扫描时间对血管造影的结果有很大的影响,且肺部血管循环快,不同医生和设备扫描得到的血管造影结果在灰度值上会有比较大的区别,因此肺动脉分割任务具备较高难度。再加之肺动脉形状不规则且宽度较窄,同时肺动脉与周边相邻的部分组织距离近且影像特征十分相似,导致传统卷积神经网络在对肺动脉进行分割时,对于肺动脉边缘的分割经常产生错误的预测。因此,现有的肺动脉分割方法的准确性较低,难以满足临床应用需求。
发明内容
本发明提供一种肺动脉分割模型构建方法、肺动脉分割方法和装置,用以解决现有技术中肺动脉边缘分割不准确的缺陷。
本发明提供一种肺动脉分割模型构建方法,包括:
基于肺动脉分割模型中共享编码器的掩膜预测分支、边缘预测分支和形状预测分支分别对样本肺部影像进行预测,得到所述样本肺部影像的肺动脉掩膜图、肺动脉边缘概率图和肺动脉形状图;
基于所述肺动脉形状图确定所述样本肺部影像的融合距离图,基于所述肺动脉边缘概率图确定肺动脉边缘图,并基于所述肺动脉边缘图与所述融合距离图之间的匹配程度,确定边缘形状匹配损失;其中,任一前景像素或任一背景像素离所述肺动脉形状图的轮廓越近,所述任一前景像素或所述任一背景像素在所述融合距离图中的值越小;
基于所述肺动脉掩膜图与标注掩膜图间的差异、所述肺动脉边缘概率图与标注边缘间的差异、所述肺动脉形状图与标注形状间的差异以及所述边缘形状匹配损失,更新所述肺动脉分割模型的参数。
根据本发明提供的一种肺动脉分割模型构建方法,所述基于所述肺动脉边缘图与所述融合距离图之间的匹配程度,确定边缘形状匹配损失,具体包括:
将所述肺动脉边缘图和所述融合距离图按像素进行点乘操作,得到边缘形状匹配图;其中,所述边缘形状匹配图中任一像素的值越大,表明所述肺动脉边缘图中的所述任一像素与所述肺动脉形状图的匹配程度越低;
将所述边缘形状匹配图按像素进行累加操作,得到所述边缘形状匹配损失。
根据本发明提供的一种肺动脉分割模型构建方法,所述基于所述肺动脉形状图确定所述样本肺部影像的融合距离图,具体包括:
对所述肺动脉形状图进行欧几里得距离变换得到前景预测距离图;
将所述肺动脉形状图取反得到背景形状图后,对所述背景形状图进行欧几里得距离变换得到背景预测距离图;
分别对所述前景预测距离图和所述背景预测距离图进行sigmoid处理后相加,得到所述融合距离图。
根据本发明提供的一种肺动脉分割模型构建方法,所述基于所述肺动脉掩膜图与标注掩膜图间的差异、所述肺动脉边缘概率图与标注边缘间的差异、所述肺动脉形状图与标注形状间的差异以及所述边缘形状匹配损失,更新所述肺动脉分割模型的参数,具体包括:
基于所述肺动脉掩膜图与标注掩膜图间的差异及其第一权重、所述肺动脉边缘概率图与标注边缘间的差异及其第二权重、所述肺动脉形状图与标注形状间的差异及其第三权重,以及所述边缘形状匹配损失及其第四权重,确定模型总损失;其中,所述第一权重、所述第二权重、所述第三权重均大于所述第四权重;
基于所述模型总损失,更新所述肺动脉分割模型的参数。
根据本发明提供的一种肺动脉分割模型构建方法,所述基于肺动脉分割模型中共享编码器的掩膜预测分支、边缘预测分支和形状预测分支分别对样本肺部影像进行预测,得到所述样本肺部影像的肺动脉掩膜图、肺动脉边缘概率图和肺动脉形状图,具体包括:
基于所述肺动脉分割模型的编码器对所述样本肺部影像进行图像特征提取,得到所述样本肺部影像的图像特征编码;
分别基于所述掩膜预测分支、所述边缘预测分支和所述形状预测分支,对所述样本肺部影像的图像特征编码进行肺动脉掩膜预测、肺动脉边缘预测和肺动脉形状预测,得到所述肺动脉掩膜图、所述肺动脉边缘概率图和所述肺动脉形状图。
根据本发明提供的一种肺动脉分割模型构建方法,所述分别基于所述掩膜预测分支、所述边缘预测分支和所述形状预测分支,对所述样本肺部影像的图像特征编码进行肺动脉掩膜预测、肺动脉边缘预测和肺动脉形状预测,得到所述肺动脉掩膜图、所述肺动脉边缘概率图和所述肺动脉形状图,具体包括:
基于所述边缘预测分支和所述形状预测分支,对所述样本肺部影像的图像特征编码进行肺动脉边缘预测和肺动脉形状预测,得到所述肺动脉边缘概率图和所述肺动脉形状图;
基于所述肺动脉边缘概率图和所述肺动脉形状图,对所述样本肺部影像的图像特征编码进行优化,强化所述样本肺部影像的图像特征编码中对应肺动脉真实边缘的像素的编码值,得到所述样本肺部影像的边缘感知特征编码;
基于所述掩膜预测分支,对所述边缘感知特征编码进行肺动脉掩膜预测,得到所述肺动脉掩膜图。
本发明还提供一种肺动脉分割方法,包括:
获取待分割肺部影像;
将所述待分割肺部影像输入至肺动脉分割模型中,基于所述肺动脉分割模型的掩膜预测分支对所述待分割肺部影像进行掩膜预测,得到所述待分割肺部影像的肺动脉掩膜图,作为所述待分割肺部影像的分割结果;
其中,所述肺动脉分割模型是基于如上述任一种所述肺动脉分割模型构建方法构建得到的。
本发明还提供一种肺动脉分割模型构建装置,包括:
样本预测单元,用于基于肺动脉分割模型中共享编码器的掩膜预测分支、边缘预测分支和形状预测分支分别对样本肺部影像进行预测,得到所述样本肺部影像的肺动脉掩膜图、肺动脉边缘概率图和肺动脉形状图;
匹配损失计算单元,用于基于所述肺动脉形状图确定所述样本肺部影像的融合距离图,基于所述肺动脉边缘概率图确定肺动脉边缘图,并基于所述肺动脉边缘图与所述融合距离图之间的匹配程度,确定边缘形状匹配损失;其中,任一前景像素或任一背景像素离所述肺动脉形状图的轮廓越近,所述任一前景像素或所述任一背景像素在所述融合距离图中的值越小;
模型参数更新单元,用于基于所述肺动脉掩膜图与标注掩膜图间的差异、所述肺动脉边缘概率图与标注边缘间的差异、所述肺动脉形状图与标注形状间的差异以及所述边缘形状匹配损失,更新所述肺动脉分割模型的参数。
本发明还提供一种肺动脉分割装置,包括:
图像获取单元,用于获取待分割肺部影像;
图像分割单元,用于将所述待分割肺部影像输入至肺动脉分割模型中,基于所述肺动脉分割模型的掩膜预测分支对所述待分割肺部影像进行掩膜预测,得到所述待分割肺部影像的肺动脉掩膜图,作为所述待分割肺部影像的分割结果;
其中,所述肺动脉分割模型是基于如上述任一种所述肺动脉分割模型构建方法构建得到的。
本发明还提供一种电子设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述程序时实现如上述任一种所述肺动脉分割模型构建方法或肺动脉分割方法。
本发明还提供一种非暂态计算机可读存储介质,其上存储有计算机程序,该计算机程序被处理器执行时实现如上述任一种所述肺动脉分割模型构建方法或肺动脉分割方法。
本发明还提供一种计算机程序产品,包括计算机程序,所述计算机程序被处理器执行时实现如上述任一种所述肺动脉分割模型构建方法或肺动脉分割方法。
本发明提供的肺动脉分割模型构建方法、肺动脉分割方法和装置,通过共享编码器的掩膜预测分支、边缘预测分支和形状预测分支分别对样本肺部影像进行预测,得到肺动脉掩膜图、肺动脉边缘概率图和肺动脉形状图,基于通过肺动脉形状图确定的融合距离图与通过肺动脉边缘概率图确定的肺动脉边缘图之间的匹配程度,确定边缘形状匹配损失,用于从肺动脉形状的轮廓以外和以内两个方向约束肺动脉边缘,使其更贴合肺动脉形状的轮廓,能够提升边缘预测分支的对于肺动脉边缘的预测能力;基于肺动脉掩膜图与标注掩膜图间的差异、肺动脉边缘概率图与标注边缘间的差异、肺动脉形状图与标注形状间的差异以及边缘形状匹配损失,更新肺动脉分割模型的参数,对于掩膜预测分支,一方面通过利用肺动脉掩膜图与标注掩膜图间的差异优化其参数提升了整个分支的性能,另一方面,通过三个分支共享的编码器,在优化边缘预测分支和形状预测分支时会优化编码器的图像特征编码能力,使得编码器能够更准确地提取肺动脉的边缘特征和形状特征,使得共享该编码器的掩膜预测分支的掩膜预测能力得以增强,从而提升整个模型对于肺动脉的分割能力。
附图说明
为了更清楚地说明本发明或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作一简单地介绍,显而易见地,下面描述中的附图是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1是本发明提供的肺动脉分割模型构建方法的流程示意图;
图2是本发明提供的边缘形状匹配损失计算方法的流程示意图;
图3是本发明提供的融合距离图确定方法的流程示意图;
图4是本发明提供的掩膜边缘形状预测方法的流程示意图;
图5是本发明提供的肺动脉分割方法的流程示意图;
图6是本发明提供的肺动脉分割模型构建装置的结构示意图;
图7是本发明提供的肺动脉分割装置的结构示意图;
图8是本发明提供的电子设备的结构示意图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面将结合本发明中的附图,对本发明中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
图1是本发明提供的肺动脉分割模型构建方法的流程示意图,如图1所示,该方法包括:
步骤110,基于肺动脉分割模型中共享编码器的掩膜预测分支、边缘预测分支和形状预测分支分别对样本肺部影像进行预测,得到所述样本肺部影像的肺动脉掩膜图、肺动脉边缘概率图和肺动脉形状图;
步骤120,基于所述肺动脉形状图确定所述样本肺部影像的融合距离图,基于所述肺动脉边缘概率图确定肺动脉边缘图,并基于所述肺动脉边缘图与所述融合距离图之间的匹配程度,确定边缘形状匹配损失;其中,任一前景像素或任一背景像素离所述肺动脉形状图的轮廓越近,所述任一前景像素或所述任一背景像素在所述融合距离图中的值越小;
步骤130,基于所述肺动脉掩膜图与标注掩膜图间的差异、所述肺动脉边缘概率图与标注边缘间的差异、所述肺动脉形状图与标注形状间的差异以及所述边缘形状匹配损失,更新所述肺动脉分割模型的参数。
具体地,肺动脉分割模型中包含三个分支,分别是掩膜预测分支、边缘预测分支和形状预测分支,上述三个分支均可以基于编码器-解码器模型构建,且三个分支可以共享编码器,以提升三个分支之间的约束关系。其中,掩膜预测分支用于对输入的肺部影像进行肺动脉的掩膜分割,其主要目的在于将整个肺动脉区域分割出来,从而输出该肺部影像对应的肺动脉掩膜图,是整个肺动脉分割模型的主要分支,其输出为输入的肺部影像对应的分割结果;边缘预测分支用于对输入的肺部影像进行肺动脉的边缘分割,其主要目的在于将肺动脉的边缘分割出来,从而输出该肺部影像对应的肺动脉边缘概率图;形状预测分支则用于对输入的肺部影像进行肺动脉的形状分割,其主要目的在于获取整个肺动脉的形状特征,从而输出该肺部影像对应的肺动脉形状图。此处,边缘预测分支和形状预测分支作为掩膜预测分支的辅助分支,可以经由三个分支共享的编码器帮助掩膜预测分支将更多注意力放在边缘预测分支输出的肺动脉边缘概率图对应的边缘以内以及区域形状与形状预测分支输出的肺动脉形状图相匹配的区域范围内,从而分割出更精确的肺动脉掩膜图。从而提升肺动脉分割的准确性。
需要说明的是,部分肺动脉分割方式中也会额外增加对于肺动脉边缘的注意力以增强肺动脉分割性能,然而,肺动脉本身形状的不规则性以及宽度较窄的特性会导致其边缘难以准确分割,即使额外增加边缘分割分支,也可能会由于分割出的边缘不够准确导致对于肺动脉掩膜分割任务的帮助微乎其微、甚至起到负面作用。因此,本发明实施例在构建肺动脉分割模型时,会在掩膜预测分支和边缘预测分支的基础上增设形状预测分支,利用肺动脉形状特征较之于边缘特征更明显、提取难度较小的特点,基于形状预测分支约束边缘预测分支,再结合形状预测分支和边缘预测分支共同约束掩膜预测分支,从而提升掩膜预测分支的分割性能,分割得到更精确的肺动脉掩膜图。
具体而言,在确定完肺动脉分割模型的结构后,可以结合该模型结构利用大量的训练样本对其进行迭代训练。在当前一轮迭代训练过程中,对于任意一个样本肺部影像,利用肺动脉分割模型的掩膜预测分支、边缘预测分支和形状预测分支分别对该样本肺部影像进行肺动脉掩膜预测、肺动脉边缘预测和肺动脉形状预测,得到该样本肺部影像的肺动脉掩膜图、肺动脉边缘概率图和肺动脉形状图。其中,肺动脉掩膜图中被预测为肺动脉像素的像素点的值为1,其余像素点的值为0,肺动脉边缘概率图中各个像素点的值为该像素点被预测为肺动脉边缘的概率,肺动脉掩膜预测和肺动脉边缘预测均为分类任务。而肺动脉形状图对应着该样本肺部影像中肺动脉的形状,被形状预测分支评估为肺动脉形状内部的像素的值为1,其余像素的值为0,肺动脉形状预测为回归任务,利用回归方式能够较准确地提取肺动脉的形状特征,从而为后续的约束能力提供基础。
为了提升边缘预测分支的性能,从而增强边缘预测分支对掩膜预测分支的辅助作用,可以在边缘预测分支和形状预测分支之间建立约束。此处,可以基于该样本肺部影像的肺动脉形状图中指示的肺动脉内部区域(即前景区域)以及肺动脉以外的区域(即背景区域),确定该样本肺部影像的融合距离图。其中,融合距离图中指示了样本肺部影像中每个像素离肺动脉形状图轮廓的距离,任一前景像素或任一背景像素离肺动脉形状图的轮廓越近,该前景像素或该背景像素在融合距离图中的值越小。
基于上述融合距离图以及基于肺动脉边缘概率图确定得到的肺动脉边缘图,可以确定二者之间的匹配程度。此处,肺动脉边缘图中被边缘预测分支预测为肺动脉边缘像素的像素对应的值为1,其余像素对应的值为0。考虑到形状特征相较于边缘特征的可靠性更高,因此将以融合距离图为基准,确定肺动脉边缘图与融合距离图之间的匹配程度,肺动脉边缘图与融合距离图之间的匹配程度越高,表明肺动脉边缘图指示的肺动脉边缘与融合距离图指示的肺动脉形状的轮廓越贴合,该肺动脉边缘图的准确性更高。需要说明的是,此处通过确定肺动脉边缘图与融合距离图之间的匹配程度判断肺动脉边缘图指示的肺动脉边缘与融合距离图指示的肺动脉形状的轮廓之间的贴合程度,由于任一前景像素或任一背景像素离肺动脉形状图的轮廓越近,该前景像素或该背景像素在融合距离图中的值越小,因此可以保证从肺动脉形状的轮廓以外(即背景区域)和以内(即前景区域)两个方向约束肺动脉边缘,使其更贴合肺动脉形状的轮廓。
随后,基于肺动脉边缘图与融合距离图之间的匹配程度,确定边缘形状匹配损失。肺动脉边缘图与融合距离图之间的匹配程度越低,边缘形状匹配损失越高,在调整肺动脉分割模型参数时,该边缘形状匹配损失对于模型参数更新的影响越强,对边缘预测分支的约束越强,从而能够提升边缘预测分支的对于肺动脉边缘的预测能力。
基于肺动脉掩膜图与标注掩膜图间的差异、肺动脉边缘概率图与标注边缘间的差异、肺动脉形状图与标注形状间的差异以及上述边缘形状匹配损失,更新肺动脉分割模型的参数,以提升肺动脉分割模型的分割能力。此处,可以使用canny算子提取样本肺部影像中标注掩膜图中肺动脉的边缘,得到标注边缘;对标注掩膜图使用欧几里得距离变化得到距离图,并对距离图进行sigmoid处理,得到标注形状。
其中,由于肺动脉分割模型存在掩膜预测分支、边缘预测分支和形状预测分支等三个分支,在优化肺动脉分割模型的参数时,可以主要利用肺动脉掩膜图与标注掩膜图间的差异优化掩膜预测分支的参数,主要利用肺动脉边缘概率图与标注边缘间的差异优化边缘预测分支的参数,主要利用肺动脉形状图与标注形状间的差异优化形状预测分支的参数。此外,对于边缘预测分支,还可以额外再利用边缘形状匹配损失对该分支添加额外约束,进一步优化边缘预测分支。对于主要分支——掩膜预测分支,一方面通过利用肺动脉掩膜图与标注掩膜图间的差异优化其参数提升了整个分支的性能,另一方面,通过三个分支共享的编码器,在优化边缘预测分支和形状预测分支时会优化编码器的图像特征编码能力,使得编码器能够更准确地提取肺动脉的边缘特征和形状特征,使得共享该编码器的掩膜预测分支的掩膜预测能力得以增强,从而提升整个模型对于肺动脉的分割能力。
本发明实施例提供的方法,通过共享编码器的掩膜预测分支、边缘预测分支和形状预测分支分别对样本肺部影像进行预测,得到肺动脉掩膜图、肺动脉边缘概率图和肺动脉形状图,基于通过肺动脉形状图确定的融合距离图与通过肺动脉边缘概率图确定的肺动脉边缘图之间的匹配程度,确定边缘形状匹配损失,用于从肺动脉形状的轮廓以外和以内两个方向约束肺动脉边缘,使其更贴合肺动脉形状的轮廓,能够提升边缘预测分支的对于肺动脉边缘的预测能力;基于肺动脉掩膜图与标注掩膜图间的差异、肺动脉边缘概率图与标注边缘间的差异、肺动脉形状图与标注形状间的差异以及边缘形状匹配损失,更新肺动脉分割模型的参数,对于掩膜预测分支,一方面通过利用肺动脉掩膜图与标注掩膜图间的差异优化其参数提升了整个分支的性能,另一方面,通过三个分支共享的编码器,在优化边缘预测分支和形状预测分支时会优化编码器的图像特征编码能力,使得编码器能够更准确地提取肺动脉的边缘特征和形状特征,使得共享该编码器的掩膜预测分支的掩膜预测能力得以增强,从而提升整个模型对于肺动脉的分割能力。
基于上述实施例,如图2所示,所述基于所述肺动脉边缘图与所述融合距离图之间的匹配程度,确定边缘形状匹配损失,具体包括:
步骤210,将所述肺动脉边缘图和所述融合距离图按像素进行点乘操作,得到边缘形状匹配图;其中,所述边缘形状匹配图中任一像素的值越大,表明所述肺动脉边缘图中的所述任一像素与所述肺动脉形状图的匹配程度越低;
步骤220,将所述边缘形状匹配图按像素进行累加操作,得到所述边缘形状匹配损失。
具体地,由于融合距离图中指示了样本肺部影像中每个像素离肺动脉形状图轮廓的距离,任一前景像素或任一背景像素离肺动脉形状图的轮廓越近,该前景像素或该背景像素在融合距离图中的值越小,因此可以将肺动脉边缘图和融合距离图按像素进行点乘操作,得到边缘形状匹配图。其中,通过按像素进行点乘操作,对于肺动脉边缘图中指示的任一边缘像素,若该边缘像素不在肺动脉形状图的轮廓上,表明该边缘像素与肺动脉形状图的轮廓不贴合,且该边缘像素离肺动脉形状图的轮廓越远,该像素在融合距离图中的值越大,点乘后的值也越大,能够体现该边缘像素与肺动脉形状图的轮廓的匹配程度越低。可见,边缘形状匹配图中任一像素的值越大,表明肺动脉边缘图中的该像素与肺动脉形状图的匹配程度越低。对于非边缘像素(即被边缘预测分支预测为非边缘的像素),则点乘后的值为0,表明该像素对于后续边缘形状匹配损失的计算不做贡献。
得到边缘形状匹配图之后,将该边缘形状匹配图按像素进行累加操作,得到边缘形状匹配损失。此处,由于边缘形状匹配图中任一像素的值越大,表明肺动脉边缘图中的该像素与肺动脉形状图的匹配程度越低,因此将边缘形状匹配图按像素进行累加之后得到的值可以代表肺动脉边缘图指示的肺动脉边缘整体与肺动脉形状图的轮廓之间的贴合程度。将边缘形状匹配图按像素进行累加之后得到的累加值越高,表明肺动脉边缘图指示的肺动脉边缘整体与肺动脉形状图的轮廓之间的贴合程度越低,因此需要对边缘预测分支进行更大幅度的调整,故将该累加值作为边缘形状匹配损失,可以有效指导边缘预测分支的优化,使其预测的肺动脉边缘图更准确。
基于上述任一实施例,如图3所示,所述基于所述肺动脉形状图确定所述样本肺部影像的融合距离图,具体包括:
步骤310,对所述肺动脉形状图进行欧几里得距离变换得到前景预测距离图;
步骤320,将所述肺动脉形状图取反得到背景形状图后,对所述背景形状图进行欧几里得距离变换得到背景预测距离图;
步骤330,分别对所述前景预测距离图和所述背景预测距离图进行sigmoid处理后相加,得到所述融合距离图。
具体地,对肺动脉形状图进行欧几里得距离变换,得到相应的前景预测距离图。其中,肺动脉形状图中被形状预测分支预测为肺动脉形状内部的像素的值为1,其余像素的值为0,在对肺动脉形状图进行欧几里得距离变换后,肺动脉形状图的轮廓内部(即前景区域)的像素在前景预测距离图中对应的距离值均为0,而肺动脉形状图的轮廓外部(即背景区域)的像素在前景预测距离图中对应的距离值随着该像素与肺动脉形状图的轮廓之间距离的增大而增大。需要说明的是,在计算肺动脉形状图与标注形状间的差异时可以通过计算该前景预测距离图与标注形状之间的差异得到。
将肺动脉形状图取反,得到背景形状图。在背景形状图中被形状预测分支预测为肺动脉形状内部的像素的值为0,其余像素的值为1。类似地,对背景形状图进行欧几里得距离变换得到背景预测距离图。在对背景形状图进行欧几里得距离变换后,肺动脉形状图的轮廓外部(即背景区域)的像素在背景预测距离图中对应的距离值均为0,而肺动脉形状图的轮廓内部(即前景区域)的像素在背景预测距离图中对应的距离值随着该像素与肺动脉形状图的轮廓之间距离的增大而增大。
分别对上述前景预测距离图和上述背景预测距离图进行sigmoid处理,以弱化离肺动脉形状图的轮廓相距较远的前景像素和背景像素在后续计算边缘形状匹配损失时造成的影响,避免由于上述相距较远的前景像素和背景像素对应的距离值过大导致计算的边缘形状匹配损失过大。再将sigmoid处理后的前景预测距离图和背景预测距离图进行叠加,得到融合距离图。
基于上述任一实施例,所述基于所述肺动脉掩膜图与标注掩膜图间的差异、所述肺动脉边缘概率图与标注边缘间的差异、所述肺动脉形状图与标注形状间的差异以及所述边缘形状匹配损失,更新所述肺动脉分割模型的参数,具体包括:
基于所述肺动脉掩膜图与标注掩膜图间的差异及其第一权重、所述肺动脉边缘概率图与标注边缘间的差异及其第二权重、所述肺动脉形状图与标注形状间的差异及其第三权重,以及所述边缘形状匹配损失及其第四权重,确定模型总损失;其中,所述第一权重、所述第二权重、所述第三权重均大于所述第四权重;
基于所述模型总损失,更新所述肺动脉分割模型的参数。
具体地,在训练肺动脉分割模型时,为了提升模型优化的效率和效果,可以将肺动脉掩膜图与标注掩膜图间的差异作为掩膜预测主损失、肺动脉边缘概率图与标注边缘间的差异作为边缘预测主损失,以及将肺动脉形状图与标注形状间的差异作为形状预测主损失,以分别引导掩膜预测分支、边缘预测分支和形状预测分支等三个分支的学习。在此基础上,对于边缘预测分支,还可以额外再利用边缘形状匹配损失对该分支添加额外约束,进一步优化边缘预测分支。具体在更新肺动脉分割模型中上述三个分支的参数时,可以为掩膜预测主损失、边缘预测主损失、形状预测主损失和边缘形状匹配损失设置相应的权重,以控制各个损失在模型参数更新过程中的影响程度。
其中,在优化肺动脉分割模型的参数时,掩膜预测主损失、边缘预测主损失和形状预测主损失可以在大方向上指引掩膜预测分支、边缘预测分支和形状预测分支朝着正确的方向(即逐渐靠近标注结果),三个分支的准确性需要由掩膜预测主损失、边缘预测主损失和形状预测主损失奠定基础。在掩膜预测分支、边缘预测分支和形状预测分支具备基础准确性之后,可以基于边缘形状匹配损失对边缘预测分支和形状预测分支进行微调,其中主要会对边缘预测分支进行进一步优化使其预测得到的肺动脉边缘图与形状预测分支预测得到的肺动脉形状图的轮廓相贴合,而微调后的边缘预测分支和形状预测分支会经由三个分支共享的编码器对掩膜预测分支产生影响。
因此,在为掩膜预测主损失、边缘预测主损失、形状预测主损失和边缘形状匹配损失设置权重时,掩膜预测主损失对应的第一权重、边缘预测主损失对应的第二权重和形状预测主损失对应的第三权重均大于边缘形状匹配损失对应的第四权重。此外,第二权重还可以大于第三权重,因为形状预测分支提取的肺动脉形状特征的稳定性更高,不容易过拟合,而边缘预测分支提取的肺动脉边缘特征的稳定性较差,因此可以将第二权重设置得较高。
基于上述任一实施例,步骤110具体包括:
基于所述肺动脉分割模型的编码器对所述样本肺部影像进行图像特征提取,得到所述样本肺部影像的图像特征编码;
分别基于所述掩膜预测分支、所述边缘预测分支和所述形状预测分支,对所述样本肺部影像的图像特征编码进行肺动脉掩膜预测、肺动脉边缘预测和肺动脉形状预测,得到所述肺动脉掩膜图、所述肺动脉边缘概率图和所述肺动脉形状图。
具体地,可以基于掩膜预测分支、边缘预测分支和形状预测分支共享的编码器对样本肺部影像进行图像特征提取,获取样本肺部影像中与肺动脉区域存在关联的图像特征,得到样本肺部影像的图像特征编码。其中,图像特征编码中包含对应影像中的肺动脉区域特征以及与肺动脉区域存在关联的其余区域的特征。
分别基于肺动脉分割模型的掩膜预测分支、边缘预测分支和形状预测分支,对样本肺部影像的图像特征编码进行肺动脉掩膜预测、肺动脉边缘预测和肺动脉形状预测,得到肺动脉掩膜图、肺动脉边缘概率图和肺动脉形状图。其中,掩膜预测分支和边缘预测分支的解码器可以基于分类方式分别进行肺动脉掩膜预测和肺动脉边缘预测,得到肺动脉掩膜图和肺动脉边缘概率图,而形状预测分支的解码器可以基于回归方式进行肺动脉形状预测,得到肺动脉形状图。
基于上述任一实施例,如图4所示,所述掩膜预测分支、所述边缘预测分支和所述形状预测分支,对所述样本肺部影像的图像特征编码进行肺动脉掩膜预测、肺动脉边缘预测和肺动脉形状预测,得到所述肺动脉掩膜图、所述肺动脉边缘概率图和所述肺动脉形状图,具体包括:
步骤410,基于所述边缘预测分支和所述形状预测分支,对所述样本肺部影像的图像特征编码进行肺动脉边缘预测和肺动脉形状预测,得到所述肺动脉边缘概率图和所述肺动脉形状图;
步骤420,基于所述肺动脉边缘概率图和所述肺动脉形状图,对所述样本肺部影像的图像特征编码进行优化,强化所述样本肺部影像的图像特征编码中对应肺动脉真实边缘的像素的编码值,得到所述样本肺部影像的边缘感知特征编码;
步骤430,基于所述掩膜预测分支,对所述边缘感知特征编码进行肺动脉掩膜预测,得到所述肺动脉掩膜图。
具体地,为了在掩膜预测分支、边缘预测分支和形状预测分支之间添加更多的约束,以使得边缘预测分支和形状预测分支可以为掩膜预测分支提供更多辅助信息,包括肺动脉的边缘信息和肺动脉的形状信息等,以提升掩膜预测分支的掩膜预测性能,可以在掩膜预测分支中的编码器之后添加编码优化模块,以提升掩膜预测分支的解码效果。
具体而言,对于边缘预测分支和形状预测分支,可以对样本肺部影像的图像特征编码进行肺动脉边缘预测和肺动脉形状预测,得到肺动脉边缘概率图和肺动脉形状图。进一步地,可以首先利用形状预测分支对样本肺部影像的图像特征编码进行肺动脉形状预测,得到肺动脉形状图,再基于注意力机制,利用肺动脉形状图对图像特征编码进行优化,基于肺动脉形状图对图像特征编码进行注意力变换,将图像特征编码中与肺动脉形状图相匹配的区域的编码值强化,得到形状感知编码。再基于边缘预测分支对形状感知编码进行肺动脉边缘预测,得到肺动脉边缘概率图,以提升肺动脉边缘概率图的准确性。
随后,对于掩膜预测分支,可以基于边缘预测分支和形状预测分支输出的肺动脉边缘概率图和肺动脉形状图,对样本肺部影像的图像特征编码进行优化,强化样本肺部影像的图像特征编码中对应肺动脉真实边缘的像素的编码值,得到样本肺部影像的边缘感知特征编码。其中,可以基于肺动脉边缘概率图和肺动脉形状图之间的相似性对肺动脉边缘概率图进行优化,强化其中与肺动脉形状图相匹配的像素的概率值,得到优化边缘概率图。再基于优化边缘概率图对样本肺部影像的图像特征编码进行优化,强化样本肺部影像的图像特征编码中对应肺动脉真实边缘的像素的编码值。此处,可以利用注意力机制,基于优化边缘概率图对图像特征编码进行注意力变换,将图像特征编码中对应肺动脉真实边缘的像素的编码值强化,得到样本肺部影像的边缘感知特征编码。再基于掩膜预测分支的解码器对边缘感知特征编码进行肺动脉掩膜预测,得到肺动脉掩膜图,以提升肺动脉掩膜图的准确性。
基于上述任一实施例,图5是本发明提供的肺动脉分割方法的流程示意图,如图5所示,该方法包括:
步骤510,获取待分割肺部影像;
步骤520,将所述待分割肺部影像输入至肺动脉分割模型中,基于所述肺动脉分割模型的掩膜预测分支对所述待分割肺部影像进行掩膜预测,得到所述待分割肺部影像的肺动脉掩膜图,作为所述待分割肺部影像的分割结果;
其中,所述肺动脉分割模型是基于如上述任一实施例提供的肺动脉分割模型构建方法构建得到的。
具体地,本实施例中采取窗宽1600、窗位-600以采集肺部CT影像。由于胸腔部分进入模型很可能会引发错误预测,故可以使用图像处理方法提取肺区域掩膜,根据肺区域掩膜生成的目标框截取肺部CT影像的单一肺部影像,随后可以将影像尺寸缩放至[128,128,128],降低内存成本。此外还可以通过调整其灰度范围以削弱影像中一些其它组织带来的影响,将CT值归一化到灰度范围(0~255),得到待分割肺部影像。此处,待分割肺部影像的处理方式与样本肺部影像的处理方式相同,即将影像输入至肺动脉分割模型之前可以进行相同的图像处理操作,以保证模型输入的统一性。
随后,将待分割肺部影像输入至肺动脉分割模型中,肺动脉分割模型会自主对待分割肺部影像进行编解码操作,得到相应的输出。此处,肺动脉分割模型中的掩膜预测分支会对待分割肺部影像进行编码解码操作,实现肺动脉掩膜预测,得到待分割肺部影像的肺动脉掩膜图,作为该待分割肺部影像的分割结果。其中,肺动脉分割模型是基于如上述任一实施例提供的肺动脉分割模型构建方法构建得到的。具体而言,可以基于肺动脉分割模型的编码器对待分割肺部影像进行图像特征提取,得到待分割肺部影像的图像特征编码,再基于掩膜预测分支的解码器,对待分割肺部影像的图像特征编码进行肺动脉掩膜预测,得到待分割肺部影像的肺动脉掩膜图,实现准确的肺动脉分割。
进一步地,还可以如上述实施例中提供的方式,基于边缘预测分支和形状预测分支,对待分割肺部影像的图像特征编码进行肺动脉边缘预测和肺动脉形状预测,得到待分割肺部影像的肺动脉边缘概率图和肺动脉形状图。基于待分割肺部影像的肺动脉边缘概率图和肺动脉形状图,对待分割肺部影像的图像特征编码进行优化,强化待分割肺部影像的图像特征编码中对应肺动脉真实边缘的像素的编码值,得到待分割肺部影像的边缘感知特征编码。再基于掩膜预测分支,对待分割肺部影像的边缘感知特征编码进行肺动脉掩膜预测,得到待分割肺部影像的肺动脉掩膜图,实现更精确的肺动脉分割。
下面对本发明提供的肺动脉分割模型构建装置进行描述,下文描述的肺动脉分割模型构建装置与上文描述的肺动脉分割模型构建方法可相互对应参照。
基于上述任一实施例,图6是本发明提供的肺动脉分割模型构建装置的结构示意图,如图6所示,该装置包括:样本预测单元610、匹配损失计算单元620和模型参数更新单元630。
其中,样本预测单元610用于基于肺动脉分割模型中共享编码器的掩膜预测分支、边缘预测分支和形状预测分支分别对样本肺部影像进行预测,得到所述样本肺部影像的肺动脉掩膜图、肺动脉边缘概率图和肺动脉形状图;
匹配损失计算单元620用于基于所述肺动脉形状图确定所述样本肺部影像的融合距离图,基于所述肺动脉边缘概率图确定肺动脉边缘图,并基于所述肺动脉边缘图与所述融合距离图之间的匹配程度,确定边缘形状匹配损失;其中,任一前景像素或任一背景像素离所述肺动脉形状图的轮廓越近,所述任一前景像素或所述任一背景像素在所述融合距离图中的值越小;
模型参数更新单元630用于基于所述肺动脉掩膜图与标注掩膜图间的差异、所述肺动脉边缘概率图与标注边缘间的差异、所述肺动脉形状图与标注形状间的差异以及所述边缘形状匹配损失,更新所述肺动脉分割模型的参数。
本发明实施例提供的装置,通过共享编码器的掩膜预测分支、边缘预测分支和形状预测分支分别对样本肺部影像进行预测,得到肺动脉掩膜图、肺动脉边缘概率图和肺动脉形状图,基于通过肺动脉形状图确定的融合距离图与通过肺动脉边缘概率图确定的肺动脉边缘图之间的匹配程度,确定边缘形状匹配损失,用于从肺动脉形状的轮廓以外和以内两个方向约束肺动脉边缘,使其更贴合肺动脉形状的轮廓,能够提升边缘预测分支的对于肺动脉边缘的预测能力;基于肺动脉掩膜图与标注掩膜图间的差异、肺动脉边缘概率图与标注边缘间的差异、肺动脉形状图与标注形状间的差异以及边缘形状匹配损失,更新肺动脉分割模型的参数,对于掩膜预测分支,一方面通过利用肺动脉掩膜图与标注掩膜图间的差异优化其参数提升了整个分支的性能,另一方面,通过三个分支共享的编码器,在优化边缘预测分支和形状预测分支时会优化编码器的图像特征编码能力,使得编码器能够更准确地提取肺动脉的边缘特征和形状特征,使得共享该编码器的掩膜预测分支的掩膜预测能力得以增强,从而提升整个模型对于肺动脉的分割能力。
下面对本发明提供的肺动脉分割装置进行描述,下文描述的肺动脉分割装置与上文描述的肺动脉分割方法可相互对应参照。
基于上述任一实施例,图7是本发明提供的肺动脉分割装置的结构示意图,如图7所示,该装置包括:图像获取单元710和图像分割单元720。
其中,图像获取单元710用于获取待分割肺部影像;
图像分割单元720用于将所述待分割肺部影像输入至肺动脉分割模型中,基于所述肺动脉分割模型的掩膜预测分支对所述待分割肺部影像进行掩膜预测,得到所述待分割肺部影像的肺动脉掩膜图,作为所述待分割肺部影像的分割结果;
其中,所述肺动脉分割模型是基于如上述任一实施例提供的肺动脉分割模型构建方法构建得到的。
图8是本发明提供的电子设备的结构示意图,如图8所示,该电子设备可以包括:处理器(processor)810、存储器(memory)820、通信接口(Communications Interface)830和通信总线840,其中,处理器810,存储器820,通信接口830通过通信总线840完成相互间的通信。处理器810可以调用存储器820中的逻辑指令,以执行肺动脉分割模型构建方法,该方法包括:基于肺动脉分割模型中共享编码器的掩膜预测分支、边缘预测分支和形状预测分支分别对样本肺部影像进行预测,得到所述样本肺部影像的肺动脉掩膜图、肺动脉边缘概率图和肺动脉形状图;基于所述肺动脉形状图确定所述样本肺部影像的融合距离图,基于所述肺动脉边缘概率图确定肺动脉边缘图,并基于所述肺动脉边缘图与所述融合距离图之间的匹配程度,确定边缘形状匹配损失;其中,任一前景像素或任一背景像素离所述肺动脉形状图的轮廓越近,所述任一前景像素或所述任一背景像素在所述融合距离图中的值越小;基于所述肺动脉掩膜图与标注掩膜图间的差异、所述肺动脉边缘概率图与标注边缘间的差异、所述肺动脉形状图与标注形状间的差异以及所述边缘形状匹配损失,更新所述肺动脉分割模型的参数。
处理器810还可以调用存储器820中的逻辑指令,以执行肺动脉分割方法,该方法包括:获取待分割肺部影像;将所述待分割肺部影像输入至肺动脉分割模型中,基于所述肺动脉分割模型的掩膜预测分支对所述待分割肺部影像进行掩膜预测,得到所述待分割肺部影像的肺动脉掩膜图,作为所述待分割肺部影像的分割结果;其中,所述肺动脉分割模型是基于如上述任一实施例提供的肺动脉分割模型构建方法构建得到的。
此外,上述的存储器820中的逻辑指令可以通过软件功能单元的形式实现并作为独立的产品销售或使用时,可以存储在一个计算机可读取存储介质中。基于这样的理解,本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来,该计算机软件产品存储在一个存储介质中,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质包括:U盘、移动硬盘、只读存储器(ROM,Read-Only Memory)、随机存取存储器(RAM,Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。
另一方面,本发明还提供一种计算机程序产品,所述计算机程序产品包括存储在非暂态计算机可读存储介质上的计算机程序,所述计算机程序包括程序指令,当所述程序指令被计算机执行时,计算机能够执行上述各方法所提供的肺动脉分割模型构建方法,该方法包括:基于肺动脉分割模型中共享编码器的掩膜预测分支、边缘预测分支和形状预测分支分别对样本肺部影像进行预测,得到所述样本肺部影像的肺动脉掩膜图、肺动脉边缘概率图和肺动脉形状图;基于所述肺动脉形状图确定所述样本肺部影像的融合距离图,基于所述肺动脉边缘概率图确定肺动脉边缘图,并基于所述肺动脉边缘图与所述融合距离图之间的匹配程度,确定边缘形状匹配损失;其中,任一前景像素或任一背景像素离所述肺动脉形状图的轮廓越近,所述任一前景像素或所述任一背景像素在所述融合距离图中的值越小;基于所述肺动脉掩膜图与标注掩膜图间的差异、所述肺动脉边缘概率图与标注边缘间的差异、所述肺动脉形状图与标注形状间的差异以及所述边缘形状匹配损失,更新所述肺动脉分割模型的参数。当所述程序指令被计算机执行时,计算机还能够执行上述各方法所提供的肺动脉分割方法,该方法包括:获取待分割肺部影像;将所述待分割肺部影像输入至肺动脉分割模型中,基于所述肺动脉分割模型的掩膜预测分支对所述待分割肺部影像进行掩膜预测,得到所述待分割肺部影像的肺动脉掩膜图,作为所述待分割肺部影像的分割结果;其中,所述肺动脉分割模型是基于如上述任一实施例提供的肺动脉分割模型构建方法构建得到的。
又一方面,本发明还提供一种非暂态计算机可读存储介质,其上存储有计算机程序,该计算机程序被处理器执行时实现以执行上述各提供的肺动脉分割模型构建方法,该方法包括:基于肺动脉分割模型中共享编码器的掩膜预测分支、边缘预测分支和形状预测分支分别对样本肺部影像进行预测,得到所述样本肺部影像的肺动脉掩膜图、肺动脉边缘概率图和肺动脉形状图;基于所述肺动脉形状图确定所述样本肺部影像的融合距离图,基于所述肺动脉边缘概率图确定肺动脉边缘图,并基于所述肺动脉边缘图与所述融合距离图之间的匹配程度,确定边缘形状匹配损失;其中,任一前景像素或任一背景像素离所述肺动脉形状图的轮廓越近,所述任一前景像素或所述任一背景像素在所述融合距离图中的值越小;基于所述肺动脉掩膜图与标注掩膜图间的差异、所述肺动脉边缘概率图与标注边缘间的差异、所述肺动脉形状图与标注形状间的差异以及所述边缘形状匹配损失,更新所述肺动脉分割模型的参数。该计算机程序被处理器执行时还可以实现以执行上述各提供的肺动脉分割方法,该方法包括:获取待分割肺部影像;将所述待分割肺部影像输入至肺动脉分割模型中,基于所述肺动脉分割模型的掩膜预测分支对所述待分割肺部影像进行掩膜预测,得到所述待分割肺部影像的肺动脉掩膜图,作为所述待分割肺部影像的分割结果;其中,所述肺动脉分割模型是基于如上述任一实施例提供的肺动脉分割模型构建方法构建得到的。
以上所描述的装置实施例仅仅是示意性的,其中所述作为分离部件说明的单元可以是或者也可以不是物理上分开的,作为单元显示的部件可以是或者也可以不是物理单元,即可以位于一个地方,或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部模块来实现本实施例方案的目的。本领域普通技术人员在不付出创造性的劳动的情况下,即可以理解并实施。
通过以上的实施方式的描述,本领域的技术人员可以清楚地了解到各实施方式可借助软件加必需的通用硬件平台的方式来实现,当然也可以通过硬件。基于这样的理解,上述技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来,该计算机软件产品可以存储在计算机可读存储介质中,如ROM/RAM、磁碟、光盘等,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行各个实施例或者实施例的某些部分所述的方法。
最后应说明的是:以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。

Claims (7)

1.一种肺动脉分割模型构建方法,其特征在于,包括:
基于肺动脉分割模型中共享编码器的掩膜预测分支、边缘预测分支和形状预测分支分别对样本肺部影像进行预测,得到所述样本肺部影像的肺动脉掩膜图、肺动脉边缘概率图和肺动脉形状图;其中,所述掩膜预测分支、边缘预测分支和形状预测分支均基于编码器-解码器模型构建得到;
基于所述肺动脉形状图确定所述样本肺部影像的融合距离图,基于所述肺动脉边缘概率图确定肺动脉边缘图,并基于所述肺动脉边缘图与所述融合距离图之间的匹配程度,确定边缘形状匹配损失;其中,任一前景像素或任一背景像素离所述肺动脉形状图的轮廓越近,所述任一前景像素或所述任一背景像素在所述融合距离图中的值越小;
基于所述肺动脉掩膜图与标注掩膜图间的差异、所述肺动脉边缘概率图与标注边缘间的差异、所述肺动脉形状图与标注形状间的差异以及所述边缘形状匹配损失,更新所述肺动脉分割模型的参数;
所述基于所述肺动脉边缘图与所述融合距离图之间的匹配程度,确定边缘形状匹配损失,具体包括:
将所述肺动脉边缘图和所述融合距离图按像素进行点乘操作,得到边缘形状匹配图;其中,所述边缘形状匹配图中任一像素的值越大,表明所述肺动脉边缘图中的所述任一像素与所述肺动脉形状图的匹配程度越低;
将所述边缘形状匹配图按像素进行累加操作,得到所述边缘形状匹配损失;
所述基于肺动脉分割模型中共享编码器的掩膜预测分支、边缘预测分支和形状预测分支分别对样本肺部影像进行预测,得到所述样本肺部影像的肺动脉掩膜图、肺动脉边缘概率图和肺动脉形状图,具体包括:
基于所述肺动脉分割模型的编码器对所述样本肺部影像进行图像特征提取,得到所述样本肺部影像的图像特征编码;
基于所述边缘预测分支和所述形状预测分支,对所述样本肺部影像的图像特征编码进行肺动脉边缘预测和肺动脉形状预测,得到所述肺动脉边缘概率图和所述肺动脉形状图;
基于所述肺动脉边缘概率图和所述肺动脉形状图之间的相似性对所述肺动脉边缘概率图进行优化,强化所述肺动脉边缘概率图中与所述肺动脉形状图相匹配的像素的概率值,得到优化边缘概率图;利用注意力机制,基于所述优化边缘概率图对所述样本肺部影像的图像特征编码进行注意力变换,将所述样本肺部影像的图像特征编码中对应肺动脉真实边缘的像素的编码值强化,得到所述样本肺部影像的边缘感知特征编码;
基于所述掩膜预测分支,对所述边缘感知特征编码进行肺动脉掩膜预测,得到所述肺动脉掩膜图。
2.根据权利要求1所述的肺动脉分割模型构建方法,其特征在于,所述基于所述肺动脉形状图确定所述样本肺部影像的融合距离图,具体包括:
对所述肺动脉形状图进行欧几里得距离变换得到前景预测距离图;
将所述肺动脉形状图取反得到背景形状图后,对所述背景形状图进行欧几里得距离变换得到背景预测距离图;
分别对所述前景预测距离图和所述背景预测距离图进行sigmoid处理后相加,得到所述融合距离图。
3.根据权利要求1所述的肺动脉分割模型构建方法,其特征在于,所述基于所述肺动脉掩膜图与标注掩膜图间的差异、所述肺动脉边缘概率图与标注边缘间的差异、所述肺动脉形状图与标注形状间的差异以及所述边缘形状匹配损失,更新所述肺动脉分割模型的参数,具体包括:
基于所述肺动脉掩膜图与标注掩膜图间的差异及其第一权重、所述肺动脉边缘概率图与标注边缘间的差异及其第二权重、所述肺动脉形状图与标注形状间的差异及其第三权重,以及所述边缘形状匹配损失及其第四权重,确定模型总损失;其中,所述第一权重、所述第二权重、所述第三权重均大于所述第四权重;
基于所述模型总损失,更新所述肺动脉分割模型的参数。
4.一种肺动脉分割方法,其特征在于,包括:
获取待分割肺部影像;
将所述待分割肺部影像输入至肺动脉分割模型中,基于所述肺动脉分割模型的掩膜预测分支对所述待分割肺部影像进行掩膜预测,得到所述待分割肺部影像的肺动脉掩膜图,作为所述待分割肺部影像的分割结果;
其中,所述肺动脉分割模型是基于如权利要求1至3任一项所述肺动脉分割模型构建方法构建得到的。
5.一种肺动脉分割模型构建装置,其特征在于,包括:
样本预测单元,用于基于肺动脉分割模型中共享编码器的掩膜预测分支、边缘预测分支和形状预测分支分别对样本肺部影像进行预测,得到所述样本肺部影像的肺动脉掩膜图、肺动脉边缘概率图和肺动脉形状图;其中,所述掩膜预测分支、边缘预测分支和形状预测分支均基于编码器-解码器模型构建得到;
匹配损失计算单元,用于基于所述肺动脉形状图确定所述样本肺部影像的融合距离图,基于所述肺动脉边缘概率图确定肺动脉边缘图,并基于所述肺动脉边缘图与所述融合距离图之间的匹配程度,确定边缘形状匹配损失;其中,任一前景像素或任一背景像素离所述肺动脉形状图的轮廓越近,所述任一前景像素或所述任一背景像素在所述融合距离图中的值越小;
模型参数更新单元,用于基于所述肺动脉掩膜图与标注掩膜图间的差异、所述肺动脉边缘概率图与标注边缘间的差异、所述肺动脉形状图与标注形状间的差异以及所述边缘形状匹配损失,更新所述肺动脉分割模型的参数;
所述基于所述肺动脉边缘图与所述融合距离图之间的匹配程度,确定边缘形状匹配损失,具体包括:
将所述肺动脉边缘图和所述融合距离图按像素进行点乘操作,得到边缘形状匹配图;其中,所述边缘形状匹配图中任一像素的值越大,表明所述肺动脉边缘图中的所述任一像素与所述肺动脉形状图的匹配程度越低;
将所述边缘形状匹配图按像素进行累加操作,得到所述边缘形状匹配损失;
所述基于肺动脉分割模型中共享编码器的掩膜预测分支、边缘预测分支和形状预测分支分别对样本肺部影像进行预测,得到所述样本肺部影像的肺动脉掩膜图、肺动脉边缘概率图和肺动脉形状图,具体包括:
基于所述肺动脉分割模型的编码器对所述样本肺部影像进行图像特征提取,得到所述样本肺部影像的图像特征编码;
基于所述边缘预测分支和所述形状预测分支,对所述样本肺部影像的图像特征编码进行肺动脉边缘预测和肺动脉形状预测,得到所述肺动脉边缘概率图和所述肺动脉形状图;
基于所述肺动脉边缘概率图和所述肺动脉形状图之间的相似性对所述肺动脉边缘概率图进行优化,强化所述肺动脉边缘概率图中与所述肺动脉形状图相匹配的像素的概率值,得到优化边缘概率图;利用注意力机制,基于所述优化边缘概率图对所述样本肺部影像的图像特征编码进行注意力变换,将所述样本肺部影像的图像特征编码中对应肺动脉真实边缘的像素的编码值强化,得到所述样本肺部影像的边缘感知特征编码;
基于所述掩膜预测分支,对所述边缘感知特征编码进行肺动脉掩膜预测,得到所述肺动脉掩膜图。
6.一种肺动脉分割装置,其特征在于,包括:
图像获取单元,用于获取待分割肺部影像;
图像分割单元,用于将所述待分割肺部影像输入至肺动脉分割模型中,基于所述肺动脉分割模型的掩膜预测分支对所述待分割肺部影像进行掩膜预测,得到所述待分割肺部影像的肺动脉掩膜图,作为所述待分割肺部影像的分割结果;
其中,所述肺动脉分割模型是基于如权利要求1至3任一项所述肺动脉分割模型构建方法构建得到的。
7.一种非暂态计算机可读存储介质,其上存储有计算机程序,其特征在于,所述计算机程序被处理器执行时实现如权利要求1至3任一项所述肺动脉分割模型构建方法或如权利要求4所述肺动脉分割方法。
CN202211372949.5A 2022-11-04 2022-11-04 肺动脉分割模型构建方法、肺动脉分割方法和装置 Active CN115423832B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211372949.5A CN115423832B (zh) 2022-11-04 2022-11-04 肺动脉分割模型构建方法、肺动脉分割方法和装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211372949.5A CN115423832B (zh) 2022-11-04 2022-11-04 肺动脉分割模型构建方法、肺动脉分割方法和装置

Publications (2)

Publication Number Publication Date
CN115423832A CN115423832A (zh) 2022-12-02
CN115423832B true CN115423832B (zh) 2023-03-03

Family

ID=84207548

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211372949.5A Active CN115423832B (zh) 2022-11-04 2022-11-04 肺动脉分割模型构建方法、肺动脉分割方法和装置

Country Status (1)

Country Link
CN (1) CN115423832B (zh)

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2020080913A (ja) * 2018-11-15 2020-06-04 株式会社医用科学研究所 非造影ct画像からの3次元メディアルアクシスモデルに基づく関心臓器画像自動セグメンテーション装置、及び自動セグメンテーション方法
CN114596275A (zh) * 2022-03-02 2022-06-07 推想医疗科技股份有限公司 肺血管分割方法、装置、存储介质及电子设备

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10157467B2 (en) * 2015-08-07 2018-12-18 Arizona Board Of Regents On Behalf Of Arizona State University System and method for detecting central pulmonary embolism in CT pulmonary angiography images
CN113569790B (zh) * 2019-07-30 2022-07-29 北京市商汤科技开发有限公司 图像处理方法及装置、处理器、电子设备及存储介质
CN114998611A (zh) * 2022-05-27 2022-09-02 电子科技大学 一种基于结构融合的目标轮廓检测方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2020080913A (ja) * 2018-11-15 2020-06-04 株式会社医用科学研究所 非造影ct画像からの3次元メディアルアクシスモデルに基づく関心臓器画像自動セグメンテーション装置、及び自動セグメンテーション方法
CN114596275A (zh) * 2022-03-02 2022-06-07 推想医疗科技股份有限公司 肺血管分割方法、装置、存储介质及电子设备

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
"A Multi-Task Network with Distance–Mask–Boundary Consistency Constraints for Building Extraction from Aerial Image";Shi, F.;《Remote Sensing》;20210706;第1-22页 *
"Lung vessel segmentation based on random forests";Zhao, B.;《Electronics Letters》;20170119;第 220-222页 *
2D/3D级联卷积在分割CT肺动脉上的应用研究;黄绍辉等;《中国数字医学》;20190515(第05期);第12-16页 *

Also Published As

Publication number Publication date
CN115423832A (zh) 2022-12-02

Similar Documents

Publication Publication Date Title
CN110310287B (zh) 基于神经网络的危及器官自动勾画方法、设备和存储介质
CN115661144B (zh) 基于可变形U-Net的自适应医学图像分割方法
CN112465834B (zh) 一种血管分割方法及装置
CN110136135B (zh) 分割方法、装置、设备以及存储介质
WO2021136368A1 (zh) 钼靶图像中胸大肌区域自动检测方法及装置
CN111899244B (zh) 图像分割、网络模型的训练方法及装置,及电子设备
CN112396605B (zh) 网络训练方法及装置、图像识别方法和电子设备
CN113554665A (zh) 一种血管分割方法及装置
CN116245832B (zh) 一种图像处理方法、装置、设备及存储介质
US8340378B2 (en) Ribcage segmentation
CN113256670A (zh) 图像处理方法及装置、网络模型的训练方法及装置
CN113870284A (zh) 用于对医学影像进行分割的方法、设备和介质
CN115546231A (zh) 一种基于半监督深度学习的自适应脑胶质瘤分割方法
CN111626379A (zh) 肺炎x光图像检测方法
CN113888566B (zh) 目标轮廓曲线确定方法、装置、电子设备以及存储介质
CN115423832B (zh) 肺动脉分割模型构建方法、肺动脉分割方法和装置
US20230334698A1 (en) Methods and systems for positioning in an medical procedure
CN113780040A (zh) 唇部关键点的定位方法及装置、存储介质、电子设备
CN117333428A (zh) 一种ct影像分割方法、装置、设备及存储介质
CN115439423B (zh) 一种基于ct图像的识别方法、装置、设备及存储介质
CN112116608B (zh) 一种导丝分割方法、装置、电子设备及存储介质
CN115861255A (zh) 用于图像处理的模型训练方法、装置、设备、介质及产品
CN114419375A (zh) 图像分类方法、训练方法、装置、电子设备以及存储介质
CN114445419A (zh) 一种基于支气管拓扑结构的肺段分割方法、装置及系统
CN115564788B (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