Discussion :

[Hoopsy]

Bonjour,

Je viens de remarquer quelque chose d'assez bizarre quand j'utilise le module Math::Trig qw(great_circle_distance deg2rad) pour calculer la distance entre deux points géographiques, alors quand je prends l'exemple suivant :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
 
    use Math::Trig qw(great_circle_distance deg2rad);
    sub NESW { deg2rad($_[0]), deg2rad(90 - $_[1])}
    my @L = NESW(5.761080000, 36.822300000);
    my @T = NESW(5.76108, 36.8223);
    my $km = great_circle_distance(@L, @T, 6378);
    print("$km\n");

Logiquement je dois trouver 0 km et bien non je trouve 0,0000950396060943604 km.
J'ai utilisé l'exemple le plus répondu sur internet à ce sujet :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
 
   use Math::Trig qw(great_circle_distance deg2rad);
   sub NESW { deg2rad($_[0]), deg2rad(90 - $_[1]) }
   my @L = NESW( -0.5, 51.3);
   my @T = NESW(139.8, 35.7);
   my $km = great_circle_distance(@L, @T, 6378); # About 9600 km.
   print("$km\n");

Moi je trouve 9605.26637021388 km et non 9600 km .

Est-ce que c'est la même chose pour vous ? Peut-être que j'utilise mal le module ou alors le calcul n'est pas très précis .

Merci de m'aider.
H.

[Beniou]

Bonjour,

Je trouve les mêmes résultats que toi : en regardant de plus prêt la fonction "great_circle_distance" je pense savoir d'où cela vient.

Voici la fonction :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
 
sub great_circle_distance {
    my ( $theta0, $phi0, $theta1, $phi1, $rho ) = @_;
 
    $rho = 1 unless defined $rho; # Default to the unit sphere.
 
    my $lat0 = pip2 - $phi0;
    my $lat1 = pip2 - $phi1;
 
    return $rho *
	acos_real(cos( $lat0 ) * cos( $lat1 ) * cos( $theta0 - $theta1 ) +
		   sin( $lat0 ) * sin( $lat1 ));
}

En essayant de récupérer les différentes valeurs dans la fonction et en faisant des prints comme suit en utilisant ton premier exemple (celui qui doit afficher 0) :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
 
my $value = cos( $lat0 ) * cos( $lat1 ) * cos( $theta0 - $theta1 ) +
		   sin( $lat0 ) * sin( $lat1 );
print "value : $value => acos_real($value) = ".acos_real($value)."\n";
print "acos_real(1) = ".acos_real(1)."\n";

on obtient :

value : 1 => acos_real(1) = 1.49011611938477e-08
acos_real(1) = 0

On remarque donc qu'il met value à 1 et le acos_real ne met pas à 0 contrairement à la ligne suivante où j'affiche le acos_real(1) = 0.

Ensuite j'ai essayé d'afficher le valeur de value avec la fonction sprintf pour avoir 15 et 20 chiffres après la virgule

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
 
 print sprintf("%.15f",$value)."\n";
 print sprintf("%.20f",$value)."\n";

on obtient respectivement:

1.000000000000000
0.99999999999999988898

Je pense donc que le problème vient de la gestion des "nombres réels" et leur approximation dans calcul de cette valeur "cos( $lat0 ) * cos( $lat1 ) * cos( $theta0 - $theta1 ) + sin( $lat0 ) * sin( $lat1 )"

[Hoopsy]

Bonjour Beniou,

Je vous remercie pour vos explications, il doit surement y avoir un problème d'approximation, c'est dommage pour les utilisateurs de ce module.

J'ai masqué ce bémol en prenant que deux chiffres après la virgule du résultat final, et là je peux enfin avoir mon 0 km , mais ceci dit le résultat n'est pas vraiment très très juste pour d'autres jeux de données.

Merciiiiiiii beaucoup.
H.