一种CT图像半自动标注方法和系统技术方案

本发明专利技术公开一种CT图像半自动标注方法和系统，属于图像处理领域。包括：获取待标注三维CT图像，对三维CT图像去背景；对去背景后三维CT图像的每个二维切片执行二维开运算，得到各个二维切片对应的若干候选连通域；响应于针对每个二维切片选择的种子点，基于区域生长算法确定目标连通域，作为二维标注结果；叠加所有二维标注结果，得到掩膜，将预处理后三维CT图像与掩膜点乘，得到掩膜后三维CT图像；对掩膜后三维CT图像执行三维开运算，得到若干候选空间域；响应于针对掩膜后三维CT图像选择的种子点，基于区域生长算法确定目标空间域，得到三维标注结果。本发明专利技术通过形态学及区域生长，综合二维及三维信息，实现对目标快速标注，同时降低假阳性标注。降低假阳性标注。降低假阳性标注。

【技术实现步骤摘要】
一种CT图像半自动标注方法和系统


[0001]本专利技术属于医学图像处理领域，更具体地，涉及一种CT图像半自动标注方法和系统。

技术介绍

[0002]各种放射方式产生的医学图像包含诸如解剖学和病理学的大量信息，这些信息可以提供大量可视信息，并与诸如基因组学、蛋白质组学、组织病理学及临床实验数据等非图像数据构成互补。在医学图像处理过程中，经常需要对感兴趣区域(ROI)的边界给出精确标注。然而，手动标注医学图像通常耗时长且开销大，而且标注结果因人而异。在此背景下，用图像分割算法的介入来节省人工与时间成本是十分有意义的。
[0003]专利CN113192099A公开一种组织提取方法，通过响应输入装置在医学图像中选择种子点的操作，自动确定种子点对应的目标组织的形态学参数，然后根据形态学参数为种子点匹配对应的区域生长策略，并根据匹配到的区域生长策略完成目标组织区域的提取，即对于形态学参数不同的种子点采用了不同的区域生长策略。
[0004]然而，该方法存在以下缺点：该方案需要根据组织半径阈值确定区域生长策略，但是该阈值需要人为设定，泛化性较弱。同时，该方案没有利用到三维图像的空间信息，容易出现误标与漏标。

技术实现思路

