polygon.h测试

最新推荐文章于 2022-06-13 09:43:27 发布
最新推荐文章于 2022-06-13 09:43:27 发布
阅读量1.2k 收藏
点赞数
分类专栏： C++、GO端游、手游开发
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/u010401391/article/details/9013703
版权
这是一个关于C++代码的测试，主要涉及polygon.h库中用于计算几何的函数，包括判断圆是否与多边形相交、线段是否与多边形相交等。代码提供了主函数main中不同场景的测试用例，如判断点在多边形内的算法InsidePolygon和CircleCrossPolygon。
摘要由CSDN通过智能技术生成

测试代码：
#include "stdafx.h"
#include "polygon.h"
#include <iostream>
using namespace std;

const int N=4;

/* 判断圆是否与多边形相交*/   
//bool CircleCrossConvexPolygon(int vcount,POINT center,float radius,POINT polygon[])
bool CircleCrossPolygon(int vcount,POINT center,float radius,POINT polygon[])
{   
 POINT P_0=POINT(center.x+radius,center.y);
 int N=10;
 for(int i=0;i<N;i++)
 {
  POINT P_i=rotate(center,i*2*PI/N,P_0);
  //bool ret=InsideConvexPolygon(vcount,polygon,P_i);
  bool ret=InsidePolygon(vcount,polygon,P_i);
  if(ret)
   return true;
 }
 P_0=rotate(center,PI/N,P_0);
 for(int i=0;i<N;i++)
 {
  POINT P_i=rotate(center,i*2*PI/N,P_0);
  //bool ret=InsideConvexPolygon(vcount,polygon,P_i);
  bool ret=InsidePolygon(vcount,polygon,P_i);
  if(ret)
   return true;
 }
 return false; 
}

/* 判断线段LINESEG laser所在直线LINE laser是否与简单多边形polygon相交*/
bool LineCrossPolygon(int vcount,LINESEG laser,POINT polygon[])
{
 LINESEG u[MAXV];
 for(int i=0;i<vcount-1;i++)
  u[i]=LINESEG(polygon[i],polygon[i+1]);
 if(vcount>2)
  u[vcount-1]=LINESEG(polygon[vcount-1],polygon[0]);
 for(int i=0;i<vcount;i++)
 {
  bool ret=intersect_l(u[i],laser);
  if(ret)
   return true;
 }
 return false;
}

