在我的软件中，我尝试计算两条线的交点。我从第一行的四个双变量（x 的两个/y 的两个）开始，第二行的四个双变量开始。数学对我来说很清楚，而且在纸面上它是有效的。为了进一步处理，我必须将交叉点保存在 CSV 文件中。但是当我打开文件时，我所拥有的只是-1.#IND，据我所知，这意味着negative infinite或Not A Number。所以这是我的代码：
double A_x=boundingbox[i][0];
double A_y=boundingbox[i][1];
double B_x=boundingbox[i+1][0];
double B_y=boundingbox[i+1][0];
double C_x = matrix[j][0];
double C_y = matrix[j][1];
double D_x = matrix[j+1][0];
double D_y = matrix[j+1][1];
所以首先我从两个向量中得到四个点的坐标。
double AB_x=B_x - A_x;
double AB_y=B_y - A_y;
double m1 = AB_y / AB_x;
double c1 = A_y - m1*A_x;
double CD_x = D_x-C_x;
double CD_y = D_y-C_y;
double m2 = CD_y/CD_x;
double c2 = C_y-m1*C_x;
有了坐标，我就可以为两条线建立方程。在下一步中，我将查看变量m1和m2. 如果m1减号m2为 0，则没有交点。
if( (m1-m2) != 0) {
double point_x = (c2 - c1) / (m1 - m2);
double point_y = m1 * point_x + c1;
intersects.push_back(std::array<double, 3>());
int q = intersects.size()-1;
intersects[q][0]=point_x;
intersects[q][1]=point_y;
intersects[q][2]=i;
}
最后，我只使用for-loop 将矢量相交保存在 CSV 文件中。
for(int i=0; i<intersects.size(); i++)
{
file << "Intersects" << ";" << intersects[i][0] << ";" << intersects[i][1] << ";" << intersects[i][2] << endl;
}
但是当我打开文件时，我刚刚提到了错误。我已经尝试过另一种数据类型，例如 float，但没有任何改变。我知道代码的第一部分和最后一部分是有效的，因为我已经将它与不同类型的数学一起使用，但是我的数学方法被错误地选择了，所以我不得不改变它。所以我认为问题必须在if-statement 的条件下。
编辑：在评论的帮助下，我现在修改了我的代码。所以这是我的新代码：
for(int a=fillingStart; a<boundingbox.size()-1; a += 2)
{
double Ax=boundingbox[a][0];
double Ay=boundingbox[a][1];
double Bx=boundingbox[a+1][0];
double By=boundingbox[a+1][1];
double B1=Bx-Ax;
double A1=By-Ay;
for(size_t j=0; j<(matrix.size()-1); j++)
{
double Cx=matrix[j][0];
double Cy=matrix[j][1];
double Dx=matrix[j+1][0];
double Dy=matrix[j+1][1];
double B2=Dx-Cx;
double A2=Dy-Cy;
double det=A1*B2-A2*B1;
if(det!=0) { // det == 0 -> Lines are parallel
double C1=A1*Ax+B1*Ay;
double C2=A2*Cx + B2*Cy;
double point_x=(B2*C1-B1*C2)/det;
double point_y=(A1*C2-A2*C1)/det;
if(IsPointInBoundingBox(Ax, Ay, Bx, By, point_x, point_y) == 1 && IsPointInBoundingBox(Cx, Cy, Dx, Dy, point_x, point_y) == 1)
{
intersects.push_back(std::array<double, 3>());
int q = intersects.size()-1;
intersects[q][0]=point_x;
intersects[q][1]=point_y;
intersects[q][2]=matrix[j][2];
}
}
}
}
int IsPointInBoundingBox(double x1, double y1, double x2, double y2, double px, double py)
{
if(min(x1, x2) <= px && px <= max(x1,x2) && min(y1,y2) <= py && py <= max(y1,y2) )
{
return 1;
}
else
return 0;
}
当我运行程序时，文件中没有任何交集。获得文件交叉点的唯一方法是修改IsPointInBoundingBox函数。为了让它工作，我必须使用这个代码：
if( (px+0.01) >= min(x1,x2) && (px-0.01) <= max(x1,x2) &&
(py+0.01) >= min(y1,y2) && (py-0.01) <= max(y1,y2) )
{
return 1;
}
else
return 0;
我的大问题是现在我必须使用不同的输入。通过添加+-0.1到我的代码中，我将不得不计算边界是否太小，或者可能更小。但我不明白，为什么我失去了所有的交叉点，我原来的功能。我已经手动计算了交叉点，并且总是在设定的边界中找到一个点。