自己写的点云点集凹包求解 ,不是凸包!!!(附C++代码)

近期遇到一个计算几何问题,需要从点集中重建一个粗略合理的几何形状。求出一定程度上反映这些散点轮廓的平面多边形,给出边的连接方式即可。如从下图的左图散点重建为右图的形状:

网上流传的版本都是求凸包,而不是求凹包,于是自己写一个版本的代码,有瑕疵,但是基本可以用;

思路:随机取一个点作为原点A,依次计算A点与其余各点的连线的正切值,以A点为原点,分割为多个正切值的区间,每个正切值的区间内取一个最远的点作为边缘点,依次连接边缘点即组成一个粗略的凹包图形。

如上图,以A点作为原点,把它的360区域按照正切值分为了4个区域,每个区域找最远的点作为边缘点。这里只分了4个区域作为示意图,实际可以分很多个区域,效果会好很多。

C++ 的代码,如下:

vector<vector <cv::Point>> contour_s;
    vector <float> threshold;
    int bin_num = 30;
    for (int i = 0; i <= bin_num; i
  • 1
    点赞
  • 3
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 0
    评论
好的,我可以为您解答这个问题。 对于凸包问题,我们可以使用分治法来解决。具体步骤如下: 1. 将所有点按照 x 坐标进行排序; 2. 将点集分成左右两组,分别求解每一组的凸包; 3. 合并左右两组的凸包。 接下来是具体的 C++ 代码实现: ```c++ #include <iostream> #include <vector> #include <algorithm> using namespace std; struct Point { int x; int y; }; // 按照 x 坐标排序 bool cmp(Point a, Point b) { return a.x < b.x; } // 叉积计算 int crossProduct(Point a, Point b, Point c) { return (b.x - a.x) * (c.y - a.y) - (c.x - a.x) * (b.y - a.y); } // 求解凸包 vector<Point> convexHull(vector<Point> points) { int n = points.size(); sort(points.begin(), points.end(), cmp); vector<Point> hull; for (int i = 0; i < n; i++) { while (hull.size() >= 2 && crossProduct(hull[hull.size() - 2], hull.back(), points[i]) < 0) { hull.pop_back(); } hull.push_back(points[i]); } int t = hull.size(); for (int i = n - 2; i >= 0; i--) { while (hull.size() > t && crossProduct(hull[hull.size() - 2], hull.back(), points[i]) < 0) { hull.pop_back(); } hull.push_back(points[i]); } hull.pop_back(); return hull; } // 合并凸包 vector<Point> mergeConvexHull(vector<Point> left, vector<Point> right) { int leftmost = 0, rightmost = 0; for (int i = 1; i < left.size(); i++) { if (left[i].x < left[leftmost].x) { leftmost = i; } } for (int i = 1; i < right.size(); i++) { if (right[i].x > right[rightmost].x) { rightmost = i; } } vector<Point> hull; hull.push_back(left[leftmost]); hull.push_back(right[rightmost]); int l = leftmost, r = rightmost; while (l < left.size() - 1 || r > 0) { if (l == left.size() - 1) { hull.push_back(right[--r]); } else if (r == 0) { hull.push_back(left[++l]); } else { int cp = crossProduct(left[l], right[r], right[r - 1]); if (cp > 0) { hull.push_back(right[--r]); } else { hull.push_back(left[++l]); } } } return hull; } // 分治求解凸包 vector<Point> divideAndConquerConvexHull(vector<Point> points) { if (points.size() <= 1) { return points; } int mid = points.size() / 2; vector<Point> left(points.begin(), points.begin() + mid); vector<Point> right(points.begin() + mid, points.end()); vector<Point> hullLeft = divideAndConquerConvexHull(left); vector<Point> hullRight = divideAndConquerConvexHull(right); return mergeConvexHull(hullLeft, hullRight); } int main() { vector<Point> points = {{0, 0}, {1, 1}, {2, 2}, {3, 3}, {2, 4}, {3, 5}, {4, 3}, {5, 2}, {4, 0}, {3, 1}}; vector<Point> hull = divideAndConquerConvexHull(points); for (int i = 0; i < hull.size(); i++) { cout << "(" << hull[i].x << ", " << hull[i].y << ")" << endl; } return 0; } ``` 以上代码实现了分治法求解凸包问题,其中 `Point` 结构体表示一个二维点,`cmp` 函数用于按照 x 坐标排序,`crossProduct` 函数用于计算叉积,`convexHull` 函数用于求解单组点集凸包,`mergeConvexHull` 函数用于合并左右两组点集凸包,`divideAndConquerConvexHull` 函数用于利用分治法求解凸包问题。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

点云-激光雷达-Slam-三维牙齿

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值