一种基于聚类融合算法的双模态脑肿瘤MRI分割方法技术

本发明专利技术涉及医学图像分析技术领域，且公开了一种基于聚类融合算法的双模态脑肿瘤MRI分割方法，包括以下步骤：步骤1输入肿瘤患者的不同模态图像，包括T2加权和Flair图像，对输入的图像进行灰度值归一化；步骤2筛选图像中脑部区域对应的像素点。该基于聚类融合算法的双模态脑肿瘤MRI分割方法，通过利用了不同MRI模态的成像特点，结合不同模态的信息和邻域信息，提高分割的准确度，在K均值算法中结合使用三种距离计算方法，改善了使用单个距离分割效果不稳定的情况，对于来自不同设备的图像具有鲁棒性，能够自动识别病变区域，本发明专利技术提出的后处理方法，有效的提升了分割的准确率，从而达到了算法实施简单、高效和运算速度块的目的。

A Bimodal Brain Tumor MRI Segmentation Method Based on Cluster Fusion Algorithms

【技术实现步骤摘要】
一种基于聚类融合算法的双模态脑肿瘤MRI分割方法
本专利技术涉及医学图像分析
，具体为一种基于聚类融合算法的双模态脑肿瘤MRI分割方法。
技术介绍
脑肿瘤列肿瘤相关性死亡的首位病因，其中胶质瘤是成人中最常见的原发性脑肿瘤，占成人恶性原发性肿瘤的70％，平均存活时间为1年，脑肿瘤具有生长速度快、复发率高和致残率高等特点，因此早期诊断和治疗能够减少对患者造成的痛苦，在所有影像学诊断方法中，核磁共振成像被认为是诊断脑肿瘤最敏感的方法，具有对比度高和无创的优点，MRI脑肿瘤分割为神经病理学分析提供了有价值的信息，可以辅助于规划治疗策略、检测疾病进展和预测患者预后。MRI脑肿瘤分割具有一定的技术难度，主要体现在以下方面：(1)胶质瘤可能会出现在大脑的任何位置，其形状、结构和大小各不相同，(2)部分高级别的胶质瘤，显弥漫、浸润性生长，在图像中表现为肿瘤部位和周围脑组织的灰度值梯度小，边界模糊，(3)脑组织的结构复杂，正常的脑组织包括灰质、白质、脑室和脑脊液等，病变组织包括肿瘤和水肿区域，其中肿瘤又包括坏死区、增强区和非增强区，(4)由于MRI采集设备和参数不同、外界噪声和患者个体差异等原因，导致个体图像差异较大，通过调整MRI图像的测量参数和注射造影剂，可以获得多种模态的MRI，常见的模态有，T1加权、T2加权、T1增强和Flair，不同的模态强调信息各有不同，结合多模态的信息有助于解决以上技术难点，图像分割方法分为有监督分割和无监督分割，目前用于MRI的有监督分割方法主要有基于卷机神经网络，和基于地图集的方法，有监督分割需要大量的人工标注数据来训练模型，但是人工标注数据和训练模型的过程耗费很多人力、资源和时间，而无监督方法具有省时省力，适用性广等优点，基于聚类的无监督方法以像素点为样本，根据其特征分为若干的类，同一个类别中的像素点彼此相似，距离小，与其他类别中的像素点相异，距离大。目前常用的聚类方法有K均值聚类、模糊C均值聚类、谱聚类、混合高斯模型和层次聚类等，有学者提出了水平集与K均值聚类结合的脑肿瘤分割方法，这种方法初始化水平集轮廓对后续结果有影响，需要调整能量函数的参数，但是脑肿瘤图像个体间差别很大，很难找到对所有图像都合适的参数，对于脑肿瘤MRI病变区分割的问题，K均值聚类相比于其他聚类方法，速度更快，算法更简单高效，然而K均值聚类不同的距离函数对图像分割的影响较大，单一的距离函数不能够满足所有类型图像的分割，现有的基于聚类的图像分割方法中，将像素点聚为若干类簇，但是如何从若干类簇中自动选择目标区域仍待解决，而且聚类结果中通常包含大量噪声点和错误分类的像素点，本专利技术提出了多种距离函数融合的K均值聚类方法，可以从分割后的图像中自动识别病变区域，以及有效去除噪声点的后处理方法。
技术实现思路
(一)解决的技术问题针对现有技术的不足，本专利技术提供了一种基于聚类融合算法的双模态脑肿瘤MRI分割方法，具备能够实现简单快速的无监督脑肿瘤图像分割等优点，解决了现有基于聚类的脑肿瘤MRI分割方法，算法复杂，速度较慢，参数较多，K均值聚类简单快速，但是效果不稳定，对不同MRI切片的分割效果参差不齐的问题。(二)技术方案为实现上述能够实现简单快速的无监督脑肿瘤图像分割的目的，本专利技术提供如下技术方案：一种基于聚类融合算法的双模态脑肿瘤MRI分割方法，包括以下步骤：步骤1输入肿瘤患者的不同模态图像，包括T2加权和Flair图像，对输入的图像进行灰度值归一化；步骤2筛选图像中脑部区域对应的像素点，提取像素点的多模态的邻域像素点灰度值作为该像素点的特征；步骤3结合三种距离的K均值聚类对像素点分类，自动调整聚类数目；步骤4从若干类中自动选出包含肿瘤的类；步骤5使用一系列形态学处理，除去错误分类的像素点，精确提取病变区，包括肿瘤和水肿。优选的，将图像灰度范围统一线性归一化为0至255的范围，去除背景区域的像素点，只计算包含脑区的部分(非零区)，如此便组成了n×p的样本空间，其中n为待聚类的像素点数，p为特征维数。优选的，对像素点聚类的三种距离具体为：平方欧氏距离、余弦距离和城市街区距离。优选的，基于聚类融合算法的算法进一步细化如下：Stpe1：输入样本空间s，样本为N个像素点，根据T2加权和Flair图像特点，病变区域在两图中均显高亮度，白质在两图中均显较低亮度，灰质和脑室在T2加权图中显高亮度而在Flair图中显低亮度，故设置k均值算法初始聚类数K＝3，从N个样本中随机选取K个类簇的聚类中心计算每个样本到聚类中心的平方欧氏距离:Stpe2：将每个样本划分到距离值最小的聚类中心的类簇中，计算每个类簇中样本的平均值，作为新的聚类中心：Stpe3：其中表示第t次迭代聚类中心的第k个类簇，表示第k个类簇中样本的个数，这里的求和是指类簇中所有元素在每列属性上的和，是一个p维向量，表示为：ck＝(ck,1,ck,h,…,,ck,h)(3)Stpe4：如此循环迭代，直至t＝T时停止迭代，得到k个最终聚类中心计算n个样本与k个聚类中心的距离，得到大小为n×k的距离矩阵Dsqe，元素表示i个样本到j个聚类中心的距离，其中i＝1,…,nj＝1,…,k，第一次聚类中k＝3，将样本分为3类，将(1)式平方欧氏距离换为余弦距离：Stpe5：重复上述算法T次迭代得到距离矩阵Dcos，再将平方欧氏距离换为城市街区距离：j为x和的第j维，重复上述算法T次迭代得到距离矩阵Dcb。Stpe6：将Dsqe，Dcos，Dcb中相应类别编号对应，将用三种距离聚类的样本定位至原图中的像素点，由此分别获得三种矩阵对应的三个聚类结果图，首先，k类像素点中在Flair图像上平均灰度值最大的一类，记为类别1，距离矩阵的第1行表示类别1的聚类中心，接下来，计算剩余聚类中心与第一类聚类中心的距离，按照距离由小到大的顺序，分别记为类别2到类别k，对应于距离矩阵的第二至第k列，最后，得到了列换序后的矩阵Dsqe*，Dcos*，Dcb*，对三个距离矩阵各个元素求和得到Dadd：Dadd＝Dsqe*+Dcos*+Dcb*(6)将样本x归为与其距离最小的聚类中心的类别。(三)有益效果与现有技术相比，本专利技术提供了一种基于聚类融合算法的双模态脑肿瘤MRI分割方法，具备以下有益效果：该基于聚类融合算法的双模态脑肿瘤MRI分割方法，通过利用了不同MRI模态的成像特点，结合不同模态的信息和邻域信息，提高分割的准确度，在K均值算法中结合使用三种距离计算方法，改善了使用单个距离分割效果不稳定的情况，对于来自不同设备的图像具有鲁棒性，能够自动识别病变区域，本专利技术提出的后处理方法，有效的提升了分割的准确率，从而达到了算法实施简单、高效和运算速度块的目的。附图说明图1是本专利技术的多模态MRI脑肿瘤图像无监督分割方法的总体流程图。图2是本专利技术改进的K均值聚类算法的流程图；图3是本专利技术实施例的示例性流程示意图；图4是本专利技术的方法中形态学操作的实施例的示例性流程示意图；图5是本专利技术在BRATS2015真实临床数据上的脑肿瘤分割结果示意图；图6是本专利技术算法在10个脑肿瘤图像的相似系数。具体实施方式本专利技术以脑肿瘤分割的应用为目标，从分割结果中提取肿瘤特征，用于制定治疗计划、分析肿瘤生长趋势和治疗效果评估，同一组织在不同本文档来自技高网...

