如何编写PostgreSQL查询，按给定经纬度由近及远排序数据

下面分两种常见场景给出解决方案，根据你的数据库环境选择合适的方式：

1. 原生PostgreSQL实现（无需额外扩展） #

如果你的数据库没有安装PostGIS扩展，可以用Haversine公式计算两点间的球面距离，再按距离排序。假设你的表名为locations，包含存储纬度的lat字段和经度的lon字段，目标经纬度为37.7749（纬度）和-122.4194（经度）：

SELECT
    id,
    name,
    lat,
    lon,
    -- 计算距离（单位：公里）
    6371 * acos(
        cos(radians(37.7749)) * cos(radians(lat)) * cos(radians(lon) - radians(-122.4194)) +
        sin(radians(37.7749)) * sin(radians(lat))
    ) AS distance_km
FROM locations
ORDER BY distance_km ASC;

关键说明： #

• 6371是地球平均半径（公里），如果需要英里单位，替换为3956即可
• radians()函数用于将度数转换为弧度，因为三角函数要求输入为弧度
• 结果按distance_km升序排列，即距离由近到远

2. 使用PostGIS扩展（推荐，适合大数据量） #

如果你的数据库安装了PostGIS空间扩展，用空间函数计算距离会更高效，尤其是数据量较大时（可通过空间索引优化查询性能）。

步骤1：确认PostGIS已启用 #

CREATE EXTENSION IF NOT EXISTS postgis;

步骤2：查询语句 #

假设你的表locations有geom字段（类型为GEOMETRY(Point, 4326)，4326是WGS84标准坐标系，对应GPS经纬度），目标点为ST_SetSRID(ST_MakePoint(-122.4194, 37.7749), 4326)：

SELECT
    id,
    name,
    ST_Y(geom) AS lat,
    ST_X(geom) AS lon,
    -- 计算距离（单位：米）
    ST_Distance(geom, ST_SetSRID(ST_MakePoint(-122.4194, 37.7749), 4326)) AS distance_m
FROM locations
ORDER BY distance_m ASC;

如果表中没有geom空间字段，也可以直接用经纬度生成点进行计算：

SELECT
    id,
    name,
    lat,
    lon,
    ST_Distance(
        ST_SetSRID(ST_MakePoint(lon, lat), 4326),
        ST_SetSRID(ST_MakePoint(-122.4194, 37.7749), 4326)
    ) AS distance_m
FROM locations
ORDER BY distance_m ASC;

关键说明： #

• ST_MakePoint的参数顺序是经度在前，纬度在后，注意不要混淆
• ST_SetSRID用于指定坐标系，4326是全球通用的GPS坐标系
• ST_Distance在地理坐标系下返回的距离单位是米，如需公里可除以1000
• 给geom字段创建空间索引可大幅提升大数据量查询速度：

  CREATE INDEX idx_locations_geom ON locations USING GIST(geom);
  

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