[OpenCV] 数字图像处理 C++ 学习——11自定义线性滤波 附完整代码

文章目录

前言

线性滤波用于去噪、平滑、边缘检测等操作。通过定义卷积核(算子)与图像进行卷积,线性滤波器可以有效地提取图像中的特定特征。本文将介绍几种常见的边缘检测算子,包括 Robert算子、Sobel 算子和拉普拉斯算子,并展示如何在 OpenCV 中实现这些滤波器的自定义操作。将提供完整的代码,以便学习使用。

1.理论基础

线性滤波是通过卷积操作实现的。在图像处理中,卷积是指使用一个小的滤波器(也称为卷积核或算子)在图像上移动,将核内的像素值进行加权平均,从而生成新的图像像素值。滤波器可以设计成不同的形状和大小,以实现不同的图像处理效果。

下面是卷积操作的动态过程,可以直观理解卷积:

(1)Roberts 算子 (Roberts Operator)

Roberts 算子是一种简单的边缘检测算子,通过计算图像的水平和垂直方向的差分来检测边缘。它使用两个 2x2 的卷积核来计算梯度的近似值,通常用于检测图像的斜边和锐利的边缘。

​ 水平方向(Roberts X) 垂直方向(Roberts Y):

(2)Sobel 算子 (Sobel Operator)

Sobel 算子是一种一阶导数滤波器,用于检测图像中的边缘。分别计算图像在水平方向和垂直方向上的梯度,通常用于边缘检测。相比 Roberts 算子,Sobel 算子更适合于检测大规模的边缘和平滑噪声。

(3)拉普拉斯算子 (Laplacian Operator)

拉普拉斯算子是一种二阶导数滤波器,用于检测图像中的边缘和细节。它通过计算图像的二阶导数来突出灰度变化较大的区域。

2.代码实现

图片下载链接实验图片链接lena.png

(1)Roberts 算子

Roberts 算子通过计算图像在水平方向和垂直方向的梯度来突出边缘特征,最终在两个窗口中显示使用 Roberts_X 和 Roberts_Y 核滤波后的结果。

cpp 复制代码
/****************************Roberts算子***********************************/
	// 定义Roberts算子核
	Mat Roberts_X = (Mat_<int>(2, 2) << 1, 0, 0, -1);
	Mat Roberts_Y = (Mat_<int>(2, 2) << 0, 1, -1, 0);
	//应用Roberts算子
	Mat Roberts_XFiltered, Roberts_YFiltered, RobertsGradient;
	filter2D(image, Roberts_XFiltered, -1, Roberts_X);
	filter2D(image, Roberts_YFiltered, -1, Roberts_Y);
	cv::imshow("Roberts_X", Roberts_XFiltered);
	cv::imshow("Roberts_Y", Roberts_YFiltered);

Roberts 算子结果:

(2)Sobel 算子

Sobel(image, sobelX, -1, 1, 0, 3);

  • ddepth=-1:输出图像的深度将与输入图像保持一致。
  • dx = 1:表示计算图像在 x 方向上的一阶导数,即水平梯度。
  • dy = 0:表示在 y 方向上不计算导数。
  • ksize = 3:使用 3x3 的 Sobel 核进行卷积操作。
cpp 复制代码
	/****************************Sobel 算子***********************************/
	// 定义Sobel X算子和Sobel Y算子
	Mat sobelX, sobelY;
	Sobel(image, sobelX, -1, 1, 0, 3); // 水平方向
	Sobel(image, sobelY, -1, 0, 1, 3); // 垂直方向
	// 显示结果
	cv::imshow("Sobel X Image", sobelX);
	cv::imshow("Sobel Y Image", sobelY);

Sobel 算子结果:

(3)拉普拉斯算子

一个 3x3 的拉普拉斯算子核,应用于图像以检测边缘。通过 filter2D 函数将该算子核与输入图像进行卷积,可以突出显示图像中的边缘和细节。

cpp 复制代码
/****************************拉普拉斯算子***********************************/
	// 定义3x3拉普拉斯算子
	Mat laplacianKernel = (Mat_<int>(3, 3) << 0, -1, 0, -1, 4, -1, 0, -1, 0);
	Mat laplacianFiltered;
	filter2D(image, laplacianFiltered, -1, laplacianKernel);
	cv::imshow("Laplacian Filtered Image", laplacianFiltered);

拉普拉斯算子结果:

3.完整代码

cpp 复制代码
#include<opencv2/opencv.hpp>
#include<highgui.hpp>

using namespace cv;
using namespace std;

void linear_filtering()
{
	cv::Mat image;
	image = imread("lena.png");
	if (image.empty()) {
		printf("could not find the image...\n");
		return;
	}
	namedWindow("input image", cv::WINDOW_AUTOSIZE);
	cv::imshow("input image", image);
	/**********************Roberts算子*******************************/
	// 定义Roberts算子核
	Mat Roberts_X = (Mat_<int>(2, 2) << 1, 0, 0, -1);
	Mat Roberts_Y = (Mat_<int>(2, 2) << 0, 1, -1, 0);
	//应用Roberts算子
	Mat Roberts_XFiltered, Roberts_YFiltered, RobertsGradient;
	filter2D(image, Roberts_XFiltered, -1, Roberts_X);
	filter2D(image, Roberts_YFiltered, -1, Roberts_Y);
	cv::imshow("Roberts_X", Roberts_XFiltered);
	cv::imshow("Roberts_Y", Roberts_YFiltered);

	/****************************Sobel算子***************************/
	// 定义Sobel X算子和Sobel Y算子
	Mat sobelX, sobelY;
	Sobel(image, sobelX, -1, 1, 0, 3); // 水平方向
	Sobel(image, sobelY, -1, 0, 1, 3); // 垂直方向
	// 显示结果
	cv::imshow("Sobel X Image", sobelX);
	cv::imshow("Sobel Y Image", sobelY);

	/************************拉普拉斯算子***************************/
	// 定义3x3拉普拉斯算子
	Mat laplacianKernel = (Mat_<int>(3, 3) << 0, -1, 0, -1, 4, -1, 0, -1, 0);
	Mat laplacianFiltered;
	filter2D(image, laplacianFiltered, -1, laplacianKernel);
	cv::imshow("Laplacian Filtered Image", laplacianFiltered);

	waitKey(0);
}
int main() 
{
	linear_filtering();
    return 0;
}
相关推荐
dreadp40 分钟前
解锁豆瓣高清海报(二) 使用 OpenCV 拼接和压缩
图像处理·python·opencv·计算机视觉·数据分析
涛ing1 小时前
32. C 语言 安全函数( _s 尾缀)
linux·c语言·c++·vscode·算法·安全·vim
独正己身2 小时前
代码随想录day4
数据结构·c++·算法
我不是代码教父4 小时前
[原创](Modern C++)现代C++的关键性概念: 流格式化
c++·字符串格式化·流格式化·cout格式化
利刃大大5 小时前
【回溯+剪枝】找出所有子集的异或总和再求和 && 全排列Ⅱ
c++·算法·深度优先·剪枝
子燕若水5 小时前
mac 手工安装OpenSSL 3.4.0
c++
*TQK*5 小时前
ZZNUOJ(C/C++)基础练习1041——1050(详解版)
c语言·c++·编程知识点
ElseWhereR6 小时前
C++ 写一个简单的加减法计算器
开发语言·c++·算法
*TQK*6 小时前
ZZNUOJ(C/C++)基础练习1031——1040(详解版)
c语言·c++·编程知识点
※DX3906※6 小时前
cpp实战项目—string类的模拟实现
开发语言·c++