与cvFindContours配合使用cvCalcPGH

Question

出于某种原因，谷歌没有使用cvCalcPGH()的好例子，也没有使用Freeman代码进行轮廓匹配的好例子。所以我试着做到：

#include "opencv2/opencv.hpp"
#include "opencv2/legacy/legacy.hpp"

int
main(int argc, char* argv[]) {

        IplImage* g_gray = cvLoadImage("lisa.png",CV_LOAD_IMAGE_GRAYSCALE);
        cvThreshold( g_gray, g_gray, 100, 255, CV_THRESH_BINARY );

        CvMemStorage* storage=cvCreateMemStorage(0);
        CvSeq* firstContour=NULL;

        int a=cvFindContours(g_gray, storage, &firstContour,sizeof(CvChain),CV_RETR_LIST, CV_CHAIN_CODE);

        CvHistogram *hist;
        int h_bins = 30, s_bins = 30;
        float h_ranges[] = { 0, 180 };
        float s_ranges[] = { 0, 255 };
        int hist_size[] = { h_bins, s_bins };
        float* ranges[] = { h_ranges, s_ranges };

        hist = cvCreateHist( 2, hist_size, CV_HIST_ARRAY, ranges, 1 );
        if( CV_SEQ_ELTYPE( firstContour ) != CV_32SC2 ) {
                printf(":(\n");
        }

        cvCalcPGH((CvSeq *)firstContour,hist);
        return 0;

}

cvFindContours()将按照文档说明运行，但cvCalcPGH()将失败：

$ ./calcpgh 
:(
OpenCV Error: Unsupported format or combination of formats (The contour is not valid or the point type is not supported) in cvCalcPGH, file /Users/tonu-samuel/OpenCV-2.4.2/modules/legacy/src/pgh.cpp, line 351
terminate called throwing an exceptionAbort trap: 6
$

assert是2.4.2，并指出OpenCV ()用于我在代码中检查的CV_SEQ_ELTYPE( contour ) != CV_32SC2。如何修复它不确定。

Answer 1

谷歌不会给你好的例子的原因是，你正在谈论的所有东西都被弃用了。

cvCalcPGH在新的opencv中不再可用，你可能已经知道了，因为你已经包含了旧库。

另外，CV_CHAIN_CODE是removed，你不能在新的opencv中使用它。

你应该找时间将你的代码移植到新的opencv上，它正在流行，并将从其他人那里获得更多帮助。

但是，正如旧文档所说，讨论您当前的代码：

CV_CHAIN_CODE -输出Freeman链码中的轮廓。所有其他方法输出多边形(顶点序列)

cvFindContours的输出始终是顶点序列，CV_CHAIN_CODE只是一个例外。这就是为什么cvCalcPGH无法理解这些代码的原因。

而且，您可能已经使用其他近似方法测试过代码，您应该已经看到cvCalcPGH可以很好地处理多边形。

您应该自己将Freeman代码转换为多边形，然后可以使用任何您想要的opencv方法。但是，如果不这样做，我还没有看到一个方法可以解释Freeman代码本身。

我不确定这个转换到底应该如何完成，但是使用下面的代码，你可以提取代码并开始自己解释它，这段代码是从Sara Maher Mohammad偷来的。

CvChainPtReader reader;

cvStartReadChainPoints((CvChain*) firstContour, &reader);
for (int i = 0; i < firstContour->total; i++) {
    CV_READ_SEQ_ELEM(reader.code, (*((CvSeqReader*)&(reader))));
    printf("%d \n", reader.code);
}

但同样，我不建议使用过时的方法和函数，您可能会找到更好的方法来做您想做的事情。

只是猜测你可能只想找出两个Freeman代码之间的距离，你可能不需要计算PGH，而只是使用一些其他的距离函数，比如讨论的here。