次要注意事项

Question

我想从matlab调用一个c函数，因为我尝试使用MEX编写一个包装函数。在编译时，我正在收集错误C2109：下标需要数组或指针类型和错误C2440：'function'：无法从'double *'转换为'double' 任何人都可以帮助我，我犯了错误？

#include "mex.h"
#include "matrix.h"
#include "CVIPtoolkit.h"
#include "CVIPtools.h"
#include "CVIPmatrix.h"
#include <stdio.h>
#include <stdlib.h>
#include <math.h>


void midd(double outdata, int type, int  height, int width){
     Image *outputImage;
     byte **output;
     int r,c;

     mexPrintf("type %d\n", type);
     mexPrintf("height %d\n", height);
     mexPrintf("width %d\n", width);

     outputImage=new_Image (PGM, GRAY_SCALE, 0, height,  width, CVIP_BYTE, REAL );
     outputImage = h_image(type, height,width);

     output = getData_Image(outputImage, 0);
     for(r=0; r < height; r++) {
            for(c=0; c < width; c++)
            {
                mexPrintf("type %d\n", type);
                mexPrintf("height %d\n", height);
                mexPrintf("width %d\n", width);

                outdata[r+height*c+height*width] =output[r][c];     /* passing data back to MATLAB variable from CVIPtools variable */
             }
         }     
 }


void mexFunction(int nlhs, mxArray *plhs[], int nrhs, const mxArray *prhs[])
{
    double *outdata;
    int type, height, width;

   // double *indata = (double *)mxGetData(prhs[0]);

   type = mxGetScalar(prhs[0]);
   height = mxGetScalar(prhs[1]);
   width = mxGetScalar(prhs[2]);

   mexPrintf("type %d\n", type);
   mexPrintf("height %d\n", height);
   mexPrintf("width %d\n", width);

   plhs[0] = mxCreateDoubleMatrix(height,width,mxREAL);

   outdata = mxGetData(plhs[0]);
   midd(outdata, type, height, width);

}

我试图调用的c函数如下：

#include "CVIPtoolkit.h"
#include <stdio.h>
#include <stdlib.h>
#include <math.h>


Image *
h_image(int type, unsigned int height, unsigned int width){

    /* type = 1, Constant
     * type = 2, Fixed mask
     * type = 3, Gaussian
     */

    unsigned int r, c, hf_w = width/2, hf_h = height/2;
    Image *outimage;
    float **outdata, sum = 0.0, sigma, tmp1, tmp2, tmp;

    if (height < 3 || width < 3) {
        fprintf(stderr, "Masksize too small, at least 3x3\n");
        return (Image *)NULL;
    }

    outimage = new_Image(PGM, GRAY_SCALE, 1, height, width, CVIP_FLOAT, REAL);
    outdata = (float **)getData_Image(outimage, 0);

    switch (type) {
        case 1:
            for (r = 0; r < height; r++) 
                for (c = 0; c < width; c++) {
                    outdata[r][c] = 1.0;
                    sum += outdata[r][c];
                }
            break;
        case 2:
            for (r = 0; r < height; r++) 
                for (c = 0; c < width; c++) {
                    outdata[r][c] = 1.0;
                    sum += outdata[r][c];
                }
            outdata[height/2][width/2] = height * width;
            sum = sum - 1.0 + outdata[height/2][width/2];
            break;
        case 3:
            c = (width + height) /4;
            r = (width + height) /2;
            sigma = sqrt(c*c / (2 * log(2) + (r - 3) * log(3)));
            sigma = 1.0 / 2.0 /sigma/sigma;
            tmp = width * height;
            for (r = 0; r < height; r++) 
                for (c = 0; c < width; c++) {
                    tmp1 = (r-hf_h)*(r-hf_h); tmp2 = (c-hf_w)*(c-hf_w);
                    outdata[r][c] = tmp*exp(- (tmp1 + tmp2) * sigma);
                    sum += outdata[r][c];
                }
            break;
        default:
            fprintf(stderr, "Incorrect mask type number: %d\n", type);
            return (Image *)NULL;
    }

    return outimage;
}

Answer 1

在您的主要功能中，outdata是指向double的指针，但您的功能midd会接收实际的double。这就是为什么你在类型中得到错误的原因。

只需更改函数声明，以便第一个输入接受指向double的指针：

void midd(double *outdata, int type, int  height, int width)
         //      ^^^^^^^^

次要注意事项

我质疑您将图像数据复制回MEX阵列：

outdata[r+height*c+height*width] =output[r][c];

您不需要height*width作为抵消。 r + height*c足以以列主要顺序访问单个通道2D矩阵。如果您有多通道图像，则只需要height*width偏移。该偏移量允许您访问其他通道中的图像数据...并且由于您只有单通道数据（看起来如此......），因此不需要此偏移。

因此，您只需要这样做：

outdata[r + height*c] = output[r][c];

如果你不这样做，我怀疑你最终会遇到分段错误，因为你最终会访问你不允许访问的部分内存。

此外，一旦完全测试了代码，就可以删除mexPrintf语句。由于您将命令提示符置于嵌套的for循环中，因此它会不必要地使用打印消息充斥您的命令提示符。我怀疑你这样做是为了调试，这很好，但我建议你将MEX函数附加到一个实际的调试器并正确调试代码而不是打印语句。

请参阅我的帖子，了解如何在此处设置该设置：Preventing a MEX file from crashing in MATLAB

从MATLAB调用c函数？

1 个答案:

次要注意事项