简介这份四川省成都市12.5米DEM数据包面向GIS、测绘、城乡规划、环境与灾害研究等从业者提供比常用30米分辨率更精细的高程栅格能更清晰呈现山脊、山谷、坡度等地形细节适合开展地形分析、坡度坡向提取、城市三维可视化与灾害风险评估。压缩包共12个文件总大小约80.13MB核心为TIFF格式高程栅格并配套shp/shx行政边界、prj投影定义、tfw坐标定位、dbf与vat属性表及ovr金字塔等文件分别用于定义坐标系、精准配准地理位置、存储属性信息和加速大图缩放可在ArcGIS、QGIS中直接加载。目前已有419人学习下载适合需快速获取成都高精度地形底图的用户使用。这份数据能省去自行下载、裁剪、拼接多源数据的流程一次拿到完整市域DEM与行政边界并支持淹没模拟、通视分析、工程选址、生态评价等进阶应用显著缩短前期数据处理时间。 在GIS数据交换群里泡久了经常会看到一类资源地名开头、分辨率居中、后缀还带个.7z。“四川省成都市12.5米DEM数据含市级范围shp文件.7z”就是典型代表。乍一看只是一个压缩包但对做过成都方向地形分析的人来说这套东西能省下大半天的预处理时间。12.5米DEM意味着每个像元对应地面12.5米见方比常见的SRTM 30米、ASTER 30米细腻不少山区沟谷、城区微地形都能看出轮廓附带市级范围shp文件相当于把裁剪边界和行政范围都给你备好了。这篇文章就围绕这套数据讲清楚数据价值、解压校验、shp边界该怎么用、裁剪怎么做以及如果数据不够用自己上哪儿重新下载还有什么坑要避开。无论你是做规划的、搞测绘的还是想用地形数据做可视化的小白都能从里面找到能直接抄走的操作。1. 12.5米DEM到底强在哪先搞清这个数据包的含金量1.1 DEM不是图片是一张高程数值表很多人第一次拿到12.5米DEM双击打开发现是一张灰蒙蒙的图以为就是普通卫片。其实DEM的每个像元记录的不是颜色而是这一小块地面的海拔高度。你可以把它理解成一张“海拔查询表”只不过这张表是栅格结构的每个格子都有明确的经度、纬度和高程值。12.5米分辨率指的是一个像元对应地面12.5米乘12.5米的实际范围。以成都市域约1.4万平方公里的体量来算一张完整的12.5米DEM格网数量至少在八九千万量级解压后动不动就是几个GB的GeoTIFF文件。这个数据量放到以前普通电脑加载都会卡现在倒是轻松很多但压缩包体积依然可观所以用7z打包几乎是标配。1.2 同为全球DEM12.5米和30米差在哪市面常见的开源DEM主要就那几款SRTM 30米、ASTER 30米、ALOS 12.5米还有商业的TanDEM-X 12米/30米。很多人只下载过30米觉得够用了但真正做过山区项目的人普遍会觉得30米在成都这种地形起伏大的地方不够看。数据源分辨率覆盖范围特点常见用途SRTM 30米约30米全球北纬60°到南纬56°覆盖广、下载方便但局部细节平滑区域宏观地形、水系提取ASTER 30米约30米全球易受云和条带影响垂直精度波动较大教学演示、大范围概览ALOS 12.5米约12.5米全球大部分陆地区域细节更突出雷达数据受天气影响小地质灾害、精细坡面分析TanDEM-X 12米约12米全球精度高但商业授权费用不低专业测绘、工程前期选择12.5米数据的核心原因是它能看出30米数据看不出来的地形细节。成都西边是龙门山一带的高山峡谷东边有龙泉山中部是平原。在山区12.5米数据能把冲沟、山脊转折、崩滑体边界这些微地貌显示得更清楚在平坝区域河道阶梯、人工台地、城中村的地形起伏也能有所体现。对地质灾害解译、洪水淹没模拟、机场净空评估、管线选线这类工作来说这种细节很关键。1.3 带shp文件为什么是加分项单独一份DEM对你来说只是一张高程图但如果配套一个成都市级范围的shp边界数据就变得好操作得多。你不必再从全国行政区划数据里手工提取边界也不用担心裁剪时把周边地市带进来。更实用的是shp文件可以做掩膜分析比如统计成都市域内的平均海拔、坡度分级面积或者生成成都范围的等高线、山体阴影图。2. 打开.7z的正确姿势解压、校验与文件清单核对2.1 GIS数据打包为什么偏爱7zGeoTIFF格式本身带有大量重复的块结构和元数据压缩空间很大7z算法能压出比zip更小的体积。更重要的是shp文件不是单个文件它必须由.shp、.shx、.dbf、.prj等多个文件共同组成。把整个shp文件夹用7z压成一个包传输、分享时就不会出现“只发了.shp别人打不开”的尴尬情况。所以“为什么我发的shp对方打不开”这类问题多半是没连其他后缀文件一起发。解开这个压缩包之前先养成一个习惯先弄清楚里面装了什么。用7z l可以不解压直接列出文件清单。这样你还没解压就能知道里面是一个完整的成都范围的tif还是按区县分幅的多个tif便于提前规划后续操作。2.2 Windows和Linux下的解压操作Windows下最简单装一个7-Zip或者Bandizip右键选择“解压到当前文件夹”即可。一个小提醒解压路径尽量不要带中文也不要把压缩包放在网盘同步目录里直接解压某些GIS老版本软件对中文路径和云盘占位文件不友好加载时会莫名其妙报错。Linux服务器上解压很多人的第一反应是unzip但7z格式需要p7zip组件# Debian/Ubuntu系 sudo apt install p7zip-full # CentOS/RHEL系 sudo yum install p7zip安装完成后解压7z x 四川省成都市12.5米DEM数据含市级范围shp文件.7z如果文件名太长或者带中文直接Tab补全即可。7za是7z的独立命令行版本功能上对解压来说差不多日常用哪个都行。2.3 解压前先做哈希校验别偷懒网盘或FTP传文件偶尔会出现下载不完整的情况。一个好几个GB的数据包如果中间坏了一个字节解压时可能只报一句“CRC Failed”然后整包报废。可靠的做法是下载完后用SHA256校验一下压缩包的哈希值和发布方提供的值比对。Linux下sha256sum 四川省成都市12.5米DEM数据含市级范围shp文件.7zWindows PowerShell下Get-FileHash .\文件名.7z -Algorithm SHA256如果发布方没有提供哈希值也可以用7z t测试压缩包完整性7z t 文件名.7z这个命令会逐个文件测试解压数据是否正确。压包测试通过再解压、再投入项目能避免很多隐藏问题。我自己的习惯是所有GIS原始数据下载完先记录哈希值再归档备份方便以后核对文件有没有被误改。3. 市级shp边界文件拿到手先确认坐标系再谈裁剪3.1 shp的一整套文件缺一个都不行解压后你会看到一堆扩展名不同的文件。.shp存几何图形.shx是空间索引.dbf存属性字段.prj记录坐标系还有可能附带.sbn、.sbx等辅助文件。很多新手只把.shp拖进GIS结果属性表是空的或者坐标系显示未知就是这个原因。用7z打包数据通常不容易缺文件但解压后如果打算把shp单独发给别人建议继续用7z打包整个shp文件夹不要只发一个.shp。这个习惯能省掉大量“你发的文件我打不开”的沟通成本。3.2 坐标系是裁剪前最需要较真的一步拿到市级shp之后第一件事不是急着裁剪而是打开图层属性看坐标系。常见情况有两种shp本身是WGS84经纬度坐标EPSG:4326或者CGCS2000地理坐标系EPSG:4490也可能是某种投影坐标系比如高斯克吕格投影或UTM。为什么要较真因为DEM的坐标系和shp不一致时裁剪结果会出现明显位移甚至裁出来一片空白。成都市域大致在东经103°到105°、北纬30°到31.5°之间这个范围如果你看到的坐标是几十万、几百万的数值那多半是投影坐标系如果是一百出头的经纬度就是地理坐标系。用图层范围快速判断比干看元数据直观得多。如果shp没有.prj文件GIS会默认按WGS84显示但实际数据可能是其他坐标系。这时候需要用“定义投影”工具手动指定正确的坐标系而不是直接做重投影。定义错了后面全部白做。3.3 从属性表确认边界范围是否覆盖全成都打开shp属性表看看里面有没有行政名称字段比如NAME、PAC或市辖区县名称。很多渠道发布的“市级范围”只包含主城区不一定覆盖全部县市。如果你做的是全市域分析先确认边界范围是否完整如果只需要中心城区反而要小心里面多出来的远郊地块。一个实用的验证方法把shp加载到GIS里叠加一份在线影像或路网数据肉眼确认边界轮廓是否符合预期。虽然土但最有效。4. 用shp裁剪成都市DEM两条主流路线实操4.1 ArcGIS操作Extract by Mask是首选工具很多人一提到裁剪第一个想到的就是“Clip”工具。栅格裁剪确实可以用Clip但做地形分析时我更推荐用“按掩膜提取”Extract by Mask。差别在于Clip是严格按矩形范围切割虽然能选“使用输入要素裁剪几何”但处理方式偏“切图”Extract by Mask则是把shp边界外的像元直接设为NoData输出结果在边界处更贴合分析需求后续做坡度、填洼时不会让边界外的高程参与计算。具体步骤打开ArcToolbox路径Spatial Analyst Tools → Extraction → Extract by Mask。Input raster选择DEMInput raster or feature mask data选择成都市级shp。在Environment设置里把Snap Raster设为DEM确保输出像元与原DEM对齐避免出现半个像元的偏移。指定输出路径运行。运行结束后叠加shp边界查看重点看边缘处是否有白边或黑边。如果边界外是大片黑色说明NoData显示设置不对如果白色通常只是符号化默认背景不代表数据错误。4.2 QGIS操作免费开源流程同样干净没有ArcGIS许可也无所谓QGIS完全能完成这个活。操作路径菜单栏选择“栅格” → “提取” → “按掩膜图层裁剪栅格”。输入图层选DEM掩膜图层选成都shp。勾选“将结果裁剪到掩膜范围”之类的选项。输出格式选GeoTIFF点击运行。QGIS还提供“构建虚拟栅格”功能如果你手头是分幅DEM可以先用这个功能把多个tif拼成一个虚拟栅格再裁剪速度快也不占磁盘空间。4.3 裁剪前必查NoData设置是否正确这是最容易踩坑的一步。12.5米DEM往往来自雷达数据处理水体、陡坡和雷达信号遮挡区域可能出现无效值。有的产品用-32768表示无数据有的用0如果加载后不检查直接把0值当成真实海拔后续生成的坡度、等高线会出现一堆离谱结果。检查方法很简单在ArcGIS里查看图层属性中的统计数据看最小值和最大值。如果最小值是一个“整数”而周围地形显然不可能有负数或零值就要警惕。用Copy Raster工具重新输出设置NoData值或者用栅格计算器里的SetNull函数把异常值剔除。命令行下可以用gdal_translate快速处理gdal_translate -a_nodata -32768 -of GTiff input.tif output_nodata.tif这一步做完再去做Extract by Mask后面所有分析才有意义。我见过不少项目做完后才发现整个盆地全变成“海拔0米”原因就是NoData没设对。4.4 裁剪结果的质量检查裁剪完成后不只是看一眼“有没有裁出来”就结束。建议再做三层检查查看输出栅格的像元数量、范围是否和shp吻合。用Zonal Statistics as Table统计成都范围内的平均高程、最大最小高程和已知地理常识比对。叠加shp边界放大边缘检查是否有锯齿、错位或异常条带。如果一切都对这份“成都市12.5米DEM”才算真正可以用了。5. DEM到手后的进阶玩法等高线、坡度坡向、三维表达5.1 一键生成成都市的等高线有了裁剪后的DEM生成等高线是再自然不过的需求。ArcGIS里直接用“等高线”Contour工具设置等距间隔成都平原地区5米或10米能画出很多细节西边山区等高线加密会非常密集出图时建议用50米或100米间隔。QGIS里同样有“等高线”功能只是工具栏位置不同。等高线生成后记得把成果转成shp或DWG方便进CAD或出专题图。配合shp边界整张图就很像正式的地形图成果了。5.2 坡度、坡向、山体阴影地形分析的三大件地形分析最常用的三个衍生数据是坡度Slope、坡向Aspect和山体阴影Hillshade。坡度能直接看出哪些区域陡峭对地质灾害排查、建设用地适宜性评价很有帮助坡向影响日照和植被分布山体阴影主要用于地貌可视化让平坦的平原也能看出立体感。实操经验是出图时把山体阴影图放在DEM图层上方并设置30%到50%的透明度底图用高程渐变着色效果会非常漂亮。这个组合我用了很多年比单纯拉伸渲染或者单一山体阴影都耐看。5.3 有了边界shp还能做渔网、KML、3D Tiles拿到市级shp它的用处不只在裁剪DEM。很多人喜欢把市域切成规则渔网然后做“百公里网格平均海拔”这类统计图用的就是GIS里的“创建渔网”Create Fishnet工具配合Zonal Statistics。原理很简单shp提供边界范围DEM提供每个格网内的高程最后输出一张按海拔着色的格网图。如果想在网页端展示地形效果可以把shp转成KML导入Google Earth或其它在线地图平台。更进阶的玩法是“shp转3D Tiles”用CesiumLab或FME之类的工具把DEM和成都边界shp处理成带地形起伏的三维瓦片之后在浏览器里就能看到整个成都市域的立体地形。做汇报演示时这种可视化效果很有冲击力。5.4 顺便说说DSM和DEM的区别很多人在网上搜“dsm生成dem”因为雷达数据经常打包提供的是DSM也就是数字表面模型它记录的是地表物体树冠、建筑、电线塔的高度而不是裸地高程。如果直接用DSM做坡度、汇水分析山区和城市区会出现与实际不符的突变。真正需要做地形分析时优先确认数据源是不是DEM。如果只有DSM可以通过滤波、局部最小值统计等算法去除地表物体的高度还原裸地地形但这个过程比较折腾。市面上所谓“12.5米分辨率”的雷达产品能不能直接当地形数据用关键就看它是不是已经做过地形校正和去除非地面点处理。6. 想自己重新下载12.5米DEM这些渠道和坑先知道6.1 常见开源渠道ASF DAAC和OpenTopography如果这套打包数据满足不了你比如需要成都周边更大的范围或者想获取不同时期的DEM做对比可以自己下载。12.5米DEM最常提到的源头是ASF DAAC搜索时选ALOS PALSAR传感器再按成都的经纬度画范围线就能看到可下载的产品标称分辨率约12.5米。OpenTopography也提供ALOS World 3D等全球DEM数据但需要注册账号下载时按提示选择范围和分辨率。自己下载的缺点是经常要按卫星轨道分块不是整整齐齐的一整幅成都范围。下载下来之后你还得自己做拼接、裁剪、重投影、NoData处理这一套流程熟练的话一两个小时不熟练的话折腾半天都有可能。所以包含市级shp的打包数据有它的价值。6.2 分幅数据回来的拼接细节如果下载到的是多个tif拼接时优先使用构建虚拟栅格Build Virtual Raster的方式。虚拟栅格不产生实际合并文件只是把多个分幅当成一个整体来读取后续裁剪、统计都直接按整幅处理方便又省空间。最终需要导出成品时再用“导出”或gdal_translate生成一个新的GeoTIFF。拼接时一定要注意各幅数据的坐标系、像元大小是否一致。如果投影不一致先统一投影如果像元大小有细微差异以较高分辨率的一幅为基准重采样避免拼接结果“东一块西一块”对不齐。6.3 精度问题12.5米不等于处处准确说句实话12.5米DEM在植被茂密的高山区、雷达阴影陡坡区、水域边界都可能出现明显误差。雷达数据虽不受云层干扰但对地形朝向敏感陡峭峡谷里会有阴影和叠掩现象。做初步选线、区域评估没问题但涉及工程勘察、用地审批这类精细决策还是得以实测数据和更高精度的测绘成果为准。另外一些产品的高程基准是大地水准面和国内常用的1985国家高程基准可能存在几米到十几米的差异。跨数据拼接时尤其要注意别用两套垂直基准的数据直接算土方量。6.4 归档备份的土办法最后分享一个我自己的归档习惯。每次拿到一套DEM数据我会在数据目录里放三个东西原始压缩包、完整的SHA256哈希值文件、一个README.txt里面写明数据来源、坐标系、高程单位、裁剪范围、处理日期。sha256sum *.7z CHECKSUM.sha256这样做的好处是半年后你翻到这份数据不用猜“这个tif到底处理过没有”一个文本文件全部说清楚。大数据时代真正值钱的不是下载速度而是数据到了手里之后的管理方式。压缩包保留一份处理后的成果放另一份哈希值核对无误下次用的时候心里特别踏实。这批成都12.5米DEM和市级shp按上面的流程走一遍解压、校验、检查NoData、确认坐标系、裁剪、分析基本不会出大问题。如果你手头正好有这份数据不妨现在就打开压缩包先7z l看看清单再决定下一步怎么处理。本文还有配套的精品资源点击获取 成都12.5米DEM数据包实用指南：解压、裁剪与地形分析
简介这份四川省成都市12.5米DEM数据包面向GIS、测绘、城乡规划、环境与灾害研究等从业者提供比常用30米分辨率更精细的高程栅格能更清晰呈现山脊、山谷、坡度等地形细节适合开展地形分析、坡度坡向提取、城市三维可视化与灾害风险评估。压缩包共12个文件总大小约80.13MB核心为TIFF格式高程栅格并配套shp/shx行政边界、prj投影定义、tfw坐标定位、dbf与vat属性表及ovr金字塔等文件分别用于定义坐标系、精准配准地理位置、存储属性信息和加速大图缩放可在ArcGIS、QGIS中直接加载。目前已有419人学习下载适合需快速获取成都高精度地形底图的用户使用。这份数据能省去自行下载、裁剪、拼接多源数据的流程一次拿到完整市域DEM与行政边界并支持淹没模拟、通视分析、工程选址、生态评价等进阶应用显著缩短前期数据处理时间。 在GIS数据交换群里泡久了经常会看到一类资源地名开头、分辨率居中、后缀还带个.7z。“四川省成都市12.5米DEM数据含市级范围shp文件.7z”就是典型代表。乍一看只是一个压缩包但对做过成都方向地形分析的人来说这套东西能省下大半天的预处理时间。12.5米DEM意味着每个像元对应地面12.5米见方比常见的SRTM 30米、ASTER 30米细腻不少山区沟谷、城区微地形都能看出轮廓附带市级范围shp文件相当于把裁剪边界和行政范围都给你备好了。这篇文章就围绕这套数据讲清楚数据价值、解压校验、shp边界该怎么用、裁剪怎么做以及如果数据不够用自己上哪儿重新下载还有什么坑要避开。无论你是做规划的、搞测绘的还是想用地形数据做可视化的小白都能从里面找到能直接抄走的操作。1. 12.5米DEM到底强在哪先搞清这个数据包的含金量1.1 DEM不是图片是一张高程数值表很多人第一次拿到12.5米DEM双击打开发现是一张灰蒙蒙的图以为就是普通卫片。其实DEM的每个像元记录的不是颜色而是这一小块地面的海拔高度。你可以把它理解成一张“海拔查询表”只不过这张表是栅格结构的每个格子都有明确的经度、纬度和高程值。12.5米分辨率指的是一个像元对应地面12.5米乘12.5米的实际范围。以成都市域约1.4万平方公里的体量来算一张完整的12.5米DEM格网数量至少在八九千万量级解压后动不动就是几个GB的GeoTIFF文件。这个数据量放到以前普通电脑加载都会卡现在倒是轻松很多但压缩包体积依然可观所以用7z打包几乎是标配。1.2 同为全球DEM12.5米和30米差在哪市面常见的开源DEM主要就那几款SRTM 30米、ASTER 30米、ALOS 12.5米还有商业的TanDEM-X 12米/30米。很多人只下载过30米觉得够用了但真正做过山区项目的人普遍会觉得30米在成都这种地形起伏大的地方不够看。数据源分辨率覆盖范围特点常见用途SRTM 30米约30米全球北纬60°到南纬56°覆盖广、下载方便但局部细节平滑区域宏观地形、水系提取ASTER 30米约30米全球易受云和条带影响垂直精度波动较大教学演示、大范围概览ALOS 12.5米约12.5米全球大部分陆地区域细节更突出雷达数据受天气影响小地质灾害、精细坡面分析TanDEM-X 12米约12米全球精度高但商业授权费用不低专业测绘、工程前期选择12.5米数据的核心原因是它能看出30米数据看不出来的地形细节。成都西边是龙门山一带的高山峡谷东边有龙泉山中部是平原。在山区12.5米数据能把冲沟、山脊转折、崩滑体边界这些微地貌显示得更清楚在平坝区域河道阶梯、人工台地、城中村的地形起伏也能有所体现。对地质灾害解译、洪水淹没模拟、机场净空评估、管线选线这类工作来说这种细节很关键。1.3 带shp文件为什么是加分项单独一份DEM对你来说只是一张高程图但如果配套一个成都市级范围的shp边界数据就变得好操作得多。你不必再从全国行政区划数据里手工提取边界也不用担心裁剪时把周边地市带进来。更实用的是shp文件可以做掩膜分析比如统计成都市域内的平均海拔、坡度分级面积或者生成成都范围的等高线、山体阴影图。2. 打开.7z的正确姿势解压、校验与文件清单核对2.1 GIS数据打包为什么偏爱7zGeoTIFF格式本身带有大量重复的块结构和元数据压缩空间很大7z算法能压出比zip更小的体积。更重要的是shp文件不是单个文件它必须由.shp、.shx、.dbf、.prj等多个文件共同组成。把整个shp文件夹用7z压成一个包传输、分享时就不会出现“只发了.shp别人打不开”的尴尬情况。所以“为什么我发的shp对方打不开”这类问题多半是没连其他后缀文件一起发。解开这个压缩包之前先养成一个习惯先弄清楚里面装了什么。用7z l可以不解压直接列出文件清单。这样你还没解压就能知道里面是一个完整的成都范围的tif还是按区县分幅的多个tif便于提前规划后续操作。2.2 Windows和Linux下的解压操作Windows下最简单装一个7-Zip或者Bandizip右键选择“解压到当前文件夹”即可。一个小提醒解压路径尽量不要带中文也不要把压缩包放在网盘同步目录里直接解压某些GIS老版本软件对中文路径和云盘占位文件不友好加载时会莫名其妙报错。Linux服务器上解压很多人的第一反应是unzip但7z格式需要p7zip组件# Debian/Ubuntu系 sudo apt install p7zip-full # CentOS/RHEL系 sudo yum install p7zip安装完成后解压7z x 四川省成都市12.5米DEM数据含市级范围shp文件.7z如果文件名太长或者带中文直接Tab补全即可。7za是7z的独立命令行版本功能上对解压来说差不多日常用哪个都行。2.3 解压前先做哈希校验别偷懒网盘或FTP传文件偶尔会出现下载不完整的情况。一个好几个GB的数据包如果中间坏了一个字节解压时可能只报一句“CRC Failed”然后整包报废。可靠的做法是下载完后用SHA256校验一下压缩包的哈希值和发布方提供的值比对。Linux下sha256sum 四川省成都市12.5米DEM数据含市级范围shp文件.7zWindows PowerShell下Get-FileHash .\文件名.7z -Algorithm SHA256如果发布方没有提供哈希值也可以用7z t测试压缩包完整性7z t 文件名.7z这个命令会逐个文件测试解压数据是否正确。压包测试通过再解压、再投入项目能避免很多隐藏问题。我自己的习惯是所有GIS原始数据下载完先记录哈希值再归档备份方便以后核对文件有没有被误改。3. 市级shp边界文件拿到手先确认坐标系再谈裁剪3.1 shp的一整套文件缺一个都不行解压后你会看到一堆扩展名不同的文件。.shp存几何图形.shx是空间索引.dbf存属性字段.prj记录坐标系还有可能附带.sbn、.sbx等辅助文件。很多新手只把.shp拖进GIS结果属性表是空的或者坐标系显示未知就是这个原因。用7z打包数据通常不容易缺文件但解压后如果打算把shp单独发给别人建议继续用7z打包整个shp文件夹不要只发一个.shp。这个习惯能省掉大量“你发的文件我打不开”的沟通成本。3.2 坐标系是裁剪前最需要较真的一步拿到市级shp之后第一件事不是急着裁剪而是打开图层属性看坐标系。常见情况有两种shp本身是WGS84经纬度坐标EPSG:4326或者CGCS2000地理坐标系EPSG:4490也可能是某种投影坐标系比如高斯克吕格投影或UTM。为什么要较真因为DEM的坐标系和shp不一致时裁剪结果会出现明显位移甚至裁出来一片空白。成都市域大致在东经103°到105°、北纬30°到31.5°之间这个范围如果你看到的坐标是几十万、几百万的数值那多半是投影坐标系如果是一百出头的经纬度就是地理坐标系。用图层范围快速判断比干看元数据直观得多。如果shp没有.prj文件GIS会默认按WGS84显示但实际数据可能是其他坐标系。这时候需要用“定义投影”工具手动指定正确的坐标系而不是直接做重投影。定义错了后面全部白做。3.3 从属性表确认边界范围是否覆盖全成都打开shp属性表看看里面有没有行政名称字段比如NAME、PAC或市辖区县名称。很多渠道发布的“市级范围”只包含主城区不一定覆盖全部县市。如果你做的是全市域分析先确认边界范围是否完整如果只需要中心城区反而要小心里面多出来的远郊地块。一个实用的验证方法把shp加载到GIS里叠加一份在线影像或路网数据肉眼确认边界轮廓是否符合预期。虽然土但最有效。4. 用shp裁剪成都市DEM两条主流路线实操4.1 ArcGIS操作Extract by Mask是首选工具很多人一提到裁剪第一个想到的就是“Clip”工具。栅格裁剪确实可以用Clip但做地形分析时我更推荐用“按掩膜提取”Extract by Mask。差别在于Clip是严格按矩形范围切割虽然能选“使用输入要素裁剪几何”但处理方式偏“切图”Extract by Mask则是把shp边界外的像元直接设为NoData输出结果在边界处更贴合分析需求后续做坡度、填洼时不会让边界外的高程参与计算。具体步骤打开ArcToolbox路径Spatial Analyst Tools → Extraction → Extract by Mask。Input raster选择DEMInput raster or feature mask data选择成都市级shp。在Environment设置里把Snap Raster设为DEM确保输出像元与原DEM对齐避免出现半个像元的偏移。指定输出路径运行。运行结束后叠加shp边界查看重点看边缘处是否有白边或黑边。如果边界外是大片黑色说明NoData显示设置不对如果白色通常只是符号化默认背景不代表数据错误。4.2 QGIS操作免费开源流程同样干净没有ArcGIS许可也无所谓QGIS完全能完成这个活。操作路径菜单栏选择“栅格” → “提取” → “按掩膜图层裁剪栅格”。输入图层选DEM掩膜图层选成都shp。勾选“将结果裁剪到掩膜范围”之类的选项。输出格式选GeoTIFF点击运行。QGIS还提供“构建虚拟栅格”功能如果你手头是分幅DEM可以先用这个功能把多个tif拼成一个虚拟栅格再裁剪速度快也不占磁盘空间。4.3 裁剪前必查NoData设置是否正确这是最容易踩坑的一步。12.5米DEM往往来自雷达数据处理水体、陡坡和雷达信号遮挡区域可能出现无效值。有的产品用-32768表示无数据有的用0如果加载后不检查直接把0值当成真实海拔后续生成的坡度、等高线会出现一堆离谱结果。检查方法很简单在ArcGIS里查看图层属性中的统计数据看最小值和最大值。如果最小值是一个“整数”而周围地形显然不可能有负数或零值就要警惕。用Copy Raster工具重新输出设置NoData值或者用栅格计算器里的SetNull函数把异常值剔除。命令行下可以用gdal_translate快速处理gdal_translate -a_nodata -32768 -of GTiff input.tif output_nodata.tif这一步做完再去做Extract by Mask后面所有分析才有意义。我见过不少项目做完后才发现整个盆地全变成“海拔0米”原因就是NoData没设对。4.4 裁剪结果的质量检查裁剪完成后不只是看一眼“有没有裁出来”就结束。建议再做三层检查查看输出栅格的像元数量、范围是否和shp吻合。用Zonal Statistics as Table统计成都范围内的平均高程、最大最小高程和已知地理常识比对。叠加shp边界放大边缘检查是否有锯齿、错位或异常条带。如果一切都对这份“成都市12.5米DEM”才算真正可以用了。5. DEM到手后的进阶玩法等高线、坡度坡向、三维表达5.1 一键生成成都市的等高线有了裁剪后的DEM生成等高线是再自然不过的需求。ArcGIS里直接用“等高线”Contour工具设置等距间隔成都平原地区5米或10米能画出很多细节西边山区等高线加密会非常密集出图时建议用50米或100米间隔。QGIS里同样有“等高线”功能只是工具栏位置不同。等高线生成后记得把成果转成shp或DWG方便进CAD或出专题图。配合shp边界整张图就很像正式的地形图成果了。5.2 坡度、坡向、山体阴影地形分析的三大件地形分析最常用的三个衍生数据是坡度Slope、坡向Aspect和山体阴影Hillshade。坡度能直接看出哪些区域陡峭对地质灾害排查、建设用地适宜性评价很有帮助坡向影响日照和植被分布山体阴影主要用于地貌可视化让平坦的平原也能看出立体感。实操经验是出图时把山体阴影图放在DEM图层上方并设置30%到50%的透明度底图用高程渐变着色效果会非常漂亮。这个组合我用了很多年比单纯拉伸渲染或者单一山体阴影都耐看。5.3 有了边界shp还能做渔网、KML、3D Tiles拿到市级shp它的用处不只在裁剪DEM。很多人喜欢把市域切成规则渔网然后做“百公里网格平均海拔”这类统计图用的就是GIS里的“创建渔网”Create Fishnet工具配合Zonal Statistics。原理很简单shp提供边界范围DEM提供每个格网内的高程最后输出一张按海拔着色的格网图。如果想在网页端展示地形效果可以把shp转成KML导入Google Earth或其它在线地图平台。更进阶的玩法是“shp转3D Tiles”用CesiumLab或FME之类的工具把DEM和成都边界shp处理成带地形起伏的三维瓦片之后在浏览器里就能看到整个成都市域的立体地形。做汇报演示时这种可视化效果很有冲击力。5.4 顺便说说DSM和DEM的区别很多人在网上搜“dsm生成dem”因为雷达数据经常打包提供的是DSM也就是数字表面模型它记录的是地表物体树冠、建筑、电线塔的高度而不是裸地高程。如果直接用DSM做坡度、汇水分析山区和城市区会出现与实际不符的突变。真正需要做地形分析时优先确认数据源是不是DEM。如果只有DSM可以通过滤波、局部最小值统计等算法去除地表物体的高度还原裸地地形但这个过程比较折腾。市面上所谓“12.5米分辨率”的雷达产品能不能直接当地形数据用关键就看它是不是已经做过地形校正和去除非地面点处理。6. 想自己重新下载12.5米DEM这些渠道和坑先知道6.1 常见开源渠道ASF DAAC和OpenTopography如果这套打包数据满足不了你比如需要成都周边更大的范围或者想获取不同时期的DEM做对比可以自己下载。12.5米DEM最常提到的源头是ASF DAAC搜索时选ALOS PALSAR传感器再按成都的经纬度画范围线就能看到可下载的产品标称分辨率约12.5米。OpenTopography也提供ALOS World 3D等全球DEM数据但需要注册账号下载时按提示选择范围和分辨率。自己下载的缺点是经常要按卫星轨道分块不是整整齐齐的一整幅成都范围。下载下来之后你还得自己做拼接、裁剪、重投影、NoData处理这一套流程熟练的话一两个小时不熟练的话折腾半天都有可能。所以包含市级shp的打包数据有它的价值。6.2 分幅数据回来的拼接细节如果下载到的是多个tif拼接时优先使用构建虚拟栅格Build Virtual Raster的方式。虚拟栅格不产生实际合并文件只是把多个分幅当成一个整体来读取后续裁剪、统计都直接按整幅处理方便又省空间。最终需要导出成品时再用“导出”或gdal_translate生成一个新的GeoTIFF。拼接时一定要注意各幅数据的坐标系、像元大小是否一致。如果投影不一致先统一投影如果像元大小有细微差异以较高分辨率的一幅为基准重采样避免拼接结果“东一块西一块”对不齐。6.3 精度问题12.5米不等于处处准确说句实话12.5米DEM在植被茂密的高山区、雷达阴影陡坡区、水域边界都可能出现明显误差。雷达数据虽不受云层干扰但对地形朝向敏感陡峭峡谷里会有阴影和叠掩现象。做初步选线、区域评估没问题但涉及工程勘察、用地审批这类精细决策还是得以实测数据和更高精度的测绘成果为准。另外一些产品的高程基准是大地水准面和国内常用的1985国家高程基准可能存在几米到十几米的差异。跨数据拼接时尤其要注意别用两套垂直基准的数据直接算土方量。6.4 归档备份的土办法最后分享一个我自己的归档习惯。每次拿到一套DEM数据我会在数据目录里放三个东西原始压缩包、完整的SHA256哈希值文件、一个README.txt里面写明数据来源、坐标系、高程单位、裁剪范围、处理日期。sha256sum *.7z CHECKSUM.sha256这样做的好处是半年后你翻到这份数据不用猜“这个tif到底处理过没有”一个文本文件全部说清楚。大数据时代真正值钱的不是下载速度而是数据到了手里之后的管理方式。压缩包保留一份处理后的成果放另一份哈希值核对无误下次用的时候心里特别踏实。这批成都12.5米DEM和市级shp按上面的流程走一遍解压、校验、检查NoData、确认坐标系、裁剪、分析基本不会出大问题。如果你手头正好有这份数据不妨现在就打开压缩包先7z l看看清单再决定下一步怎么处理。本文还有配套的精品资源点击获取