首页 > 专利 > 杭州电子科技大学 > 一种减少图像标签融合冗余计算的方法专利详情

一种减少图像标签融合冗余计算的方法   0    0

有效专利 查看PDF
专利申请流程有哪些步骤?
专利申请流程图
申请
申请号:指国家知识产权局受理一件专利申请时给予该专利申请的一个标示号码。唯一性原则。
申请日:提出专利申请之日。
2017-10-16
申请公布
申请公布指发明专利申请经初步审查合格后,自申请日(或优先权日)起18个月期满时的公布或根据申请人的请求提前进行的公布。
申请公布号:专利申请过程中,在尚未取得专利授权之前,国家专利局《专利公报》公开专利时的编号。
申请公布日:申请公开的日期,即在专利公报上予以公开的日期。
2018-04-06
授权
授权指对发明专利申请经实质审查没有发现驳回理由,授予发明专利权;或对实用新型或外观设计专利申请经初步审查没有发现驳回理由,授予实用新型专利权或外观设计专利权。
2020-09-22
预估到期
发明专利权的期限为二十年,实用新型专利权期限为十年,外观设计专利权期限为十五年,均自申请日起计算。专利届满后法律终止保护。
2037-10-16
基本信息
有效性 有效专利 专利类型 发明专利
申请号 CN201710962059.2 申请日 2017-10-16
公开/公告号 CN107798680B 公开/公告日 2020-09-22
授权日 2020-09-22 预估到期日 2037-10-16
申请年 2017年 公开/公告年 2020年
缴费截止日
分类号 G06T7/10 主分类号 G06T7/10
是否联合申请 独立申请 文献类型号 B
独权数量 1 从权数量 2
权利要求数量 3 非专利引证数量 1
引用专利数量 0 被引证专利数量 0
非专利引证 1、CN 107093190 A,2017.08.25郭其淼.多图谱医学图像分割方法研究及应用《.中国优秀硕士学位论文全文数据库 信息科技辑》.2015,(第2期),;
引用专利 被引证专利
专利权维持 5 专利申请国编码 CN
专利事件 事务标签 公开、实质审查、授权
申请人信息
申请人 第一申请人
专利权人 杭州电子科技大学 当前专利权人 杭州电子科技大学
发明人 颜成钢、赵崇宇、黄海亮、付祖贤、张腾 第一发明人 颜成钢
地址 浙江省杭州市下沙高教园区2号大街 邮编 310018
申请人数量 1 发明人数量 5
申请人所在省 浙江省 申请人所在市 浙江省杭州市
代理人信息
代理机构
专利代理机构是经省专利管理局审核,国家知识产权局批准设立,可以接受委托人的委托,在委托权限范围内以委托人的名义办理专利申请或其他专利事务的服务机构。
杭州君度专利代理事务所 代理人
专利代理师是代理他人进行专利申请和办理其他专利事务,取得一定资格的人。
朱月芬
摘要
本发明公开了一种减少图像标签融合冗余计算的方法。本发明的方法是利用GPU多线程并行计算,整体上提取目标图像和多图谱图像的灰度特征,再建立一个从目标图像的灰度特征到多图谱图像的灰度特征的映射地图,通过映射地图寻找每个点匹配到的多图谱点集和对应的标签集,接着用欧氏距离计算匹配到的点集处的权重值,然后统计可能匹配的每个标签值的权重,找到目标图像每个点处对应权重值总和最大的标签值,最后得到带标签的分割图像。本发明所述的方法能够保证在不改变原有算法的分割精度的条件下,根据GPU多线程并行计算的特点,通过改为适用于GPU多线程的计算方式,来减少标签融合时的冗余计算,提高医学图像的分割效率。
  • 摘要附图
    一种减少图像标签融合冗余计算的方法
