opencv矩阵运算(2)

最新推荐文章于 2023-10-06 16:29:31 发布
最新推荐文章于 2023-10-06 16:29:31 发布
阅读量2.7k 收藏 4
点赞数
分类专栏： 图形图像 opencv 文章标签： opencv 图像处理
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/u011630458/article/details/46840705
版权

简介

　　本篇承接上一篇，继续opencv下矩阵计算的函数使用。

计算矩阵的逆

　　注意：矩阵A是可逆矩阵的充分必要条件是行列式detA不等于0。

具体代码

double x[3][3] = {{1, 2, 3}, {2, 2, 1}, {3, 4, 3}};
double y[3][3] = {{1, 0, 0}, {0, 2, 0}, {0, 0, 3}};
 
void showMatdate(Mat tmpMat){
	int i, j;
	CvScalar s1;
	Width = tmpMat.rows;
	Height = tmpMat.cols;
	IplImage tmp;
 
	tmp = tmpMat;
 
	for(i=0; i< Width; i++){
		for(j=0; j<Height; j++){
			s1 = cvGet2D(&tmp, i, j);
			printf("%0.1lf  ", s1.val[0]);
		}
		printf("\n");
	}
	printf("\n");
}
 
int main(int argc, char *argv[]){
 
/*************初始化矩阵*****************************/
	mat1 = Mat(3, 3, CV_64FC1, x);
	src1 = mat1;
	mat2 = Mat(3, 3, CV_64FC1, y);
	src2 = mat2;
/*************显示矩阵数据***************************/
	printf("mat1:\n");
	showMatdate(mat1);
 
/*****************矩阵的逆**********************/
	mat3 = Mat(3, 3, CV_64FC1);
	src3 = mat3;
	cvInvert(&src1, &src3);
	showMatdate(mat3);
	return 0;
}

结果显示

矩阵元素自然对数

具体代码

/*****************矩阵元素自然对数**********************/
    mat3 = Mat(3, 3, CV_64FC1);
    src3 = mat3;
    cvLog(&src1, &src3);
    showMatdate(mat3);

结果显示

矩阵查找表转换

具体代码

/*****************矩阵查找表转换**×******************/
    uchar lut[256];
 
    mat3 = Mat(3, 3, CV_8UC1);
    src3 = mat3;
    mat1.convertTo(mat1, CV_8UC1);
    src1 = mat1;
    mat2 = cvCreateMatHeader(1, 256, CV_8UC1);
    src2 = mat2;
    for (int i = 0; i < 256; i++) {
        lut[i] = 255 - i;
    }
    cvSetData(&src2, lut, 0);
    cvLUT(&src1, &src3, &src2);
    printf("cvLUT(mat1):\n");
    showMatdate(mat3);
    注意：mat3 = src2[mat1].（如果mat1格式为CV_8U）
          mat3 = src2[mat1 + 128].（如果mat1格式为CV_8S）

结果显示

计算向量间马氏距离

具体代码

/*************显示矩阵数据***************************/
    printf("mat1:\n");
    showMatdate(mat1);
    printf("mat2:\n");
    showMatdate(mat2);
 
/*****************计算向量间马氏距离**********************/                                                                          
    mat3 = Mat(3, 3, CV_64FC1, z);
    src3 = mat3;
    printf("mat3:\n");
    showMatdate(mat3);
 
    tmp = cvMahalanobis(&src1, &src2, &src3);
    printf("cvMahalanobis(mat1, mat2, mat3): %.1lf\n", tmp);
    马氏距离的定义，参考如下：http://blog.csdn.net/jmy5945hh/article/details/20536929

结果显示

获得矩阵元素间最大值

具体代码

/*****************计算矩阵参数间最大值**********************/
    mat3 = Mat(3, 3, CV_64FC1);
    src3 = mat3;
    cvMax(&src1, &src2, &src3);
    printf("cvMax(mat1, mat2):\n");
    showMatdate(mat3);
   类似的还有cvMaxS:计算矩阵元素和参数的最大值。
             cvAvg：计算矩阵元素的平均值。
             cvAvgSdv：计算矩阵元素的平均值和标准差。
             cvMin：计算矩阵参数间最小值。
             cvMinS：计算矩阵元素和参数的最小值。

结果显示

单通道合成多通道矩阵

具体代码

    mat3 = Mat(1, 3, CV_8UC3);                                                                                                       
    src3 = mat3;
    cvMerge(&src1, &src2, 0, 0, &src3);
    printf("cvMax(mat1, mat2):\n");
    showMatdate(mat3);
  对应的函数为：split()（将多通道分离为单通道矩阵）。

