1. Toy模板网首页
  2. > Toy博客

（数字图像处理MATLAB+Python）第九章图像形态学运算-第三节：二值图像的形态学处理

10月前 作者：快乐江湖 分类：Toy博客 阅读(66) 违法举报

这篇具有很好参考价值的文章主要介绍了（数字图像处理MATLAB+Python）第九章图像形态学运算-第三节：二值图像的形态学处理。希望对大家有所帮助。如果存在错误或未考虑完全的地方，请大家不吝赐教，您也可以点击"举报违法"按钮提交疑问。

一：形态滤波

（1）概述

形态滤波：是一种在数字图像处理中常用的图像处理技术，用于改善图像的质量、提取图像的特定特征或去除图像中的噪声。形态滤波主要基于形态学运算，通过结构元素（也称为模板）对图像进行局部区域的操作，从而改变图像的形状和结构。选择不同形状(如各向同性的圆、十字架、矩形、不同朝向的有向线段等)、不同尺寸的结构元素可以提取图像的不同特征

（2）程序

如下：实现对二值图像的特征提取

（数字图像处理MATLAB+Python）第九章图像形态学运算-第三节：二值图像的形态学处理

matlab实现：

I=imread(‘pattern.jpg’);
BW1=im2bw(I,h);
BW1=1-BW1;
se=strel('square',3); 
BW2=1-imopen(BW1,se);
figure;imshow(BW2);title('矩形块提取');
se45=strel('line',25,45);     
BW3=1-imopen(BW1,se45);
figure;imshow(BW3);title('线段提取');

python实现：

import cv2
import numpy as np
import matplotlib.pyplot as plt

# 读取图像
I = cv2.imread('pattern.jpg', cv2.IMREAD_GRAYSCALE)

# 二值化处理
_, BW1 = cv2.threshold(I, 0, 255, cv2.THRESH_BINARY + cv2.THRESH_OTSU)
BW1 = 255 - BW1

# 矩形块提取
se = cv2.getStructuringElement(cv2.MORPH_RECT, (3, 3))
BW2 = cv2.morphologyEx(BW1, cv2.MORPH_OPEN, se)
BW2 = 255 - BW2

# 显示矩形块提取结果
plt.imshow(BW2, cmap='gray')
plt.title('矩形块提取')
plt.show()

# 线段提取
se45 = cv2.getStructuringElement(cv2.MORPH_RECT, (25, 25))
se45 = cv2.rotate(se45, cv2.ROTATE_90_CLOCKWISE)
BW3 = cv2.morphologyEx(BW1, cv2.MORPH_OPEN, se45)
BW3 = 255 - BW3

# 显示线段提取结果
plt.imshow(BW3, cmap='gray')
plt.title('线段提取')
plt.show()

二：图像的平滑处理

（1）概述

图像的平缓处理：是数字图像处理中的一项重要任务，它有助于减少图像中的噪声和细节，并使图像更加平滑和连续。通过结合开运算和闭运算，可以实现图像的更加全面的平滑处理。首先，使用开运算可以消除细小的噪声和不连续区域，平滑边缘，同时保留较大的结构。接着，使用闭运算可以填补空洞，连接区域，消除缺陷和断裂，使图像更加连续和平滑。对二值图像平滑处理的形态学变换为

Y = ( X ∘ S ) ⋅ S Y = ( X ⋅ S ) ∘ S \begin{array}{l}\boldsymbol{Y}=(\boldsymbol{X} \circ \boldsymbol{S}) \cdot \boldsymbol{S} \\\boldsymbol{Y}=(\boldsymbol{X} \cdot \boldsymbol{S}) \circ \boldsymbol{S}\end{array} Y=(X∘S)⋅SY=(X⋅S)∘S​

例如下图，结构元素 S S S大于所有噪声和噪声块的尺寸

（数字图像处理MATLAB+Python）第九章图像形态学运算-第三节：二值图像的形态学处理

（2）程序

如下

（数字图像处理MATLAB+Python）第九章图像形态学运算-第三节：二值图像的形态学处理

matlab实现：

Image=imread('A.bmp');
BW=im2bw(Image);
SE=strel('square',3);
result1=imclose(imopen(BW,SE),SE); 
figure,imshow(result1);
result2=imopen(imclose(BW,SE),SE);               
figure,imshow(result2);

Python实现：

import cv2
import numpy as np
import matplotlib.pyplot as plt

# 读取图像
Image = cv2.imread('A.bmp', cv2.IMREAD_GRAYSCALE)