法律状态
序号 法律状态公告日 法律状态 法律状态信息
1 2020-09-22 授权
2 2019-01-04 著录事项变更 发明人由赵崇宇 颜成钢 黄海亮 付祖贤张腾变更为颜成钢 赵崇宇 黄海亮 付祖贤张腾
3 2018-04-06 实质审查的生效 IPC(主分类): G06T 7/10 专利申请号: 201710962059.2 申请日: 2017.10.16
4 2018-03-13 公开
权利要求
权利要求书是申请文件最核心的部分,是申请人向国家申请保护他的发明创造及划定保护范围的文件。
1.一种减少图像标签融合冗余计算的方法,其特征在于包括如下步骤:
步骤(1):读取待分割目标图像的灰度图、多个模板图谱的灰度图以及相对应已分割的标签图;根据模板图谱的标签信息,建立一张标记目标图像待处理点位置信息的三维地图,并通过搜索所有目标图像待处理点附近的领域,统计并建立一张标记模板图谱的待处理点位置信息的三维地图;
步骤(2):建立两个用于存储位置信息的数组,分别存储目标图像待处理点集和模板图谱的待处理点集的三维坐标,且建立一张模板图谱待处理点三维坐标与位置数组下标序号的映射地图Ⅰ和一张模板图谱待处理点集的位置数组下标序号与对应标签值的映射地图Ⅱ;对目标图像每个待处理点,以该待处理点为中心搜索一个三维领域,并根据刚才建立的映射地图Ⅰ,将目标图像的每个待处理点与模板图谱点集进行匹配,并获取匹配后模板图谱点集对应的数组下标序号集合,得到一个从目标图像的待处理点到模板图谱点集对应的数组下标序号集合的映射关系,即映射地图Ⅲ;
步骤(3):将目标图像的灰度图、n个模板图谱、位置数组、映射地图Ⅰ、映射地图Ⅱ和映射地图Ⅲ拷贝到GPU显存中,然后在GPU中进行分割计算,GPU开启目标图像待处理点数目个线程,进行大数据处理:
根据目标图像待处理点和模板图谱待处理点的位置数组,得到一个目标图像待处理点集的所有特征向量集合和每个模板图谱待处理点集的所有特征向量集合;考虑到每个模板图谱的不同特征和分割特性,需要对每个模板图谱的映射地图Ⅲ进行匹配数据的筛选,得到一张新的目标图像待处理点与每个模板图谱下标序号集的映射地图IV;接着计算目标图像每个待处理点处对应的权重集和标签集;统计目标图像待处理点可能匹配的每个标签值的权重分布情况,将最大权重对应的标签值作为目标图像待处理点的分割结果,最终得到了目标图像的分割;
步骤1具体如下:
1-
1.从N个模板图谱中选取一个模板图谱的灰度图作为待分割目标图像;
1-
2.再从剩余的N-1个模板图谱中选取选出n个模板图谱用作模板;
1-
3.将待分割目标图像的灰度图、选出的n个模板图谱读入到计算机内存中;
1-
4.根据选出来的n个模板图谱,根据它们标签值的位置来建立一张目标图像待处理点位置信息的三维地图;在目标图像待处理点的三维地图中,考虑到搜索越界问题,对每一个待处理点,以该待处理点为中心搜索一个大小为(3,3,3)三维领域,然后统计并建立模板图谱待处理点位置信息的三维地图。

2.根据权利要求1所述的一种减少图像标签融合冗余计算的方法,其特征在于步骤2具体如下:
根据目标图像待处理点位置信息的三维地图和模板图谱的待处理点位置信息的三维地图,建立两个位置数组,分别存储目标图像待处理点集的三维坐标和模板图谱的待处理点集的三维坐标,并且建立一张模板图谱待处理点三维坐标与位置数组下标序号的映射地图Ⅰ和一张模板图谱待处理点集的位置数组下标序号与对应标签值的映射地图Ⅱ;对目标图像每个待处理点,以该点为中心搜索一个大小为(3,3,3)三维领域,并根据刚才建立的映射地图Ⅰ,得到目标图像的每个待处理点匹配模板图谱点集所对应的数组下标序号集合,即得到一张目标图像待处理点与模板图谱待处理点在数组中下标序号集的映射地图Ⅲ。

