美文网首页
灰度形态学的基本应用

灰度形态学的基本应用

作者: 此间不留白 | 来源:发表于2020-08-06 22:10 被阅读0次

    上海交通大学 医学图像处理技术

    形态学边缘检测(Edge Detection)

    用灰度形态学实现的边缘检测,按照其实现方式的不同,可以分为以下三种类型,分别是:

    • Standard Edge Detection
      Edge_S(F,K) = (F \oplus K)-(F \ominus K)
    • External Edge Detection
      Edge_E(F,K) = (F \oplus K) -F
    • Internal Edge Detection
      Edge_I(F,K) = F -(F \ominus K)

    如下图所示,形象地说明了三种边缘检测算法的基本工作原理。


    如上图所示,有图像F,和半径为r的圆形SE,首先对用SE对原始图像进行膨胀,原始图像在原来的尺寸上,半径增加r/2,再用SE对原始图像进行腐蚀,原始图像的的半径减少r/2,经过膨胀和腐蚀得到的图像相减,得到Strandard Edge; 膨胀后的图像减原始图像,得到External Edge;而原始图像,减去腐蚀后的图像,得到Internal Edge.

    Egde Detection 的应用实例 MRI图像进行Edge Detection 的应用实例

    形态学梯度(Morphological Gradient)

    针对以上三种不同的Edge Detection,有三种不同的Morphological Gradient,用数学公式的描述,如下所示:

    • Gradient_S(F) = \frac{1}{2}[(F \oplus K)-(F \ominus K)]
    • Gradient_I(F) = \frac{1}{2}[ (F \oplus K) -F]
    • Gradient_E(F) = \frac{1}{2} [F -(F \ominus K)]

    形态学梯度的计算公式和形态学边缘检测的计算公式非常类似,对于形态学梯度的理解,可以用以下图片表示:



    如上图所示,用一条曲线表示原始图像的像素变化,可以看出,膨胀是对原始图像做了一个取上界的平滑,而腐蚀是对原始图像,做了一个取下限的平滑,膨胀后的图像减去腐蚀后的图像就会得到一个梯度的变换的曲线,可以看出,当图像的曲线在某一特定范围内剧烈变换时,梯度曲线的变化也越剧烈,相当于对原始图像实现了一个求导的操作。

    形态学梯度的应用,如下图所示,可以看出,利用形态学梯度,可以很好地区分图像不同区域的灰度变化,所以,在图像分割领域,形态学梯度有很好的应用。

    形态学平滑(Morphological Smoothing)

    利用形态学开运算和闭运算可以实现图像的形态学平滑(也可以称之为图像除噪)处理,想象一下,假如有一块高低不平的田地,为了使这块平整这块土地,你需要挖掉高处的土,填到田地的低洼处,而闭运算就相当于挖掉高处的土,开运算就相当于将在高处挖掉的图填至低处。所以,按照这个原理,利用开运算和闭运算就可以实现图像的平滑,其计算公式如下:
    MSmooth(F) = (F \circ_g K) \bullet_g K
    如下图所示,是图像平滑的一个具体示意图


    如下图所示,是形态学平滑的效果图,可以看出,形态学平滑消除图像中的一些噪点。

    利用开运算和闭运算实现的形态学平滑,具有计算量大,复杂度较高等缺点,俄日了应对这一问题,可以直接使用膨胀和腐蚀来实现这一运算,其实现公式如下所示:

    DSmooth(F) = \frac{1}{2} = [(F \oplus_g K)-(F \ominus_g K)]

    利用膨胀和腐蚀实现的形态学smooth,如下图所示:


    Top-hat Transform

    Top-hat Transform算法在中文里被称为“高帽算法”,简称为TT,Top-hat算法。提取不光滑背景中的图像特征时,Top-hat算法是一种非常有效的工具。按照Top-hat算法的实现方式,TT算法可以分为以下两类:

    • White Top-hat Transform(WTT)
      T_i = F-F \circ_gK_i

    • Black Top-hat Transform(BTT)
      T_i = F \bullet_gr_iK-F
      TT 算法在三维拓扑中的直观描述如下图所示:


      如下图所示的信号曲线,通过灰度开运算,SE从下往上进行fit时,有一些尖端无法被fit到,所以会被抹去,得到原始信号经由灰度开运算得到的信号,原始信号通过和灰度开运算得到的信号相减,得到原始信号的Top-hat

    TT算法的一些应用

    如下图所示,通过BTT和WTT两种算法,得到效果图,可以看出来,通过BTT算法,得到了图像的一些“白区”,而WTT算法得到了图像的一些“纹路”

    如下图所示,是一张卫星图片,通过WTT,二值化,区域重建等算法等算法的组合,可以求取图像中飞机跑道的长度。


    2.JPG

    如下图所示的一张化石图片中,需要检测化石图片中线虫的位置,可以通过BTT算法,除噪,二值化,Skeleton,消除图像中较短的结构,图像重建等一系列算法的组合,完成化石中线虫的检测。


    TT算法的缺点

    TT算法虽然被在形态学上有广泛的应用,但有以下两个显著的缺陷:

    1. 难以区分图像中复杂的ROI
    2. TT算法在应用时,很难确定SE的大小。

    Difference Transform Top-hat(DTT)

    针对TT算法的缺陷,科研人员提出了一种新的TT算法,该算法即是DTT算法,DTT算法,采用了一系列形状相同,尺寸不同的SE,求取不同尺寸大小的SE下TT算法的差值,其数学公式如下所
    示:
    F_i = |T_i-T_{i-1}|_B+F^{'}_{i-1}, \ \ \ \ \ \ \ F^{'}_{i}= \bigcup_{1<j<i}F_j;F_i=\oslash

    DTT算法的工作原理是:

    • 找到一系列的TT(通过尺寸不同,形状相同的SE计算得到)
    • 计算不同TT之间的差值
    • 通过阈值的方式将求取同一阈值下的TT
    • 通过或运算,将同一阈值条件下的TT相加

    DTT算法的优点:
    与TT算法相比,DTT算法具有以下优点:

    1. 能够有效求取梯度变换较快的区域
    2. 能够将图像中的物体按照其尺寸进行排序
    DTT算法的原理示意图

    如下图所示,TT和DTT算法效果比较,右上角的曲线反映了图像的梯度变化,可以看出,DTT能够消除掉图像中梯度较小的区域,能够更加有效的提取图像中梯度变换更快的区域。


    DTT算法的应用

    如下图所示的杂志封面处理中,可以看出,图像中间杂志的封面,文字和杂志中的图片易于区分,非常适合处理,而其他杂志封面中文字和图片混在一起,处理难度较大,通过DTT算法可以有效解决这一问题



    Utlimate Erosion(UE)

    想想一个剥洋葱的过程,通过将一个洋葱一层一层剥下去,会得到一个芯,同剥洋葱一样,通过对一副图像进行不间断的Erosion,直至图像中的一些点消失之前停止,会得到图像的一些种子点,如下图所示,有了种子点之后,也很容易地可以实现对图像的重建(比如通过膨胀算法等)。



    通过Utlimate Erosion可以很容易地实现对图像中一些相连的部分进行分割,如下图所示,


    Geodesic Influence

    Geodesic Influence(GI)算法与UE算法相对应,只不过是一个Recursive Dilation的过程,这种膨胀算法也称之为是受限制的膨胀算法(条件膨胀算法)。

    GI算法的主要工作原理是,有一些种子点,通过GI算法,会得到经由这些种子点膨胀后的图像,这些以不同种子点为中心不断膨胀的过程中,种子附近的像素分配原则是“先到先得”,即边界处的像素点,按照不同种子膨胀后的边界的接触顺序进行分配。最后,通过GI算法,可以非常有效地实现对复杂图像的区域分割,如下图所示:

    相关文章

      网友评论

          本文标题:灰度形态学的基本应用

          本文链接:https://www.haomeiwen.com/subject/sakgrktx.html