int main(int argc,char* argv[])
{
 不能构成凸多边形，但是简单多边形
 //CPoint mathPointSet[4]={CPoint(0,0),CPoint(100,200),CPoint(300,400),CPoint(-200,600)};
 不能构成凸多边形，但是简单多边形
 //CPoint mathPointSet[4]={CPoint(0,0),CPoint(300,400),CPoint(100,200),CPoint(-200,600)};
 //能构成凸多边形，且是凸多边形
 //CPoint mathPointSet[4]={CPoint(0,0),CPoint(400,200),CPoint(300,400),CPoint(-200,600)};
 //能构成凸多边形，但不是凸多边形,不是简单多边形
 CPoint mathPointSet[4]={CPoint(0,0),CPoint(300,400),CPoint(400,200),CPoint(-200,600)};

 //bool ret1=InsideConvexPolygon(sizeof(mathPointSet)/sizeof(CPoint),mathPointSet,CPoint(150,290));
 bool ret1=InsidePolygon(sizeof(mathPointSet)/sizeof(CPoint),mathPointSet,CPoint(150,290));
 //bool ret1=InsidePolygon(sizeof(mathPointSet)/sizeof(CPoint),mathPointSet,CPoint(230,290));
 cout<<"ret1="<<ret1<<endl;

 bool ret=CircleCrossPolygon(sizeof(mathPointSet)/sizeof(CPoint),CPoint(500,300),150,mathPointSet);
 //bool ret=CircleCrossPolygon(sizeof(mathPointSet)/sizeof(CPoint),CPoint(500,300),130,mathPointSet);
 cout<<"ret="<<ret<<endl;

 system("pause");
 return 0;

 float p[N][2];
 for(int i=0;i<N;i++)
 {
  p[i][0]=mathPointSet[i].x;
  p[i][1]=mathPointSet[i].y;
 }
 if (isconvexEx(p,N))
  cout<<"能";
 else
  cout<<"不能";

 CPoint mathPointSet2[4]={};

 //测试顶点顺序
 int len=0;
 cout<<endl;
 cout<<"排序前";
 cout<<endl;
 for(int i=0;i<N;i++)
 {
  cout<<" 第"<<i<<"个点坐标:";
  cout<<mathPointSet[i].x<<","<<mathPointSet[i].y<<endl;
 }
 if (issimple(N,mathPointSet))
  cout<<"无序的多边形顶点集合pt,issimple=true";
 else
  cout<<"无序的多边形顶点集合pt,issimple=false";
 cout<<endl;
 if (isconvex(N,mathPointSet))
  cout<<"是凸多边形";
 else
  cout<<"不是凸多边形";
 Graham_scan(mathPointSet,mathPointSet2,N,len);
 cout<<endl;
 cout<<"排序后";
 cout<<endl;
 for(int i=0;i<N;i++)
 {
  cout<<" 第"<<i<<"个点坐标:";
  cout<<mathPointSet[i].x<<","<<mathPointSet[i].y<<endl;
 }
 if (issimple(N,mathPointSet))
  cout<<"排序成逆时针顺序后的顶点集合pt，issimple=true";
 else
  cout<<"排序成逆时针顺序后的顶点集合pt，issimple=false";
 cout<<endl;
 if (isconvex(N,mathPointSet))
  cout<<"是凸多边形";
 else
  cout<<"不是凸多边形";
 cout<<endl;
 system("pause");
 return 0;
}
//int main(int argc,char* argv[])
//{
//   float p[N][2];
//   char ch='y';
//   while (ch=='y'||ch=='Y')
//   {
//    cout<<"请输入这"<<N<<"个点坐标:"<<endl;
//    for(int i=0;i<N;i++)
//    {
//     cout<<" 第"<<i<<"个点坐标:";
//    cin>>p[i][0]>>p[i][1];
//    //cin>>**p>>*(*p+1);
//    cout<<p[i][0]<<","<<p[i][1]<<endl;
//    }
//    if (isconvexEx(p,N))
//     cout<<"能";
//    else
//     cout<<"不能";
//    POINT pt[N],pt2[N];
//    for(int i=0;i<N;i++)
//    {
//     pt[i].x=(float)p[i][0];
//     pt[i].y=(float)p[i][1];
//    }
//    //测试顶点顺序
//    int len=0;
//    cout<<endl;
//    cout<<"排序前";
//    for(int i=0;i<N;i++)
//    {
//     cout<<" 第"<<i<<"个点坐标:";
//     cout<<pt[i].x<<","<<pt[i].y<<endl;
//    }
//    if (issimple(N,pt))
//     cout<<"无序的多边形顶点集合pt,issimple=true";
//    else
//     cout<<"无序的多边形顶点集合pt,issimple=false";
//    cout<<endl;
//    if (isconvex(N,pt))
//     cout<<"是凸多边形";
//    else
//     cout<<"不是凸多边形";
//    Graham_scan(pt,pt2,N,len);
//    cout<<endl;
//    cout<<"排序后";
//    for(int i=0;i<N;i++)
//    {
//     cout<<" 第"<<i<<"个点坐标:";
//     cout<<pt[i].x<<","<<pt[i].y<<endl;
//    }
//    if (issimple(N,pt))
//     cout<<"排序成逆时针顺序后的顶点集合pt，issimple=true";
//    else
//     cout<<"排序成逆时针顺序后的顶点集合pt，issimple=false";
//    cout<<endl;
//    if (isconvex(N,pt))
//     cout<<"是凸多边形";
//    else
//     cout<<"不是凸多边形";
//    cout<<endl<<"是否继续判断?(继续请输入y&Y)";
//    cin>>ch;
//   }
// system("pause");
// return 0;
//}


