[MySQL] MySQL 版本不支持 ST_Distance_Sphere替代方案和解决方案

如果你的 MySQL 版本不支持 ST_Distance_Sphere（常见于 MySQL 5.7.6 以下版本），或者因配置问题无法调用该函数，以下是替代方案和解决方案：

检查 MySQL 版本及功能支持

确认 MySQL 版本：

SELECT VERSION();

若版本低于 5.7.6，则 ST_Distance_Sphere 不可用。

检查空间扩展是否启用：

SHOW VARIABLES LIKE 'have_geometry';

若结果为 YES，表示支持基础空间函数（如 ST_Distance），但不一定包含球面计算函数。

替代方案：使用 ST_Distance（平面距离计算）

如果无需高精度球面距离，可用 ST_Distance 计算平面坐标系下的最小距离（适用于小范围区域，单位与坐标系一致，如度、米等）：

SELECT 
    id,
    name,
    ST_Distance(geom, @pt) AS planar_distance
FROM polygons
ORDER BY planar_distance;

示例结果：
id name planar_distance

1 Polygon1 0
2 Polygon2 ~11.31（度或米）

手动实现球面距离计算（Haversine 公式）

若需精确地理距离，可手动实现球面距离公式。假设多边形和点的坐标是经纬度（WGS84），以下是通过 SQL 函数计算两点间距离的示例：

步骤 1：创建自定义球面距离函数

DELIMITER //
CREATE FUNCTION haversine_distance(
    lat1 FLOAT, lon1 FLOAT,
    lat2 FLOAT, lon2 FLOAT
) RETURNS FLOAT
BEGIN
    DECLARE R FLOAT DEFAULT 6371000; -- 地球半径（米）
    DECLARE dlat FLOAT;
    DECLARE dlon FLOAT;
    DECLARE a FLOAT;
    DECLARE c FLOAT;
    
    SET dlat = RADIANS(lat2 - lat1);
    SET dlon = RADIANS(lon2 - lon1);
    
    SET a = 
        SIN(dlat/2) * SIN(dlat/2) +
        COS(RADIANS(lat1))  COS(RADIANS(lat2)) 
        SIN(dlon/2) * SIN(dlon/2);
    
    SET c = 2 * ATAN2(SQRT(a), SQRT(1-a));
    
    RETURN R * c;
END //
DELIMITER ;

步骤 2：提取多边形顶点坐标并计算最小距离

需先提取多边形的顶点坐标，逐一计算点到顶点的距离，再取最小值：
-- 示例：计算点 (3,3) 到所有多边形顶点的最小球面距离

SELECT 
    p.id,
    p.name,
    MIN(
        haversine_distance(
            3, 3,  -- 查询点的经纬度
            ST_Y(point), ST_X(point)  -- 多边形顶点的经纬度
        )
    ) AS min_distance_meters
FROM polygons p
JOIN (
    -- 提取多边形顶点（需解析GEOMETRY字段）
    SELECT 
        id,
        ST_X(ST_PointN(geom, n)) AS lon,
        ST_Y(ST_PointN(geom, n)) AS lat
    FROM polygons
    CROSS JOIN (SELECT 1 AS n UNION SELECT 2 UNION SELECT 3 UNION SELECT 4) AS numbers
    WHERE n <= ST_NumPoints(geom)
) AS vertices ON p.id = vertices.id
GROUP BY p.id, p.name
ORDER BY min_distance_meters;

升级 MySQL 或迁移到兼容工具

升级 MySQL：升级到 5.7.6+ 或 8.0+ 以原生支持 ST_Distance_Sphere。

使用 PostGIS：若需高性能地理计算，建议迁移到 PostgreSQL + PostGIS，支持更丰富的空间函数（如 ST_DistanceSpheroid）。

注意事项

坐标系一致性：确保所有几何数据使用 WGS84（SRID 4326）。

性能问题：手动解析多边形顶点和计算距离可能对大型数据集性能较差。

精度取舍：ST_Distance 的平面计算在经纬度坐标系下误差较大，适合小范围；Haversine 公式适用于大范围但计算复杂。

通过以上方案，即使 ST_Distance_Sphere 不可用，仍能实现点到多边形的距离计算。