ST_Distance不能用于半径内分组，因其返回连续距离值而非离散分组键，导致每行距离几乎唯一，无法有效归并；且GROUP BY作用于标量计算结果，无法表达“某点周围r米内所有点”这一空间邻域关系。
直接说结论：ST_Distance 不能用于 GROUP BY 分组聚合的“半径内分组”，它只能算两点距离；真要按地理半径聚类，得用 ST_Within、ST_DWithin 或空间索引辅助的邻近查询，而不是靠 GROUP BY + ST_Distance。
为什么 ST_Distance 不能直接用于半径分组聚合
很多人写类似 GROUP BY ST_Distance(geom, 'POINT(x y)') 或试图在 HAVING 里写 ST_Distance(...) ，结果要么报错，要么逻辑完全不对——因为 <code>ST_Distance 返回的是浮点数（单位是坐标系单位，比如度），不是可分组的地理区域标识；而且它不带空间索引优化，暴力计算全表两两距离，性能崩得极快。
常见错误现象：
-
GROUP BY ST_Distance(a.geom, b.geom)报错：PostGIS 不允许在 GROUP BY 中使用非标量、非确定性或未索引的空间函数 - 写成
SELECT *, ST_Distance(geom, 'POINT(116.4 39.9)') AS d FROM poi HAVING d ：语法错误（HAVING 必须配合 GROUP BY），且没用上空间索引 - 即使加了
WHERE ST_Distance(geom, ...) ，也因无法走 <code>GIST索引而全表扫描，10 万条数据查 3 秒以上
正确做法：用 ST_DWithin 实现“某点半径内聚合”
ST_DWithin 是专为“半径内筛选”设计的函数，支持空间索引加速，语义清晰，且能自然嵌套进聚合流程。它判断两个几何是否在指定距离内（单位由 SRID 决定，如 EPSG:4326 下用度，EPSG:3857 或 EPSG:32650 下用米）。
典型场景：统计北京国贸 3 公里内的餐厅数量、平均评分
- 确保 geom 列有
GIST索引：CREATE INDEX idx_poi_geom ON poi USING GIST (geom); - 明确 SRID：如果用 WGS84（
EPSG:4326），距离单位是“度”，3 公里 ≈ 0.027 度（粗略值）；更推荐先转成米制坐标系（如EPSG:32650），再用整数米值 - 聚合写法示例（PostgreSQL + PostGIS）：
SELECT COUNT(*) AS cnt, AVG(rating) AS avg_rating FROM poi WHERE ST_DWithin( geom, ST_Transform(ST_SetSRID(ST_MakePoint(116.434, 39.909), 4326), 32650), 3000 -- 单位：米 );
注意：ST_DWithin 的第三个参数必须与前两个几何的 SRID 单位一致；强制转坐标系比硬凑度值更可靠。
想对多个中心点分别做半径聚合？用 LATERAL 或子查询
如果要对一批门店（centers 表）各自统计 500 米内用户数，不能简单 JOIN ... ON ST_Distance(...) ，必须用 <code>ST_DWithin + LATERAL 才能触发索引下推：
- 错误写法：
JOIN users u ON ST_Distance(c.geom, u.geom) → 全表笛卡尔积 + 暴力算距 - 正确写法（PostgreSQL 9.3+）：
SELECT c.id, c.name, COALESCE(t.cnt, 0) AS nearby_users FROM centers c LEFT JOIN LATERAL ( SELECT COUNT(*) AS cnt FROM users u WHERE ST_DWithin(u.geom, c.geom, 500) ) t ON true;
关键点：LATERAL 让子查询能引用外层 c.geom，且 PostGIS 会为每次迭代自动利用 u.geom 上的 GIST 索引快速找邻近点，效率提升百倍以上。
容易被忽略的坐标系和单位陷阱
最常踩的坑不是函数不会用，而是单位和 SRID 混乱：
- 用
ST_Distance或ST_DWithin前，务必确认两个几何的 SRID 一致；否则结果无意义，甚至报错 -
ST_Distance('POINT(0 0)', 'POINT(0 1)')在EPSG:4326下返回约 1.0（度），不代表 1 公里——实际约 111 公里 - 批量处理时别手写“0.01 度 ≈ 1 公里”换算，用
ST_Transform转到投影坐标系（如 UTM）再操作，数值直观、索引有效、精度可控 - PostGIS 3.0+ 支持
geography类型，可直接用米制距离（如ST_DWithin(geom::geography, 'POINT(...)::geography', 5000)），但要求原始列是geography或显式转换，且索引需建在geography列上
真正难的不是写对一行 SQL，而是让那一行在百万数据上跑得动、结果对得上、换个城市还能复用——坐标系选错，后面全白搭。
