PostGIS计算棱锥体积的实现方法及报错解决方案

你遇到的XX000错误核心有两个问题：

• 构造的POLYHEDRALSURFACE面顶点环绕方向不一致，所有面的法向量没有统一朝向实体外部，导致ST_MakeSolid将其识别为未连通的面集合，无法生成有效实体
• 原始SQL未将几何字符串显式转换为geometry类型，存在隐式类型转换的兼容性风险

调整每个面的顶点顺序保证法向量统一朝外，加上显式几何类型转换，修改后的可运行SQL如下：

SELECT 
    ST_Volume(ST_MakeSolid(geom)) AS pyramid_volume
FROM (
    SELECT ST_GeomFromText('POLYHEDRALSURFACE(
        ((0 0 0, 1 0 0, 0 1 0, 0 0 0)),
        ((0 0 0, 0 1 0, 0 0 1, 0 0 0)),
        ((0 0 0, 0 0 1, 1 0 0, 0 0 0)),
        ((1 0 0, 0 1 0, 0 0 1, 1 0 0))
    )') As f(geom);

执行后返回的结果为0.166666666666667，也就是该三棱锥的体积1/6，符合预期。

• 构造用于生成实体的多面曲面时，所有面的顶点顺序需要遵循「从实体外部观察为逆时针排列」的规则，保证法向量统一朝外
• 计算体积前可先用ST_IsSolid(ST_MakeSolid(geom))验证生成的实体是否有效，提前排查拓扑问题
• 复杂多面体可调用ST_OrientateSurface函数统一调整所有面的法向方向，无需手动调整每个面的顶点顺序

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