【OpenCV C++20 学习笔记】拉普拉斯(Laplace)二阶求导-边缘检测

拉普拉斯二阶求导

  • 原理
    • [拉普拉斯算子(Laplacian Operator)](#拉普拉斯算子(Laplacian Operator))
  • API
  • 实例

原理

在OpenCV中,Sobel算法可以对图片中的值求一阶导数,从而计算出图片中的边缘线。其原理如下面的示意图:

那么,如果再求一次导数的,即求二阶导数,其实也可以找出这个颜色值显著变化的分界点:

可以看到,现在颜色值显著变化的位置,其导数值为0.

但是这有一个问题,就是二阶导数为0的也可以是一些无意义的值。所以,必须要进行一些过滤。

拉普拉斯算子(Laplacian Operator)

拉普拉斯算子的算法公式定义如下:
L a p l a c e ( f ) = ∂ 2 f ∂ x 2 + ∂ 2 f ∂ y 2 Laplace(f) = \frac{\partial^2f}{\partial x^2} + \frac{\partial^2f}{\partial y^2} Laplace(f)=∂x2∂2f+∂y2∂2f

可以看到拉普拉斯算法可以同时对两个维度进行求导,这是它相对于Sobel算法的优势。但是由于拉普拉斯算法还是要求斜率,所以其内部仍然调用了Sobel算法。

API

在OpenCV中,使用Laplacian()函数来进行拉普拉斯计算,其函数原型如下:

cpp 复制代码
void cv::Laplacian(	InputArray	src,							//输入图
					OutputArray	dst,							//输出图
					int			ddepth,							//输出的数据类型,-1表示与输入图一致
					int			ksize = 1,						//卷积核尺寸,必须是正奇数
					double		scale =1,						//计算结果的放大系数,默认为1,即不放大
					double		delta = 0,						//计算结果的偏移值,默认为0,即不偏移
					int			borderType = BORDER_DEFAULT)	//图像边缘的扩充方式,默认为镜像复制
  • ksize = 1时,使用一个 3 × 3 3 \times 3 3×3的卷积核,如下:

    0 1 0 1 − 4 1 0 1 0 \] \\begin{bmatrix} 0 \& 1 \& 0 \\\\ 1 \& -4 \& 1 \\\\ 0 \& 1 \& 0 \\end{bmatrix} 0101−41010

实例

在进行拉普拉斯求导之前也要进行滤波和灰度化,以去除噪音。

这里我们将拉普拉斯计算的结果中的数据类型定义为CV_16S,是为了防止溢出。接着又通过convertScaleAbs()函数转换回了CV_8U类型。

完整代码如下:

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

using namespace cv;

int main() {
	Mat src{ imread("lena.jpg") };

	//高斯滤波
	Mat blured;
	GaussianBlur(src, blured, Size(3, 3), 0, 0, BORDER_DEFAULT);

	//灰度化
	Mat gray;
	cvtColor(blured, gray, COLOR_BGR2GRAY);

	//拉普拉斯
	Mat dst;
	Laplacian(gray, dst, CV_16S, 3, 1, 0, BORDER_DEFAULT);

	//转换为CV_8U
	Mat abs_dst;
	convertScaleAbs(dst, abs_dst);

	imshow("原图", src);
	imshow("Laplace", abs_dst);
	waitKey(0);
}

运行结果如下:

相关推荐
DKPT2 小时前
Java桥接模式实现方式与测试方法
java·笔记·学习·设计模式·桥接模式
好好研究5 小时前
学习栈和队列的插入和删除操作
数据结构·学习
彭祥.5 小时前
Jetson边缘计算主板:Ubuntu 环境配置 CUDA 与 cudNN 推理环境 + OpenCV 与 C++ 进行目标分类
c++·opencv·分类
新中地GIS开发老师6 小时前
新发布:26考研院校和专业大纲
学习·考研·arcgis·大学生·遥感·gis开发·地理信息科学
SH11HF6 小时前
小菜狗的云计算之旅,学习了解rsync+sersync实现数据实时同步(详细操作步骤)
学习·云计算
Tony沈哲6 小时前
macOS 上为 Compose Desktop 构建跨架构图像处理 dylib:OpenCV + libraw + libheif 实践指南
opencv·算法
Frank学习路上6 小时前
【IOS】XCode创建firstapp并运行(成为IOS开发者)
开发语言·学习·ios·cocoa·xcode
Chef_Chen8 小时前
从0开始学习计算机视觉--Day07--神经网络
神经网络·学习·计算机视觉
X_StarX10 小时前
【Unity笔记02】订阅事件-自动开门
笔记·学习·unity·游戏引擎·游戏开发·大学生
MingYue_SSS10 小时前
开关电源抄板学习
经验分享·笔记·嵌入式硬件·学习