PostGIS空间汇总函数如何实现区域数据聚合？关键参数与优化技巧详解（附：实战代码）

wowwwai GIS研习社 · 工具流程与项目排障

在实际项目中，PostGIS空间汇总函数如何实现区域数据聚合？关键参数与优化技巧详解（附：实战代码）这个问题通常出现在“点数据很多、行政区很多、统计结果很慢或不准”的场景里。比如把门店、人口点、污染监测点、订单位置按街道、区县或网格汇总，输出每个区域内的数量、总和、均值、最大值等指标。

本文以PostGIS空间汇总函数为主线，讲清楚区域数据聚合的常用写法、关键空间关系函数、索引优化方法，以及如何检查结果是否可靠。示例使用行政区面数据和业务点数据，适合 GIS 工程师、空间数据分析师和 WebGIS 后端开发者参考。

引言：PostGIS空间汇总函数解决什么问题

所谓PostGIS空间汇总函数，并不是某一个单独函数，而是一组空间关系函数、空间计算函数与 PostgreSQL 聚合函数组合起来的查询模式。常见组合包括：

• ST_Intersects：判断两个几何是否相交，常用于点落入面、线穿过面、面与面叠加。
• ST_Contains：判断一个面是否包含另一个几何，常用于行政区包含点。
• ST_Within：判断一个几何是否位于另一个几何内部，语义上与 ST_Contains 方向相反。
• COUNT、SUM、AVG、MIN、MAX：PostgreSQL 原生聚合函数，用于数量、求和、均值等统计。
• GROUP BY：按区域编号、区域名称或网格 ID 分组。

一个典型问题是：已有一张区县面表和一张企业点表，如何统计每个区县内企业数量、注册资本总额和平均注册资本？这就是典型的PostGIS区域数据聚合。

背景：为什么区域数据聚合容易出错

很多人写 PostGIS 空间汇总 SQL 时，能跑出结果，但结果经常存在三个问题：慢、不准、重复。原因通常不是 PostGIS 本身，而是数据、坐标系、空间关系选择和索引使用不当。

常见业务场景

• 按行政区统计 POI 数量。
• 按街道汇总人口点、企业点、案件点。
• 按网格统计订单量、设备数量、风险点数量。
• 按流域或保护区汇总监测站点指标。
• 按缓冲区统计周边设施数量。

常见错误表现

• 查询执行十几分钟甚至更久。
• 边界上的点没有被统计，或被多个区域重复统计。
• 同一个点因为行政区面重叠被重复计算。
• 经纬度数据直接做面积或距离，结果单位不符合预期。
• SQL 能运行，但没有使用空间索引。

区域数据聚合的核心不是“把两个表 JOIN 一下”，而是先明确空间关系，再保证坐标系一致、几何有效、索引可用，最后才是聚合统计。

原理：PostGIS区域数据聚合的核心逻辑

PostGIS空间汇总函数的基本逻辑可以拆成三步：先建立空间匹配关系，再按区域字段分组，最后使用聚合函数计算指标。

基础 SQL 结构

SELECT
    a.region_id,
    a.region_name,
    COUNT(p.id) AS point_count,
    SUM(p.amount) AS amount_sum,
    AVG(p.amount) AS amount_avg
FROM admin_region a
LEFT JOIN business_point p
    ON ST_Intersects(a.geom, p.geom)
GROUP BY
    a.region_id,
    a.region_name
ORDER BY
    a.region_id;

这里的 admin_region 是区域面表，business_point 是业务点表，geom 是几何字段。LEFT JOIN 可以保留没有点落入的区域，统计结果中的数量为 0 或金额为空。

ST_Intersects、ST_Contains、ST_Within怎么选

函数	适用场景	注意事项
ST_Intersects(a.geom, p.geom)	点、线、面只要有交集就匹配	边界点也可能被认为相交，适合多数空间连接
ST_Contains(a.geom, p.geom)	区域面包含点或小面	边界上的点通常不算严格包含
ST_Within(p.geom, a.geom)	点位于区域面内部	语义清晰，但同样要注意边界问题
ST_Covers(a.geom, p.geom)	区域覆盖点，包括边界	处理边界点时比 ST_Contains 更稳妥

如果你的业务要求“边界上的点也算入该区域”，可以优先考虑 ST_Covers 或 ST_Intersects。如果区域之间存在重叠，则需要先处理面数据，否则一个点可能被多个区域统计。

空间索引为什么重要

PostGIS 的空间索引通常使用 GiST 索引。没有索引时，数据库可能需要把每个区域和每个点逐一比较，数据量稍大就会非常慢。有了空间索引后，PostGIS 可以先用外接矩形快速过滤候选对象，再做精确空间判断。

