大家好,我是小鸭酱,博客地址为:http://www.cnblogs.com/xiaoyajiang  
此篇博客实现了判定平面一点是否在给定多边形内部的功能。精确,性能优良,因为只包含加法和乘法运算,效率极高。
实现的C++源码如下,注意点在源码的相应注释中。
有错误和不好的地方还望各位批评指正。
 
#include<iostream>
using namespace std;
 
struct Point2f  // 用以进行简单测试的结构。opencv中也有同名的数据类型,注意其x和y与我们想的不一样。
{
    double x;
    double y;
};
 
bool onPolygonTest(Point2f polygon[], int num, Point2f test_point);  // 用以测试在多边形边上的退化情况
bool onLineTest(Point2f p, Point2f q, Point2f s);  // 用以测试在边上的退化情况
bool inPolygonTest(Point2f polygon[], int num, Point2f test_point);  // 一般情况
bool toLeftTest(Point2f p, Point2f q, Point2f s);  // 核心技术
double Area2(Point2f p, Point2f q, Point2f s);  // 核心技术实现
/*----------------------------------------------------------------------*/
/*   函数bool isInner(Point2f [], int num, Point2f, bool)               */
/*   输入                                                                                 */
/*         参数1 多边形顶点: Point2f polygon[]                             */
/*         参数2 多边形顶点个数                                                    */
/*         参数3 测试点:Point2f test_point                                   */
/*         参数4 是否包含边上的点:bool isIncludeLine                    */
/*   输出                                                                                */
/*         true:在多边形内;                                                      */
/*         false:不在多边形内                                                     */
/*  作者:小鸭酱的书签                                                            */ 
/*  日期:2016年9月28日                                                        */ 
/*--------------------------------------------------------------------*/
bool isInner(Point2f polygon[], int num, Point2f test_point, bool isIncludeLine)  
{
    if(num < 3)
        return false;
    if (isIncludeLine)  // 首先判断退化情况
        return (onPolygonTest(polygon, num, test_point) || inPolygonTest(polygon, num, test_point));
    else  // 一般情况
        return inPolygonTest(polygon, num, test_point);
}
 
bool onPolygonTest(Point2f polygon[], int num, Point2f test_point)  // 只要点在任意一条边上,都算作在内部
{
    int sum = 0;
    int i = 0;
    for(; i < num - 1; ++i)
    {
        sum += onLineTest(polygon[i], polygon[i+1], test_point);
        if (sum > 0)
            return true;
    }
    sum += onLineTest(polygon[i], polygon[0], test_point);
    return sum > 0;
}
bool onLineTest(Point2f p, Point2f q, Point2f s)  // 处理退化情况,判断点是否在线上;这里只处理多边形是没有旋转的矩形的情况。后期更新任意线段
{
    if(p.x == q.x && q.x == s.x && (s.y - p.y) * (s.y - q.y) <= 0 )
        return true;
    else if(q.y == p.y && p.y == s.y && (s.x - p.x) * (s.x - q.x) <= 0 )
        return true;
    else
        return false;
}
 
bool inPolygonTest(Point2f polygon[], int num, Point2f test_point)  // 根据点在每条边的left或者right的pattern来判断是否在内部
{
    int sum = 0;  // 用以记录有多少个left的pattern
    int i = 0;
    for(; i < num - 1; ++i)
    {
        sum += toLeftTest(polygon[i], polygon[i+1], test_point);
    }
    sum += toLeftTest(polygon[i], polygon[0], test_point);
 
    return (sum == 0 || sum == num) ? true : false;  // 很显然,当点在多边形的每条有向边的right或者left,它才会在内部
}
 
bool toLeftTest(Point2f p, Point2f q, Point2f s)  // 主要技术。有向面积值大于0则在有向边的左边
{
    return Area2(p, q, s) > 0;
}
 
double Area2(Point2f p, Point2f q, Point2f s)  // 可见其只包含加法和乘法运算,计算给定点和多边形的一条有向边所组成的三角形的2倍面积,我们是要使用其符号来判断在左边还是右边
{
 
    return
    p.x * q.y - p.y * q.x
    + q.x * s.y - q.y * s.x
    + s.x * p.y - s.y * p.x;
 
}
 
int main()                  // test
{
    Point2f polygon[4];
    polygon[0].x = 0;
    polygon[0].y = 0;
    polygon[1].x = 2;
    polygon[1].y = 0;
    polygon[2].x = 2;
    polygon[2].y = 1;
    polygon[3].x = 0;
    polygon[3].y = 1;
    Point2f test_point;
 
    test_point.x = 0;
    test_point.y = 2;
    int n = sizeof(polygon) / sizeof(Point2f);
    cout << "in = " << inPolygonTest(polygon, n, test_point) << endl;
    cout << "online = " << onPolygonTest(polygon, n, test_point) << endl;
    cout << "output = " << isInner(polygon, n, test_point, false) << endl;
    cout << n << endl << endl;
 
    return 0;
}

