ST_Distance_Sphereのデフォルトの地球半径はどこから来ますか？

MySQLのドキュメントでは ST_Distance_Sphere

計算では、球状の地球と構成可能な半径を使用します。オプションの半径引数はメートル単位で指定する必要があります。省略した場合、デフォルトの半径は6,370,986メートルです。radius引数が存在するが正ではない場合、ER_WRONG_ARGUMENTSエラーが発生します。

PostGISはのドキュメントで述べてST_Distance_Sphereいます（ドキュメントはもう正確ではありませんが）

球状の地球と6370986メートルの半径を使用します。

どこからデフォルトの6,370,986メートルを取得しましたか？ WGS84によると、長軸半径は6,378,137.0 mです。現在、平均半径を使用しているPostGISは、基本的に6371008を使用しています。

コードを見る

#define WGS84_MAJOR_AXIS 6378137.0
#define WGS84_INVERSE_FLATTENING 298.257223563
#define WGS84_MINOR_AXIS (WGS84_MAJOR_AXIS - WGS84_MAJOR_AXIS / WGS84_INVERSE_FLATTENING)
#define WGS84_RADIUS ((2.0 * WGS84_MAJOR_AXIS + WGS84_MINOR_AXIS ) / 3.0)

つまり

-- SELECT 6378137.0 - 6378137.0 / 298.257223563;
WGS84_MINOR_AXIS = 6356752.314245179498
-- SELECT ( 2.0 * 6378137.0 + ( 6378137.0 - 6378137.0 / 298.257223563) ) / 3.0;
WGS84_RADIUS = 6371008.771415059833

新しいバージョンははるかに効率が悪く、より複雑で、Pro4jを使用していますが、同じことをしているようです。

それでも6370986はどこから来たのですか？

OK、これはhilarriuusssです。これを追跡しました。lwgeom/lwgeom_spheroid.cPostGIS 1.0.0rc4の古いコピーでは、これを見ることができます、

/*
 * This algorithm was taken from the geo_distance function of the 
 * earthdistance package contributed by Bruno Wolff III.
 * It was altered to accept GEOMETRY objects and return results in
 * meters.
 */
PG_FUNCTION_INFO_V1(LWGEOM_distance_sphere);
Datum LWGEOM_distance_sphere(PG_FUNCTION_ARGS)
{
        const double EARTH_RADIUS = 6370986.884258304;

のドキュメントに進むとearthdistance、次のことがわかります。

モジュールのキューブベースの部分とは異なり、ユニットはここで固定されていることに注意してください。earth()関数を変更しても、この演算子の結果には影響しません。

そしてそのハードワイヤード番号：EARTH_RADIUSここで見ることができます

/* Earth's radius is in statute miles. */
static const double EARTH_RADIUS = 3958.747716;

だから、簡単にできます。

EARTH_RADIUS * MILES_TO_METERS = EARTH_RADIUS_IN_METERS
 3958.747716 * 1609.344        = 6370986.884258304

そして、あなたはあなたの6370986.884258304。もちろん、それを切り捨ててに格納してlongください。

だから、本質的には、MySQLでの半径がたのPostGISからレイジーコピージョブから持ち上げメートルにマイルで半径を変換からあいまいな定数からランダムに20歳のPostgreSQLモジュール。

earth_distanceBruce MomjianによるPre-PostGISモジュールです。私はここで6370986をBmomjian定数と宣言します。これは、MySQLを満足させるためのメートル単位の地球の適切な近似です。たぶん長い間ではありませんが。