PostgreSQL经纬度距离查询实现及radians函数报错解决

PostgreSQL 抛出 function radians(character varying) does not exist 的核心原因是：kr_user 表中的 latitude、longitude 字段存储类型为变长字符串（character varying），而内置的 radians() 函数仅接受数值类型（双精度浮点数、数值型等）的输入参数，无法直接匹配字符串类型的调用，因此触发函数签名不匹配的报错。

给传入 radians() 的经纬度字符串添加显式类型转换，将其转为可被函数识别的双精度浮点数即可，修正后可正常运行的SQL如下：

SELECT id, 
( 3959 * acos( 
    cos( radians(14.685327) ) 
    * cos( radians( CAST(latitude AS DOUBLE PRECISION) ) ) 
    * cos( radians( CAST(longitude AS DOUBLE PRECISION) ) - radians(80.054889) ) 
    + sin( radians(14.685327) ) * sin(radians( CAST(latitude AS DOUBLE PRECISION) )) 
) ) AS distance 
FROM kr_user
HAVING distance < 50
ORDER BY distance

注意：公式中的3959是地球半径的英里制数值，如果需要以公里为单位计算距离，将这个值替换为6371即可。
长期优化建议：尽量在表设计阶段就将经纬度字段设置为数值类型（如 double precision）存储，避免每次查询都做类型转换的性能损耗，也能减少空值、非数字字符串导致的转换报错。

手写球面距离公式的方案仅适合小数据量的临时查询，存在两个明显短板：一是逐行计算三角函数无法走普通索引，数据量超过10万行后查询延迟会明显上升；二是浮点精度问题容易导致acos函数的输入值超出[-1,1]的合法范围，触发数学域计算错误。
更成熟的实现路径包括：

• 启用PostgreSQL原生的PostGIS地理空间扩展，这是数据库层面专门处理空间数据的工具，内置高精度的经纬度距离计算函数，支持空间GiST索引，百万级数据的距离检索也能保持毫秒级响应
• 适配Laravel开发场景：可以通过查询构造器直接封装距离计算逻辑，不需要硬拼接原生SQL；如果引入PostGIS扩展，可通过Eloquent作用域封装通用的按距离检索方法，在业务代码中直接复用查询逻辑

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