如何使用单条SQL语句查找距离最近的坐标对？

核心实现注意点 #

要满足单SQL、无变量辅助的要求，必须先规避两个常见的坑：

• 杜绝重复配对：自连接生成坐标对时，必须加a.id < b.id的约束，避免同一个坐标对被正反计算两次（比如同时出现(1,2)和(2,1)两组无效重复配对），既减少计算量也避免结果重复。
• 降低计算冗余：找最近点的过程本质是比较距离大小，短距离近似场景下不需要对距离做开平方运算，直接对比距离平方值的大小即可，平方值最小对应的真实距离一定最小，能大幅减少计算开销。

可直接运行的单条SQL实现 #

适配示例的三列坐标表结构（id、lat纬度、long经度），完整SQL不需要临时变量、不需要分段执行：

SELECT
  a.id AS first_point_id,
  a.lat AS first_point_lat,
  a.long AS first_point_long,
  b.id AS second_point_id,
  b.lat AS second_point_lat,
  b.long AS second_point_long,
  POWER(a.lat - b.lat, 2) + POWER(a.long - b.long, 2) AS distance_squared
FROM coordinates a
-- 核心去重连接条件
JOIN coordinates b
  ON a.id < b.id
-- 按距离从小到大排序，取第一行就是最近的坐标对
ORDER BY distance_squared ASC
LIMIT 1;

示例数据计算验证 #

针对给出的三行示例坐标：

id lat long
1  -12 35 
2  -1  100
3  -9  3

所有不重复配对的计算结果：

• 配对(1,2)：距离平方值为(-12+1)^2 + (35-100)^2 = 4346
• 配对(1,3)：距离平方值为(-12+9)^2 + (35-3)^2 = 1033
• 配对(2,3)：距离平方值为(-1+9)^2 + (100-3)^2 = 9473
  SQL最终返回的最小距离配对就是id为1和3的坐标对，和手工计算结果完全一致。

精度升级方式 #

如果需要计算真实的球面最短距离，不需要修改整体SQL结构，只需要把原有距离平方的计算表达式替换成内嵌的Haversine公式即可，全程依然是单条SQL、无外部变量：

替换后的距离计算表达式（结果单位为公里）：
2 * 6371 * ASIN(SQRT(POWER(SIN(RADIANS(b.lat - a.lat)/2), 2) + COS(RADIANS(a.lat)) * COS(RADIANS(b.lat)) * POWER(SIN(RADIANS(b.long - a.long)/2), 2)))
替换完成后保留原有的排序、取第一行逻辑，即可得到球面距离下的最近坐标对。

内容的提问来源于stack exchange，提问作者pp182