使用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();
}