简介一套面向渭河流域研究与应用的空间数据包涵盖流域边界、高程、水系网络及省市县三级行政区划覆盖甘肃、陕西、宁夏3省区及20个地级市、112个县适合GIS分析、水文建模、生态评估与流域管理研究。数据基于高分辨率DEM、权威水系数据库和行政边界资料制作统一采用WGS84坐标系统经ArcGIS环境整理与一致性校核附带可编辑MXD工程、标准Shapefile矢量文件和高清TIF成图MXD已预设图层分组、图例、比例尺和注记样式TIF为300dpi高清输出。资源共57个文件以shp/shx/dbf等矢量要素、mxd制图文档、tif底图及辅助文档为主压缩包约202MB目录结构清晰便于直接调用。已有471人浏览学习适合需要快速获取渭河流域空间基础数据的学生、规划研究者与工程技术人员可直接用于流系提取、水源管理、生态连通分析和区域协调发展等工作。1. 三件套到手先别急出图mxd、shapefile、TIF 各管渭河流域数据的哪一段搞引水工程前期分析同事拷来一个数据包标题写着“渭河流域【矢量边界高程水系覆盖省市县】-可编辑mxd文件标准shape文件标准成图TIF”第一反应是双击 mxd 出图。翻车往往就在这第一步ArcMap 里的 mxd 只记图层路径和符号坐标数据本体全在 shapefile 和 TIF 里路径一断就只剩感叹号。这套打包方式本身是标准做法但“标准”二字只代表对方按规范存了不代表你打开就能套合。下面按我跑流域项目的流程来怎么验收这套数据、怎么把流域 DEM 裁出来以及 Clip 和掩膜提取的关键区别。这篇笔记既适合刚入行做水保、GIS 数据处理的新手也适合常年被坐标偏移问题折腾的熟手属于拿到现成数据包后必做的那套体检流程。2. 标准 shapefile 与 TIF 的验收清单从文件后缀到坐标系一步步确认数据没缺腿要理解“可编辑 mxd 标准 shapefile 标准成图 TIF”这套组合得先分清三者的职责边界。mxd 是 ArcMap 的地图文档记录图层列表、渲染符号、标注、比例尺等图面整饰信息它保存的是“这个数据包怎么展示”的说明书而不是坐标数据本身。你双击 mxdArcMap 会按文档里记录的数据源路径去磁盘上找对应 shp 或 tif找到就正常显示找不到就变成红感叹号。这也解释了为什么标题要强调“可编辑”——对方的意思是图层数据源都是原始文件不是导出成图片的静态图你还能继续改属性、调符号、做分析。2.1 mxd 存的是“地图的皮”矢量数据本体全在 shapefile 里我经手过的数据包问题恰恰常出在“可编辑”三个字被误读。有人以为 mxd 是万能容器拿过去改了两笔又另存结果原目录里的 shp 被删了mxd 变成空壳地图。正规的通行做法是把 mxd 和原始 shp、tif 放在同一个根目录里用相对路径指向数据源。ArcMap 默认存绝对路径你的工程目录和对方不一样所以打开第一件事就是看图层有没有红感叹号。如果有右键图层 → 属性 → 源 → 设置数据源逐个指回即可。对渭河流域这种跨甘肃、宁夏、陕西等多个省区的数据包shapefile 内部往往还按专题拆成多个图层流域边界、河流水系、湖库、高程点、行政分区等。mxd 就是把它们叠好、配色配好你打开就能当流域底图用。但如果 mxd 里引用的 shp 缺了 .prj坐标系信息丢失哪怕路径正常边界和水系叠到 DEM 上也是错位的。判断 mxd 是否真的“可编辑”有一个笨办法右键图层看“数据源”路径指向的是 shp 或 tif 而不是图片说明还能继续编辑指向 jpg、png 就只是背景图。2.2 标准 shapefile 不是单个文件后缀缺一不可.prj 决定套合成败“标准 shapefile”在 ArcGIS 里的硬标准是至少由 .shp、.shx、.dbf 三个基础文件组成。.shp 存几何坐标.shx 是几何索引.dbf 存属性字段。更完整的发行包还会带 .prj投影定义、.cpg属性编码、.sbn/.sbx空间索引、.xml元数据。文件少了ArcMap 会报“无法打开要素类”或“未找到对象”。传输时最常见的丢件是 .shx它很不起眼但一丢整个图层就打不开。网盘传几次、聊天软件发几次系统经常把同名的 .shp 和 .shx 当成重复文件过滤这是真事。拿到包先列目录比在 ArcMap 里挨个加载快得多。.prj 文件的重要性要单独提。渭河流域大约在东经 103°—110°、北纬 33°—37° 之间跨度大既有用经纬度的 WGS84、CGCS2000 地理坐标系也有用投影坐标的高斯克吕格 3 度带。一个数据包里边界和水系可能是 CGCS2000 3 度带 108E 投影DEM 却是 WGS84 经纬度叠在一起就是几十米到几百米的偏移这几乎是数据包验收时第一项要查的。检查方法是在图层属性 → 源里看“坐标系”或者右键图层 → 数据 → 属性 → XY 坐标系。更稳的是用一个小脚本读取每个要素类的 spatialReference后面第 6 章给现成代码。2.3 成图 TIF 与数据 TIF 的差别带地理编码的 GeoTIFF 才是成图标题里的“标准成图 TIF”一般指两种东西千万别混。第一种是成果制图把流域边界、水系、地形渲染结果按比例尺排版后导出成 TIF这种图往往没有坐标信息是给人看的图不能叠进 ArcMap 继续分析。第二种是带空间参考的 GeoTIFF像元里存着高程值或渲染值加载后能直接做坡度、淹没分析。判断方法很简单把 TIF 拖进 ArcMap右键属性看“空间参考”和“像元大小”是否都有值再拖动图层看能不能跟边界自动套合如果能说明带地理编码如果孤零零落在左下角基本是普通成图。对 DEM 类 TIF 还要查存储位深16 位有符号整数还是 32 位浮点。渭河河谷高程大约从 300 多米到 3000 多米的秦岭北麓16 位整型够用但做坡度计算和高程分级前建议转成 32 位浮点避免高程值被取整的精度问题。图层属性 → 源里可以看到像素类型。如果发现是 8 位 RGB 的渲染图那它不是真正的 DEM 数据只是地形阴影图标题里的“高程”就要打个问号得回头找原始栅格。按这个顺序验收数文件后缀 → 开 mxd 查红叹号 → 逐个图层看坐标系 → 把 TIF 和 shp 同时拖进去看是否重合四步走完数据包有没有缺腿就基本清楚了。3. 验证矢量边界、高程与水系把“覆盖省市县”从标题翻译成属性表检查3.1 用属性表验证省市县覆盖字段名与行政区代码是硬指标标题里的“覆盖省市县”落在 shapefile 属性表里通常是几列编码字段省名称、市名称、县名称外加对应的行政区划代码。拿到后先打开属性表看字段结构再做字段统计看有哪些省市县。渭河流域以甘肃省东部、宁夏南部、陕西省关中及陕北部分区域为主体属性表里应当能查到天水、平凉、定西、庆阳、固原、宝鸡、咸阳、西安、渭南、延安、铜川等地名。重点不是有没有这些名字而是行政区划代码字段是否完整——后续做县界统计、分市汇总靠的就是权威代码字段去关联别的业务表。这里有个高频踩坑点字段名乱码或属性显示为问号。标准 shapefile 的属性表如果是在中文 Windows 系统里建的编码通常是 GBK如果数据包后来被第三方工具重写过可能变成 UTF-8。打开的 dbf 如果不带 .cpg 文件ArcMap 默认按系统区域解释中文属性容易变成乱码。验证方法是在字段列表里看列名是否正常在属性表里看“陕西省”“甘肃省”这类值有没有显示成“????”。如果乱码补一个与 .shp 同名、放在同目录的 .cpg 文件内容写“UTF-8”或“OEM 936”重新加载图层往往能救回来。3.2 边界闭合与拓扑检查面图层最大的隐性坑矢量边界如果做的是面要素最基础的要求是“闭合”且“无自相交”。但野外采集或矢量化出来的流域边界常带尖角、重复节点和微小缝隙ArcMap 加载时看不出问题等你把它当作裁剪范围去裁 DEM出来的结果就缺一块或带锯齿。检查方式用 Data Management Tools → Features → Check Geometry它可以扫出空几何、自相交、重复顶点等病灶发现问题后用同目录下的 Repair Geometry 一键修。但要注意Repair Geometry 能修复拓扑错误解决不了“面不封闭”这种根本问题。如果检查发现要素类型是 Polyline 或面不封闭常见做法是先用 Feature To Polygon要素转面重建一次面再重新做拓扑检查。渭河流域这种边界图我一般把“县域面”和“流域面”分开查两者精度来源不同混在一起很难定位问题。3.3 高程和坐标系一致性河流、边界、DEM 必须在一个坐标系里数据包拆开看矢量边界、水系、DEM 高程、TIF 成图可能来自不同时期、不同采集批次坐标系不一致几乎是常态。ArcMap 10.x 在加载多个不同坐标系图层时会默认开启动态投影表面上看着重合但导出数据或做叠加分析时系统按各自坐标系参与计算结果就可能出现偏移或投影变形。最稳的做法是统一到 CGCS2000 或 WGS84。如果只做流域级分析常见做法是把所有数据都投影到 CGCS2000 3 度带 108 度中央经线的高斯投影这样面积和距离量测都贴近真实地面。矢量转换用 Data Management Tools → Projections and Transformations → Feature → Project栅格 DEM 换坐标系用 Project Raster 工具重采样方法选双线性或三次卷积不要用最邻近——最邻近在山地区域会损失高程细节生成阶梯状边缘。转换完成后在 ArcMap 里把边界、水系、DEM 三层全打开留意状态栏右下角的坐标读数看悬停光标时坐标是否在同一量级。如果只是“看着重合”但坐标读数差出几个数量级说明动态投影失效必须先解决坐标系问题再谈后续分析。3.4 数据完整性shp 缺文件、空要素、重复要素标题写了“标准 shapefile”数据包的完整性就得按标准抠。右键要素类属性在常规页看要素计数计数为 0图层就是空壳。渭河流域水系这种线图层正常应该有成百上千条河段如果统计出来只有几十条多半是对方只导出了主河道支流被遗漏。重复要素是另一类隐蔽问题两条完全或近乎相同的线段叠在一起出图看不出但计算水系总长度会直接翻倍。用 ArcToolbox 里的 Find Identical 工具把几何字段设为 shape查出来的重复记录删掉即可。数据包是别人打包的加工过程有没有产生重复顶点你无法预判体检查一次不亏。注意拿到数据包当天把原始文件整个复制一份到只读存档目录之后所有操作都在工作副本上做。原始 shp 一旦被 Repair Geometry 或字段编辑器改过返回去基本没有后悔药。4. 用面图层裁剪 DEM TIFArcMap 的 Clip 与掩膜提取到底该选谁4.1 两种工具的本质差异范围裁剪 vs 掩膜提取“arcmap 中依靠面图层裁剪 dem 栅格 tif 文件”和“依靠面图层掩码提取”是两类高频操作很多人混着用结果在边界精度和运行效率上吃亏。Clip裁剪在栅格里属于 Data Management Tools → Raster → Raster Processing → Clip。它默认按面图层的“外接矩形”去切除非勾选“Use Input Features for Clipping Geometry”使用输入要素裁剪几何才会按面边界切。不勾选时输出 DEM 是矩形四个角延伸到流域外拿去做面积统计会多算。Extract by Mask按掩膜提取在 Spatial Analyst Tools → Extraction → Extract by Mask它严格按面要素几何提取像元边界外像元变成 NoData输出栅格可以设置像元大小和重采样方式缺点是要启用 Spatial Analyst 扩展许可处理大栅格时内存占用高速度明显比 Clip 慢。给个判断题你要的是“渭河流域范围内的高程栅格边界外不要”优先用 Extract by Mask你要的是“拿流域面把大区域 DEM 快速切出一块后面再细算”用 Clip 加裁剪几何更快。实际项目里我通常先用 Clip 粗裁再用 Extract by Mask 精切两次配合既省时间又保精度。原因是预先裁剪一次后掩膜提取只需处理小范围栅格内存压力小出错率也低。4.2 实际操作步骤ArcMap 工具箱面板流程以 ArcMap 10.x 为例假设数据包里边界图层是Hydrological_Basin.shpDEM 是DEM_30m.tif操作步骤如下打开 ArcToolbox → Data Management Tools → Raster → Raster Processing → Clip。输入栅格选 DEM_30m.tifOutput Extent 下拉选择“要素类”并指定 Hydrological_Basin.shp勾选“Use Input Features for Clipping Geometry”。输出位置选工作目录文件名写 DEM_basin_clip.tif格式选 TIFF。执行 Clip 后再打开 Spatial Analyst Tools → Extraction → Extract by Mask。输入栅格选 DEM_basin_clip.tif输入掩膜数据选 Hydrological_Basin.shp输出栅格写 DEM_basin_mask.tif。参数上要注意Clip 工具的 NoData 值项如果留空输出栅格会沿用输入的 NoData 设置如果手动填 0会把数据里真实存在的高程 0 值也当成 NoData 丢掉。像元大小同理如果之前做过投影转换像元大小可能不是整数比如 29.9994 米不要手动改成 30除非你明确知道自己在做什么。在掩膜提取里输出像元大小默认取输入栅格但如果你希望输出与某个参考栅格对齐比如后续要与土地利用栅格叠加计算可以在环境设置里把栅格分析的“像元大小”指定为那个参考栅格。这两个细节直接决定两个栅格能不能严格像元对齐栅格计算前一定要先检查。4.3 Python 批处理脚本arcpy 一键输出流域 DEM 和水系叠加图在大范围流域分析中手工点工具箱要重复很多次尤其当你有多个子流域边界文件时。常见做法是把上面流程写成 arcpy 脚本。以下代码在 ArcMap 的 Python 窗口或独立 Python 环境里都能跑# -*- coding: utf-8 -*- import arcpy from arcpy.sa import ExtractByMask # 设置工作空间和许可 arcpy.env.workspace rD:\weishui arcpy.CheckOutExtension(Spatial) # 输入流域边界、原始 DEM、输出目录 basin rD:\weishui\Hydrological_Basin.shp dem rD:\weishui\DEM_30m.tif out_dir rD:\weishui\working # 第一步用面图层粗裁剪 DEM快速切出范围 arcpy.Clip_management( in_rasterdem, rectangle#, # 不手动指定矩形范围 out_rasterout_dir \\dem_clip.tif, in_template_datasetbasin, # 用边界要素类作为范围 nodata_value#, # 不覆盖原始 NoData 值 clipping_geometryClippingGeometry, # 按要素边界裁剪 maintain_clipping_extentFalse ) print(Clip 粗裁剪完成) # 第二步按掩膜精确提取边界外像元设为 NoData out_mask ExtractByMask(out_dir \\dem_clip.tif, basin) out_mask.save(out_dir \\dem_basin.tif) print(掩膜提取完成) arcpy.CheckInExtension(Spatial)脚本关键点Clip_management 的 clipping_geometry 参数在 ArcMap 10.x 里写 ClippingGeometry输出才会是面要素形状只写 NONE 就退化成外接矩形。ExtractByMask 来自 arcpy.sa要求脚本开头显式导入并 CheckOutExtension(Spatial)。第一次跑建议先选一个小范围边界测试确认输出栅格的像元大小、NoData 值符合预期再全量跑。脚本执行报错时优先看环境设置的输出坐标系——如果输入 DEM 是 WGS84 经纬度边界是 CGCS2000 投影动态投影会生效但缺少地理变换参数时输出照样偏移。解决办法是脚本开头显式设置 arcpy.env.outputCoordinateSystem别让坐标系问题变成黑匣子。5. 避坑记录渭河流域数据落到 ArcMap 最常见的 5 个翻车现场5.1 mxd 打开后一片红色感叹号相对路径没设置现象双击“可编辑 mxd 文件”加载后所有图层显示红色感叹号图层内容一片空白。 原因mxd 默认记录的是打包电脑上的绝对路径比如 E:\gisdata\shp\basin.shp你机器上没有这个盘符或目录ArcMap 找不到数据源。 解决先确认原始数据目录还在右键图层属性把数据源指回来只能解决当前会话。想根治打开菜单栏“文件 → 地图文档属性 → 数据源选项”勾选“存储相对路径名”再另存一次 mxd。以后整个目录任意移动只要 mxd 与数据目录的相对结构不变都不会断链。我见过不少项目团队把 mxd 和 shp 分别拷进“出图”和“原始数据”两个文件夹目录结构一变就全部断链这就是没有一开始约定存储结构的结果。5.2 边界叠加 DEM 整体偏移坐标系混用现象边界和水系相互能对上但一叠加 DEM所有矢量整体平移几十米到几百米。 原因矢量是 CGCS2000 高斯投影DEM 是 WGS84 地理坐标ArcMap 的动态投影虽然能在屏幕上勉强显示但叠加分析输出时坐标定义混乱。 解决先把所有图层统一到一个坐标系不要依赖动态投影做分析。矢量用 Project 工具栅格用 Project Raster 工具两者都必须使用同一套地理变换参数。对渭河流域我一般用 CGCS2000 3 度带 108 中央经线或 Albers 等积圆锥投影前者适合距离量测后者适合面积统计。重投影后开启数据框属性里的“变换”提示验证中心点坐标是否一致光看屏幕重合不算数。5.3 属性表中文全变问号缺少 .cpg 文件现象打开县界、水系图层的属性表省份、河流名称字段显示为“????”或乱码。 原因shapefile 的 dbf 属性字段编码可能是 GBK 或 UTF-8但目录里没有 .cpg 文件说明编码ArcMap 按系统默认区域解释解读方式不对就乱码。 解决在 shp 同目录下用记事本创建一个同名 .cpg 文件内容写“UTF-8”或“OEM 936”保存后重新加载图层。注意 .cpg 主文件名必须与 .shp 完全一致。这个小文件是中文 GIS 数据里最容易被忽略的后补药。如果两种编码都试过仍乱码说明 dbf 二进制内容被第三方工具改写需要用支持编码切换的 dbf 工具打开检查原始字段字节。5.4 面图层变线图层要素类型与文件名对不上现象边界.shp 加载进 ArcMap 后图标是线状用来做掩膜提取时提示“输入掩膜必须为栅格或面要素”。 原因原始边界可能是面但对方在后期编辑里把面转成了线保存或者整个“边界.shp”发布时本来就是线要素。要素类型写在 .shp 文件头里ArcMap 读的是类型定义不看文件名。 解决右键图层属性 → 源看数据格式里的要素类型。如果已变成线用 Feature To Polygon 重建面重建后面要素会丢失原属性需要用空间连接把属性按空间位置挂回来。我的习惯是把流域边界和水系分开存放边界只保留面水系只保留线从源头避免这种问题。5.5 TIF 成图背景一大片黑色NoData 与黑边处理不当现象把成图 TIF 叠加到边界图层上大片黑色背景把矢量全盖住。 原因GeoTIFF 里背景值常见为 0 或 255理论上应作为 NoData但 ArcMap 默认按黑色渲染或者对方导出时没有设置透明色。 解决右键图层 → 属性 → 符号系统把 NoData 颜色设为“无颜色”。如果是多波段 RGB 成图背景不是 NoData 而是 RGB(0,0,0)做法是用“按属性设置”把 0 值映射为透明或整体调整透明度。黑色背景问题在山区地形阴影图里尤其常见因为阴影本身偏暗不放大看很难判断是背景还是渲染内容建议切换到全图范围查看四角。6. 进阶做法写一段 arcpy 验收脚本给数据包做体检并导出报告这几步做完数据包基本能用了。但想让数据包以后不再“开盲盒”最好把验收动作固定成一个脚本每次拿到新数据都跑一遍。下面脚本在 ArcMap Python 窗口运行逐个检查要素类的要素类型、坐标系、要素数量和空几何输出一份文本报告# -*- coding: utf-8 -*- import arcpy arcpy.env.workspace rD:\weishui tol rD:\weishui\result.txt fc_list arcpy.ListFeatureClasses() with open(tol, w) as f: for fc in fc_list: desc arcpy.Describe(fc) sr desc.spatialReference count int(arcpy.GetCount_management(fc).getOutput(0)) f.write(图层: {} 类型: {} 数量: {}\n.format(fc, desc.shapeType, count)) f.write(坐标系: {}\n.format(sr.name)) with arcpy.da.SearchCursor(fc, [SHAPE]) as cur: empty 0 for row in cur: if row[0] is None or row[0].isNull: empty 1 f.write(空几何数量: {}\n\n.format(empty)) print(验收报告已生成打开查看, tol)如果在旧版 ArcMap 10.x 的 Python 2.7 环境下运行注意脚本里不要用 f-string上面统一用了 format路径含中文时要在字符串前加 u 前缀。这个脚本与其说是自动化验收不如说是一份“留痕”。项目周期拉长到半年以后回看当初用的是哪个版本的边界、坐标系是什么凭记忆很容易剪不断理还乱报告文件反而是最可靠的项目档案。我还习惯在验收后做一次子流域拆分测试随机挑渭河流域其中一段界跑一遍第 4 章的 Clip 和掩膜提取记下运行用时和输出像元大小判断全流域能不能在当天算完。说一个我自己的教训以前拿到这类数据包第一反应都是直接往工程里加等出图错了才回头查坐标、补文件。后来凡经手的流域数据无论来源多正规一律先跑这套体检脚本再叠 DEM 验证前后不到二十分钟省下来的是后面加班的数天。做 GIS 数据处理数据源越是不透明的下载包越该在进门时把户口查清楚。希望帮到你。本文还有配套的精品资源点击获取