企业官网建设流程全解析

简介：宜昌市DEM数字高程数据以30米分辨率刻画地表起伏，适合GIS专业学生、地理信息行业从业者以及使用ArcGIS、QGIS等软件进行空间分析的工程师。资源覆盖宜昌市行政边界，并外扩部分周边区域，可直接用于地形渲染、坡度坡向提取、流域边界划分、可视域分析以及城市规划中的三维场景构建。压缩包内共12个文件，既有shp格式的市级边界矢量数据及其dbf属性表、prj投影文件，也有tif格式的高程栅格、tfw地理参考文件和ovr金字塔预览文件，还包含用于索引的sbn、sbx、shx及xml元数据，整体约63.67MB，文件组织规范，便于不同软件读取。目前已有824人学习下载。借助这些数据，使用者可快速完成高程建模与地形指标计算，也可将边界shp与DEM叠加开展洪涝风险、光照遮挡等专题分析，为教学实验、科研入门和区域评估提供可靠的基础底图与练习素材。

1. 一份宜昌市30米DEM数据包：从解压到坡度分析的完整落地思路

做宜昌市域的地形分析时，最烦人的往往不是算法，而是找不到一张能把整个市完整框住的DEM。网上下来的数据要么是省域大图裁不干净，要么精度低到连地名都看不出来。这份资源是一份覆盖湖北省宜昌市的DEM数字高程数据，分辨率30米，压缩包里还带了本市级范围的shp文件。有了这两个东西，你可以在ArcGIS里直接把宜昌的地形底图拉出来，做坡度、坡向、山体阴影、等高线甚至是简单的汇水分析。适合做国土空间规划、流域治理、选址评估的GIS从业者，也适合想拿真实地形数据练手的测绘学生。下面按实际使用顺序，把这份数据的关键用法和坑都过一遍。

2. DEM和shp的组合逻辑：先搞清文件构成与坐标系再动手

很多人拿到zip包第一反应是找到后缀为 .tif 的主文件，双击试试能不能打开，这没错，但很容易漏掉坐标系问题。30米DEM的数据组织比普通图片要严谨得多，主文件之外配套的元数据文件、投影文件，甚至shp边界的编码格式，都会影响后续分析是否正确。这一章先把包内文件和坐标系讲透，因为后面所有操作都建立在这两个基础上。

2.1 一个标准的30米DEM压缩包，里面应该有什么

这类高程数据常见的主格式是GeoTIFF，也可能打包成IMG或DAT，但无论是哪种，核心信息都在同名辅助文件里。打开zip先看文件结构，不要只认tif。通常一个完整的宜昌市DEM包应该包含下面这些类别，我列一个清单：

在ArcGIS中读取shp时，如果只是把.shp拖进地图，系统也会去找同名的.shx和.dbf。缺任何文件都会报“缺少Shapefile组件”的错误。不少人栽在这个细节上，以为下载不完整，其实只是漏了小文件。另外，包里可能有Readme.txt，别删，里面偶尔会写数据坐标系和采集时间等关键信息。

30米DEM的栅格主文件通常是一整块覆盖全市的大tif，但也可能按1:10000标准分幅切成若干小块，特别是从地理空间数据云下载的分幅产品。如果看到一堆带行列号的文件，比如srtm_67_05.tif，那说明需要先镶嵌再裁剪。记住这个判断逻辑，后面第3章的Mosaic操作会用到。

2.2 30米数据算什么精度：它适合哪些事，不适合哪些事

30米分辨率意味着一个像元在地面上对应大约30×30米的正方形区域，900平方米的高程取一个平均值。这个精度在宏观分析里很香，但并非所有任务都适合。做市域尺度的地形类型划分、径流方向通量计算、山体阴影渲染，它完全够用；要是拿它去找某个村庄附近的小冲沟，那结果只能说“仅供参考”。

给一个参照系，方便判断手里项目该不该用这一层数据：

