opencv自定义线性滤波、边缘处理、Sobel算子(6)

本文详细介绍了使用OpenCV进行图像处理的各种技术,包括自定义线性滤波、边缘处理、Sobel算子和拉普拉斯算子的应用。通过具体代码示例,展示了如何实现图像的边缘检测和特征提取。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

1、自定义线性滤波

#include<iostream>
#include<opencv2/opencv.hpp>
using namespace std;
using namespace cv;
int main()
{
	Mat srcImg = imread("C:/Users/admin/Desktop/1.JPG");
	if (!srcImg.data)
	{
		printf("没有找到!!!");
		return -1;
	}
	Mat dstImg, dstImg2, dstImg3,dstImg4;
	Mat kernel_x = (Mat_<int>(3, 3) << -1, 0, 1, -2, 0, 2, -1, 0, 1);//sobel
	Mat kernel_y = (Mat_<int>(3, 3) << -1, -2, -1, 0, 0, 0, 1, 2, 1);
	Mat lap = (Mat_<int>(3, 3) << 0, -1, 0, -1, 4, -1, 0, -1, 0);//拉普拉斯算子

	int c = 0;
	int index = 0;
	int ksiz

e = 3;
while (true)
{
c = waitKey(500);
if ((char)c==27)//27键盘上esc
{
break;
}
ksize = 4 + (index %20) * 2 + 1;
Mat kernel = Mat::ones(Size(ksize, ksize), CV_32F) / (float)(ksize*ksize);
filter2D(srcImg, dstImg4, -1, kernel, Point(-1, -1));
index++;
imshow(“dstImg4 windows”, dstImg4);
}

filter2D(srcImg, dstImg3, -1, lap, Point(-1, -1));
namedWindow("dstImg3 windows", CV_WINDOW_AUTOSIZE);
imshow("dstImg3 windows",dstImg3);
filter2D(srcImg, dstImg, -1, kernel_x, Point(-1, -1));
namedWindow("dstImg windows", CV_WINDOW_AUTOSIZE);
filter2D(srcImg, dstImg2, -1, kernel_y, Point(-1, -1));
namedWindow("dstImg2 windows", CV_WINDOW_AUTOSIZE);
imshow("dstImg2 windows", dstImg2);
namedWindow("dstImg windows", CV_WINDOW_AUTOSIZE);
imshow("dstImg windows", dstImg);
namedWindow("srcImg windows", CV_WINDOW_AUTOSIZE);
imshow("srcImg windows", srcImg);
waitKey(0);
return 0;

}

2、边缘处理

#include<iostream>
#include<opencv2/opencv.hpp>
using namespace std;
using namespace cv;
int main()
{
	Mat srcImg = imread("C:/Users/admin/Desktop/1.JPG");
	if (!srcImg.data)
	{
		printf("没有找到!!!");
		return -1;
	}

	int top = (int)(0.05*srcImg.rows);
	int bottom = (int)(0.05*srcImg.rows);
	int left = (int)(0.05*srcImg.cols);
	int right = (int)(0.05*srcImg.cols);

	RNG rng(12345);
int c = 0;

int borderType = BORDER_DEFAULT;//默认处理方式
Mat dstImg;
while (true)
{
	c = waitKey(500);
	if ((char)c==27)
	{
		break;
	}
	if ((char)c=='r')
	{
		borderType = BORDER_REPLICATE;//用已知的边缘像素值填充
	}
	else if ((char)c=='w')
	{
		borderType = BORDER_WRAP;//用另一边的像素来补偿填充
	}
	else if ((char)c=='c')
	{
		borderType = BORDER_CONSTANT;//用指定像素来填充
	}
	Scalar color = Scalar(rng.uniform(0, 255), rng.uniform(0, 255), rng.uniform(0, 255));
	copyMakeBorder(srcImg, dstImg, top, bottom, left, right, borderType, color);
	namedWindow("dstImg windows", CV_WINDOW_AUTOSIZE);
	imshow("dstImg windows", dstImg);
}
/*namedWindow("dstImg windows", CV_WINDOW_AUTOSIZE);
imshow("dstImg windows", dstImg);*/
namedWindow("srcImg windows", CV_WINDOW_AUTOSIZE);
imshow("srcImg windows", srcImg);
waitKey(0);

return 0;

}

