C++ 如何计算两条线段在何处相交?
我从中实现了线段交点公式。它非常有效,但我还需要知道这两条线段的交点。我如何修改代码才能做到这一点?(C++) 首先是Point类(实际上只是一个2d向量) 最后是检查交叉点的函数。但它发生在哪里C++ 如何计算两条线段在何处相交?,c++,intersection,C++,Intersection,我从中实现了线段交点公式。它非常有效,但我还需要知道这两条线段的交点。我如何修改代码才能做到这一点?(C++) 首先是Point类(实际上只是一个2d向量) 最后是检查交叉点的函数。但它发生在哪里 bool doIntersect(Point p1, Point q1, Point p2, Point q2) { // Find the four orientations needed for general and // special cases int o1 = orienta
bool doIntersect(Point p1, Point q1, Point p2, Point q2)
{
// Find the four orientations needed for general and
// special cases
int o1 = orientation(p1, q1, p2);
int o2 = orientation(p1, q1, q2);
int o3 = orientation(p2, q2, p1);
int o4 = orientation(p2, q2, q1);
// General case
if (o1 != o2 && o3 != o4)
return true;
// Special Cases
// p1, q1 and p2 are colinear and p2 lies on segment p1q1
if (o1 == 0 && onSegment(p1, p2, q1)) return true;
// p1, q1 and q2 are colinear and q2 lies on segment p1q1
if (o2 == 0 && onSegment(p1, q2, q1)) return true;
// p2, q2 and p1 are colinear and p1 lies on segment p2q2
if (o3 == 0 && onSegment(p2, p1, q2)) return true;
// p2, q2 and q1 are colinear and q1 lies on segment p2q2
if (o4 == 0 && onSegment(p2, q1, q2)) return true;
return false; // Doesn't fall in any of the above cases
}
你只需要实现其中的一个等式 在下面的代码中,我选择实现第一个等式
// your previous code
Point CalcIntersection(Point p1, Point q1, Point p2, Point q2){
struct Point pInt;
if(doIntersect(p1, q1, p2, q2)){
pInt.x= ((p1.x * q1.y - p1.y * q1.x)*(p2.x-q2.x) - (p2.x * q2.y - p2.y * q2.x)*(p1.x-q1.x)) / ((p1.x-q1.x) * (p2.y-q2.y) -(p1.y-q1.y) * (p2.x-q2.x));
pInt.y = ((p1.x * q1.y - p1.y * q1.x)*(p2.y-q2.y) - (p2.x * q2.y - p2.y * q2.x)*(p1.y-q1.y)) / ((p1.x-q1.x) * (p2.y-q2.y) -(p1.y-q1.y) * (p2.x-q2.x));
}
return pInt;
}
Rq:我认为你的x和y应该是一个浮点数,而不是一个整数。请显示你的代码。你只是从网站上复制了代码。显示你的代码不是复制的。只需要一个谷歌就可以得到答案。甚至维基百科也有一个相当可读的条目。是的,这是我的代码。正如我在文章中所说,我从网站上复制了它。我能找到的只是直线交点,而不是线段,维基上的答案对我来说有点太专业了,因为数学不是我的强项。这不是直线吗?我检查过了,稍后有一个直线段的公式,我想:pFloat=Point(x1+t(x2-x1),y1+t(y2-y1))。我要试试那个工具,看看能不能用。
// To find orientation of ordered triplet (p, q, r).
// The function returns following values
// 0 --> p, q and r are colinear
// 1 --> Clockwise
// 2 --> Counterclockwise
int orientation(Point p, Point q, Point r)
{
// See https://www.geeksforgeeks.org/orientation-3-ordered-points/
// for details of below formula.
int val = (q.y - p.y) * (r.x - q.x) -
(q.x - p.x) * (r.y - q.y);
if (val == 0) return 0; // colinear
return (val > 0)? 1: 2; // clock or counterclock wise
}
// The main function that returns true if line segment 'p1q1'
// and 'p2q2' intersect.
bool doIntersect(Point p1, Point q1, Point p2, Point q2)
{
// Find the four orientations needed for general and
// special cases
int o1 = orientation(p1, q1, p2);
int o2 = orientation(p1, q1, q2);
int o3 = orientation(p2, q2, p1);
int o4 = orientation(p2, q2, q1);
// General case
if (o1 != o2 && o3 != o4)
return true;
// Special Cases
// p1, q1 and p2 are colinear and p2 lies on segment p1q1
if (o1 == 0 && onSegment(p1, p2, q1)) return true;
// p1, q1 and q2 are colinear and q2 lies on segment p1q1
if (o2 == 0 && onSegment(p1, q2, q1)) return true;
// p2, q2 and p1 are colinear and p1 lies on segment p2q2
if (o3 == 0 && onSegment(p2, p1, q2)) return true;
// p2, q2 and q1 are colinear and q1 lies on segment p2q2
if (o4 == 0 && onSegment(p2, q1, q2)) return true;
return false; // Doesn't fall in any of the above cases
}
// your previous code
Point CalcIntersection(Point p1, Point q1, Point p2, Point q2){
struct Point pInt;
if(doIntersect(p1, q1, p2, q2)){
pInt.x= ((p1.x * q1.y - p1.y * q1.x)*(p2.x-q2.x) - (p2.x * q2.y - p2.y * q2.x)*(p1.x-q1.x)) / ((p1.x-q1.x) * (p2.y-q2.y) -(p1.y-q1.y) * (p2.x-q2.x));
pInt.y = ((p1.x * q1.y - p1.y * q1.x)*(p2.y-q2.y) - (p2.x * q2.y - p2.y * q2.x)*(p1.y-q1.y)) / ((p1.x-q1.x) * (p2.y-q2.y) -(p1.y-q1.y) * (p2.x-q2.x));
}
return pInt;
}