OpenCV图像形态学操作

import cv2
import numpy as np

# 腐蚀函数：cv2.erode(src, kernel, dst, anchor, iterations, borderType, borderValue)
# 参数说明：
# src: 输入图像
# kernel: 结构元素，若为Mat()则默认3×3矩形
# iterations: 迭代次数，次数越多腐蚀效果越明显

sun = cv2.imread('sun.png')
cv2.imshow('src', sun)
cv2.waitKey(0)

kernel = np.ones((4, 4), np.uint8)  # 定义4×4结构元素
erosion_1 = cv2.erode(sun, kernel, iterations=2)  # 迭代2次腐蚀
cv2.imshow('erosion_1', erosion_1)
cv2.waitKey(0)

src
erosion_1

2.图像膨胀

原理：

膨胀与腐蚀相反，会使前景物体边界向外扩张。核心逻辑是：用结构元素遍历图像，当结构元素与前景区域有任何重叠（即结构元素至少有一个像素落在前景中）时，就将结构元素中心位置的像素设为前景。

效果：可填补小缝隙、连接断开的物体、增大物体尺寸。

代码实现：

python 复制代码

# 膨胀函数：cv2.dilate(img, kernel, iteration)
# 参数说明：
# img: 目标图片
# kernel: 操作内核，默认3×3矩阵
# iterations: 膨胀次数，默认1

wenzi = cv2.imread('wenzi.png')
cv2.imshow('src1', wenzi)
cv2.waitKey(0)

kernel = np.ones((4, 4), np.uint8)  # 定义4×4结构元素
wenzi_new = cv2.dilate(wenzi, kernel, iterations=2)  # 迭代2次膨胀
cv2.imshow('wenzi_new', wenzi_new)
cv2.waitKey(0)

src1
wenzi_new

3.开运算与闭运算

开运算：先腐蚀后膨胀，用于平滑轮廓、断开窄颈、消除细突起
闭运算：先膨胀后腐蚀，用于弥合间断、消除孔洞、填补轮廓断裂

python 复制代码

# 开运算示例（指纹去噪）
zhiwen = cv2.imread('zhiwen.png')
cv2.imshow('src2', zhiwen)
cv2.waitKey(0)

kernel = np.ones((2, 2), np.uint8)  # 2×2结构元素
zhiwen_open = cv2.morphologyEx(zhiwen, cv2.MORPH_OPEN, kernel)  # 开运算
cv2.imshow('zhiwen_open', zhiwen_open)
cv2.waitKey(0)

# 闭运算示例（指纹补全）
zhiwen_duan = cv2.imread('zhiwen_duan.png')
cv2.imshow('src3', zhiwen_duan)
cv2.waitKey(0)

kernel = np.ones((4, 4), np.uint8)  # 4×4结构元素
zhiwen_close = cv2.morphologyEx(zhiwen_duan, cv2.MORPH_CLOSE, kernel)  # 闭运算
cv2.imshow('zhiwen_close', zhiwen_close)
cv2.waitKey(0)

src2
zhiwen_open
src3
zhiwen_close

4.梯度运算

梯度运算 = 膨胀结果 - 腐蚀结果，用于突出图像中强度变化剧烈的区域（提取边缘）。

python 复制代码

wenzi = cv2.imread('wenzi.png')
cv2.imshow('wenzi_original', wenzi)
cv2.waitKey(0)

kernel = np.ones((2, 2), np.uint8)  # 2×2结构元素

# 单独展示膨胀和腐蚀结果
pz_wenzi = cv2.dilate(wenzi, kernel, iterations=1)
cv2.imshow('dilated', pz_wenzi)
cv2.waitKey(0)

fs_wenzi = cv2.erode(wenzi, kernel, iterations=1)
cv2.imshow('eroded', fs_wenzi)
cv2.waitKey(0)

# 梯度运算提取边缘
bianyuan = cv2.morphologyEx(wenzi, cv2.MORPH_GRADIENT, kernel)
cv2.imshow('edge_gradient', bianyuan)
cv2.waitKey(0)

wenzi_original
dilated
eroded
edge_gradient

5.顶帽和黑帽

顶帽：原始图像 - 开运算结果，用于提取图像中比周围亮的细节
黑帽：闭运算结果 - 原始图像，用于提取图像中比周围暗的细节

python 复制代码

sun = cv2.imread('sun.png')
cv2.imshow('sun_original', sun)
cv2.waitKey(0)

kernel = np.ones((2, 2), np.uint8)  # 2×2结构元素

# 开运算结果
open_sun = cv2.morphologyEx(sun, cv2.MORPH_OPEN, kernel)
cv2.imshow('open_sun', open_sun)
cv2.waitKey(0)

# 顶帽运算
tophat = cv2.morphologyEx(sun, cv2.MORPH_TOPHAT, kernel)
cv2.imshow('TOPHAT', tophat)
cv2.waitKey(0)

# 闭运算结果
close_sun = cv2.morphologyEx(sun, cv2.MORPH_CLOSE, kernel)
cv2.imshow('close_sun', close_sun)
cv2.waitKey(0)

