尧图网站设计 尧图网站设计YAOTU DESIGN
ARTICLE DETAIL

资讯详情

深耕网站设计与一线实操的经验洞察。

全球DEM数据shp文件实测:从解压到ArcGIS/QGIS地形分析全流程

全球DEM数据shp文件实测:从解压到ArcGIS/QGIS地形分析全流程 简介一份覆盖全球的DEM数字高程模型数据包面向使用ArcGIS进行地形分析与制图的GIS学习者与从业者解决了快速获取全球范围基础高程数据的痛点。压缩包共4个文件核心为全球DEM.tif栅格影像配套全球DEM.tif.ovr金字塔文件以提升显示与分析速度另含两个xml元数据文件记录投影与空间参考信息整体仅18.82MB便于下载与存储。该资源已有1615人学习使用数据经预处理并已注明含投影可直接在ArcGIS中加载用于坡度坡向提取、流域划分、地形剖面等基础分析。对于需要开展全球尺度的教学演示、初步科研或可视化项目这套轻量级DEM数据提供了经济实用的解决方案借助配套文件可减少重复处理步骤快速产出分析结果。 做GIS的人应该都经历过这种时刻满世界找数据不是下载链接失效就是文件大得离谱解压要等半天。所以我看到“全球DEM数据(shp文件不大很好用).rar”这个资源的时候第一反应是“靠谱吗”。实测下来这份数据确实对得起“不大很好用”这几个字文件体积控制得不错解压后是标准的shp格式在ArcGIS、QGIS里都能直接拖进去用尤其适合做小范围分析、出图底图、快速查看地形起伏。这篇文章我就从实战角度把这份数据的解压、打开、转换、应用整个流程拆开讲清楚顺便把过程中最容易踩的坑一并列出来。1. 这份数据包里到底有什么1.1 为什么DEM会用shp格式很多刚接触GIS的朋友一听到DEM就默认是栅格数据比如tif、img、grid这些。这个理解没错标准DEM确实是一张灰度图每个像素值就是高程。但网上下载到的“DEM数据”不一定是纯栅格尤其是标注了shp文件的资源通常是下面几种情况之一等高线shp把DEM栅格提取成矢量等高线每根线带一个高程属性字段适合制图和快速查询。高程点shp在DEM上均匀或按特征点采样的点数据每个点有高程值适合做插值或地形分析。DEM范围/图幅边界shp这个更像是一个轮廓告诉你这个DEM覆盖了哪里本身没有高程信息。这份全球DEM数据包里我解压后主要看到的是等高线和范围边界两类shp同时附带一个说明文档。对于只想看地形走势、做规划分析、生成坡度坡向图的朋友来说等高线shp反而比原始栅格更直观而且文件体积小得多打开不卡。1.2 解压后的文件结构用解压工具打开rar后里面不是乱糟糟一堆文件而是按大洲或区域分好了子文件夹每个文件夹里有对应的shp文件组。这里要提醒一下shp格式不是单个文件而是一组文件的总称至少包含.shp几何信息也就是点线面的坐标.shx形状索引负责几何位置索引.dbf属性表存高程值、名称等字段.prj投影坐标系信息这个特别重要没有它软件里经常显示不了正确位置我拿到的这份数据每个区域都完整带上了这四个文件所以加载基本没问题。1.3 坐标系和精度要注意什么打开数据后第一件事不是急着出图而是看一眼右下角的坐标系信息。这份数据用的是WGS84地理坐标系单位是度也就是经纬度。适合做宏观分析比如看山脉走向、大区域地形起伏。但如果你要做面积计算、坡度分析这种依赖精确距离的活儿建议先投影到适合你所在区域的高斯投影或UTM投影否则算出来的坡度和长度会有偏差。精度方面这份数据相当于中等分辨率能看清大的山谷山脊但看不清小沟小坎。做小范围工程分析的时候最好还是去下载更高精度的DEM推荐几个免费渠道放在后面章节这里先不展开。2. rar文件解压和处理全流程2.1 下载前先看这几件事下载rar文件之前我习惯先确认三件事文件大小、来源是否可信、有没有附带说明文档。这份“全球DEM数据(shp文件不大很好用).rar”体积不大几十MB级别说明数据经过了合理的简化和压缩不是那种塞满原始栅格的巨型包。下载到本地后我建议先核对一下文件能不能正常打开再决定是否解压。有时下载断断续续文件损坏一解压就报错白等半天。2.2 解压工具怎么选rar格式的解压工具市面上主流的有这么几个WinRAR最老牌支持rar压缩和解压界面老派但稳定7-Zip免费开源解压速度飞快支持格式多360压缩、Bandizip国内用得多操作简单双击就能解压我平时主力用7-Zip因为它免费、没广告、解压rar也没问题。如果你遇到解压时提示需要密码先别急着找破解工具这通常有两种情况一是文件加密了发布者会给密码二是文件头损坏工具误报。优先回去看下载页有没有附带密码说明。2.3 解压报错排查三步走解压过程中最常见的报错是“文件损坏”或“不可预料的压缩文件末端”。我的排查顺序是重新下载一次尽量用浏览器自带下载别用那些多线程下载工具有时反而会把文件下坏。用WinRAR打开rar文件菜单里选“工具-修复压缩文件”让它尝试修复。如果还不行看看是不是下载的文件名被系统改了扩展名手动改成.rar再试。这套流程能解决九成以上的解压问题。3. 在GIS软件里加载shp文件3.1 ArcGIS/ArcMap中打开shp的正确姿势ArcMap和ArcGIS Pro加载shp文件很多人习惯用“添加数据”按钮去一层层找文件。其实有个更快的方法打开目录窗口Catalog直接定位到你的解压文件夹把shp文件拖到地图视图里就行。有一点要特别注意如果你的电脑装了多个GIS软件shp文件的默认打开方式可能会被占用。比如双击shp文件结果电脑用记事本打开了看到的是一堆乱码。这不是文件坏了是系统关联错了。正确做法是打开GIS软件用拖拽的方式加载。在ArcMap里加载后如果发现图形显示不全双击图层打开图层属性切到“符号系统”把等高线的分级设置调整一下就行。这份数据的属性表里有一个高程字段我通常用“分级色彩”来渲染一眼就能看清低海拔到高海拔的过渡。3.2 QGIS导出shp最简单方法QGIS相比ArcGIS更轻量加载shp也更直接。打开QGIS直接把shp文件拖进“图层”面板或者用“图层-添加图层-添加矢量图层”也行。如果你想从QGIS里导出一个新的shp最简单的操作是在图层上右键选“导出-保存要素为”格式选ESRI Shapefile设置好文件名和坐标系确定就完事。整个过程不超过10秒比ArcGIS里还要干脆。QGIS还有一个优势可以边预览边调整符号。比如把等高线按高程值用不同颜色显示在“图层属性-符号化”里选“已分类”字段选高程字段分类数量设成8-10档点“分类”就自动生成非常方便。3.3 把shp转成DEM栅格如果你想从等高线shp得到一份真正的DEM栅格也就是每个像元都有高程值那就需要插值了。分两种情况点是高程点用反距离权重法IDW或克里金插值在ArcGIS里是“3D分析-栅格插值-IDW”QGIS里是“插值-IDW”。线是等高线要先转成TIN再把TIN转栅格。ArcGIS里依次用“3D分析-创建TIN”和“TIN转栅格”。插值的结果取决于原始数据质量。这份数据做出来的DEM是平滑趋势型的细节不多但做区域地形分析、坡度坡向分级、剖面分析都够用。4. DEM数据的实际应用场景与操作4.1 用DEM提取坡度坡向拿到DEM栅格后最常用的两个操作就是提坡度和坡向。ArcGIS里路径是“3D分析-坡度/坡向”QGIS里是“栅格-分析-坡度/坡向”。提取完坡度之后很多人会按角度直接分级制图比如0-5度是平缓区5-15度是缓坡15-25度是陡坡25度以上是峻坡。实际项目中还要结合项目需求细分。如果你习惯用Python批量处理这里给一段常用代码import rasterio import numpy as np # 读取DEM with rasterio.open(dem.tif) as src: dem src.read(1) transform src.transform # 计算坡度简单差分法 x, y np.gradient(dem) slope np.arctan(np.sqrt(x**2 y**2)) * (180 / np.pi) # 计算坡向 aspect np.arctan2(-x, y) * (180 / np.pi) aspect np.where(aspect 0, 360 aspect, aspect) print(坡度范围, slope.min(), -, slope.max()) print(坡向范围, aspect.min(), -, aspect.max())这段代码适合快速预览精度要求高的场景建议还是用GIS软件内置算法毕竟处理了投影变形等细节。4.2 多幅DEM镶嵌有时候你下载的是分幅的DEM比如按经纬度1度一格切好的要用的时候得拼接起来。ArcGIS里的工具叫“镶嵌至新栅格”在“数据管理工具-栅格-栅格数据集”里。实际操作时有个细节把多幅DEM放进去后注意选好“像元类型”和“波段数”不然输出结果会出现黑边或值域不对。我一般设定像元类型为16位有符号整型波段数1这样高程值不会被截断。QGIS里对应的是“栅格-杂项-合并”。两种软件出来的结果都行选自己顺手的。4.3 加载到图新地球或三维地球软件shp格式的DEM数据还能直接加载到三维地球软件里比如图新地球。这个操作对做项目汇报、看实景关系特别有用。图新地球里点击“添加-添加矢量图层”选择你的shp文件再配合在线影像底图就能把地形等高线叠加到三维地形上。如果数据带高程属性还可以在三维场景里做拉伸显示效果相当直观。需要注意的是图新地球对shp的坐标系比较敏感如果你的数据不是WGS84加载后可能跑到海里。解决办法是在GIS软件里先把坐标系转换成WGS84再导出成新的shp。5. 常见问题与避坑指南5.1 属性表打开是空的这种情况通常不是文件问题而是shp的dbf文件编码不兼容。ArcGIS打开中文属性表乱码QGIS里也可能显示问号。最简单的办法用记事本打开扩展名为.cpg的文件如果没有就新建一个文件名和shp一致内容改成UTF-8保存后重新加载。5.2 shp坐标偏移了坐标偏移最常见的原因就是.prj缺失。你用GIS软件打开shp软件默认用了软件自带的坐标系自然就偏了。解决方法是找到原数据来源确认坐标系然后在图层属性里手动指定。如果实在找不到坐标系有一个笨办法把全球矢量边界也加载进来手动对位置然后用“空间校正”工具纠正。5.3 数据量“不大”就一定不够用吗未必。这份数据“不大”是因为它保存的是矢量化后的关键地形要素文件量小、加载快、操作流畅。对于大范围制图、前期踏勘、课堂演示它比几十GB的原始DEM更方便。但如果你需要做水文分析、淹没分析、汇水区提取这类对地形细节极其敏感的工作建议换用高精度DEM免费渠道有NASA SRTM最高30米分辨率全球覆盖ALOS AW3D30也是30米地理空间数据云国内访问快注册后就能下载USGS EarthExplorer找数据全但需要一点学习成本5.4 rar带密码怎么办压缩包带密码大概率是发布者为了保护数据或防止二次传播。先回去翻下载页面或说明文档看有没有附带密码。实在找不到用软件暴力破解不现实耗时且容易触发安全软件报毒。我的建议是直接放弃这包数据换其他渠道下载。市面上的免费DEM源很多没必要跟一个密码较劲。6. 实操经验杂谈最后分享几个我在处理这类数据时总结的小习惯。第一拿到任何shp数据第一步永远是打开属性表看字段而不是急着出图。属性表能告诉你这份数据真正包含什么、精度如何、字段命名是否规范。第二做任何分析前先复制一份原始数据扔到备份文件夹所有操作在副本上进行。这个习惯救过我很多次尤其是插值和镶嵌这类会写入新文件的工具。第三处理完的数据建议统一转成GeoPackage或FileGDB格式比散装shp文件组好管理太多一个文件打包搞定发给同事也不怕漏文件。数据这东西很多时候不是越贵越好越全越好适合自己的场景才是最好的。像这份“不大很好用”的全球DEM数据我拿到手后几十秒就能出图、就能跑分析省下来的时间够你多试好几套方案了。希望这篇内容能帮你少走点弯路把时间花在真正重要的分析上。本文还有配套的精品资源点击获取
返回列表