# 二值化处理
_, BW = cv2.threshold(Image, 0, 255, cv2.THRESH_BINARY + cv2.THRESH_OTSU)

# 创建结构元素
SE = cv2.getStructuringElement(cv2.MORPH_RECT, (3, 3))

# 先开运算再闭运算
result1 = cv2.morphologyEx(cv2.morphologyEx(BW, cv2.MORPH_OPEN, SE), cv2.MORPH_CLOSE, SE)

# 显示结果1
plt.imshow(result1, cmap='gray')
plt.show()

# 先闭运算再开运算
result2 = cv2.morphologyEx(cv2.morphologyEx(BW, cv2.MORPH_CLOSE, SE), cv2.MORPH_OPEN, SE)

# 显示结果2
plt.imshow(result2, cmap='gray')
plt.show()

三：图像的边缘提取

（1）概述

图像的边缘提取：是数字图像处理中的一项重要任务，它有助于凸显图像中的边缘和轮廓信息，以便于进一步的分析和处理。结合开运算和闭运算可以实现图像的边缘提取。首先，应用闭运算可以填补图像中的空洞和断裂，连接边缘区域，使边缘更加连续。然后，应用开运算可以平滑边缘、去除细小的孤立点和细线条，突出边缘特征。提取物体的轮廓边缘的形态学变换

  • 内边界： Y = X − ( X ⊖ S ) \boldsymbol{Y}=\boldsymbol{X}-(\boldsymbol{X} \ominus \boldsymbol{S}) Y=X−(X⊖S)
  • 外边界： Y = ( X ⊕ S ) − X \boldsymbol{Y}=(\boldsymbol{X} \oplus \boldsymbol{S})-\boldsymbol{X} Y=(X⊕S)−X
  • 形态学梯度： Y = ( X ⊕ S ) − ( X ⊖ S ) \boldsymbol{Y}=(\boldsymbol{X} \oplus \boldsymbol{S})-(\boldsymbol{X} \ominus \boldsymbol{S}) Y=(X⊕S)−(X⊖S)

例如下图，使用结构元素 S S S对图像 X X X进行边缘提取

（数字图像处理MATLAB+Python）第九章图像形态学运算-第三节：二值图像的形态学处理

（数字图像处理MATLAB+Python）第九章图像形态学运算-第三节：二值图像的形态学处理

（2）程序

如下

（数字图像处理MATLAB+Python）第九章图像形态学运算-第三节：二值图像的形态学处理

matlab实现：

Image=imread('menu.bmp');
BW=im2bw(Image);
SE=strel('square',3);
result1=BW-imerode(BW,SE); 
result2=imdilate(BW,SE)-BW; 
result3=imdilate(BW,SE)-imerode(BW,SE);     
figure,imshow(result1);
figure,imshow(result2);
figure,imshow(result3);

python实现：

import cv2
import numpy as np
import matplotlib.pyplot as plt

# 读取图像
Image = cv2.imread('menu.bmp', cv2.IMREAD_GRAYSCALE)

# 二值化处理
_, BW = cv2.threshold(Image, 0, 255, cv2.THRESH_BINARY + cv2.THRESH_OTSU)

# 创建结构元素
SE = cv2.getStructuringElement(cv2.MORPH_RECT, (3, 3))

# 边缘提取结果1
result1 = BW - cv2.erode(BW, SE)

# 显示结果1
plt.imshow(result1, cmap='gray')
plt.show()

# 边缘提取结果2
result2 = cv2.dilate(BW, SE) - BW

# 显示结果2
plt.imshow(result2, cmap='gray')
plt.show()

# 边缘提取结果3
result3 = cv2.dilate(BW, SE) - cv2.erode(BW, SE)

# 显示结果3
plt.imshow(result3, cmap='gray')
plt.show()

四：区域填充

（1）概述

概述：是数字图像处理中的一项常见任务，它有助于填补图像中的空洞或断裂，并连接相邻的区域，使图像中的区域变得连续和完整。结合开运算和闭运算可以实现图像的区域填充。首先，通过应用开运算，可以去除图像中的小型空洞、断裂和孤立点，使得区域边缘更加连续。接着，再应用闭运算，可以填补空洞、连接区域，使得图像中的区域变得完整和连续

X k = ( X k − 1 ⨁ S ) ∩ A C , k = 1 , 2 , 3 , ⋯ X_{k}=\left(X_{k-1} \bigoplus S\right) \cap A^{C}, \quad k=1,2,3, \cdots Xk​=(Xk−1​⨁S)∩AC,k=1,2,3,⋯