宜昌地形差异很大，西部是山区，东部往江汉平原过渡，沿江一带地势平缓。30米的数据能把这种宏观格局表达清楚，但到了平坦区域，比如枝江一带，地面起伏只有几米，栅格化的高程噪声会掩盖真实微地形。所以后面第4章会专门讲平坦区域的处理。

2.3 坐标系和投影：坡度算错的最大元凶

这部分必须单独说，因为绝大多数坡度分析翻车，都栽在坐标系上。DEM的高程值单位是米，而地理坐标系（也就是经纬度）的水平单位是度。直接把WGS1984经纬度下的DEM扔进坡度工具，ArcGIS会警告水平距离与垂直距离单位不一致，实际算出来的坡度要么接近90度，要么全部是0度。

正确的做法是先把它变成投影坐标系。投影坐标系下的水平单位通常是米，这样坡度公式里的水平距离和垂直距离才统一。宜昌市经度范围约在110°E到112°E之间，UTM Zone 49N（中央经线111°E）或者CGCS2000 3度带高斯投影都很常用。原包里的.prj文件如果是GCS_WGS_1984，那就要在分析前执行投影栅格转换。

shp文件也一样。本市级范围shp可能自带投影坐标，也可能是经纬度边界。在ArcGIS里用“投影”工具把shp转换到DEM所在坐标系，保证两者坐标基准一致。如果原始数据是WGS84，而你项目的其他数据用CGCS2000，别嫌麻烦，必须统一，不然叠加后会有肉眼可见的偏移。我的习惯是第一步就打开.prj和shp的属性表源信息，先记录两个坐标系，再动手操作，这能省掉后面整块重来的时间。

3. 解压、加载、裁剪：把宜昌DEM用起来的三个必做步骤

这一章直接进入实操。从解压zip开始，到生成第一张坡度图为止，每步都有可复现的命令或工具参数。我会同时给出ArcGIS菜单操作和GDAL命令行两种路径，日常用哪个顺手就选哪个。

3.1 解压与文件检查：别把路径搞出中文

从Windows资源管理器里右键解压是最简单的方式，但遇到数据文件名带中文，或者路径层级太深时，ArcGIS偶尔会不识别。更稳的办法是用命令行解压，可控编码。如果数据包文件名是英文，直接解压；如果包含中文目录名，建议用Git Bash或WSL执行：

unzip -O gbk 宜昌DEM30m.zip -d ./yichang_dem

-O gbk指定文件名按GBK编码解压，解决Windows中文文件名在Linux环境下解压乱码的问题。-d ./yichang_dem把输出目录固定为当前目录下的yichang_dem，方便后面对照路径。解压完成后进入目录，用ls -lh查看文件大小，确保tif不是0字节或几十K的破损文件。

如果你的环境没有unzip，用7-Zip也可以，但注意解压后把顶层文件夹名改成纯英文，例如Yichang_DEM_30m。中文路径在ArcGIS的历史版本里出现过“栅格数据集无法创建”的报错，英文路径能绕开这个坑。

3.2 在ArcGIS中加载并检查图层：先看范围再看值域

打开ArcMap或ArcGIS Pro，把tif拖进地图，再把shp拖进去。第一件事不是缩放，而是右键图层，点击“属性→源”，看空间参考。确认DEM和shp的坐标系是否一致。如果一致，继续；如果不一致，先记录。然后右键DEM图层选择“缩放至图层”，用目测对比shp是否在DEM范围内。

如果DEM是多幅小tif，ArcGIS加载进来会各自为政，此时需要用“镶嵌至新栅格”工具：搜索该工具，输入所有tif，像素类型选择“16 bit Signed”（SRTM类数据常见）或“32 bit Float”，波段数填1，镶嵌运算符选BLEND。输出后用“复制栅格”工具再生成一份干净的单波段文件。这一步能消除接缝和重叠区域的颜色跳变。

