矩阵运算提速——玩转opencv::Mat

介绍:用Eigen或opencv::Mat进行矩阵的运算,比用cpp的vector或vector进行矩阵运算要快吗?

使用 Eigen 或 OpenCV 的 cv::Mat 进行矩阵运算通常比使用 std::vector<int>std::vector<double> 更快。这主要有以下几个原因:

  1. 优化的底层实现:Eigen 和 OpenCV 都是为高效的矩阵运算而设计的库,内部实现了许多优化,例如使用 SIMD(单指令多数据)指令、并行计算等。这些优化可以显著提高计算速度。

  2. 内存布局:Eigen 和 OpenCV 在内存布局上进行了优化,通常使用连续的内存块来存储矩阵数据,这样可以提高缓存命中率,从而加快访问速度。

  3. 专用的算法:这些库提供了许多专用的算法和函数,能够高效地执行常见的矩阵运算(如乘法、求逆、特征值分解等),而使用 std::vector 进行这些操作时,通常需要手动实现算法,效率较低。

  4. 类型安全和表达能力:Eigen 提供了强类型的矩阵和向量类,能够在编译时进行更多的检查,减少运行时错误,同时也使得代码更易于阅读和维护。

因此,如果你的应用涉及大量的矩阵运算,建议使用 Eigen 或 OpenCV 这样的专用库,而不是直接使用 std::vector

构造

link1

索引

访问Mat的通道数

#include <opencv2/opencv.hpp>
#include <iostream>int main() {// 创建一个3通道的彩色图像(例如,BGR格式)cv::Mat colorImage = cv::Mat::zeros(100, 100, CV_8UC3);// 创建一个单通道的灰度图像cv::Mat grayImage = cv::Mat::zeros(100, 100, CV_8UC1);// 获取通道数int colorChannels = colorImage.channels();int grayChannels = grayImage.channels();// 输出通道数std::cout << "彩色图像的通道数: " << colorChannels << std::endl; // 应该输出3std::cout << "灰度图像的通道数: " << grayChannels << std::endl;   // 应该输出1return 0;
}

定义4维Mat

#include <opencv2/opencv.hpp>
#include <iostream>int main() {// 定义4维矩阵,大小为2x3x4x5,数据类型为CV_32F(32位浮点数)cv::Mat mat4D(2, new int[4]{3, 4, 5}, CV_32F);// 填充矩阵for (int i = 0; i < 2; ++i) {for (int j = 0; j < 3; ++j) {for (int k = 0; k < 4; ++k) {for (int l = 0; l < 5; ++l) {mat4D.at<float>(i, j, k, l) = static_cast<float>(i * 1000 + j * 100 + k * 10 + l);}}}}// 输出矩阵的形状和内容std::cout << "4维矩阵的大小: " << mat4D.size << std::endl;std::cout << "4维矩阵的内容:" << std::endl;for (int i = 0; i < 2; ++i) {for (int j = 0; j < 3; ++j) {for (int k = 0; k < 4; ++k) {for (int l = 0; l < 5; ++l) {std::cout << mat4D.at<float>(i, j, k, l) << " ";}std::cout << std::endl;}std::cout << std::endl;}}// 释放动态分配的内存delete[] mat4D.size;return 0;
}

取某一行

#include <opencv2/opencv.hpp>
#include <iostream>int main() {// 创建一个 3x3 的矩阵cv::Mat mat = (cv::Mat_<float>(3, 3) << 1, 2, 3,4, 5, 6,7, 8, 9);// 获取第 3 行(索引为 2)的所有元素cv::Mat thirdRow = mat.row(2); // 行索引从 0 开始// 输出结果std::cout << "第三行的元素是:" << std::endl;std::cout << thirdRow << std::endl;return 0;
}

提取块

#include <opencv2/opencv.hpp>
#include <iostream>int main() {// 创建一个 5x5 的矩阵cv::Mat mat = (cv::Mat_<float>(5, 5) << 1, 2, 3, 4, 5,6, 7, 8, 9, 10,11, 12, 13, 14, 15,16, 17, 18, 19, 20,21, 22, 23, 24, 25);// 定义要提取的块的起始位置和大小int startRow = 1; // 起始行索引int startCol = 1; // 起始列索引int blockRows = 3; // 块的行数int blockCols = 3; // 块的列数// 提取块cv::Mat block = mat(cv::Range(startRow, startRow + blockRows), cv::Range(startCol, startCol + blockCols));// 输出结果std::cout << "提取的块是:" << std::endl;std::cout << block << std::endl;return 0;
}

