数字图像处理第九章——形态学图像处理

数字图像处理—形态学图像处理

同样的，暂时对书上已经写得很清楚的知识点不再重复赘述，主要做一些总结，思考以及知识点的梳理和扩展。

（一）预备知识

介绍一下形态学中的一些基本概念。

1. 用数学形态学（也称图像代数）表示以形态为基础对图像进行分析的数学工具
2. 基本思想是用具有一定形态的结构元素去度量和提取图像中的对应形状以达到对图像分析识别的目的
3. 形态学图像处理的数学基础和所用语言是集合论
4. 形态学图像处理的应用可以简化图像数据， 保持它们基本的形状特性，并除去不相干的结 构
5. 形态学图像处理的基本运算有4个：膨胀、 腐蚀、开操作和闭操作

1.1 集合理论中的基本概念

介绍一下比较陌生的几个概念，其他的看书就好：

1. 所有像素坐标的集合均不属于集合A，记为$A^c$，由下式给出:
   
   这个集合称为集合A的补集。

2. 集合B的反射，定义为：
   
   即关于原集合原点对称 .

3. 集合A平移到点z=(z1,z2),表示为(A)z，定义为:
   

1.2 二值图像、集合及逻辑算子

二值图像：

二值图像(Binary Image)，按名字来理解只有两个值，0和1，0代表黑，1代表白，或者说0表示背景，而1表示前景。其保存也相对简单，每个像素只需要1Bit就可以完整存储信息。如果把每个像素看成随机变量，一共有N个像素，那么二值图有2的N次方种变化，而8位灰度图有255的N次方种变化，8为三通道RGB图像有255255255的N次方种变化。也就是说同样尺寸的图像，二值图像保存的信息更少。二值图像（binary image），即图像上的每一个像素只有两种可能的取值或灰度等级状态，人们经常用黑白、B&W、单色图像表示二值图像。

二值图像集合：

如果A和B是二值图像，那么C=A∪B仍是二值图像。这里，如 果 A 和B中相应的像素不是前景像素就是背景像素，那么 C中的这个像素就是前景像素。以第一种观点，函数 C由下式给出:

另一方面，运用集合的观点，C由下式给出：

集合运算：

1. A为图像集合，B为结构元素（集合）。
2. 数学形态学运算时B对A进行操作。
3. 结构元素要有1个原点（即结构元素参与形态学运算的参考点），可以是中心像素，原则上可选任何像素。
   注意：原点可以包含在结构元素中，也可以不包含在结构元素中，但运算的结果常不相同。

编码：

f = imread('D:\数字图像处理\第九章学习\Fig0903(a).tif');
g = imread('D:\数字图像处理\第九章学习\Fig0903(b).tif');
subplot(2,3,1), imshow(f);title('(a)二值图像 A:');
subplot(2,3,2), imshow(g);title('(b)二值图像 B:');
subplot(2,3,3), imshow(~f);title('(c)A的补集~A:');
subplot(2,3,4), imshow(f|g);title('(d) A和B的并集 A|B:');
subplot(2,3,5), imshow(f&g);title('(e)A和B的交集 A & B:');
subplot(2,3,6), imshow(f&~g);title('(f)A和B的差集 A&~B');

代码运行效果如下：

分析：

图像(d)是 “ UTK”和 “ GT” 图像的并集，包括来自两幅图像的所有前景像素。相反，两幅图像的交集(图(e))显示了字母 “ UTK”和 “ GT”中重叠的像素。最后，集合的差集图像(图(f))显示了 “ UTK”中除去 “ GT” 像素后的字母。

（二）膨胀和腐蚀

2.1 膨胀

膨胀：膨胀是在二值图像中“加长”或“变粗”的操作。这种特殊的方式和变粗的程度由一个称为结构元素的集合控制。（实际就是将结构元素的原点与二值图像中的1重叠，将二值图像中重叠部分不是1的值变为1，完成膨胀）。

公式：

A和B是两个集合，A被B膨胀定义为：

公式解释：

1. B的反射进行平移与A的交集不为空。
2. B的反射：相对于自身原点的映象。
3. B的平移：对B的反射进行位移

图解：

      

(a)集合A    (b)结构元素B (黑色为原点所在)

      

(c)结构元素B的映像    (d)图中两种阴影部分（深色为扩大的部分）合起来为A+B

注意：

1. 膨胀运算只要求结构元素的原点在目标图像的内部平移，换句话说，当结构元素在目标图像上平移时，允许结构元素中的非原点像素超出目标图像的范围。
2. 膨胀运算具有扩大图像和填充图像中比结果元素小的成分的作用，因此在实际应用中可以利用膨胀运算连接相邻物体和填充图像中的小孔和狭窄的缝隙。

膨胀举例：

膨胀函数：

D = imdilate(A,B)

图像膨胀的应用：桥接文字裂缝

编码：

A = imread('D:\数字图像处理\第九章学习\Fig0906(a).tif');
B = [0 1 0; 1 1 1; 0 1 0];   %指定结构元素由0和1组成的矩阵
A2 = imdilate(A, B);    %二值图像
subplot(1,2,1), imshow(A);title('(a)包括断开文本的输入图像:');
subplot(1,2,2), imshow(A2);title('(b)膨胀后图像:');

图片中字体的加粗,且填充了字母中的小孔和狭窄的缝隙。

2.2 结构元的分解

公式：

公式理解：

B膨胀A等同于B1先膨胀A，再用B2膨胀之前的结果。

举例：

下面是由1组成的5x5数组的膨胀：

这个结构元能够分解为值为 1 的 5 元素行矩阵和值为 1 的 5 元素列矩阵:

分析：

在原结构元中，元素个数为 25; 但在行列分解后，总元素数目仅为 10。这意味着首先用 行结构元膨胀，再用列结构元膨胀，能够比 5x5 的数组膨胀快 2.5 倍。在实践中，速度的增长稍微慢一些，因为在每个膨胀运算中总有些其他开销。然而，由分解执行获得的速度方面的增 长仍然有很大意义