如下图，边界点用灰色表示，赋值为1；所有非边界点是白色部分，赋值为0

（数字图像处理MATLAB+Python）第九章图像形态学运算-第三节：二值图像的形态学处理

（2）程序

如下

（数字图像处理MATLAB+Python）第九章图像形态学运算-第三节：二值图像的形态学处理

matlab实现：

Image=imread('coin.bmp');
BW=im2bw(Image);
imshow(BW); title('二值图像');
result1=imfill(BW,'holes');
figure,imshow(result1);title('二值图像的区域填充');

python实现：

import cv2
import numpy as np
import matplotlib.pyplot as plt

# 读取图像
Image = cv2.imread('coin.bmp', cv2.IMREAD_GRAYSCALE)

# 二值化处理
_, BW = cv2.threshold(Image, 0, 255, cv2.THRESH_BINARY + cv2.THRESH_OTSU)

# 显示二值图像
plt.imshow(BW, cmap='gray')
plt.title('二值图像')
plt.show()

# 区域填充
result1 = cv2.fillHoles(BW)

# 显示区域填充结果
plt.imshow(result1, cmap='gray')
plt.title('二值图像的区域填充')
plt.show()

五：目标探测（击中与否变换）

（1）概述

目标探测（击中与否变换）：击中与否变换的原理利用腐蚀特性——腐蚀的过程相当于对可以填入结构元素的位置作标记的过程。目标检测既要检测目标内部，也要检测外部，即在一次运算中要同时捕获内外标记。因此，进行目标探测，需要采用两个结构元素 S 1 S_{1} S1​、 S 2 S_{2} S2​，构成一个结构元素对 S = ( S 1 , S 2 ) S=(S_{1},S_{2}) S=(S1​,S2​)

  • S 1 S_{1} S1​：探测目标内部
  • S 2 S_{2} S2​：探测目标外部

则图像 X X X用结构元素 S S S进行击中与否变换，记为

X ∗ S = ( X ⊖ S 1 ) ∩ ( X C ⊖ S 2 ) X ∗ S = ( X ⊖ S 1 ) ∩ ( X C ⨁ S ^ 2 ) X ∗ S = { x ∣ S 1 + x ⊆ X  and  S 2 + x ⊆ X C } \begin{array}{l}\boldsymbol{X} * \boldsymbol{S}=\left(\boldsymbol{X} \ominus \boldsymbol{S}_{\mathbf{1}}\right) \cap\left(\boldsymbol{X}^{\boldsymbol{C}} \ominus \boldsymbol{S}_{\mathbf{2}}\right) \\\boldsymbol{X} * \boldsymbol{S}=\left(\boldsymbol{X} \ominus \boldsymbol{S}_{\mathbf{1}}\right) \cap\left(\boldsymbol{X}^{\boldsymbol{C}} \bigoplus \widehat{\boldsymbol{S}}_{\mathbf{2}}\right) \\\boldsymbol{X} * \boldsymbol{S}=\left\{\boldsymbol{x} \mid \boldsymbol{S}_{\mathbf{1}}+\boldsymbol{x} \subseteq \boldsymbol{X} \quad \text { and } \quad \boldsymbol{S}_{\mathbf{2}}+\boldsymbol{x} \subseteq \boldsymbol{X}^{\boldsymbol{C}}\right\}\end{array} X∗S=(X⊖S1​)∩(XC⊖S2​)X∗S=(X⊖S1​)∩(XC⨁S 2​)X∗S={x∣S1​+x⊆X and S2​+x⊆XC}​

当且仅当结构元素 S 1 S_{1} S1​平移到某一点可填入集合 X X X的内部，结构元素 S 2 S_{2} S2​平移到该点可填入集合 的外部时，该点才在击中击不中变换的输出中

例如下图，(a)为由四个物体：矩形、小方形、大方形、带有小凸出部分的大方形组成的图像 X X X；(b)为结构元素对 S = ( S 1 , S 2 ) S=(S_{1},S_{2}) S=(S1​,S2​)。要求通过击中与否运算，能正确识别方形

（数字图像处理MATLAB+Python）第九章图像形态学运算-第三节：二值图像的形态学处理

效果如下

（数字图像处理MATLAB+Python）第九章图像形态学运算-第三节：二值图像的形态学处理

（2）程序

如下

