形态学操作—边界提取

  图像形态学是一种用于图像处理的数学理论,它涉及对图像中形状和结构的分析与处理。图像形态学操作中的边界提取是一种常见的处理方式,可以用来检测图像中物体的轮廓或边缘。

原理

  边界提取的主要原理是使用形态学操作中的膨胀(Dilation)和腐蚀(Erosion)。膨胀操作会使物体区域膨胀,扩大物体的边界;腐蚀操作会使物体区域收缩,减小物体的边界。通过将膨胀后的图像与腐蚀后的图像进行差异,可以得到物体的边界。

作用和适用场景

  边界提取常用于图像分割、对象检测和特征提取等领域。它可以帮助定位图像中的物体边界,用于识别形状、计算物体的面积和周长等应用场景。

数学公式

  在图像形态学中,膨胀和腐蚀操作可以使用结构元素(structuring element)对图像进行操作。结构元素通常是一个小的二值图像(例如矩形、圆形或十字形),它在进行膨胀或腐蚀时与原始图像进行卷积操作。
  膨胀操作 Dilate ( I , S E ) \text{Dilate}(I, SE) Dilate(I,SE) 可以表示为:
( Dilate ( I , S E ) ) ( x , y ) = max ⁡ ( i , j ) ∈ S E I ( x − i , y − j ) (\text{Dilate}(I, SE))(x, y) = \max_{(i,j) \in SE} I(x - i, y - j) (Dilate(I,SE))(x,y)=(i,j)SEmaxI(xi,yj)
  腐蚀操作 Erode ( I , S E ) \text{Erode}(I, SE) Erode(I,SE) 可以表示为:
( Erode ( I , S E ) ) ( x , y ) = min ⁡ ( i , j ) ∈ S E I ( x + i , y + j ) (\text{Erode}(I, SE))(x, y) = \min_{(i,j) \in SE} I(x + i, y + j) (Erode(I,SE))(x,y)=(i,j)SEminI(x+i,y+j)
  其中, I I I 是输入的二值图像, S E SE SE 是结构元素。

代码示例(使用Python和OpenCV)

import cv2
import numpy as np

def show_images(image):
    cv2.namedWindow('image',cv2.WINDOW_KEEPRATIO)
    cv2.imshow('image',image)
    cv2.waitKey()
    cv2.destroyAllWindows()

def Boundary_Extracted(image):
    # 设定结构元素
    kernel = np.ones((5, 5), np.uint8)
    # 膨胀操作
    dilated = cv2.dilate(image, kernel, iterations=1)
    # 腐蚀操作
    eroded = cv2.erode(image, kernel, iterations=1)
    # 边界提取
    boundary = cv2.absdiff(dilated, eroded)
    return boundary

if __name__ == '__main__':
    # 读取图像
    img = cv2.imread('cat-dog.png', flags=0)
    re_img=Boundary_Extracted(img)
    # top_row = np.hstack((img, re_img[0]))
    # bottom_row = np.hstack((re_img[1], re_img[2])) #水平
    # combined_img = np.vstack((img, re_img))# 垂直
    combined_img=np.hstack((img,re_img))
    show_images(combined_img)
  • 0
    点赞
  • 2
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 0
    评论
好的,我可以在前面的代码基础上稍作修改来实现提取边界的功能。以下是代码示例: ```python import cv2 # 读取图像 img = cv2.imread('image/coins.png') # 灰度化处理 gray = cv2.cvtColor(img, cv2.COLOR_BGR2GRAY) # 边缘检测 edges = cv2.Canny(gray, 100, 200) # 膨胀操作 kernel = cv2.getStructuringElement(cv2.MORPH_RECT, (3, 3)) dilated = cv2.dilate(edges, kernel) # 腐蚀操作 eroded = cv2.erode(edges, kernel) # 形态学梯度 gradient = cv2.morphologyEx(edges, cv2.MORPH_GRADIENT, kernel) # 提取边界 inner_contours, _ = cv2.findContours(eroded, cv2.RETR_CCOMP, cv2.CHAIN_APPROX_SIMPLE) inner_boundary = cv2.drawContours(img, inner_contours, -1, (0, 0, 255), 2) # 提取边界 outer_contours, _ = cv2.findContours(dilated, cv2.RETR_CCOMP, cv2.CHAIN_APPROX_SIMPLE) outer_boundary = cv2.drawContours(img, outer_contours, -1, (0, 255, 0), 2) # 提取形态学梯度边界 gradient_contours, _ = cv2.findContours(gradient, cv2.RETR_CCOMP, cv2.CHAIN_APPROX_SIMPLE) gradient_boundary = cv2.drawContours(img, gradient_contours, -1, (255, 0, 0), 2) # 显示图像 cv2.imshow('Original Image', img) cv2.imshow('Inner Boundary', inner_boundary) cv2.imshow('Outer Boundary', outer_boundary) cv2.imshow('Gradient Boundary', gradient_boundary) cv2.waitKey(0) cv2.destroyAllWindows() ``` 运行以上代码,我们可以得到三个边界提取结果:内边界、外边界形态学梯度边界。其中,内边界可以通过先膨胀再腐蚀的操作得到;外边界可以通过先腐蚀再膨胀的操作得到;形态学梯度边界可以通过使用`cv2.drawContours()`函数和`cv2.findContours()`函数来实现。在`cv2.drawContours()`函数中,我们需要指定要绘制轮廓的图像,找到的轮廓,轮廓的索引(-1表示绘制所有轮廓),线条颜色和线条宽度等参数。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

Make_magic

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值