最近在看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
结果
在线预览: