opencv处理矩阵系列子程序

// cvMatTest.cpp : 定义控制台应用程序的入口点。
//


#include "stdafx.h"
#include <stdio.h>
#include "cv.h"




int _tmain(int argc, _TCHAR* argv[])
{






/*
//用固定数据创建矩阵
CvMat rotmat;
int sizes[2];
float array[] = {0,1,2,3,4,5,6,7,8,9,10,11,12,13,14};
//注意，此句更改为 cvInitMatHeader(&rotmat,3,5,CV_32FC1,NULL);时，cvGetDimSize不能起作用
cvInitMatHeader(&rotmat,3,5,CV_32FC1,array);
int type = cvGetElemType(&rotmat);
cvGetDims(&rotmat,sizes);
printf("type=%d\n",type);
printf("height=%d,width=%d\n",sizes[0],sizes[1]); 
printf("height=%d,width=%d\n",cvGetDimSize(&rotmat, 0),cvGetDimSize(&rotmat, 1));*/










/*
//矩阵数据存取简单方法------CV_MAT_ELEM_PTR
CvMat* mat = cvCreateMat(3,3,CV_32FC1);
cvZero(mat);//将矩阵置0
//为矩阵元素赋值
CV_MAT_ELEM( *mat, float, 0, 0 ) = 1.f;
CV_MAT_ELEM( *mat, float, 0, 1 ) = 2.f;
CV_MAT_ELEM( *mat, float, 0, 2 ) = 3.f;
CV_MAT_ELEM( *mat, float, 1, 0 ) = 4.f;
CV_MAT_ELEM( *mat, float, 1, 1 ) = 5.f;
CV_MAT_ELEM( *mat, float, 1, 2 ) = 6.f;
CV_MAT_ELEM( *mat, float, 2, 0 ) = 7.f;
CV_MAT_ELEM( *mat, float, 2, 1 ) = 8.f;
CV_MAT_ELEM( *mat, float, 2, 2 ) = 9.f;
//获得矩阵元素的值
float element1 = CV_MAT_ELEM(*mat,float,2,2);
float element2=*((float*)CV_MAT_ELEM_PTR(*mat,2,2));
printf("element1: %f\n",element1);
printf("element2: %f\n",element2);*/








/*
//矩阵数据存取困难方法------cvPtr2D
CvMat* mat = cvCreateMat(3,3,CV_32FC1);
cvZero(mat);//将矩阵置0
//为矩阵元素赋值
CV_MAT_ELEM( *mat, float, 0, 0 ) = 1.f;
CV_MAT_ELEM( *mat, float, 0, 1 ) = 2.f;
CV_MAT_ELEM( *mat, float, 0, 2 ) = 3.f;
CV_MAT_ELEM( *mat, float, 1, 0 ) = 4.f;
CV_MAT_ELEM( *mat, float, 1, 1 ) = 5.f;
CV_MAT_ELEM( *mat, float, 1, 2 ) = 6.f;
CV_MAT_ELEM( *mat, float, 2, 0 ) = 7.f;
CV_MAT_ELEM( *mat, float, 2, 1 ) = 8.f;
CV_MAT_ELEM( *mat, float, 2, 2 ) = 9.f;
//获得矩阵元素(0,2)的值
float *p = (float*)cvPtr2D(mat, 0, 2);
printf("矩阵元素(0,2)的值: %f\n",*p);*/


/*
//一维矩阵数据遍历--------cvPtr2D
CvMat* mat = cvCreateMat(3,3,CV_32FC1);
cvZero(mat);//将矩阵置0


//获得矩阵元素(0,0)的指针
float *p = (float*)cvPtr2D(mat, 0, 0);
//为矩阵赋值
for(int i = 0; i < 9; i++)
{
*p = (float)i;
p++;
}


//打印矩阵的值
p = (float*)cvPtr2D(mat, 0, 0);


for(int i = 0; i < 9; i++)
{
printf("%f\t",*p);
p++;
if((i+1) % 3 == 0)
printf("\n");
}*/




/*
//矩阵元素为三通道浮点数存取与遍历-------cvPtr2D
CvMat* mat = cvCreateMat(3,3,CV_32FC3);//注意此处用CV_32FC3表示RGB数值
cvZero(mat);//将矩阵置0
//为矩阵元素赋值


//获得矩阵元素(0,0)的指针
float *p = (float*)cvPtr2D(mat, 0, 0);
//为矩阵赋值
for(int i = 0; i < 9; i++)
{
//为每个通道赋值
*p = (float)i*10; 
p++;
*p = (float)i*10+1;
p++;
*p = (float)i*10+2;
p++;
}


//打印矩阵的值
p = (float*)cvPtr2D(mat, 0, 0);


for(int i = 0; i < 9; i++)
{
printf("%2.1f,%2.1f,%2.1f\t",*p,*(p+1),*(p+2));
p+=3;
if((i+1) % 3 == 0)
printf("\n");
}*/







/*
//采用cvGetReal2D
//矩阵元素为1通道浮点型数据
CvMat*mat=cvCreateMat(3,3,CV_32FC1 );
cvZero(mat);//将矩阵置0
//为矩阵元素赋值


//获得矩阵元素(0,0)的指针
float *p=(float*)cvPtr2D(mat,0,0);
//为矩阵赋值
for(int i=0;i<9;i++)
{
//为每个通道赋值
*p=(float)i*10; 
p++;
}


for(int i=0;i<3;i++)
for(int j=0;j<3;j++)
{
printf("%lf\t",cvGetReal2D(mat,i,j));
if ((i*3+j+1)%3==0)
printf("\n");
}*/








//采用cvSet*D给三通道矩阵复制3*4*3
//矩阵元素为三通道8位浮点数
/*
CvMat *mat=cvCreateMat(3,4,CV_32FC3 );
cvZero(mat);//将矩阵置0
//为矩阵元素赋值


for(int i = 0; i < 3; i++)
for(int j = 0; j < 4; j++)
cvSet2D( mat, i, j, cvScalar(i*10,j*10,i*j*10) );


for(int i=0;i<3;i++)
for(int j=0;j<4;j++)
{
printf("%lf,%lf,%lf\t",cvGet2D( mat, i, j ).val[0],cvGet2D( mat, i, j ).val[1],cvGet2D( mat, i, j ).val[2]);
if ((i*4+j+1)%3==0)
printf("\n");
}*/


//累加一个三通道矩阵的所有元素
float sum( const CvMat* newmat );
CvMat *mat=cvCreateMat(3,4,CV_32FC3 );
cvZero(mat);//将矩阵置0
//为矩阵元素赋值


for(int i = 0; i < 3; i++)
for(int j = 0; j < 4; j++)
cvSet2D( mat, i, j, cvScalar(i*10,j*10,i*j*10) );
for(int i=0;i<3;i++)
for(int j=0;j<4;j++)
{
printf("%lf,%lf,%lf\t",cvGet2D( mat, i, j ).val[0],cvGet2D( mat, i, j ).val[1],cvGet2D( mat, i, j ).val[2]);
if ((i*4+j+1)%3==0)
printf("\n");
}
printf("\n");
double s = sum(mat);
printf("Sum of the 3 channel is %f\n",s);


getchar();
return 0;
}


float sum(const CvMat* mat)
{
float s=0.0f;
for(int row=0;row<mat->rows;row++)
{
float* prt=(float*)(mat->data.ptr+row*mat->step);
for (int col=0;col<mat->cols;col++)
{
s+=prt[3*col]+prt[3*col+1]+prt[3*col+2];
}
}
return (s);
}