简介覆盖湖北省宜昌市全域及周边部分区域这套30米分辨率数字高程模型DEM数据包附带宜昌市行政边界Shapefile适合GIS初学者与专业人士用于地形分析、环境模拟和地图制图。压缩包共12个文件总大小约63.67MB核心包含宜昌市dem.tif高程栅格、宜昌市范围.shp矢量边界并配套prj投影文件、tfw世界文件、xml元数据及sbn/sbx空间索引可被ArcGIS、QGIS等软件直接加载兼顾快速预览与坐标配准。数据同时提供dbf属性文件与XML说明便于核对地理参考和属性字段。目前已有825人学习下载该数据经过实际应用验证可直接用于坡度计算、流域划分、可视性和淹没分析等教学或科研场景为城市规划、环境研究和灾害风险评估提供可靠的高精度地形底图。1. 宜昌市30米DEM数据拿到手先别急着出图宜昌的地形不简单西部山区沟谷交错中部丘陵过渡东边才是沿江平原。做地形分析最耗时间的环节往往不是算法本身而是数据准备。这份宜昌市DEM数字高程数据30m含本市级范围shp文件压缩包把全市30米分辨率DEM栅格和本市级行政范围SHP文件放在一起省去了拼接分幅、自己画行政边界这一步。30米分辨率在市级尺度上看坡度、坡向、山体阴影够用文件体积也适中ArcMap打开不卡。适合GIS从业者、测绘相关专业学生以及正在做宜昌市域地形研究但被数据预处理卡住的人。拿到手后别急着直接拉伸渲染先做三件事确认坐标系、确认裁剪范围、确认有效高程值。2. DEM数据选型与宜昌市数据构成30米分辨率、SHP面层与坐标系的三个判断2.1 30米分辨率对宜昌山区地形的适用边界先算一笔账。30米分辨率意味着一个像元对应地面上约900平方米。宜昌西部山区的冲沟、陡坎在30米像元里几乎都会被平滑掉但长江河谷、主要支流、山体走向这种百米级地形特征还能保留下来。所以这份数据最适合的场景是市级尺度的栅格分析土地利用的坡度分级、山洪淹没范围的初步圈定、城市开发的适建性评价。如果拿到村一级或者要给某一条具体线路做土方估算30米是明显不够的至少要到10米甚至5米以下。我见过一些新手拿30米DEM直接做小流域汇水分析结果集水面积和实测对不上问题不在地貌而在分辨率。一个30米像元的误差放到坡度算法里会被放大成好几度的偏差。反过来如果一上来就追求5米、2米高程数据成本和时间都会成倍增加而且对全市域范围来说数据处理慢到让人失去耐心。顺便补充一个现实选择现在用BigEMap这类工具可以下载到5米精度高程但那是商业在线数据源下载量受权限和网络限制。对于宜昌全市域范围的预研分析30米DEM是成本和精度的平衡点。分析场景建议分辨率典型用途宜昌市域/县域宏观分析30m坡度分级、山体阴影、淹没模拟初判乡镇/小流域分析10m-12.5m汇水分析、灾评、生态敏感性线路/场地级设计1m-2m土方量、边坡、工程可行性2.2 数据包里的双层结构DEM栅格与本市级SHP的配合关系这个压缩包的核心是两部分一个DEM栅格TIFF文件一个宜昌市本市级范围的SHP面文件。SHP不是用来做属性关联的它的唯一使命是当裁剪掩膜把DEM从大范围数据里切到宜昌市行政边界内。打开压缩包后建议先列一下目录结构。常见做法是用Python快速看一眼栅格和矢量各自的信息import arcpy arcpy.env.workspace rD:\yichang_dem print(arcpy.ListRasters()) print(arcpy.ListFeatureClasses())逻辑很简单确认当前目录里有哪些栅格和要素类。列出之后再检查SHP的空间范围和字段desc arcpy.Describe(yichang_city.shp) print(范围:, desc.extent) fields arcpy.ListFields(yichang_city.shp) for f in fields: print(f.name, f.type)这段代码会打印SHP的XMin、YMin、XMax、YMax和字段列表。我一般会先记住四个范围值后面裁完DEM再对比一次。SHP通常带字段比如地级市名称或行政区划代码但这里用不上裁剪只认几何不认属性。2.3 坐标系统一CGCS2000在ArcMap中的识别与转换拿到数据后第一个坑往往不是裁剪而是坐标系。DEM可能来自不同数据源栅格本身可能没有投影信息只有经纬度SHP可能是CGCS2000高斯投影也可能是WGS84。宜昌的经纬度大致在110°15′E到112°04′E、29°56′N到31°34′N之间这类区域用UTM 49N或者CGCS2000 3度分带都合适。我在处理这类数据时的习惯是先右键图层看属性源标签里看空间参考。如果DEM显示的是GCS_WGS_1984或GCS_CGCS_2000这样的地理坐标系就得先投影再裁剪。否则后续算坡度时会被经纬度坐标的尺度问题带偏。sr_proj arcpy.SpatialReference(WGS 1984 UTM Zone 49N) arcpy.ProjectRaster_management( yichang_dem_orig.tif, yichang_dem_utm49.tif, sr_proj, BILINEAR )参数解释第一参数是输入栅格第二是输出路径第三是目标投影坐标系第四是重采样方式。DEM是连续表面数据重采样选BILINEAR或CUBIC不要选NEAREST否则会在地形上产生阶梯状断层。有一种情况要区分如果栅格显示已经有投影坐标系但SHP却没有.prj文件ArcMap会直接弹警告说缺少空间参考信息。这时候不要用Project工具强行转换而应该先用Define Projection把SHP的坐标系定义成和DEM一致再投影。定义和投影是两个完全不同的操作前者只是给数据一个身份标签后者才是真正改变坐标数值。提示裁剪前花两分钟做坐标系检查比裁剪后发现图对不上再返工划算得多。3. ArcMap裁剪DEM的两种路线按掩膜提取与面裁剪的结果差异3.1 按掩膜提取Extract by Mask的操作路径与参数ArcMap里最常用的是Spatial Analyst工具下的Extract by Mask对应热搜词里那句“依靠面图层掩码提取”。操作路径是ArcToolbox - Spatial Analyst Tools - Extraction - Extract by Mask。工具界面里只有三个关键参数输入栅格选DEM输入栅格或要素掩膜数据选宜昌市SHP输出栅格给一个路径。这里有个容易忽略的点环境设置里的“捕捉栅格”。如果环境里没有指定输出栅格的像元对齐方式默认跟随输入DEM也就是说裁剪结果不会因为SHP边界偏移而重建网格这是好事保证了后续分析的一致性。用ArcPy的命令长这样from arcpy.sa import * dem yichang_dem_utm49.tif mask yichang_city.shp out ExtractByMask(dem, mask) out.save(yichang_dem_mask.tif)逻辑是遍历SHP范围内每个DEM像元像元中心点落在面内的保留落在面外的置为NoData。所以裁剪边界不会是绝对平滑的行政边界线而是由像元中心点决定的锯齿边界。这是正常行为不是数据错误。若想边界更贴合SHP可以先把SHP向外或向内做微小的Buffer但这会改变分析范围要谨慎。3.2 用SHP面图层直接裁剪Clip的设定与边界条件另一条路线是Data Management工具下的Clip。很多用户习惯点这个因为名字直观。但Clip有个关键开关Use Input Features for Clipping Geometry。勾选时输出按SHP的实际几何裁剪不勾选时输出只是SHP总外接矩形的范围周围会留下大量NoData区域。arcpy.Clip_management( yichang_dem_utm49.tif, #, yichang_dem_clip.tif, yichang_city.shp, 0, ClippingGeometry, NO_MAINTAIN_EXTENT )参数逐个说第二参数rectangle传#表示不手动指定外接矩形完全交给SHP第四参数是裁剪要素第五参数nodata_value传0表示把原来无效值写成0实际使用时可以传空字符串让它保留原有NoData第六参数ClippingGeometry是关键不写它或写成NONE输出就是外接矩形第七参数控制是否保持原栅格范围我一般用NO_MAINTAIN_EXTENT让输出范围贴合裁剪结果。从性能和结果上比Clip通常比Extract by Mask快文件体积也更小因为它直接按几何取像元不需要重新计算掩膜内部关系。但Clip处理边缘像元的规则与Extract by Mask相近都是按像元中心点判断所以边界差别不大。3.3 两种方式的实际差异边缘像元、NoData与文件体积不少初学者会问掩膜提取和面裁剪到底有啥区别用一张表可以讲清楚。维度按掩膜提取面裁剪勾选要素几何面裁剪不勾选输出形状与SHP几何一致与SHP几何一致SHP外接矩形SHP范围外全部NoData取决于NoData设置保留NoData或0边缘像元像元中心判断像元中心判断矩形边界截断速度较慢较快最快典型用途面积统计、后续栅格分析快速出图、文件瘦身临时查看实际项目里我处理这份宜昌市DEM时更倾向于用Extract by Mask因为它对面积的统计更干净。而如果只是想和图斑叠个底图用Clip更快。无论选哪条都要注意一个暗坑SHP里如果存在多个面要素比如本市级范围SHP包含主城区、飞地、岛屿多个面两种工具都会按所有面的并集处理。如果数据包里SHP还附带市级以下细分面裁剪时不会自动按细分区县分别输出必须用Split by Attributes或循环迭代。4. DEM裁剪避坑与排查五个最常见翻车现场4.1 现象裁出来的DEM大半个屏幕是黑色NoData占比极高原因SHP范围超出了DEM有效数据范围或者DEM数据本身在边界处有无效值掩膜提取后无效区域被原样保留成了NoData。这个现象在分幅拼接的DEM里非常常见因为每幅图边缘都有拼接缝。解决先对比SHP范围与DEM范围如果两者差异很大说明SHP不是从DEM上切下来的需要换范围更贴合的数据源。如果只是边缘一圈NoData可以在环境设置里打开“范围”并选择SHP把处理范围锁定。对剩余NoData区域可以用栅格计算器里的Con函数填充相邻值但要注意这会人为修改地形数据除非必要不建议做。out_nodata Con(IsNull(yichang_dem_clip.tif), 0, yichang_dem_clip.tif)Con函数的逻辑是第一个参数是条件第二个是条件成立时的值第三个是条件不成立时的值。这里把NoData像元填成0能让表面看起来完整但对坡度分析极不友好0米高程会和周边真实高程形成假陡坎。所以这个操作只适合出图预览不适合定量计算。4.2 现象DEM在ArcMap里渲染成一整片刺眼的彩色马赛克原因符号化方式用了分类渲染默认分5类把高程值切成几个区间山区地形细节全被抹平了。这是新手最容易犯的错不是数据问题是显示问题。解决在图层面板里把符号系统改成“拉伸”类型选“标准偏差”或“最值”再配一个适合高程的色带。宜昌这种山地地形我习惯把拉伸类型选为“最值”色带用从深绿到棕再到白的渐变。拉伸显示的细节远比分类渲染丰富。如果还觉得不够叠加一层山体阴影作为透明底图把DEM主图透明度调成50%左右地形立体感马上出来。4.3 现象SHP和DEM在ArcMap里显示位置错开一个在东一个在西原因没有对齐坐标系。最常见的是SHP有.prj文件但DEM没有或者两个坐标系基准面不一致。宜昌范围内WGS84与CGCS2000的差值在平面上看只有几十厘米到一米量级但如果一个是投影坐标一个是地理坐标画出来就是十万八千里。解决先分别看两个图层的空间参考。把SHP投影到和DEM相同的坐标系或者反过来。如果SHP缺少.prj用Define Projection指定不要用Project。方向错了数据就废了。4.4 现象高程值比真实地面普遍高或低几十米原因DEM生产时采用的高程基准和当地常用高程基准不一致。国内常见的DEM有基于EGM96大地水准面的也有基于1985国家高程基准的两者在某些山区会差十几米到几十米。对地形形态分析比如坡度、坡向这个整体偏移不影响结果但如果你要和高程控制点比对或者计算水位淹没范围就必须校正。解决找宜昌市内几个已知高程点统计DEM提取值与实测值的差值。如果差值比较恒定在栅格计算器里直接加一个常数就行。4.5 现象ArcMap在加载或裁剪DEM时直接闪退原因ArcMap默认是32位进程单个栅格过大时内存不够用。30米分辨率的全市范围栅格本不该这么大但如果没有建金字塔加载时会尝试把整个栅格读入内存很容易爆。解决加载前先在目录里给TIFF生成金字塔或者在ArcToolbox里用Build Pyramids工具扫一遍。还有一招是把TIFF转成CRF格式CRF对栅格运算的支持更好读取也更快。arcpy.BuildPyramids_management(yichang_dem_utm49.tif, PYRAMIDS, , BILINEAR)参数说明第二个参数PYRAMIDS表示创建金字塔第三、第四参数控制压缩方式和重采样方法。BILINEAR保证缩小显示时地形平滑。5. 把DEM变成实用地形产品山体阴影、等高线与坡度一次生成5.1 山体阴影Hillshade的参数与显示效果山体阴影是DEM最直观的地形可视化方式。工具默认的方位角是315度也就是阳光从西北方向照过来高度角45度。宜昌西部山区沟谷多为了突出峡谷细节我通常会把高度角降到30度到35度之间。高度角越低阴影越长地形越立体但太低了会丢失沟谷底部的细节。hillshade Hillshade(yichang_dem_clip.tif, 315, 35, 1) hillshade.save(yichang_hillshade.tif)参数解释315是太阳方位角35是太阳高度角最后的1是Z因子也就是垂直方向缩放系数。Z因子不是随意填的——当DEM坐标单位和高程单位都是米时Z因子为1如果坐标系是经纬度Z因子要填约0.000009但这就是治标不治本正确做法是先转投影坐标系再做分析。我给宜昌市区做过一个地形制图把山体阴影透明化处理压在坡度图下面效果比纯色渲染好很多。这个叠加顺序很关键山体阴影放上层透明度调到60%到70%下面的高程或坡度图就能保持色彩信息又有立体感。5.2 等高线生成间隔选择与平滑处理从DEM生成等高线在ArcToolbox里对应Contour工具。间隔怎么选取决于展示尺度市级总体规划图我通常用50米间隔县级或重点片区用20米河谷平坝区域用10米。宜昌市区从江边高程约50米到周边山区1000米以上20米间隔生成的等高线数量已经相当可观。contour Contour( yichang_dem_clip.tif, yichang_contour_20m.shp, 20, 0 )参数解析输入栅格为DEM输出是线要素类20是等高线间隔0是起始等高线高程。输出线要素的高程字段会自动命名为CONTOUR。如果等高线毛刺太多不要急着用Generalize工具把线平滑到面目全非那是自欺欺人。30米DEM本身就有随机误差毛刺是真实反映。想减少毛刺可以先对DEM做一次3×3的Focal Statistics低通滤波把高频噪声压下去再生成等高线。滤波必然损失真实细节鱼与熊掌要自己权衡。5.3 坡度分析与宜昌复杂地形的适配坡度工具有两个单位度DEGREE和百分比PERCENT。宜昌西部山区坡度陡坡面很多超过25度用度展示直观但如果做水土流失评价标准阈值通常用百分比比如15%、25%、35%。这个细节很多人不看直接用了默认的度结果和标准对不上。slope Slope(yichang_dem_clip.tif, DEGREE, 1) slope.save(yichang_slope_deg.tif)参数说明DEGREE表示输出坡度单位为度1依然是Z因子。生成之后按经验阈值做重分类0到5度为平地5到15度为缓坡15到25度为中等坡25到35度为陡坡大于35度为极陡坡。宜昌山区在35度以上的区域不在少数特别是三峡沿岸。最后提醒一句坡度分析对DEM噪声高度敏感。如果最终坡度图上有明显条带或方块效应先回头检查DEM原始数据是不是有拼接痕迹而不是换更复杂的算法。6. 裁完DEM之后用像元数验证数据完整性的一个习惯数据裁剪完不是立刻做分析先算一笔账。把有效像元数乘以像元面积得到的结果应该和SHP面积接近。这个过程我每次都会走一遍花不了两分钟但能避免灾难性的返工。desc arcpy.Describe(yichang_dem_clip.tif) cell desc.meanCellWidth total int(arcpy.GetRasterProperties_management( yichang_dem_clip.tif, CELL_COUNT ).getOutput(0)) nodata int(arcpy.GetRasterProperties_management( yichang_dem_clip.tif, NODATA_CELL_COUNT ).getOutput(0)) area_km2 (total - nodata) * cell * cell / 1e6 print(有效面积约:, round(area_km2, 1), km²)这段代码先从栅格属性里拿像元尺寸再拿总像元数和NoData像元数两者相减就是有效像元数乘以像元面积换算成平方公里。值得注意meanCellWidth只有在投影坐标系下才是直接用米表示的像元宽度如果是经纬度坐标这个值会是度算出来的面积就没有意义了。然后取SHP的SHAPE_AREA字段做对比with arcpy.da.SearchCursor(yichang_city.shp, [SHAPE_AREA]) as cursor: for row in cursor: shp_area_km2 row[0] / 1e6 print(SHP面积:, round(shp_area_km2, 1), km²)相对误差在3%以内说明裁剪结果可信。边界像元的中心点取舍会让面积有小幅偏差这是正常的。但如果误差超过3%就要回头查SHP范围、DEM范围、投影是否统一这几件事。这个习惯是我在给三峡库区周边做山地灾害分析时被逼出来的。当时裁完的DEM在等高线图上少了一条完整山脊线数据表面看着正常等算到汇水面积才发现偏差有8%。从那以后我每裁一个DEM都会强制走一遍这个像元数验证流程。数据是别人处理的但最后为结果负责的是自己。希望这个验证流程也能帮到你。本文还有配套的精品资源点击获取
