判断点在凸四边形内

#include <stdio.h>
#include <stdlib.h>

struct pt
{
	float x,y;
};

struct quat
{
	pt points[4];
};

float crossmulti2d( float x1, float y1, float x2, float y2)
{
	return x1 * y2 - x2 * y1;
}

bool inquat( quat _q, pt _pt)
{
	pt vec1, vec2;
	vec1.x = _q.points[1].x - _q.points[0].x;
	vec1.y = _q.points[1].y - _q.points[0].y;
	vec2.x = _pt.x - _q.points[0].x;
	vec2.y = _pt.y - _q.points[0].y;
	if( crossmulti2d( vec2.x, vec2.y, vec1.x, vec1.y ) < 0 )
	{
		return false;
	}
	vec1.x = _q.points[2].x - _q.points[1].x;
	vec1.y = _q.points[2].y - _q.points[1].y;
	vec2.x = _pt.x - _q.points[1].x;
	vec2.y = _pt.y - _q.points[1].y;
	if( crossmulti2d(vec2.x, vec2.y, vec1.x, vec1.y) < 0 )
	{
		return false;
	}
	vec1.x = _q.points[3].x - _q.points[2].x;
	vec1.y = _q.points[3].y - _q.points[2].y;
	vec2.x = _pt.x - _q.points[2].x;
	vec2.y = _pt.y - _q.points[2].y;
	if( crossmulti2d(vec2.x, vec2.y, vec1.x, vec1.y) < 0 )
	{
		return false;
	}
	vec1.x = _q.points[0].x - _q.points[3].x;
	vec1.y = _q.points[0].y - _q.points[3].y;
	vec2.x = _pt.x - _q.points[3].x;
	vec2.y = _pt.y - _q.points[3].y;
	if( crossmulti2d(vec2.x, vec2.y, vec1.x, vec1.y ) < 0 )
	{
		return false;
	}
	return true;
}

int main()
{
	quat shape = {
		{ {1, 1},{2, 1}, {1, 0},{0, 0} }
	};

	pt tests[] = {
		{0.0, 1.0},
		{1.0, 0.5},
		{1.5, 0},
		{1.5, 0.5},
		{1.5, 0.52},
		{1.5, 0.49}
		};
	for(int i =0; i<sizeof(tests)/sizeof(pt); ++i)
	{
		if( inquat(shape, tests[i]) == true )
		{
			printf("pass\n");
		}
		else
		{
			printf("failed\n");
		}
	}	
	return 0;
}

网站名称:判断点在凸四边形内
转载源于:https://www.cdcxhl.com/article10/jjddgo.html

成都网站建设公司_创新互联,为您提供搜索引擎优化营销型网站建设企业建站动态网站定制网站

广告

声明:本网站发布的内容(图片、视频和文字)以用户投稿、用户转载内容为主,如果涉及侵权请尽快告知,我们将会在第一时间删除。文章观点不代表本网站立场,如需处理请联系客服。电话:028-86922220;邮箱:631063699@qq.com。内容未经允许不得转载,或转载时需注明来源: 创新互联

网站建设网站维护公司