MySQL GIS Lat/Lon 到 X/Y：哪个是哪个？

Question

MySQL 中的 GIS 点是否应存储为 POINT($latitude $longitude) 或 POINT($longitude $latitude) 是否存在约定？

经度对应于笛卡尔地图上的 X，在视觉上更有意义，北指向上，但常见的说法是“纬度和经度”。

Answer 1

在 MySQL 中，您可能会使用 GeomFromText() 函数在空间字段中插入数据。此函数使用 WKT（众所周知的文本）格式来定义几何图形，并在POINT 情况下，它被定义为：

POINT ($longitude $latitude)

Answer 2

对于在 MySQL 8+ 中使用 GPS 坐标，接受的答案是不正确的，并且它会遇到麻烦（尚未使用以前版本的 MySQL 进行测试）。

TL;DR; 在 MySQL 8+ 中使用“POINT($lat $long)”作为 WKT 字符串，但使用 POINT($long, $lat) 与 POINT() 函数。

完整答案：

在使用 SRID 4326（您应该用于 GPS 坐标系的坐标系）时使用 WKT 表示法作为“POINT($longitude $latitude)”会导致不正确的距离计算，即使在整个过程中始终使用应用程序。请继续阅读以了解详细信息。

例如，让我们考虑一下多伦多的加拿大国家电视塔和纽约的世界贸易中心一号大楼之间的直接距离约为。根据 Google 地图，549,18 公里。

GPS 坐标：

• 加拿大国家电视塔：43.64386666880877，-79.38670551139633
• 世界贸易中心一号大楼：40.689321781458446，-74.04415571126154

预期距离：549.18km

以下查询将产生正确的结果结果：

SELECT 
    ST_DISTANCE(
        ST_GEOMFROMTEXT('POINT(40.689321781458446 -74.04415571126154)', 4326),
        ST_GEOMFROMTEXT('POINT(43.64386666880877 -79.38670551139633)', 4326),
        'metre'
    )
FROM DUAL;

    -- results in 549902.432032006 meters which is around 549.9km (CORRECT)

但是，如果您在 WKT 中首先提供经度（按照接受的答案中的建议），则会计算出错误的距离：

SELECT 
    ST_DISTANCE(
        ST_GEOMFROMTEXT('POINT(-74.04415571126154 40.689321781458446)', 4326),
        ST_GEOMFROMTEXT('POINT(-79.38670551139633 43.64386666880877)', 4326),
        'metre'
    )
FROM DUAL;

    -- results in 601012.8595500318 which is around 601km (WRONG)

正如您所看到的 POINT($long $lat) WKT字符串方法不正确，大约有偏差。与 POINT($lat $long) 方法相比，距离为 51 公里，误差几乎为 10%。实际上，你走得越远，情况就会变得更糟。

解释：

这似乎是因为当 MySQL 在 GPS 坐标上下文中考虑 WKT 字符串时，它会将第一个参数视为纬度，将第二个参数视为经度。尝试运行以下查询：

SELECT 
    ST_Latitude(ST_GEOMFROMTEXT('POINT(40.689321781458446 -74.04415571126154)',4326)) as latitude,
    ST_Longitude(ST_GEOMFROMTEXT('POINT(40.689321781458446 -74.04415571126154)',4326)) as longitude
FROM dual;

-- results in
latitude, longitude
40.689321781458446,-74.04415571126154

请注意，当使用 POINT(x, y) 函数而不是 WKT 字符串时，情况恰恰相反！

例子：

SELECT 
    ST_DISTANCE(
        ST_SRID(POINT(-74.04415571126154, 40.689321781458446), 4326),
        ST_SRID(POINT(-79.38670551139633, 43.64386666880877), 4326),
        'metre'
    )
FROM DUAL;

-- results in 549902.432032006 meters which is around 549.9km (CORRECT)