基于傅里叶变换的医学图像分割算法
【技术领域】
[0001] 本发明属于图像处理领域，具体涉及一种基于傅里叶变换的医学图像分割算法。
【背景技术】
[0002] 手动圈定病灶和使用图像分析软件分割图像是现有的两大类医学图像分割方法。 手动方法大多用在比较复杂的场合，严重依赖操作者的病灶判断经验，因此并不适合大部 分缺乏相关专业领域经验的操作人员去使用。并且，手动方法分割图像的效率也很低，不适 合对大批量图像的分割作业，因此其普及和推广受到限制。使用图像分析软件分割图像，存 在对图像分割的精确度和准确性不足的问题，易产生误判，而误判对于病人及家属来说是 不可接受的，不利于病人康复，也是大部分医患矛盾的根源。
[0003] 图像分割是指利用图像信息中某些特征，从图像中提取用户感兴趣目标。对待分 割图像进行信息分析，对重要的信息进行特征提取。在进行初始化分割之前，先对图像中各 像素提取特征，特征提取后，可以得到一个相应的特征图像，这个特征图像包含三个通道， 每个通道对应着一组特征值。但是对不同场景的图像，其特征图像具有不同的特性，例如具 有强烈纹理特性的图像，其方向特征会较其他两个特征有更大的对比度；而对于自然场景 图像，其颜色和方向的对比度会比较大。概括地说，给定某种一致性属性准则（度量）P，将 图像X正确地划分成为互不交迭的区域集丨Si，的过程称之为分割，正确分割应 满足下列条件：
[0008] 用作图像分割的度量准则不是唯一的，它与应用场景图像及应用目的有关，用于 图像分割的场景图像特征信息有亮度、色彩、纹理、结构、温度、频谱、运动、形状、位置、梯度 和模型等。
[0009] 傅里叶变换是一种正交变换，其是将时域信号分解为不同频率的正弦信号或者余 弦函数叠加之和，在一维信号处理中得到了广泛的应用，目前很多学者将该方法应用到图 像处理领域，提出图像特征提取分割方法，可以较好地对图像中的相位、幅度、高低频率以 及图像分割增强等进行处理。但这类图像特征提取分割难以无法完成对图像关键特征信息 的适度提取，同时其分割算法复杂，造成分割的图像的分辨率低，清晰度不高。
【发明内容】
[0010] 本发明所要解决的技术问题是：针对现有技术的不足，提供一种对医学图像分割 精度高的基于傅里叶变换的医学图像分割算法。
[0011] 本发明解决上述技术问题所采用的技术方案为：基于傅里叶变换的医学图像分割 算法，包括以下步骤：
[0012] (1)获取一幅来源于X射线机、CT或核磁共振诊断设备的待分割的医学扫描图像， 然后将该待分割的医学扫描图像转换成灰度图像，再采用傅里叶变换算法对灰度图像进行 分析，得到初始感兴趣的轮廓{A1，A2，A3, 一,Anh其中A1，A2，A3, 一,An定义为构成初始感 兴趣的轮廓的所有轮廓点；
[0013] ⑵根据初始感兴趣的轮廓{A1，A2, A3, ···，!!}经包络得到灰度图像中的感兴趣 区域，根据该感兴趣区域的外形选取一个能够将该感兴趣区域包含在内的方形区域，定义 该方形区域的宽度为W、高度为H，如果WXH能够被uXu整除，则将该方形区域定义为当前 灰度图像，然后直接将当前灰度图像划分@
个互不重叠的尺寸大小为uXu的子块； 如果WXH不能够被uXu整除，则扩展该方形区域使其尺寸大小能够被uXu整除，将扩展 后的方形区域定义为当前灰度图像，然后将当前灰度图像划分成
个互不重叠的尺 寸大小为uXu的子块，其中，W'和H'对应表示扩展后的方形区域的宽度和高度，W' == W 且 Η' >H 或 W' >W 且 H' = = H 或 W' >W 且 H' >H，u 取 4、5 或 6 ;
[0014] (3)采用区域生长法对当前灰度图像进行分割，得到当前灰度图像中的多个初步 目标区域；
[0015] (4)将当前灰度图像中当前待处理的初步目标区域定义为当前初步目标区域；
[0016] (5)从当前灰度图像中提取出与当前初步目标区域对应的所有尺寸大小为uXu 的子块，按序对当前初步目标区域对应的每个尺寸大小为uXu的子块进行处理，将当前待 处理的尺寸大小为uXu的子块定义为当前子块；
[0017] (6)将当前子块中的所有像素点各自的像素值作为输入参数输入非线性优化模型 中进行优化，该非线性优化模型为：
[0019] 其中：^表示非线性优化模型输出的信噪比，A为固定周期信号的幅度，f。为固定 周期信号的频率，t为布朗粒子的运动时间，Φ为固定周期信号的初相位，m、η均为双稳态 势皇实参数，g(t)表示非线性优化模型的输入参数，s为布朗粒子的运动坐标，α为噪声强 度，ξ (t)为均值为0的高斯白噪声；
[0020] 将当前子块中的所有像素点各自的像素值作为输入参数输入非线性优化模型中 进行优化后，非线性优化模型输出当前子块中的所有像素点各自的信噪比，如果当前子块 中的每个像素点的信噪比与事先存储于数据库中的医学扫描图像的灰度图像中对应像素 点的信噪比的误差小于10%，则确定当前子块优化成功，其中，事先存储于数据库中的医学 扫描图像的灰度图像的尺寸大小与当前灰度图像的尺寸大小相同；
[0021] (7)将下一个待处理的尺寸大小为uXu的子块作为当前子块，然后返回步骤（6) 继续执行，直至当前初步目标区域对应的所有尺寸大小为uXu的子块处理完毕，得到对应 的最终目标区域；
[0022] (8)将当前灰度图像中下一个待处理的初步目标区域作为当前初步目标区域，然 后返回步骤（5)继续执行，直至当前灰度图像中的所有初步目标区域处理完毕，得到当前 灰度图像中的多个最终目标区域，至此完成医学扫描图像的分割。
[0023] 作为优选，步骤（6)中事先存储于数据库中的医学扫描图像的灰度图像中的每个 像素点的信噪比的具体获取过程为：
[0024] (6. 1)选取一幅来源于X射线机、CT或核磁共振诊断设备的具有代表性的医学扫 描图像，然后将该医学扫描图像转换为灰度图像；
[0025] (6. 2)定义该灰度图像的宽度为W、高度为H，如果WXH能够被uXu整除，则将该 灰度图像作为待处理灰度图像，然后直接将待处理灰度图像划分试
个互不重叠的尺 寸大小为uXu的子块；如果WXH不能够被uXu整除，则扩展该灰度图像使其尺寸大小能 够被uXu整除，将扩展后的灰度图像作为待处理灰度图像，然后将待处理灰度图像划分成
t互不重叠的尺寸大小为uXu的子块，其中，W'和H'对应表示扩展后的灰度图像 的宽度和高度，W' ==W且H'>H或W'>W且Η' ==H或W'>W且H'>H，u取4、5或6;
[0026] (6. 3)手动圈定待处理灰度图像中的每个目标区域；
[0027] (6.4)计算待处理灰度图像中的每个目标区域对应的每个尺寸大小为uXu的子 块中的每个像素点的信噪比。
[0028] 与现有技术相比，本发明的优点在于：本发明公开的基于傅里叶变换的医学图像 分割算法，采用傅里叶变换算法对灰度图像进行分析，获得感兴趣的灰度图像区域，并将感 兴趣的灰度图像区域划分为互不重叠的子块，采用区域生长法对当前灰度图像进行分割， 综合考虑了所有像素，提高了像素分割的准确性，避免了现有分割方法的不足；采用本发明 分割算法对医学图像进行分割，分割效果明显，有利于提高对医学图像及病例判断的准确 性，尤其适用于噪声污染较严重和边缘化较模糊的医学图像边缘的准确确定。
【附图说明】
[0029] 图1为经本发明分割算法分割前的鼻咽癌MRI图像；
[0030] 图2为经本发明分割算法分割后的鼻咽癌MRI图像。
【具体实施方式】
[0031] 以下结合附图实施例对本发明作进一步详细描述。
[0032] 对于图1所示的鼻咽癌MRI图像，采用本发明的基于傅里叶变换的医学图像分割 算法进行分割，包括以下步骤：
[0033] (1)获取图1所示的鼻咽癌MRI图像的扫描图像，然后将该扫描图像转换成灰度 图像，再采用傅里叶变换算法对灰度图像进行分析，得到初始感兴趣的轮廓{A1，A2，A3,… ，An}，其中Al, A2, A3,…，An定义为构成初始感兴趣的轮廓的所有轮廓点；
[0034] (2)根据初始感兴趣的轮廓{A1，A2, A3, ···，!!}经包络得到灰度图像中的感兴趣 区域，根据该感兴趣区域的外形选取一个能够将该感兴趣区域包含在内的方形区域，定义 该方形区域的宽度为W、高度为H，如果WXH能够被5 X 5整除，则将该方形区域定义为当前 灰度图像，然后直接将当前灰度图像划分姑
个互不重叠的尺寸大小为5X5的子块； 如果WXH不能够被5 X 5整除，则扩展该方形区域使其尺寸大小能够被5 X 5整除，将扩展 后的方形区域定义为当前灰度图像，然后将当前灰度图像划分试
t互不重叠的尺 寸大小为5X5的子块，其中，W'和H'对应表示扩展后的方形区域的宽度和高度，W' == W 且 Η' >H 或 W' >W 且 H' = = H 或 W' >W 且 Η' >H ;
[0035] (3)采用区域生长法对当前灰度图像进行分割，得到当前灰度图像中的多个初步 目标区域；
[0036] (4)将当前灰度图像