使用空间截面投影的脑肿瘤MRI图像分割方法技术

本发明专利技术公开了一种使用空间截面投影的脑肿瘤MRI图像分割方法，采用基于截面投影的三维Otsu法对当前灰度图像进行初步分割，综合考虑了所有像素，提高了像素分割的准确性，避免了现有分割方法的不足；采用本发明专利技术分割方法对脑肿瘤MRI图像进行分割，可获得最好的分割效果，有利于提高对脑肿瘤病例判断的准确性；本发明专利技术分割方法对含不同噪声类型的脑肿瘤MRI图像也能取得较好的分割效果，适用性好；本发明专利技术分割方法的分割精度高、抗噪性强、分割效果好、分割效率高，可应用于噪声干扰严重及时间方面要求较高的系统，表现出重要的研究价值和广阔的应用前景。

【技术实现步骤摘要】

本专利技术属于图像处理领域，具体涉及一种使用空间截面投影的脑肿瘤 MRI (Magnetic Resonance Imaging，磁共振成像）图像分割算法。
技术介绍
图像分割是计算机视觉领域内的一项关键技术，在图像处理与分析间起着重要的 桥梁作用。在众多的分割方法中，阈值分割法因其简单、有效而被广泛应用于图像处理领 域。常用的阈值法主要有最大类间方差法（简称Otsu法）、最大熵法、最小交叉熵法等。 其中日本学者大津展之提出的Otsu法因有模式识别的相关理论为基础且图像分割性能良 好，而逐渐成为一种最为流行的阈值分割技术。由于一维Otsu法仅考虑了像素的灰度信息 而没考虑像素间的空间相关性，故当图像包含噪声时，一维Otsu法的分割效果不佳。对此， 刘健庄等提出了二维Otsu法，该方法同时考虑了像素的灰度信息和其邻域的空间相关性， 改善了一维Otsu法的抗噪性。景晓军等在二维Otsu法的基础上进一步引入了邻域中值灰 度，并由此提出了三维Otsu法，该算法能更好地对低信噪比图像进行分割。 传统三维Otsu法对图像进行分割的过程： 对于尺寸为MXN，灰度级为L的图像，设像素（x，y)的灰度值为f(x，y)，则该像素 KXK邻域的均值灰度g(x，y)和中值灰度h(x，y)可分别定义如下： h (x, y) = med {f (x+m, y+n), m = _k/2,…，k/2 ;n = _k/2,…，k/2} (2) 由g(x，y)及h(x，y)的定义可知，像素的邻域均值灰度和邻域中值灰度的灰度级 也为L。定义f (X，y)、g(x, y)、h(x, y)形成的三元组（i, j, k)为一个三维直方图，则该直方 图定义在一大小为LXLXL的正方体区域内，其3个坐标分别表示像素的灰度值、邻域均值 灰度值和邻域中值灰度值。直方图中任一点的值记为P l jk，Pl jk定义为 其中，Cljk为（i，j，k)出现的频数，0彡i，j，k彡L-1，且 传统三维OtSU法将直方图划分成如图1所示的8个区域。因目标和背景内部的 像素间存在很强的相关性，故像素的灰度值、邻域均值灰度、邻域中值灰度三者非常接近； 而位于目标与背景边界附近的像素，上述3个数值存在明显差异。基于上述认识，传统三维 Otsu法分别将区域0、区域1视为目标和背景，而将区域2~7视为边缘和噪声。由于边界 附近的像素数通常很少，故传统Otsu法假设区域2~7上的所有p 1]k~0。传统三维Otsu 法采用如下的距离测度函数确定最佳阈值 CN 105139410 A ~P 2/6 页 其中，w〇(s, t, q)、 U Ti、U Tj、U Tk、 μ i (s, t, q)、μ .j (s, t, q)、μ k (s, t, q)的含义及计 算方式见文献。 使得式（4)取得最大值的阈值（s。，t。，q。）即为最佳阈值，BP 传统三维Otsu法虽然增强了二维Otsu法的抗噪性，但是仍然存在以下几点不足： 1)与二维Otsu法类似，仅考虑了主对角线附近〇、1两个区域，而假设区域2~7中的概率 为〇,这与实际不符，往往造成分割不准确；2)阈值维数的增加虽能改善分割效果，但却严 重影响了分割效率；3)由式（4)可知，三维Otsu法的计算公式较为复杂，难以扩展到多阈 值的情况；4)三维Otsu法的抗噪性仍不足，有待进一步增强。尽管传统三维Otsu法忽略了 区域2~7中的噪声，但是区域0、1中通常也包含噪声，而传统三维Otsu法并未对其做进 一步处理，从而影响了分割效果。可见，传统三维Otsu法虽然同时考虑了邻域均值和灰度 中值灰度，但是分割效率低、抗噪性不足，故有必要对其进行改进以提高分割效率、增强抗 噪性；此外，由于未对目标和背景区域的噪声进行任何处理，当噪声干扰严重时，三维Otsu 的分割效果仍不够理想。 范九伦等提出的快速递推法虽提高了分割效率，但要以消耗大量的额外空间为代 价。Puthipon提出的等效三维Otsu法虽提高了分割效率且无需额外的空间消耗，但因按照 传统的三维直方图划分方式，仅考虑对角线附近的2个长方体区域而忽略了其他区域，往 往造成分割不准确。基于截面投影的三维Otsu法是现有的对传统三维Otsu法的改进方法， 以截面投影直方图为基础建立，虽然该方法考虑了整个直方图区域，避免了传统三维Otsu 法近似假设造成的不合理性，但是该方法的分割精度和分割效果仍有待提高。
技术实现思路
本专利技术所要解决的技术问题是：针对现有技术的不足，提供一种对图像分割精度 高、抗噪性强、分割效果好、分割效率高的使用空间截面投影的脑肿瘤MRI图像分割算法。 本专利技术解决上述技术问题所采用的技术方案为：使用空间截面投影的脑肿瘤MRI 图像分割算法，包括以下步骤： (1)获取一幅来源于脑部的MRI扫描设备的待分割的脑肿瘤MRI扫描图像，然后将 该待分割的脑肿瘤MRI扫描图像转换成灰度图像，再采用基于截面投影的三维Otsu法对灰 度图像进行分析，得到初始感兴趣的轮廓{Al，A2, A3,…，An}，其中Al, A2, A3,…，An定义为 构成初始感兴趣的轮廓的所有轮廓点； (2)根据初始感兴趣的轮廓{A1，A2, A3,…，An}经包络得到灰度图像中的感兴趣 区域，根据该感兴趣区域的外形选取一个能够将该感兴趣区域包含在内的方形区域，定义 该方形区域的宽度为W、高度为H，如果WXH能够被uXu整除，则将该方形区域定义为当前 灰度图像，然后直接将当前灰度图像划分成^个互不重叠的尺寸大小为uXu的子块； UXU 如果WXH不能够被uXu整除，则扩展该方形区域使其尺寸大小能够被uXu整除，将扩展 W'x IP 后的方形区域定义为当前灰度图像，然后将当前灰度图像划分成个互不重叠的尺寸 u.龍： 大小为uXu的子块，其中，W'和H'对应表示扩展后的方形区域的宽度和高度，W' == W且 Η' >H 或 W' >W 且 H' = = H 或 W' >W 且 H' >H，u 取 4、5 或 6 ; (3)采用区域生长法对当前灰度图像进行分割，得到当前灰度图像中的多个初步 病灶区域； (4)将当前灰度图像中当前待处理的初步病灶区域定义为当前初步病灶区域； (5)从当前灰度图像中提取出与当前初步病灶区域对应的所有尺寸大小为uXu 的子块，按序对当前初步病灶区域对应的每个尺寸大小为uXu的子块进行处理，将当前待 处理的尺寸大小为uXu的子块定义为当前子块； (6)将当前子块中的所有像素点各自的像素值作为输入参数输入非线性优化模型 中进行优化，该非线性优化模型为： 其中：^表示非线性优化模型输出的信噪比，A为固定周期信号的幅度，f。为固定 at 周期信号的频率，t为布朗粒子的运动时间，Φ为固定周期信号的初相位，m、η均为双稳态 势皇实参数，g(t)表示非线性优化模型的输入参数，s为布朗粒子的运动坐标，α为噪声强 度，ξ (t)为均值为0的高斯白噪声； 将当前子块中的所有像素点各自的像素值作为输入参数输入非线性优化模型中 进行优化后，非线性优化模型输出当前子块中的所有像素点各自的信噪比，如果当前子块 中的每个像素点的信噪比与事先存储于数据库中的脑肿瘤MRI本文档来自技高网...