这种检查流程别跳过。我有一次拿到某省30米DEM，加载后看起来和shp完全套合，结果投影一看，DEM是GCS_WGS_1984，shp是CGCS2000 3度带，两者在平原地区偏移了将近80米。因为屏幕缩放比例小看不出来，最后在公里网格上才暴露。所以源头信息比眼睛可靠。

3.3 用shp裁剪DEM：提取宜昌市域范围

加载无误后，需要把DEM裁剪到shp边界内，这样能剔除周边区域数据，文件体积变小，分析速度也快。在ArcGIS里搜索“按掩膜提取（Extract by Mask）”，这是空间分析模块的标准工具。参数设置如下：

• 输入栅格：宜昌DEM（如果是多幅，用镶嵌后的tif）
• 输入掩膜：宜昌市shp
• 输出栅格：指定路径，命名如Yichang_dem_cut.tif

环境设置里有一个关键点：在“地理处理选项”的环境页签中，把“处理范围”设为“与掩膜相同”，输出的像元大小设为与输入DEM一致（30米）。否则裁剪结果可能四边多出半个像元，或者被重采样成别的分辨率。

不习惯用ArcGIS工具箱的话，用GDAL一行命令也能完成，而且确定性更强：

gdalwarp -cutline Yichang.shp -crop_to_cutline -dstnodata -9999 -of GTiff Yichang_dem.tif Yichang_dem_cut.tif

-cutline指定边界shp，-crop_to_cutline表示仅保留边界内部，-dstnodata -9999把边界外的高程设为空值，-of GTiff指定输出格式。这条命令在QGIS内置终端里也能跑。注意shp的坐标系必须与DEM相同，否则gdalwarp会先做投影转换，可能产生偏移。

3.4 生成坡度和坡向图：一条标准生产流

坡度是DEM最常用的派生数据。在ArcGIS中找“坡向”旁边的“坡度”工具（Spatial Analyst → 表面分析 → 坡度），输入前一步裁剪好的tif。输出参数里，有两种单位可选：度和百分比。一般分析用度，水土保持计算里常用百分比。宜昌山地坡度普遍较大，度更直观。Z因子保持1，因为投影坐标下水平和垂直都是米，不需要缩放。如果输入的是经纬度坐标系，ArcGIS会提示使用0.000009左右甚至更大的Z因子，看到这个提示就别继续了，先重投影。

生成坡度图后，默认灰色符号系统很难看。点击符号系统，选择“拉伸”类型，色带选“高程配色”或自己搭一版黄绿棕分带。进一步把连续坡度栅格重分类成标准坡度级别：≤5°、5-15°、15-25°、25-35°、>35°，用“重分类”工具保存成整型栅格，方便后续出图或面积统计。

4. 30米DEM的精度边界：拿到手先做这三件验证

数据能显示不等于数据可靠。30米DEM在生产过程中会有无效值、内陆凹陷、平坦区域噪声等常见瑕疵。刚拿到宜昌数据时别急着做分析，先花十分钟做三件验证。

4.1 高程值域检查：区分无效值与真实低地

用“获取栅格属性”工具查看DEM的最小值、最大值。宜昌市海拔范围大约从沿江的几十米到西部山区的两千米上下，如果最小值出现-32768、-9999这类极端负值，说明该区域是NoData。在后续坡度计算时，NoData区域会沿边界扩散。

处理方式是先把NoData统一成固定空值。ArcGIS里有“设置空值”工具，可以把所有非有效值（比如小于某个阈值，或等于-32768）改为空。也可以用栅格计算器写条件式：

Con(IsNull("dem") | ("dem" < 0), -9999, "dem")

逻辑说明：Con是ArcGIS栅格计算器里的条件函数，第一段是判断条件，第二段是满足该条件时输出的值，第三段是其他情况沿用的原值。这里把NoData或小于0的像元统一赋成-9999，便于后续统一处理。注意宜昌沿江部分地区海拔确实能低到几十米，所以阈值不要设成“低于50米就置空”，只处理负大值和NoData。