inPolygonTest学习和C++实现的更多相关文章

  1. 从直播编程到直播教育:LiveEdu.tv开启多元化的在线学习直播时代

    2015年9月,一个叫Livecoding.tv的网站在互联网上引起了编程界的注意.缘于Pingwest品玩的一位编辑在上网时无意中发现了这个网站,并写了一篇文章<一个比直播睡觉更奇怪的网站:直 ...

  2. Angular2学习笔记(1)

    Angular2学习笔记(1) 1. 写在前面 之前基于Electron写过一个Markdown编辑器.就其功能而言,主要功能已经实现,一些小的不影响使用的功能由于时间关系还没有完成:但就代码而言,之 ...

  3. ABP入门系列(1)——学习Abp框架之实操演练

    作为.Net工地搬砖长工一名,一直致力于挖坑(Bug)填坑(Debug),但技术却不见长进.也曾热情于新技术的学习,憧憬过成为技术大拿.从前端到后端,从bootstrap到javascript,从py ...

  4. 消息队列——RabbitMQ学习笔记

    消息队列--RabbitMQ学习笔记 1. 写在前面 昨天简单学习了一个消息队列项目--RabbitMQ,今天趁热打铁,将学到的东西记录下来. 学习的资料主要是官网给出的6个基本的消息发送/接收模型, ...

  5. js学习笔记:webpack基础入门(一)

    之前听说过webpack,今天想正式的接触一下,先跟着webpack的官方用户指南走: 在这里有: 如何安装webpack 如何使用webpack 如何使用loader 如何使用webpack的开发者 ...

  6. Unity3d学习 制作地形

    这周学习了如何在unity中制作地形,就是在一个Terrain的对象上盖几座小山,在山底种几棵树,那就讲一下如何完成上述内容. 1.在新键得项目的游戏的Hierarchy目录中新键一个Terrain对 ...

  7. 《Django By Example》第四章 中文 翻译 (个人学习,渣翻)

    书籍出处:https://www.packtpub.com/web-development/django-example 原作者:Antonio Melé (译者注:祝大家新年快乐,这次带来<D ...

  8. 菜鸟Python学习笔记第一天:关于一些函数库的使用

    2017年1月3日 星期二 大一学习一门新的计算机语言真的很难,有时候连函数拼写出错查错都能查半天,没办法,谁让我英语太渣. 关于计算机语言的学习我想还是从C语言学习开始为好,Python有很多语言的 ...

  9. 多线程爬坑之路-学习多线程需要来了解哪些东西?(concurrent并发包的数据结构和线程池,Locks锁,Atomic原子类)

    前言:刚学习了一段机器学习,最近需要重构一个java项目,又赶过来看java.大多是线程代码,没办法,那时候总觉得多线程是个很难的部分很少用到,所以一直没下决定去啃,那些年留下的坑,总是得自己跳进去填 ...

随机推荐

  1. Django 反向生成 从数据库生成Model

    Django 反向生成 从数据库生成Model 使用Django生成Model python manage.py inspectdb或python manage.py inspectdb > m ...

  2. Codeforces 358D Dima and Hares

    http://codeforces.com/contest/358/problem/D 题意:给出n个数,每个数取走的贡献与相邻的数有关,如果取这个数的时候,左右的数都还没被取,那么权值为a,如果左右 ...

  3. JAVA Stack栈和Heap堆的区别(转)

          首先分清楚Stack,Heap的中文翻译:Stack—栈,Heap—堆.         在中文里,Stack可以翻译为“堆栈”,所以我直接查找了计算机术语里面堆和栈开头的词语:      ...

  4. 5.1.1 读取Redis 数据

    Redis 服务器是Logstash 推荐的Broker选择,Broker 角色就意味会同时存在输入和输出两个插件. 5.1.1 读取Redis 数据 LogStash::Input::Redis 支 ...

  5. sourceforge软件下载方式

    访问http://www.mirrorservice.org/sites/download.sourceforge.net/pub/sourceforge/ 根据软件名称在文件列表中查找

  6. 面向切面编程(AOP)的理解

    在传统的编写业务逻辑处理代码时,我们通常会习惯性地做几件事情:日志记录.事务控制及权限控制等,然后才是编写核心的业务逻辑处理代码.当代码编写完成回头再看时,不禁发现,扬扬洒洒上百行代码中,真正用于核心 ...

  7. 机房收费系统中的Grid++Report报表设计器的应用

    在进行账单查询功能的时候我应用了Grid++Report报表设计器,下面我就为大家介绍一下,还望大家多多指点. 首先,在Grid++Report报表设计器中进行报表界面的设置.在属性编辑窗口中这里对报 ...

  8. 【软件技巧】Sublime Text为不同语法定义不同高亮

    Sublime Text默认的语法高亮已经非常美丽了,可是对于个别语言还是有些不爽. 默认高亮规则叫Monokai,能够从Preferences->Settings - Default中看到: ...

  9. 初步swift语言学习笔记9(OC与Swift杂)

    笔者:fengsh998 原文地址:http://blog.csdn.net/fengsh998/article/details/34440159 转载请注明出处 假设认为文章对你有所帮助,请通过留言 ...

  10. [跟我学spring学习笔记][IoC]

    IoC基础 什么是IoC Ioc—Inversion of Control,即“控制反转”,不是什么技术,而是一种设计思想. ioc做什么 IoC容器帮对象找相应的依赖对象并注入,而不是由对象主动去找 ...