判断线段相交程序

    最近写程序用到，先上一段找来的别人的程序：

/*
如何判断两条线段是否相交？

也许看到这个题目，大多数人会想到利用计算直线是否相交的方法，即通过解直线方程，来求得两直线的交点（如果存在的话），然后再通过判断该交点是否在两条线段上的方法来解决。嗯，这样是可以，但是有没有更加快捷的方法呢？其实利用向量叉乘的性质（乘积的正负性可以表示两个向量之间的相对转动方向，如a×b为正，则b通过顺时针旋转可以与a共线，否则通过逆时针旋转可以与a共线）就可以很方便的解决这个问题。 

快速排斥实验：包含两条线段的最小矩形是否相交，不相交则肯定不行，这样可省掉一些计算。 
跨立试验：如果两线段相交，则两线段必然相互跨立对方。若a1b1跨立a2b2，矢量(a1-a2)和(b1-a2)位于矢量(a2-b2) 的两侧，即   
(a1-a2)×(a1-b1)*(a1-b1)×(a1-b2)>=0 同理,a2b2跨立a1b1   
(a2-a1)×(a2-b2)*(a2-b2)×(a2-b1)>=0
*/

typedef struct SDL_Point {
  Sint16 x;
  Sint16 y;
} SDL_Point;


int max(int x, int y)

{

    //比较两个数的大小，返回大的数

    if(x > y) return x;

    else return y; 

}

 

int min(int x, int y)

{

    //比较两个数的大小，返回小的数

    if(x < y) return x;

    else return y; 

}


int multi(SDL_Point p1, SDL_Point p2, SDL_Point p0)

{

    //求矢量[p0, p1], [p0, p2]的叉积

    //p0是顶点

    return (p1.x - p0.x) * (p2.y - p0.y) - (p2.x - p0.x) * (p1.y - p0.y);

    //若结果等于0，则这三点共线

    //若结果大于0，则p0p2在p0p1的逆时针方向

    //若结果小于0，则p0p2在p0p1的顺时针方向

}

 

bool isIntersected(SDL_Point s1, SDL_Point e1, SDL_Point s2, SDL_Point e2)

{

    //判断线段是否相交

    //1.快速排斥试验判断以两条线段为对角线的两个矩形是否相交

    //2.跨立试验

    if(

              (max(s1.x, e1.x) > min(s2.x, e2.x)) &&

              (max(s2.x, e2.x) > min(s1.x, e1.x)) &&

              (max(s1.y, e1.y) > min(s2.y, e2.y)) &&

              (max(s2.y, e2.y) > min(s1.y, e1.y)) &&

              (multi(s2, e1, s1) * multi(e1, e2, s1) >= 0) &&

              (multi(s1, e2, s2) * multi(e2, e1, s2) >= 0)

              )  return true;

    

    return false;    

}




              







    大致思路是分两步：



1. 一条线段大的坐标要比别一条小的坐标大



2. 一条线段的两个端点在另一线段的两边（通过比较叉乘得到）







    按照段代码，写到程序里，一般运行没什么问题，但总是有时候会有些问题，相交判断不准确。于是开始怀疑这个算法的正确性，经手算验证后发现这段代码没问题。经过反复测试，才发问题确实出在这段代码上。最后才发现问题所在：溢出。



    对于一个16位表示的int，其表示范围在–32,768 到 32,767。在我的实际程序中


的坐标范围是800*480，坐标略大的话，Multi之后的结果就超 出32767。



    于是改用32位长的int，这样能够正确表示Multi的结果，但是对于下面这段代码：



              (multi(s2, e1, s1) * multi(e1, e2, s1) >= 0) && (multi(s1, e2, s2) * multi(e2, e1, s2) >= 0)



    multi的结果相乘后仍将溢出，可以改为逻辑判断解决。修改后使用如下的代码：



int MAX(a, b)
{
	if (a > b)
	{
		return a;
	}
	else
	{
		return b;
	}
}

int MIN(a, b)
{
	if (a < b)
	{
		return a;
	}
	else
	{
		return b;
	}
}

int MULTI(int x1, int y1, int x2, int y2, int x0, int y0)
{
	return ((x1 - x0) * (y2 - y0) - (x2 - x0) * (y1 - y0));
}

int isIntersected(int x1, int y1, int x2, int y2,
				  int x3, int y3, int x4, int y4)
{
	int t1, t2;

	if ((MAX(x1, x2) < MIN(x3, x4)) || (MAX(x3, x4) < MIN(x1, x2)) ||
		(MAX(y1, y2) < MIN(y3, y4)) || (MAX(y3, y4) < MIN(y1, y2)))
	{
		return 0;
	}

	t1 = MULTI(x3, y3, x2, y2, x1, y1);
	t2 = MULTI(x2, y2, x4, y4, x1, y1);

	if (t1 >= 0 && t2 < 0 || t1 < 0 && t2 >= 0)
	{
		return 0;
	}

	t1 = MULTI(x1, y1, x4, y4, x3, y3);
	t2 = MULTI(x4, y4, x2, y2, x3, y3);

	if (t1 >= 0 && t2 < 0 || t1 < 0 && t2 >= 0)
	{
		return 0;
	}
	printf("%d, %d, %d, %d, %d, %d, %d, %d, ",
				  x1, y1, x2, y2,
				  x3, y3, x4, y4);
	printf("isInter? %d, %d, %d, %d, %d, %d, %d, %d, %d, %d/n",
	(MAX(x1, x2) > MIN(x3, x4)),
		(MAX(x3, x4) > MIN(x1, x2)),
		(MAX(y1, y2) > MIN(y3, y4)),
		(MAX(y3, y4) > MIN(y1, y2)),
		MULTI(x3, y3, x2, y2, x1, y1) * MULTI(x2, y2, x4, y4, x1, y1),
		MULTI(x1, y1, x4, y4, x3, y3) * MULTI(x4, y4, x2, y2, x3, y3),
		MULTI(x3, y3, x2, y2, x1, y1) , MULTI(x2, y2, x4, y4, x1, y1),
		MULTI(x1, y1, x4, y4, x3, y3) , MULTI(x4, y4, x2, y2, x3, y3));
	return 1;
}




 



    终于，运行正常~