3. 空间数据基础与坐标系统

断面位置、汇水边界、河网拓扑这三类对象，在GIS里往往来自不同渠道：断面坐标可能来自手持GPS实测点，矢量边界来自水利部门下发的行政或流域界线，水系数据来自测绘部门的基础地理信息库，DEM则可能是30米或12.5米分辨率的栅格产品。它们各自的坐标系、高程基准、精度水平都不相同。如果不在分析前把它们统一到同一个空间基准上，后续的汇水范围提取、断面上下游关系判断、面积量算都会出现系统性偏差——这种偏差不是随机误差，而是可以大到几百米甚至跨越整个流域的错位。

四类基础数据的角色分工

先把这四类数据在分析链条中的位置理清楚，才能判断统一基准时该以谁为准。

• DEM（数字高程模型）：提供连续的地面高程场，是汇水范围提取、流向计算、河网生成的唯一地形输入。它的分辨率决定了汇水边界的最小可分辨尺度，它的高程基准决定了水流方向的正确性。
• 矢量边界：包括流域边界、行政区边界、监测断面控制单元边界。它既是分析的裁剪范围，也是成果的统计单元，还常被用作汇水范围提取结果的校核参照。
• 水系数据：以线要素表达河流中心线，附带河流名称、等级、流向等属性。它用于确定断面在河网中的位置、判断上下游拓扑关系，也用于对DEM提取的河网进行比对修正。
• 断面点位：通常是点要素，带断面编号、名称、监测项目等属性。它是整个分析的锚点——所有汇水范围都是相对于某个断面向上游追溯得到的。

这四类数据中，DEM是栅格，其余三类是矢量。栅格与矢量的坐标转换规则不同，这是统一基准时最容易出问题的地方。

投影与坐标转换的关键判断

坐标系统包含两个层面：地理坐标系（如CGCS2000经纬度）和投影坐标系（如CGCS2000 3度带高斯-克吕格投影）。地理坐标系用角度单位，适合存储和交换；投影坐标系用长度单位，适合量算面积、距离和进行栅格运算。

汇水范围分析必须使用投影坐标系，原因很直接：流向计算依赖栅格的行列间距，如果间距以度为单位，不同纬度上同样的角度差对应的实际距离不同，坡度计算和流向判断都会失真。因此第一步是把所有数据转换到同一个投影坐标系下。以我国常用的情况为例，若研究区位于东经114°附近，可选用CGCS2000 3度带第38带（中央经线114°），投影坐标单位为米。

转换时要注意三个易错点：

1. 基准面转换与投影转换要分开处理。如果源数据是西安80或北京54坐标系，转到CGCS2000需要先做基准面转换（七参数或三参数），再做投影。直接把投影参数改掉而不做基准面转换，会引入几十到上百米的偏移。
2. 矢量与栅格的转换路径不同。矢量数据用投影变换工具直接转换坐标值；栅格数据除了投影变换，还需要重采样。重采样会改变像元值，必须选择合适的方法：高程用双线性或三次卷积，土地利用类型必须用最邻近法，否则会产生不存在的地类值。
3. 转换后必须做残差检查。选取至少3个分布均匀的控制点，比较转换前后在同一已知点上的坐标差。残差在厘米级可接受，若达到米级，说明参数或基准面设置有问题。

高程基准与栅格填洼

DEM的高程基准同样需要统一。我国常用1985国家高程基准，而部分历史数据可能使用1956黄海高程系，两者在多数地区相差几厘米到十几厘米，对平原河网地区的流向判断可能有影响。更关键的是DEM本身的误差：植被、建筑物造成的虚假洼地会让流向计算在中途“断流”，因此汇水范围提取前必须做填洼处理。填洼阈值要结合研究区地形设定，平原区阈值过大会抹平真实微地形，山区阈值过小则填不干净。

下面的图示展示了四类数据从各自原始状态汇聚到统一空间基准，再进入断面与汇水范围分析的处理关系。注意图中标注的转换动作与检查点，它们对应上面提到的关键判断。

一个可操作的基准统一检查清单

把上面的判断落到实际操作上，可以在正式提取汇水范围之前，用下面这份清单逐项确认。清单中的参数值属于示例设定，实际使用时需按研究区情况调整。

需要强调的是，统一基准不是一次性的格式转换，而是一个带验证的质量控制环节。转换完成后，最有效的验证方式是把断面点位、水系线和DEM晕渲图叠加在一起目视检查：断面点是否落在河道中心线上，水系线是否沿DEM中的谷底线延伸，矢量边界与提取出的汇水范围是否大体吻合。如果断面点明显偏离河道，或者水系线与地形谷底系统性错开，问题多半出在基准面转换或投影参数上，而不是后续的流向算法。

完成这一步之后，四类数据才真正处在同一个可比较、可运算的空间框架里。后续的填洼、流向计算、汇水范围提取，其结果的可靠性都建立在这个基准之上；基准不统一，再精细的算法也无法弥补输入层面的错位。