用C/C++读取MATLAB中mat文件中的Cell数据

一、引言

最近由于研究需要，要用C语言调用MATLAB中的cell数据。于是，百度了一下，很容易找到了C++读写.mat文件。

但发现该文中说的是如何用C/C++调用mat文件中的矩阵，而我要调用的是mat文件中的cell数据，于是，又开始在网上狂找资料。不过这次就没那么幸运了，最后终于找到了一根”救命稻草“——mxGetCell()，之后在C++读写.mat文件的基础上自己摸索了一下，顺利解决了问题！！详见第二节。

二、C/C++调用MATLAB中的cell数据

1、mxGetCell原型

/*
 * Get a pointer to the specified cell element. 
 */ 
EXTERN_C mxArray *mxGetCell(const mxArray *pa, mwIndex i);

参数说明：

输入：

pa：指向mxArray类型的指针；

i:元胞中元素的标号；按列优先，依次为0,1,2…..。

输出：

指向mxArray类型的指针。

2、实例详解

a、这里要读入的是以下形式的元胞数组：

b、其中元胞a的主要内容如下：

c、C代码详解：

#include "mat.h"
int main()
{

    MATFile *pmat;//指向mat文件的指针
    const char **dir;//元素名列表
    const char *file;//要打开的mat文件名
    int         ndir;//mat文件中的元素(矩阵、元胞)个数
    mxArray *cell1;//指向要读取的元胞数据的指针
    mxArray *mat1;//指向元胞中某个元素的指针
    double *a;//指向元胞中某个元素第一个数据的指针
    int cellM,cellN,M,N;
    int count;
    int i,j;

    file="cellMat.mat";  //要打开的mat文件名
    pmat = matOpen(file, "r");//打开文件，返回指向文件指针
    if (pmat == NULL)
    {
        printf("打开mat文件失败！");
        return 0;
    }
    //读取mat文件中的矩阵列表(返回mat中包含的元素的名字)
    dir = (const char **)matGetDir(pmat, &ndir);
    if (dir == NULL)
    {
        printf("读取mat文件失败！");
        return 0;
    }   

    //输出cellMat.mat中的元素数量
    printf("%s文件中一共有%d个元素\n",file,ndir);

    //从打开的mat文件中(pamt)读取名字为dir[0]的元胞
    //(返回指向该矩阵的指针，注意：MATLAB中的矩阵是按列优先存储到)
    cell1 = matGetVariable(pmat,dir[0]);//指向名字为dir[0]的元胞

    cellM = (int)mxGetM(cell1);//获取元胞的行数
    cellN = (int)mxGetN(cell1);//获取元胞的列数
    //按列优先依次读取元胞中的数据
    for(count=0;count<cellM*cellM;count++)
    {
        mat1 = mxGetCell(cell1,count);//按列优先，指向元胞中的第count+1个元素

        a = (double*)mxGetData(mat1);//指向mat1中的第一个数据

        M = (int)mxGetM(mat1);//获取mat1的行数
        N = (int)mxGetN(mat1);//获取mat1的列数

        printf("元素%d的维数为:(%d,%d)\n",count,M,N);

        for (i=0;i<M;i++)
        {
            for (j=0;j<N;j++)
                printf("%-5.4f ",a[j*M+i]);
            printf("\n");
        }
    }
    matClose(pmat);//关闭文件
    return 0;
}

d、代码运行结果：

从运行结果来看和上述MATLAB中显示的数据是一致的。

三、进一步探讨

由于上述读取的mat文件中只包含一个元胞数组，过于简单，大家可以参考上述代码，尝试读一下包含多个cell的mat文件以及多个cell与矩阵混合的mat文件，以此来加深理解。

最后，再次感谢C++读写.mat文件的作者！

说明：

对于上述代码中表述不清或不准确的地方欢迎大家批评指正！！