做实验需要将matlab实现的meanshift的结果中的region的Iabels矩阵,需要把labels.mat读入VS2010中,实现功能,在此把实现过程记录下来。
C++读取mat文件的步骤如下。
1. vs2010的配置
新建一个工程,在属性管理器中选择 工程->属性->VC++目录。
设置包含目录为:D:\MATLAB\R2010a\extern\include(我的matlab装在D盘);
设置库目录为:D:\MATLAB\R2010a\extern\lib\win32\microsoft。
2.设置环境变量path
计算机->属性->高级系统设置->高级->环境变量
在Path中已经存在D:\MATLAB_2011a\R2011a\bin\,但是D:\MATLAB_2011a\R2011a\bin\win32下有个重要的libmat.dll。如果不把d:\MATLAB\R2010a\bin\改为D:\MATLAB_2011a\R2011a\bin\win32程序将无法运行。
理论上D:\MATLAB_2011a\R2011a\bin\应该会向深一步检索到libmat.dll。坑爹的windows这一步并没有去做。这步设置完后要重启(环境配置完毕)。
3.使用MATLAB 提供的MAT 文件接口函数
在C/C++程序中有两种方式可以读取MAT文件数据。利用MATLAB提供的有关MAT文件的编程接口函数。MATLAB的库函数中包含了MAT文件接口函数库,其中有各种对MAT文件进行读写的函数,都是以mat开头的函数。如表1所示。
表1 C语言中的MAT文件读写函数
MAT 函 数 功 能
matOpen打开 MAT 文件
matClose关闭 MAT 文件
matGetDir从 MAT 文件中获得 MATLAB 阵列的列表
matGetFp获得一个指向MAT 文件的ANSI C 文件指针
matGetVariable从MAT文件中读取MATLAB阵列
matPutVariable写MATLAB阵列到MAT文件
matGetNextVariable从MAT文件中读取下一个MATLAB阵列
matDeleteVariable从MAT文件中删去下一个MATLAB阵列
matPutVariableAsGlobal从MATLAB阵列写入到MAT文件中
matGetVariableInfo从MAT文件中读取MATLAB阵列头信息
matGetNextVariableInfo从MAT文件中读取下一个MATLAB阵列头信息
4.程序举例
读取.mat文件,并把mat内容存储在图像中并显示。
注意:在Studio中一定要注意,默认读取.mat的方式是按列读取,不是按行读取。这是matlab数据存储的特色!
#include "stdafx.h" int loadMat(Mat& regId) { MATFile *pmat; const char **dir; const char *file; const char *name; int ndir; mxArray *pa; file="D:/labels3_ms.mat";//双反斜杠防止转义 pmat=matOpen(file, "r");//打开文件,返回指向文件指针 if (pmat == NULL) { cout<<"Error opening file"<<file<<endl; return(1); } dir = (const char **)matGetDir(pmat, &ndir); //ndir 表示mat文件中含有矩阵数目 if (dir == NULL) { std::cout<<"Error reading directory of file"<<file<<endl; return(1); } else { cout<<"Error opening file"<<file<<endl; for (int i=0; i < ndir; i++) cout<<"Mat NUM: "<<dir[i]<<endl;//输出所含矩阵数目 } pmat=matOpen(file, "r");//重新打开文件; pa = matGetNextVariable(pmat, &name); //返回指向文件头文件信息的指针,指针类型为*mxArray; //name 是矩阵的名字; cout<<name<<endl; pmat=matOpen(file, "r"); int ii=mxGetM(pa);//矩阵行数400 int jj=mxGetN(pa);//矩阵列数300 int *pm=(int*)mxGetPr(pa); //获取矩阵数值,返回指向矩阵第一数值的指针; regId=Mat::zeros(jj,ii,CV_8UC1); map<<span>int,int>regSize; int m=0; for (int r=0;r { int n=0; m++; uchar *id=regId.ptr(r); for (int c=0;c { id[c]=(*pm); pm++; regSize[*pm]++; } } cout<<"regSize:"<<regSize.size()<<endl; // for(map::iterator iter=regSize.begin();iter!=regSize.end();iter++) // { // cout<<iter->first<<" size: "<<iter->second<<endl; // } mxFree(dir); // //normalize(regId,regId,0,1,CV_MINMAX,CV_32FC1); regId=regId.t(); //imshow("show",regId); //waitKey(); //imwrite("D:/regId1.jpg",regId); //return 1; return regSize.size(); }