简介面向地理信息、制图与区域研究人员的黄土高原空间分布及地形条件数据包以标准Shapefile矢量边界为核心可直接读取黄土高原面积、周长、经纬度等属性信息解决区域底图制作与地形分析的基础数据需求。包内共19个文件压缩包大小约186.45MB主要包含shp、dbf、prj等矢量组件、可编辑mxd制图工程、标准成图TIF以及全国高程TIF数据并附有中国行政区等辅助zip覆盖从矢量编辑到栅格成图的常用需求。mxd工程便于按个人偏好调整版式标准TIF成图可即拿即用全国高程数据则支持绘制全国及黄土高原区域的高程专题图适合教学、科研和日常制图场景。目前已有83人学习下载是快速获取权威空间分布与地形底图的实用选择。1. 拿到“黄土高原空间分布及地形条件”资源包后先别急着双击那个 MXD第一次接触到这种资源包的人目光往往会被 .mxd 文件吸住感觉双击打开就能看到一幅现成的黄土高原地形图。实际干过一轮就会发现这套包里真正决定成图质量的是那几个“标准 shape 文件”和那张“全国高程文件 TIF”而 MXD 更像一个把所有素材按制图规范绑到一起的框架。资源包解决的痛点是你不用再花一周时间去搜集黄土高原边界、拼接全国 DEM也不用从零去调色带、排图例。可它也不是傻瓜包常见翻车点是数据源路径失效以及把“裁剪”和“掩码提取”当成一回事。这篇文章我会按常规落地流程拆开讲文件结构、ArcMap 中依靠面图层裁剪 DEM 栅格 TIF 的两种做法、标准成图 TIF 的导出设置、以及进入生产前最值得看的几个坑。适合需要做地貌制图、水土保持或流域分析的人参考也适合 GIS 新手把这套包当成练习模板。2. 拆开数据包可编辑 MXD、标准 shape 文件和全国高程 TIF 各是干什么的整个资源包可以按“数据 → 制图 → 输出”三层理解。标准 shape 文件提供黄土高原的空间分布边界全国高程 TIF 提供地形条件的基础栅格可编辑 MXD 负责把这两类数据组合成一张带有图例、比例尺、指北针的专题图最后的“标准成图 TIF”则是面向甲方或报告交付的成品影像。四者不是并列关系MXD 依赖前两者成图 TIF 又依赖 MXD 的布局与符号设置。很多人拿到后直接在 ArcMap 里双击 MXD内容列表一片红色感叹号于是断定资源包有问题。其实 MXD 不保存栅格和矢量数据它只保存数据源路径和显示规则。路径一旦变化链接必然断。理解这一点后下面拆开每个文件类型。2.1 可编辑 MXD制图模板没数据也能复用可编辑 MXD 的价值在于“制图模板化”。黄土高原沟壑纵横地形起伏大如果 DEM 用默认的全蓝色带或 green-to-brown 色带很难体现黄土高原的干沟万壑。有经验的制图员会在符号系统里把海拔分成 58 级用浅黄、土黄、棕褐这一套颜色来表达黄土阶地和高原面。MXD 把这些配色记录了下来你拿到后不需要重新做一遍符号工程。同时MXD 的布局视图里已经排好了图例、指北针、比例尺、图名和图框。你只要把数据源修好图面排布几乎不用动。但也正因为“可编辑”它很容易被改乱。我的建议是先复制一份原始 MXD命名为“原版模板.mxd”存档再在副本上做实验。这样哪怕把自己的配色调翻了车也有后悔药可吃。如果你打开 MXD 后发现某个图层前面有红色感叹号那就是数据源路径失效了。右键图层 → “数据” → “修复数据源”把路径指到对应文件就好。对于一次要修复几十个图层的情况后面第 5 章会给出批量方案。2.2 标准 shape 文件边界规范直接影响裁剪范围“标准 shape 文件”重点在“标准”二字。Shapefile 不是一个单独文件而是 .shp、.dbf、.shx、.prj 四个文件组成有时还带 .cpg。标准 shapefile 意味着坐标系定义完整、字段命名一致、几何没有明显拓扑错误。这个资源包里的 shape 文件通常至少有一个面图层代表黄土高原的分布边界。它是后续裁剪 DEM 的几何依据。如果面图层的边界有缝隙、自相交或者坐标范围异常裁剪出来的栅格边缘就会出现条带状 NoData影响坡度、坡向等地形分析结果。我拿到 shapefile 后有三步固定动作第一步在 ArcCatalog 中查看属性确认 .prj 里写的是 WGS84 还是 CGCS2000第二步加载到 ArcMap打开“要素类属性”里的“范围”选项卡看看经纬度或投影坐标是否合理第三步用“地理处理 → 修复几何”工具跑一次把所有潜在的多边形错误修掉。这一步看起来多余但能省下后面不少排查时间。2.3 全国高程文件 TIF 与坐标系的取舍全国高程文件 TIF 一般是覆盖全国的 DEM像元大小可能到 30 米。直接拖进 ArcMap缩放手感会非常滞涩内存占用也高。这不是 TIF 坏了而是没有建金字塔或者没有裁剪。为了后续产线顺畅我的做法是先把全国 TIF 建金字塔然后在 MXD 数据框里只保留黄土高原书签范围。有人会问既然只是做黄土高原为什么不直接提供一个局部 TIF因为“空间分布及地形条件”这个题目往往还要分析流域连通性或者对比整个黄河中游的地形特征。局部 TIF 在边缘处会截断河网全国 TIF 保留了更完整的分析弹性。但保留全国 TIF 不意味着全图参与计算。在做任何坡度、坡向分析前都需要用标准 shape 文件把 DEM 按边界裁出来。裁剪前还要确定坐标系如果 DEM 是地理坐标系WGS84 经纬度计算坡度得到的是“度”如果要算面积、周长或坡度百分比则应该先把 DEM 投影到 Albers 等积圆锥或高斯克吕格投影。常规做法是在 ArcToolbox 里用“投影栅格”重采样方法选双线性像元大小保持 30 米输出坐标系选“Asia North Albers Equal Area Conic”。注意不要选默认的“邻近”重采样否则坡面会出阶梯状锯齿看起来像像素画。3. ArcMap 中依靠面图层裁剪 DEM 栅格 TIF 文件裁剪和掩码提取的区别这一章进入真正动手的环节。标题里写的“黄土高原空间分布及地形条件”能不能落成可用栅格靠的就是这一步。ArcMap 里有两处工具都叫“裁剪”一个是“数据管理工具 → 栅格 → 栅格处理 → 裁剪”另一个是“空间分析工具 → 提取分析 → 按掩膜提取”。检索时大家常把后者说成“掩码提取”其实它对应的英文是 Extract by Mask。两者算法不同结果面差很大。3.1 裁剪只是“框选”掩膜提取才“扣边”先说底层逻辑。裁剪工具Clip本质上是拿面要素的外接矩形去切栅格输出结果还是一个矩形只是矩形的范围缩到了面要素的外接矩形四角全部被 NoData 填充。比如黄土高原边界是一个非常不规则的肺叶状轮廓用裁剪工具得到的栅格仍是一个规整矩形文件体积一点没少因为那些 NoData 区域也占着像元位置。按掩膜提取Extract by Mask则不同它会以面要素的几何形状为掩膜逐像元判断是否落于面内部落在面外的像元直接设为 NoData。输出的栅格不再是矩形而是贴合边界的不规则形状像元矩阵紧凑后续做面积统计时不会把外部的噪音算进去。所以选型逻辑很直接如果只是临时看一眼或者最终成果要转成 ASCII 和反演算法对接用裁剪工具没关系但如果是做标准成图 TIF、面积统计或地形因子计算必须用按掩膜提取。用错工具后最诡异的现象是明明能看到黄土高原轮廓统计面积却偏大原因就是 NoData 区域被当成有效区域参与了某些映射计算。3.2 手工操作“按掩膜提取”的最短步骤在 ArcMap 10.x 系列版本里最短操作步骤如下在 ArcMap 内容列表加载“全国高程文件 TIF”。加载“标准 shape 文件”中的黄土高原边界面图层。点击“地理处理”下拉菜单搜索“按掩膜提取”。在工具对话框中输入栅格选全国 TIF输入栅格数据或要素掩膜数据选面图层输出栅格路径写到你自己的工作空间。点击“环境”按钮展开“栅格分析”输出坐标系选“与输入相同”处理范围选“与掩膜图层相同”像元大小选“与输入 DEM 相同”。第 5 步是很多结果与预期不符的根源。特别是“处理范围”如果保持默认的“全图”工具会把掩膜提取后的大范围 NoData 像元全部保留输出栅格的行列数和全国 DEM 一样大文件大小也就降不下来。这个操作看起来跑完了其实等于没做。跑完后在内容列表把输出层的符号系统改成“拉伸”拉伸类型选“标准偏差”色带选从浅黄到深棕的渐变。你会看到黄土高原边界内部是连续海拔渐变边界外是透明的 NoData这个效果就是掩码提取和裁剪的最大区别。3.3 用 arcpy 脚本批量裁剪多区域一次跑完手工 UI 适合单次工作。如果要把黄土高原按流域拆成十几个子区域每次走一遍“地理处理 → 按掩膜提取”会非常繁琐。我一般会把流程写成 Python 脚本放进 ArcMap 的 Python 窗口跑。下面是一段可复用代码# -*- coding: utf-8 -*- import arcpy from arcpy import env from arcpy.sa import * # 必须先开启空间分析扩展否则工具会报错 arcpy.CheckOutExtension(Spatial) env.workspace rD:\loess_data\project.gdb env.overwriteOutput True # 输入数据全国 DEM TIF 黄土高原边界 shapefile dem_tif rD:\source\dem_wgs84.tif mask_shp rD:\source\loess_boundary.shp # 让分析范围严格等于掩膜要素范围避免多出一圈 NoData mask_desc arcpy.Describe(mask_shp) env.extent mask_desc.extent out_raster rD:\loess_data\loess_dem_masked.tif out ExtractByMask(dem_tif, mask_shp) out.save(out_raster) # 顺手生成坡度验证掩膜后的地形连续性 slope Slope(out, DEGREE) slope.save(rD:\loess_data\loess_slope_deg.tif)代码的逻辑分为四块。第一块“CheckOutExtension”是硬门槛没有它空间分析工具不会被授权脚本会在 ExtractByMask 这一行直接终止。第二块“env.extent”是从掩膜 shapefile 读取范围并写进环境变量这一步效仿了手工操作里“处理范围”的设置。第三块是执行提取并保存注意“ExtractByMask”的第二个参数既可以是点、线、面要素也可以是栅格传 shapefile 时它内部会先转成掩膜。第四块用坡度数据做质检如果掩膜边缘有缺陷坡度结果在边界处会出现异常栅格。这段脚本里最容易被忽略的是arcpy.Describe(mask_shp).extent。直接给env.extent赋值一个 shapefile 路径是不被接受的必须先通过 Describe 拿到 Extent 对象。另外建议把env.snapRaster也指向 DEM这样输出像元会与 DEM 严格对齐避免后续叠加时出现半个像元的错位。批量化时在外层套一个for循环逐个读取要素输出文件加流域名后缀即可。4. 标准成图 TIF在可编辑 MXD 里安排布局再导出参数不能拍脑袋填经过掩膜提取数据层面已经就绪。下一步是制作“标准成图 TIF”。这个 TIF 不是简单把 DEM 导出图片而是把黄土高原的高程分布、坡度分级、行政边界和图例组成一张可打印的专题图。ArcMap 里一般通过可编辑 MXD 的“布局视图”来完成。4.1 把掩膜结果和 shapefile 装进 MXD 布局打开原始 MXD如果内容列表已经预置了图层组那最好。没有的话按照“边界 地形 辅助要素”三层结构添加。第一层放标准 shapefile 的边界面第二层放按掩膜提取后的 DEM 栅格第三层放等高线或水系等辅助数据。图层的顺序不能乱边界层通常放在 DEM 上方并设置为空心填充避免盖住地形效果。在内容列表里右键 DEM 层的符号系统选择“分类”而不是“拉伸”分类数设在 58 类。黄土高原这种高程差相对集中的区域8 类足够把低海拔河谷和高原面区分开。图例上要明确“高程米”。如果 MXD 里已有图例需要右键图例查看它引用的图层名一旦图层名变了图例会变成空白。然后切换到“布局视图”。图纸上会有一个数据框数据框里显示的就是前面内容列表的数据。要在保证比例尺整数化的前提下把数据框范围调整到黄土高原边界附近。常见做法是先用“书签管理器”添加一个目标范围书签再在数据框属性 → “数据框范围”里选择“其他”并引用该书签。这样输出范围不会因为鼠标拖动而偏离。4.2 导出 TIF 的参数设置与常见误解布局完成后选“文件 → 导出地图”保存类型选 TIFF。导出设置里几个关键点依次说明分辨率打印输出用 300 dpi仅在屏幕或 PDF 里看150 dpi 足够。A3 幅面 300 dpi 大约产生 4500×6200 像素的 TIF再高除了拖慢加载没有更多收益。压缩选 LZW。LZW 是无损压缩对 DEM 这类渐变栅格效果很好JPEG 压缩会在边界处产生块状瑕疵不适合标准成图。颜色模式工程协作默认 RGB出印刷品再选 CMYK。背景色透明如果在 ArcMap 成图里需要把数据框外的白色区域去掉勾选“背景色透明”并把“背景颜色”设置为一个不常用的纯色如 255,0,255。否则输出 TIF 在其他软件里会有白色边框盖住底图。很多第一次导出的人以为“成图 TIF”就是把这个 TIF 再导出一遍误会就大了。MXD 里引用的 DEM 是数值栅格而导出的成图 TIF 是 RGB 影像是配色后的结果不再保留高程值。也就是说成图 TIF 只用来展示汇报面积统计和坡度计算还得用第三章的输出。4.3 数据刷新的复现工作流因为 MXD 可编辑这套工作流可以反复执行。典型场景当拿到新版 DEM 或更新后的边界 shapefile不必重建 MXD只需在内容列表右键旧栅格 → “移除”再“添加数据”选新掩膜结果之后重新导出 TIF 即可。但要注意掩膜提取输出和布局导出不要写成循环嵌套即不要在脚本里把“掩膜提取结果 .tif”当作 MXD 数据源时又对这个 TIF 执行“按掩膜提取”那样会把 NoData 区域当作黑边反复叠加最后成图边界一圈黑框。为了可追溯我会把 MXD 与数据源放在同一个工作目录并保持相对路径。如果你的数据包是拷贝到另一台电脑打开“文件 → 地图文档属性”勾选“存储相对路径名”。这一步能避免路径变化带来的红感叹号属于一劳永逸的小技巧。5. 五条避坑记录从 MXD 数据源到掩膜提取的边界问题下面这些坑都是实际项目里反复出现过的按“现象 → 原因 → 解决”记下来方便对号入座。5.1 MXD 换电脑后全是红色感叹号现象把 MXD 连带数据包拷贝到另一台电脑打开后每个图层前面都是红色感叹号双击图层没有反应。 原因MXD 默认存储绝对路径。原电脑上的路径是D:\gisdata\...新电脑没有这个目录。 解决打开 MXD 后右键任意图层 → “数据” → “修复数据源”手动定位到数据包新路径。如果图层很多干脆用 arcpy 脚本把所有图层的数据源一次重定位。做法是先遍历arcpy.mapping.ListLayers(doc)用findAndReplaceWorkspacePath()把根路径统一替换。5.2 “裁剪”和“掩码提取”结果范围一样统计结果却差一截现象有人在“数据管理 → 裁剪”工具下选了“使用输入要素裁剪几何”结果看起来和掩膜提取差不多但用栅格统计出的总面积明显偏大。 原因裁剪工具即使勾选几何也只是把被裁剪区域之外的像元设为 NoData但栅格矩阵仍是矩形掩膜提取输出的是紧凑矩阵。后续统计面积时很多统计函数会跳过 NoData但由于范围不同、像元对齐不同结果自然会偏差。 解决需要统计面积或写入 ASCII 的成果一律走“按掩膜提取”用裁剪工具出来的数据只做快速浏览。识别方法很简单打开栅格属性看行数和列数。如果行列数与全国 DEM 一样说明范围和范围没缩到位掩膜提取失败或环境范围没设置。5.3 掩膜提取后 TIF 还是几个 GB文件体积没变现象明明做了按掩膜提取输出 TIF 的文件大小和原全国 DEM 差不多磁盘占用一点没少。 原因环境设置里的“处理范围”没有设置。ArcMap 默认处理范围是输入栅格全图掩膜提取只在掩膜范围外写 NoData但矩阵尺寸仍然是全国尺寸。 解决在按掩膜提取对话框点“环境”把“栅格分析 → 处理范围”改成“与图层相同”然后指定黄土高原边界图层。脚本里用arcpy.env.extent arcpy.Describe(mask_shp).extent。如果已经生成再用“数据管理 → 栅格 → 裁剪”裁掉 NoData 外框文件也会正常瘦身。5.4 全国高程 TIF 加载后 MXD 卡死转圈十几分钟现象把全国 TIF 拖进 MXD 后缩放和平移都非常卡内存占用居高不下甚至提示脚本停止响应。 原因TIF 没有构建金字塔ArcMap 每次重绘都要读取全分辨率像元。 解决在 ArcCatalog 中右键 TIF → “属性” → “金字塔” → 构建金字塔。更自动化一点用 Pythonimport arcpy arcpy.BuildPyramids_management( rD:\source\dem_wgs84.tif, PYRAMIDS, 2, BILINEAR )参数说明PYRAMIDS表示构建内部金字塔文件2是重采样级别取 1 或 2 即可BILINEAR是重采样方法DEM 用双线性比最近邻法平滑。构建完后再回到 MXD你会发现缩放流畅多了。另外把数据框坐标系固定成投影坐标系也能减少动态投影的计算压力。5.5 成图 TIF 边界出现黑框或白边现象导出的成图 TIF 在 PhotoShop 或报告里打开黄土高原边界外是一圈黑色或白色块叠在其他图上会遮挡底图。 原因导出时没有处理背景色透明。如果数据框外空白区域默认是白色黑框则通常在数据框内 NoData 区与图例透明区域冲突时出现。 解决导出地图时勾选“背景色透明”并将背景颜色设置为罕见色同时在图层属性里把 DEM 的 NoData 颜色设置为“无颜色”。之前有次我忘了勾选成图 TIF 叠加在卫星底图上边界外一大块白被客户一眼看穿这就是经验所在。6. 验证成果用 Python 快速核对 TIF 与 shapefile 的范围和投影到了最后一个环节。不管是用裁剪还是掩膜提取最终都要验证输出产品没有范围错位和投影丢失。肉眼判断在地形复杂区域很容易漏掉问题我习惯用一段短脚本在交付前检查栅格与矢量的一致性。from osgeo import gdal, ogr dem gdal.Open(rD:\loess_data\loess_dem_masked.tif) gt dem.GetGeoTransform() x_res abs(gt[1]) y_res abs(gt[5]) cols dem.RasterXSize rows dem.RasterYSize left gt[0] top gt[3] right left cols * gt[1] bottom top rows * gt[5] print(像元大小:, x_res, y_res) print(栅格范围:, left, bottom, right, top) shp ogr.Open(rD:\source\loess_boundary.shp) lyr shp.GetLayer(0) ext lyr.GetExtent() print(矢量范围:, ext[0], ext[2], ext[1], ext[3])这段代码做的验证很有意思读取栅格的仿射变换算出四个角坐标跟 shapefile 的范围对比。如果掩膜提取时的环境范围没有设置到边界自身你会看到栅格范围远大于矢量范围证明输出里面还有大块 NoData。像元大小也可以顺带检查如果原本是 30 米出来却变成 0.000333 度说明 DEM 还在地理坐标系要再确认是否满足你的分析需求。此外还要看投影字符串是否一致dem.GetProjection()和lyr.GetSpatialRef().ExportToWkt()可以打印出来对比。遇到坐标系不一致最简单的做法是在 ArcMap 里对输出的栅格做“投影栅格”转成与 shapefile 一致的投影。不要相信 MXD 数据框的动态投影那只改变显示不改变数据。我的个人记性是每跑完一步就把成果另存为带版本号的 TIF比如loess_dem_mask_v2.tif这样既能回退也能在交付前保留一条完整的可追溯链路。这种习惯在反复调试制图模板时尤其有用毕竟地理处理的“黑匣子”有时候不是报错而是悄悄让你在边界上多一条缝。希望这套流程能帮你把这个资源包真正用起来做到黄土高原空间分布与地形条件的数据可以干净、快速地落到图上。本文还有配套的精品资源点击获取 
长河君 · 2026/10/8 16:58:13 · 阅读 1.2k ArcMap中裁剪DEM做黄土高原地形图：掩膜提取与成图TIF避坑指南
简介面向地理信息、制图与区域研究人员的黄土高原空间分布及地形条件数据包以标准Shapefile矢量边界为核心可直接读取黄土高原面积、周长、经纬度等属性信息解决区域底图制作与地形分析的基础数据需求。包内共19个文件压缩包大小约186.45MB主要包含shp、dbf、prj等矢量组件、可编辑mxd制图工程、标准成图TIF以及全国高程TIF数据并附有中国行政区等辅助zip覆盖从矢量编辑到栅格成图的常用需求。mxd工程便于按个人偏好调整版式标准TIF成图可即拿即用全国高程数据则支持绘制全国及黄土高原区域的高程专题图适合教学、科研和日常制图场景。目前已有83人学习下载是快速获取权威空间分布与地形底图的实用选择。1. 拿到“黄土高原空间分布及地形条件”资源包后先别急着双击那个 MXD第一次接触到这种资源包的人目光往往会被 .mxd 文件吸住感觉双击打开就能看到一幅现成的黄土高原地形图。实际干过一轮就会发现这套包里真正决定成图质量的是那几个“标准 shape 文件”和那张“全国高程文件 TIF”而 MXD 更像一个把所有素材按制图规范绑到一起的框架。资源包解决的痛点是你不用再花一周时间去搜集黄土高原边界、拼接全国 DEM也不用从零去调色带、排图例。可它也不是傻瓜包常见翻车点是数据源路径失效以及把“裁剪”和“掩码提取”当成一回事。这篇文章我会按常规落地流程拆开讲文件结构、ArcMap 中依靠面图层裁剪 DEM 栅格 TIF 的两种做法、标准成图 TIF 的导出设置、以及进入生产前最值得看的几个坑。适合需要做地貌制图、水土保持或流域分析的人参考也适合 GIS 新手把这套包当成练习模板。2. 拆开数据包可编辑 MXD、标准 shape 文件和全国高程 TIF 各是干什么的整个资源包可以按“数据 → 制图 → 输出”三层理解。标准 shape 文件提供黄土高原的空间分布边界全国高程 TIF 提供地形条件的基础栅格可编辑 MXD 负责把这两类数据组合成一张带有图例、比例尺、指北针的专题图最后的“标准成图 TIF”则是面向甲方或报告交付的成品影像。四者不是并列关系MXD 依赖前两者成图 TIF 又依赖 MXD 的布局与符号设置。很多人拿到后直接在 ArcMap 里双击 MXD内容列表一片红色感叹号于是断定资源包有问题。其实 MXD 不保存栅格和矢量数据它只保存数据源路径和显示规则。路径一旦变化链接必然断。理解这一点后下面拆开每个文件类型。2.1 可编辑 MXD制图模板没数据也能复用可编辑 MXD 的价值在于“制图模板化”。黄土高原沟壑纵横地形起伏大如果 DEM 用默认的全蓝色带或 green-to-brown 色带很难体现黄土高原的干沟万壑。有经验的制图员会在符号系统里把海拔分成 58 级用浅黄、土黄、棕褐这一套颜色来表达黄土阶地和高原面。MXD 把这些配色记录了下来你拿到后不需要重新做一遍符号工程。同时MXD 的布局视图里已经排好了图例、指北针、比例尺、图名和图框。你只要把数据源修好图面排布几乎不用动。但也正因为“可编辑”它很容易被改乱。我的建议是先复制一份原始 MXD命名为“原版模板.mxd”存档再在副本上做实验。这样哪怕把自己的配色调翻了车也有后悔药可吃。如果你打开 MXD 后发现某个图层前面有红色感叹号那就是数据源路径失效了。右键图层 → “数据” → “修复数据源”把路径指到对应文件就好。对于一次要修复几十个图层的情况后面第 5 章会给出批量方案。2.2 标准 shape 文件边界规范直接影响裁剪范围“标准 shape 文件”重点在“标准”二字。Shapefile 不是一个单独文件而是 .shp、.dbf、.shx、.prj 四个文件组成有时还带 .cpg。标准 shapefile 意味着坐标系定义完整、字段命名一致、几何没有明显拓扑错误。这个资源包里的 shape 文件通常至少有一个面图层代表黄土高原的分布边界。它是后续裁剪 DEM 的几何依据。如果面图层的边界有缝隙、自相交或者坐标范围异常裁剪出来的栅格边缘就会出现条带状 NoData影响坡度、坡向等地形分析结果。我拿到 shapefile 后有三步固定动作第一步在 ArcCatalog 中查看属性确认 .prj 里写的是 WGS84 还是 CGCS2000第二步加载到 ArcMap打开“要素类属性”里的“范围”选项卡看看经纬度或投影坐标是否合理第三步用“地理处理 → 修复几何”工具跑一次把所有潜在的多边形错误修掉。这一步看起来多余但能省下后面不少排查时间。2.3 全国高程文件 TIF 与坐标系的取舍全国高程文件 TIF 一般是覆盖全国的 DEM像元大小可能到 30 米。直接拖进 ArcMap缩放手感会非常滞涩内存占用也高。这不是 TIF 坏了而是没有建金字塔或者没有裁剪。为了后续产线顺畅我的做法是先把全国 TIF 建金字塔然后在 MXD 数据框里只保留黄土高原书签范围。有人会问既然只是做黄土高原为什么不直接提供一个局部 TIF因为“空间分布及地形条件”这个题目往往还要分析流域连通性或者对比整个黄河中游的地形特征。局部 TIF 在边缘处会截断河网全国 TIF 保留了更完整的分析弹性。但保留全国 TIF 不意味着全图参与计算。在做任何坡度、坡向分析前都需要用标准 shape 文件把 DEM 按边界裁出来。裁剪前还要确定坐标系如果 DEM 是地理坐标系WGS84 经纬度计算坡度得到的是“度”如果要算面积、周长或坡度百分比则应该先把 DEM 投影到 Albers 等积圆锥或高斯克吕格投影。常规做法是在 ArcToolbox 里用“投影栅格”重采样方法选双线性像元大小保持 30 米输出坐标系选“Asia North Albers Equal Area Conic”。注意不要选默认的“邻近”重采样否则坡面会出阶梯状锯齿看起来像像素画。3. ArcMap 中依靠面图层裁剪 DEM 栅格 TIF 文件裁剪和掩码提取的区别这一章进入真正动手的环节。标题里写的“黄土高原空间分布及地形条件”能不能落成可用栅格靠的就是这一步。ArcMap 里有两处工具都叫“裁剪”一个是“数据管理工具 → 栅格 → 栅格处理 → 裁剪”另一个是“空间分析工具 → 提取分析 → 按掩膜提取”。检索时大家常把后者说成“掩码提取”其实它对应的英文是 Extract by Mask。两者算法不同结果面差很大。3.1 裁剪只是“框选”掩膜提取才“扣边”先说底层逻辑。裁剪工具Clip本质上是拿面要素的外接矩形去切栅格输出结果还是一个矩形只是矩形的范围缩到了面要素的外接矩形四角全部被 NoData 填充。比如黄土高原边界是一个非常不规则的肺叶状轮廓用裁剪工具得到的栅格仍是一个规整矩形文件体积一点没少因为那些 NoData 区域也占着像元位置。按掩膜提取Extract by Mask则不同它会以面要素的几何形状为掩膜逐像元判断是否落于面内部落在面外的像元直接设为 NoData。输出的栅格不再是矩形而是贴合边界的不规则形状像元矩阵紧凑后续做面积统计时不会把外部的噪音算进去。所以选型逻辑很直接如果只是临时看一眼或者最终成果要转成 ASCII 和反演算法对接用裁剪工具没关系但如果是做标准成图 TIF、面积统计或地形因子计算必须用按掩膜提取。用错工具后最诡异的现象是明明能看到黄土高原轮廓统计面积却偏大原因就是 NoData 区域被当成有效区域参与了某些映射计算。3.2 手工操作“按掩膜提取”的最短步骤在 ArcMap 10.x 系列版本里最短操作步骤如下在 ArcMap 内容列表加载“全国高程文件 TIF”。加载“标准 shape 文件”中的黄土高原边界面图层。点击“地理处理”下拉菜单搜索“按掩膜提取”。在工具对话框中输入栅格选全国 TIF输入栅格数据或要素掩膜数据选面图层输出栅格路径写到你自己的工作空间。点击“环境”按钮展开“栅格分析”输出坐标系选“与输入相同”处理范围选“与掩膜图层相同”像元大小选“与输入 DEM 相同”。第 5 步是很多结果与预期不符的根源。特别是“处理范围”如果保持默认的“全图”工具会把掩膜提取后的大范围 NoData 像元全部保留输出栅格的行列数和全国 DEM 一样大文件大小也就降不下来。这个操作看起来跑完了其实等于没做。跑完后在内容列表把输出层的符号系统改成“拉伸”拉伸类型选“标准偏差”色带选从浅黄到深棕的渐变。你会看到黄土高原边界内部是连续海拔渐变边界外是透明的 NoData这个效果就是掩码提取和裁剪的最大区别。3.3 用 arcpy 脚本批量裁剪多区域一次跑完手工 UI 适合单次工作。如果要把黄土高原按流域拆成十几个子区域每次走一遍“地理处理 → 按掩膜提取”会非常繁琐。我一般会把流程写成 Python 脚本放进 ArcMap 的 Python 窗口跑。下面是一段可复用代码# -*- coding: utf-8 -*- import arcpy from arcpy import env from arcpy.sa import * # 必须先开启空间分析扩展否则工具会报错 arcpy.CheckOutExtension(Spatial) env.workspace rD:\loess_data\project.gdb env.overwriteOutput True # 输入数据全国 DEM TIF 黄土高原边界 shapefile dem_tif rD:\source\dem_wgs84.tif mask_shp rD:\source\loess_boundary.shp # 让分析范围严格等于掩膜要素范围避免多出一圈 NoData mask_desc arcpy.Describe(mask_shp) env.extent mask_desc.extent out_raster rD:\loess_data\loess_dem_masked.tif out ExtractByMask(dem_tif, mask_shp) out.save(out_raster) # 顺手生成坡度验证掩膜后的地形连续性 slope Slope(out, DEGREE) slope.save(rD:\loess_data\loess_slope_deg.tif)代码的逻辑分为四块。第一块“CheckOutExtension”是硬门槛没有它空间分析工具不会被授权脚本会在 ExtractByMask 这一行直接终止。第二块“env.extent”是从掩膜 shapefile 读取范围并写进环境变量这一步效仿了手工操作里“处理范围”的设置。第三块是执行提取并保存注意“ExtractByMask”的第二个参数既可以是点、线、面要素也可以是栅格传 shapefile 时它内部会先转成掩膜。第四块用坡度数据做质检如果掩膜边缘有缺陷坡度结果在边界处会出现异常栅格。这段脚本里最容易被忽略的是arcpy.Describe(mask_shp).extent。直接给env.extent赋值一个 shapefile 路径是不被接受的必须先通过 Describe 拿到 Extent 对象。另外建议把env.snapRaster也指向 DEM这样输出像元会与 DEM 严格对齐避免后续叠加时出现半个像元的错位。批量化时在外层套一个for循环逐个读取要素输出文件加流域名后缀即可。4. 标准成图 TIF在可编辑 MXD 里安排布局再导出参数不能拍脑袋填经过掩膜提取数据层面已经就绪。下一步是制作“标准成图 TIF”。这个 TIF 不是简单把 DEM 导出图片而是把黄土高原的高程分布、坡度分级、行政边界和图例组成一张可打印的专题图。ArcMap 里一般通过可编辑 MXD 的“布局视图”来完成。4.1 把掩膜结果和 shapefile 装进 MXD 布局打开原始 MXD如果内容列表已经预置了图层组那最好。没有的话按照“边界 地形 辅助要素”三层结构添加。第一层放标准 shapefile 的边界面第二层放按掩膜提取后的 DEM 栅格第三层放等高线或水系等辅助数据。图层的顺序不能乱边界层通常放在 DEM 上方并设置为空心填充避免盖住地形效果。在内容列表里右键 DEM 层的符号系统选择“分类”而不是“拉伸”分类数设在 58 类。黄土高原这种高程差相对集中的区域8 类足够把低海拔河谷和高原面区分开。图例上要明确“高程米”。如果 MXD 里已有图例需要右键图例查看它引用的图层名一旦图层名变了图例会变成空白。然后切换到“布局视图”。图纸上会有一个数据框数据框里显示的就是前面内容列表的数据。要在保证比例尺整数化的前提下把数据框范围调整到黄土高原边界附近。常见做法是先用“书签管理器”添加一个目标范围书签再在数据框属性 → “数据框范围”里选择“其他”并引用该书签。这样输出范围不会因为鼠标拖动而偏离。4.2 导出 TIF 的参数设置与常见误解布局完成后选“文件 → 导出地图”保存类型选 TIFF。导出设置里几个关键点依次说明分辨率打印输出用 300 dpi仅在屏幕或 PDF 里看150 dpi 足够。A3 幅面 300 dpi 大约产生 4500×6200 像素的 TIF再高除了拖慢加载没有更多收益。压缩选 LZW。LZW 是无损压缩对 DEM 这类渐变栅格效果很好JPEG 压缩会在边界处产生块状瑕疵不适合标准成图。颜色模式工程协作默认 RGB出印刷品再选 CMYK。背景色透明如果在 ArcMap 成图里需要把数据框外的白色区域去掉勾选“背景色透明”并把“背景颜色”设置为一个不常用的纯色如 255,0,255。否则输出 TIF 在其他软件里会有白色边框盖住底图。很多第一次导出的人以为“成图 TIF”就是把这个 TIF 再导出一遍误会就大了。MXD 里引用的 DEM 是数值栅格而导出的成图 TIF 是 RGB 影像是配色后的结果不再保留高程值。也就是说成图 TIF 只用来展示汇报面积统计和坡度计算还得用第三章的输出。4.3 数据刷新的复现工作流因为 MXD 可编辑这套工作流可以反复执行。典型场景当拿到新版 DEM 或更新后的边界 shapefile不必重建 MXD只需在内容列表右键旧栅格 → “移除”再“添加数据”选新掩膜结果之后重新导出 TIF 即可。但要注意掩膜提取输出和布局导出不要写成循环嵌套即不要在脚本里把“掩膜提取结果 .tif”当作 MXD 数据源时又对这个 TIF 执行“按掩膜提取”那样会把 NoData 区域当作黑边反复叠加最后成图边界一圈黑框。为了可追溯我会把 MXD 与数据源放在同一个工作目录并保持相对路径。如果你的数据包是拷贝到另一台电脑打开“文件 → 地图文档属性”勾选“存储相对路径名”。这一步能避免路径变化带来的红感叹号属于一劳永逸的小技巧。5. 五条避坑记录从 MXD 数据源到掩膜提取的边界问题下面这些坑都是实际项目里反复出现过的按“现象 → 原因 → 解决”记下来方便对号入座。5.1 MXD 换电脑后全是红色感叹号现象把 MXD 连带数据包拷贝到另一台电脑打开后每个图层前面都是红色感叹号双击图层没有反应。 原因MXD 默认存储绝对路径。原电脑上的路径是D:\gisdata\...新电脑没有这个目录。 解决打开 MXD 后右键任意图层 → “数据” → “修复数据源”手动定位到数据包新路径。如果图层很多干脆用 arcpy 脚本把所有图层的数据源一次重定位。做法是先遍历arcpy.mapping.ListLayers(doc)用findAndReplaceWorkspacePath()把根路径统一替换。5.2 “裁剪”和“掩码提取”结果范围一样统计结果却差一截现象有人在“数据管理 → 裁剪”工具下选了“使用输入要素裁剪几何”结果看起来和掩膜提取差不多但用栅格统计出的总面积明显偏大。 原因裁剪工具即使勾选几何也只是把被裁剪区域之外的像元设为 NoData但栅格矩阵仍是矩形掩膜提取输出的是紧凑矩阵。后续统计面积时很多统计函数会跳过 NoData但由于范围不同、像元对齐不同结果自然会偏差。 解决需要统计面积或写入 ASCII 的成果一律走“按掩膜提取”用裁剪工具出来的数据只做快速浏览。识别方法很简单打开栅格属性看行数和列数。如果行列数与全国 DEM 一样说明范围和范围没缩到位掩膜提取失败或环境范围没设置。5.3 掩膜提取后 TIF 还是几个 GB文件体积没变现象明明做了按掩膜提取输出 TIF 的文件大小和原全国 DEM 差不多磁盘占用一点没少。 原因环境设置里的“处理范围”没有设置。ArcMap 默认处理范围是输入栅格全图掩膜提取只在掩膜范围外写 NoData但矩阵尺寸仍然是全国尺寸。 解决在按掩膜提取对话框点“环境”把“栅格分析 → 处理范围”改成“与图层相同”然后指定黄土高原边界图层。脚本里用arcpy.env.extent arcpy.Describe(mask_shp).extent。如果已经生成再用“数据管理 → 栅格 → 裁剪”裁掉 NoData 外框文件也会正常瘦身。5.4 全国高程 TIF 加载后 MXD 卡死转圈十几分钟现象把全国 TIF 拖进 MXD 后缩放和平移都非常卡内存占用居高不下甚至提示脚本停止响应。 原因TIF 没有构建金字塔ArcMap 每次重绘都要读取全分辨率像元。 解决在 ArcCatalog 中右键 TIF → “属性” → “金字塔” → 构建金字塔。更自动化一点用 Pythonimport arcpy arcpy.BuildPyramids_management( rD:\source\dem_wgs84.tif, PYRAMIDS, 2, BILINEAR )参数说明PYRAMIDS表示构建内部金字塔文件2是重采样级别取 1 或 2 即可BILINEAR是重采样方法DEM 用双线性比最近邻法平滑。构建完后再回到 MXD你会发现缩放流畅多了。另外把数据框坐标系固定成投影坐标系也能减少动态投影的计算压力。5.5 成图 TIF 边界出现黑框或白边现象导出的成图 TIF 在 PhotoShop 或报告里打开黄土高原边界外是一圈黑色或白色块叠在其他图上会遮挡底图。 原因导出时没有处理背景色透明。如果数据框外空白区域默认是白色黑框则通常在数据框内 NoData 区与图例透明区域冲突时出现。 解决导出地图时勾选“背景色透明”并将背景颜色设置为罕见色同时在图层属性里把 DEM 的 NoData 颜色设置为“无颜色”。之前有次我忘了勾选成图 TIF 叠加在卫星底图上边界外一大块白被客户一眼看穿这就是经验所在。6. 验证成果用 Python 快速核对 TIF 与 shapefile 的范围和投影到了最后一个环节。不管是用裁剪还是掩膜提取最终都要验证输出产品没有范围错位和投影丢失。肉眼判断在地形复杂区域很容易漏掉问题我习惯用一段短脚本在交付前检查栅格与矢量的一致性。from osgeo import gdal, ogr dem gdal.Open(rD:\loess_data\loess_dem_masked.tif) gt dem.GetGeoTransform() x_res abs(gt[1]) y_res abs(gt[5]) cols dem.RasterXSize rows dem.RasterYSize left gt[0] top gt[3] right left cols * gt[1] bottom top rows * gt[5] print(像元大小:, x_res, y_res) print(栅格范围:, left, bottom, right, top) shp ogr.Open(rD:\source\loess_boundary.shp) lyr shp.GetLayer(0) ext lyr.GetExtent() print(矢量范围:, ext[0], ext[2], ext[1], ext[3])这段代码做的验证很有意思读取栅格的仿射变换算出四个角坐标跟 shapefile 的范围对比。如果掩膜提取时的环境范围没有设置到边界自身你会看到栅格范围远大于矢量范围证明输出里面还有大块 NoData。像元大小也可以顺带检查如果原本是 30 米出来却变成 0.000333 度说明 DEM 还在地理坐标系要再确认是否满足你的分析需求。此外还要看投影字符串是否一致dem.GetProjection()和lyr.GetSpatialRef().ExportToWkt()可以打印出来对比。遇到坐标系不一致最简单的做法是在 ArcMap 里对输出的栅格做“投影栅格”转成与 shapefile 一致的投影。不要相信 MXD 数据框的动态投影那只改变显示不改变数据。我的个人记性是每跑完一步就把成果另存为带版本号的 TIF比如loess_dem_mask_v2.tif这样既能回退也能在交付前保留一条完整的可追溯链路。这种习惯在反复调试制图模板时尤其有用毕竟地理处理的“黑匣子”有时候不是报错而是悄悄让你在边界上多一条缝。希望这套流程能帮你把这个资源包真正用起来做到黄土高原空间分布与地形条件的数据可以干净、快速地落到图上。本文还有配套的精品资源点击获取