C++-ObjC OpenCV 约束 Delaunay

Question

我在 OpenCV 2.3.1 中成功实现了轮廓的 Delaunay 三角剖分。

使用 cvPointPolygonTest 我可以获得凸包中的所有三角形，然后我尝试对三角形质心执行另一个 cvPointPolygonTest 以了解它们是否在主轮廓中，这样我就可以进行约束三角剖分轮廓。

但是，它效果不佳，因为某些三角形（例如，两条腿分开的步行者）“越过”一个洞。

有谁知道执行约束三角测量的方法。我考虑过凸面缺陷，但无法理解如何开始。

提前致谢 ！

Ben

---

事实上，这不是凸包缺陷问题，而是三角剖分问题。此图片将向您展示问题：

特别是在三角外壳的底部，您可以看到三角测量位于轮廓内和轮廓外，因为 OpenCV 正在对凸包进行三角测量。我想找到一种方法来对轮廓本身进行三角测量。

我发现了一些关于在轮廓本身中添加斯坦纳点的想法，但找不到从 OpenCV 开始的地方。

我的想法是：

• 测试三角形是否在轮廓内和轮廓外；
• 如果 true ：获取交点；
• 并将其添加到 cvSubdiv2D 中。

我这样说对吗？

感谢您的耐心和解答！

Answer 1

最后我在我的程序中实现了poly2tri库。

这是结果（OpenCV 获取轮廓，然后 poly2tri 执行三角测量）：

// -------------------- poly2tri --------------------


NSMutableArray* temp = [[NSMutableArray alloc] init];

vector<p2t::Triangle*> triangles;
vector< vector<p2t::Point*> > polylines;
vector<p2t::Point*> polyline;

for(int i = 0; i < contour32->total; i++ ) {
    CvPoint2D32f* pt32 = CV_GET_SEQ_ELEM(CvPoint2D32f, contour32, i);
    polyline.push_back(new p2t::Point((double)pt32->x, (double)pt32->y));
}

polylines.push_back(polyline);


p2t::CDT* cdt = new p2t::CDT(polyline);



// TODO -> holes with CV_RETR_TREE

// Triangulation !
cdt->Triangulate();

// On exporte nos triangles
triangles = cdt->GetTriangles();

for (int i = 0; i < triangles.size(); i++) {

    p2t::Triangle& t = *triangles[i];
    p2t::Point& a = *t.GetPoint(0);
    p2t::Point& b = *t.GetPoint(1);
    p2t::Point& c = *t.GetPoint(2);

    double x1 = (width / rWidth * (double)a.x) - (width / 2.f);
    double y1 = (height / rHeight * (double)a.y) - (height / 2.f);                   
    double x2 = (width / rWidth * (double)b.x) - (width / 2.f);
    double y2 = (height / rHeight * (double)b.y) - (height / 2.f);
    double x3 = (width / rWidth * (double)c.x) - (width / 2.f);
    double y3 = (height / rHeight * (double)c.y) - (height / 2.f);

    [temp addObject:[[NSArray arrayWithObjects:
                      [NSArray arrayWithObjects:
                       [NSNumber numberWithDouble:x1],
                       [NSNumber numberWithDouble:y1],
                       [NSNumber numberWithDouble:0.],
                       nil], 
                      [NSArray arrayWithObjects:
                       [NSNumber numberWithDouble:x2],
                       [NSNumber numberWithDouble:y2],
                       [NSNumber numberWithDouble:0.],
                       nil], 
                      [NSArray arrayWithObjects:
                       [NSNumber numberWithDouble:x3],
                       [NSNumber numberWithDouble:y3],
                       [NSNumber numberWithDouble:0.],
                       nil], 
                      nil] autorelease]];

}

[outDelaunay addObject:temp];                       

其中contour32 是使用 cvFindContours 找到的 OpenCV 轮廓，然后转换为 float32。

Answer 2

您必须使用：

    CvSeq* cvConvexityDefects(
const CvArr* contour,
const CvArr* convexhull,
CvMemStorage* storage = NULL
);

这将返回 CvConvexityDefect 的序列（CvSeq），而不是 CvPoint 或您可能习惯的任何其他类似 CvSeq 的内容。

您可以在循环中遍历缺陷，如下所示：

CvSeq* defects =  cvConvexityDefects(.....);
for (i = 0; i < defects->total; i++)
{
CvConvexityDefect* def = (CvConvexityDefect*)cvGetSeqElem(defects, i);
//do something with them
}

缺陷具有以下结构：

    typedef struct CvConvexityDefect {
// point of the contour where the defect begins
CvPoint* start;
// point of the contour where the defect ends
CvPoint* end;
// point within the defect farthest from the convex hull
CvPoint* depth_point;
// distance between the farthest point and the convex hull
float depth;
} CvConvexityDefect;

因此，在获取每个缺陷后，您可以从中构建点的 CvSeq 并在三角形的质心上使用 cvPointPolygonTest看看它们是否在里面。
如果它们位于缺陷内部，则意味着它们位于主轮廓外部。

希望这是您所需要的并且对您有所帮助。