Вычисление точек пересечения двух линий возвращает -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 к моему коду я должен был бы вычислить, являются ли границы маленькими или могут быть меньшими. Но я не понимаю, почему я теряю все пересечения с моей первоначальной функцией. Я уже рассчитал пересечения вручную и всегда находил точку в заданных границах.

-3

Решение

Проверьте деление на ноль. AB_x, CD_x может быть 0.

1

Другие решения

Других решений пока нет …

По вопросам рекламы [email protected]