CGAL求最小外包矩形

有两种所谓的最小外包矩形，第一种通过求所有节点的最小与最大xy来求的，这种叫与坐标轴平行的最小外包矩形；另外一种则是本文说的这种，与范围的形状与走势有关的，叫非坐标轴平行的最小外包矩形，效果如下图所示:

CGAL求最小外包矩形 > image2018-8-21_10-33-22.png

/// 这里做测试，只选中LWPOLYLINLE，也就是CAD中的多段线实体

struct resbuf* rb = ads_buildlist(RTDXF0, _T("LWPOLYLINE"), RTNONE);

ads_name ent, ssLines;

if (RTNORM != ads_ssget(NULL, NULL, NULL, rb, ssLines))

{

       ads_relrb(rb);

      return;

}

ads_relrb(rb);

AcDbObjectId id;

Adesk::Int32 len = 0, l = 0;

ads_sslength(ssLines, &len);

std::vector<Point_2> nodesn,nodes_hull;

////////这里遍历选择集获取所有实体的节点

for (l = 0;l < len;l++)

{

      ads_ssname(ssLines, l, ent);

      if (Acad::eOk != acdbGetObjectId(id, ent))

      {

            continue;

      }

      AcDbEntity* pEnt = NULL;

      if (acdbOpenAcDbEntity(pEnt, id, AcDb::kForRead) != Acad::eOk)

      {

            continue;

      }

      if (!pEnt->isKindOf(AcDbPolyline::desc()))

      {

            continue;

      }

      AcDbPolyline* pLine = AcDbPolyline::cast(pEnt);

      AcGePoint3d ptS, ptE;

      for (int idx = 0;idx < pLine->numVerts();idx++)

      {

            AcGePoint2d pt;

            pLine->getPointAt(idx, pt);

            nodesn.push_back(Point_2(pt.x, pt.y));

            acutPrintf(_T("Random coordinate:%d: %.3lf %.3lf\n"), idx+1,pt.x, pt.y);

      }

      pEnt->close();

      break;

}

ads_ssfree(ssLines);

///发现在求最小外接矩形的时候传入的点不能有坐标一直的点，而且点的方向必须是逆时针排列

///所以这里先求节点的凸包，然后通过Polygon_2来把非逆时针排列的节点纠正成逆时针

CGAL::convex_hull_2(nodesn.begin(), nodesn.end(), std::back_inserter(nodes_hull));

Polygon_2 p;

std::vector<Point_2>::iterator itn;

for (itn = nodes_hull.begin();itn != nodes_hull.end();itn++)

{

      p.push_back(*itn);

}

AcGePoint2dArray nodes;

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

{

      double x = CGAL::to_double(p.vertex(i).x());

      double y = CGAL::to_double(p.vertex(i).y());

      nodes.append(AcGePoint2d(x, y));

}

Common::DrawPolyline2D(nodes, true);

if (p.is_clockwise_oriented())

{

      p.reverse_orientation();

}

Polygon_2 m_p;

CGAL::min_rectangle_2(p.vertices_begin(), p.vertices_end(), std::back_inserter(m_p));

Polygon_2::Vertex_const_iterator it;

nodes.removeAll();

for (it = m_p.vertices_begin(); it != m_p.vertices_end(); it++)

{

      double x = CGAL::to_double((*it).x());

      double y = CGAL::to_double((*it).y());

      nodes.append(AcGePoint2d(x, y));

}

Common::DrawPolyline2D(nodes, true);

巴特西

CGAL求最小外包矩形

最新文章

热门文章