DBで緯度経度から距離を計算する【ST_Distance_Sphere】

FUNCTION ST_Distance_Sphere does not exist

DBから緯度経度で距離を測りたかったので調べたら、
ST_Distance_Sphereというのが出てきたので、使ってみると

FUNCTION ST_Distance_Sphere does not exist

とでてきました。
調べてみると、関数作ってやってるのがちらほら見つかったので、
それを試してみました！

ちなみに、今回はGoogle Mapと比べたときに近いものを正とすることにしました！

ST_Distance_Sphere対応策1

CREATE FUNCTION st_distance_sphere(pt1 POINT, pt2 POINT)
RETURNS double(10,2)

RETURN 6371 * 2 * ASIN(
    SQRT(
        POWER(SIN((ST_Y(pt2) - ST_Y(pt1)) * pi()/180 / 2), 2) +
        COS(ST_Y(pt1) * pi()/180 ) *
        COS(ST_Y(pt2) * pi()/180) *
        POWER(SIN((ST_X(pt2) - ST_X(pt1)) * pi()/180 / 2), 2)
    )
);

https://stackoverflow.com/questions/59370894/mariadb-no-equivalent-function-for-st-distance-sphere

SELECT  st_distance_sphere(ST_GeomFromText('POINT(35.41675663857636 139.4444246043602)', 4326), ST_GeomFromText('POINT(35.4273603630893 139.49807795163338)', 4326));

ST_Distance_Sphere対応策1の結果

近しい数値がでていますが、なんか納得できないですね、原因もわからないですし

ST_Distance_Sphere対応策2

CREATE FUNCTION st_distance_sphere2(latitude_1 double, longitude_1 double, latitude_2 double, longitude_2 double)
RETURNS double
return 6371 * acos(
    cos(radians(latitude_2))
    * cos(radians(latitude_1))
    * cos(radians(longitude_1) - radians(longitude_2))
    + sin(radians(latitude_2))
    * sin(radians(latitude_1))
);

https://blogenist.jp/2019/09/02/9037/#2

式を見ても正直わかりませんが、さっきのよりわかる人にはわかる感じですね。
おそらく球面距離計算がこれなんですかね？

SELECT st_distance_sphere(35.41675663857636, 139.4444246043602, 35.4273603630893, 139.49807795163338);

ST_Distance_Sphere対応策2の結果

ほぼ、同じですね、
ということは、地球は楕円型でスフィアではなくスフィロイドというらしいのですが
Googleでもこんな感じでやってるのかな？

6371ってなに？

地球の半径だそうです！
https://www.google.com/search?q=地球の半径は何キロですか

4326ってなに？

ちなみに、4326というのは「空間参照系」「WGS84・緯度経度座標」というものらしいです。

https://lonlat.info/spatialreferencesystem/