本帖最后由 破晓 于 2015-1-20 10:41 编辑
最近在看delaunay三角剖分,其中第一步就是求点集的凸包。凸包,即最外层的点连接起来构成的凸多边型。 查了下资料 凸包算法很容易。 步骤为: 输入的点集合必须没有相同的点 这是前提。 1.找到点集x最小点 x相同的情况 找y最小点。 2.将其余点 和 此点连线, 求直线与y轴夹角 逆时针排序。
3.然后用一个list, 先push入最左边的点(x坐标最小的点),再把点 0- 7依次push入,最后再一次把最左边的点push入。 用一个新list result, push 入 最左边的点, 点 0 点 1. 然后从点2开始 检查 向量01 12是否为左转,即检查直线 01 12 两者形成的角度是否向左。(可以理解角度朝向最左边的点)
向量(x1, y1) (x2, y2) 左转检测 即两者叉积 x1y2 - x2y1 < 0 如果 012为左转 result push入点2 继续下个点 检查 123 是否左转 直到 list遍历完毕。 这里实际上123 非左转, 所以 我们删掉 点2 然后 连接 点 013 看是否左转。如果点013也不左转 则删掉 点 1。 重复 直到 list遍历完毕。
所谓左转 即:
如图 01 12 为左转, 12 23非左转。
检测的核心代码: - var t:int = -1;
- result[++t] = list[0];
- result[++t] = list[1];
- result[++t] = list[2];
- for(i = 3; i<len; i++)
- {
- while(!left(result[t-1], result[t], list[i] )) t--;
- result[++t] = list[i];
- }
复制代码
上面 i = 3 为 图中点2 结果
在线预览:(用鼠标点地两分割线中间的空白区域绘制图形)
|