//矩阵的传入与传出
//20151021矩阵[1 3 50 100;3 5 50 100;5 7 50 100;7 9 50 100;9 11 50 100]出错，得进行调试
//opencv 的结果为[5 7 50 100]
//[2 4 50 100;4 6 50 100;6 6 50 100;8 10 50 100;10 12 50 100] 结果 6 8 50 100

#include "mex.h" 
#include <iostream>
#include <vector>
#include <opencv2/opencv.hpp>
using namespace cv;
using namespace std;
    void mexFunction(int nlhs, mxArray *plhs[], int nrhs, const mxArray *prhs[]){
        //输入，输入进来的矩阵一定是ｎ×４，不能为空
        double *inMatrix; 
        inMatrix = mxGetPr(prhs[0]); //得到输入矩阵的第一个元素的指针
        int mrows = mxGetM(prhs[0]);   //获得矩阵的行
        int ncols = mxGetN(prhs[0]);   //获得矩阵的列，

        vector<Rect> dst;
        dst.clear();
        Rect tempt;
        int dst_size;
        printf("%d_%d\n", mrows, ncols);  //打印行和列

        for(int i = 0; i < mrows; i++)
        {
            tempt.x=inMatrix [0 * mrows + i];//[当前列*行数+当前行]
            tempt.y=inMatrix [1* mrows + i];
            tempt.width=inMatrix [2* mrows + i];
            tempt.height=inMatrix [3 * mrows + i];
            dst.push_back(tempt);
        }
    //    printf("oriSize: %d\n",dst.size());  
        //cout<<"oriSize:"<<dst.size()<<endl;
        groupRectangles(dst, 1, 0.2);
        dst_size=dst.size();//当前行数
          printf("curSize: %d\n",dst.size());  
        //cout<<"curSize:"<<dst_size<<endl;

        //输出
        plhs[0] = mxCreateDoubleMatrix(dst_size,4, mxREAL);  //输出，每一列为一个窗口，行数=维数，列数=窗口数
        double *outMatrix;
        outMatrix = mxGetPr(plhs[0]);
        for(int i = 0; i < dst_size; i++)
        {
            outMatrix[0*dst_size+i ]=dst[i].x;
            outMatrix[1*dst_size+i ]=dst[i].y;
            outMatrix[2*dst_size+i ]=dst[i].width;
            outMatrix[3*dst_size+i ]=dst[i].height;
        }

}