基于简化带权无向图的统计平均模型构建方法技术

本发明专利技术公开一种基于简化带权无向图的统计平均模型构建方法，能避免选择特定模板图像带来的偏差以及无模板groupwise配准收敛速度慢，效率低等问题。方法包括：将刚性对齐后的图像集作为输入，构成图模型的顶点集合，求解阈值，构造图像之间的边，建立无向图模型；在图像仿射传播聚类结果的基础上，删除类内与类间冗余的边，对无向图模型的简化；在简化图模型的基础上，根据每幅图像与周围相连图像的差异，设置边的权重；通过配准计算每幅图像与与其连接图像之间的形变场；将形变场依据对应图模型上的权重进行复合；将复合后的形变场作用于图像，使图像之间的差异逐步减小，简化带权无向图模型逐渐收缩；求取图像集的平均图像，从而获得图谱。

【技术实现步骤摘要】

本专利技术涉及图像配准
，具体涉及一种基于简化带权无向图的统计平均模型构建方法。
技术介绍
医学图像配准是指将不同条件下获取的图像变换到同一坐标空间，实现多种不同信息的融合。通过将不同模态、不同时间或者不同病人的数据配准，从而使得不同数据在对应点上达到空间的一致，为医生诊断疾病提供更加全面的信息。其在计算机辅助治疗领域拥有广泛的应用前景。随着医学图像配准技术的发展，大量精确的配准算法被提出。2002年，Shen提出HAMMER弹性配准算法，其利用边界类型、图像灰度值和几何不变矩组成属性向量对每个像素点进行属性描述，根据两幅图像像素点特征之间的差异，求解最优形变场，最终实现两幅图像之间的配准。2009年，Vercauteren提出的Diffeomorphic demons方法，它在demons配准算法的基础上，引入对称的微分同胚形变模型，从而实现精确地非参数配准。2012年，Kim提出基于机器学习的配准算法，其根据训练集中图像之间的差异，通过支持向量回归，为浮动图像学习出初始形变场，根据初始形变场，求解出中间模板图像，而中间模板图像与目标图像之间的差异较小，然后再将中间模板与目标图像配准，从而实现精确、快速地配准。上述提到的方法，其特点在于需要为每幅图像选择特定的模板，即固定图像，一经选定，在配准过程中就不会改变。该类算法在处理两幅图像之间的pairwise匹配时，可以获得较好的配准效果。但是，当需要处理大量数据，以分析它们之间的相似与差异时，选取特定模板过程中的主观因素将会导致最终配准结果的有偏性。因此，为
了解决该问题，用于分析群组数据的groupwise配准算法被提出。基于动态模板图像的groupwise配准算法的提出，试图去解决pairwise配准算法的问题。所谓动态模板图像，并不选择一个固定的图像作为模板，而是在每次配准时都会根据当前图像间相似关系，选取不同图像作为本次迭代的模板。Joshi等人提出的方法是现在常用的groupwise配准算法，该方法在每次迭代前，首先计算整体图像分布流型上的karcher均值图像，将其作为groupwise配准的模板图像，引导所有图像向其配准。Sabuncu等人则利用ICluster将图像集聚集成几个子类，将类中心图像选为该类的模板引导类内图像配准，最终所有图像配准到同一坐标空间。将均值图像作为模板，避免了偏差的产生。Sharpmean算法利用平均图像作为模板引导配准，每幅图像不再以等权重的方式构成平均图像，而是根据上次迭代的结果与图像之间的差异来决定该图像的权重。通过不同的权重，让差异较大的一些图像构成较精细的平均图像，从而引导其他图像进行配准。而对于第一次迭代来说，没有上一次迭代的结果，则选取与所有图像差异之和最小的图像作为平均图像。Shen提出了ABSORB算法，该算法的核心在于为每幅图像选取邻域图像，邻域图像即引导其配准的中间模板图像。邻域则根据该图像周围局部分布信息确立。对于无模板的groupwise配准算法，不需要预先选取模板。它们依据图像之间的相似性关系，建立配准框架，在每次迭代的过程中，通过形变场复合，逐步将图像中所有图像配准到同一坐标空间。该类groupwise配准算法中，常用的是HUGS方法。在该方法中，计算任意两幅图像之间的差异，通过线性搜索可以确定一个最优阈值。根据该阈值，建立任意两幅之间的联系(即判断阈值与两幅图像之间相似性程度的关系，以确定它们之间是否存在边)，从而构造图模型。依据图模型中边上切向力的引导(即配准)，最终将所有图像对
齐。这样一来，避免了选取模板可能带来的误差。ABSORB、HUGS和其它基于树的配准算法它们都利用了图论的相关知识，它们的核心都是利用图对图像集分布建模，根据图模型的引导，从而最终实现将大量图像配准到同一坐标空间。但是不同在于ABSORB算法的框架类似于按住模型中的一点，其余点向其移动，最终实现图模型的收缩。而HUGS建立的图中，任意两点连接在一起，且沿着边的切向力最终使所有点收缩到一个点。显而易见的是，ABSORB算法收敛的速度较快，但是结果是有偏差的。而HUGS算法避免了偏差，但是其建立的无向图收敛速度较慢，从而影响了算法的效率。
技术实现思路
有鉴于此，本专利技术实施例提供一种基于简化带权无向图的统计平均模型构建方法，能够避免选择特定模板图像带来的偏差，同时也能避免无模板groupwise配准收敛速度慢，效率低等问题。本专利技术实施例提出一种基于简化带权无向图的统计平均模型构建方法，包括：S1、将刚性配准后的图像集作为输入，构成图模型的顶点集合，求解阈值，并通过图像间相似度与阈值的比对，构造图像之间的边，从而建立无向图模型；S2、在图像集进行仿射传播聚类，在仿射传播聚类结果的基础上，根据图模型将边分为类内和类间两类，确定并删除类内与类间冗余的边，实现对无向图模型的简化，得到简化图模型；S3、在简化图模型的基础上，根据每幅图像与周围相连图像的差异，设置边的权重，该权重将用于求解图像的形变场；S4、通过配准计算每幅图像与与其连接图像之间的形变场，将形变场依据对应图模型上的权重进行复合；S5、通过将复合后的形变场作用于图像，将图像变形；S6、计算所有图像间的差异，判断所述差异是否小于预设的第一阈值，若所述差异小于所述第一阈值，则执行步骤S7，或者若所述差异不小于所述第一阈值，则将最新得到的图像作为输入重复执行步骤S1至S6，直至所有图像间的差异小于所述第一阈值，执行步骤S7；S7、求取图像集的平均图像，从而获得图谱。本专利技术实施例提供的基于简化带权无向图的统计平均模型构建方法，通过分布建模与图模型收缩，利用图模型实现快速精确的图谱求解，在求解的过程中，通过建立无向图模型对图像集分布建模，避免了在选择特定模板图像过程中的主观因素带来的偏差；综合图像集全局与局部分布信息以简化图模型，提高了groupwise配准框架的运算效率；通过设置形变场复合权重改善了无模板groupwise配准收敛速度慢等问题。附图说明图1为本专利技术基于简化带权无向图的统计平均模型构建方法一实施例的流程示意图。具体实施方式为使本专利技术实施例的目的、技术方案和优点更加清楚，下面将结合本专利技术实施例中的附图，对本专利技术实施例中的技术方案进行清楚地描述，显然，所描述的实施例是本专利技术一部分实施例，而不是全部的实施例。基于本专利技术中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本专利技术保护的范围。如图1所示，本实施例公开一种基于简化带权无向图的统计平均模型构建方法，包括：S1、将刚性配准后的图像集作为输入，构成图模型的顶点集合，求解阈值，并通过图像间相似度与阈值的比对，构造图像之间的边，
从而建立无向图模型；S2、在图像集进行仿射传播聚类，在仿射传播聚类结果的基础上，根据图模型将边分为类内和类间两类，确定并删除类内与类间冗余的边，实现对无向图模型的简化，得到简化图模型；S3、在简化图模型的基础上，根据每幅图像与周围相连图像的差异，设置边的权重，该权重将用于求解图像的形变场；S4、通过配准计算每幅图像与与其连接图像之间的形变场，将形变场依据对应图模型上的权重进行复合；S5、通过将复合后的形本文档来自技高网...

