使用VS2010调用matlab的mat格式文件

时间:2023-03-09 19:35:42
使用VS2010调用matlab的mat格式文件

做实验需要将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();
} dir = (const char **)matGetDir(pmat, &ndir); //ndir 表示mat文件中含有矩阵数目
if (dir == NULL) {
std::cout<<"Error reading directory of file"<<file<<endl;
return();
}
else
{
cout<<"Error opening file"<<file<<endl;
for (int i=; 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=;
for (int r=;r
{
int n=;
m++;
uchar *id=regId.ptr(r);
for (int c=;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();
}