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语言对归并排序算法的使用,归并排序算法的平均事件复杂度为(n\log n),需要的朋友可以参考下2016-05-05
这篇文章主要介绍了二分查找算法在C语言程序中的使用示例,文中最后提到了使用二分查找法一个需要注意的地方,需要的朋友可以参考下2021-08-08
我们今天要给大家讲的底层的IO看上去简单，实则抽象。并且在它之上衍生出了语言层面用于实战的技术，比如我们熟悉的java语言中的NIO或者像Netty这样的框架2021-06-06
这篇文章主要为大家详细介绍了C++实现教职工信息管理系统，文中示例代码介绍的非常详细，具有一定的参考价值，感兴趣的小伙伴们可以参考一下2022-03-03
线程池，简单来说就是有一堆已经创建好的线程，接下来通过本文给大家介绍C/C++ 原生API实现线程池的方法，感兴趣的朋友跟随小编一起看看吧2021-11-11
八分音符酱作为一款声音控制类游戏，当时还是很受大家的喜爱的。本文将用Matlab制作一款简易版的八分音符酱游戏，感兴趣的可以学习一下2022-02-02
这篇文章主要介绍了如何让Dev-C++支持auto关键字问题,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不吝赐教2023-12-12
程序在运行过程中，有对也就有错，正确那么就不用说了，但是如果错误，那么我们如何快速的定位到错误的位置，以及知道发生了什么错误。当一个函数发现自己无法处理的异常，就会抛出一个异常，让函数调用者直接或者间接的处理这个错误。本文将详解介绍c++中的异常2021-06-06
这篇文章主要为大家详细介绍了C语言实现密码本小项目，文中示例代码介绍的非常详细，具有一定的参考价值，感兴趣的小伙伴们可以参考一下2020-02-02
这篇文章主要介绍了C++基于栈实现铁轨问题,实例分析了C++使用栈实现铁轨问题的思路与解决方法,具有一定参考借鉴价值,需要的朋友可以参考下2015-04-04
最新评论