访问某行某列的元素

在C++中,使用OpenCV库的cv::Mat类来表示图像或矩阵。要访问cv::Mat中的特定行和列的元素,可以使用at(row, col)方法,其中type是元素的数据类型。

#include <iostream>
#include <opencv2/opencv.hpp>int main() {// 创建一个3x3的矩阵,类型为CV_8UC1(单通道8位无符号整数)cv::Mat mat = (cv::Mat_<uchar>(3, 3) << 1, 2, 3,4, 5, 6,7, 8, 9);// 输出整个矩阵std::cout << "矩阵内容:\n" << mat << std::endl;// 访问特定行和列的元素int row = 1; // 第二行(索引从0开始)int col = 2; // 第三列(索引从0开始)// 使用at方法访问元素uchar value = mat.at<uchar>(row, col);std::cout << "元素在 (" << row << ", " << col << ") 的值: " << (int)value << std::endl;// 修改特定行和列的元素mat.at<uchar>(row, col) = 10;std::cout << "修改后的矩阵内容:\n" << mat << std::endl;return 0;
}

示例输出:

矩阵内容:
[1, 2, 3;
4, 5, 6;
7, 8, 9]
元素在 (1, 2) 的值: 6
修改后的矩阵内容:
[1, 2, 3;
4, 5, 10;
7, 8, 9]

运算

矩阵乘法和Hamornoid积

#include <iostream>
#include <opencv2/opencv.hpp>int main() {// 创建两个相同维度的矩阵cv::Mat mat1 = (cv::Mat_<float>(2, 2) << 1, 2,3, 4);cv::Mat mat2 = (cv::Mat_<float>(2, 2) << 5, 6,7, 8);// 输出原始矩阵std::cout << "矩阵1:\n" << mat1 << std::endl;std::cout << "矩阵2:\n" << mat2 << std::endl;// 方法1:使用 cv::multiplycv::Mat result1;cv::multiply(mat1, mat2, result1);std::cout << "矩阵乘法 (使用 cv::multiply):\n" << result1 << std::endl;// 方法2:使用 cv::mulcv::Mat result2 = mat1.mul(mat2); // 也可以使用 result2 = mat1 * mat2;std::cout << "矩阵乘法 (使用运算符 *):\n" << result2 << std::endl;// 方法3:使用运算符 *cv::Mat result3 = mat1*mat2; // 也可以使用 result2 = mat1 * mat2;std::cout << "逐元素相乘结果 (使用运算符 *):\n" << result3 << std::endl;return 0;
}

类型一致

  • opencv::Mat之间的加减乘除必须保证类型相同
  • 比如cv::divide()运算
cv::divide(new_cen, cen_index[cen_id].size(), new_cen);
  • opencv::Mat的类型有如下几种:
    在这里插入图片描述* 可在初始化时指定
	Mat new_cen=Mat::zeros(1,128,CV_32FC1);

opencv::Mat是继承Eigen的实现吗?

cv::Mat并不是继承自Eigen库的实现。OpenCV和Eigen是两个独立的库,它们各自有自己的数据结构和实现方式。

OpenCV的cv::Mat

  • cv::Mat是OpenCV中用于表示图像和矩阵的主要数据结构。它提供了多种功能,包括图像处理、矩阵运算、图像读取和写入等。
  • cv::Mat支持多维数组,能够处理不同的数据类型(如8位无符号整数、浮点数等),并且可以表示单通道和多通道图像(如RGB图像)。
  • OpenCV的设计目标是高效处理图像和视频数据,特别是在计算机视觉和图像处理领域。

Eigen库:

  • Eigen是一个用于线性代数的C++模板库,提供了矩阵和向量的高效操作。它专注于数学运算,特别是线性代数、几何变换和数值计算。
  • Eigen的设计目标是提供高性能的数学运算,同时保持易用性和灵活性。

结论:

