鱼眼相机模型-MEI

embedded/2024/11/26 14:19:19/

参考文献: Single View Point Omnidirectional Camera Calibration from Planar Grids

1. 相机模型如下:

   // 相机坐标系下的点投影到畸变图像// 输入:相机坐标系点坐标cam  输出: 畸变图像素点坐标disPtvoid FisheyeCamAdapter::cam2DistImg(cv::Point3f &cam, cv::Point2f &disPt){double r = cv::norm(cam);double dx_u = 0, dy_u = 0;if (r != 0){cv::Point3f ps = cam / r;double x = ps.x / (ps.z + camInt.fXi);double y = ps.y / (ps.z + camInt.fXi);distortion(x, y, &dx_u, &dy_u);x += dx_u;y += dy_u;disPt.x = x * camInt.fGammaX + camInt.fCx;disPt.y = y * camInt.fGammaY + camInt.fCy;}else{disPt.x = camInt.fCx;disPt.y = camInt.fCy;}}// 无畸变图像到畸变图像void FisheyeCamAdapter::distortion(double mx_u, double my_u, double *dx_u, double *dy_u){double mx2_u = 0., my2_u = 0., mxy_u = 0., rho2_u = 0., rad_dist_u = 0.;double k1 = camInt.distortCoeff[0];double k2 = camInt.distortCoeff[1];double p1 = camInt.distortCoeff[2];double p2 = camInt.distortCoeff[3];double k5 = camInt.distortCoeff[4];mx2_u = mx_u * mx_u;my2_u = my_u * my_u;mxy_u = mx_u * my_u;rho2_u = mx2_u + my2_u;rad_dist_u = k1 * rho2_u + k2 * rho2_u * rho2_u + k5 * rho2_u * rho2_u * rho2_u;*dx_u = mx_u * rad_dist_u + 2 * p1 * mxy_u + p2 * (rho2_u + 2 * mx2_u);*dy_u = my_u * rad_dist_u + 2 * p2 * mxy_u + p1 * (rho2_u + 2 * my2_u);}

   反投影过程:畸变图中的像素坐标计算相机坐标系下坐标

    cv::Point3f FisheyeCamAdapter::pointDis2Camera(const cv::Point2f &disPoint){double mx_d, my_d, mx_u, my_u;double lambda;double xi = camInt.fXi;// Lift points to normalised planefloat inv_K11 = 1 / camInt.fGammaX;float inv_K13 = -camInt.fCx / camInt.fGammaX;float inv_K22 = 1 / camInt.fGammaY;float inv_K23 = -camInt.fCy / camInt.fGammaY;mx_d = inv_K11 * (disPoint.x) + inv_K13;my_d = inv_K22 * (disPoint.y) + inv_K23;undistortGN(mx_d, my_d, &mx_u, &my_u, 100); // 去畸变坐标// Lift normalised points to the sphere (inv_hslash)cv::Point3f camera;if (xi == 1){lambda = 2 / (mx_u * mx_u + my_u * my_u + 1);camera.x = lambda * mx_u;camera.y = lambda * my_u;camera.z = lambda - 1;}else{double sqrt_i = 1.0 + (1.0 - xi * xi) * (mx_u * mx_u + my_u * my_u);if (sqrt_i < 0){camera.x = -1000000;camera.y = -1000000;camera.z = 1;}else{lambda = (xi + sqrt(sqrt_i)) / (1.0 + mx_u * mx_u + my_u * my_u);camera.x = lambda * mx_u;camera.y = lambda * my_u;camera.z = lambda - xi;}}return camera;}


http://www.ppmy.cn/embedded/140640.html

相关文章

设计模式之结构型模式

设计模式是软件开发中常见的解决方案&#xff0c;它们提供了一种在特定情况下解决常见问题的模板或框架。设计模式可以分为三大类&#xff1a;创建型模式、结构型模式和行为型模式。本文将重点介绍结构型模式&#xff08;Structural Design Patterns&#xff09;&#xff0c;并…

【6.10】位运算-实现四则运算

前言 从我们开始上学起就了解到&#xff0c;若要进行加减乘除就要使用对应的运算符号…… 甚至在如今的计算机中也是如此。我们通常只知道如何使用这些运算符号&#xff0c;却很少去关注其底层是如何实现的。假如某天在面试中遇到一道题&#xff0c;要求不使用 “&#xff0c;-…

昇腾CANN环境下Whisper.cpp安装指南

前置检查 确认昇腾AI处理器已经安装妥当 lspci | grep Processing accelerators ❕务必确认操作系统架构及版本、Python版本满足要求 软件 版本 操作系统 openEuler20.03/22.03, Ubuntu 20.04/22.04 Python 3.8, 3.9, 3.10 uname -m && cat /etc/*release py…

非标自动化项目管理如何做

非标自动化项目管理的关键在于&#xff1a;深入理解客户需求、制定详细的项目计划、有效的资源调配、严格的质量控制、持续的风险管理、高效的沟通协调、灵活应对变更、项目总结与持续改进。深入理解客户需求是项目成功的基础。通过与客户的深入沟通&#xff0c;全面了解其生产…

[面试]-golang基础面试题总结

文章目录 panic 和 recover**注意事项**使用 pprof、trace 和 race 进行性能调试。**Go Module**&#xff1a;Go中new和make的区别 Channel什么是 Channel 的方向性&#xff1f;如何对 Channel 进行方向限制&#xff1f;Channel 的缓冲区大小对于 Channel 和 Goroutine 的通信有…

网络安全设备系列--安全隔离网闸

0x00 定义: 网闸&#xff08;GAP&#xff09;全称安全隔离网闸。安全隔离网闸是一种由带有多种控制功能专用硬件在电路上切断网络之间的链路层连接&#xff0c;并能够在网络间进行安全适度的应用数据交换的网络安全设备。 网闸在两个不同安全域之间&#xff0c;通过协议转换的…

Websocket++ 框架

概述 基于 C 标准的 WebSocket 协议实现&#xff0c;适合使用 C11 或更高版本开发的项目&#xff0c;支持构建 WebSocket 客户端和服务器&#xff0c;功能灵活&#xff0c;设计模块化&#xff0c;非常适合高性能的网络通信场景 主要特点 支持 WebSocket 客户端和服务器默认支持…

lanqiaoOJ 3745:餐厅排队 ← 数组模拟队列

【题目来源】https://www.lanqiao.cn/problems/3745/learning/【题目描述】 在蓝桥学院的新餐厅&#xff0c;学生们在取餐窗囗形成了一条长队。小蓝&#xff0c;餐厅的经理&#xff0c;希望能够实时了解队伍最前面和最后面的学生编号。你需要执行以下三种操作&#xff1a; 1.学…