C语言实现的PNPoly算法代码例子

 更新时间:2014年07月14日 09:21:22   投稿:junjie  
这篇文章主要介绍了C语言实现的PNPoly算法代码例子,PNPoly算法j是判断一个坐标点是否在不规则多边形内部的算法,需要的朋友可以参考下

写C语言的实验用到的一个算法,判断一个点是否在多边形的内部。C的代码如下:

int pnpoly(int nvert, float *vertx, float *verty, float testx, float testy)
{
  int i, j, c = 0;
  for (i = 0, j = nvert-1; i < nvert; j = i++) {
    if ( ((verty[i]>testy) != (verty[j]>testy)) && 
      (testx < (vertx[j]-vertx[i]) * (testy-verty[i]) / 
      (verty[j]-verty[i]) + vertx[i]) )
      c = !c;
    }
  return c;
}


其中nvert是多边形顶点的个数,vertx和verty分别是多边形顶点横、纵坐标的数组,textx和testy是待测点的坐标。这个算法是由W. Randolph Franklin提出的,根据Jordan curve theorem,多边形将平面分为内外两个区域,假设待测点在多边形内部,从待测点引出一条射线必然会与多边形有至少一个交点。该射线与多边形第一次相交时将“冲出”多边形,第二次相交将“进入”多边形,依此类推,若射线与多边形有奇数个交点,则该点在多边形内部,反之则在外部。

PNPoly算法正是从待测点引出一条水平向右的射线,并计算与多边形的交点个数。解释一下这段代码:for (i = 0, j = nvert-1; i < nvert; j = i++)循环的含义就是始终让j = i – 1,如果i = 0那么j = nvert – 1,从而依次检验多边形的每条边。接下来的重点就是条件语句,(verty[i]>testy) != (verty[j]>testy)很好理解,就是一条边上的两个顶点分别在待测点的上方和下方,通过这条语句可以知道从待测点向右引出的射线有可能与该条边相交(只要待测点在边的左侧即可)。

但具体判断相交就要交给解析几何了。建系写出该条边所在直线的方程:

变形一下:

代入待测点坐标,根据图形关系得到这个不等式:

也就是语句testx < vertx[j]-vertx[i]) * (testy-verty[i]) / (verty[j]-verty[i]) + vertx[i]了.
 

相关文章

  • C语言实现骑士飞行棋小游戏

    C语言实现骑士飞行棋小游戏

    这篇文章主要为大家详细介绍了C语言实现骑士飞行棋小游戏,文中示例代码介绍的非常详细,具有一定的参考价值,感兴趣的小伙伴们可以参考一下
    2020-02-02
  • C++重载运算符实现分数加减乘除

    C++重载运算符实现分数加减乘除

    这篇文章主要为大家详细介绍了C++重载运算符实现分数加减乘除,文中示例代码介绍的非常详细,具有一定的参考价值,感兴趣的小伙伴们可以参考一下
    2021-06-06
  • C++进程共享数据封装成类实例

    C++进程共享数据封装成类实例

    这篇文章主要介绍了C++进程共享数据封装成类的方法,以实例形式讲述了其封装代码与具体用法,具有一定的实用价值,需要的朋友可以参考下
    2014-10-10
  • C++你可能不知道地方小结

    C++你可能不知道地方小结

    c++中编译器替我们完成了许多事情,我们可能不知道,但也可能习以为常
    2013-01-01
  • C语言如何在指针中隐藏数据详解

    C语言如何在指针中隐藏数据详解

    这篇文章主要给大家介绍了关于C语言如何在指针中隐藏数据的相关资料,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的参考学习价值,需要的朋友们下面来一起看看吧
    2018-12-12
  • C语言版实现三子棋游戏

    C语言版实现三子棋游戏

    这篇文章主要为大家详细介绍了C语言版实现三子棋游戏,文中示例代码介绍的非常详细,具有一定的参考价值,感兴趣的小伙伴们可以参考一下
    2021-07-07
  • C/C++实现快速排序算法的思路及原理解析

    C/C++实现快速排序算法的思路及原理解析

    这篇文章主要介绍了C/C++实现快速排序算法的思路及原理解析,本文给大家介绍的非常详细,对大家的学习或工作具有一定的参考借鉴价值,需要的朋友可以参考下
    2021-01-01
  • map插入自定义对象总结

    map插入自定义对象总结

    黑树在插入节点时,必须依照大小比对之后在一个合适的位置上执行插入动作。所以作为关键字,起码必须有“<”这个比较操作符
    2013-09-09
  • 详解C++成员函数的override和final说明符的用法

    详解C++成员函数的override和final说明符的用法

    这篇文章主要介绍了详解C++成员函数的override和final说明符的用法,分别用于重写和禁止继承类,要的朋友可以参考下
    2016-01-01
  • C++空类详解

    C++空类详解

    以下是对C++中的空类进行了详细的介绍,需要的朋友可以过来参考下
    2013-09-09

最新评论