VC++ 平面坐标中点集的判断

于 2019-12-30 16:54:16 发布
阅读量663 收藏 2
点赞数
分类专栏： VC++（日积月累篇） VC++（UI设计篇）
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/u012156872/article/details/103769075
版权

之前做雷达扫射图的时候用到的一个接口类，本接口实现以下功能：

1、已知点P(x,y)和多边形Poly，判断点P(x,y)是否在多边形内部。
bool IsPointInPolygon(std::vector<POINT> poly,POINT pt);

2、已知两条线段P1P2和Q1Q2，判断P1P2和Q1Q2是否相交，若相交，求出交点
bool GetCrossPOINT(const POINT &p1,const POINT &p2,const POINT &q1,const POINT &q2,long &x,long &y);

3、多边形交集-已知两个多边形Poly1和Poly2，分别由点集C1={P1,P2,...,Pm}和C2={Q1,Q2,...,Qn}表示，求这两个多边形的交集。
bool PolygonClip(const vector<POINT> &poly1,const vector<POINT> &poly2, std::vector<POINT> &interPoly);

4、判断点在圆弧内
bool IsPointInAngleArc(const POINT &pCenter, const float nRadius, const int nRadian, const int nAngle, const POINT &pPoint);

头文件声明（MathFormula.h）：

#pragma once
#include <vector>
using namespace std;

#define PI	3.1415926
inline int sgn(float d)
{
	if(d<0) return -1;
	else if (d==0) return 0;
	else return 1;
}

class CMathFormula
{
public:
	CMathFormula(void);
	~CMathFormula(void);

private:
	// 多边形点集排序-若点a大于点b,即点a在点b顺时针方向,返回true,否则返回false
	bool POINTCmp(const POINT &a,const POINT &b,const POINT &center);

	// 线段相交-排斥实验
	bool IsRectCross(const POINT &p1,const POINT &p2,const POINT &q1,const POINT &q2);

	// 线段相交-跨立判断
	bool IsLineSegmentCross(const POINT &pFirst1,const POINT &pFirst2,const POINT &pSecond1,const POINT &pSecond2);
protected:
	// 多边形点集排序
	// 已知多边形点集C={P1,P2,...,PN}，其排列顺序是杂乱，依次连接这N个点，无法形成确定的多边形，需要对点集C进行排序后，再绘制多边形。
	void ClockwiseSortPOINTs(std::vector<POINT> &vPOINTs);

public:
	// 点在多边形内部
	// 已知点P(x,y)和多边形Poly，判断点P(x,y)是否在多边形内部。
	bool IsPointInPolygon(std::vector<POINT> poly,POINT pt);

	// 线段相交
	// 已知两条线段P1P2和Q1Q2，判断P1P2和Q1Q2是否相交，若相交，求出交点
	bool GetCrossPOINT(const POINT &p1,const POINT &p2,const POINT &q1,const POINT &q2,long &x,long &y);

	// 多边形交集-已知两个多边形Poly1和Poly2，分别由点集C1={P1,P2,...,Pm}和C2={Q1,Q2,...,Qn}表示，求这两个多边形的交集。
	bool PolygonClip(const vector<POINT> &poly1,const vector<POINT> &poly2, std::vector<POINT> &interPoly);

	// 判断点在圆弧内
	bool IsPointInAngleArc(const POINT &pCenter, const float nRadius, const int nRadian, const int nAngle, const POINT &pPoint);
};

 源码实现（MathFormula.cpp）：

#include "stdafx.h"
#include "MathFormula.h"