获取矩阵最大最小元素

具体代码

    printf("mat1:\n");
    showMatdate(mat1);
    cvMinMaxLoc(&src1, &min, &max, &min_p1, &max_p2);
    printf("min:%lf, min_p1.x:%d, min_p1.y:%d\n", min, min_p1.x, min_p1.y);
    printf("max:%lf, max_p2.x:%d, max_p2.y:%d\n", max, max_p2.x, max_p2.y);
     获取到最大最小值：max，min。以及它们的对应位置坐标：min_p1, max_p2。

结果显示

两个矩阵傅里叶频谱相乘

具体代码

    printf("mat1:\n");
    showMatdate(mat1);
    printf("mat2:\n");
    showMatdate(mat2);
 
    mat3 = Mat(3, 3, CV_64FC1);                                                                                                      
    src3 = mat3;
 
    cvMulSpectrums(&src1, &src2, &src3, DFT_ROWS);
    printf("mat3:\n");
    showMatdate(mat3);

结果显示

矩阵乘法

具体代码

    printf("mat1:\n");
    showMatdate(mat1);
    printf("mat2:\n");
    showMatdate(mat2);
 
    mat3 = Mat(3, 3, CV_64FC1);
    src3 = mat3;
 
    cvMul(&src1, &src2, &src3, 3);
    printf("cvMul(mat1 * mat2 * 3):\n");                                                                                             
    showMatdate(mat3);

结果显示

矩阵和转置的乘积

具体代码

    printf("mat1:\n");
    showMatdate(mat1);
    printf("mat2:\n");                                                                                                               
    showMatdate(mat2);
 
    mat3 = Mat(3, 3, CV_64FC1);
    src3 = mat3;
 
    cvMulTransposed(&src1, &src3, 0, &src2);
    printf("cvMulTransposed(mat1):\n");
    showMatdate(mat3);
    void cvMulTransposed( const CvArr* src, CvArr* dst, int order, const CvArr* delta=NULL );
    src:输入矩阵
    dst:目标矩阵
    order:乘法顺序
    delta:一个可选数组， 在乘法之前从 src 中减去该数组。函数 cvMulTransposed 计算 src 和它的转置的乘积。
函数求值公式：
        如果 order=0
              dst=(src-delta)*(src-delta)T
        否则
              dst=(src-delta)T*(src-delta)

结果显示

矩阵绝对差等

具体代码

    printf("mat1:\n");
    showMatdate(mat1);
    printf("mat2:\n");
    showMatdate(mat2);
 
    mat3 = Mat(3, 3, CV_64FC1);
    src3 = mat3;
 
    tmp = cvNorm(&src1, &src2, NORM_L1);                                                                                             
    printf("cvNorm(mat1, mat2, NORM_INF):%lf\n", tmp);
  double cvNorm(const CvArr* arr1, const CvArr* arr2=NULL, int norm_type=CV_L2, const CvArr* mask=NULL )
  如果arr2 == NULL
则：
    
否则：
      
   或者

结果显示

极性坐标转换到笛卡尔坐标

具体代码

    printf("mat1:\n");
    showMatdate(mat1);
    printf("mat2:\n");
    showMatdate(mat2);
 
    mat3 = Mat(3, 1, CV_64FC1);
    src3 = mat3;
    mat4 = Mat(3, 1, CV_64FC1);
    src4 = mat4;
 
    cvPolarToCart(&src1, &src2, &src3, &src4, true);
    printf("cvPolarToCart(mat1, mat2)--x:\n");
    showMatdate(mat3);
    printf("cvPolarToCart(mat1, mat2)--y:\n");
    showMatdate(mat4);
  void cvPolarToCart(const CvArr* magnitude, const CvArr* angle, CvArr* x, CvArr* y, int angle_in_degrees=0)

  magnitude:极坐标的长度。
  angle:极坐标的角度。
  x：笛卡尔X坐标。
  y：笛卡尔Y坐标。
  angle_in_degrees：若为true，表示输入的是角度，否则表示输入的是弧度。

结果显示

矩阵元素求幂

具体代码

    printf("mat1:\n");
    showMatdate(mat1);
 
    mat3 = Mat(3, 1, CV_64FC1);
    src3 = mat3;
 
    cvPow(&src1, &src3, 2);
    printf("cvPow(mat1, 2):\n");                                                                                                     
    showMatdate(mat3);

结果显示

矩阵简化为向量