【技术保护点】
1.一种基于聚类融合算法的双模态脑肿瘤MRI分割方法，其特征在于，包括以下步骤：步骤1输入肿瘤患者的不同模态图像，包括T2加权和Flair图像，对输入的图像进行灰度值归一化；步骤2筛选图像中脑部区域对应的像素点，提取像素点的多模态的邻域像素点灰度值作为该像素点的特征；步骤3结合三种距离的K均值聚类对像素点分类，自动调整聚类数目；步骤4从若干类中自动选出包含肿瘤的类；步骤5使用一系列形态学处理，除去错误分类的像素点，精确提取病变区，包括肿瘤和水肿。

【技术特征摘要】
1.一种基于聚类融合算法的双模态脑肿瘤MRI分割方法，其特征在于，包括以下步骤：步骤1输入肿瘤患者的不同模态图像，包括T2加权和Flair图像，对输入的图像进行灰度值归一化；步骤2筛选图像中脑部区域对应的像素点，提取像素点的多模态的邻域像素点灰度值作为该像素点的特征；步骤3结合三种距离的K均值聚类对像素点分类，自动调整聚类数目；步骤4从若干类中自动选出包含肿瘤的类；步骤5使用一系列形态学处理，除去错误分类的像素点，精确提取病变区，包括肿瘤和水肿。2.如权利要求1所述的基于聚类融合算法的双模态脑肿瘤MRI分割方法，其特征在于，将图像灰度范围统一线性归一化为0至255的范围，去除背景区域的像素点，只计算包含脑区的部分(非零区)，如此便组成了n×p的样本空间，其中n为待聚类的像素点数，p为特征维数。3.如权利要求1所述的基于聚类融合算法的双模态脑肿瘤MRI分割方法，其特征在于，对像素点聚类的三种距离具体为：平方欧氏距离、余弦距离和城市街区距离。4.如权利要求1所述的基于聚类融合算法的双模态脑肿瘤MRI分割方法，其特征在于，基于聚类融合算法的算法进一步细化如下：Stpe1：输入样本空间s，样本为N个像素点，根据T2加权和Flair图像特点，病变区域在两图中均显高亮度，白质在两图中均显较低亮度，灰质和脑室在T2加权图中显高亮度而在Flair图中显低亮度，故设置k均值算法初始聚类数K＝3，从N个样本中随机选取K个类簇的聚类中心计算每个样本到聚类中心的平方欧氏距离:Stpe2：将...

【专利技术属性】
技术研发人员：施建宇，张安琪，牛晨，李奎，卢慧甍，尚晓娅，李静迪，
申请(专利权)人：西北工业大学，
类型：发明
国别省市：陕西,61