虽然OpenCV和Eigen都涉及矩阵和线性代数的概念,但它们是独立的库,cv::Mat并不继承自Eigen的实现。实际上,OpenCV和Eigen可以一起使用,开发者可以在OpenCV中处理图像数据,然后使用Eigen进行更复杂的数学运算,或者反之亦然。两者之间的结合可以利用各自的优势。

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.rhkb.cn/news/500748.html

如若内容造成侵权/违法违规/事实不符,请联系长河编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

mac m2 安装 docker

文章目录 安装1.下载安装包2.在downloads中打开3.在启动台打开打开终端验证 修改国内镜像地址小结 安装 1.下载安装包 到官网下载适配的安装包&#xff1a;https://www.docker.com/products/docker-desktop/ 2.在downloads中打开 拖过去 3.在启动台打开 选择推荐设置 …

redis的集群模式与ELK基础

一、redis的集群模式 1.主从复制 &#xff08;1&#xff09;概述 主从模式&#xff1a;这是redis高可用的基础&#xff0c;哨兵和集群都是建立在此基础之上。 主从模式和数据库的主从模式是一样的&#xff0c;主负责写入&#xff0c;然后把写入的数据同步到从服务器&#xff…

建立一个Macos载入image的实例含界面

前言 为了方便ios程序的开发&#xff0c;有时候需要先用的Macos平台进行一些功能性的程序开发。 作为对比和参考。 1、创建一个MacOS的App 2、主界面控件的增加 添加的控件方法与ios相同&#xff0c;也是再用commandshiftL&#xff08;CtrlShiftL&#xff09;,就会弹出控件…

《机器学习》从入门到实战——逻辑回归

目录 一、简介 二、逻辑回归的原理 1、线性回归部分 2、逻辑函数&#xff08;Sigmoid函数&#xff09; 3、分类决策 4、转换为概率的形式使用似然函数求解 5、对数似然函数 ​编辑 6、转换为梯度下降任务 三、逻辑回归拓展知识 1、数据标准化 &#xff08;1&#xf…

实践:事件循环

实践&#xff1a;事件循环 代码示例 console.log(1); setTimeout(() > console.log(2), 0); Promise.resolve(3).then(res > console.log(res)); console.log(4);上述的代码的输出结果是什么 1和4肯定优先输出&#xff0c;因为他们会立即方式堆栈的执行上下文中执行&am…

【机器学习】工业 4.0 下机器学习如何驱动智能制造升级

我的个人主页 我的领域&#xff1a;人工智能篇&#xff0c;希望能帮助到大家&#xff01;&#xff01;&#xff01;&#x1f44d;点赞 收藏❤ 随着科技的飞速发展&#xff0c;工业 4.0 浪潮正席卷全球制造业&#xff0c;而机器学习作为这一变革中的关键技术&#xff0c;正以前…

自从学会Git,感觉打开了一扇新大门

“同事让我用 Git 提交代码&#xff0c;我居然直接把项目文件压缩发过去了……”相信很多初学者都经历过类似的窘境。而当你真正掌握 Git 时&#xff0c;才会发现它就像一本魔法书&#xff0c;轻松解决代码管理的种种难题。 为什么 Git 能成为程序员的标配工具&#xff1f;它究…

Mono里运行C#脚本21—mono_image_init_name_cache

前面分析了怎么样加载mscorlib.dll文件,然后把文件数据读取到内存。 接着下来,就会遇到加载整个C#的类型系统,比如System. Object,大体类型如下图所示: 在对CIL编译之前,需要把这些类型全部加载到内存里,以便快捷地访问它们。 mono_image_init_name_cache函数就是完成…

【Triton-ONNX】如何使用 ONNX 模型服务与 Triton 通信执行推理任务上-Triton快速开始

模型部署系列文章 前置-docker 理解:【 0 基础 Docker 极速入门】镜像、容器、常用命令总结前置-http/gRPC 的理解: 【HTTP和gRPC的区别】协议类型/传输效率 /性能等对比【保姆级教程附代码】Pytorch (.pth) 到 TensorRT (.plan) 模型转化全流程【保姆级教程附代码(二)】Pytor…

win32汇编环境,对话框中显示bmp图像文件

;运行效果 ;win32汇编环境&#xff0c;对话框中显示bmp图像文件 ;显示的是一张尺寸267*400的bmp位图,及一张缩小为原来三分之一的位图 ;将代码复制进radasm软件里&#xff0c;直接编译就可以运行了 ;下面为asm文件 ;>>>>>>>>>>>>>>&…