4.2 叠加矢量边界检查位置套合度

验证位置精度最直接的办法，是把DEM和高分辨率影像或者已有的河流shp叠在一张图里。在ArcGIS里加载宜昌的河流水系shp（很多公开地理数据平台都有），叠加到DEM上放大到1:10万左右，观察河流线是否沿着山谷的最低处走。如果河流线往山坡上跑，或者明显偏离沟谷中线，说明笛卡儿坐标基准有问题。

这个偏移可能来自投影转换参数。比如WGS84测定的坐标与CGCS2000之间存在几十厘米到几米的系统差，在30米分辨率下肉眼几乎不可见，但沿山脊线对比时能发现。如果偏差大于一个像元（30米），直接把DEM和shp都转换到同一个CGCS2000投影下再重新裁剪。位置验证这步不复杂，但能避免你把一套错位数据用到底。

4.3 平坦区域的噪声陷阱

城区和水田是30米DEM最容易出噪声的地方。宜昌主城区在现代化建设中，许多地块被推平，但DEM像元里还保留着改造前的起伏痕迹；农田则因为灌溉季积水，高程值跳跃明显。这些噪声在坡度分析中会表现为短促的、毫无规律的坡度突变，比如平地上突然冒出一个20度的像元。

针对这种情况，我一般会先做一个“焦点统计（均值）滤波”，把每个像元周围3×3邻域的高程平均值赋给中心格，然后再算坡度。在ArcGIS里：搜索"焦点统计"，邻域设为矩形3×3，统计类型MEAN，输出后替代原DEM。注意不要重复多次，两次以上会过度平滑，陡坡的真实转角会被抹掉。

另外，填洼（Fill）也是必然步骤，尤其做汇水分析时。SRTM和ASTER类数据在河谷区域经常出现伪洼地，不填掉的话流向计算会出现断层。ArcGIS的“填洼”工具只需输入DEM，主要参数“填充阈值”保持默认即可；如果生成的流向图还有断头，再手动提高阈值，但要谨慎，每次只加10米。

5. 避坑指南：DEM和shp一起用最容易翻车的五个场景

这一章写我遇过的高频故障，全部都是“现象→原因→解决”结构。对照排查，能省下不少刷论坛的时间。

5.1 现象：DEM加载后一整片黑，什么都看不到

原因：栅格图层的默认拉伸方式可能把有效高程全映射到了黑色端；或者NoData值被当成有效值参与渲染。常见于原始包里的.aux.xml统计信息缺失，ArcGIS临时计算了错误的直方图。

解决：右键图层→属性→符号系统→拉伸→拉伸类型改为“百分比截断”（或“标准差”），数据裁剪到10%-90%区间，色带重选。如果还是黑，就重新生成统计信息：在目录窗口中找到tif，右键→属性→源→计算统计信息。注意检查值域，正常30米DEM的高程单位是米，不是厘米；如果tif是浮点型且值域是两三千米，拉伸百分比选15%-85%一般能正常显示。

5.2 现象：shp拖进ArcGIS提示“没有可显示的数据源”，或只显示一部分

原因：shp配套文件缺失，或者shp里的.prj定义与内容不匹配。很多下载站打包时漏了.shx或.dbf，ArcGIS找不到坐标和属性索引就打不开shp。

解决：检查压缩包中是否有.shp、.shx、.dbf、.prj、.cpg这五个文件。缺shx，用其他工具重建（如QGIS里的“修复数据源”），但更省事的是回到下载源重新获取或者在末端补充属性信息。另外，路径中不能有中文和空格，把整个目录复制到纯英文路径下重试。

5.3 现象：坡度图看起来全是0度和90度，没有中间过渡

原因：这是典型的坐标系问题。直接用经纬度坐标系DEM算坡度，ArcGIS内部把水平距离当作度，垂直距离当米，计算出的人工比值就变成极端值。

