mysql 5.7+ 优先用 st_distance_sphere，注意 point(经度,纬度) 顺序、坐标系一致性和 mbrcontains 配合索引；5.6- 需手写 haversine 并确保弧度转换正确；存储过程须用参数占位符，分页必加 order by distance_m, id。
MySQL 5.7+ 直接用 ST_Distance_Sphere 最省事
如果你的 MySQL 版本 ≥ 5.7 且启用了空间索引（POINT 字段类型），ST_Distance_Sphere 是首选——它内置球面距离计算，精度高、语法短、无需手写 Haversine 公式。
常见错误是传参顺序反了：POINT(longitude, latitude)，不是 (lat, lng)。MySQL 的 POINT 要求先经度后纬度，和多数 API（如百度、高德）习惯相反，一错就返回 0 或报错 NULL。
示例（查离某点 5km 内的所有门店）：
SELECT id, name,
ST_Distance_Sphere(POINT(lng, lat), POINT(-73.9857, 40.7484)) AS distance_m
FROM stores
WHERE ST_Distance_Sphere(POINT(lng, lat), POINT(-73.9857, 40.7484))
-
lng和lat必须是数值型字段（DECIMAL(10,8)或DOUBLE），不能是字符串 - 如果表里存的是百度坐标系（BD-09），
ST_Distance_Sphere会按 WGS84 算，结果偏差可达百米级，需先转坐标系 - 该函数在
WHERE中无法走空间索引，纯计算过滤；要加速得配合MBRContains做初步框选
MySQL 5.6 及以下只能手写 Haversine，注意单位和括号嵌套
5.6 不支持 ST_Distance_Sphere，必须用 Haversine 公式。网上抄来的 SQL 经常漏掉 PI()/180 转弧度，或 ASIN 里开根号没包全，导致 NULL 或负距离。
关键参数含义：
-
6371：地球平均半径（km），用6371000得到米 -
RADIANS(x)比x * PI() / 180更安全，避免浮点误差 -
ROUND(..., 1)建议保留一位小数，整数取整会抹掉 300 米内的差异
推荐写法（兼容 5.6）：
SELECT id, name,
ROUND(6371000 * 2 * ASIN(SQRT(
POW(SIN((RADIANS(lat) - RADIANS(40.7484)) / 2), 2) +
COS(RADIANS(40.7484)) * COS(RADIANS(lat)) *
POW(SIN((RADIANS(lng) - RADIANS(-73.9857)) / 2), 2)
)), 1) AS distance_m
FROM stores
WHERE ...
封装成存储过程时，别把经纬度当字符串拼进 SQL
用 CREATE PROCEDURE 封装距离查询很常见，但新手常犯的错是把输入参数直接拼进字符串再 PREPARE 执行——这既危险（SQL 注入），又低效（每次编译）。
正确做法是用参数占位符：
DELIMITER $$
CREATE PROCEDURE GetNearbyStores(IN in_lng DOUBLE, IN in_lat DOUBLE, IN max_dist_m INT)
BEGIN
SELECT id, name,
ST_Distance_Sphere(POINT(lng, lat), POINT(in_lng, in_lat)) AS distance_m
FROM stores
WHERE ST_Distance_Sphere(POINT(lng, lat), POINT(in_lng, in_lat))
- 所有输入参数都用
IN声明，直接参与计算，不拼 SQL - 如果要用动态表名或列名（极少见），才考虑
CONCAT+PREPARE，但必须严格校验输入 - 存储过程里调
ST_Distance_Sphere和普通查询性能一致，无额外开销
距离排序 + 分页时，LIMIT 前必须 ORDER BY，否则结果不稳定
这是线上最容易被忽略的坑：只写 LIMIT 20 不加 ORDER BY distance_m，MySQL 可能每次返回不同记录，尤其数据有并发更新时。
更隐蔽的问题是“距离相等时排序未定义”：两个点距离都是 1234.5 米，MySQL 可能随机选一个排前面，翻页时同一条记录可能重复出现或丢失。
稳妥写法：
ORDER BY distance_m ASC, id ASC
- 主序按距离升序，次序按主键（或其他唯一字段）保稳定
- 如果用
ST_Distance_Sphere计算，注意它返回的是 double，比较时存在浮点误差，ORDER BY本身已足够，不用额外ROUND - 分页偏移大（如
OFFSET 10000）时，建议改用游标分页（last_id + distance_m 组合条件），避免全表扫
地理距离计算真正难的不是公式，而是坐标系一致性、索引优化和浮点边界处理。只要记住 POINT(lng, lat)、用 RADIANS、加二级排序，大部分场景就稳了。