CREATE INDEX idx_admin_region_geom
ON admin_region
USING GIST (geom);

CREATE INDEX idx_business_point_geom
ON business_point
USING GIST (geom);

ANALYZE admin_region;
ANALYZE business_point;

ANALYZE 用于更新统计信息，让 PostgreSQL 查询优化器更准确地选择执行计划。很多PostGIS空间汇总函数查询慢的问题，实际是索引建了但统计信息没有更新。

步骤：用PostGIS空间汇总函数实现区域数据聚合

步骤一：检查表结构与坐标系

先确认两个表都有几何字段，并且空间参考系统一致。

SELECT
    Find_SRID('public', 'admin_region', 'geom') AS region_srid,
    Find_SRID('public', 'business_point', 'geom') AS point_srid;

如果一个是 4326，另一个是 3857，直接空间汇总可能出现匹配错误。应先统一坐标系。

ALTER TABLE business_point
ADD COLUMN geom_4326 geometry(Point, 4326);

UPDATE business_point
SET geom_4326 = ST_Transform(geom, 4326);

如果原始数据本来就是经纬度，但 SRID 缺失，不能盲目使用 ST_Transform，应先使用 ST_SetSRID 正确声明坐标系。

UPDATE business_point
SET geom = ST_SetSRID(geom, 4326)
WHERE ST_SRID(geom) = 0;

步骤二：检查几何有效性

区域面如果存在自相交、空几何或非法环，空间判断可能不稳定。建议先检查。

SELECT
    region_id,
    region_name,
    ST_IsValid(geom) AS is_valid,
    ST_IsValidReason(geom) AS valid_reason
FROM admin_region
WHERE NOT ST_IsValid(geom);

对于可修复的面数据，可以使用 ST_MakeValid 生成有效几何。但生产环境中建议先备份原表。

UPDATE admin_region
SET geom = ST_MakeValid(geom)
WHERE NOT ST_IsValid(geom);

步骤三：建立空间索引

对参与空间关系判断的几何字段建立 GiST 索引。

CREATE INDEX IF NOT EXISTS idx_admin_region_geom
ON admin_region
USING GIST (geom);

CREATE INDEX IF NOT EXISTS idx_business_point_geom
ON business_point
USING GIST (geom);

ANALYZE admin_region;
ANALYZE business_point;

如果使用的是转换后的几何字段，例如 geom_4326，索引也必须建在实际参与查询的字段上。

CREATE INDEX IF NOT EXISTS idx_business_point_geom_4326
ON business_point
USING GIST (geom_4326);

步骤四：编写基础区域统计 SQL

以下示例统计每个行政区内的企业数量、注册资本总额和平均值。

SELECT
    a.region_id,
    a.region_name,
    COUNT(p.id) AS company_count,
    COALESCE(SUM(p.capital), 0) AS capital_sum,
    ROUND(AVG(p.capital)::numeric, 2) AS capital_avg
FROM admin_region a
LEFT JOIN business_point p
    ON ST_Covers(a.geom, p.geom)
GROUP BY
    a.region_id,
    a.region_name
ORDER BY
    company_count DESC;

这里使用 ST_Covers，原因是它对边界点更友好。如果业务明确要求点必须严格位于区域内部，可以改为 ST_Contains。

步骤五：把汇总结果写入新表

如果统计结果需要给 WebGIS 或报表系统使用，建议生成一张结果表。

DROP TABLE IF EXISTS region_company_summary;

CREATE TABLE region_company_summary AS
SELECT
    a.region_id,
    a.region_name,
    a.geom,
    COUNT(p.id) AS company_count,
    COALESCE(SUM(p.capital), 0) AS capital_sum,
    ROUND(AVG(p.capital)::numeric, 2) AS capital_avg
FROM admin_region a
LEFT JOIN business_point p
    ON ST_Covers(a.geom, p.geom)
GROUP BY
    a.region_id,
    a.region_name,
    a.geom;

结果表保留区域几何字段后，可以直接用于 QGIS、ArcGIS Pro 或 GeoServer 发布专题图。

CREATE INDEX idx_region_company_summary_geom
ON region_company_summary
USING GIST (geom);

ANALYZE region_company_summary;

步骤六：使用EXPLAIN检查是否走索引

优化PostGIS空间汇总函数时，不要只看 SQL 写法，还要看执行计划。

EXPLAIN ANALYZE
SELECT
    a.region_id,
    COUNT(p.id) AS point_count
