Distancia entre dos puntos

[[escafandra]]

Empezado por DanielMC

Alguien pudiera traducir el código que puso el Sr.Neftali a c++!!!.

Mira a ver si esta traducción a C++ Builder te sirve.

Código PHP:

#include <Math.hpp>
#include <Math.h>

double RoundTo(double F, int d=0)
{
   double p = pow10(d);
   return floor(F*p+0.5)/p;
}

double WGS84Distance(double plon1, double plat1, double plon2, double plat2, BOOL Almostnul = false)
{
  double d, s, c;
  double lat1, lon1, lat2, lon2, almost1;
  int R = 6378137; // Medium earth radius in meter

  // d = R x arcos [ sin(lat1) x sin(lat2) + cos(lat1) x cos(lat2) x cos(lon2-lon1) ]
  lat1 = DegToRad(plat1);
  lon1 = DegToRad(plon1);
  lat2 = DegToRad(plat2);
  lon2 = DegToRad(plon2);
  if (lon1 == lon2 && lat1 == lat2)
    d = 0;
  else{
    s = sin(lat1) * sin(lat2);
    c = cos(lat1) * cos(lat2) * cos(lon2 - lon1);

    if(Almostnul)
      almost1 = RoundTo((s + c) * 100000000) / 100000000;
    else
      almost1 = (s + c);

    if(almost1 !=1)
      d = R * ArcCos(s + c);
    else
      d = 0;
  }
  return d;
} 


Saludos.

[[Delphius]]

Empezado por DanielMC

Gracias a ambos por responder, mirando el enlace que puso el Sr. Casimiro me aclara unas cuantas cosas, pero la solución final la dan en PHP, la cual probé y funciona perfectamente. Lo de restar las longitudes y latitudes, elevarla y luego sacarle la raíz cuadrada, lo probare a ver qué tal su precisión
Alguien pudiera traducir el código que puso el Sr.Neftali a c++!!!.

Esto va a funcionar siempre y cuando estés hablando en el sistema euclideano.
La función de la distancia, de nuevo, en el sistema euclideano, no es más que ¡el Teorema de Pitágoras aplicado en dos dimensiones! Ahora bien si vas a hacerlo por esta vía tal función debe ser revisada.
Teóricamente es como dices: Sqrt((x2-x1)^2 + (y2-y1)^2)) pero... en el mundo de las computadoras donde aplica la aritemética de punto flotante las cosas deben encararse un pelín diferente.
Tanto para números pequeños como grandes esa fórmula dará underflow y/o overflow respectivamente. En su lugar se emplea la función Hypot, cuya implementación previene ese problema. Está demás decir que la mayoría de los lenguajes ya cuentan con esta función implementada.

Ahora bien, tengo entendido que lo ideal sería que emplearas el desarrollo por medio de la geometría esférica (una geometría NO euclideana). Desconozco si es que el desarrollo propuesto por Neftali y el que te recomendaron que revises en el link con código en PHP aplica esto.

Lo que si me permito hacer un buen llamado de atención. Existen varios sistemas de coordenadas... tienes las geodésicas (como Gauss Kruger), y las geográficas (Longitud/Latitud). Ahora bien... Hay diferentes sistemas de medición... y de donde ponen el punto de origen... Un sistema pone el 0 en UTM, mientras que otros se manejan con el falso este y falso sur... Es decir no cuentan con S, O. Simplemente tienen N y E.

En diferentes partes del mundo emplean diferentes sistemas. ¡Yo que vos acudo a un agrimensor!

Saludos,