3.根据权利要求2所述的一种减少图像标签融合冗余计算的方法,其特征在于步骤3具体如下:
3-
1.根据目标图像待处理点和模板图谱待处理点的位置数组,对目标图像的每个待处理点和每个模板图谱的每个待处理点,搜索一个(5,5,5)三维领域的灰度值取为长度为125的特征向量,即得到一个目标图像待处理点集的所有特征向量集合和每个模板图谱待处理点集的所有特征向量集合;
3-
2.由于每个图谱的不同特征和分割特性,需要进行预选择筛选操作,即对每个模板图谱的映射地图Ⅲ进行匹配数据的筛选,得到一张新的目标图像待处理点与每个模板图谱下标序号集的映射地图IV;
3-
3.对目标图像所有特征向量和每个模板图谱所有特征向量进行正则化处理;
3-
4.根据映射地图IV,用欧氏距离计算两个特征之间的差异,得到目标图像每个待处理点处对应的权重集和标签集;
3-
5.对目标图像每个待处理点,统计该待处理点可能匹配的每个标签值的权重分布情况;然后在这个权重分布中,找到最大权重对应的标签值,并将该标签值作为目标图像在该待处理点处的标签值。
说明书

技术领域

[0001] 本发明属于医学图像分割技术领域,具体涉及对图像进行标签融合时,减少其中冗余计算的方法。技术背景
[0002] 随着生物成像技术的快速发展以及现代计算机存储能力和计算能力的不断提升,近年来,医学影像分析已经逐渐成为一个非常重要的研究方向。各种各样的成像技术,如计算机断层扫描(CT)、磁共振成像(MRI)等,已广泛地应用于病理检测,脑组织分析,临床诊断等各个方面,为人类的健康做出了巨大的贡献,而其中医学图像的分割,作为一项关键技术,扮演着非常重要的角色。
[0003] 图像分割就是把图像分成若干个特定的、具有独特性质的区域或者说分离出感兴趣的区域的技术和过程。在临床上,通过观察分割出来的各个功能区域图像,可以更好地来诊断或预测疾病的发生。
[0004] 近些年来,基于多图谱的图像分割方法在实际应用中取得了很大的成功。它在原有的单幅图像的基础上,发展形成了利用多幅图像配准并进行标签传播的基本框架。简单的说,图谱就是一个模板,包括图像与图像的标签。多图谱分割过程:(1)将每一幅模板图谱匹配到目标图像,并相应的将模板的标签传播到目标图像作为它的一个分割;(2)将得到的所有分割以一定的方法进行组合进而得到最终的标签。
[0005] 现有的医学图像分割方法主要集中在目标图像上每一点与多图谱上点集的映射关系。如果单单从算法本身出发,对目标图像上的一个处理点,就要进行一次多图谱的匹配点集的搜索,再用程序来实现,这是可行的,但是计算中存在很多冗余,增加了CPU的负担和分割的时间。随着CUDA并行计算框架的推出,让用GPU来处理海量数据并行计算成为了可能。因此,利用GPU来并行处理医学图像的分割将会更加高效,将减少原本程序设计中的冗余计算,提高了分割效率。

发明内容