# 黑帽运算
blackhat = cv2.morphologyEx(sun, cv2.MORPH_BLACKHAT, kernel)
cv2.imshow('BLACKHAT', blackhat)
cv2.waitKey(0)

cv2.destroyAllWindows()

sun_original
open_sun
TOPHAT
close_sun
BLACKHAT

三、形态学操作的扩展应用与参数调整

1.结构元素的多样化选择

除了代码中使用的矩形结构元素，OpenCV 还支持自定义形状：

python 复制代码

# 矩形结构元素（默认）
rect_kernel = cv2.getStructuringElement(cv2.MORPH_RECT, (5, 5))
# 十字形结构元素（对线性特征敏感）
cross_kernel = cv2.getStructuringElement(cv2.MORPH_CROSS, (5, 5))
# 椭圆形结构元素（平滑处理曲线）
ellipse_kernel = cv2.getStructuringElement(cv2.MORPH_ELLIPSE, (5, 5))

print("矩形核:\n", rect_kernel)
print("十字核:\n", cross_kernel)
print("椭圆核:\n", ellipse_kernel)

python 复制代码

矩形核:
 [[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]]
十字核:
 [[0 0 1 0 0]
 [0 0 1 0 0]
 [1 1 1 1 1]
 [0 0 1 0 0]
 [0 0 1 0 0]]
椭圆核:
 [[0 0 1 0 0]
 [1 1 1 1 1]
 [1 1 1 1 1]
 [1 1 1 1 1]
 [0 0 1 0 0]]

2.形态学操作的链式应用

组合多种操作处理复杂图像，例如先去噪再提取边缘：

python 复制代码

def add_peppersalt_noise(image, n=10000):
    result = image.copy()
    h, w = image.shape[:2]  # 获取图片的高和宽
    for i in range(n):  # 生成n个椒盐噪声
        x = np.random.randint(low=1, high=h)
        y = np.random.randint(low=1, high=w)
        if np.random.randint(low=0, high=2) == 0:
            result[x, y] = 0
        else:
            result[x, y] = 255
    return result

noise_image = cv2.imread('../opencv/sun.png')
cv2.imwrite('noise.png', add_peppersalt_noise(noise_image))
# 读取带噪声的图像
noise_img = cv2.imread('noise.png', 0)  # 灰度图读取
cv2.imshow('带噪声原图', noise_img)
cv2.imwrite('noise_img.png', noise_img)
cv2.waitKey(0)

# 1. 开运算去噪
kernel = cv2.getStructuringElement(cv2.MORPH_ELLIPSE, (3, 3))
denoised = cv2.morphologyEx(noise_img, cv2.MORPH_OPEN, kernel)
cv2.imshow('去噪后', denoised)
cv2.imwrite('denoised.png', denoised)
cv2.waitKey(0)
 
# 2. 梯度运算提取边缘
edge = cv2.morphologyEx(denoised, cv2.MORPH_GRADIENT, kernel)
cv2.imshow('提取的边缘', edge)
cv2.imwrite('edge.png', edge)
cv2.waitKey(0)
cv2.destroyAllWindows()

带噪声原图
去噪后
提取的边缘

3.迭代次数对结果的影响

迭代次数直接影响操作强度：

python 复制代码

img = cv2.imread('../opencv/sun.png', 0)
kernel = np.ones((3, 3), np.uint8)

# 不同迭代次数的腐蚀效果
erosion_1 = cv2.erode(img, kernel, iterations=1)
erosion_7 = cv2.erode(img, kernel, iterations=7)
erosion_15 = cv2.erode(img, kernel, iterations=15)

cv2.imshow('迭代1次', erosion_1)
cv2.imwrite('erosion1.jpg', erosion_1)
cv2.imshow('迭代7次', erosion_7)
cv2.imwrite('erosion7.jpg', erosion_7)
cv2.imshow('迭代15次', erosion_15)
cv2.imwrite('erosion15.jpg', erosion_15)
cv2.waitKey(0)
cv2.destroyAllWindows()

迭代1次
迭代7次
迭代15次

现象：迭代次数越多，前景收缩越明显（腐蚀）或扩张越显著（膨胀）。

4.二值图像与灰度图像的处理差异

二值图像（0 和 255）：直接改变前景与背景分布，效果最明显
灰度图像：基于灰度值操作（腐蚀取最小值，膨胀取最大值）

python 复制代码

# 灰度图膨胀示例
gray_img = cv2.imread('gray_image.png', 0)
kernel = np.ones((3, 3), np.uint8)
gray_dilate = cv2.dilate(gray_img, kernel, iterations=1)

cv2.imshow('灰度原图', gray_img)
cv2.imshow('灰度膨胀后', gray_dilate)
cv2.waitKey(0)
cv2.destroyAllWindows()

灰度原图
灰度膨胀后

四、总结与实践建议

参数调试：从 3×3 结构元素和 1-5 次迭代开始尝试，避免过度处理
场景匹配 ：
- 去噪 → 开运算
- 补全轮廓 → 闭运算
- 提取边缘 → 梯度运算
- 提取亮 / 暗细节 → 顶帽 / 黑帽
技术结合 ：常与阈值分割（cv2.threshold()）配合，先二值化再处理可提升效果