mysqlの緯度と経度を使用して距離を計算します

Question

人々は緯度と経度を入力してplacesテーブルから最も近い場所を表示できるスクリプトを作成しようとしています。 "like"または "="を使用すると、完全一致の結果が取得されます。しかし、私は正確な一致がない場合、最も近い場所を表示できるようにしたい。 MySQLで距離を計算する方法はありますか？ 助けてください。ありがとう

Answer 1

ビンセント式https://en.wikipedia.org/wiki/Vincenty%27s_formulaeを使用できます。

MySQLではこのコードとの距離を計算する関数を作成することができます。

DELIMITER $$ 
DROP FUNCTION IF EXISTS vincenty_distance$$ 
CREATE FUNCTION vincenty_distance(lat1 DOUBLE, lon1 DOUBLE, lat2 DOUBLE, lon2 DOUBLE) RETURNS DOUBLE 
BEGIN 
    DECLARE TO_RAD DOUBLE; 
    DECLARE a INT; 
    DECLARE b DOUBLE; 
    DECLARE f DOUBLE; 
    DECLARE L DOUBLE; 
    DECLARE U1 DOUBLE; 
    DECLARE U2 DOUBLE; 
    DECLARE sinU1 DOUBLE; 
    DECLARE cosU1 DOUBLE; 
    DECLARE sinU2 DOUBLE; 
    DECLARE cosU2 DOUBLE; 
    DECLARE lambda DOUBLE; 
    DECLARE lambdaP DOUBLE; 
    DECLARE iterLimit INT; 
    DECLARE sinLambda DOUBLE; 
    DECLARE cosLambda DOUBLE; 
    DECLARE sinSigma DOUBLE; 
    DECLARE cosSigma DOUBLE; 
    DECLARE sigma DOUBLE; 
    DECLARE sinAlpha DOUBLE; 
    DECLARE cosSqAlpha DOUBLE; 
    DECLARE cos2SigmaM DOUBLE; 
    DECLARE C DOUBLE; 
    DECLARE D DOUBLE; 
    DECLARE E DOUBLE; 
    DECLARE uSq DOUBLE; 
    DECLARE deltaSigma DOUBLE; 
    DECLARE s DOUBLE; 
    SET TO_RAD = pi()/180; /*converts degree to raians*/ 
    SET a = 6378137; 
    SET b = 6356752.3142; 
    SET f = 1/298.257223563; /* WGS-84 ellipsoid params*/ 
    SET L = (lon2-lon1) * TO_RAD; 
    SET U1 = atan((1 - f) * tan(lat1 * TO_RAD)); 
    SET U2 = atan((1 - f) * tan(lat2 * TO_RAD)); 
    SET sinU1 = sin(U1); 
    SET cosU1 = cos(U1); 
    SET sinU2 = sin(U2); 
    SET cosU2 = cos(U2); 
    SET lambda = L; 
    SET iterLimit = 100; 
    REPEAT 
     SET sinLambda = sin(lambda); 
     SET cosLambda = cos(lambda); 
     SET sinSigma = sqrt((cosU2 * sinLambda) * (cosU2 * sinLambda) + (cosU1 * sinU2 - sinU1 * cosU2 * cosLambda) * (cosU1 * sinU2 - sinU1 * cosU2 * cosLambda)); 

     IF 0 = sinSigma THEN 
      RETURN 0; /* co-incident points*/ 
     END IF; 
     SET cosSigma = sinU1 * sinU2 + cosU1 * cosU2 * cosLambda; 
     SET sigma = atan2(sinSigma, cosSigma); 
     SET sinAlpha = cosU1 * cosU2 * sinLambda/sinSigma; 
     SET cosSqAlpha = 1 - sinAlpha * sinAlpha; 
     IF (cosSqAlpha = 0) THEN 
      SET cos2SigmaM = 0; /* equatorial line: cosSqAlpha = 0 (§6)*/ 
     ELSE 
      SET cos2SigmaM = cosSigma - 2 * sinU1 * sinU2/cosSqAlpha; 
      SET C = f/16 * cosSqAlpha * (4 + f * (4 - 3 * cosSqAlpha)); 
     END IF; 
     SET lambdaP = lambda; 
     SET lambda = L + (1 - C) * f * sinAlpha * (sigma + C * sinSigma * (cos2SigmaM + C * cosSigma * (-1 + 2 * cos2SigmaM * cos2SigmaM))); 
     SET iterLimit = iterLimit - 1; 
    UNTIL (abs(lambda - lambdaP) <= 0.0000000001 && iterLimit = 0) END REPEAT; 
    SET uSq = cosSqAlpha * (a * a - b * b)/(b * b); 
    SET D = 1 + uSq/16384 * (4096 + uSq * (-768 + uSq * (320 - 175 * uSq))); 
    SET E = uSq/1024 * (256 + uSq * (-128 + uSq * (74 - 47 * uSq))); 
    SET deltaSigma = E * sinSigma * (cos2SigmaM + E/4 * (cosSigma * (-1 + 2 * cos2SigmaM * cos2SigmaM) - E/6 * cos2SigmaM * (-3 + 4 * sinSigma * sinSigma) * (-3 + 4 * cos2SigmaM * cos2SigmaM))); 
    SET s = b * D * (sigma - deltaSigma); 

    RETURN round(s, 3); /* round to 1mm precision*/ 
END$$ 
DELIMITER ; 

使用するには、単に緯度と経度パラメータで機能vincenty_distanceを呼び出す必要があり、例えば：

SELECT vincenty_distance(47.6593,10.97647,46.2512010,10.069972); 

Vincenty Formulaは距離を計算する最も正確な方法の1つです。なぜなら、それは地球の姿を偏平な回転楕円体として仮定しているからです。

アルゴリズムは距離をメートルで返します。精度は1mmです。

そうでなければあなたは（わずかな距離で有効な）速いが、不正確な計算を好むならば、あなたの使用ピタゴラスの定理： Calculating distance (pythagoras) and running count in sql query

距離を計算するための他の式の多くがありますが、あなたがに応じていずれかを選択する必要があります精度とあなたが必要とするパフォーマンス：

https://en.wikipedia.org/wiki/Geographical_distance

MySQL Great Circle Distance (Haversine formula)