[0006] 本发明提出一种利用GPU多线程并行计算来减少图像标签融合冗余计算的方法。该方法可以提高医学图像分割的效率。
[0007] 本发明的方法是利用GPU多线程并行计算,整体上提取目标图像和多图谱图像的灰度特征,再建立一个从目标图像的灰度特征到多图谱图像的灰度特征的映射地图,通过映射地图寻找每个点匹配到的多图谱点集和对应的标签集,接着用欧氏距离计算匹配到的点集处的权重值,然后统计可能匹配的每个标签值的权重,找到目标图像每个点处对应权重值总和最大的标签值,最后得到带标签的分割图像。
[0008] 本发明解决其技术问题所采用的技术方案具体包括以下步骤:
[0009] 步骤(1):读取待分割目标图像的灰度图、多个模板图谱的灰度图以及相对应已分割的标签图;根据模板图谱的标签信息,建立一张标记目标图像待处理点位置信息的三维地图,并通过搜索所有目标图像待处理点附近的领域,统计并建立一张标记模板图谱的待处理点位置信息的三维地图。
[0010] 具体的:
[0011] 1-1.从N个模板图谱中选取一个模板图谱的灰度图作为待分割目标图像;
[0012] 1-2.再从剩余的N-1个模板图谱中选取选出n个模板图谱用作模板;
[0013] 1-3.将待分割目标图像的灰度图、选出的n个模板图谱读入到计算机内存中;
[0014] 1-4.根据选出来的n个模板图谱,根据它们标签值的位置来建立一张目标图像待处理点位置信息的三维地图;在目标图像待处理点的三维地图中,考虑到搜索越界问题,对每一个待处理点,以该待处理点为中心搜索一个大小为(3,3,3)三维领域,然后统计并建立模板图谱待处理点位置信息的三维地图。
[0015] 步骤(2):建立两个用于存储位置信息的数组,分别存储目标图像待处理点集和模板图谱的待处理点集的三维坐标,且建立一张模板图谱待处理点三维坐标与位置数组下标序号的映射地图Ⅰ和一张模板图谱待处理点集的位置数组下标序号与对应标签值的映射地图Ⅱ;对目标图像每个待处理点,以该待处理点为中心搜索一个三维领域,并根据刚才建立的映射地图Ⅰ,将目标图像的每个待处理点与模板图谱点集进行匹配,并获取匹配后模板图谱点集对应的数组下标序号集合,得到一个从目标图像的待处理点到模板图谱点集对应的数组下标序号集合的映射关系,即映射地图Ⅲ;
[0016] 具体的:
[0017] 根据目标图像待处理点位置信息的三维地图和模板图谱的待处理点位置信息的三维地图,建立两个位置数组,分别存储目标图像待处理点集的三维坐标和模板图谱的待处理点集的三维坐标,并且建立一张模板图谱待处理点三维坐标与位置数组下标序号的映射地图Ⅰ和一张模板图谱待处理点集的位置数组下标序号与对应标签值的映射地图Ⅱ;对目标图像每个待处理点,以该点为中心搜索一个大小为(3,3,3)三维领域,并根据刚才建立的映射地图Ⅰ,得到目标图像的每个待处理点匹配模板图谱点集所对应的数组下标序号集合,即得到一张目标图像待处理点与模板图谱待处理点在数组中下标序号集的映射地图Ⅲ。
[0018] 步骤(3):将目标图像的灰度图、n个模板图谱、位置数组、映射地图Ⅰ、映射地图Ⅱ和映射地图Ⅲ拷贝到GPU显存中,然后在GPU中进行分割计算,GPU开启目标图像待处理点数目个线程,进行大数据处理:
[0019] 根据目标图像待处理点和模板图谱待处理点的位置数组,得到一个目标图像待处理点集的所有特征向量集合和每个模板图谱待处理点集的所有特征向量集合;考虑到每个模板图谱的不同特征和分割特性,需要对每个模板图谱的映射地图Ⅲ进行匹配数据的筛选,得到一张新的目标图像待处理点与每个模板图谱下标序号集的映射地图IV;接着计算目标图像每个待处理点处对应的权重集和标签集;统计目标图像待处理点可能匹配的每个标签值的权重分布情况,将最大权重对应的标签值作为目标图像待处理点的分割结果,最终得到了目标图像的分割。具体的:
[0020] 3-1.根据目标图像待处理点和模板图谱待处理点的位置数组,对目标图像的每个待处理点和每个模板图谱的每个待处理点,搜索一个(5,5,5)三维领域的灰度值取为长度为125的特征向量,即得到一个目标图像待处理点集的所有特征向量集合和每个模板图谱待处理点集的所有特征向量集合;
[0021] 3-2.由于每个图谱的不同特征和分割特性,需要进行预选择筛选操作,即对每个模板图谱的映射地图Ⅲ进行匹配数据的筛选,得到一张新的目标图像待处理点与每个模板图谱下标序号集的映射地图IV;
[0022] 3-3.对目标图像所有特征向量和每个模板图谱所有特征向量进行正则化处理;
[0023] 3-4.根据映射地图IV,用欧氏距离计算两个特征之间的差异,得到目标图像每个待处理点处对应的权重集和标签集;
[0024] 3-5.对目标图像每个待处理点,统计该待处理点可能匹配的每个标签值的权重分布情况;然后在这个权重分布中,找到最大权重对应的标签值,并将该标签值作为目标图像在该待处理点处的标签值。
[0025] 本发明的有益效果是:
[0026] 本发明所述的方法能够保证在不改变原有算法的分割精度的条件下,根据GPU多线程并行计算的特点(即能很好地用于处理大数据),通过改为适用于GPU多线程的计算方式,来减少标签融合时的冗余计算,提高医学图像的分割效率。

