原文:http://blog.sina.com.cn/s/blog_4b0020f30101075w.html
CvMat矩阵数据结构是OpenCV的基础数据类型,对于图像处理这种密级型运算,经常需要访问,修改,设置其元素的值。OpenCV提供了很多优良的函数,能够很简单的实现上述功能。在《学习OpenCV》一书中,作者分别就简单的方法,麻烦的方法,和恰当的方法对相关函数进行了讲解,讲得比较清晰。本文主要讲解通过指针高效访问CvMat元素的方法z及平时容易忽视的一些小问题。
1.关于元素数据类型
CvMat中数据类型由几个部分构成 CV_<bit_depth>(S|U|F)channels,S表示有符号的,U表示无符号的,F表示浮点数;比如CV_32F1,表示32位1通道浮点数;CV_8U3,表示8位无符号3通道整形;数据类型重要的原因在于
A.它决定了CvMat数据的分布,比如,若元素类型CV_8UC1(常用于灰度图像 ),那么CvMat的数据排列是每行按照ggggggg(g表示一个像素的灰度值)的格式排列;若元素类型为CV_8UC3,则可以表示彩色图像,其行排列成为bgrbgrbgr(分别表示蓝绿红三个通道值,三个通道值表征1个像素)的形式;
B.在访问其数据类型时,如何正确转换成对应的数据类型指针;后面的例子会说到;
2.访问CvMat中的元素
简单的通过CV_MAT_ELEM宏,cvGetRealXD()函数即可实现,但是图像处理是计算密集型操作,这些函数虽然简单易用,但是效率比较低。因此最常用的是采用指针来访问CvMat中的元素。
CvMat结构中data结构对于指针访问其元素非常重要.
其data成员为
union{
uchar* ptr;
short* s;
int* i;
float* fl;
double* db;
}data;
由于是联合体,因此在访问时,指针可以在这几种类型的指针之间转换。当然这还归功于CvMat中每行的字节长度是固定的,成员step记录了CvMat每行的字节数。下面的代码说明如何高效访问矩阵元素。
CvMat* mat=cvCreateMat(5,3,CV_32FC1);//建立一个5行3列的矩阵,元素类型为32位单通道浮点数;数据类型对于后面使用指针访问矩阵元素非常重要。
cvZero(mat);
int row,col;
//下面的代码给矩阵的每个元素赋值
for (row=0;row<mat->height;row++)
{
float* pData=(float*)(mat->data.ptr+row*mat->step);//获取第row行的行首指针,因为数据类型为浮点型,因此,通过data.ptr与step获得的字节指针需要转换为float* 这样的指针
for (col=0;col<mat->width;col++)
{
*pData=(row+col);
pData++;//因为,指针后移一位,也即是指向下一个浮点数
}
}
//下面的3段代码功能是一样的,都是在控制台显示各个元素的值
///code1
for (row=0;row<mat->height;row++)
{
float* pData=(float*)(mat->data.ptr+row*mat->step);//我们通过转换成float*这样的指针
for (col=0;col<mat->width;col++)
{
std::cout<<*pData<<" ";
pData++;//指向下一个元素
}
}
///code2
for (row=0;row<mat->height;row++)
{
float* pData=mat->data.fl+row*mat->step/4;//这儿就用到了union中的fl指针,由于step是以1字节来计算的,因此实际步长应该为mat->step/4(浮点数的长度为4字节)
for (col=0;col<mat->width;col++)
{
std::cout<<*pData<<" ";
pData++;
}
std::cout<<std::endl;
}
///code 3
for (row=0;row<mat->height;row++)
{
float* pData=mat->data.fl+row*mat->step/4;
for (col=0;col<mat->width;col++)
{
std::cout<<*(pData+j)<<" ";//自己计算列指针的位置
}
std::cout<<std::endl;
}
///
用指针除了顺序访问CvMat中的元素外,还可以访问任意位置的元素,当然前提是需要自己计算指针。比如:
for (row=0;row<mat->height;row++)
{
float* pData=mat->data.fl+row*mat->step/4;
for (col=0;col<mat->width;col++)
{
if (row>0)
{
std::cout<< *(pData+col-mat->step/4)<<" ";//访问mat(row-1,col)元素
if (col>0)
std::cout<<*(pData+col-1-mat->step/4)<<"";//访问mat(row-1,col-1)元素
}
}
std::cout<<std::endl;
}
cvReleaseMat(&mat);//用完释放相关资源
因此,对CvMat中的元素不要拘泥于书上提供的几种方式,在程序效率很重要的情况下,可以合理使用指针结合step完美的访问CvMat中的元素,当然,使用指针,也有缺点,出现错误不容易发现,自己曾经遭过道。