c语言实现点在多边形内部,C语言中实现 点在多边形内 算法
来源:
天极网
作者:
若水
2008-05-15/01:29
本文是采用射线法判断点是否在多边形内的C语言程序。多年前,我自己实现了这样一个算法。但是随着时间的推移,我决定重写这个代码。参考周培德的《计算几何》一书,结合我的实践和经验,我相信,在这个算法的实现上,这是你迄今为止遇到的最优的代码。
这是个C语言的小算法的实现程序,本来不想放到这里。可是,当我自己要实现这样一个算法的时候,想在网上找个现成的,考察下来竟然一个符合需要的也没有。我对自己大学读书时写的代码没有信心,所以,决定重新写一个,并把它放到这里,以飨读者。也增加一下BLOG的点击量。
首先定义点结构如下:
以下是引用片段:
/* Vertex structure */
typedef struct
{
double x, y;
} vertex_t;
本算法里所指的多边形,是指由一系列点序列组成的封闭简单多边形。它的首尾点可以是或不是同一个点(不强制要求首尾点是同一个点)。这样的多边形可以是任意形状的,包括多条边在一条绝对直线上。因此,定义多边形结构如下:
以下是引用片段:
/* Vertex list structure – polygon */
typedef struct
{
int num_vertices; /* Number of vertices in list */
vertex_t *vertex; /* Vertex array pointer */
} vertexlist_t;
为加快判别速度,首先计算多边形的外包矩形(rect_t),判断点是否落在外包矩形内,只有满足落在外包矩形内的条件的点,才进入下一步的计算。为此,引入外包矩形结构rect_t和求点集合的外包矩形内的方法vertices_get_extent,代码如下:
以下是引用片段:
/* bounding rectangle type */
typedef struct
{
double min_x, min_y, max_x, max_y;
} rect_t;
/* gets extent of vertices */
void vertices_get_extent (const vertex_t* vl, int np, /* in vertices */
rect_t* rc /* out extent*/ )
{
int i;
if (np > 0){
rc->min_x = rc->max_x = vl[0].x; rc->min_y = rc->max_y = vl[0].y;
}else{
rc->min_x = rc->min_y = rc->max_x = rc->max_y = 0; /* =0 ? no vertices at all */
}
for(i=1; i
{
if(vl[i].x min_x) rc->min_x = vl[i].x;
if(vl[i].y min_y) rc->min_y = vl[i].y;
if(vl[i].x > rc->max_x) rc->max_x = vl[i].x;
if(vl[i].y > rc->max_y) rc->max_y = vl[i].y;
}
}
当点满足落在多边形外包矩形内的条件,要进一步判断点(v)是否在多边形(vl:np)内。本程序采用射线法,由待测试点(v)水平引出一条射线B(v,w),计算B与vl边线的交点数目,记为c,根据奇内偶外原则(c为奇数说明v在vl内,否则v不在vl内)判断点是否在多边形内。
具体原理就不多说。为计算线段间是否存在交点,引入下面的函数:
(1)is_same判断2(p、q)个点是(1)否(0)在直线l(l_start,l_end)的同侧;
(2)is_intersect用来判断2条线段(不是直线)s1、s2是(1)否(0)相交;
#p#分页标题#e#以下是引用片段:
/* p, q is on the same of line l */
static int is_same(const vertex_t* l_start, const vertex_t* l_end, /* line l */
const vertex_t* p,
const vertex_t* q)
{
double dx = l_end->x - l_start->x;
double dy = l_end->y - l_start->y;
double dx1= p->x - l_start->x;
double dy1= p->y - l_start->y;
double dx2= q->x - l_end->x;
double dy2= q->y - l_end->y;
return ((dx*dy1-dy*dx1)*(dx*dy2-dy*dx2) > 0? 1 : 0);
}
/* 2 line segments (s1, s2) are intersect? */
static int is_intersect(const vertex_t* s1_start, const vertex_t* s1_end,
const vertex_t* s2_start, const vertex_t* s2_end)
{
return (is_same(s1_start, s1_end, s2_start, s2_end)==0 &&
is_same(s2_start, s2_end, s1_start, s1_end)==0)? 1: 0;
}
下面的函数pt_in_poly就是判断点(v)是(1)否(0)在多边形(vl:np)内的程序:
以下是引用片段:
int pt_in_poly ( const vertex_t* vl, int np, /* polygon vl with np vertices */
const vertex_t* v)
{
int i, j, k1, k2, c;
rect_t rc;
vertex_t w;
if (np
return 0;
vertices_get_extent(vl, np, &rc);
if (v->x x > rc.max_x || v->y y > rc.max_y)
return 0;
/* Set a horizontal beam l(*v, w) from v to the ultra right */
w.x = rc.max_x + DBL_EPSILON;
w.y = v->y;
c = 0; /* Intersection points counter */
for(i=0; i
{
j = (i+1) % np;
if(is_intersect(vl+i, vl+j, v, &w))
{
c++;
}
else if(vl[i].y==w.y)
{
k1 = (np+i-1)%np;
while(k1!=i && vl[k1].y==w.y)
k1 = (np+k1-1)%np;
k2 = (i+1)%np;
while(k2!=i && vl[k2].y==w.y)
k2 = (k2+1)%np;
if(k1 != k2 && is_same(v, &w, vl+k1, vl+k2)==0)
c++;
if(k2 <= i)
break;
i = k2;
}
}
return c%2;
}
本想配些插图说明问题,但是,CSDN的文章里放图片我还没用过。以后再试吧!实践证明,本程序算法的适应性极强。但是,对于点正好落在多边形边上的极端情形,有可能得出2种不同的结果。
c语言实现点在多边形内部,C语言中实现 点在多边形内 算法相关推荐
- c语言实现点在多边形内部,求教这个判断点是否在多边形内的C语言算法用的什么原理....
该楼层疑似违规已被系统折叠 隐藏此楼查看此楼 #include //#define N const int N=4; int tt(double p[][2],int a,int b) //a为下标, ...
- 点在多边形内算法,C#判断一个点是否在一个复杂多边形的内部
判断一点是否在不规则图像的内部算法,如下图是由一个个点组成的不规则图像,判断某一点是否在不规则矩形内部,先上效果图 算法实现如下,算法简单,亲试有效 public class PositionAlgo ...
- 向量叉乘判断点是否在多边形内部
向量叉乘: 二维向量的叉乘 ( x1 , y1 ) * ( x2 , y2 ) = x1y2 - y1x2 如果值大于0 , 则表明 ( x2 , y2 ) 在 ( x1 , y1 )左边,反之在右边 ...
- 求任意多边形内部水平方向似最大矩形算法实现
背景说明 前段时间有个求点是否在多边形内部的需求,折腾了不少时间,现截取其中的的重点部分--求任意多边形内部水平方向似最大矩形--来搞篇博客. 求点是否在多边形内部这个算法很容易搞,一搜一大把,但数据 ...
- R语言ggplot2可视化散点图、并使用geom_encircle函数自定义多边形圈定可视化图像中的指定区域、使用geom_smooth函数基于loess方法拟合数据点之间的趋势关系曲线
R语言ggplot2可视化散点图.并使用geom_encircle函数自定义多边形圈定可视化图像中的指定区域(Scatterplot With Encircling).使用geom_smooth函数基 ...
- iOS开发之APP内部切换语言
2019独角兽企业重金招聘Python工程师标准>>> 前言 iOS开发中,随着APP越来越完善,很多APP都做了国际化,也实现了APP内部切换语言. 原理 国际化都会走到NSBun ...
- 判定点是否在不规则多边形内部的问题
2019独角兽企业重金招聘Python工程师标准>>> 问题如下: 话说在平面内有一个任意的不规则的封闭多边形,另外在这个平面内还有一个点,问题:如何高效的判定这个点是在这个多边形内 ...
- 【计算几何】点在多边形内部
问题描述:已知点P(x,y)和多边形Poly,判断点P(x,y)是否在多边形内部. 基本方法:射线法 以点P为端点,向左方作射线L,由于多边形是有界的,所以射线L的左端一定在多边形外部,考虑沿着L从无 ...
- [转]如何判断一个点是否在一个多边形内部
原文地址:http://hi.baidu.com/wayright/item/ad18e4c0c5446b2dee4665c9 在多边性的存储中,每一个多边形都是由一系列连续的点组成,例如保存为数组P ...
最新文章
- android 下的网络图片加载
- PySpark安装和测试
- 老男孩Linux运维第41期20171016第六周学习重点课堂记录
- 成功解决numpy.linalg.LinAlgError: singular matrix
- ai外呼营销系统_郑州专业智能营销系统开发,外呼电话线路好用吗-有效!-「线路稳定 不封号」...
- python和javascript哪个好_Python
- 微信公众号怎么快速导出一个月的文章数据
- Windows/Linux如何查看动态库的依赖关系?
- 计算机三级网络技术总结,计算机三级网络技术总结
- JavaScript数据类型之Number
- 计算机必懂的54个英文单词和缩写
- 【气动学】基于matlab GUI外弹道仿真系统【含Matlab源码 1044期】
- esp ghost引导_UEFI引导修复教程和工具
- centos安装VirtualBox增强包VBoxGuestAdditions
- netty实现gmssl_gmssl国密总结
- 我的世界java甘蔗机_我的世界甘蔗机怎么做? MineCraft甘蔗机图文教程
- arkit 人脸捕捉_iPhone X上的ARKit人脸追踪
- 因特尔显卡自定义分辨率_为什么从最新的英特尔显卡控制面板中移除了自定义分辨率?...
- 日语学习网站分类汇总
- 第十二章、GDK学习
热门文章
- 基于AVR和MT8870的远程家电控制系统设计
- 关于webcontrols的TreeView中转义符的处理问题
- 用神经网络做分子模型:乙烯和乙炔的实验数据
- 【Paper】2019_Consensus Control of Multiple AUVs Recovery System Under Switching Topologies and Time D
- 【Python】调用百度云API驾驶行为分析 Driver Behavior
- 2.10 m 个样本的梯度下降-深度学习-Stanford吴恩达教授
- TCL with SNPS file existsfile rename/copy redirect lsearch split
- Android + kernel USB host 如何判断插入的设备是高速还是低速的设备
- CentOS8更换源后失败问题
- redis 高级特性一