PostGIS求交集怎么写?空间分析函数详解?

“为什么我的PostGIS交集查询返回空？”——你不是一个人

上周一位在国土规划院实习的研究生私信我：“Dr. Gis，我用ST_Intersection查两个面图层的重叠区域，结果全是空，是不是数据坏了？”——其实90%的情况，问题不在数据，而在你没搞懂“空间谓词”和“空间操作”的区别。今天我就带你从原理到实战，彻底吃透PostGIS的空间交集运算。

空间交集 ≠ 简单重叠：先搞懂“拓扑关系”这层窗户纸

很多初学者以为“求交集”就是找两个图形重叠的部分，这没错，但PostGIS的底层逻辑比你想得更严谨。它遵循的是DE-9IM（Dimensionally Extended 9-Intersection Model）模型——听名字很吓人，其实你可以把它想象成“地理版的布尔代数”。

我在参与某市生态红线划定项目时，就吃过这个亏：两个面看起来完全重合，但因为坐标精度差了0.0001度，ST_Intersection直接返回空。后来改用ST_Intersects先做预筛选，再缓冲0.1米处理边界，才解决问题。

简单类比：求交集就像用两把剪刀同时剪一张纸——只有真正重合的那部分纸片才会留下来。如果两把剪刀根本没碰到一起（不相交），或者只擦了个边（边界接触但无面积重叠），结果自然就是空。

三大核心函数实战：ST_Intersection、ST_Intersects、ST_Overlaps怎么选？

PostGIS里和“交集”相关的函数有好几个，新手最容易混淆。下面这张表帮你理清它们的核心差异：

记住口诀：先问“有没有”，再算“是什么”。也就是说，先用ST_Intersects过滤出真正相交的要素，再对它们执行ST_Intersection，能极大提升查询效率。

手把手实战：从建表到输出GeoJSON的完整流程

假设我们有两个面图层：一个是“城市建成区”(urban_area)，另一个是“生态保护区”(eco_zone)。目标是找出所有与保护区重叠的建成区，并计算重叠面积。

-- 步骤1：创建示例表（实际项目中你已有数据）
CREATE TABLE urban_area (
    id SERIAL PRIMARY KEY,
    name VARCHAR(50),
    geom GEOMETRY(Polygon, 4326)
);

CREATE TABLE eco_zone (
    id SERIAL PRIMARY KEY,
    zone_name VARCHAR(50),
    geom GEOMETRY(Polygon, 4326)
);

-- 步骤2：插入模拟数据（略，假设已存在）

-- 步骤3：核心查询 - 求交集并计算面积
SELECT 
    u.name AS urban_name,
    e.zone_name,
    ST_Area(ST_Intersection(u.geom, e.geom)::geography) / 1000000 AS overlap_sq_km, -- 转为平方公里
    ST_AsGeoJSON(ST_Intersection(u.geom, e.geom)) AS overlap_geom
FROM urban_area u
JOIN eco_zone e 
ON ST_Intersects(u.geom, e.geom)  -- 先用谓词加速
WHERE ST_Intersection(u.geom, e.geom) IS NOT NULL  -- 双保险防空几何
AND ST_GeometryType(ST_Intersection(u.geom, e.geom)) = 'ST_Polygon'; -- 确保结果是面（排除线/点）

关键细节解释：

• ST_Area(...::geography)：将几何转为地理类型计算真实面积（单位平方米），再除以10⁶得平方公里。
• ST_AsGeoJSON()：直接输出标准GeoJSON，方便前端可视化。
• ST_GeometryType() = 'ST_Polygon'：过滤掉退化成线或点的无效交集（比如两个面仅边界相切）。

避坑指南：三个高频报错的终极解决方案

1. 错误：ERROR: BOOM! Could not generate outside point!
   原因：几何自相交或拓扑异常。解决：用ST_MakeValid(geom)修复几何。

2. 结果为空但肉眼可见重叠？
   原因：坐标系不一致或精度误差。解决：统一SRID + 使用ST_Buffer(geom, 0.00001)微调边界。

3. 性能极慢？
   原因：未建空间索引。解决：对geom字段创建GIST索引：CREATE INDEX idx_urban_geom ON urban_area USING GIST(geom);

总结：交集运算的“三要三不要”

✅ 要先用ST_Intersects做预筛选
✅ 要检查几何有效性（ST_IsValid）
✅ 要统一坐标系并建空间索引
❌ 不要直接对大数据集跑ST_Intersection
❌ 不要忽略结果几何类型（可能是点/线/面）
❌ 不要忘记面积单位转换（degree² ≠ m²）

空间分析的魅力就在于：看似简单的“求交集”，背后藏着拓扑学、计算几何和工程实践的层层智慧。你现在踩过的坑，都是未来简历上的亮点。

你在项目中用PostGIS求交集时遇到过什么奇葩问题？评论区留下你的故事，点赞最高的送《PostGIS实战手册》电子版！