计算一个线段与一个矩形的相交面积

最新推荐文章于 2022-04-12 10:15:35 发布
最新推荐文章于 2022-04-12 10:15:35 发布
阅读量405 收藏 1
点赞数 1
分类专栏： C++ 算法 文章标签： 算法 c++
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/xuyin1204/article/details/106335957
版权
C++ 同时被 2 个专栏收录
14 篇文章 4 订阅
订阅专栏
算法
10 篇文章 0 订阅
订阅专栏 
#include <iostream>
#include <cmath>
#include <cstring>
#include <algorithm>
#include <vector>
// #include <cstdlib>
#include "stdlib.h"

using namespace std;
const double PI = 3.141592654;

const int maxn = 300;
const double eps = 1e-6;
struct Point
{
    double x, y;
};
Point p1[maxn], p2[maxn];
int n1, n2;

// 计算两个线段的角度
double normals(double x1, double y1, double x2, double y2)
{
	double n = x1 * x2 + y1 * y2;
	double m = sqrt(x1 * x1 + y1 * y1) * sqrt(x2 * x2 + y2 * y2);
	if (m == 0) {
		cout << "参数错误：";
		return -1;
	}
	return acos(n / m) * (180 / PI);
}

// 计算一个点（x2, y2）在 矢量线段[(x0, y0), (x1, y1)]的哪侧，即在线段的顺时针还是逆时针位置，还是线段
double pointLocation(double x0, double y0, double x1, double y1, double x2, double y2)
{
	return (x0 - x2) * (y1 - y2) - (y0 - y2) * (x1 - x2);
}

// 计算一个点在一个矢量图的那侧，正数在顺时针位置，负数为逆时针位置，零为在线段上或者延长线上
double pointLocation2(Point p0, Point p1, Point p2)
{
	return (p0.x - p2.x) * (p1.y - p2.y) - (p0.y - p2.y) * (p1.x - p2.x);
}

Point intersection(Point a,Point b,Point c,Point d)
{
    Point p = a;
    double t =((a.x-c.x)*(c.y-d.y)-(a.y-c.y)*(c.x-d.x))/((a.x-b.x)*(c.y-d.y)-(a.y-b.y)*(c.x-d.x));
    p.x +=(b.x-a.x)*t;
    p.y +=(b.y-a.y)*t;
    cout << "intersection p.x=" << p.x << ", p.y=" << p.y << endl;
    return p;
}

/*
Point getCrossPoint(Point a,Point b, Point a1,Point b1)
{
	Point base = b1 - b;
    double d1= abs(cross(base,a-b));
    double  d2= abs(cross(base,a1-b));
    double t=d1/(d1+d2);
    Point temp=(a1-a)*t;
    return a+temp;
}
*/


// Calculate the cross product
float cross(Point a, Point b, Point c)
{
    return (a.x - c.x) * (b.y - c.y) - (b.x - c.x) * (a.y - c.y);
}

//Calculate the intersection of two line segments
Point segmentIntersection(Point a, Point b, Point c, Point d)
{
	Point temp={-1, -1};
    if ( min(a.x, b.x) <= max(c.x, d.x) && max(a.x, b.x) >= min(c.x, d.x)
            && min(a.y, b.y) <= max(c.y, d.y) && max(a.y, b.y) >= min(c.y, d.y)
            && cross(d, a, c) * cross(b, d, c) >= 0 && cross(b, c, a) * cross(d, b, a) >= 0) {
        if (cross(d, a, c) * cross(b, d, c) == 0 && cross(b, c, a) * cross(d, b, a) == 0) {
            cout << "The two lines segment are parallel." << endl;
            return temp;
        }
        Point p = a;
        float t = ((a.x - c.x) * (c.y - d.y) - (a.y - c.y) * (c.x - d.x)) / ((a.x - b.x) * (c.y - d.y)-(a.y - b.y) * (c.x - d.x));
        p.x += (b.x - a.x) * t;
        p.y += (b.y - a.y) * t;
        cout << "intersection p.x= " << p.x << ", p.y= " << p.y << endl;
        return p;
    }
    cout << "not intersection."<< endl;
    return temp;
}

