本发明专利技术涉及一种基于粒子群算法的自适应脊柱CT图像分割方法，该方法包括如下步骤：(1)输入原始图像；(2)初始化第一代种群；(3)计算个体最优和全局最优；(4)产生新的个体；(5)计算新的个体最优和全局最优；(6)判断是否满足最大迭代次数，如果满足则转到步骤(7)，否则，返回步骤(4)；(7)进行初步的图像分割，作为初始的分割结果；(8)采取拓扑学操作，对脊柱进一步准确的分割，输出图像。本发明专利技术优化了搜索粒度，由粗到细控制算法收敛的速度，大大加快了算法的收敛速度，提高了分割精度，并利用待分割区域的先验知识，采用拓扑学对脊柱进一步进行准确的分割，大大提升了图像分割的效率和质量。
【技术实现步骤摘要】
本专利技术涉及一般的图像处理，具体涉及图像分析，尤其是涉及通过图像属性的分 析来分割图像。
技术介绍
随着计算机科学技术的不断发展，图像处理技术也得到了一定的发展，逐步形成 了自己的科学理论体系。图像分割是图像处理的重要技术，在实际应用和理论研究中受到 越来越多人的更加关注。图像分割技术是把图像划分为若干个有意义的区域，把感兴趣的目标提取出来， 以便后续工作的有效进行，是图像处理到图像分析的关键步骤。现有的图像分割方法可分 以下几类：基于区域的分割方法，基于边缘的分割方法，基于聚类的方法。在基于聚类的方 法中，是将聚类的理论应道到图像分割上，将图像分割问题，转化为一个求解目标函数的最 优值问题。现有的聚类算法解决图像分割问题主要包括，蚁群算法，神经网络算法，遗传算 法和粒子群算法，其中粒子群优化算法的目标是解决现实世界中的各种最优化问题，该特 点正好和图像分割思想符合，因此粒子群优化算法广泛被学者应用到图像分割中去。但是 粒子群算法虽然实现简单，参数少，但是容易陷入局部最优，参数的选取的合适与否直接影 响到图像分割的速度和精度。 章慧等人在"ImageThresholdSegmentationMethodBasedonImproved ParticleSwarmOptimization"（ComputerScience39 (2012) 289-301)中提出一种基于 混沌粒子群算法的图像阈值分割方法，该方法引入混沌思想，对全局最优值进行混动扰动， 因此增强了算法的全局搜索能力，能够在一定程度上避免陷入局部最优的状态，但是算法 的分割精度和时间复杂度度提高了。
技术实现思路
本专利技术所要解决的技术问题的提供一种基于粒子群算法的自适应脊柱CT图像分 割方法，该方法可改善算法的自适应性，提高求解速度和精度，加速算法收敛。 本专利技术解决上述问题的技术方案如下： -种基于粒子群算法的自适应脊柱CT图像分割方法，该方法包括如下步骤： (1)输入原始图像，计算出图像的最小灰度值和最大灰度值； (2)初始化第一代种群：在图像的最小灰度值至最大灰度值之间分别随机产生一 个整数来初始化第一代种群任一个体的每一维，再以同样的方法初始化第一代种群下一个 体的每一维，依次完成第一代种群每一个体的初始化； (3)按以下方法计算个体最优值和全局最优值： 3a)利用0STU法得到图像的类间距方差，将类间距方差作为适应度函数； 3b)将种群个体带入到适应度函数计算出种群个体的适应度函数值，然后，先挑选 种群个体中适应度最大的作为个体得到个体最优，再从个体最优中挑选适应度最大的个体 作为全局最优； (4)按以下方法产生新的个体： 4a)先按下述公式公式（I)计算出k是自适应权重：k=w· (l-(t/n)1/r) (I)， 式（I)中，t表示当前的迭代次数，n表示最大迭代次数，r是大于等于1的常 数，《为惯性因子； 4b)再将自适应权重按下述公式（II)产生新一代的个体： xid(t+l) =k·xid(t)+C! ·randiO· (pid-xid (t))+c2 ·rand2 () · (pgd-Xid(t)) (II), 式（II)中，(^和c2为基本粒子群算法中的加速常数，randi()和rand2〇表示两 个在范围里变化的随机值，Pld是个体最优值，pgd是全局最优值，Xld(t)为第t代的 第i个个体的第d维，xld(t+l)为第t+Ι代的第i个个体的第d维； (5)按以下方法产生新的个体最优和全局最优： 将新个体带入适应度函数，计算新的种群新个体的适应度函数值；当新个体的适 应度值比个体最优大，那么则用新个体代替个体最优，否则个体最优不变；当新个体的适应 度值比全局最优的适应度函数值大则用个体最优代替全局最优，否则全局最优不变； (6)判断是否满足最大迭代次数，如果满足，得到全局最优值，否则返回步骤（4); (7)以全局最优的每一维数据作为阈值，对得到的图像矩阵进行分割；(8)采取拓扑学操作，对脊柱进一步准确的分割，输出图像。 上述方案中，步骤3a)所述的0STU法是：设图像被{t」，t_2，…，t_MJ个阈值分割 为Μ个部分，则最佳分割阈值{V」，t%…，t%J需要满足如下条件： {t*-!，t*-2,…，t*-MJ=argmax{ 〇 2 (t!，t2,…，tM !)} (III) 其中，σ2(怀t2,···，tMD为图像灰度值的方差值，"Ck为分割图像，k= 1，2,…，M，PRh⑴/N为各阈值的 计算概率，iV= 为图像的总像素数，其中，L为所述图像的灰度值，h(i)为灰度值为 i的像素数； 分割图像(^对应阈值0,…，t」，分割图像(：2对应阈值t」+l，…，t_2，...，分割图像 CM对应阈值t-M !，…，L，0 <t-t-…<t-M ^L。 本专利技术与现有的技术相比具有以下优点： 本专利技术基于以下思路：假设空间上有一个起始点和一个目标点，在算法的刚开始 阶段，粒子距离目标较远，可以加大搜索幅度，大步接近目标点，随着迭代次数的增加，粒子 距离目标越来越近时，粒子搜索幅度变小，搜索的幅度更加精细，局部搜索能力变强，小幅 度的逼近目标，最终到达目标点。也就是说，在给定相同的迭代次数的情况下，本专利技术和其 他粒子群算法的区别在于，对于每一次迭代，本专利技术给分配给每次迭代的搜索幅度是不同 的，本专利技术搜索幅度随着迭代次数的增加自适应的减少，而其他粒子群算法，每次迭代的 搜索幅度都是不变的。 粒子群算法的速度的更新幅度会直接影响到算法的搜索性能，而公式（I)中的 自适应权重k关系到搜索的幅度大小，因此本专利技术通过改进权重k，通过控制权重的大小来 控制搜索幅度，进而将当前迭代次数和总迭代次数考虑进去，自适应改变搜索幅度，在算法 刚开始时，搜索幅度较大，搜索面积较广，可以较快的接近全局最优，经过不断的迭代，算法 搜索面积减小，搜索精度提尚，最终更加逼近最优解，从而具有更快的收敛速率，大大提尚 了图像分割精度。 另外，脊柱的分割容易受到肾脏，肋骨和一些病灶点的干扰，本专利技术利用一系列拓 扑学的操作，通过剔除非连通区域，以及通过检测脊柱的宽度变化去除肾脏，肋骨的影响， 以达到对脊柱的精确分割。【附图说明】 图1是本专利技术的流程图。 图2是公式（I)的函数图像。 图3是脊柱的原始CT图像。 图4是初步分割图，其中，（a)图是本专利技术所述方法的初步分割图，（b)图为传统方 法的初步分割图。 图5是本专利技术所述方法的最终分割的脊柱区域图像。 图6是本专利技术所述方法的最终分割的脊柱区域映射到原始图像的边界图。【具体实施方式】 下面结合附图对本专利技术做进一步的描述。 参照图1，本专利技术的具体实现步骤如下： 步骤1，输入图3所示的原始图像，求出图像的最小灰度值和最大灰度值。 步骤2,初始化第一代种群；最小灰度值至最大灰度值之间分别随机产生一个整 数来初始化种群的个体的每一维；本专利技术实施例中将维数取为1。 步骤3,计算个体最优值和全局最优。 设图像被{t」，t_2,…，t_MJ个阈值分割为Μ个部分，则最佳分割阈值{V」，t%… ，t%J需要满足如下条件： {t*-!，t*-2,…，t*-M=argmax{σ2 (t!，t2,…，本文档来自技高网...
【技术保护点】
一种基于粒子群算法的自适应脊柱CT图像分割方法，该方法包括如下步骤：(1)输入原始图像，计算出图像的最小灰度值和最大灰度值；(2)初始化第一代种群：在图像的最小灰度值至最大灰度值之间分别随机产生一个整数来初始化第一代种群任一个体的每一维，再以同样的方法初始化第一代种群下一个体的每一维，依次完成第一代种群每一个体的初始化；(3)按以下方法计算个体最优值和全局最优值：3a)利用OSTU法得到图像的类间距方差，将类间距方差作为适应度函数；3b)将种群个体带入到适应度函数计算出种群个体的适应度函数值，然后，先挑选种群个体中适应度最大的作为个体得到个体最优，再从个体最优中挑选适应度最大的个体作为全局最优；(4)按以下方法产生新的个体：4a)先按下述公式公式(Ⅰ)计算出k是自适应权重：k＝w·(1‑(t/n)1/r) (Ⅰ)，式(Ⅰ)中，t表示当前的迭代次数，n表示最大迭代次数，r是大于等于1的常数，w为惯性因子；4b)再将自适应权重按下述公式(Ⅱ)产生新一代的个体：xid(t+1)＝k·xid(t)+c1·rand1()·(pid‑xid(t))+c2·rand2()·(pgd‑xid(t)) (Ⅱ)，式(Ⅱ)中，c1和c2为基本粒子群算法中的加速常数，rand1()和rand2()表示两个在[0，1]范围里变化的随机值，pid是个体最优值，pgd是全局最优值，xid(t)为第t代的第i个个体的第d维，xid(t+1)为第t+1代的第i个个体的第d维；(5)按以下方法产生新的个体最优和全局最优：将新个体带入适应度函数，计算新的种群新个体的适应度函数值；当新个体的适应度值比个体最优大，那么则用新个体代替个体最优，否则个体最优不变；当新个体的适应度值比全局最优的适应度函数值大则用个体最优代替全局最优，否则全局最优不变；(6)判断是否满足最大迭代次数，如果满足，得到全局最优值，否则返回步骤(4)；(7)以全局最优的每一维数据作为阈值，对得到的图像矩阵进行分割；(8)采取拓扑学操作，对脊柱进一步准确的分割，输出图像。...
【技术特征摘要】
【专利技术属性】
技术研发人员：李彭军，严静东，杨小燕，朱旭阳，
申请(专利权)人：南方医科大学，广州南方宜信信息科技有限公司，
类型：发明
国别省市：广东;44
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。