/*
	引用：http://www.cnblogs.com/dwdxdy/p/3230485.html
	问题描述：已知两条线段P1P2和Q1Q2，判断P1P2和Q1Q2是否相交，若相交，求出交点。
	两条线段的位置关系可以分为三类：有重合部分、无重合部分但有交点、无交点。

	算法的步骤如下：
	1.快速排斥实验。
	设以线段P1P2为对角线的矩形为R，设以线段Q1Q2为对角线的矩形为T，如果R和T不相交，则两线段不相交。
	
	2.跨立实验。
	如果两线段相交，则两线段必然相互跨立对方。
	若P1P2跨立Q1Q2，则矢量(P1-Q1)和(P2-Q1)位于矢量(Q2-Q1)的两侧，即( P1 - Q1 ) × ( Q2 - Q1 ) * ( P2 - Q1 ) × ( Q2 - Q1 ) < 0。
	若Q1Q2跨立P1P2，则矢量(Q1-P1)和(Q2-P1)位于矢量(P2-P1)的两侧，即( Q1 - P1 ) × ( P2 - P1 ) * ( Q2 - P1 ) × ( P2 - P1 ) < 0。
*/

//排斥实验
bool CMathFormula::IsRectCross(const POINT &p1,const POINT &p2,const POINT &q1,const POINT &q2)
{
	bool ret = min(p1.x,p2.x) <= max(q1.x,q2.x)    &&
		min(q1.x,q2.x) <= max(p1.x,p2.x) &&
		min(p1.y,p2.y) <= max(q1.y,q2.y) &&
		min(q1.y,q2.y) <= max(p1.y,p2.y);
	return ret;
}

//跨立判断
bool CMathFormula::IsLineSegmentCross(const POINT &pFirst1,const POINT &pFirst2,const POINT &pSecond1,const POINT &pSecond2)
{
	long line1,line2;
	line1 = pFirst1.x * (pSecond1.y - pFirst2.y) +
		pFirst2.x * (pFirst1.y - pSecond1.y) +
		pSecond1.x * (pFirst2.y - pFirst1.y);
	line2 = pFirst1.x * (pSecond2.y - pFirst2.y) +
		pFirst2.x * (pFirst1.y - pSecond2.y) + 
		pSecond2.x * (pFirst2.y - pFirst1.y);
	if (((line1 ^ line2) >= 0) && !(line1 == 0 && line2 == 0))
		return false;

	line1 = pSecond1.x * (pFirst1.y - pSecond2.y) +
		pSecond2.x * (pSecond1.y - pFirst1.y) +
		pFirst1.x * (pSecond2.y - pSecond1.y);
	line2 = pSecond1.x * (pFirst2.y - pSecond2.y) + 
		pSecond2.x * (pSecond1.y - pFirst2.y) +
		pFirst2.x * (pSecond2.y - pSecond1.y);
	if (((line1 ^ line2) >= 0) && !(line1 == 0 && line2 == 0))
		return false;
	return true;
}

bool CMathFormula::GetCrossPOINT(const POINT &p1,const POINT &p2,const POINT &q1,const POINT &q2,long &x,long &y)
{
	if(IsRectCross(p1,p2,q1,q2))
	{
		if (IsLineSegmentCross(p1,p2,q1,q2))
		{
			//求交点
			long tmpLeft,tmpRight;
			tmpLeft = (q2.x - q1.x) * (p1.y - p2.y) - (p2.x - p1.x) * (q1.y - q2.y);
			tmpRight = (p1.y - q1.y) * (p2.x - p1.x) * (q2.x - q1.x) + q1.x * (q2.y - q1.y) * (p2.x - p1.x) - p1.x * (p2.y - p1.y) * (q2.x - q1.x);

			x = (int)((double)tmpRight/(double)tmpLeft);

			tmpLeft = (p1.x - p2.x) * (q2.y - q1.y) - (p2.y - p1.y) * (q1.x - q2.x);
			tmpRight = p2.y * (p1.x - p2.x) * (q2.y - q1.y) + (q2.x- p2.x) * (q2.y - q1.y) * (p1.y - p2.y) - q2.y * (q1.x - q2.x) * (p2.y - p1.y); 
			y = (int)((double)tmpRight/(double)tmpLeft);
			return true;
		}
	}
	return false;
}

