opecv c++计算图像的曲率

公式

κ = z x x ⋅ z y 2 − 2 ⋅ z x ⋅ z y ⋅ z x y + z y y ⋅ z x 2 ( z x 2 + z y 2 + 1 ) 3 / 2 \kappa = \frac{z_{xx} \cdot z_y^2 - 2 \cdot z_x \cdot z_y \cdot z_{xy} + z_{yy} \cdot z_x^2}{(z_x^2 + z_y^2 + 1)^{3/2}}\newline κ=(zx2+zy2+1)3/2zxx⋅zy2−2⋅zx⋅zy⋅zxy+zyy⋅zx2

其中:
z x = ∂ z ∂ x z_x = \frac{\partial z}{\partial x}\newline zx=∂x∂z
z y = ∂ z ∂ y z_y = \frac{\partial z}{\partial y}\newline zy=∂y∂z
z x x , z y y , z x y z_{xx}, z_{yy}, z_{xy} zxx,zyy,zxy 分别为 ( z ) 关于 ( x ) 和 ( y ) 的二阶导数和交叉导数。

代码

cpp 复制代码
void computeCurvature(const cv::Mat& depth_map, cv::Mat& curvature_map)
{
	cv::Mat gradient_x, gradient_y;
	cv::Mat gradient_xx, gradient_yy, gradient_xy;

	// 计算一阶导数
	Sobel(depth_map, gradient_x, CV_64F, 1, 0, 3);
	Sobel(depth_map, gradient_y, CV_64F, 0, 1, 3);

	// 计算二阶导数
	Sobel(gradient_x, gradient_xx, CV_64F, 1, 0, 3);
	Sobel(gradient_y, gradient_yy, CV_64F, 0, 1, 3);
	Sobel(gradient_x, gradient_xy, CV_64F, 0, 1, 3);

	// 计算曲率
	cv::Mat denominator = gradient_x.mul(gradient_x) + gradient_y.mul(gradient_y) + 
		cv::Mat::ones(depth_map.size(), CV_64F);
	cv::pow(denominator, 1.5, denominator);

	curvature_map = (gradient_xx.mul(gradient_y.mul(gradient_y)) - 2.0 * gradient_x.mul(gradient_y.mul(gradient_xy)) + gradient_yy.mul(gradient_x.mul(gradient_x))) / denominator;
}
相关推荐
做人不要太理性16 分钟前
【C++】深入哈希表核心:从改造到封装,解锁 unordered_set 与 unordered_map 的终极奥义!
c++·哈希算法·散列表·unordered_map·unordered_set
程序员-King.25 分钟前
2、桥接模式
c++·桥接模式
chnming198729 分钟前
STL关联式容器之map
开发语言·c++
程序伍六七42 分钟前
day16
开发语言·c++
小陈phd1 小时前
Vscode LinuxC++环境配置
linux·c++·vscode
华清远见IT开放实验室1 小时前
【每天学点AI】实战图像增强技术在人工智能图像处理中的应用
图像处理·人工智能·python·opencv·计算机视觉
火山口车神丶1 小时前
某车企ASW面试笔试题
c++·matlab
只怕自己不够好1 小时前
《OpenCV 图像缩放、翻转与变换全攻略:从基础操作到高级应用实战》
人工智能·opencv·计算机视觉
是阿建吖!2 小时前
【优选算法】二分查找
c++·算法
Ajiang28247353044 小时前
对于C++中stack和queue的认识以及priority_queue的模拟实现
开发语言·c++