Calcule distance sur POINT

**bannik** · 05/03/2013, 11h26

Bonjour

Je recherche une méthode de calcul de la distance entre deux point enregistrés en base de donnée (colonnes de type POINT)

Je connais la méthode a partir des latitudes longitudes:

Code :

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
17
18
19
20
21
 
 
CREATE FUNCTION `Distance`(lat1 DOUBLE, lng1 DOUBLE, lat2 DOUBLE, lng2 DOUBLE) RETURNS double
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 (6378.137 * 2 * ATAN2(SQRT(a), SQRT(1 - a)));
END

Je suppose que vu que ce sont des données géo ça devrai être plus simple.

Quelqu'un connais la méthode?

Merci

**bannik** · 05/03/2013, 13h27

Bon, je pense avoir la solution:

Code :

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

1
2
 
SELECT (GLength(LineStringFromWKB(LineString(point1,point2)))) * 111.18957696  as dist

Pour un retour en Km

Quelqu'un pour confirmer?

Merci

**ericd69** · 05/03/2013, 18h03

salut,

la première est une approximation juste la deuxième semble pouvoir engendrer des erreurs assez grossière partout où on en parle...

après c'est à tester en terme de performance et d'écart de valeur...

ta première fonction est très mal optimisée et en fait il existe 2 formule à peu près équivalentes...
moi j'ai choisi l'autre et introduite le support de l'altitude (surtout pour des points espacés de quelques km sinon ça prend évidemment un facteur d'erreur de plus en plus important):

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
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$$

radians() ne fait que multiplier une constante à ton nombre donc tu peux mutualiser c'est transformation dans certains cas et ainsi limiter aussi le nombre de variables allouées...
exemple si tu as:

Code :

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

1
2
x=radians(a),y=radians(b)
z=cos(x-y)

ça équivaut à:

Code :

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

z=cos(radians(a-b))

si a et b ne sont utilisé que pour faire a-b alors tu peux éviter nombre d'opérations comme ça, tu le vois, ou en stockant radians(a-b) s'il apparait à divers endroits par exemple...