图像处理：Python使用OpenCV进行图像锐化 （非锐化掩模、拉普拉斯滤波器）

非锐化掩模 (Unsharp Masking)

步骤：

1. 模糊图像：使用高斯模糊滤波器对原图像进行模糊处理，得到模糊图像。
2. 计算细节层：通过从原图像中减去模糊图像，得到细节层。
3. 增强图像：将细节层乘以一个增益系数后加回到原图像，得到增强后的图像。

公式：
设原图像为 ( I )，模糊图像为 (I blur )，细节层为 ( D )，增益系数为 ( k )，最终的锐化图像 ( I’ ) 计算如下：

$$D=I-I_{\text{blur}}$$

$$I^{\prime }=I+k\cdot D$$

代码示例：

import cv2
import numpy as npdef unsharp_mask(image, k=1.5):# 高斯模糊图像blurred = cv2.GaussianBlur(image, (9, 9), 10.0)# 计算细节层detail = image - blurred# 增强图像sharpened = image + k * detailreturn np.clip(sharpened, 0, 255).astype(np.uint8)image = cv2.imread('Task3.jpg')
sharpened_image = unsharp_mask(image)
cv2.imshow('Original Image', image)
cv2.imshow('Unsharp Masked Image', sharpened_image)
cv2.waitKey(0)
cv2.destroyAllWindows()

拉普拉斯滤波器 (Laplacian Filter)

步骤：

1. 计算拉普拉斯图像：使用拉普拉斯算子计算图像的二阶导数，得到拉普拉斯图像。
2. 增强图像：将拉普拉斯图像加回到原图像中，得到锐化后的图像。

公式：
设原图像为 ( I )，拉普拉斯图像为 ( L )，最终的锐化图像 ( I’ ) 计算如下：

$$L=\Delta I=\frac{{\partial }^{2}I}{\partial x^2}+\frac{{\partial }^{2}I}{\partial y^2}$$

$$I^{\prime }=I+k\cdot L$$

代码示例：

import cv2
import numpy as npdef laplacian_sharpen(image, k=1.0):# 计算拉普拉斯图像laplacian = cv2.Laplacian(image, cv2.CV_64F)laplacian = np.uint8(np.absolute(laplacian))# 增强图像sharpened = cv2.addWeighted(image, 1, laplacian, k, 0)return sharpenedimage = cv2.imread('path_to_your_image.jpg')
sharpened_image = laplacian_sharpen(image)
cv2.imshow('Original Image', image)
cv2.imshow('Laplacian Sharpened Image', sharpened_image)
cv2.waitKey(0)
cv2.destroyAllWindows()

效果对比

了将非锐化掩模（Unsharp Masking）和拉普拉斯滤波器（Laplacian Filter）的方法整合到一个代码中，并对比展示效果，将两个锐化方法的结果放在同一个窗口中进行展示

import cv2
import numpy as np
import matplotlib.pyplot as pltdef unsharp_mask(image, k=1.5):# 高斯模糊图像blurred = cv2.GaussianBlur(image, (9, 9), 10.0)# 计算细节层detail = image - blurred# 增强图像sharpened = image + k * detailreturn np.clip(sharpened, 0, 255).astype(np.uint8)def laplacian_sharpen(image, k=1.0):# 计算拉普拉斯图像laplacian = cv2.Laplacian(image, cv2.CV_64F)laplacian = np.uint8(np.absolute(laplacian))# 增强图像sharpened = cv2.addWeighted(image, 1, laplacian, k, 0)return sharpeneddef display_images(original, unsharp, laplacian):titles = ['Original Image', 'Unsharp Masked ', 'Laplacian Sharpened ']images = [original, unsharp, laplacian]plt.figure(figsize=(10, 10)) for i in range(3):plt.subplot(1, 3, i + 1)plt.imshow(cv2.cvtColor(images[i], cv2.COLOR_BGR2RGB))plt.title(titles[i])plt.xticks([]), plt.yticks([])plt.show()def main():image_path = 'Task3.jpg'  # 请替换为你的图像路径image = cv2.imread(image_path)if image is None:print(f"Error: Unable to load image at {image_path}")returnunsharp_image = unsharp_mask(image)laplacian_image = laplacian_sharpen(image)display_images(image, unsharp_image, laplacian_image)if __name__ == "__main__":main()

具体效果对比如下：不同的图片的效果可能不同

总结

这两种锐化方法各有优缺点，要根据具体需求选择合适的方法：

• 非锐化掩模：

  • 优点：能够灵活控制图像的锐化程度，通过调整增益系数和模糊程度，可以获得较为自然的锐化效果。
  • 缺点：在处理带有高噪声的图像时，容易放大噪声。

• 拉普拉斯滤波器：

  • 优点：计算简单，能够快速增强图像边缘和细节。
  • 缺点：容易引入噪声和伪影，对噪声不敏感的图像效果更好。

通过应用这些方法，可以有效增强图像的边缘和细节，使图像看起来更加清晰和锐利。