Version imprimable
un petit problème pour les calculs des doubles.
j'utilise evc4, voici mon code:
double lat1=45.8541;
double lat2=45.854;
double deltaLat=lat1-lat2; //cette ligne donne bien 0.0001
deltaLat=deltaLat/2; //et là 5.00000000016
vous savez qu'est-ce qu'il s'est passé ??
Je ne voit pas de pb chez moi...
j'obtiens 5.0000000001659828e-005
Tu n'aurrais pas oublié de regarder le e-005 ? :mouarf:
oui, effectivement 8O
mais problème n'est pas résolu, en fait, j'ai fait une fonction qui calcul la distance entre 2 points latitude/longitude, mais le résultat est totalement faux.
#include <math.h>
double Distance(double Lat1,double Lng1,double Lat2,double Lng2)
{
double rayonTerre=6371;
double deltaLat=Lat1-Lat2;
double deltaLng=Lng1-Lng2; //jusqu'à ici c'est juste
double a=sin(deltaLat/2)*sin(deltaLat/2)+
cos(Lat1)*cos(Lat2)*sin(deltaLng/2)*sin(deltaLng/2);
double c=2*atan2(sqrt(a),sqrt(1-a));
double distance=rayonTerre*c;
return distance; //le résultat est complètement faux, du gerne 0.7 au lieu de 0.02
}
Quelles sont tes données d'entrées ?
Les données entrées:
Lat1=45.854 Lng1=4.7319
Lat2=45.8541 Lng2=4.7321
J'attend comme résultat 0.01907 (en théorie) mais j'ai 0.740652288866354 :cry:
L'incertitude (740 mètres par rapport à 20 mètres) est quand même inacceptable...
Je pense que tu devrais convertir tes angles en radians...
:roll:
Bonjour,
Je ne prendrai pas la peine de lire ton code.
Une telle erreur vient forcément de l'algorithme, car si c'était le compilateur ou/et le processeur, ça se saurait depuis longtemps, n'est-ce pas ?
Donc revois ton calcul, en faisant quand même attention aux approximations inévitables quand on calcule avec des réels (par exemple, faire la différence entre deux valeurs très proches peut donner un résultat loin de la valeur attendue, idem diviser par une valeur très faible...).
Je suis d'accord avec thewho
2 choses me choquent dans ton algorithme :
- Tu utilises apparemment des degrés plutôt que des radians
- Tu renvois des km plutôt que des m
Place toi toujours dans un système universel ;)
j'utilise le formulaire d'Haversine, les lat/lng sont sous forme de degré décimale.
R = earth’s radius (6371km)
Δlat = lat2− lat1
Δlong = long2− long1
a = sin²(Δlat/2) + cos(lat1).cos(lat2).sin²(Δlong/2)
c = 2.atan2(√a, √(1−a))
d = R.c
donc l'algorithme n'a aucun erreur.
sauf que les fonction sin, atan2, ... prennent des radians et non des degrés comme la formule le laisse à penser. Fais le test au lieu de dire que ce n'est pas ça.
Ouiiiiii, c'est bien ça le problème ;)
j'ai converti les degrés en radians, et ça donne des bonnes résultats.
Merci beaucoup les gars !!