Кратчайшее расстояние до точки

Question

У меня есть 2 вектора, один (вектор1 структур (Point)) заполняется X количеством точек, а другой (vector2 of structs (PrimeTemplate)) заполняется Y количеством точек. Я хочу найти все значения ниже порога, и я чувствую, что мой код просто этого не делает. На данный момент я просто проигнорирую, если одна точка отображает более одного другого. Что мне не хватает? Я генерирую только несколько моментов, и я знаю, что я должен получать больше.

struct Template{ 
    int tempX; 
    int tempY; 
}; 

struct PrimeTemplate{ 
    double tempX; 
    double tempY; 
}; 

int matches = 0; 
for (int outerLoop = 0; outerLoop < vector1 .size(); outerLoop++) 
    { 
     for (int innerLoop = 0; innerLoop < vector2.size(); innerLoop++) 
     { 
      double tempEuclidianX = std::pow(abs(vector1 [outerLoop].tempX - vector2[innerLoop].tempX), 2.0); 
      double tempEuclidianY = std::pow(abs(vector1 [outerLoop].tempY - vector2[innerLoop].tempY), 2.0); 
      double Euclidian = sqrt(tempEuclidianX + tempEuclidianY); 
      if (Euclidian <= 5) //less than threshold 
      { 
       matches++; 
      } 
     } 
    } 

Пример ввода из файла будет выглядеть следующим образом (два разных файла, случайные числа) (не беспокойтесь о получении данных, все это есть)

245 21 
452 54 
124 68 
485 78 
111 29 
97 75 
78 113 
300 124 
411 101 

Answer 1

Это был метод, который я использовал для вычисления кратчайших расстояний между парой. Он перемещается по текстовому файлу и загружает векторы, которые вы видите. Оказалось, что вопрос был в моей реализации до этого кода, который был нормализацией биометрических точек.

for (int outerLoop = 0; outerLoop < Tvector.size(); outerLoop++) 
     { 
      for (int innerLoop = 0; innerLoop < QPrimeVector.size(); innerLoop++) 
      { 
       double tempEuclidianX = std::pow((QPrimeVector[innerLoop].tempX - Tvector[outerLoop].tempX), 2.0); 
       double tempEuclidianY = std::pow((QPrimeVector[innerLoop].tempY - Tvector[outerLoop].tempY), 2.0); 
       double Euclidian = sqrt(tempEuclidianX + tempEuclidianY); 
       if (Euclidian <= THRESHOLD) //less than threshold and not taken already 
       { 
        if (Euclidian < minEuclidian) 
        { 
         minEuclidian = Euclidian; 
         if (!Tvector[outerLoop].marked) 
         { 
          matched = innerLoop; 
         } 
        } 
       } 
       if (matched != -1) 
       { 
        matches++; 
       } 
       matched = -1; 
       minEuclidian = 10; 
      } 

      if (matches > masterMatchCount) 
      { 
       masterMatchCount = matches; 
       deltaThetaMaster = deltaTheta; 
       deltaXMaster = deltaX; 
       deltaYMaster = deltaY; 
      } 
     } 
     for (int reset = 0; reset < Tvector.size(); reset++) 
     { 
      Tvector[reset].marked = false; //reset all the matches 
     } 
     QPrimeVector.clear(); 
    } 

Answer 2

Что с ваш код является то, что вы используете abs() перед возведением в квадрат. Нет необходимости брать абсолютное значение вообще, прежде чем возводить в квадрат, конечно, но если вы собираетесь, то вы хотите использовать fabs, как только абс принимает и возвращает целое число. Возможно, это дополнительное округление может быть причиной неправильного ответа.