FROM admin_region a
LEFT JOIN business_point p
    ON ST_Intersects(a.geom, p.geom)
GROUP BY a.region_id;

如果执行计划中出现 Index Scan、Bitmap Index Scan 或与 GiST 索引相关的信息，通常说明空间索引被使用了。如果看到大量 Seq Scan，需要检查索引、统计信息、字段类型和 SQL 条件。

常见坑：PostGIS空间汇总函数结果不准的原因

坑一：坐标系不一致

空间关系判断要求几何坐标在同一空间参考下。如果区域是 EPSG:4326，点是 EPSG:3857，坐标数值完全不在同一体系里，结果可能为空或严重错误。

• 用 ST_SRID 或 Find_SRID 检查 SRID。
• 用 ST_Transform 做真实坐标转换。
• 用 ST_SetSRID 只声明坐标系，不改变坐标值。

坑二：边界点统计规则不清

点刚好落在行政区边界上时，ST_Contains 可能不统计它。对于行政区边界、网格边界、缓冲区边界，必须提前定义业务规则。

• 边界点要统计：考虑 ST_Covers 或 ST_Intersects。
• 边界点不能重复统计：需要行政区拓扑无重叠，并增加去重规则。
• 点落在多个区域：需要按优先级、面积占比或最近中心点规则处理。

坑三：面数据重叠导致重复汇总

如果行政区面之间有重叠，一个点可能同时匹配多个区域，导致数量和金额被重复计算。可以先检查区域之间是否重叠。

SELECT
    a.region_id AS region_a,
    b.region_id AS region_b
FROM admin_region a
JOIN admin_region b
    ON a.region_id < b.region_id
   AND ST_Intersects(a.geom, b.geom)
   AND ST_Area(ST_Intersection(a.geom, b.geom)) > 0;

如果存在重叠，需要回到数据生产环节修复边界，或者在统计逻辑中明确主归属区域。

坑四：COUNT字段选择不当

COUNT(*) 在 LEFT JOIN 场景下可能让没有点的区域也统计为 1。区域聚合时通常应使用 COUNT(p.id)。

SELECT
    a.region_id,
    COUNT(*) AS wrong_count,
    COUNT(p.id) AS correct_count
FROM admin_region a
LEFT JOIN business_point p
    ON ST_Covers(a.geom, p.geom)
GROUP BY a.region_id;

坑五：在JOIN条件里直接ST_Transform导致索引失效

下面这种写法虽然方便，但可能导致已有空间索引无法有效使用。

SELECT
    a.region_id,
    COUNT(p.id)
FROM admin_region a
JOIN business_point p
    ON ST_Intersects(a.geom, ST_Transform(p.geom, 4326))
GROUP BY a.region_id;

更推荐提前生成统一坐标系的几何字段，并对该字段建立空间索引。

ALTER TABLE business_point
ADD COLUMN IF NOT EXISTS geom_4326 geometry(Point, 4326);

UPDATE business_point
SET geom_4326 = ST_Transform(geom, 4326)
WHERE geom_4326 IS NULL;

CREATE INDEX IF NOT EXISTS idx_business_point_geom_4326
ON business_point
USING GIST (geom_4326);

方法比较：不同PostGIS区域数据聚合写法怎么选

方法	示例	优点	适用场景
点面空间连接后 GROUP BY	ST_Covers(a.geom, p.geom)	直观、通用、易维护	按行政区、街道、网格统计点数据
相交后按面积加权	ST_Intersection + ST_Area	适合面与面叠加统计	土地利用、人口栅格转面、规划区叠加
先过滤外接矩形再精确判断	&& + ST_Intersects	可增强索引过滤意识	大数据量空间查询优化
预计算归属区域 ID	点表增加 region_id	后续统计非常快	数据更新频率低、查询频率高的业务系统
物化视图	CREATE MATERIALIZED VIEW	适合报表与看板	定时更新的区域统计结果

面与面叠加汇总示例

如果不是点落入面，而是“地类图斑按行政区汇总面积”，就不能简单 COUNT。通常需要相交切割后按面积统计。

SELECT
    a.region_id,
    a.region_name,
    l.land_type,
    SUM(
        ST_Area(
            ST_Intersection(a.geom, l.geom)
        )
    ) AS area_sum
FROM admin_region a
JOIN landuse_polygon l
    ON ST_Intersects(a.geom, l.geom)
GROUP BY
    a.region_id,
    a.region_name,
    l.land_type;

这里要特别注意坐标系。如果使用 EPSG:4326 经纬度坐标，ST_Area 计算结果不是常见的平方米。应转换到合适的投影坐标系后再计算面积。

