Вычисление точек пересечения двух линий возвращает -1.#IND в CSV-файле
В моем программном обеспечении я пытаюсь вычислить точки пересечения двух линий. Я начинаю с четырех двойных переменных для первой строки (две для х / двух для у) и четырех двойных переменных для второй строки. Математика мне понятна, и на бумаге это работает. Для дальнейшей обработки я должен сохранить точки пересечения в файле 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минус m20 - нет точки пересечения.
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петля для сохранения вектора пересекается в CSV-файле.
for(int i=0; i<intersects.size(); i++)
{
file << "Intersects" << ";" << intersects[i][0] << ";" << intersects[i][1] << ";" << intersects[i][2] << endl;
}
Но когда я открывал файл, у меня только что была упомянутая ошибка. Я уже пробовал другой тип данных, такой как float, но ничего не изменилось. Я знаю, что первая и последняя части кода работают, потому что я уже использовал их с другим видом математики, однако мой математический подход был выбран неправильно, поэтому мне пришлось его изменить. Поэтому я думаю, что проблема должна быть в состоянии для if-заявление.
Изменить: с помощью комментариев, я теперь адаптировал свой код. Итак, это мой новый код:
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 к моему коду я должен был бы вычислить, являются ли границы маленькими или могут быть меньшими. Но я не понимаю, почему я теряю все пересечения с моей первоначальной функцией. Я уже рассчитал пересечения вручную и всегда находил точку в заданных границах.