（数字图像处理MATLAB+Python）第九章图像形态学运算-第三节：二值图像的形态学处理

matlab实现：

Image=imread('test.bmp');
BW=im2bw(Image);                 
interval=[-1  -1  -1  -1  -1
                  -1  -1  -1  -1  -1
                  -1  -1  1   1   1
                  -1  -1  1   1   1
                  -1  -1  1   1   1];      %定义结构元素对
result=bwhitmiss(BW,interval);   %击中击不中
figure,imshow(result); title('击中与否变换结果');

python实现：

import cv2
import numpy as np
import matplotlib.pyplot as plt

# 读取图像
Image = cv2.imread('coin.bmp', cv2.IMREAD_GRAYSCALE)

# 二值化处理
_, BW = cv2.threshold(Image, 0, 255, cv2.THRESH_BINARY + cv2.THRESH_OTSU)

# 显示二值图像
plt.imshow(BW, cmap='gray')
plt.title('二值图像')
plt.show()

# 区域填充
result1 = cv2.fillHoles(BW)

# 显示区域填充结果
plt.imshow(result1, cmap='gray')
plt.title('二值图像的区域填充')
plt.show()

六：细化

（1）概述

细化：骨架化结构是目标图像的重要拓扑描述。对目标图像进行细化处理，就是求图像的中央骨架的过程，是将图像上的文字、曲线、直线等几何元素的线条沿着其中心轴线将其细化成一个像素宽的线条的处理过程。基于数学形态学变换的细化算法为

X ⊙ S = X − ( X ∗ S ) \boldsymbol{X} \odot \boldsymbol{S}=\boldsymbol{X}-(\boldsymbol{X} * \boldsymbol{S}) X⊙S=X−(X∗S)

可见，细化实际上为从集合 X X X中去掉被结构元素 S S S击中的结果。具体采用的细化方法为

X 1 = X ⊙ S , X 2 = X 1 ⊙ S , ⋯ ⋯ X n = X n − 1 ⊙ S \boldsymbol{X}_{\mathbf{1}}=\boldsymbol{X} \odot \boldsymbol{S}, \boldsymbol{X}_{\mathbf{2}}=\boldsymbol{X}_{\mathbf{1}} \odot \boldsymbol{S}, \cdots \cdots \boldsymbol{X}_{n}=\boldsymbol{X}_{n-1} \odot \boldsymbol{S} X1​=X⊙S,X2​=X1​⊙S,⋯⋯Xn​=Xn−1​⊙S

如下示例

（数字图像处理MATLAB+Python）第九章图像形态学运算-第三节：二值图像的形态学处理

（数字图像处理MATLAB+Python）第九章图像形态学运算-第三节：二值图像的形态学处理

（2）程序

如下

（数字图像处理MATLAB+Python）第九章图像形态学运算-第三节：二值图像的形态学处理

matlab实现：

Image=imread('menu.bmp');
BW=im2bw(Image); 
result1=bwmorph(BW,'thin',1); 
result2=bwmorph(BW,'thin',Inf); 
figure,imshow(result1);title('细化一次');
figure,imshow(result2);title('细化至只有一个像素宽');

python实现：文章来源地址https://www.toymoban.com/news/detail-475615.html

Image=imread('menu.bmp');
BW=im2bw(Image); 
result1=bwmorph(BW,'thin',1); 
result2=bwmorph(BW,'thin',Inf); 
figure,imshow(result1);title('细化一次');
figure,imshow(result2);title('细化至只有一个像素宽');

到了这里，关于（数字图像处理MATLAB+Python）第九章图像形态学运算-第三节：二值图像的形态学处理的文章就介绍完了。如果您还想了解更多内容，请在右上角搜索TOY模板网以前的文章或继续浏览下面的相关文章，希望大家以后多多支持TOY模板网！

本文来自互联网用户投稿，该文观点仅代表作者本人，不代表本站立场。本站仅提供信息存储空间服务，不拥有所有权，不承担相关法律责任。如若转载，请注明出处： 如若内容造成侵权/违法违规/事实不符，请点击违法举报进行投诉反馈，一经查实，立即删除！

分享到：
领支付宝红包 赞助服务器费用

相关文章

觉得文章有用就打赏一下文章作者

支付宝扫一扫打赏

博客赞助

微信扫一扫打赏

请作者喝杯咖啡吧~博客赞助

支付宝扫一扫领取红包，优惠每天领

二维码1

领取红包

二维码2

领红包