SELECT
    a.region_id,
    l.land_type,
    SUM(
        ST_Area(
            ST_Intersection(
                ST_Transform(a.geom, 4547),
                ST_Transform(l.geom, 4547)
            )
        )
    ) AS area_m2
FROM admin_region a
JOIN landuse_polygon l
    ON ST_Intersects(a.geom, l.geom)
GROUP BY
    a.region_id,
    l.land_type;

如果数据量很大，不建议在高频查询中反复 ST_Transform 和 ST_Intersection，可以提前生成投影字段或中间结果表。

检查清单：上线前如何验证PostGIS空间汇总结果

• 确认区域表和业务表的 SRID 是否一致。
• 确认几何字段是否为空，是否存在无效几何。
• 确认空间索引是否建立在实际查询字段上。
• 使用 EXPLAIN ANALYZE 检查查询计划。
• 确认边界点采用 ST_Contains、ST_Covers 还是 ST_Intersects。
• 确认区域面之间是否存在重叠。
• 确认 COUNT(p.id) 是否比 COUNT(*) 更符合业务含义。
• 抽取几个区域在 QGIS 或 ArcGIS Pro 中人工核对。
• 对于金额、人口等数值字段，检查 NULL 值是否需要用 COALESCE 处理。
• 对于面积统计，确认是否使用了合适的投影坐标系。

推荐的结果核对 SQL

可以先核对全表点数量与各区域汇总数量是否一致。前提是区域不重叠，且每个点最多落入一个区域。

SELECT COUNT(*) AS total_points
FROM business_point;

SELECT SUM(company_count) AS summary_points
FROM region_company_summary;

如果两个结果差异明显，需要检查是否有点落在区域外，或者区域之间存在重叠。

SELECT COUNT(p.id) AS outside_points
FROM business_point p
LEFT JOIN admin_region a
    ON ST_Covers(a.geom, p.geom)
WHERE a.region_id IS NULL;

FAQ：PostGIS空间汇总函数常见问题

PostGIS空间汇总函数是不是只有ST_SummaryStats？

不是。ST_SummaryStats 主要用于栅格统计，而本文讨论的PostGIS空间汇总函数是矢量区域聚合场景，通常由 ST_Intersects、ST_Covers、GROUP BY、COUNT、SUM 等组合实现。

PostGIS区域数据聚合应该用ST_Contains还是ST_Intersects？

点面统计中，如果边界点也要算入区域，优先考虑 ST_Covers 或 ST_Intersects。如果要求点必须严格在面内部，可以使用 ST_Contains。实际项目中要先确认业务口径。

为什么PostGIS空间汇总查询很慢？

常见原因包括没有 GiST 空间索引、索引建在错误字段上、查询中直接使用 ST_Transform 导致索引难以利用、几何过于复杂、统计信息未更新。建议先建索引，再执行 ANALYZE，最后用 EXPLAIN ANALYZE 检查。

没有点的区域如何显示为0？

使用 LEFT JOIN 保留所有区域，并对数值汇总结果使用 COALESCE。数量字段建议使用 COUNT(p.id)。

SELECT
    a.region_id,
    COUNT(p.id) AS point_count,
    COALESCE(SUM(p.amount), 0) AS amount_sum
FROM admin_region a
LEFT JOIN business_point p
    ON ST_Covers(a.geom, p.geom)
GROUP BY a.region_id;

PostGIS面与面区域汇总能直接SUM面积吗？

不能简单直接汇总原图斑面积。面与面叠加时，应先用 ST_Intersection 得到相交部分，再用 ST_Area 计算相交面积，并确保坐标系适合面积计算。

区域统计结果如何用于QGIS或WebGIS？

可以把汇总结果写入包含几何字段的新表，或者创建物化视图。QGIS 可以直接连接 PostGIS 表进行分级设色；GeoServer、MapServer 或后端 API 也可以基于结果表发布服务。

结论：稳定的PostGIS区域聚合要同时关注SQL、数据和索引

PostGIS空间汇总函数实现区域数据聚合的核心写法并不复杂：区域表与业务表通过空间关系函数连接，再用 GROUP BY 和聚合函数统计。但要让结果准确、查询高效，必须同时处理坐标系、几何有效性、边界规则、重复统计和空间索引。

在生产项目中，建议按“检查坐标系、修复几何、建立索引、选择空间关系、编写汇总 SQL、检查执行计划、抽样核对结果”的顺序执行。这样写出来的PostGIS区域数据聚合不仅能跑通，也更容易维护、复用和上线。