ST_Distance_Sphere中的默认地球半径从何而来？

MySQL在文档中说 ST_Distance_Sphere

计算使用球形地球和可配置的半径。可选的radius参数应以米为单位。如果省略，则默认半径为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仍然来自哪里？

好的，这是希拉里乌斯。我追踪到了。在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。当然，只需将其截断并将其存储在a中，long因为为什么不这样做。

因此，在本质上，在MySQL的半径从懒惰复制工作从PostGIS的提升是转换的半径英里米，从一个不起眼的常数从随机20岁的PostgreSQL模块。

earth_distance是Bruce Momjian编写的PostGIS之前的模块。我在此宣布6370986为Bmomjian常数：以米为单位的近似地球，可以满足MySQL的要求。虽然可能不会很久。