2D凸包算法(四):Quick Hull Algorithm

Quick Hull Algorithm 图示

在这里插入图片描述
先确定两个距离最大的点,连接后构成线。然后寻找离线的最远的点,构成三角形。以此类推,向外快速扩展,直到所有点都在凸包内。

  • 时间复杂度: O(N²) ,是最有效率的演算法。
    (预先剔除凸包内部大部分的点,而且不必预先排序所有点。)
具体思想:
  1. 找出最左点与最右点,连线,所有点分为上半部与下半部。(上半部与下半部分开求解)
  2. 处理上半部,找出上半凸包:
  3. 找到距离底线最远的点,形成一个三角形区域。(如果有许多点,就找最靠近左端点的那一点。避免凸包上多点共线。)
  4. 运用叉积运算,找出三角形外部的点,分为左、右两份。
  5. 至于三角形内部、三角形上的点,不会是凸包顶点,尽数剔除之。
  6. 左、右两份分别递归求解,直到找出上半凸包。(三角形的两边,分别做为左、右两份的底线。)
  7. 处理下半部,找出下半凸包,重复3-6步。

在这里插入图片描述

Quick Hull Algorithm 代码

目前测出当点数偏多时容易出现问题,code 仅供参考。

void quickhull(const vpdd &input)
{

    int min_x = 0, max_x = 0;
    for (int i = 0; i < input.size(); i++)
    {
        if (input[i].first < input[min_x].first)
        {
            min_x = i;
        }
        if (input[i].first > input[max_x].first)
        {
            max_x = i;
        }
    }

    hull_quick.push_back(input[min_x]);
    hull_quick.push_back(input[max_x]);

    vpdd left, right;
    for (int i = 0, temp; i < input.size(); i++)
    {
        if (i == min_x || i == max_x)
        {
            continue;
        }
        temp = getsign(input[min_x], input[max_x], input[i]);
        if (temp == 0)
        {
            continue;
        }
        if (temp < 0)
        {
            left.push_back(input[i]);
        }
        else
        {
            right.push_back(input[i]);
        }
    }
    qhullrecur(right, input[min_x], input[max_x]);
    qhullrecur(left, input[max_x], input[min_x]);
}

void qhullrecur(const vpdd &input, pdd a, pdd b)
{

    if (input.size() == 0)
    {
        return;
    }
    vpdd right_p;
    pid max_ind = make_pair(-1, 0);

    for (int i = 0; i < input.size(); i++)
    {
        if (get_right(a, b, input[i]))
        {
            if (max_ind.second < dist(a, b, input[i]))
            {
                max_ind.first = i;
                max_ind.second = dist(a, b, input[i]);
            }
        }
    }

    pdd c = input[max_ind.first];
    hull_quick.push_back(c);

    vpdd acr, bcr;
    for (int i = 0; i < input.size(); i++)
    {
        if (get_right(a, c, input[i]))
        {
            acr.push_back(input[i]);
        }
        if (get_right(c, b, input[i]))
        {
            bcr.push_back(input[i]);
        }
    }
    qhullrecur(acr, a, c);
    qhullrecur(bcr, c, b);
}
int getsign(pdd a, pdd b, pdd c)
{
    int d = (c.first - a.first) * (b.second - a.second) - (c.second - a.second) * (b.first - a.first);
    return d;
}

bool get_right(pdd a, pdd b, pdd c)
{
    return ((c.first - a.first) * (b.second - a.second) - (c.second - a.second) * (b.first - a.first)) > 0;
}

double dist(pdd a, pdd b, pdd c)
{
    return abs((c.first - a.first) * (b.second - a.second) - (c.second - a.second) * (b.first - a.first));
}

  • 4
    点赞
  • 24
    收藏
    觉得还不错? 一键收藏
  • 3
    评论
1 2/3维图像分割工具箱 2 PSORT粒子群优化工具箱 3 matlab计量工具箱Lesage 4 MatCont7p1 5 matlab模糊逻辑工具箱函数 6 医学图像处理工具箱 7 人工蜂群工具箱 8 MPT3安装包 9 drEEM toolbox 10 DOMFluor Toolbox v1.7 11 Matlab数学建模工具箱 12 马尔可夫决策过程(MDP)工具箱MDPtoolbox 13 国立SVM工具箱 14 模式识别与机器学习工具箱 15 ttsbox1.1语音合成工具箱 16 分数阶傅里叶变换的程序FRFT 17 魔方模拟器与规划求解 18 隐马尔可夫模型工具箱 HMM 19 图理论工具箱GrTheory 20 自由曲线拟合工具箱ezyfit 21 分形维数计算工具箱FracLab 2.2 22 For-Each 23 PlotPub 24 Sheffield大学最新遗传算法工具箱 25 Camera Calibration 像机标定工具箱 26 Qhull(二维三维三角分解、泰森图)凸包工具箱 2019版 27 jplv7 28 MatlabFns 29 张量工具箱Tensor Toolbox 30 海洋要素计算工具箱seawater 31 地图工具箱m_map 32 othercolor配色工具包 33 Matlab数学建模工具箱 34 元胞自动机 35 量子波函数演示工具箱 36 图像局域特征匹配工具箱 37 图像分割graphcut工具箱 38 NSGA-II工具箱 39 chinamap中国地图数据工具箱(大陆地区) 40 2D GaussFit高斯拟合工具箱 41 dijkstra最小成本路径算法 42 多维数据快速矩阵乘法 43 约束粒子群优化算法 44 脑MRI肿瘤的检测与分类 45 Matlab数值分析算法程序 46 matlab车牌识别完整程序 47 机器人工具箱robot-10.3.1 48 cvx凸优化处理工具箱 49 hctsa时间序列分析工具箱 50 神经科学工具箱Psychtoolbox-3-PTB 51 地震数据处理工具CREWES1990版 52 经济最优化工具箱CompEcon 53 基于约束的重构分析工具箱Cobratoolbox 54 Schwarz-Christoffel Toolbox 55 Gibbs-SeaWater (GSW)海洋学工具箱 56 光声仿真工具箱K-Wave-toolbox-1.2.1 57 语音处理工具箱Sap-Voicebox 58 贝叶斯网工具箱Bayes Net Toolbox(BNT) 59 计算机视觉工具箱VFfeat-0.9.21 60 全向相机校准工具箱OCamCalib_v3.0 61 心理物理学数据分析工具箱Palamedes1_10_3 62 生理学研究工具箱EEGLAB 63 磁共振成像处理工具箱CONN 18b 64 matlab 复杂网络工具箱 65 聚类分析工具箱FuzzyClusteringToolbox 66 遗传规划matlab工具箱 67 粒子群优化工具箱 68 数字图像处理工具箱DIPUM Toolbax V1.1.3 69 遗传算法工具箱 70 鱼群算法工具箱OptimizedAFSAr 71 蚁群算法工具箱 72 matlab优化工具箱 73 数据包络分析工具箱 74 图像分割质量评估工具包 75 相关向量机工具箱 76 音频处理工具箱 77 nurbs工具箱 78 Nurbs-surface工具箱 79 grabit数据提取工具箱 80 量子信息工具箱QLib 81 DYNAMO工具箱 82 NEDC循环的整车油耗量 83 PlotHub工具箱 84 MvCAT_Ver02.01 85 Regularization Tools Version 4.1 86 MatrixVB 4.5(含注册) 87 空间几何工具箱 matGeom-1.2.2 88 大数计算工具箱 VariablePrecisionIntegers 89 晶体织构分析工具包 mtex-5.7.0 90 Minimal Paths 2工具箱 91 Matlab数学建模工具箱

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

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

余额充值