叠加分析从几何交集生成新关系

叠加分析把两组空间对象放在同一坐标空间中，通过相交、差集或联合生成新几何，并组合相关属性。它不仅回答“哪些对象相交”，还可以产生相交部分本身，用于计算受影响面积、切分管理单元或形成新的空间分区。

选择相交对象与生成交集不同

空间谓词只返回关系真假。查询与保护区相交的地块，可以保留完整地块几何；几何交集则返回地块落在保护区内部的部分。两者用于不同统计。

若用完整地块面积统计受保护面积，会高估跨边界地块。若只需列出受影响地块，生成交集又会增加计算和碎片。输出目标应在选择函数前确定。

线与面相交可能返回线，面与面相交可能返回面、线或点，取决于实际接触。下游不能假设结果始终与输入维度相同。

交集保留两个几何共同部分，差集保留 A 中不属于 B 的部分，联合形成两者覆盖的组合。对称差保留只属于其中一方的区域。裁剪通常用一个范围几何限制另一数据集。

操作方向很重要。A 减 B 与 B 减 A 不同；属性传播也取决于哪个数据集作为主体。工作流应明确输出记录代表哪一类业务对象。

“union”在不同软件中可能指几何合并，也可能指保留两层所有分区并组合属性。使用函数名称前要确认其几何与表结构语义。

参与叠加的数据必须处于可比较参考系。临时投影转换应选择适合区域和面积、距离需求的操作，并记录精度。错误标注 CRS 会产生看似正常但位置错误的结果。

两套边界来自不同尺度或年代时，理论上重合的线可能略有偏移，生成狭长缝隙和细碎多边形。它不是运算错误，而是输入几何不一致的显现。

精度网格或坐标吸附可以减少数值噪声，但会移动顶点并改变面积。容差应依据数据精度设定，应用前后要统计变化，不能任意四舍五入坐标。

自相交多边形、未闭合环或错误嵌套可能让叠加失败或产生异常结果。处理前应检查有效性，并记录修复方法。自动修复可能把一个对象拆成多个部分或删除退化部分。

空几何与 NULL 不同。两个对象没有共同部分时，结果可能为空；输入缺失则表示没有可计算几何。统计时要分别处理。

GeometryCollection 和 MultiPolygon 需要下游支持。只提取面会丢掉边界接触形成的线和点，这可能符合面积任务，但必须是明确选择。

一个地块与多个区域相交时，会形成多条输出记录，每条携带双方属性。面积或人口等原对象属性若原样复制，再直接求和会重复计算。

需要按交集面积分摊时，应先确认变量能否按面积线性分配。人口、产量和风险不一定均匀分布，面积比例只是模型假设。

字段同名、类型不同和编码域冲突也要处理。输出 schema 应明确来源和单位，而不是让软件自动加后缀后无人理解。

复杂边界叠加会产生大量小碎片。部分是真实空间分区，部分来自边界不一致或浮点误差。按面积阈值删除碎片可能同时删除真实小对象。

可以先统一边界、简化到符合任务的尺度或按网格分块处理，但简化会改变拓扑，分块需要在接缝处去重。任何优化都要保持业务可接受的几何误差。

数据库中应先用空间索引筛选候选，再计算精确交集。将全部对象两两求交会产生巨大的无效组合。执行计划、候选数和输出顶点数都应监控。

栅格叠加要求网格对齐或先重采样。坐标参考、分辨率、原点不同的栅格即使覆盖同一区域，也不能直接逐像元运算。

分类栅格组合产生类别交叉表，连续栅格则可能执行代数计算。无数据传播规则必须明确，否则一个输入缺失可能被错误当作零。

矢量与栅格叠加需要决定按像元中心、覆盖比例还是触碰关系选取。不同规则会影响边界像元和小对象统计。

可用面积守恒检查结果：在无重叠、无空洞的理想分区中，交集部分面积之和应接近原对象面积。差异需要解释投影、精度和无效几何处理。

还应统计空结果、多部件、碎片数量、最小面积和每个输入的匹配数，并在边界复杂区域抽样制图。仅检查程序没有报错不足以证明结果正确。

叠加分析把几何关系转化为新数据。可靠输出取决于明确的集合操作、输入精度、属性传播和统计假设，而不只是调用一个交集函数。

结果还应保留输入数据版本、处理参数和生成时间。上游边界更新后，同名统计可能发生变化；没有可追溯记录，就无法判断差异来自现实变化、数据修订还是叠加规则调整。

PostGIS Project：PostGIS Official Manual，https://www.postgis.net/documentation/manual/
Open Geospatial Consortium：OGC Simple Feature Access，https://www.ogc.org/standards/sfa/