4. 空间叠加分析基础与GIS操作框架
生态保护红线与敏感区的空间叠加，本质上是把两类或多类空间数据放到同一坐标参考下，通过几何求交、求并、求差或栅格代数运算，得到新的空间关系与属性组合。要保证叠加结果可信，先要弄清三件事：用哪种叠加类型、依据什么空间关系谓词判定、以及矢量与栅格各自适合什么条件。本节围绕这三点建立可操作的GIS框架。
一、叠加分析的四种基本类型
矢量叠加按运算语义可分为四类，它们决定了输出要素的几何与属性继承方式：
- 相交（Intersection）：保留两个图层公共部分，输出几何为交集，属性同时携带两图层的字段。用于提取“红线内且属于敏感区”的重叠斑块，是冲突识别最常用的运算。
- 联合（Union）：保留两图层全部范围，重叠处拆分并同时保留双方属性，非重叠处另一方属性为空。适合一次性生成全域分区底图，再按属性筛选。
- 擦除（Erase）：从输入图层中扣除擦除图层的覆盖范围。用于“红线扣除已建合法用地”或“敏感区扣除红线后剩余可开发空间”。
- 标识（Identity）：以输入图层为范围，叠加标识图层的属性但不裁剪输入几何。适合在保持红线斑块完整的前提下，标注每个斑块落入的敏感区类型。
选择依据是“输出要保留谁的几何、继承谁的属性”。冲突面积统计用相交；分区制图用联合；扣除合法用地用擦除；保持红线完整性做属性挂接用标识。
二、空间关系谓词：判定“是否叠加”的规则
叠加前必须先明确判定条件。OGC简单要素规范中的谓词是最通用的表达方式：Intersects（相交，含边界接触）、Within（完全位于内部）、Contains（包含）、Overlaps（部分重叠但不完全包含）、Touches（仅边界接触）、Crosses（穿越）、Disjoint（相离）。生态红线冲突识别通常用Intersects做初筛，再用面积阈值过滤边界接触产生的零面积或微小面积要素；若要求“红线斑块整体落入敏感区”，则改用Within。
谓词判定依赖几何有效性。实际数据常见自相交、悬挂节点、重复点等问题，需先执行修复几何（如ArcGIS的Repair Geometry、QGIS的“修复几何”），否则Intersects可能漏判或报错。此外，两图层必须统一到同一投影坐标系，跨带或经纬度数据直接叠加会产生面积失真。
三、矢量与栅格的适用条件
矢量叠加精度高、边界清晰、属性表灵活，适合红线边界、行政区、地块等具有明确边界的对象，但大数据量下求交计算量大，且叠加后易产生大量碎小斑块，需要按面积阈值合并。栅格叠加以像元为单元做代数运算，适合连续型变量（坡度、植被指数、生境适宜性）和多图层加权评价，运算快、易并行，但结果受分辨率制约，边界呈锯齿状，面积统计存在像元化误差。
判断原则：以“边界精确性”为首要需求选矢量；以“连续变量综合”或“大范围快速筛查”为首要需求选栅格；两者结合时，可先用栅格做适宜性初筛，再对高值区用矢量精确求交。
四、可执行的操作流程
以“红线与敏感区冲突识别”为例，给出一个可复现的矢量叠加流程（示例参数为教学假设，非实测结果）：
- 数据准备：红线图层
redline.shp、敏感区图层sensitive.shp，均转换为CGCS2000 3度带高斯投影，确保面积单位为平方米。 - 几何修复：对两图层执行修复几何，删除空几何与自相交要素，记录修复前后要素数。
- 相交运算：以红线为输入、敏感区为叠加层执行Intersection，输出
overlap.shp，属性表同时含红线编号与敏感区类型。 - 面积计算与过滤：新增字段
area_m2，计算几何面积；按area_m2 >= 100过滤，剔除边界接触产生的微小碎斑。 - 汇总统计：按敏感区类型分组求和面积，得到各类敏感区被红线覆盖的面积清单。
- 质量检查：核对叠加前后总面积守恒关系（红线面积 = 重叠面积 + 未重叠面积），偏差超过0.5%时回查几何与投影。
掌握叠加类型、谓词与矢量/栅格适用条件后，后续的红线冲突识别、敏感区准入判断与分区制图，都可以在这套框架上按需组合运算，而不必每次重新设计流程。