PostgreSQL中实现Shoelace formula计算多边形面积（无PostGIS）

鞋带公式核心逻辑是通过累加相邻顶点的叉积总和，取绝对值除以2得到多边形面积，是有序顶点计算平面多边形面积的最优方案。

前提假设 #

你的有序顶点存储在表polygon_vertices中，表结构如下：

• polygon_id：多边形唯一ID，用于区分不同多边形
• point_order：顶点的排序序号，需按顺时针/逆时针连续排列
• x：转换完成的笛卡尔坐标x值
• y：转换完成的笛卡尔坐标y值

实现代码 #

WITH closed_vertices AS (
    SELECT
        polygon_id,
        x,
        y,
        -- 匹配当前顶点的下一个顶点，最后一个顶点的下一个取第一个顶点，完成多边形闭合
        LEAD(x, 1, FIRST_VALUE(x) OVER (PARTITION BY polygon_id ORDER BY point_order)) 
            OVER (PARTITION BY polygon_id ORDER BY point_order) AS next_x,
        LEAD(y, 1, FIRST_VALUE(y) OVER (PARTITION BY polygon_id ORDER BY point_order)) 
            OVER (PARTITION BY polygon_id ORDER BY point_order) AS next_y
    FROM polygon_vertices
),
cross_calculate AS (
    SELECT
        polygon_id,
        SUM(x * next_y - next_x * y) AS cross_total
    FROM closed_vertices
    GROUP BY polygon_id
)
SELECT
    polygon_id,
    -- 坐标单位为米时用下方语句，结果为平方公里
    ABS(cross_total) / 2 / 1000000 AS area_km2
    -- 坐标单位为公里时用下方语句，结果为平方公里
    -- ABS(cross_total) / 2 AS area_km2
FROM cross_calculate
-- 过滤顶点数不足3的无效多边形
WHERE polygon_id IN (
    SELECT polygon_id 
    FROM polygon_vertices 
    GROUP BY polygon_id 
    HAVING COUNT(point_order) >= 3
);

注意事项 #

• 顶点的point_order必须严格按多边形的顺时针/逆时针顺序连续编号，顺序错误会直接导致计算结果完全偏差
• 上述代码为标准SQL语法，兼容所有支持窗口函数的数据库（MySQL 8.0+、PostgreSQL、SQL Server等）
• 若你的数据库不支持窗口函数，可改用自关联的方式配对相邻顶点，计算逻辑完全一致

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