【技术保护点】
使用空间截面投影的脑肿瘤MRI图像分割算法，其特征在于，包括以下步骤：(1)获取一幅来源于脑部的MRI扫描设备的待分割的脑肿瘤MRI扫描图像，然后将该待分割的脑肿瘤MRI扫描图像转换成灰度图像，再采用基于截面投影的三维Otsu法对灰度图像进行分析，得到初始感兴趣的轮廓{A1,A2,A3,…,An}，其中A1,A2,A3,…,An定义为构成初始感兴趣的轮廓的所有轮廓点；(2)根据初始感兴趣的轮廓{A1,A2,A3,…,An}经包络得到灰度图像中的感兴趣区域，根据该感兴趣区域的外形选取一个能够将该感兴趣区域包含在内的方形区域，定义该方形区域的宽度为W、高度为H，如果W×H能够被u×u整除，则将该方形区域定义为当前灰度图像，然后直接将当前灰度图像划分成个互不重叠的尺寸大小为u×u的子块；如果W×H不能够被u×u整除，则扩展该方形区域使其尺寸大小能够被u×u整除，将扩展后的方形区域定义为当前灰度图像，然后将当前灰度图像划分成个互不重叠的尺寸大小为u×u的子块，其中，W'和H'对应表示扩展后的方形区域的宽度和高度，W'＝＝W且H'>H或W'>W且H'＝＝H或W'>W且H'>;H，u取4、5或6；(3)采用区域生长法对当前灰度图像进行分割，得到当前灰度图像中的多个初步病灶区域；(4)将当前灰度图像中当前待处理的初步病灶区域定义为当前初步病灶区域；(5)从当前灰度图像中提取出与当前初步病灶区域对应的所有尺寸大小为u×u的子块，按序对当前初步病灶区域对应的每个尺寸大小为u×u的子块进行处理，将当前待处理的尺寸大小为u×u的子块定义为当前子块；(6)将当前子块中的所有像素点各自的像素值作为输入参数输入非线性优化模型中进行优化，该非线性优化模型为：dsdt=A×sin(2π×f0×t+ψ)+m×g(t)+n×s-m×s3+2α×ξ(t)]]>其中：表示非线性优化模型输出的信噪比，A为固定周期信号的幅度，f0为固定周期信号的频率，t为布朗粒子的运动时间，ψ为固定周期信号的初相位，m、n均为双稳态势垒实参数，g(t)表示非线性优化模型的输入参数，s为布朗粒子的运动坐标，α为噪声强度，ξ(t)为均值为0的高斯白噪声；将当前子块中的所有像素点各自的像素值作为输入参数输入非线性优化模型中进行优化后，非线性优化模型输出当前子块中的所有像素点各自的信噪比，如果当前子块中的每个像素点的信噪比与事先存储于数据库中的脑肿瘤MRI扫描图像的灰度图像中对应像素点的信噪比的误差小于10％，则确定当前子块优化成功，其中，事先存储于数据库中的脑肿瘤MRI扫描图像的灰度图像的尺寸大小与当前灰度图像的尺寸大小相同；(7)将下一个待处理的尺寸大小为u×u的子块作为当前子块，然后返回步骤(6)继续执行，直至当前初步病灶区域对应的所有尺寸大小为u×u的子块处理完毕，得到对应的最终病灶区域；(8)将当前灰度图像中下一个待处理的初步病灶区域作为当前初步病灶区域，然后返回步骤(5)继续执行，直至当前灰度图像中的所有初步病灶区域处理完毕，得到当前灰度图像中的多个最终病灶区域，至此完成脑肿瘤MRI扫描图像的分割。...

【技术特征摘要】

【专利技术属性】
技术研发人员：汤旭翔，傅均，陈赛，陈柳柳，曹海洋，
申请(专利权)人：浙江工商大学，
类型：发明
国别省市：浙江;33