具体代

    printf("mat1:\n");
    showMatdate(mat1);
 
    mat3 = Mat(3, 1, CV_64FC1);
    src3 = mat3;
 
    cvReduce(&src1, &src3, 1, CV_REDUCE_MAX);
    printf("cvReduce(mat1, 1 , CV_REDUCE_MAX):\n");
    showMatdate(mat3);
  void cvReduce(const CvArr* src, CvArr* dst, int dim=-1, int op=CV_REDUCE_SUM)
  src:待简化的矩阵。
  dst：生成的向量。
  dim：0意味着矩阵被处理成一行,1意味着矩阵被处理成为一列,-1时维数将根据输出向量的大小自动选择.
  op:
       CV_REDUCE_SUM-输出是矩阵的所有行/列的和.
       CV_REDUCE_AVG-输出是矩阵的所有行/列的平均向量.
       CV_REDUCE_MAX-输出是矩阵的所有行/列的最大值.
       CV_REDUCE_MIN-输出是矩阵的所有行/列的最小值.

结果显示
羽凌寒
关注
专栏目录
矩阵伪逆的opencv实现
qq_38231713的博客
12-06 1000
1、矩阵的逆 定义： 设A是数域上的一个n阶方阵，若在相同数域上存在另一个n阶矩阵B，使得： AB=BA=I。 则我们称B是A的逆矩阵，而A则被称为可逆矩阵。 可逆条件： A是可逆矩阵的充分必要条件是，即可逆矩阵就是非奇异矩阵。（当 时，A称为奇异矩阵） 性质： -矩阵A可逆的充要条件是A的行列式不等于0。 -可逆矩阵一定是方阵。 -如果矩阵A是可逆的，A的逆矩阵是唯一的。 -可逆矩阵也被称为非奇...
openCV入门之矩阵的运算
qq_47896251的博客
09-23 487
矩阵的运算 加法 1.mat加法 重载‘+’方法（mat+mat） 直接Mat m = mat1 + mat2; 注： 两边类型需相同 若值超出范围，则截断为最大值 例：uchar类型最大为255 ，结果123+145=268>255 则结果直接为255 重载‘+’方法（mat+value） 直接 Mat mat = mat1+value; value可以是任意类型的值，返回的mat的类型与mat1相同 add(mat1,mat2,dst, Mat(),CV_64FC1); 注： mat1,
OpenCV - Operations on Arrays 对数组(矩阵)的一些操作
weixin_34060299的博客
11-20 251
2019独角兽企业重金招聘Python工程师标准>>> ...
opencv中的矩阵操作
Rachel Zhang的专栏
03-23 2万+
有很多函数有mask，代表掩码，如果某位mask是0,那么对应的src的那一位就不计算，mask要和矩阵/ROI/的大小相等  大多数函数支持ROI，如果图像ROI被设置，那么只处理ROI部分  少部分函数支持COI，如果COI设置，只处理感兴趣的通道 矩阵逻辑运算  void cvAnd(const CvArr* src1, const CvArr* src2, CvArr* dst,
OpenCV矩阵函数操作大全
架构设计
07-10 205
转自：http://blog.csdn.net/freeboy1015/article/details/7256110 OpenCv矩阵操作 有很多函数有mask，代表掩码，如果某位mask是,那么对应的src的那一位就不计算，mask要和矩阵/ROI/的大小相等 大多数函数支持ROI，如果图像ROI被设置，那么只处理ROI部分 少部分函数支持COI，如果COI设置，...
OpenCV矩阵运算
merryken的专栏
03-06 906
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 26 27 28 29 30 31 32 33 34 35
opencv矩阵运算
08-18
下面是一些常见的OpenCV矩阵运算操作： 1. 创建矩阵： 使用`cv::Mat`类可以创建一个矩阵对象，并指定其大小和数据类型。例如，创建一个3x3的浮点型矩阵： ```cpp cv::Mat mat(3, 3, CV_32F); ``` 2. 访问像素...
opencv矩阵运算加速cuda
07-17
对于OpenCV中的矩阵运算加速，可以使用CUDA（Compute Unified Device Architecture）来利用GPU并行计算的能力。CUDA是NVIDIA提供的一种并行计算平台和编程模型，可以通过利用GPU的多个核心同时进行计算，从而加速...
opencv 平面法向量_OpenCV矩阵运算之顶点法向量计算
weixin_33931886的博客
02-08 847
OpenCV进行顶点法向量计算先介绍一下我接触此题的背景。我需要将一个.obj文件中的vn即顶点法向量计算出来。现在有的数据是f(三角面片索引)，v(顶点坐标)，但是网上找的代码不是很多，逻辑也不是很清楚，因此自己写了代码，并与大家分享。有错误敬请指出！顶点法向量的计算方法这里大概摘自这篇文章最简单也是最直接的求顶点的法向量的方法是遍历此顶点所关联的所有三角形，将这些三角形上的法向量单位化后叠加，...
OpenCV矩阵（向量）运算
SarKerson的博客
08-14 7021
Opencv 矩阵（向量）运算四则运算⎛⎝⎜⎜100010001⎞⎠⎟⎟+⎛⎝⎜⎜100010001⎞⎠⎟⎟=⎛⎝⎜⎜200020002⎞⎠⎟⎟ \left( \begin{array}{ccc} 1 & 0 & 0 \\ 0 & 1 & 0 \\ 0 & 0 & 1 \end{array} \right)+ \left( \begin{array}{ccc} 1 & 0 & 0 \\ 0
OpenCv 关于矩阵的相关计算函数
pekingFloater的博客
03-31 1995
GEMM通用矩阵乘法void cvGEMM( const CvArr* src1, const CvArr* src2, double alpha, const CvArr* src3, double beta, CvArr* dst, int tABC=0 ); #define cvMatMulAdd( src1, src2, src3, dst ) cvGEMM( src1, src2,...
opencv笔记 | 矩阵操作函数
【黑键】
05-21 178
矩阵操作函数
opencv求矩阵的逆_OpenCv矩阵操作函数大全
weixin_39873325的博客
12-19 486
有很多函数有mask，代表掩码，如果某位mask是0,那么对应的src的那一位就不计算，mask要和矩阵/ROI/的大小相等大多数函数支持ROI，如果图像ROI被设置，那么只处理ROI部分少部分函数支持COI，如果COI设置，只处理感兴趣的通道矩阵逻辑运算void cvAnd(const CvArr* src1,const CvArr* src2, CvArr* dst, const CvArr*...
OpenCV中对Mat矩阵加、减、乘、除、转置等操作的总结
热门推荐
进击的Explorer
11-08 3万+
1、矩阵乘 1） 矩阵乘A*B 是以数学运算中矩阵相乘的方式实现的，即Mat矩阵A和B被当做纯粹的矩阵做乘法运算，这就要求A的列数等于B的行数时，才能定义两个矩阵相乘。如A是m×n矩阵，B是n×p矩阵，它们的乘积AB是一个m×p矩阵。 例如： Mat A=Mat::ones(2,3,CV_32FC1); Mat B=Mat::ones(3,2,CV_32FC1); ... AB=A*B; 2） ...
opencv如何把一个矩阵不同列分离开_利于OPENCV实现分割算法的问题
weixin_39613951的博客
12-22 244
该楼层疑似违规已被系统折叠隐藏此楼查看此楼// meanshift_segmentation.cpp : 定义控制台应用程序的入口点。//#include "stdafx.h"#include #include #include #include #includeusing namespace cv;using namespace std;Mat src,dst;int spatialRad=1...
OpenCV----对图像进行幂次变换
裸奔人生
08-05 5725
uchar T[256]; for(int i = 0;i  {  T[i] = (uchar) (255 * pow(i,r) / pow(255.0f,r)); //幂次变换 } CvMat *T_mat = cvCreateMatHeader(1,256,CV_8UC1); //定义单通道，类型为8位无符号字符型的数组 cvSetData(T_mat,T,0);//指派用户数据给数组头 cv
OpenCV-除法运算cv::divide
翟天保的博客
08-08 1万+
介绍了OpenCV中除法运算函数cv::divide
【opencv】图像数字化——矩阵的运算（ 6 其他运算）
qq_40297054的博客
04-23 208
6 其他运算 6.1 指数和对数运算 针对矩阵中每一个数值进行相应的运算 exp函数 log函数：以e为底 以上两个函数的输入矩阵的数据类型只能是CV_32F或者CV_64F #include <opencv2/core/core.hpp> #include<iostream> using namespace std; using namespace cv; int main(){ Mat src1 = (Mat_<double>(2, 3) <
OpenCV学习笔记——图像非线性变换（幂律、对数）
最新发布
mimanchi99的博客
10-06 1161
幂律变换又称伽马变换，和对数变换的原理差不多，只是参数多了一个，可变宽带的输入像素值范围可选，但是把低值带拉伸还是高值带拉伸取决于伽马的设定。对数变换可以拉伸范围较窄的低灰度值，同时压缩范围较宽的高灰度值，可以用来扩展图像中的暗像素值，同时压缩亮像素值。·当r1,把输入的窄带暗值映射到宽带输出亮值，提高图像中亮区域的对比度。·幂律变换通过幂次变换曲线中的r值把输入的窄带宽值映射到宽带输出值。
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值