基于字典学习的医学图像交互式联合分割制造技术

本发明专利技术公开了一种基于字典学习的医学图像交互式联合分割方法，主要解决现有交互式分割腹部图像计算量大的问题。其实现过程是：输入腹部CT图像序列，用分水岭方法进行初分割并提取每个区域的灰度梯度特征；然后从中选出一幅图像由用户交互式标记目标区域和背景区域；再基于区域间的最大相似性准则对未标记的区域进行合并，直到整幅图像被分为目标和背景两部分为止；然后对得到的目标和背景区域进行训练生成目标字典和背景字典，再用这些字典去逼近序列中待分割的图像并找到逼近误差最小的区域作为标记目标和标记背景，最后用区域合并方法分割出目标。本发明专利技术具有对医学图像分割效果好、省时省力的优点，可用于对腹部CT序列图像的分割。

【技术实现步骤摘要】

本专利技术属于图像处理
，特别是医学图像的处理，主要可用于腹部CT图像的胃癌病变区域分割。
技术介绍
医学图像分割是现代医学图像处理中的一个重要的研究领域，也是病变区域提取、检测、跟踪等后续医学操作的基础。因为这些都需要有一个准确的分割结果为前提，然后才能对各组织进行准确地定位和计算，所以医学图像分割结果的准确性对于医生对病情的诊断并做出正确治疗方案具有非常重要的意义。医学图像分割的方法大致能分为自动分割，手动分割和交互式分割。手动分割费时费力，其结果完全取决于操作者的经验，并且分割结果无法重复。由于医学图像有复杂的背景，自动分割运算量大，在没有人工指导的情况下准确率很难保证。交互式分割属于半自 动分割，它允许人工参与提供一些先验信息，但又不完全依赖于人工参与。由于交互式分割在允许人工参与的情况下保证了准确率。因此，交互式分割方法是现在应用最热门和最广泛的处理医学图像的方法。它的具体做法是由用户标记出图像中感兴趣区域的大致范围，从而提供部分先验信息，然后分割算法再根据这些信息对图像进行分割。交互式分割只需要人工提供少量信息，就能适应不同的图像特性和需求，提高分割准确率，同时可以有效降低计算复杂度。虽然交互式分割有以上优点，但它也存在缺点，那就是一次只能分割一幅图像，用户必须对每幅待分割的图像标记出目标和背景，然后一幅一幅的分割。对于CT图像，一个病人的序列图有上百幅，如果单幅分割，则极大的增加了工作量。
技术实现思路
本专利技术的目的在于针对上述交互式图像分割必须单幅处理的耗时耗力缺陷，提出一种基于字典学习的医学图像交互式联合分割方法，实现对一个人的CT序列中的多幅图像的分割，避免单幅交互式分割的繁琐性，节省时间与精力，减小工作量。实现本专利技术上述目的技术方案，包括如下步骤(I)输入腹部CT序列图像I = U1, 12，…，IJ，n为从序列图中选取图像的数量，η最大值为100，并用分水岭方法对腹部CT序列图像I进行初始分割，得到初分割子图像Γ ;(2)统计每幅图初分割后的区域的数目L，并计算每一个区域的灰度梯度特征向量;(3)从序列图像I中随机选取一幅图像Ii, I彡i彡η，由用户在其对应初分割子图I, i上标记部分目标区域Mtj和部分背景区域Mb，其中Mtj包含初分割子图I' ^勺匕个区域，Mb包含初分割子图I'个区域，剩余的未标记区域N包含初分割子图I',的s个区域，s的大小为L-kfh，该目标区域代表腹部CT图像中的胃癌病变区域；(4)基于最大相似性准则对上述标记后的图像进行区域合并，即先将未标记区域N中与标记背景Mb中相似性最大的区域合并得到新的标记背景区域Μ, B，重复合并剩下的未标记区域与M' B，直到s = O合并完成后，得到了图像Ii的全部目标区域和背景区域；(5)用K-SVD算法对步骤⑷中得到的目标区域所对应的特征向量训练生成一个目标字典％，并将步骤(4)中得到的背景区域用k均值方法聚为4类，然后分别对每类训练生成背景字典Db = {DB1，Db2, Db3, Db4}，其中Db1，Db2, Db3, Db4分别为4类背景区域对应的背景字典；(6)由目标字典Dq去逼近一幅待分割的图像Γ」的特征向量X，1彡j彡n，j关i，计算出逼近误差最小的区域作为标记目标区域由4个背景字典DB1，Db2, Db3, Db4分别去逼近待分割图像I, j的特征向量X，计算出的4个逼近误差最小的区域一起作为标记背景区域Mb，再次按照步骤(4)进行区域合并，得到图像L的目标区域；(7)重复执行步骤(6)，对输入的一组序列图像I中还未分割的图像进行目标区域提取，直到所有序列图像都分割出目标区域为止。本专利技术与现有的技术相比具有以下优点I、本专利技术由于只需要对一组CT序列图像中的一幅进行简单的划线操作来人机交互，其余的图像不再需要人工参与便可分割出目标，因而不仅使多幅效果和单幅交互式分割一样，而且极大的节省了时间并减少了工作量；2、本专利技术针对CT图像自身的特点，用灰度梯度作为区域特征，更好的表示了区域特性，并提高了分割效果。附图说明图I是本专利技术的流程图；图2是本专利技术所使用的一个人的3幅腹部CT图像序列，其中2a为序列图的第78幅图，2b为第89幅图，2c为第96幅图；图3是本专利技术对序列图进行初始分割后的结果图像，其中图3a、3b和3c分别是对图2a、2b和2c初始分割后的图像；图4是本专利技术对图3a手动划线标记后的图像；图5是本专利技术对图4进行区域合并后的结果图像；图6是本专利技术通过字典匹配后区域合并的分割结果图。具体实施例方式参照图1，本专利技术基于字典学习的医学图像交互式联合分割方法给出以下实施例步骤I :输入一组腹部CT图像序列，如图2a，2b和2c。步骤2 :用分水岭算法对原图2a，2b和2c分别进行初始分割得到图3a、3b和3c，所述的分水岭算法是参考自文献“王忠义，“分水岭变换在医学图像分割中的应用”，电脑知识与技术，V01. 5，No. 12，April 2009，pp. 3221-3223”，具体分割步骤如下2a)对图2a中每幅图的每个像素的灰度级进行由低到高排序；2b)再对排序后的灰度级从低到高实现淹没，淹没过程中对每一个局部极小值在高度为h灰度级的影响域采用先进先出FIFO进行判断及标注，即可得到分割为多个区域的图3a，其中h的范围为0-255，根据CT图像特性，本实例选取16 ；2c)对于图2b和2c，重复执行步骤2a)和2b)，得到初分割后的图3b和图3c。步骤3 :分别统计图3a、图3b和图3c中所有区域的数目L及每个区域中的像素点数，提取每一个区域中每个像素点的灰度梯度特征，然后将每个区域内所有像素点的特征向量均值作为这个区域的特征向量。步骤4 :选取图2a做交互式分割，针对要分割的目标在图3a上进行标记。4a)根据医生经验，胃部淋巴结主要出现在胃壁周围5cm内，把胃壁周围5cm以内的病变区域作为要分割的目标； 4b)由用户在图3a上划不同的线段分别对目标和背景进行标记，如图4所示，例如用绿色线段标记目标区域Mtj,用蓝色线段标记部分背景区域Mb，剩余的未标记区域为N，其中Mq包含图3a的kQ个区域，Mb包含图3a的kB个区域，剩余的未标记区域N包括了图3a的s个区域，s的大小为L- ，该目标区域代表腹部CT图像中的胃癌病变区域。步骤5 :基于最大相似性准则对图4进行区域合并，即先将未标记区域N中与标记背景Mb中相似性最大的区域合并后得到新的标记背景区域M' b中，重复合并剩下的未标记区域与M' B，直到s = O合并完成后，如图5得到了图像Ii的全部目标区域和背景区域，所述的区域合并方法是参考自文献“Jifeng Ning等人，“Interactive image segmentationby maximal similarity based region merging”， Pattern Recognition,43,2010,pp. 445-456”。步骤6 :用K-SVD算法对图5的目标区域的特征向量训练生成一个目标字典％，用k均值将图5的背景区域聚为4类，然后分别对每类中各区域的特征向量训练生成背景字典Db — {DB1，Db2，Db3，Db4I， 其中 Dbi，Db2 Db本文档来自技高网...

【技术保护点】

【技术特征摘要】

【专利技术属性】
技术研发人员：缑水平，马丽敏，周治国，刘芳，张晓鹏，唐磊，王之龙，王云利，
申请(专利权)人：西安电子科技大学，
类型：发明
国别省市：