基于多尺度空间模糊聚类的无监督mri图像分割方法
【技术领域】
[0001] 本发明涉及一种图像分割方法，尤其是基于多尺度空间模糊聚类的无监督MRI图 像分割方法。
【背景技术】
[0002] 多尺度分析是正确认识事物和现象的重要方法之一，目前广泛地用于医学图像分 析中。但是目前已有的方法普遍存在计算方法太简单化，在一定条件下容易受到图像对比 度与亮度变化的影响，如果使用人员要求得最佳阈值，往往需要进行较大规模的遍历计算 并计算出方差，当计算量大时效率会很低。同时，在实际图像中，由于图像本身灰度分布以 及噪声干扰等因素的影响，目前以后的方法并不能使图像分割得到满意的结果，虽在一定 程度上可以消除噪声的影响，但该方法计算量相当大，难以应用到实时系统。
[0003] 由于磁共振成像（Magnetic Resonance Imaging，简称MRI)设备获取的人脑核磁 图像受噪声、射频场不均匀性、脑不同组织之间的差异性以及部分容积效应等因素影响，造 成人脑核磁图像的均匀性变差，因此仅依靠图像的灰度密度信息给精确的脑图像分类带来 了很大困难，如果要得到正确的脑组织分类，首先必须对灰度进行校正。
[0004] 近年来，一些基于模糊聚类的多尺度空间的图像数据处理技术被陆续提出，该类 方法的过程通常是在尺度空间图像序列的每一个尺度层次上直接进行模糊聚类方法，上一 尺度层次的分割结果经常作为下一尺度层次的初始状态，由于没有引进有效的不同的尺度 层次之间的约束来优化当前尺度层次内的分割，所以分割结果只在最后分割的尺度层次得 到优化，因此，这些方法对于退化图像的鲁棒性不是很好。
【发明内容】
[0005] 本发明所要解决的技术问题是提供一种图像分割结果较为准确且有效地去除图 像中的噪声、提高了退化图像的鲁棒性的基于多尺度空间模糊聚类的无监督MRI图像分割 方法。
[0006] 本发明解决上述技术问题所采用的技术方案为：基于多尺度空间模糊聚类的无监 督MRI图像分割方法，包括以下步骤：
[0007] ①从医院的MRI医学影像设备中获取一幅待分割的MRI扫描图像，然后将该待分 割的MRI扫描图像转换成灰度图像，再采用模糊聚类算法对灰度图像进行分析，得到初始 感兴趣的轮廓认1，厶2,43，*"^11}，其中41、4233、*"、411表示构成初始感兴趣的轮廓 的所有边缘点；
[0008] ②根据初始感兴趣的轮廓{Al，A2, A3，…，An}，经过包络获得灰度图像中的感兴 趣区域，根据该感兴趣区域的外形选取一个能够将该感兴趣区域包含在内的方形区域，定 义该方形区域的宽度为W、高度为H，如果该方形区域的尺寸大小WXH能够被uXu整除，则 H 7 X /7 将该方形区域定义为当前灰度图像，然后直接将当前灰度图像划分成1--个互不重叠的 U X U 尺寸大小为uXu的子块；如果WXH不能够被uXu整除，则扩展该方形区域直至其尺寸大 小W' XH'恰好能够被uXu整除，其中，W'表示经过扩展的方形区域的宽度，H'表示经过 扩展的方形区域的高度，再将经过扩展的方形区域定义为当前灰度图像，然后将当前灰度 图像划分试
个互不重叠的尺寸大小为uXu的子块，W' = = W且H'>H，或W'>W且H' ==H，或 W' >W 且 H' >H，u 取 4 或 5 或 6 ;
[0009] ③采用标记分水岭分割算法对当前灰度图像进行分割，得到经过分割后的多块初 步目标区域；
[0010] ④对获得的所有初步目标区域中当前待处理的一块初步目标区域定义为当前初 步目标区域；
[0011] ⑤从当前灰度图像中提取出与当前初步目标区域对应的所有尺寸大小为uXu的 子块，对提取出的所有子块中当前待处理的一个子块定义为当前子块；
[0012] ⑥将当前子块中的每个像素点的像素值分别作为输入参数输入非线性优化模型 中进行优化，该非线性优化模型为：
[0014] 其中表示非线性优化模型输出的信噪比，A为固定周期信号的幅度，f。为固定 Llt 周期信号的频率，t为布朗粒子的运动时间，Φ为固定周期信号的初相位，m、η均为双稳态 势皇实参数，g(t)表示非线性优化模型的输入参数，s为布朗粒子的运动坐标，α为噪声强 度，ξ (t)为均值为0的高斯白噪声；
[0015] 将当前子块中的每个像素点的像素值分别作为输入参数输入非线性优化模型中 进行优化后，由非线性优化模型输出当前子块中的每个像素点的信噪比，如果当前子块中 的每个像素点的信噪比与事先存储于数据库中的MRI扫描图像的灰度图像中对应的像素 点的信噪比相比，误差小于10%，则判定当前子块优化成功，将优化成功的当前子块定义为 保留子块，其中，事先存储于数据库中的MRI扫描图像的灰度图像的尺寸大小与当前灰度 图像的尺寸大小相同；
[0016] ⑦判断提取出的所有子块是否全部处理完毕，若全部处理完毕，则得到由所有保 留子块形成的最终目标区域，并执行步骤⑧；若有未处理的子块，则将提取出的所有子块中 下一个待处理的子块定义为当前子块，返回步骤⑥；
[0017] ⑧判断获得的所有初步目标区域是否全部处理完毕，若全部处理完毕，则得到当 前灰度图像中经过分割后形成的多个最终目标区域，并执行步骤⑨；若有未处理的初步目 标区域，则将获得的所有初步目标区域中下一个待处理的初步目标区域定义为当前初步目 标区域，返回步骤⑤；
[0018] ⑨完成对MRI扫描图像的分割过程。
[0019] 在步骤②中，u = 5。
[0020] 在步骤②中，根据初始感兴趣的轮廓{A1，A2,A3，…，An}，经过包络获得感兴趣的 灰度图像区域的具体过程为：首先根据初始感兴趣的轮廓的外形选取一个最大范围的方形 区域，使其能够被uXu整除，再在感兴趣区域内除已选取的方形区域之外的剩余区域中选 取一个最大范围的方形区域，使其能够被uXu整除，重复以上步骤直至所有选取出的方形 区域均能够被uXu整除，最后将无法被uXu整除的区域通过空白像素填充直至恰好能够 被uXu整除。
[0021] 步骤⑥中事先存储于数据库中的MRI扫描图像的灰度图像中的每个像素点的信 噪比的具体获取过程为：
[0022] ⑥-1选取一幅具有代表性的MRI扫描图像，该具有代表性的MRI扫描图像的宽度 为W，高度为H，然后将该MRI扫描图像转换为灰度图像；
[0023] ⑥-2如果WXH能够被uXu整除，则将该灰度图像作为待处理灰度图像，然后直 接将待处理灰度图像划分成
个互不重叠的尺寸大小为uXu的子块；如果WXH不能 够被uXu整除，则扩展该灰度图像直至其尺寸大小恰好能够被uXu整除，将扩展后的灰度 图像作为待处理灰度图像，然后将待处理灰度图像划分成
个互不重叠的尺寸大小 为uXu的子块，其中，W'表示扩展后的灰度图像的宽度，H'表示扩展后的灰度图像的高度， W' = = W 且 Η' >H，或 W' >W 且 H' = = H，或 W' >W 且 Η' >H ;
[0024] ⑥-3手动圈定待处理灰度图像中的每个目标区域；
[0025] ⑥-4计算得到待处理灰度图像中的每个目标区域对应的每个尺寸大小为uXu的 子块中的每个像素点的信噪比。
[0026] 与现有技术相比，本发明的优点在于在实际图像中，首先将待分割的MRI扫描图 像转换成灰度图像，使图像分割的结果不容易受到MRI扫描图像本身灰度分布以及噪声干 扰等因素的影响，提高了图像分割结果的准确度；采用基于多尺度空间模糊聚类算法对灰 度图像进行分析，获得感兴趣的灰度图像区域，并将感兴趣的灰度图像区域划分为互不重 叠的子块，采用区域生长法对当前灰度图像进行分割，综合考虑了所有像素，提高了像素分 割的准确性，避免了现有分割方法的不足；通过将参数输入非线性优化模型进行优化的非 线性扩散技术，可以有效地去除图像中的噪声，最大限度地保护感兴趣的目标；本发明引进 了有效的不同的尺度层次之间的约束，优化了当前尺度层次内的分割，分割结果在当前尺 度层次就得到优化，因此，提高了退化图像的鲁棒性。
【附图说明】
[0027] 图1为本发明的原理框图。
【具体实施方式】
[0028] 以下结合附图实施例对本发明作进一步详细描述。
[0029] 实施例一：基于多尺度空间模糊聚类的无监督MRI图像分割方法，包括以下步骤：
[0030] ①从医院的MRI医学影像设备中获取一幅待分割的MRI扫描图像，然后将该待分 割的MRI扫描图像转换成灰度图像，再采用模糊聚类算法对灰度图像进行分析，得到初始 感兴趣的轮廓认1，厶2,43，*"^11}，其中41、4233、*"、411表示构成初始感兴趣的轮廓 的所有边缘点（参考《基于模糊聚类的多尺度空间数据挖掘模型及其矿山应用》，上海交通 大学学报，第42卷第2期，2008年2月出版）；
[0031] ②根据初始感兴趣的轮廓{Al，A2, A3，···，An}，经过包络获得灰度图像中的感兴 趣区域，根据该感兴趣区域的外形选取一个能够将该感兴趣区域包含在内的方形区域，定 义该方形区域的宽度为W、高度为H，如果该方形区域的尺寸大小WXH能够被uXu整除，则 将该方形区域定义为当前灰度图像，然后直接将当前灰度图像划分姑
个互不重叠的 尺寸大小为uXu的子块；如果WXH不能够被uXu整除，则扩展该方形区域直至其尺寸大 小W' XH'恰好能够被uXu整除，其中，W'表示经过扩展的方形区域的宽度，H'表示经过扩 展的方形区域的高度，再将经过扩展的方形区域定义为当前灰度图像，然后将当前灰度图 像划分试
个互不重叠的尺寸大小为uXu的子块，W' == W且H'>H，或W'>W且H' ==H，或 W' >W 且 H' >H，u 取 4 或 5 或 6。
[0032] 根据初始感兴趣的轮廓{Al，A2, A3，…，An}，经过包络获得感兴趣的灰度图像区 域的具体过程为：首先根据初始感兴趣的轮廓{A1，A2，A3, 一,An}的外形选取一个最大范 围