没问题，我来帮你搞定这个空间查询！既然你已经有了带经纬度和价格的房屋表，咱们结合PostGIS的特性来实现这个需求，步骤很清晰：

1. 确认PostGIS扩展已启用 #

首先得确保你的PostgreSQL数据库已经装上了PostGIS扩展，没装的话先执行这条命令：

CREATE EXTENSION IF NOT EXISTS postgis;

如果已经装过，这条命令也不会报错，放心执行就行。

2. 优化表结构（可选但推荐） #

假设你的表名叫houses，字段是id, latitude（纬度）, longitude（经度）, price（价格）。为了空间查询的性能，建议给表添加一个地理类型字段（GEOGRAPHY），用来存储经纬度对应的空间点：

ALTER TABLE houses ADD COLUMN geom GEOGRAPHY(POINT, 4326);
UPDATE houses SET geom = ST_SetSRID(ST_MakePoint(longitude, latitude), 4326);

这里4326是WGS84坐标系的SRID，也就是咱们常用的GPS经纬度坐标系。

3. 创建空间索引（关键！） #

3万条数据不算少，为了让查询速度飞起来，一定要给地理字段创建空间索引：

CREATE INDEX idx_houses_geom ON houses USING GIST(geom);

GIST索引是PostGIS里专门针对空间数据优化的索引类型，能大幅缩小查询范围。

4. 编写查询语句 #

现在就可以写核心的查询了！比如你要找经度116.403874、纬度39.914885周围5公里内的所有房屋，语句如下：

SELECT id, latitude, longitude, price
FROM houses
WHERE ST_DWithin(
    geom,
    ST_SetSRID(ST_MakePoint(116.403874, 39.914885), 4326)::GEOGRAPHY,
    5 * 1000  -- 5公里转成米，因为GEOGRAPHY类型的单位是米
);

语句解释： #

• ST_MakePoint(longitude, latitude)：把经纬度转换成PostGIS的点对象（注意顺序是经度在前，纬度在后，别搞反了！）
• ST_SetSRID(..., 4326)：给点对象指定坐标系为WGS84
• ::GEOGRAPHY：把点对象转换成地理类型（如果你的geom字段已经是GEOGRAPHY，这一步可以省略）
• ST_DWithin(a, b, distance)：判断空间对象a和b之间的距离是否小于等于指定的distance，这里用米做单位，所以X公里就乘以1000。

如果不想添加geom字段，直接用经纬度查询： #

如果你暂时不想修改表结构，也可以直接在查询里转换经纬度，不过性能会差一些（因为没法用索引）：

SELECT id, latitude, longitude, price
FROM houses
WHERE ST_DWithin(
    ST_SetSRID(ST_MakePoint(longitude, latitude), 4326)::GEOGRAPHY,
    ST_SetSRID(ST_MakePoint(116.403874, 39.914885), 4326)::GEOGRAPHY,
    X * 1000  -- 把X换成你要的公里数
);

还是建议先加geom字段和索引，毕竟3万条数据，带索引的查询速度会快很多。

额外提示 #

• 如果你用的是**几何类型（GEOMETRY）**而不是地理类型，那计算距离的时候需要考虑投影坐标系，单位可能是米或者度，不如GEOGRAPHY直接用米方便，所以优先推荐用GEOGRAPHY。
• 可以用ST_Distance函数查看具体的距离值，比如在SELECT里加上ST_Distance(geom, ST_SetSRID(ST_MakePoint(116.403874, 39.914885), 4326)::GEOGRAPHY)/1000 AS distance_km，就能看到每条记录离指定点的公里数。

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