/*
int dcmp(double x)
{
    if(x > eps) return 1;
    return x < -eps ? -1 : 0;
}

//计算多边形面积
double PolygonArea(Point p[], int n)
{
    if(n < 3) return 0.0;
    double s = p[0].y * (p[n - 1].x - p[1].x);
    for(int i = 1; i < n - 1; ++ i) {
        s += p[i].y * (p[i - 1].x - p[i + 1].x);
        // cout << "p[i-1].x =" << p[i-1].x << ", p[i-1].y=" << p[i-1].y << endl;
        // cout << "p[i].x =" << p[i].x << ", p[i].y=" << p[i].y << endl;
        // cout << "p[i+1].x =" << p[i+1].x << ", p[i+1].y=" << p[i+1].y << endl;
    }
    s += p[n - 1].y * (p[n - 2].x - p[0].x);
    cout << "s =" << s << endl;
    return fabs(s * 0.5);
}
double CPIA(Point a[], Point b[], int na, int nb) //ConvexPolygonIntersectArea
{
    Point p[20], tmp[20];
    int tn, sflag, eflag;
    memcpy(p,b,sizeof(Point)*(nb));
    for(int i = 0; i < na && nb > 2; i++)
    {
    	if (i == na - 1) {
    		sflag = dcmp(cross(a[0], p[0],a[i]));
    	} else {
    		sflag = dcmp(cross(a[i + 1], p[0],a[i]));
    	}
        for(int j = tn = 0; j < nb; j++, sflag = eflag)
        {
            if(sflag>=0) {
            	tmp[tn++] = p[j];
            }
            if (i == na - 1) {
            	if (j == nb -1) {
            		eflag = dcmp(cross(a[0], p[0], a[i]));
				} else {
					eflag = dcmp(cross(a[0], p[j + 1], a[i]));
				}
			} else {
				if (j == nb -1) {
					eflag = dcmp(cross(a[i + 1], p[0], a[i]));
				} else {
					eflag = dcmp(cross(a[i + 1], p[j + 1], a[i]));
				}
			}
            if((sflag ^ eflag) == -2){
            	if (i == na - 1) {
            		if (j == nb -1) {
            			tmp[tn++] = intersection(a[i], a[0], p[j], p[0]); //求交点
            		} else {
            			tmp[tn++] = intersection(a[i], a[0], p[j], p[j + 1]);
            		}
				} else {
					if (j == nb -1) {
						tmp[tn++] = intersection(a[i], a[i + 1], p[j], p[0]);
					} else {
						tmp[tn++] = intersection(a[i], a[i + 1], p[j], p[j + 1]);
					}
				}
            }
        }
        memcpy(p, tmp, sizeof(Point) * tn);
        nb = tn, p[nb] = p[0];
    }
    if(nb < 3) return 0.0;
    return PolygonArea(p, nb);
}
double SPIA(Point a[], Point b[], int na, int nb) //SimplePolygonIntersectArea 调用此函数
{
    int i, j;
    Point t1[na], t2[nb];
    double res = 0, num1, num2;
    t1[0] = a[0], t2[0] = b[0];
    for(i = 2; i < na; i++)
    {
        t1[1] = a[i-1], t1[2] = a[i];
        num1 = dcmp(cross(t1[1], t1[2],t1[0]));
        if(num1 < 0) swap(t1[1], t1[2]);
        for(j = 2; j < nb; j++)
        {
            t2[1] = b[j - 1], t2[2] = b[j];
            num2 = dcmp(cross(t2[1], t2[2],t2[0]));
            if(num2 < 0) swap(t2[1], t2[2]);
            res += CPIA(t1, t2, 3, 3) * num1 * num2;
        }
    }
    cout << "Sum::res=" <<res << endl;
    return res;
}
*/