解决：先对DEM执行“投影栅格”工具，输出坐标系选择UTM Zone 49N或CGCS2000 3度带高斯投影，再跑一次坡度工具。投影栅格参数里“重采样技术”选双线性（Bilinear），“像元大小”保持原样（30）。别选最近邻，虽然快，但会让台阶状地形更明显。

5.4 现象：DEM与shp边界错位，相差几十到几百米

原因：两个数据的坐标系基准不一致。DEM是WGS84，shp是CGCS2000；或者DEM是经纬度，shp是投影坐标。还有一种隐蔽情况：shp的.prj写的是北京54，但实际坐标是西安80系，如果是第三方转换过的边界文件，建议复查。

解决：右键shp→属性→源，先看坐标。然后用“投影”工具把shp转换到DEM的坐标系。如果两个都是地理坐标但基准椭球不同，就需要用“创建自定义地理变换”或者直接重新下载套合过的边界shp。对30米分辨率来说，WGS84和CGCS2000的差异在毫厘级，肉眼看不出；大幅错位一定是因为一个用了投影坐标一个用了经纬度。

5.5 现象：多幅DEM拼接后出现明显的接缝色差或断崖

原因：不同图幅之间的采集时间、传感器标定或后处理流程不一致，导致接缝位置的高程坡度突变。拼好后直接做坡度分析时，接缝处会出现一条带状“伪山脊”。

解决：“镶嵌至新栅格”工具的镶嵌运算符选BLEND，系统会在重叠区域做渐变色加权平均。如果接缝还明显，把拼接结果再做一次3×3中值滤波，强度不要太大。更专业的做法是用“按颜色镶嵌”，但这只针对色彩影像，DEM按高程值混合即可。经验是：先镶嵌再滤波，滤波接缝区域比滤镜整个区域好得多。

6. 把30米DEM榨干：山体阴影、等高线与地形剖面

到这章，数据已经能正常用了，但30米分辨率的价值不止坡度坡向那么简单。三个进阶操作能让它产出更接近成果图的素材。

山体阴影（Hillshade）是制图里最常用的底图。ArcGIS里搜索“山体阴影”，输入裁剪后的DEM，方位角设为315°，太阳高度角设为45°，Z因子默认1。输出之后把它放在原始DEM图层上面，在图层属性里把透明度调到40%～60%，再叠上宜昌市界shp和道路线，就是一张能直接出制图效果的地形底图。投影坐标系下的山体阴影不需要改Z因子，但如果你在WGS84经纬度下硬算，阴影方向会错乱。

等高线提取方面，在ArcGIS里用3D Analyst的“等高线”工具。输入DEM后，等距设置需要根据地形决定。宜昌西部山区地势陡峭，设20米到50米能得到合理线间距；东部平原设5米到10米，但30米DEM的折线会出现大量锯齿。所以我习惯先提取原始等高线，再加一个“平滑”工具，方法选PAEK算法，容差设30米。平滑后的等高线再裁到shp范围内，导出为dxf，就能在CAD里继续编辑。需要提醒的是，别试图从30米DEM里提取2米等高距，那是在考验数据分辨率，结果只会是一团麻线。

地形剖面图是检验数据是否能用的一招妙棋。在ArcGIS Pro中，用“创建剖面图”工具，在DEM上画一条穿越宜昌市南北的线，比如从长阳山区延伸到枝江平原，系统会立刻生成一条高程曲线。如果曲线平滑过渡，说明数据质量尚可；如果出现频繁相邻落差超50米的锯齿，那就要怀疑原始数据是不是碰到了数据空洞，回到前面第4章补预处理。

这让我想到一个习惯：每次拿到新的栅格数据，我都强制走一遍“坐标系查看→值域检查→范围核对→噪声滤波”的全套流程，哪怕只是用来出一张示意图也绝不跳过。等习惯成自然，你会发现很多“数据坏了”的问题，其实是加载顺序和数据源没理清。希望这些积累能帮你在宜昌地形数据上少走几步弯路。

本文还有配套的精品资源，点击获取