实施方案

[0031] 下面结合具体实施方式对本发明进行详细的说明。
[0032] 本发明提出的一种减少图像标签融合冗余计算的方法,是以标签融合算法中的非局部块加权方法(PBM)为实例来展开的,实施流程如图1所示。本发明所述方法包括以下步骤:
[0033] 步骤(1):实验中,我们有40个模板图谱(灰度图和标签图),从中选取一个图谱的灰度图作为目标灰度图像(待分割目标图像),再从剩余的39个模板图谱中选取选出10个模板图谱用作模板。读取待分割目标图像的灰度图和多个模板图谱的灰度图以及相对应已分割的标签图到计算机内存中。根据选出来的10个模板图谱,统计它们标签值的位置来建立一张目标图像待处理点位置信息的三维地图。在目标图像待处理点的三维地图中,考虑到搜索越界问题,对每一个待处理点,以该待处理点为中心搜索一个大小为(3,3,3)三维领域,然后统计并建立模板图谱待处理点位置信息的三维地图。如图2所示,这让目标图像相邻待处理点搜索到的模板图谱的待处理点集不重复,减少冗余搜索和计算。
[0034] 步骤(2):如图3和图4,根据目标图像待处理点位置信息的三维地图和模板图谱的待处理点位置信息的三维地图,建立两个位置数组,分别存储目标图像待处理点集的三维坐标和模板图谱的待处理点集的三维坐标,并且建立一张模板图谱待处理点三维坐标与位置数组下标序号的映射地图Ⅰ和一张模板图谱待处理点集的位置数组下标序号与对应标签值的映射地图Ⅱ;对目标图像每个待处理点,以该点为中心搜索一个大小为(3,3,3)三维领域,并根据刚才建立的映射地图Ⅰ,得到目标图像的每个待处理点匹配模板图谱点集所对应的数组下标序号集合,即得到一张目标图像待处理点与模板图谱待处理点在数组中下标序号集的映射地图Ⅲ。
[0035] 步骤(3):将目标图像的灰度图、10个模板图谱(灰度图和标签图)、目标图像和模板图谱待处理点集的位置数组、映射地图Ⅰ、映射地图Ⅱ和映射地图Ⅲ数据拷贝到GPU显存中。然后在GPU中进行分割计算,实例中GPU开启近150万个线程进行大数据处理:
[0036] 1)根据目标图像待处理点和模板图谱待处理点的位置数组,对目标图像的每个待处理点和每个模板图谱的每个待处理点,搜索一个(5,5,5)三维领域的灰度值取为长度为125的特征向量,即得到一个目标图像待处理点集的所有特征向量集合和每个模板图谱待处理点集的所有特征向量集合。
[0037] 2)由于每个图谱的不同特征和分割特性,需要进行预选择筛选操作,即对每个模板图谱的映射地图Ⅲ进行匹配数据的筛选,得到一张新的目标图像待处理点与每个模板图谱下标序号集的映射地图IV。
[0038] 3)对目标图像所有特征向量和每个模板图谱所有特征向量进行正则化处理。
[0039] 4)根据映射地图IV,用欧氏距离计算两个特征之间的差异,得到目标图像每个待处理点处对应的权重集和标签集。
[0040] 5)对目标图像每个待处理点,统计该待处理点可能匹配的每个标签值的权重分布情况;然后在这个权重分布中,找到最大权重对应的标签值,并将该标签值作为目标图像在该待处理点处的标签值。
[0041] GPU的每个线程处理一个目标图像待处理点,最终使得目标图像所有点都拥有标签值,即得到了目标图像的分割标签图。

附图说明

[0027] 图1是减少图像标签融合冗余计算的方法实现流程图
[0028] 图2是建立模板图谱的待处理点三维地图的过程
[0029] 图3是目标图像待处理点与目标图像的位置数组、模板图谱待处理点在数组中下标序号集的映射关系图
[0030] 图4是模板图谱待处理点与模板图谱的位置数组、标签数组、模板图谱待处理点位置-下标序号的映射关系图
版权所有:盲专网 ©2023 zlpt.xyz  蜀ICP备2023003576号