C++计算几何库polygon.h的代码：
/* 
计算几何目录----已测试算法模块  
㈠ 点的基本运算  
1. 平面上两点之间距离 1  
2. 判断两点是否重合 1  
3. 矢量叉乘 1  
4. 矢量点乘 2  
5. 判断点是否在线段上 2  
6. 求一点饶某点旋转后的坐标 2  
7. 求矢量夹角 2  

㈡ 线段及直线的基本运算  
1. 点与线段的关系 3  
2. 求点到线段所在直线垂线的垂足 4  
3. 点到线段的最近点 4  
4. 点到线段所在直线的距离 4  
5. 点到折线集的最近距离 4  
6. 判断圆是否在多边形内 5  
7. 求矢量夹角余弦 5  
8. 求线段之间的夹角 5  
9. 判断线段是否相交 6  
10.判断线段是否相交但不交在端点处 6  
11.求线段所在直线的方程 6  
12.求直线的斜率 7  
13.求直线的倾斜角 7  
14.求点关于某直线的对称点 7  
15.判断两条直线是否相交及求直线交点 7  
16.判断线段是否相交，如果相交返回交点 7  

㈢ 多边形常用算法模块  
1. 判断多边形是否简单多边形----issimple（POINT polygon[]理解为一个确定的多边形对象或有序的顶点集，结合Graham_scan函数使用） 
2. 检查多边形顶点的凸凹性 9  
3. 判断多边形是否凸多边形----isconvex（POINT polygon[]理解为一个确定的多边形对象或有序的顶点集，结合Graham_scan函数使用），isconvexEx（POINT polygon[]理解为一个确定的简单多边形对象或无序的顶点集）  
4. 求多边形面积 9  
5. 判断多边形顶点的排列方向，方法一 10  
6. 判断多边形顶点的排列方向，方法二 10  
7. 射线法判断点是否在多边形内 10  
8. 判断点是否在凸多边形内----InsideConvexPolygon
8a.判断点是否在平面多边形(可以是非简单多边形，可以是非凸多边形)内----InsidePolygon
9. 寻找点集的graham算法----Graham_scan
10.寻找点集凸包的卷包裹法 13  
11.判断线段是否在多边形内 14  
12.求简单多边形的重心 15  
13.求凸多边形的重心 17  
14.求肯定在给定多边形内的一个点 17  
15.求从多边形外一点出发到该多边形的切线 18  
16.判断多边形的核是否存在 19  

㈣ 圆的基本运算  
1 .点是否在圆内 20  
2 .求不共线的三点所确定的圆 21  

㈤ 矩形的基本运算  
1.已知矩形三点坐标，求第4点坐标 22  

㈥ 常用算法的描述 22  

㈦ 补充  
1．两圆关系： 24  
2．判断圆是否在矩形内： 24  
3．点到平面的距离： 25  
4．点是否在直线同侧： 25  
5．镜面反射线： 25  
6．矩形包含： 26  
7．两圆交点： 27  
8．两圆公共面积： 28  
9. 圆和直线关系： 29  
10. 内切圆： 30  
11. 求切点： 31  
12. 线段的左右旋： 31  
13．公式： 32  
*/
#pragma once
#include<cmath>
#include<xutility>
using namespace std;

/* 常用的常量定义 */   
const float    INF     = 1E20;//无穷大      
const float    EP      = 1E-10;//无穷小   
const int       MAXV    = 300;   
const float    PI      = 3.14159265; 

//float max(float x, float y)   
//{   
// return (x > y ? x : y);   
//}   
//
//float min(float x, float y)   
//{   
// return (x < y ? x : y);   
//}
  
/* 基本几何结构 */   
//struct POINT   
//{   
//    float x;   
//    float y;   
//    POINT(float a=0, float b=0) { x=a; y=b;} //constructor    
//};
struct qcPoint
{
 float x, y;
 qcPoint() : x(0), y(0) {}
 qcPoint(float _x, float _y) : x(_x), y(_y) {}
 qcPoint& operator= (const qcPoint& point) {x = point.x; y = point.y; return *this;}
};

typedef qcPoint TPoint;
typedef qcPoint POINT;
typedef qcPoint CPoint;

struct LINESEG   
{   
    POINT s;   
    POINT e;   
    LINESEG(POINT a, POINT b) { s=a; e=b;}   
    LINESEG() { }   
};   
struct LINE           // 直线的解析方程 a*x+b*y+c=0  为统一表示，约定 a >= 0    
{   
   float a;   
   float b;   
   float c;   
   LINE(float d1=1, float d2=-1, float d3=0) {a=d1; b=d2; c=d3;}   
};   
  
/********************** 
 *                    *  
 *   点的基本运算     *  
 *                    *  
 **********************/   
  
