【opencv】Mat 基本运算

mat 是为了与matlab矩阵运算对应而新加的opencv 类，有一些基本的操作，在下边的代码中列出，代码的资源见

http://download.youkuaiyun.com/detail/u010141025/8269695

//注意Mat的行列号是从0开始的
	//定义矩阵a,b,c
 	Mat a,b,c; 
	//生成三行四列的全一矩阵 CV_64F表示精度
 	a=Mat::ones(3,4,CV_64F);
	//a=mat::zeros(3,4,CV_64F);为生成全0
	//把矩阵a复制给矩阵b  注意不能用b=a
 	b=a.clone();
	//矩阵a每一个元素乘以2
	a=a.mul(2);
	//矩阵b每一个元素乘以4
 	b=b.mul(4);
	//矩阵a点乘矩阵b
 	c=a.mul(b);
   	cout<<"a"<<a<<endl;
   	cout<<"b"<<b<<endl;
   	cout<<"c"<<c<<endl;

运行结果：

       

       

 //常用的矩阵运算
 	Mat d,e,f,g;
	//带权重的矩阵加法  d=a*1+b*1+4,4为常数，1,1,4均可以改变值
 	addWeighted(a,1,b,1,4,d);
	//矩阵a除以a，结果为e
 	divide(d,a,e);
	//常数6除以矩阵e中的每一个元素，结果为f
 	divide(6,e,f);
	//常数2加矩阵f的每一个元素，结果为g，其中2可以为矩阵
 	add(2,f,g);
	cout<<"d"<<d<<endl;
	cout<<"e"<<e<<endl;
	cout<<"f"<<f<<endl;
 	cout<<"g"<<g<<endl;

运行结果：

 

//IplImage转为Mat
	IplImage *pImg=cvLoadImage("2.bmp",1);
	Mat mat(pImg);
	//若IplImage为3通道，则此时Mat也为三通道，可以将3通道转为单通道
	Mat matSingle;
	cvtColor(mat,matSingle,CV_RGB2GRAY);
	cout<<"mat"<<mat<<endl;
	cout<<"matSingle"<<matSingle<<endl;

运行结果：

     

//对矩阵的局部进行赋值改变
 	Mat onesMat=Mat::ones(3,6,CV_64F);
	//用常数初始化矩阵
	onesMat.setTo(4);
	//矩阵onesMat乘以2后赋值给matCount
 	Mat matCount=onesMat.mul(2); 
	//拷贝一整行元素
 	onesMat.row(1).copyTo(matCount.row(1));
 	//拷贝一整列元素
	onesMat.col(3).copyTo(matCount.col(3));
	//拷贝某一个元素
	onesMat.row(2).col(5).copyTo(matCount.row(2).col(5));
	cout<<matCount<<endl;

运行结果：
      

 

//访问矩阵中的元素
	for (int i=0;i<onesMat.cols;i++)
	{
		if (i==1||i==2)
		{
			for (int j=0;j<onesMat.rows;j++)
			{
				onesMat.row(j).col(i).setTo(-2);
			}
		}
	}
	cout<<"onesMat"<<onesMat<<endl;
	
	
	//用迭代器的方法访问矩阵中的每一个元素
	Mat matt(6,8,CV_64F);
	matt.setTo(4);
	double dTemp=0;
	MatIterator_<double> iterCount=matt.begin<double>();  //起始迭代器
	MatIterator_<double> iterEnd=matt.end<double>();      //终止迭代器
	for (;iterCount!=iterEnd;iterCount++)
	{
		dTemp=*iterCount; //取出元素
		dTemp=sqrt(dTemp);//对元素进行运算
		*iterCount=dTemp; //将运算后的结果返回
	}
	cout<<"matt"<<matt<<endl;

运行结果：

//将矩阵写入txt
	ofstream ofs("mat.txt");
	if (!ofs)
	{
		return 0;


	}
	int nCount=0;


	MatIterator_<double> iterCount1=matt.begin<double>();
	MatIterator_<double> iterEnd1=matt.end<double>();
	for (;iterCount1!=iterEnd1;iterCount1++)
	{
		ofs<<*iterCount1<<" ";
		nCount++;
		if (nCount==matt.cols)  //对矩阵进行换行
		{
			ofs<<endl;
			nCount=0;
		}
	}
	ofs.close();