Discussion :

[xavier_dcf]

Bonjour je dois creer une fonction pour calculer des distances entre deux points ayant des coordonnées géographiques .
J'ai trouvé ces formules :
-d=R*acos(sin(lat1)*sin(lat2)+cos(lat1)*cos(lat2)*cos(lon1-lon2))
d=R*2*asin(sqrt((sin((lat1-lat2)/2))^2 + cos(lat1)*cos(lat2)*(sin((lon1-lon2)/2))^2))
avec des latitudes longitudes en radian . Est ce que qqun peut confirmer ces formules déja?

Ensuite j'ai trouvé une fonction pour faire cela :
CREATE FUNCTION get_distance_metres (lat1 DOUBLE, lng1 DOUBLE, lat2 DOUBLE, lng2 DOUBLE) RETURNS DOUBLE
DETERMINISTIC
NO SQL
BEGIN
DECLARE rlo1 DOUBLE;
DECLARE rla1 DOUBLE;
DECLARE rlo2 DOUBLE;
DECLARE rla2 DOUBLE;
DECLARE dlo DOUBLE;
DECLARE dla DOUBLE;
DECLARE a DOUBLE;

SET rlo1 = RADIANS(lng1);
SET rla1 = RADIANS(lat1);
SET rlo2 = RADIANS(lng2);
SET rla2 = RADIANS(lat2);
SET dlo = (rlo2 - rlo1) / 2;
SET dla = (rla2 - rla1) / 2;
SET a = SIN(dla) * SIN(dla) + COS(rla1) * COS(rla2) * SIN(dlo) * SIN(dlo);
RETURN (6378137 * 2 * ATAN2(SQRT(a), SQRT(1 - a)));
END|
DELIMITER ;

qui calcule la deuxième formule (sauf le return mais ca doit etre une formule équivalente). Cependant mysql dit qu'il y a des erreurs au niveau des declare et return!
qu'est ce qui ne va pas??
Merci

[Fred_34]

Il te manque juste DELIMITER | au début.


La première formule est probablement suffisante (et surtout plus simple):

Code :

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

1
2
3
4
5
6
7
8
9
DELIMITER |
DROP FUNCTION IF EXISTS get_distance_metres|
CREATE FUNCTION get_distance_metres (lat1 DOUBLE, lng1 DOUBLE, lat2 DOUBLE, lng2 DOUBLE) RETURNS DOUBLE
DETERMINISTIC
NO SQL
BEGIN
RETURN 6366000 * ACOS(SIN(RADIANS(lat1)) * SIN(RADIANS(lat2)) + COS(RADIANS(lat1)) * COS(RADIANS(lat2))* COS(RADIANS(lng1)-RADIANS(lng2)));
END|
DELIMITER ;

Le résultat a l'air assez correct puisque en faisant
select get_distance_metres (43.592631,3.849873,43.593435,3.849313)J'ai obtenu 100,05. ( J'ai pris sur google map des coordonnées d'un terrain de rugby.

En prenant 6378137 comme rayon de la terre ça donne 100,24.

[xavier_dcf]

merci!
ca marche bien

[ericd69]

salut,

pour info, la même fonction tenant compte de la hauteur des points et optimisée en terme de code...

Code sql :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
delimiter $$
drop function if exists distance$$
create function distance(la1 double, lo1 double, alt1 double, la2 double, lo2 double, alt2 double)
returns double
deterministic
begin
	declare rla1,rla2,d,h double;
	set rla1=radians(la1),rla2=radians(la2);
	set d=6378137*acos(sin(rla1)*sin(rla2)+cos(rla1)*cos(rla2)*cos(radians(lo2-lo1))),h=alt2-alt1;
	if h=0 then
		return d;
	else
		return sqrt(d*d+h*h);
	end if;
end$$
delimiter ;

en factorisant les conversion et minimisant les variables nécessaires