Da dove viene il raggio terrestre predefinito in ST_Distance_Sphere?

MySQL dice nei documenti per ST_Distance_Sphere

I calcoli utilizzano una terra sferica e un raggio configurabile. L'argomento facoltativo del raggio dovrebbe essere espresso in metri. Se omesso, il raggio predefinito è 6.370.986 metri. Se l'argomento raggio è presente ma non positivo, si ER_WRONG_ARGUMENTSverifica un errore.

PostGIS dice nei documenti di ST_Distance_Sphere(anche se i documenti non sono più precisi )

Utilizza una terra sferica e un raggio di 6370986 metri.

Da dove hanno preso i 6.370.986 metri di default? WGS84 afferma che il raggio dell'asse maggiore è 6.378.137,0 m. PostGIS che ora utilizza un raggio medio utilizza essenzialmente 6371008.

Guardando il codice

#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)

questo significa

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

Le versioni più recenti sono molto meno efficienti, più complesse e usano Pro4j ma sembrano fare la stessa cosa.

Da dove viene 6370986?

Ok, questo è hilarriuusss . L'ho rintracciato. In una vecchia copia di lwgeom/lwgeom_spheroid.cPostGIS 1.0.0rc4 puoi vedere questo,

/*
 * 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;

Passando ai documenti di earthdistance, troverai questo:

Si noti che, diversamente dalla parte del modulo basata su cubi, le unità sono cablate qui: la modifica della earth()funzione non influirà sui risultati di questo operatore.

E quel numero cablato: EARTH_RADIUSpuò essere visto qui

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

Quindi puoi fare un semplice.

EARTH_RADIUS * MILES_TO_METERS = EARTH_RADIUS_IN_METERS
 3958.747716 * 1609.344        = 6370986.884258304

E tu hai il tuo 6370986.884258304. Naturalmente, basta troncarlo e archiviarlo in un longperché perché no.

Quindi, in sostanza, il raggio di MySQL è stato sollevato da un lavoro pigro di PostGIS che ha convertito un raggio in miglia in metri da un'oscura costante da un modulo PostgreSQL di 20 anni casuale .

earth_distanceè un modulo pre-PostGIS di Bruce Momjian. Con la presente proclamo 6370986 la costante di Bmomjian: una buona approssimazione della Terra in metri per soddisfare MySQL. Anche se forse non per molto.