/*
	引用：http://www.cnblogs.com/dwdxdy/p/3230647.html
	问题描述：已知点P(x,y)和多边形Poly，判断点P(x,y)是否在多边形内部。
	基本方法：射线法
	以点P为端点，向左方作射线L，由于多边形是有界的，所以射线L的左端一定在多边形外部，考虑沿着L从无究远处开始自左向右移动。
	遇到和多边形的第一个交点的时候，进入到了多边形的内部，遇到第二个交点的时候，离开了多边形...
	因而当L和多边形的交点数目C是奇数的时候，P在多边形内，是偶数，则P在多边形外。
	特殊情况分析，如图下图(a),(b),(c),(d)所示。

	图(a)中，L和多边形的顶点相交，交点只能计算一个。
	图(b)中，L和多边形顶点的交点不应被计算。
	图(c)和(d)中，L和多边形的一条边重合，这条边应该被忽略不计。

*/
//  The function will return YES if the POINT x,y is inside the polygon, or
//  NO if it is not.  If the POINT is exactly on the edge of the polygon,
//  then the function may return YES or NO.
bool CMathFormula::IsPointInPolygon(std::vector<POINT> poly,POINT pt)
{
	int i,j;
	bool c = false;
	for (i = 0,j = poly.size() - 1;i < poly.size();j = i++)
	{
		if ((((poly[i].y <= pt.y) && (pt.y < poly[j].y)) ||
			((poly[j].y <= pt.y) && (pt.y < poly[i].y)))
			&& (pt.x < (poly[j].x - poly[i].x) * (pt.y - poly[i].y)/(poly[j].y - poly[i].y) + poly[i].x))
		{
			c = !c;
		}
	}
	return c;
}

/*
	问题描述：已知多边形点集C={P1,P2,...,PN}，其排列顺序是杂乱，依次连接这N个点，无法形成确定的多边形，需要对点集C进行排序后，再绘制多边形。
	点集排序过程中，关键在于如何定义点的大小关系。
	以按逆时针排序为例，算法步骤如下：

	定义：点A在点B的逆时针方向，则点A大于点B
	1.计算点集的重心O，以重心作为逆时针旋转的中心点。
	2.计算点之间的大小关系。

	大小关系的计算，可由两种方法进行计算。

	方法1：
	以重心O作一条平行于X轴的单位向量OX，然后依次计算OPi和OX的夹角，根据夹角的大小，确定点之间的大小关系。
	OPi和OX夹角越大，说明点Pi越小。

	方法2：
	根据向量叉积的定义，向量OPi和OPj的叉积大于0，则向量OPj在向量OPi的逆时针方向，即点Pj小于点Pi。
*/
//若点a大于点b,即点a在点b顺时针方向,返回true,否则返回false
bool CMathFormula::POINTCmp(const POINT &a,const POINT &b,const POINT &center)
{
	if (a.x >= 0 && b.x < 0)
		return true;
	if (a.x == 0 && b.x == 0)
		return a.y > b.y;
	//向量OA和向量OB的叉积
	int det = (a.x - center.x) * (b.y - center.y) - (b.x - center.x) * (a.y - center.y);
	if (det < 0)
		return true;
	if (det > 0)
		return false;
	//向量OA和向量OB共线，以距离判断大小
	int d1 = (a.x - center.x) * (a.x - center.x) + (a.y - center.y) * (a.y - center.y);
	int d2 = (b.x - center.x) * (b.x - center.y) + (b.y - center.y) * (b.y - center.y);
	return d1 > d2;
}

void CMathFormula::ClockwiseSortPOINTs(std::vector<POINT> &vPOINTs)
{
	//计算重心
	POINT center;
	double x = 0,y = 0;
	for (int i = 0;i < vPOINTs.size();i++)
	{
		x += vPOINTs[i].x;
		y += vPOINTs[i].y;
	}
	center.x = (int)x/vPOINTs.size();
	center.y = (int)y/vPOINTs.size();

	//冒泡排序
	for(int i = 0;i < vPOINTs.size() - 1;i++)
	{
		for (int j = 0;j < vPOINTs.size() - i - 1;j++)
		{
			if (POINTCmp(vPOINTs[j],vPOINTs[j+1],center))
			{
				POINT tmp = vPOINTs[j];
				vPOINTs[j] = vPOINTs[j + 1];
				vPOINTs[j + 1] = tmp;
			}
		}
	}
}

