• 博客园logo
  • 会员
  • 周边
  • 众包
  • 新闻
  • 博问
  • 闪存
  • 赞助商
  • Chat2DB
    • 搜索
      所有博客
    • 搜索
      当前博客
  • 写随笔 我的博客 短消息 简洁模式
    用户头像
    我的博客 我的园子 账号设置 会员中心 简洁模式 ... 退出登录
    注册 登录
FightingForWorldFinal
博客园    首页    新随笔    联系   管理    订阅  订阅

计算几何之极角排序

1.利用叉积的正负来作cmp.(即是按逆时针排序).

 bool cmp(const point &a, const point &b)//逆时针排序 
 {   
      point origin;
    origin.x = origin.y = 0;
    return cross(origin,b,origin,a) < 0;
}

2.利用complex的内建函数。

 #include<complex>
 bool cmp(const Point& p1, const Point& p2)//360度范围逆时针排序
 {
    complex<double> c1(p1.x,p1.y),c2(p2.x,p2.y);
     return arg(c1) < arg(c2);
 }

3.利用arctan计算极角大小。(范围『-180,180』)

bool cmp(const Point& p1, const Point& p2)
{     
    return atan2(p1.y, p1.x) < atan2(p2.y, p2.x);
}

4.利用象限加上极角,叉积。

 bool cmp2(const Point &a, const Point &b)//先按象限排序,再按极角排序,再按远近排序 
{
    if (a.y == 0 && b.y == 0 && a.x*b.x <= 0)return a.x>b.x;
    if (a.y == 0 && a.x >= 0 && b.y != 0)return true;
    if (b.y == 0 && b.x >= 0 && a.y != 0)return false;
    if (b.y*a.y <= 0)return a.y>b.y;
    Point one;
    one.y = one.x = 0;
    return crossDet(one,a,b) > 0 || (crossDet(one,a,b) == 0 && a.y > b.y);
}

顺便推荐个网站:http://www.csie.ntnu.edu.tw/~u91029/PointLinePlane2.html(上面讲的很详细)

posted @ 2013-10-24 20:09  Sky-J  阅读(713)  评论(0)    收藏  举报
刷新页面返回顶部
博客园  ©  2004-2025
浙公网安备 33010602011771号 浙ICP备2021040463号-3