Postgres空间数据操作：如何查找1公里范围内的最大点簇

前置依赖 #

确保你的PostgreSQL数据库已经安装并启用PostGIS空间扩展：

CREATE EXTENSION IF NOT EXISTS postgis;

方案1：使用DBSCAN密度聚类（推荐，PostGIS 2.3及以上版本支持） #

DBSCAN是基于密度的聚类算法，刚好适配「1公里范围内的点簇」需求，可直接指定距离阈值识别符合要求的点簇，统计后取点数最多的即可。
示例代码如下：

WITH clustered_points AS (
    SELECT 
        对象id,
        坐标字段,
        -- eps设置为1000即1000米/1公里，minpoints为簇最少包含的点数，可按需调整
        ST_ClusterDBSCAN(坐标字段::geography, eps := 1000, minpoints := 2) OVER () AS cluster_id
    FROM 你的数据表名
),
cluster_stat AS (
    SELECT
        cluster_id,
        COUNT(*) AS point_total,
        ST_Collect(坐标字段) AS cluster_geometry
    FROM clustered_points
    WHERE cluster_id IS NOT NULL
    GROUP BY cluster_id
)
-- 按簇包含的点数倒序，取第一条即为最大点簇
SELECT * FROM cluster_stat ORDER BY point_total DESC LIMIT 1;

注意事项 #

• 上述代码中直接将geometry转geography类型计算球面距离，适配默认WGS84（EPSG:4326）坐标系的坐标，结果精度更高
• 如果你的坐标字段已经是单位为米的投影坐标系（比如UTM、Web墨卡托EPSG:3857），可以去掉::geography强转，直接传入坐标字段即可
• minpoints参数可根据业务调整：如果要求点簇至少包含3个点，就将参数值改为3，小于该数量的点集合不会被识别为有效簇

方案2：使用ST_ClusterWithin（适配低版本PostGIS） #

如果你使用的PostGIS版本低于2.3，没有ST_ClusterDBSCAN函数，可以用ST_ClusterWithin实现：

WITH clusters AS (
    -- 1000为距离阈值，单位米，需对应坐标系单位
    SELECT ST_ClusterWithin(坐标字段::geography, 1000) AS cluster_list
    FROM 你的数据表名
),
cluster_unnest AS (
    SELECT
        unnest(cluster_list) AS cluster_geometry,
        ST_NumGeometries(unnest(cluster_list)) AS point_total
    FROM clusters
)
SELECT * FROM cluster_unnest ORDER BY point_total DESC LIMIT 1;

内容的提问来源于stack exchange，提问作者liza.sap