/*
	引用：http://www.cnblogs.com/dwdxdy/archive/2013/08/02/3232110.html
	问题描述：已知两个多边形Poly1和Poly2，分别由点集C1={P1,P2,...,Pm}和C2={Q1,Q2,...,Qn}表示，求这两个多边形的交集。
	算法思想：
	两个多边形相交后，其顶点要么是两个多边形边的交点，要么是在多边形内部的点。

	算法步骤：
	1.计算两个多边形每条边之间的交点。
	2.计算包含在多边形内部的点。
	3.将交点和多边形内部的点，按逆时针(或顺时针)排序，得出最终的点集。
*/
bool CMathFormula::PolygonClip(const vector<POINT> &poly1,const vector<POINT> &poly2, std::vector<POINT> &interPoly)
{
	if (poly1.size() < 3 || poly2.size() < 3)
	{
		return false;
	}

	long x,y;
	//计算多边形交点
	for (int i = 0;i < poly1.size();i++)
	{
		int poly1_next_idx = (i + 1) % poly1.size();
		for (int j = 0;j < poly2.size();j++)
		{
			int poly2_next_idx = (j + 1) % poly2.size();
			if (GetCrossPOINT(poly1[i],poly1[poly1_next_idx],
				poly2[j],poly2[poly2_next_idx],
				x,y))
			{
				POINT pt;
				pt.x = x;
				pt.y = y;
				interPoly.push_back(pt);
			}
		}
	}

	//计算多边形内部点
	for(int i = 0;i < poly1.size();i++)
	{
		if (IsPointInPolygon(poly2,poly1[i]))
		{
			interPoly.push_back(poly1[i]);
		}
	}
	for (int i = 0;i < poly2.size();i++)
	{
		if (IsPointInPolygon(poly1,poly2[i]))
		{
			interPoly.push_back(poly2[i]);
		}
	}

	if(interPoly.size() <= 0)
		return false;

	//点集排序 
	ClockwiseSortPOINTs(interPoly);
	return true;
}

// 判断点在圆弧内
/*
	pCenter-圆心
	nRadius-半径
	Radian-弧度
	pPoint-点坐标
*/

/*
			/
           / |
	    C /  | B
		 /   |
		/    |
	   ------
	      A

	△ABC中，∠C=90度，A，B，C的对边分别为a，b，c，那么有以下等式： a²+b²=c²（勾股定理）
	 sinA=cosB=a/c cosA=sinB=b/c tanA=cotB=a/b cotA=tanB=b/a sinA²+cosA²=1 tanA×cotA=1
	 tanB=b/a tanA=a/b
	 ∠A=atan(a/b) ∠B=atan(b/a)

	 函数里面的是弧度
	 弧度->角度	弧度*180/PI
	 角度->弧度	角度*PI/180
*/
bool CMathFormula::IsPointInAngleArc(const POINT &pCenter, const float nRadius, const int nRadian, const int nAngle, const POINT &pPoint)
{
	float x = pPoint.x - pCenter.x;
	float y = pCenter.y - pPoint.y;

	float a = abs((float)(x));
	float b = abs((float)(y));
	float fDistance = sqrt((float)(pow((float)a, 2)) + (float)(pow((float)b, 2)));
	if (fDistance > nRadius)
		return false;	// 距离大于半径
	
	// 求∠B大小，与fRadian进行比较
	float fV = b/a;
	float fArcB = atan(fV) * 180 / PI;
	if (0 == sgn(x))
	{
		if (0 == sgn(y))
			fArcB = 0;
		else if (sgn(y) > 0)
			fArcB = 90;
		else 
			fArcB = 270;
	}
	else if (0 == sgn(y))
	{
		if (sgn(x) < 0)
			fArcB = 180;
		else 
			fArcB = 0;
	}
	else if (sgn(x)>0 && sgn(y)>0)	// 1象限
	{
		fArcB = atan(fV) * 180 / PI;
	}
	else if (sgn(x)<0 && sgn(y)>0)	// 2象限
	{
		fArcB = 180 - atan(fV) * 180 / PI;
	}
	else if (sgn(x)<0 && sgn(y)<0)	// 3象限
	{
		fArcB = 180 + atan(fV) * 180 / PI;
	}
	else if (sgn(x)>0 && sgn(y)<0)	// 4象限
	{
		fArcB = 360 - atan(fV) * 180 / PI;
	}

	float fMinAngleArc = nAngle - nRadian/2;
	float fMaxAngleArc = nAngle + nRadian/2;
	if (fMinAngleArc>=0)
	{
		if (fArcB>fMaxAngleArc || fArcB<fMinAngleArc)
			return FALSE;
	}
	else
	{
		if (fArcB>fMaxAngleArc && fArcB<(360+fMinAngleArc))
			return FALSE;
	}
	return true;
}