[0005]针对现有技术的缺陷和改进需求，本专利技术提供了一种CT图像半自动标注方法和系统，其目的在于降低人工标注成本，减少标注过程的超参数设定，提高标注方法的繁华能力以及精度。
[0006]为实现上述目的，按照本专利技术的第一方面，提供了一种CT图像半自动标注方法，该方法包括：
[0007]S1.获取待标注的三维CT图像，对三维CT图像执行去背景操作；
[0008]S2.对去背景后三维CT图像的每个二维切片执行二维开运算，得到各个二维切片对应的若干候选连通域；
[0009]S3.响应于针对每个二维切片选择的种子点，基于区域生长算法确定目标连通域，作为二维标注结果；
[0010]S4.叠加所有二维标注结果，得到掩膜，将预处理后的三维CT图像与掩膜进行点乘，得到掩膜后的三维CT图像；
[0011]S5.对掩膜后的三维CT图像执行三维开运算，得到若干候选空间域；
[0012]S6.响应于针对掩膜后的三维CT图像选择的种子点，基于区域生长算法确定目标空间域，得到三维标注结果。
[0013]优选地，步骤S1中，所述对三维CT图像执行去背景操作，具体如下：
[0014](1)获取待标注对象的CT值范围；
[0015](2)对三维CT图像中不在所述CT值范围内的全部像素点置零。
[0016]有益效果：针对人工标注方法难以区分目标与非目标组织的问题，本专利技术通过二值化处理，基于不同组织在CT图像中的特有范围来去除背景区域，实现目标与非目标组织的初步分离。
[0017]优选地，步骤S2中，所述二维开运算的算子为圆形，其直径通过以下方式确定：
[0018](1)遍历二维切片矩阵中每一行，统计每一行中连续的非0像素长度，遍历二维切片矩阵中每一列，统计每一列中连续的非0像素长度，所有连续的非0像素长度构成一个序列；
[0019](2)保留该序列中所有小于5的长度，计算均值作为圆形形态学算子的直径。
[0020]有益效果：针对形态学开运算中算子大小难以确定的问题，本专利技术通过提取图像矩阵的连续性特征，自适应确定算子大小，在断开细小连接的同时确保目标区域的完整性。
[0021]优选地，步骤S3中，所述区域生长算法为4邻域区域生长算法或8邻域区域生长算法。
[0022]有益效果：针对人工标注方法对目标的边界难以精确到单个像素点的问题，本专利技术通过区域生长算法，对单个像素点进行分析，实现对目标边界的精确标注。
[0023]优选地，步骤S6中，所述区域生长算法包括：
[0024]S61.获取初始种子点；
[0025]S62.基于当前种子点遍历其所处二维切片周围8邻域，以及下一张二维切片中对应位置的像素点，将其中CT值不为0且不为种子点的像素点标记为新种子点；
[0026]S63.重复S62，直至没有新种子点出现，所有种子点集合为标注结果G
f
；
[0027]S64.将所有的种子点清空，保留初始种子点；
[0028]S65.基于当前种子点遍历其所处二维切片周围8邻域，以及上一张二维切片中对应位置的像素点，将其中CT值不为0且不为种子点的像素点标记为新种子点；
[0029]S66.重复S65，直至没有新种子点出现，所有种子点集合为标注结果G
b
；
[0030]S67.叠加G
f
与G
b
，得到最终的标注结果。
[0031]有益效果：针对现有26邻接三维区域生长算法，算法复杂度以及时间成本高的问题，本专利技术通过双向9邻接搜索的方法，增加算法迭代次数，降低单次遍历范围，在保证算法精度的同时，算法整体复杂度大大降低，提高召回率。
[0032]优选地，该方法还包括：
[0033]S7.对三维标注结果进行高斯平滑处理。
[0034]有益效果：针对区域生长算法结果存在细小空洞的问题，本专利技术通过高斯平滑算法，补全空洞，提高标注精度。
[0035]为实现上述目的，按照本专利技术的第二方面，提供了一种CT图像半自动标注系统，包括：计算机可读存储介质和处理器；
[0036]所述计算机可读存储介质用于存储可执行指令；
[0037]所述处理器用于读取所述计算机可读存储介质中存储的可执行指令，执行第一方面所述的CT图像半自动标注方法。
[0038]总体而言，通过本专利技术所构思的以上技术方案，能够取得以下有益效果：
[0039]针对现在人工标注方法人力与时间成本较高，半自动标注方法参数设定无法自适
应导致泛化能力弱等问题，本专利技术通过提取图像连续性特征，自适应设定形态学算子大小，同时综合二维以及三维信息，实现对目标快速地标注，同时降低假阳性标注结果并提高标注结果的召回率。
附图说明
[0040]图1为本专利技术实施例中一种面向CT图像的半自动标注方法流程图；
[0041]图2为本实施例提供的步骤S3过程示意图；
[0042]图3为本实施例提供的基于二维标注结果生成掩膜示意图；
[0043]图4为本实施例提供的双向9邻接区域生长算法的搜索方式示意图。
具体实施方式
[0044]为了使本专利技术的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本专利技术进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本专利技术，并不用于限定本专利技术。此外，下面所描述的本专利技术各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。
[0045]CT图像通过不同的CT值来表现不同的组织，同一组织具有相同的CT值范围。利用CT图像的这一特点，本专利技术基于区域生长算法提出一种高精度的半自动的标注方法，节省本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种CT图像半自动标注方法，其特征在于，该方法包括：S1.获取待标注的三维CT图像，对三维CT图像执行去背景操作；S2.对去背景后三维CT图像的每个二维切片执行二维开运算，得到各个二维切片对应的若干候选连通域；S3.响应于针对每个二维切片选择的种子点，基于区域生长算法确定目标连通域，作为二维标注结果；S4.叠加所有二维标注结果，得到掩膜，将预处理后的三维CT图像与掩膜进行点乘，得到掩膜后的三维CT图像；S5.对掩膜后的三维CT图像执行三维开运算，得到若干候选空间域；S6.响应于针对掩膜后的三维CT图像选择的种子点，基于区域生长算法确定目标空间域，得到三维标注结果。2.如权利要求1所述的方法，其特征在于，步骤S1中，所述对三维CT图像执行去背景操作，具体如下：(1)获取待标注对象的CT值范围；(2)对三维CT图像中不在所述CT值范围内的全部像素点置零。3.如权利要求1所述的方法，其特征在于，步骤S2中，所述二维开运算的算子为圆形，其直径通过以下方式确定：(1)遍历二维切片矩阵中每一行，统计每一行中连续的非0像素长度，遍历二维切片矩阵中每一列，统计每一列中连续的非0像素长度，所有连续的非0像素长度构成一个序列；(2)保留该序列中所有小于5的长度，计算均值作为圆形形态学算子的直径。4.如权利要求1所述的方法，其特征在于，步骤S...

【专利技术属性】
技术研发人员：喻莉，匡卓，
申请(专利权)人：华中科技大学，
类型：发明
国别省市：