【技术保护点】
一种基于简化带权无向图的统计平均模型构建方法，其特征在于，包括：S1、将刚性配准后的图像集作为输入，构成图模型的顶点集合，求解阈值，并通过图像间相似度与阈值的比对，构造图像之间的边，从而建立无向图模型；S2、在图像集进行仿射传播聚类，在仿射传播聚类结果的基础上，根据图模型将边分为类内和类间两类，确定并删除类内与类间冗余的边，实现对无向图模型的简化，得到简化图模型；S3、在简化图模型的基础上，根据每幅图像与周围相连图像的差异，设置边的权重，该权重将用于求解图像的形变场；S4、通过配准计算每幅图像与与其连接图像之间的形变场，将形变场依据对应图模型上的权重进行复合；S5、通过将复合后的形变场作用于图像，将图像变形；S6、计算所有图像间的差异，判断所述差异是否小于预设的第一阈值，若所述差异小于所述第一阈值，则执行步骤S7，或者若所述差异不小于所述第一阈值，则将最新得到的图像作为输入重复执行步骤S1至S6，直至所有图像间的差异小于所述第一阈值，执行步骤S7；S7、求取图像集的平均图像，从而获得图谱。

【技术特征摘要】
1.一种基于简化带权无向图的统计平均模型构建方法，其特征在于，包括：S1、将刚性配准后的图像集作为输入，构成图模型的顶点集合，求解阈值，并通过图像间相似度与阈值的比对，构造图像之间的边，从而建立无向图模型；S2、在图像集进行仿射传播聚类，在仿射传播聚类结果的基础上，根据图模型将边分为类内和类间两类，确定并删除类内与类间冗余的边，实现对无向图模型的简化，得到简化图模型；S3、在简化图模型的基础上，根据每幅图像与周围相连图像的差异，设置边的权重，该权重将用于求解图像的形变场；S4、通过配准计算每幅图像与与其连接图像之间的形变场，将形变场依据对应图模型上的权重进行复合；S5、通过将复合后的形变场作用于图像，将图像变形；S6、计算所有图像间的差异，判断所述差异是否小于预设的第一阈值，若所述差异小于所述第一阈值，则执行步骤S7，或者若所述差异不小于所述第一阈值，则将最新得到的图像作为输入重复执行步骤S1至S6，直至所有图像间的差异小于所述第一阈值，执行步骤S7；S7、求取图像集的平均图像，从而获得图谱。2.根据权利要求1所述的基于简化带权无向图的统计平均模型构建方法，其特征在于，在所述...

【专利技术属性】
技术研发人员：杨健，付天宇，王涌天，李勤，范敬凡，丛伟建，
申请(专利权)人：北京理工大学，
类型：发明
国别省市：北京;11