CMathFormula::CMathFormula(void)
{
}

CMathFormula::~CMathFormula(void)
{
}

 

老狼主
关注
  • 0
    点赞
  • 踩
  • 2
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
博客
TightVNC H264编解码（一）
12-31 1597
经过尽二个月的不断摸索，TightVNC终于支持H264编解码了，前期真正编写H264编解码器只花了一周左右时间，但是测试发现效果并不是太理想（帧率太低，延迟大），这大概是所有远程控制软件要处理的通病吧。作为一位经验丰富的程序员，造成这种问题的原因，有以下几点：1、抓屏？2、H264编码？3、TCP网络发送？4、H264解码？5、视频渲染？为了排除到底哪个环节出了问题，看来不摸透底层事务流程是行不通的，不得已摸透代码画出了下面这个简陋的流程图。Server端事务处理流程图.
博客
QT HttpServer
07-31 862
资源从网上收集，包含httpserver、sslserver、http-parser三个模块，解决了编译过程中由于路径引用错误问题。下载链接：https://download.csdn.net/download/u012156872/20664799，下载编译即可应用到自己项目中，简单快速，不用造车。...
博客
QT 界面设计篇（水波纹进度条QProgressBarWater）
11-30 1688
效果图：头文件声明（QProgressBarWater.h）：#ifndef QPROGRESSBARWATER_H#define QPROGRESSBARWATER_H#include <QWidget>#include <QPainter>#include <QPaintEvent>class QProgressBarWater : public QWidget{ Q_OBJECTpublic: explic.
博客
EasyPR编译指南
08-02 3227
1、下载源码https://github.com/liuruoze/EasyPRDownloadZIP或者git（https://github.com/liuruoze/EasyPR.git）将代码解压到硬盘目录下。例如F:\车牌识别\EasyPR（车牌识别）\EasyPR2、双击解决方案文件EasyPR.sln，用VSIDE（这里使用VS2019）打开，提示点击确定。3、选择libeasypr项目，编译出现如下错误这是因为没有配置OpenCV环境所致。4、...
博客
C++ Json中文乱码问题
07-14 242
但是不知道从哪个版本开始，Json内部会自动将中文变为Unicode编码格式，也就是\u打头的字符串，这个字符串其实就是将中文转换为了Unicode格式，如果再次将strJson转换为utf-8编码发送给对方，这个时候对方接收显示的时候肯定就是乱码了。改为以后，之前cJson赋值咋样，我们得到的strJson就是咋样！
博客
C++ linux下FFmpeg转发摄像头视频流(续)
06-19 168
本篇文章介绍使用v4l2接口捕获摄像头图像数据进行ffmpeg转发。
博客
C++ linux下FFmpeg转发摄像头视频流
06-19 267
在liunx下，捕获摄像头视频流有2种方案，一种是使用ffmpeg捕获流，另一种是使用v4l2接口捕获视频流。前者使用方法简单，但是有后遗症，并不是所有型号摄像头都能够获取成功，这就是我写这篇文章的目的。老规矩，直接粘贴代码(这里虽然是使用Qt代码格式，稍加变换，即可支持其他平台)。源文件申明（PushVideoThread.cpp）
博客
QT MySQL驱动
05-14 154
在某些场合QT自带的QSQLDatabase驱动有效，但是open的时候，程序直接崩溃了，原因未知。不得已将上次使用过的。根据QT接口进行改动，进行无缝衔接。
博客
QT GDAL加载显示tif图像
05-11 365
2、加载tiff图像文件。在此之前记得引用头文件。3、GDAL库资源释放。1、GDAL库初始化。
博客
QT vlc音视频开发踩坑
05-11 201
为此我检查错误耗时将近1整天时间，甚至于都一度怀疑人生（之前在VC++下开发一切顺利）。SDK初始化一切顺利，最终在调用如下代码的时候始终失败，返回NULL值。问题解决，原来如此简单啊！
博客
QT 在线安装失败
05-07 287
QT在线安装使用的国外服务器，所以导致安装极其缓慢或者失败，唯一能想到的解决方法就是切换到国内的镜像源（除非翻墙）。
博客
QT 无标题栏窗口的拖动和拉伸
12-14 523
如下2种方案，2种方案都不完美（拉伸窗口都会抖动），但是如果一定要使用，强烈推荐第1种方案（可减小拉伸时窗口抖动）
博客
C++ libVNC开源库
12-11 654
VNC具体用途，不用过多介绍了。如果需要实时更新服务器端画面，while循环就可以了。
博客
QT QCefView自编译
12-09 598
6、回到QCefView目录，执行批处理generate-win-x86_64.bat即可，出现如下信息表示生成成功。5、设置QT版本，修改配置文件QCefView/cmake/QtConfig.cmake，如下。4、请根据对应的版本修改编译文件CefViewCore/CMakeLists.txt，如下。2、下载QCefViewCore，放到QCefViewCore目录下，地址。请注意最新版本兼容性可能存在问题，我下载的版本是。1、下载QCefView，地址。3、下载Cef库文件，地址。
博客
Ubuntu 配置共享文件夹
10-26 286
安装虚拟机经常会遇到已经共享了文件夹，但是就是找不到目录文件，搜集网上资料，汇总一下指令，留作以后备用吧：vmware-hgfsclient在文件最后一行添加如下内容：
博客
C++ OleDB数据库驱动（MySql驱动姊妹篇）
10-11 251
这里只介绍OelDB调用数据的方法。OleDB因为兼容ODBC和ADO，所以受大众追捧。大家也都非常熟悉了，这里就直接上代码吧。
博客
Visual Assist X插件安装失败总结
09-27 1149
2、安装成功以后，却显示插件无法启用，或者是插件启动了，但是实则没有启动，看不到插件选项。1、如果是运行VA_X_SETUP.exe无反应，无法启动，请右键属性，修改兼容模式。3、此时重启应该就正常了。
博客
常用时间日期格式符号定义
09-27 2759
编程过程中经常会用到时间日期格式化，这里是一些常见的日期时间格式字符定义的列表：年份：yyyyyy月份：MMMMMM天数：ddd小时：hhh分钟：mmm秒数：sss毫秒：zzzz星期几：dddddddAP时区：t。
博客
C++ linux下使用X11实现屏幕截图（续）
09-18 315
使用libpng库将XGetImage获取到的内容保存为PNG格式的文件。
博客
C++ MFC中如何根据不同文件编码格式读取和显示文件内容
09-17 837
这里假设VC++的字符集编码格式是Unicode，要显示的文本内容是strResult。上一篇，介绍了写文件，这一篇我要介绍的如何读取文件内容并显示出来。3、读取unicode-le编码格式的文本文件。4、读取unicode-be编码格式的文本文件。2、读取utf-8编码格式的文本文件。1、读取ANSI编码格式的文本文件。
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值