3、Sobel算子

#include<iostream>
#include<opencv2/opencv.hpp>
using namespace std;
using namespace cv;
int main()
{
	Mat srcImg = imread("C:/Users/admin/Desktop/1.JPG");
	if (!srcImg.data)
	{
		printf("没有找到!!!");
		return -1;  
	}

	Mat gray_img,dstImg;
	//进行高斯模糊
	GaussianBlur(srcImg, dstImg, Size(3, 3), 0, 0);
	//转为灰度图像
	cvtColor(dstImg, gray_img, CV_BGR2GRAY);
	imshow("gray_img window", gray_img);

	//sobel算子
Mat x_grad, y_grad;
Scharr(gray_img, x_grad, CV_16S, 1, 0);
Scharr(gray_img, y_grad, CV_16S, 0, 1);
//Sobel(gray_img, x_grad, CV_16S, 1, 0, 3); //CV_16S 深度
//Sobel(gray_img, y_grad, CV_16S, 0, 1, 3);
convertScaleAbs(x_grad, x_grad);
convertScaleAbs(y_grad, y_grad);

Mat xy_grad=Mat(x_grad.size(),x_grad.type());

//addWeighted(x_grad, 0.5, y_grad, 0.5, 0, xy_grad);//图像混合
int width = xy_grad.cols;
int height = xy_grad.rows;
//printf("%d", x_grad.type());
for (int row = 0; row < height; row++)
{
	for (int col = 0; col < width; col++)
	{
		int xg = x_grad.at<uchar>(row, col);
		int yg = y_grad.at<uchar>(row, col);
		
		xy_grad.at<uchar>(row, col) = saturate_cast<uchar>( xg + yg);
	}
}


imshow("xy_grad window", xy_grad);
imshow("x_grad window", x_grad);
imshow("y_grad window", y_grad);

waitKey(0);

return 0;

}

4、拉普拉斯算子

#include<iostream>
#include<opencv2/opencv.hpp>
using namespace std;
using namespace cv;

//拉普拉斯算子 原理:二阶微分,最大变化处的值为零,即边缘是零值。通过二阶导数计算;
/*
1、高斯模糊-去噪声GaussianBlur
2、转化为灰度图像 cvtColor()
3、拉普拉斯二阶导数计算
4、取绝对值 convertScaleAbs()
5、显示结果
*/
int main()
{

	Mat srcImg = imread("C:/Users/admin/Desktop/1.JPG");
	if (!srcImg.data)
	{
		printf("没有找到!!!");
		return -1;
	}
	//1、高斯模糊 - 去噪声GaussianBlur
	Mat gaussianBlur_image;
GaussianBlur(srcImg, gaussianBlur_image, Size(3, 3), 0, 0);
//2、转化为灰度图像 cvtColor()
Mat gray_img;
cvtColor(gaussianBlur_image, gray_img, CV_BGR2GRAY);
//3、拉普拉斯二阶导数计算
Mat laplance_img;
Laplacian(gray_img, laplance_img, CV_16S, 3);
//4、取绝对值 convertScaleAbs()
convertScaleAbs(laplance_img, laplance_img);
//5、显示结果
imshow("laplance window", laplance_img);

threshold(laplance_img, laplance_img, 0, 255, THRESH_OTSU | THRESH_BINARY);
imshow("laplance1 window", laplance_img);
namedWindow("srcImg windows", CV_WINDOW_AUTOSIZE);
imshow("srcImg windows", srcImg);
waitKey(0);
return 0;

}

