如何在PostGIS中统计线段表的唯一几何图形数量？

针对3200万行的线段表，以下是几种高效统计唯一几何数量的方法，按需选择：

基础精确统计（不考虑线段方向） #

如果仅需严格匹配几何数据（包括线段方向），通过几何的二进制表示分组统计，性能优于直接使用COUNT(DISTINCT geom)：

SELECT COUNT(*)
FROM (
    SELECT 1
    FROM your_table_name
    GROUP BY ST_AsBinary(geom)
) AS unique_geoms;

统计时忽略线段方向 #

若认为反向线段（如LINESTRING(0 0,1 1)和LINESTRING(1 1,0 0)）属于重复，可先统一线段方向再分组：

SELECT COUNT(*)
FROM (
    SELECT 1
    FROM your_table_name
    GROUP BY 
        CASE 
            WHEN ST_StartPoint(geom) < ST_EndPoint(geom) 
            THEN geom 
            ELSE ST_Reverse(geom) 
        END
) AS unique_geoms;

性能优化建议 #

1. 临时创建几何索引：若表中无geom字段的GiST索引，可临时创建以加速分组查询（注意：3200万行建索引需消耗一定时间和磁盘空间）：

CREATE INDEX idx_temp_your_table_geom ON your_table_name USING GIST(geom);

统计完成后可删除索引：

DROP INDEX idx_temp_your_table_geom;

2. 启用并行查询：确保PostgreSQL配置中max_parallel_workers_per_gather参数设置合理（如设为4或更高），让数据库利用多核加速大表查询。

3. 哈希分组替代：若追求更快的分组速度，可使用几何二进制的MD5哈希值分组（存在极小碰撞概率，需权衡）：

SELECT COUNT(*)
FROM (
    SELECT 1
    FROM your_table_name
    GROUP BY md5(ST_AsBinary(geom))
) AS unique_geoms;

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