Function mexFunction

maxflow_matlab/automaxflowmex.cpp:16–119

[flow label]=automaxflowmex(I,fgProb, bgProb,lambda,sigma)

// [flow label]=automaxflowmex(I,fgProb, bgProb,lambda,sigma)
void mexFunction(int			nlhs, 		/* number of expected outputs */
                 mxArray		*plhs[],	/* mxArray output pointer array */
                 int			nrhs, 		/* number of inputs */
                 const mxArray	*prhs[]		/* mxArray input pointer array */)
{
    // input checks
    if (nrhs != 5 )
    {
        mexErrMsgTxt ("USAGE: [flow label]=automaxflowmex(I,fgProb, bgProb,lambda,sigma);");
    }
    const mxArray *I = prhs[0];
    const mxArray *fgProb = prhs[1];
    const mxArray *bgProb = prhs[2];
    double lamda= * mxGetPr(prhs[3]);
    double sigma= * mxGetPr(prhs[4]);
    
    unsigned char * IPr=(unsigned char *)mxGetPr(I);
    double * fgProbPr=mxGetPr(fgProb);
    double * bgProbPr=mxGetPr(bgProb);
    // size of image
    mwSize m = mxGetM(I);//height
    mwSize n = mxGetN(I);//width
    
    //construct graph
    typedef Graph<float,float,float> GraphType;
    GraphType *g = new GraphType(/*estimated # of nodes*/ m*n, /*estimated # of edges*/ 2*m*n);
    g->add_node(m*n);
    
    float maxWeight=-10000;
    for(int x=0;x<n;x++)
    {
        for(int y=0;y<m;y++)
        {
            //n-link
            float pValue=(float)*(IPr+x*m+y);
            int uperPointx=x;
            int uperPointy=y-1;
            int LeftPointx=x-1;
            int LeftPointy=y;
            float n_weight=0;
            if(uperPointy>=0 && uperPointy<m)
            {
                float qValue=(float)*(IPr+uperPointx*m+uperPointy);
                n_weight=lamda*exp(-(pValue-qValue)*(pValue-qValue)/(2*sigma*sigma));
                int pIndex=x*m+y;
                int qIndex=uperPointx*m+uperPointy;
                
                g->add_edge(qIndex,pIndex,n_weight,n_weight);
            }
            if(n_weight>maxWeight)
            {
                maxWeight=n_weight;
            }
            
            if(LeftPointx>=0 && LeftPointx<n)
            {
                float qValue=(float)*(IPr+LeftPointx*m+LeftPointy);
                n_weight=lamda*exp(-(pValue-qValue)*(pValue-qValue)/(2*sigma*sigma));