compareHist.cpp

Go to the documentation of this file.

#include "mexopencv.hpp"
using namespace std;
using namespace cv;

namespace {
const ConstMap<string,int> HistComp = ConstMap<string,int>
    ("Correlation",     cv::HISTCMP_CORREL)
    ("ChiSquare",       cv::HISTCMP_CHISQR)
    ("Intersection",    cv::HISTCMP_INTERSECT)
    ("Bhattacharyya",   cv::HISTCMP_BHATTACHARYYA)
    ("Hellinger",       cv::HISTCMP_HELLINGER)
    ("AltChiSquare",    cv::HISTCMP_CHISQR_ALT)
    ("KullbackLeibler", cv::HISTCMP_KL_DIV);
}

void mexFunction(int nlhs, mxArray *plhs[], int nrhs, const mxArray *prhs[])
{
    // Check the number of arguments
    nargchk(nrhs>=2 && (nrhs%2)==0 && nlhs<=1);

    // Argument vector
    vector<MxArray> rhs(prhs, prhs+nrhs);

    // Option processing
    int method = cv::HISTCMP_CORREL;
    for (int i=2; i<nrhs; i+=2) {
        string key(rhs[i].toString());
        if (key=="Method")
            method = (rhs[i+1].isChar()) ?
                HistComp[rhs[i+1].toString()] : rhs[i+1].toInt();
        else
            mexErrMsgIdAndTxt("mexopencv:error","Unrecognized option");
    }

    // Process
    double d = 0;
    if (rhs[0].isSparse() && rhs[1].isSparse()) {
        SparseMat H1(rhs[0].toSparseMat(CV_32F)),
                  H2(rhs[1].toSparseMat(CV_32F));
        d = compareHist(H1, H2, method);
    }
    else {
        MatND H1(rhs[0].toMatND(CV_32F)),
              H2(rhs[1].toMatND(CV_32F));
        d = compareHist(H1, H2, method);
    }
    plhs[0] = MxArray(d);
}