图像处理|闭运算

闭运算（Closing）是形态学操作中的另一种基本操作，它与开运算（Opening）类似，但执行的步骤相反。闭运算结合了膨胀和腐蚀操作，顺序为 先膨胀后腐蚀。这种操作通常用于填补图像中的小空洞，同时去除细小的黑色噪声。

闭运算具体过程：

膨胀操作：将图像中的白色区域扩展，填补小的黑色区域或间隙。
腐蚀操作：膨胀之后，对图像进行腐蚀，减少过多的边界，保持物体的形状。

闭运算的效果是填补图像中的小空洞、消除小黑点噪声等。

闭运算的作用

填补小的黑色空洞：闭运算可以填补物体内部的小空洞和裂缝，使物体更加连贯。
去除小的白色噪声：可以去除图像中的一些小白噪声点。
平滑图像边缘：闭运算会略微平滑物体的边缘，同时保持物体的整体形态。

在 OpenCV 中的实现

OpenCV提供了cv2.morphologyEx函数来实现形态学操作，包括闭运算。我们可以通过参数 cv2.MORPH_CLOSE 来进行闭运算。

主要函数

cv2.morphologyEx(src, op, kernel)

src：输入图像（通常是二值图像）。
op：形态学操作的类型，例如 cv2.MORPH_CLOSE 代表闭运算。
kernel：结构元素，定义了操作的范围和形状。

代码示例

以下代码演示了如何使用 OpenCV 实现闭运算去除白色方块中的黑色噪点。

import cv2
import numpy as np

image = cv2.imread('./noise.png')
image = cv2.cvtColor(image, cv2.COLOR_BGR2GRAY)
_, threshed = cv2.threshold(image, 127, 255, cv2.THRESH_BINARY)
print(threshed.shape)
k = cv2.getStructuringElement(cv2.MORPH_RECT, (5, 5))
closed = cv2.morphologyEx(threshed, cv2.MORPH_CLOSE, k)
stack_img = np.hstack((image, closed))
cv2.imshow('img', stack_img)
cv2.waitKey(0)
cv2.destroyAllWindows()

解释

原图像：包含小的黑色空洞或噪声点。
闭运算结果：通过先膨胀再腐蚀的方式，填补了物体内部的小空洞，并去除了小的噪声点。

结构元素的选择

在闭运算中，结构元素（kernel）的形状和大小会影响最终的结果。常见的结构元素有：

矩形结构元素：cv2.MORPH_RECT
椭圆形结构元素：cv2.MORPH_ELLIPSE
交叉形结构元素：cv2.MORPH_CROSS

可以根据实际需要选择不同形状的结构元素：

kernel = cv2.getStructuringElement(cv2.MORPH_ELLIPSE, (5, 5))

应用场景

填补物体中的小空洞：闭运算常用于图像中物体的细节处理，尤其是填补物体中间的小空洞或裂缝，使物体看起来更加连贯。
去除小型噪声点：通过先膨胀后腐蚀，能够有效去除小的白色噪声点（在黑色背景上的小白点）。
图像的边界平滑：在图像的边缘进行平滑处理，同时保持整体形状。
增强图像特征：通过闭运算增强图像中的结构，尤其是在物体识别或图像分割前的预处理阶段。

常见问题与优化

选择合适的结构元素大小：
- 结构元素的大小决定了闭运算操作的强度。若选择太小的结构元素，可能无法有效填补空洞；若选择过大的结构元素，可能会导致边界的过度扩展。
```
# 创建一个较大的结构元素，适合更大区域的填补 
kernel = cv2.getStructuringElement(cv2.MORPH_RECT, (15, 15))
```

处理彩色图像：

对彩色图像进行闭运算时，可以对每个通道（例如RGB三个通道）分别进行操作，或者将图像转换为灰度图后进行处理。

import cv2
import numpy as np

image = cv2.imread("frame.jpg")
kernel = cv2.getStructuringElement(cv2.MORPH_RECT, (5, 5))
channels = cv2.split(image)
processed_channels = [cv2.morphologyEx(ch, cv2.MORPH_CLOSE, kernel) for ch in channels]
result = cv2.merge(processed_channels)

stack = np.hstack((image, result))
cv2.imshow('img', stack)
cv2.waitKey(0)
cv2.destroyAllWindows()