### 线性滤波的概念 线性滤波是一种常见的图像处理技术,其中目标像素的新值是其本身及周围像素值的加权和。这一过程可以通过卷积操作实现,在此过程中使用的矩阵被称为核、掩模或模板。对于每一个位置(x, y),新的像素值由原图中对应位置及其邻域内像素值乘上相应权重后的总和决定[^2]。 ```python import cv2 import numpy as np def apply_linear_filter(image, kernel): """ 应用线性滤波器到给定图片 参数: image (numpy.ndarray): 输入图像数组. kernel (numpy.ndarray): 卷积核/滤波器. 返回: filtered_image (numpy.ndarray): 过滤后的图像数据. """ # 使用OpenCV库执行二维卷积运算 filtered_image = cv2.filter2D(src=image, ddepth=-1, kernel=kernel) return filtered_image ``` ### Sobel算子的工作机制 Sobel算子作为一种离散微分算子,主要用于边缘检测。该方法不仅考虑了局部差异还引入了一定量的平滑效果来减少噪声干扰。具体来说,Sobel算子分别沿X轴和Y轴方向计算一阶梯度幅值,并据此判断是否存在显著变化即可能存在的边界处。最终输出的结果反映了各点附近亮度变化的程度,从而帮助识别物体轮廓或其他特征结构[^4]。 ```python from matplotlib import pyplot as plt def sobel_operator(image): """ 对输入图像应用Sobel算子进行边缘检测 参数: image (numpy.ndarray): 输入灰度图像. 返回: edges_x (numpy.ndarray), edges_y (numpy.ndarray): X,Y 方向上的梯度幅度. """ # 定义Sobel算子的两个方向上的核 sobel_kernel_x = np.array([[-1, 0, 1], [-2, 0, 2], [-1, 0, 1]], dtype=np.float32) sobel_kernel_y = np.array([[1, 2, 1], [0, 0, 0], [-1, -2, -1]], dtype=np.float32) # 分别计算水平和垂直方向的一阶偏导数 edges_x = cv2.Sobel(image, cv2.CV_64F, 1, 0, ksize=3) edges_y = cv2.Sobel(image, cv2.CV_64F, 0, 1, ksize=3) return abs(edges_x).astype(np.uint8), abs(edges_y).astype(np.uint8) # 显示结果 fig, axs = plt.subplots(1, 3, figsize=(15, 5)) axs[0].imshow(cv2.cvtColor(original_img, cv2.COLOR_BGR2RGB)) axs[0].set_title('Original Image') edges_x, edges_y = sobel_operator(gray_scale_img) axs[1].imshow(edges_x, cmap='gray') axs[1].set_title('Edges detected along X axis') axs[2].imshow(edges_y, cmap='gray') axs[2].set_title('Edges detected along Y axis') plt.show() ``` ### Laplacian算子的作用方式 Laplacian算子是一个二阶微分算子,能够突出显示图像中的快速强度转换部分——也就是所谓的零交叉点(zero-crossings),这些地方往往标志着不同对象之间的交界或是表面属性的变化之处。当应用于图像时,Laplacian算子会强调那些具有较大曲率的位置;而在实际应用中,则经常用来做为一种简单的锐化工具或者是作为更复杂的多尺度分析框架的一部分[^3]。 ```python def laplacian_operator(image): """ 对输入图像应用拉普拉斯算子以增强细节或查找边缘 参数: image (numpy.ndarray): 输入灰度图像. 返回: sharpened_image (numpy.ndarray): 锐化的图像. """ # 创建一个自定义大小为3x3的拉普拉斯核 laplacian_kernel = np.array([ [0, 1, 0], [1,-4, 1], [0, 1, 0] ], dtype=float) # 执行卷积并返回结果 result = cv2.filter2D(image, -1, laplacian_kernel) return result sharpened = laplacian_operator(gray_scale_img) cv2.imshow("Sharpened", sharpened) cv2.waitKey(0); cv2.destroyAllWindows(); ```
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

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

余额充值