float dist(POINT p1,POINT p2)                // 返回两点之间欧氏距离    
{   
    return( sqrt( (p1.x-p2.x)*(p1.x-p2.x)+(p1.y-p2.y)*(p1.y-p2.y) ) );   
}   
bool equal_point(POINT p1,POINT p2)           // 判断两个点是否重合     
{   
    return ( (abs(p1.x-p2.x)<EP)&&(abs(p1.y-p2.y)<EP) );   
}   
/******************************************************************************  
r=multiply(sp,ep,op),得到(sp-op)和(ep-op)的叉积  
r>0：ep在矢量opsp的逆时针方向；  
r=0：opspep三点共线；  
r<0：ep在矢量opsp的顺时针方向  
*******************************************************************************/   
float multiply(POINT sp,POINT ep,POINT op)   
{   
    return((sp.x-op.x)*(ep.y-op.y)-(ep.x-op.x)*(sp.y-op.y));   
}   
/*  
r=dotmultiply(p1,p2,op),得到矢量(p1-op)和(p2-op)的点积，如果两个矢量都非零矢量  
r<0：两矢量夹角为锐角； 
r=0：两矢量夹角为直角； 
r>0：两矢量夹角为钝角  
*******************************************************************************/   
float dotmultiply(POINT p1,POINT p2,POINT p0)   
{   
    return ((p1.x-p0.x)*(p2.x-p0.x)+(p1.y-p0.y)*(p2.y-p0.y));   
}   
/******************************************************************************  
判断点p是否在线段l上 
条件：(p在线段l所在的直线上) && (点p在以线段l为对角线的矩形内) 
*******************************************************************************/   
bool online(LINESEG l,POINT p)   
{   
    return( (multiply(l.e,p,l.s)==0) &&( ( (p.x-l.s.x)*(p.x-l.e.x)<=0 )&&( (p.y-l.s.y)*(p.y-l.e.y)<=0 ) ) );   
}   
// 返回点p以点o为圆心逆时针旋转alpha(单位：弧度)后所在的位置    
POINT rotate(POINT o,float alpha,POINT p)   
{   
    POINT tp;   
    p.x-=o.x;   
    p.y-=o.y;   
    tp.x=p.x*cos(alpha)-p.y*sin(alpha)+o.x;   
    tp.y=p.y*cos(alpha)+p.x*sin(alpha)+o.y;   
    return tp;   
}   
/* 返回顶角在o点，起始边为os，终止边为oe的夹角(单位：弧度)  
    角度小于pi，返回正值  
    角

最低0.47元/天 解锁文章
华仔Ivan
关注
专栏目录
boost::polygon模块voronoi相关的测试程序
希望我的博客，能帮上你解决学习中工作中所遇到的问题
06-19 453
boost::polygon模块voronoi相关的测试程序实现功能C++实现代码 实现功能 boost::polygon模块voronoi相关的测试程序 C++实现代码 #include <cmath> #include <cstdio> #include <ctime> #include <string> #include <ieee754.h> typedef long double fpt80; #include <boost/ran
WIN7&Ubuntu配置OpenGL详解（附.h, .lib, .dll 文件）
xl928471061的博客
05-08 537
最近要用OpenGL跑人脸模型，于是开始捣鼓怎么配置openGL。但网上的各种博客，虽说讲的很详细，但给出的步骤总是不全，让我这样的小白很是苦恼，倒腾了不下十个博客才解决。所以，为了不让之后的小白苦恼，也为了自己能更熟悉，在此进行步骤详解。
polygon中心点四叉树计算测试
地理信息技术杂谈 博客，欢迎技术交流
06-13 263
https://mp.weixin.qq.com/s?__biz=MzU2ODYzNzc4OQ==&mid=2247486129&idx=1&sn=306b0827c6fbe369fcdcfa400ec8cc5a&chksm=fc8baa83cbfc23958f23100387194f13a25c217439b4198cd8b3723272b5685f0573b1bea49e&token=463482986&lang=zh_CN#rd
Google Maps API 进级： GPolygon对象2
weixin_34319999的博客
12-18 81
转自：http://hi.baidu.com/xfm_zhr/blog/item/bcfc4d8b8444e715c9fc7a2d.html 1.       GPolygon对象 多边形对象。在地图上显示一段指定颜色及位置的多边形。可用来标记指定区域。 注意：使用该对象，必须在HTML标签包含命名空间：v="urn:schemas-microsoft-com:vml"。示例：&lt;htm...
如何将选中的点集转换成Polygon
swfcsunboy
05-02 3236
本例要实现的功能是根据选中的Points创建一个Polygon，并且保存到Polygon类型的FeatureLayer中，要求被选择的Points最少为3个。l 要点根据选择的点创建一个Polygon，首先要判断生成的Polygon是否是Simple,这里用到接口ITopologicalOperator2的属性IsSimple。如果不是,则要对做Polygon排序等处理。此外还用到了接口I
ROS移植机器人小车：问题集(2)
tanmx219的博客
03-08 4042
小问题比较多，根据自己编译的情况罗列，不一定是最佳解决办法。 ------------------------------------------------------------------------- CMake Error at /usr/share/cmake-3.16/Modules/FindQt4.cmake:1314 (message): Found unsuitable Qt version "5.12.8" from /usr/bin/qmake, this code re..
判断点在多边形内的算法
windydreams的专栏
11-07 2584
最近这几天在做一个算啊，主要是判断点是否在多边形内，在
light_occur_polygon.rar_多面体 相交_多面体相交
最新发布
09-21
"light_occur_polygon.rar_多面体 相交_多面体相交"这个压缩包文件包含了相关代码和资源，帮助我们理解并实现光线与多面体相交的算法。 1. **光线与多面体相交的基本原理**： 光线在3D空间中通常表示为从一个点...
OpenGLES—API(gl2.h)
GameCocos2dx的专栏
06-10 3274
#ifndef __gl2_h_ #define __gl2_h_ /* $Revision: 10602 $ on $Date:: 2010-03-04 22:35:34 -0800 #$ */ #include #ifdef __cplusplus extern "C" { #endif /* * This document is licensed under the SGI Fr
point-in-polygon:确定点是否在多边形内
05-05
多边形点 确定点是否在多边形内部。 该模块基于从查询点投射半无限射线并计算交叉点。 如果您需要一个数值健壮的解决方案并愿意为此牺牲一些性能，请使用 。 例子 var pointInPolygon = require ( 'point-in-polygon' ) ; var polygon = [ [ 1 , 1 ] , [ 1 , 2 ] , [ 2 , 2 ] , [ 2 , 1 ] ] ; console . log ( pointInPolygon ( [ 1.5 , 1.5 ] , polygon ) ) ; // true console . log ( pointInPolygon ( [ 4.9 , 1.2 ] , polygon ) ) ; // false console . log ( pointInPolygon ( [ 1.8 , 1.1 ] , polygon
单链表的运用&运算符重载——多项式相加
qq_43432875的博客
04-04 1980
编译环境：VS2017 #poly.h #pragma once #include<iostream> using namespace std; class Node { public: int coe;//系数 int ind;//指数 Node *next; }; class Poly { private: Node *head; public: Poly();//...
编译时报警fatal:jsk_rviz_plugins/PictogramArray.h: No such file or directory
jiaostyle的专栏
01-12 995
一、问题描述： 在编译工程时，报警： fatal:jsk_rviz_plugins/PictogramArray.h: No such file or directory 二、解决： sudo apt-get install ros-kinetic-jsk-rviz-plugins
点到多边形的切线
大鹏的blog
02-28 2113
前言 暴力算法 二分搜索
【C++ 程序】 tvj::Polygon类 （多边形类）
Teddy van Jerry的博客
01-04 2526
代码 Polygon.h /* * Project: Polygon Class * File: Polygon.h * ------------------ * * @author: Teddy van Jerry * * @version: 2021/01/03 * - initial version */ /////////////////////////////////////////////////// // NOTE: // The class name Poly
多边形自相交判断issimple
reecho_cd的专栏
03-10 3692
 最近正在做多边形自相交的检查，多边形是从GDB里直接读取的。按照网上搜来的方法，首先把iskonwnsimple设为了false，再调用issimple的。但是issimple并没有返回false，而是直接返回了true。不知道是不是ITopologicalOperator2接口存在bug（我用的是engine9.2）的原因，感觉就是issimple没有真正执行一样。能够获得到的IPolygon
