[发明专利]一种基于数学形态学的医学图像的分割算法

说明书

本发明公开了医学图像处理技术领域的一种基于数学形态学的医学图像的分割算法，其特征在于：该基于数学形态学的医学图像的分割算法包括如下步骤：S1：医学图像形态重建；S2：多尺度医学图像形态梯度计算；S3：医学图像形态梯度修正；S4：医学图像分水岭变换及区域合并，该种基于数学形态学的医学图像的分割算法，设计合理，实用性强，能够有效的除去噪音图像分隔的干扰，而且可以消除过分隔区域，同时具有很强的区域轮廓定位能力，大大提高了图像分割效果，使得医生或者医学研究人员对医学图像上显示的器官和组织的更好的辨认。

技术领域

本发明涉及医学图像处理技术领域，具体为一种基于数学形态学的医学图像的分割算法。

背景技术

医学图像分割是现代医学图像处理中的关键技术，它是正常组织和病变组织的三维重建、定量分析等后续操作的基础，也是临床医学应用的瓶颈，医学CT图像分割的主要目的是提取感兴趣的器官、组织等部分以用于临床诊断和病理学研究等诸多方面，近年来，将一些特定的新兴理论应用到医学图像分割中是医学图像处理的一个研究热点，数学形态学是一门新兴理论，其典型的图像分割算法是形态学分水岭变换，但是传统的形态学分水岭变换存在严重的过分割现象，而且受到外部噪音的影响过大，导致医学图像的医学图像分割效果降低，影响医生对医学图像上显示的器官和组织的辨认，导致医生产生误诊等情况，为此，我们提出了一种基于数学形态学的医学图像的分割算法。

发明内容

本发明的目的在于提供一种基于数学形态学的医学图像的分割算法，以解决上述背景技术中提出的传统的形态学分水岭变换存在严重的过分割现象，而且受到外部噪音的影响过大，导致医学图像的医学图像分割效果降低，影响医生对医学图像上显示的器官和组织的辨认，导致医生产生误诊等情况的问题。

为实现上述目的，本发明提供如下技术方案：一种基于数学形态学的医学图像的分割算法，该基于数学形态学的医学图像的分割算法包括如下步骤：

S1：医学图像形态重建：首先对医学图像进行一次开重建，具体公式为：(Y^～_N(f)＝Y^(rec)(Y_N(f)，f))，当对医学图像开重建完成后，再对开重建后的医学图像进行一次闭重建，具体公式为(T^～_N(f)＝Y^(rec)(T_N(f)，f))，式中Y，T分别表示形态学中的开运算和闭运算，f表示要重建的图像，(rce)表示重建运算，N表示运算使用的结构元素的尺寸；

S2：多尺度医学图像形态梯度计算：假设B_i(0≤i≤n)为一组正方形的结构元素，B_i的大小为(2i+1)×(2i+1)个像素点，则多尺度梯度的计算公式为式中f为原始图像，n为运算使用的结构元素的尺寸；

S3：医学图像形态梯度修正：假设梯度等级为l时，相应的结构半径r(l)为：0≤l≤max(l),0<a<1，式中r₀为最大半径，max(l)为梯度最大值，a为调节因子，设对应的层面经修正后转变为g_l'，则：g_l'＝T(l,r(l)),0≤l≤max(l)，将每一层面的输出进行合并，即确定修正后的梯度g'：