int main() {
	//double angle;
	//double x0, y0, x1, y1, x2, y2;
	//cout << "请输入三个点的坐标值："<<endl;
	//cin >> x0 >> y0 >> x1 >> y1 >> x2 >> y2;
	//angle = normals(x1-x0, y1 - y0, x2 - x0, y2 - y0);
	//cout << "夹角为 ：" << angle << endl;
	//double s = pointLocation(x0, y0, x1, y1, x2, y2);
	// cout << "s = " << s << endl;
	/*
	while(cin>>n1)
	{
		for(int i = 0; i < n1; i++) scanf("%lf%lf", &p1[i].x, &p1[i].y);
		double A1 = pointLocation2(p1[0], p1[1], p1[2]);
		cout << "A1 =" << A1 << endl;
	}
	*/
	Point temp0= {-1, -1};
	Point temp1= {-1, -1};
	Point temp2= {-1, -1};
	Point temp3= {-1, -1};
	Point cenPoint= {-1, -1};
	vector<Point> pointsVec;
	Point tar_points_center[4];

	while(cin>>n1>>n2)
	{
		for(int i = 0; i < n1; i++) scanf("%lf%lf", &p1[i].x, &p1[i].y);
		for(int i = 0; i < n2; i++) scanf("%lf%lf", &p2[i].x, &p2[i].y);
		// double A1 = pointLocation2(p1[0], p1[1], p1[2]);
		// cout << "A1 =" << A1 << endl;

		tar_points_center[0].x = (p2[0].x + p2[1].x) / 2.0;
		tar_points_center[0].y = (p2[0].y + p2[1].y) / 2.0;
		tar_points_center[1].x = (p2[1].x + p2[2].x) / 2.0;
		tar_points_center[1].y = (p2[1].y + p2[2].y) / 2.0;
		tar_points_center[2].x = (p2[2].x + p2[3].x) / 2.0;
		tar_points_center[2].y = (p2[2].y + p2[3].y) / 2.0;
		tar_points_center[3].x = (p2[3].x + p2[0].x) / 2.0;
		tar_points_center[3].y = (p2[3].y + p2[0].y) / 2.0;
		/*
		temp0 = segmentIntersection(p1[0], p1[1], p2[0], p2[1]);
		cout << "temp0.x = " << temp0.x << ", temp0.y=" << temp0.y << endl;
		temp1 = segmentIntersection(p1[0], p1[1], p2[1], p2[2]);
		cout << "temp1.x = " << temp1.x << ", temp1.y=" << temp1.y << endl;
		temp2 = segmentIntersection(p1[0], p1[1], p2[2], p2[3]);
		cout << "temp2.x = " << temp2.x << ", temp2.y=" << temp2.y << endl;
		temp3 = segmentIntersection(p1[0], p1[1], p2[3], p2[0]);
		cout << "temp3.x = " << temp3.x << ", temp3.y=" << temp3.y << endl;
		cout << endl;

		if (temp0.x != -1 && temp0.y != -1) {
			pointsVec.push_back(temp0);
		} else if (pointLocation2(p1[0], p1[1], tar_points_center[0]) > 0){
			pointsVec.push_back(p2[0]);
						pointsVec.push_back(p2[1]);
		}
		for (unsigned int i = 0; i < pointsVec.size(); i++) {
			cout << "temp0: x = " << pointsVec[i].x << ", y=" << pointsVec[i].y << endl;
		}
		cout << endl;

		if (temp1.x != -1 && temp1.y != -1) {
			pointsVec.push_back(temp1);
		} else if(pointLocation2(p1[0], p1[1], tar_points_center[1]) > 0) {
			pointsVec.push_back(p2[1]);
			pointsVec.push_back(p2[2]);
		}
		for (unsigned int i = 0; i < pointsVec.size(); i++) {
			cout << "temp1: x = " << pointsVec[i].x << ", y=" << pointsVec[i].y << endl;
		}
		cout << endl;

		if (temp2.x != -1 && temp2.y != -1) {
			pointsVec.push_back(temp2);
		} else if(pointLocation2(p1[0], p1[1], tar_points_center[2]) > 0) {
			pointsVec.push_back(p2[2]);
			pointsVec.push_back(p2[3]);
		}
		for (unsigned int i = 0; i < pointsVec.size(); i++) {
			cout << "temp2: x = " << pointsVec[i].x << ", y=" << pointsVec[i].y << endl;
		}
		cout << endl;

		if (temp3.x != -1 && temp3.y != -1) {
			pointsVec.push_back(temp3);
		} else if(pointLocation2(p1[0], p1[1], tar_points_center[3]) > 0) {
			pointsVec.push_back(p2[3]);
			pointsVec.push_back(p2[0]);
		}
		for (unsigned int i = 0; i < pointsVec.size(); i++) {
			cout << "temp3: x = " << pointsVec[i].x << ", y=" << pointsVec[i].y << endl;
		}
		cout << endl;


		for (unsigned int i = 0; i < pointsVec.size(); i++) {
			cout << "x = " << pointsVec[i].x << ", y=" << pointsVec[i].y << endl;
		}

		if (pointsVec.size() >= 2) {
			for (unsigned int i = 0; i < pointsVec.size() - 1; i++) {
				for (unsigned int j = i + 1; j < pointsVec.size(); j++)
					 if(pointsVec[j].x == pointsVec[i].x && pointsVec[j].y == pointsVec[i].y)
						 pointsVec.erase(pointsVec.begin() + j);
			}
		}
		cout << "after erase ..." << endl;
		for (unsigned int i = 0; i < pointsVec.size(); i++) {
			cout << "x = " << pointsVec[i].x << ", y=" << pointsVec[i].y << endl;
		}
		int flag = 0;
		cenPoint = intersection(p2[0], p2[2], p2[1], p2[3]);
		int loc = pointLocation2(p1[0], p1[1], cenPoint);
		if (0 == pointsVec.size()) {
			if (loc > 0) {
				flag = 1;
			} else if (loc < 0) {
				flag = -1;
			}
		} else {

		}
		cout << "flag = " << flag << endl;
		*/

		temp0 = intersection(p1[0], p1[1], p2[0], p2[1]);
		cout << "temp0.x = " << temp0.x << ", temp0.y=" << temp0.y << endl;
		temp1 = intersection(p1[0], p1[1], p2[1], p2[2]);
		cout << "temp1.x = " << temp1.x << ", temp1.y=" << temp1.y << endl;
		temp2 = intersection(p1[0], p1[1], p2[2], p2[3]);
		cout << "temp2.x = " << temp2.x << ", temp2.y=" << temp2.y << endl;
		temp3 = intersection(p1[0], p1[1], p2[3], p2[0]);
		cout << "temp3.x = " << temp3.x << ", temp3.y=" << temp3.y << endl;

	}
 	return 0;
}
xuyin1204
关注
专栏目录
POJ 1410 || Intersection(线段矩形相交
野生
04-01 603
判断线段与矩形相交有两种相交方式： 1）线段与矩形的四条边的其中一条至少有一个交点（X积判断 2）线段再矩形内（用坐标大小判断 十分钟一A一计算几何水题，感觉自己越来越棒了。 513KB,0MS #include #include #include using namespace std; const double eps = 1e-8; int dcmp( double x ) {
计算几何——Intersection(线段与矩形相交)
mumei的博客
07-14 4301
题目链接http://poj.org/problem?id=1410 You are to write a program that has to decide whether a given line segment intersects a given rectangle. An example: line: start point: (4,9) end point: (11,2) rect...
C++计算俩矩形的重叠面积
10-15
C++根据计算几何算法，充分考虑矩形重叠的各种情况，可以准确计算矩形的重叠率
线段 矩形 矩形与矩形 相交
luyuncsd123的专栏
03-30 2117
[csharp] view plaincopyprint? /*     功能：判断线段和矩形是否相交　     先判断线段的俩个端点是否在矩形的内部，在就必然相交           其次判断线段的包围盒是否和矩形相交，不相交的话线段和矩形肯定也不相交              最后判断，矩形的四个顶点是否位于线段的两侧，是则必然相交，否则就不相交 
求线段交点
dong2022214184的博客
09-08 431
德布是一个很大的可乐饮酒者，但他无法帮助它。 为了少喝可乐，他设计了以下装置： 一种长度和高度几乎无限的水箱，底部有无限量的可乐。 在水箱中间有两个木板，木板位于距离水箱边缘无限远的地方，木板和可乐之间的距离是无限的，两个木板可能相交， Debu将在中午将罐子倒置（稍后完成），可乐将垂直下降。可乐的速度很快，可乐到达底部片刻，所以一些 木板可以保留可乐（因为两块木板可能相交），而德布会喝木板保留的...
如何求两个矩形相交面积
qq_36238670的博客
08-23 6809
[转载]https://blog.csdn.net/jinzhichaoshuiping/article/details/51192108 如何求两个相交矩形的交集的面积? 就是交接的公共部分面积? class Rectangle { Point min; Point max; } Rectangle rect1, rect2; Rectangle rect; re...
python判断两个矩形是否相交_使用Python判断线段是否与矩形相交
weixin_39869959的博客
12-04 1586
直线相交的判定我们在初等几何中学过如何判断两条直线是否相交。在欧几里得平面上，两条直线要么平行，要么相交，要么重合。这是欧几里得第五公设的推论。相交的两条直线恰好有一个交点，而重合的两条直线有无数个交点。判断两条直线位置关系的代数方法是：联立两条之直线方程，无解，则两条直线平行，存在唯一解则两条直线相交，存在两个及以上的解，则两直线重合。联立直线方程线段相交的判定：分为矩形判定和叉积判定两步ste...
java判断线与矩形相交_判断任意多边形与矩形的相交(线段与矩形相交或线段与线段相交)...
weixin_32515349的博客
02-17 929
任意多边形与矩形的相交，其实就是判断多条线段是否与这个矩形相交，再简单点就是判断线段是否与矩形的每一条边相交了。那现在，我们先来看看判断一条线段与矩形的其中一条线段的相交的情况(上方水平线)：(图形中的a/b=d/c如果看不明白的，我也无语了，回去翻下几何图形的那中学课本。)图中已知的是红色的线段，就是要检测的线段，其起点就是star，重点是end，其坐标是已知的。图中的蓝色线是辅助理解的线。而图...
Python判断直线和矩形是否相交的方法
09-21
在Python编程中，判断直线和矩形是否相交是一个常见的几何问题，特别是在计算机图形学、游戏开发或物理模拟等场景中。本节将详细介绍如何利用Python实现这一功能，重点涉及二维坐标系下的数学运算。 首先，我们需要...
线段与矩形的交叉长度：如果线段在矩形外则返回0-matlab开发
05-30
cross_perimeter.m用于segment_rectange.m的辅助函数，它定义线段的交叉点和矩形的周长 zz_segment_rectange_test.m 是脚本测试，它通过使用随机段绘制它并打印段附近的长度来测试 segment_rectange.m。 segment_...
Python计算两个矩形重合面积代码实例
09-18
主要介绍了Python 实现两个矩形重合面积代码实例,文中通过示例代码介绍的非常详细，对大家的学习或者工作具有一定的参考学习价值,需要的朋友可以参考下
线段与矩形 是否相交 矩形与矩形
I&U
05-11 2668
/*    功能：判断线段和矩形是否相交　    先判断线段的俩个端点是否在矩形的内部，在就必然相交         其次判断线段的包围盒是否和矩形相交，不相交的话线段和矩形肯定也不相交            最后判断，矩形的四个顶点是否位于线段的两侧，是则必然相交，否则就不相交 */ using System; using System.Collections.Generic; using Sy
求两个矩阵相交部分面积
Hello_Rainy的博客
04-12 951
给你 二维 平面上两个 由直线构成且边与坐标轴平行/垂直 的矩形，请你计算并返回两个矩形覆盖的总面积。 每个矩形由其 左下 顶点和 右上 顶点坐标表示： 第一个矩形由其左下顶点 (ax1, ay1) 和右上顶点 (ax2, ay2) 定义。 第二个矩形由其左下顶点 (bx1, by1) 和右上顶点 (bx2, by2) 定义。 class Solution { public: int computeArea(int ax1, int ay1, int ax2, int ay2, int bx
线段和矩形相交 POJ 1410
weixin_30624825的博客
09-09 128
1 // 线段和矩形相交 POJ 1410 2 3 // #include <bits/stdc++.h> 4 #include <iostream> 5 #include <cstdio> 6 #include <cstdlib> 7 #include <algorithm> 8 #i...
算法：求矩形相交区域的面积
qq_36647492的博客
11-08 5597
题目描述：给定两个矩形，求这两个矩形相交区域的面积，矩形是通过处于对角线的两个点A(x1,y1)和点B(x2,y2)定义的，点A和点B的关系是下图中的一种： 输入：[-3，0]，[3，4]，[0，-1]，[9，2] 输出：6 说明： 分析： 首先分析一波，图形相交有什么规律 根据所能想象到的情况，如上图所示，不难发现有一个规律，就是两个矩形如果不相交，那么它俩的横纵坐标在坐标系上都是相邻的。 图3：矩形A的上边框，矩形A的下边框，矩形B的上边框，矩形B的下边框 图...
基础算法题——矩形面积相交
Zhengyangxinn的博客
03-01 3527
基础算法题——矩形面积相交 本人为一名普通二本学校自动化专业的大二学生，对编程有着少许兴趣。致力将算法写得更加通俗易懂。 做题心得 该题主要考查了判断、线段交，只要能够将矩形之间的关系弄清楚题目就会变得很简单。 矩形面积相交题目 题目分析 题目“对于每个矩形，我们给出它的一对相对顶点的坐标”，这也告诉了我们它的特殊情况会有很多，如果直接将数据记录下来然后就开始计算面积，无疑难度会很大，那我们有什...
【蓝桥杯】矩形面积交（C++解决）
学无止境
02-28 2348
题目描述 平面上有两个矩形，它们的边平行于直角坐标系的X轴或Y轴。对于每个矩形，我们给出它的一对相对顶点的坐标，请你编程算出两个矩形的交的面积。 输入 输入仅包含两行，每行描述一个矩形。 在每行中，给出矩形的一对相对顶点的坐标，每个点的坐标都用两个绝对值不超过10^7的实数表示。 输出 输出仅包含一个实数，为交的面积，保留到小数后两位。 样例输入 1 1 3 3 2 2 4 4 样例输...
uva11722——直线与矩形相交
mumei的博客
09-30 1288
题目链接：https://vjudge.net/contest/325833#problem/M Sidney以相同的概率密度在[t1,t2]到达不高山；Gandtom也以相同的概率密度在[s1,s2]到达。每个人都会停留W时间，问发生事情的概率。/捂脸 Input 第一行组数T(T < 500)。 接下来T行每行包含五个整数t1,t2,s1,s2,w(360≤t1 <t2 &l...
c++判断线段与矩形相交
最新发布
05-26
如果线段与矩形的任意一条边相交，则线段与矩形相交。 下面是一段示例代码，演示了如何判断线段与矩形相交： ```c++ #include using namespace std; struct Point { double x, y; }; struct Segment { Point ...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值