编译 pycaffe时报错：fatal error: numpy/arrayobject.h没有那个文件或目录
大数据超哥的博客
10-10 2114
解决方法命令行输入进pythonimport numpy as np np.get_include() 得到：/usr/local/lib/python2.7/dist-packages/numpy/core/include在Makefile.config找到PYTHON_INCLUDE，发现有点不同：PYTHON_INCLUDE := /usr/include/python2.7 \
几何算法大集合
jazywoo_在路上
08-01 3816
#ifndef _H_GEOCOMP #define _H_GEOCOMP /* 计算几何 目录 ㈠ 点的基本运算 1. 平面上两点之间距离 1 2. 判断两点是否重合 1 3. 矢量叉乘 1 4. 矢量点乘 2 5. 判断点是否在线段上 2 6. 求一点饶某点旋转后的坐标 2 7. 求矢量夹角 2 ㈡ 线段及直线的基本运算 1. 点与线段的关系 3 2. 求点到线段所在直线垂线的垂足
polygon函数用法
热门推荐
zhouqian88423的专栏
12-06 2万+
polygon  添加摘要目录 [隐藏] 1 函数功能 2 函数原型 3 主要参数 polygon-函数功能 该函数画一个由直线相闻的两个以上顶点组成的多边形，用当前画笔画多边形轮廓，用当前画刷和多边形填充模式填充多边形。polygon-函数原型 BOOL Polygon(HDC hdc, C
面向对象设计：多边形与点的数据绑定机制
"该资源是一份关于期末考核题的说明文档，主要涵盖了两个知识点：面向对象设计中的多边形Polygon类实现以及Point类的数据绑定机制。文档提到了具体的设计要求和实现方法，包括构造函数、拷贝构造函数、复制操作符、...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值