GAN对抗生成网络(二)——算法及Python实现

1 算法步骤 上一篇提到的GAN的最优化问题是&#xff0c;本文记录如何求解这一问题。 首先为了表示方便&#xff0c;记&#xff0c;这里让最大的可视作常量。 第一步&#xff0c;给定初始的&#xff0c;使用梯度上升找到 ,最大化。关于梯度下降&#xff0c;可以参考笔者另一篇…

多线程访问FFmpegFrameGrabber.start方法阻塞问题

一、背景 项目集成网络摄像头实现直播功能需要用到ffmpeg处理rtmp视频流进行web端播放 通过网上资源找到大神的springboot项目实现了rtmp视频流转为http请求进行视频中转功能&#xff0c;其底层利用javacv的FFmpegFrameGrabber进行拉流、推流&#xff0c;进而实现了视频中转。 …

X86、X64、64位、32位归纳总结

梳理一下位数区别和命名规范。 操作系统的位数决定了内存寻址空间的大小 X86_32的简称是X86。32位&#xff0c;最多只能识别4GB的内存。 X86_64的简称是X64。64位&#xff0c;最多能识别数十TB内存。 由于历史发展原因&#xff0c;形成了这样的简称&#xff0c;简称很关键。…

redis的学习(二)

4 哈希表 哈希类型中的映射关系通常称为field-value&#xff0c;⽤于区分Redis整体的键值对&#xff08;key-value&#xff09;&#xff0c; 注意这⾥的value是指field对应的值&#xff0c;不是键&#xff08;key&#xff09;对应的值&#xff0c; 4.1 操作命令 hset&#xff…

前端编码技巧与规范

当我们完成项目的构建&#xff0c;进入开发阶段的时候&#xff0c;除了你需要了解框架本身的知识点外&#xff0c;我们还需要提前掌握一些项目的编码技巧与规范&#xff0c;在根源上解决之后因编码缺陷而导致的项目维护困难、性能下降等常见问题&#xff0c;为项目多人开发提供…

谷歌广告关键词出价根据什么来判断?

投放广告的目的是为了盈利&#xff0c;而关键字的出价直接关系到广告费用的支出。因此&#xff0c;设定出价上限时&#xff0c;不仅要参考关键字规划师的建议&#xff0c;还需结合广告的盈利表现来合理判断。 可以在谷歌广告账户的后台查看“平均每次点击费用”和“每次点击的…

《我在技术交流群算命》(二):QGraphicsItem怎么写自定义信号啊(QObject多继承顺序问题)

某位群友突然无征兆的抛出以下问题&#xff1a; QGraphicsItem怎么写自定义信号啊 看到这个问题的时候我是比较疑惑的&#xff0c;按鄙人对 Qt 的了解&#xff0c;自定义信号只需: 继承QObject类中加入Q_OBJECT宏声明一个信号并使用 但该群友毕竟也不是一个Qt新手&#xff0…

filebeat采集应用程序日志和多行匹配

1 filebeat采集nginx json日志 01 修改nginx的日志为json格式 elk93节点安装nginx&#xff0c;注释掉默认的nginx日志格式&#xff1a;# access_log /var/log/nginx/access.log;&#xff0c;在下方增加以下配置。然后重启nginx log_format wzy_nginx_json {"timestamp&…

大语言模型提示技巧(二)-给模型时间思考

在与大语言模型交互的时候&#xff0c;如果模型给出了错误的结论&#xff0c;不要着急否定大模型的能力&#xff0c;我们应当尝试重新构建查询&#xff0c;请求模型在提供它的最终答案之前进行一系列相关的推理。也就是说&#xff0c;如果给模型一个在短时间或用少量文字无法完…

深入剖析MySQL数据库架构:核心组件、存储引擎与优化策略(一)

sql语句分为两大类&#xff1a;查询&#xff08;select&#xff09;、增删改----修改&#xff08;update&#xff09; select语句的执行流程 执行sql语句的流程&#xff1a;连接数据库、缓存查询、解析器、优化器、执行器、存储引擎操作数据 客户端&#xff1a;图形界面工具…