执行st_union计算无重叠线长时出现GEOS错误如何解决？

解决方案 #

方案1：实现报错自动回退逻辑 #

SQL原生语法不支持直接捕获运行时错误，你可以通过自定义PL/pgSQL函数捕获st_union抛出的拓扑异常，异常触发时自动回退到求和逻辑：

CREATE OR REPLACE FUNCTION calc_non_overlap_length(name_val text)
RETURNS float8 AS $$
DECLARE
    result_len float8;
BEGIN
    -- 优先执行去重统计逻辑
    SELECT st_length(st_union(t.geom)) INTO result_len
    FROM mytable t
    WHERE t.name = name_val AND st_isvalid(t.geom);
    RETURN result_len;
EXCEPTION
    -- 捕获拓扑异常等所有运行时错误，回退到普通求和
    WHEN OTHERS THEN
        SELECT sum(st_length(t.geom)) INTO result_len
        FROM mytable t
        WHERE t.name = name_val AND st_isvalid(t.geom);
        RETURN result_len;
END;
$$ LANGUAGE plpgsql STABLE;

函数创建完成后直接调用即可生成视图：

CREATE OR REPLACE VIEW line_total_length AS
SELECT name, calc_non_overlap_length(name) AS total_length
FROM mytable
GROUP BY name;

方案2：优先修复几何精度问题避免报错（更推荐） #

你遇到的拓扑异常90%以上是浮点精度误差导致的微小交叉未被结点化，不用回退逻辑，只要先对齐坐标精度就可以正常执行st_union，结果也更准确：

SELECT 
    name,
    st_length(st_union(ST_SnapToGrid(t.geom, 0.001))) AS total_length
FROM mytable t 
WHERE st_isvalid(t.geom) 
GROUP BY name;

上面语句里的0.001是坐标对齐精度，单位和你的空间参考单位一致，可根据你的数据精度要求调整（比如米制单位下0.001就是对齐到毫米，完全不影响统计精度），绝大多数场景下都可以解决拓扑异常问题。

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