Question

I have already passed an image to my mexFunction but now I need to pass an array of images and I am struggling to get the thing right. This is my code to get the simple Image. This works perfectly but when I go into 3D I don t understand how the information is arranged in the mxArray.

void mexFunction(int nlhs, mxArray *plhs[], int nrhs, mxArray *prhs[]) 
{
    mxArray *matrixIn = prhs[0];
    inputImage=(double *)mxGetPr(matrixIn);
    int x = int(dims[0]);
    int y = int(dims[1]);
    volume3D image(inputImage, x, y, 1);
}

volume3D::volume3D(double* image, int x, int y, int z) 
{
    allocateVolume( x, y, z); 
    for(int i=0; i<xSize; i++)
        for(int j=0; j<ySize; j++) {
            volume[i][j][0] = double(image[(i)*x+j]);
        }
}

我这样做是为了绕过它。

mwSize mrows,ncols; mrows = mxGetM(prhs[0]); ncols = mxGetN(prhs[0]);

plhs[0] = mxCreateNumericMatrix(mrows, ncols, mxDOUBLE_CLASS, mxREAL);
double *matlabTumorMap = mxGetPr(plhs[0]);
const int * dims = mxGetDimensions( plhs[0]);
int x = int(dims[0]);
int y = int(dims[1]);
int z = int(dims[2]);



mwIndex subs[3];
mexPrintf("x %i
",x);
mexPrintf("y %i
",y);
mexPrintf("z %i
",z);
mxArray *matrixTumor = plhs[0];

for(subs[0]=0; subs[0]<x; subs[0]++)
    for(subs[1]=0; subs[1]<y; subs[1]++)
        for(subs[2]=0; subs[2]<z; subs[2]++)
        {
            mwIndex x = mxCalcSingleSubscript( matrixTumor,3,subs);
            matlabTumorMap[x] = tumorMap.getVoxel(subs[0],subs[1],subs[2]);
        }

Answer 1

http://www.mathwork.de/help/techdoc/apiref/bqoqnz0.html” rel=“nofollow” http://www.mathwork.de/help/techdoc/apiref/bqoqnz0.html, http://www.mathwork.de/help/techdoc/apiref/mxcalcsinglesubscript.html” rel=“nofollow”mxCalcSingle Subscript 。这有助于你计算这些数据。

Something like

    mxArray *matrixIn = prhs[0];
    volume3D image(matrixIn);
}

volume3D::volume3D(MxArray* matrixIn) 
{
    double * inputImage=(double *)mxGetPr(matrixIn);
    assert(mxGetNumberOfDimensions(matrixIn) >= 3)
    mwSize * dims = mxGetDimensions(matrixIn);
    int x = int(dims[0]);
    int y = int(dims[1]);
    int z = int(dims[2]);

    double * image = mxGetPr(matrixIn);

    mwIndex subs[3];

    allocateVolume( x, y, z); 
    for(subs[0]=0; subs[0]<x; subs[0]++)
        for(subs[1]=0; subs[1]<y; subs[1]++)
            for(subs[2]=0; subs[2]<z; subs[2]++) {
                mwIndex x = mxCalcSingleSubscript(matrixIn, 3, subs);
                /* <unsure> */volume[subs[0]][subs[1]][subs[2]]  /* </unsure> */ = image[x];
            }

BTW:如果把C和C++混为一谈,则会由于名字命名而导致更多的头痛。

Answer 2

你们就是这样做的。

我认为,唯一的问题是你的指数化。请写:

 volume[i][j][0] = double(image[i+j*x]);

以及你打算写:

mwSize* dims = mxGetDimensions(matrixIn);

友情链接