做地理信息这么多年,“新疆地形数据下载”是我被问得最多的问题之一,尤其是“30米、12.5米、5米DEM”这三个分辨率到底去哪下、怎么下、下完怎么处理,很多人卡在第一步。新疆面积大、地形变化剧烈,从准噶尔盆地到塔里木盆地,从天山到昆仑山、阿尔泰山,一套靠谱的DEM数据往往是项目能否开工的前提。更关键的是,这三个分辨率对应着截然不同的使用场景:30米适合区域宏观分析,12.5米适合县级单元和中小流域,5米则直接关系到工程级应用。这篇文章我从数据源选型、账号注册、区域框选、批量下载到后期处理、常见坑位,把整条链路完整盘一遍,帮你省掉自己摸索的时间。
1. 先理思路:30米、12.5米、5米分别解决什么问题
1.1 三档分辨率到底怎么选
很多人一上来就问“哪个分辨率最好”,这是典型的选型误区。分辨率本质上是尺度的选择,不是越高越好,而是跟你的应用场景匹配。
30米DEM适合什么?省级、地州级的宏观分析,比如全疆范围的地貌分区、宏观坡度分级、流域脉络提取。这种尺度下你关心的是趋势,不是单条冲沟的细节。用30米数据做全疆的坡度图、坡向图,数据量可控,处理速度快,叠加其他30米分辨率遥感影像也正好匹配。常用的有SRTM 1弧秒全球数据、AW3D30、ASTER GDEM。
12.5米DEM适合什么?县级行政区、中小流域、新能源场址踏勘,比如一个县域的光伏选址、一条中小河流的汇水区划分。12.5米在新疆这种地形起伏大的区域,已经能看出山脊线、河谷阶地的轮廓,比30米细不少。目前最容易获取的是ALOS PALSAR数据,实际产品间距约12.5米,覆盖全球,在ASF DAAC平台可以免费下载。
5米DEM适合什么?工程前期、地质灾害评价、单体边坡分析、光伏阵列微观选址。5米分辨率可以支撑1:1万甚至更大比例尺的基本地形分析,能识别出小型滑坡体、采坑边界、冲沟摆动。但要注意,目前没有全球免费公开的5米DEM,新疆区域的5米数据主要来自商业卫星立体像对、无人机航测或机载LiDAR。后面我会讲两条实际可行的落地路径。
1.2 数据源选择的总体原则
我的经验是,选数据源先问五个问题:区域范围多大、需要的精度等级、是否有云量或时效性要求、坐标系是否好转换、版权或获取门槛是否可接受。
对新疆来说,还有一个现实因素:区域跨度大、条带多。比如下载全疆30米SRTM1数据,涉及多个经纬度分幅,不是单幅能搞定的。从数据源角度,优先选择标准分幅规则、命名清晰、能按坐标框选的平台。否则后期拼接时会因为文件名混乱而头大。
另外,我强烈建议下载时顺手把原始数据说明文档一并保存。很多项目做着做着要写报告,需要说明DEM来源、分辨率、采集时间、垂直精度,这些信息都藏在数据集的Metadata里。别等出报告时再回头找,那时平台版本可能已经更新,旧数据的元数据就找不回来了。
2. 主流DEM数据源盘点与选型对比
2.1 30米级:SRTM、AW3D30、ASTER GDEM怎么选
30米级数据源最常用的就是三个:SRTM、AW3D30、ASTER GDEM,它们在新疆区域的适用性有明显差异。
SRTM(Shuttle Radar Topography Mission)分为1弧秒(约30米)和3弧秒(约90米)两种。其中SRTM1覆盖北纬60°到南纬56°,新疆完全在内。它的优势是数据成熟、空洞相对少、元数据清晰,很多论文和工程报告默认用它。要注意:SRTM在极陡山区(比如天山、昆仑山局部)会有一定的空洞和条带噪声,但整体质量稳定。
AW3D30是日本提供的全球30米数字表面模型,由ALOS卫星立体像对生成。它的特点是保留了地表覆盖信息,比如建筑物和树冠高度,所以严格说起来更像DSM。在新疆的城市区域和绿洲农业区,AW3D30能看到更细腻的地表起伏,但裸岩山区反而可能因为影像匹配误差出现纹理噪声。适合做城市三维分析和植被覆盖区地形修正,不太适合做高精度山体切割分析。
ASTER GDEM同样是用卫星立体像对生成的,覆盖范围广,但在新疆局部地区水体和云区容易产生异常坑洞。我的建议是:优先SRTM1,如果某个区域空洞明显,再用AW3D30或者ASTER补洞。实际处理中,我通常把三者都下载同一区域的切片,用于互相校验和空洞填补。
2.2 12.5米级:ALOS PALSAR为什么能拿到
12.5米分辨率在全球公开DEM里属于“稀缺档位”,目前最容易获取的就是ALOS PALSAR系列卫星数据,由日本宇宙航空研究开发机构提供。在ASF DAAC平台上发布的ALOS DEM产品,横向间距约为12.5米,覆盖全球陆地,无需特殊授权,注册账号即可下载。
这个数据来自L波段SAR干涉测量,对地表具有一定穿透性,尤其在干旱区效果不错。新疆大部分区域属于干旱半干旱环境,地表植被稀疏,正好是SAR数据的擅长场景。实测下来,ALOS DEM在戈壁、荒漠、低植被区的表现比光学影像生成的DEM要稳定得多,很少出现大面积空洞。
需要提醒的是,ALOS DEM是DEM,不是DSM,它反映的是裸地表高程,对建筑、树木等人工地物进行了处理。这意味着在城区用12.5米数据看房子高度是不行的,但做地形分析、坡度提取、洪水淹没模拟非常合适。
2.3 5米级:公开还是自产
关于5米DEM,先把话说清楚:目前没有任何平台提供全球或全新疆免费的5米DEM。新疆区域能拿到的5米数据通常来自三类渠道:
第一类是商业卫星立体像对生成,比如0.5米到1米分辨率的立体影像,通过专业软件做密集匹配,生成5米格网DEM。优点是覆盖范围大,一次性可获取几千平方公里;缺点是成本高、处理周期长,而且对软件操作有门槛。
第二类是无人机航测或机载LiDAR,适合中小范围项目,比如一个水库库区、一段公路走廊带、一个矿山排土场。机载LiDAR穿透植被能力强,可以直接生成高精度DEM,5米格网只是最低要求,实际上能轻松做到1米甚至更细。
第三类是通过公开的DSM数据反算DEM。比如用高分辨率立体影像生成的DSM,配合滤波算法去掉人工地物和植被,再重采样到5米格网。这个方法适合预算有限但有一定处理能力的团队,后面第5章我详细讲。
3. 实操流程一:30米DEM下载全程(以SRTM1为例)
3.1 注册与定位新疆区域
在新疆做项目,我推荐用USGS EarthExplorer下载SRTM1,覆盖全疆、操作流畅、文件格式标准。第一次使用需要注册账号,流程不复杂,按要求填写邮箱和机构信息即可。注意密码有复杂度要求,建议用自己常用的强密码并保存好,后面下载12.5米ALOS数据时还要用类似的账号体系。
登录后,地图窗口就是你的选区工具。新疆东西跨度大,想下载全疆范围,最省事的方式是用坐标框选。可以在地图上手动拖动矩形框,也可以直接在左侧“Coordinates”标签下输入四角经纬度。比如要覆盖新疆中东部,可以输入北纬49°、南纬34°、东经96°、西经73°这样的范围。框选范围越大,命中的分幅就越多,下载总量越大。建议先做小范围试验,确认流程能跑通再放全疆。
3.2 数据筛选与下载
框选好区域后,点击“Data Sets”,在数据集列表里找到“SRTM”分类,勾选“SRTM 1-Arc Second Global”。这里有个小技巧:搜索结果里会出现多个SRTM版本,包括1弧秒和3弧秒。3弧秒是90米产品,不符合要求时别勾错。另外还有“SRTM Void Filled”版本,这是经过空洞填充的,优先选它,能省很多后期补洞的功夫。
选好数据集后点击“Results”,系统会列出所有落在你框选范围内的分幅文件。SRTM1按1度×1度标准分幅,每个文件对应一个小方格。这时你可以按文件名判断覆盖位置,也可以通过“Show footprint”在地图中查看覆盖范围。批量下载时,逐个勾选后点击“Download”选择GeoTIFF格式即可。
3.3 批量下载与文件说明
SRTM1分幅文件命名规则一般是“N40E075”这种格式,前两位表示起始纬度,后三位表示起始经度。下载后你会发现每个压缩包里有多个文件,包括一个GeoTIFF高程文件、一个数字高程模型说明文件和多个元数据文件。实际使用时只需要那个后缀为.tif的高程文件,其他文件可以保留备查但不必全部解压。
这里提醒一个下载顺序的问题:新疆全疆涉及上百个分幅,如果手动一个个点下载,很容易漏。建议在结果列表中先把所有分幅加入“Bulk Download”,然后生成下载清单,用浏览器自带的下载管理排队下载。下载过程中不要随意中断网络,SRTM1压缩包单个不大但数量多,中断后恢复比较麻烦,最好用支持断点续传的下载工具。
4. 实操流程二:12.5米ALOS DEM下载流程
4.1 ASF DAAC区域搜索
12.5米ALOS DEM从ASF DAAC平台下载,网址是search.asf.alaska.edu。这个平台由阿拉斯加卫星设施维护,界面和EarthExplorer风格不同但逻辑类似。首次使用同样需要注册,填写邮箱设置密码后就能登录。
进入搜索界面后,地图左侧有一个“Search”区域,可以通过导入GeoJSON、KML或者手动绘制多边形来限定范围。如果你想下载一个县域的12.5米数据,最简单的办法是导入该县界KML。没有KML的话,直接用矩形框选也行,但会多下载一些边界外数据。
4.2 筛选RTC产品并下载
在搜索条件里找到数据源或产品类型,选择ALOS,再选中ALOS DEM相关产品。需要说明的是,ASF平台上有ALOS PALSAR的多种产品,包括原始雷达影像、RTC产品、DEM产品。我们下载DEM时明确选择ALOS DEM,而不是原始雷达影像,否则下载下来的数据还是复数影像格式,普通GIS软件根本打不开。
设定好区域和产品类型后,点击搜索,平台会返回所有覆盖该区域的产品记录。每条记录的“Granule”字段是唯一标识,右侧的“Download”按钮可以直接下载单个文件。12.5米数据多以GeoTIFF格式提供,单景覆盖范围比SRTM1大一些,但受轨道姿态影响,新疆区域可能需要几景拼接才能覆盖完整研究区。
4.3 12.5米数据落地后的文件说明
下载下来的ALOS DEM文件命名通常包含轨道号、日期、版本号等一串字符。使用前最好在GIS软件里先加载一下,查看其空间参考和值域范围。ALOS DEM的坐标基准通常是WGS84或UTM投影,高程单位是米。如果发现高程值出现负值,先别慌,可能是局部区域低于海平面,比如吐鲁番盆地,这是正常的。
这里有一个我踩过的坑:ALOS DEM在新疆山区的细节比SRTM1多,但个别版本在陡峭峡谷处会产生异常的“仰角”或“锯齿”现象,看起来像长条形的错误凸起。遇到这种情况,加载山体阴影图逐景检查,如果异常区域面积不大,可以在后期处理中通过平滑或插值修复;如果面积大,就换一景日期不同的数据试试。
5. 实操流程三:5米DEM的两条落地路径
5.1 路径一:用高分辨率立体影像生成DSM再转为DEM
这条路径的核心思路是利用商业或免费高分辨率影像的立体像对,生成DSM数字表面模型,再做滤波处理得到DEM。常规流程是:导入立体像对、进行空三加密、生成点云、内插DSM、滤波去除非地面点、重采样为5米DEM。
具体操作中,空三加密是第一步,也是最依赖经验的一步。你需要量测控制点、设置连接点匹配参数。如果区域内有已知高程控制点,一定要优先利用,能显著提高绝对高程精度。没有控制点时,可以借助12.5米或30米DEM做基准约束,但这种方式的绝对精度会弱一些。
生成DSM后不要直接当DEM用。DSM里包含了树冠、房屋、输电塔等地物,直接用于工程分析会产生很大误差。滤波时可以采用渐进加密三角网、形态学开运算或布料模拟滤波算法。处理后需做质量检查:生成等高线与原始影像叠加,重点看山谷是否存在负异常、房屋边界是否残留凸包。
这种方法通常能稳定输出5米或更细DEM。在新疆的大范围戈壁区,地表变化平缓,重采样到5米几乎没有问题;但在山区,点云密度往往不够,输出5米格网会有一些小孔洞,需要插值填补。
5.2 路径二:机载LiDAR点云生成5米DEM
如果项目区不超过几百平方公里,机载LiDAR是最省心的5米DEM获取方式。激光雷达直接记录三维点云,经过航带平差、噪点剔除、地面点分类后,非地面点自动去掉,生产DEM就是水到渠成的事。
实际操作中,地面点分类是最关键的环节。新疆草原、荒漠区地面点分类相对容易,但遇到密集灌木或高差大的陡崖,算法容易把低矮植被点误分为地面点。建议把点云分类后的结果做成剖面图逐段检查。点云密度越高,5米格网内的点数量越多,DEM也就越平滑。对于一个5米格网,至少保证每个格网内有5个以上的地面点,否则生成的DEM会有明显格网效应。
生成DEM时,可以采用反距离加权、自然邻域或克吕金插值。在工程项目里,我喜欢用克吕金,它能根据半变异函数拟合空间相关性,在起伏地形上效果好一些。但要注意克吕金参数多、计算慢,如果项目区很大,先用普通克里金或反距离加权跑一版,再局部修正。
5.3 路径三:直接采购航测成果
实践中有不少项目是直接招标或采购当地测绘单位已有的5米DEM。新疆近年开展了大量基础测绘项目,部分地区已有现成的高精度DEM成果,甚至还有数字正射影像DOM产品。你在启动新项目前,可以先问当地自然资源或测绘资料馆,是否能提供已有成果,这往往比从原始数据做起快得多。
这种方式需要注意数据使用授权和数据时相问题。已有成果的拍摄季节、坐标系、精度指标都可能跟你的项目不完全一致,拿到手后要检查Metadata或质检报告,确认高程中误差、平面中误差是否满足要求。对绝大多数地形分析任务来说,5米DEM够用且富余。
6. 下载后的常用处理:格式、拼接、裁剪与空洞修复
6.1 DEM文件格式和坐标基准问题
实际工作中你会遇到GeoTIFF、IMG、HGT、DTED等各种DEM文件格式。SRTM1常见的HGT格式是一个没有压缩的裸二进制文件,ESRI的ArcGIS和QGIS都支持直接读取;ALOS DEM和高分成果则以GeoTIFF为主。GeoTIFF的优势是自带地理参考信息,拖进软件里就能正确叠合。
新疆使用的是CGCS2000坐标系,而很多标清DEM原始成果是WGS84经纬度坐标。两者在新疆的差值不大,但对1:1万以上大比例尺工程来说不能忽略。建议在项目开始前把DEM统一转换到CGCS2000 / 高斯-克吕格投影,比如UTM Zone 43N或带号对应的3度带投影。转换时使用ArcGIS的Project Raster或QGIS的Warp工具,注意选择合适的重采样方法,高程数据推荐双线性或三次卷积,不要用最邻近法,否则会出现阶梯状地形。
6.2 DEM拼接、裁剪与分割
下载得到的分幅DEM必须经过拼接才能使用。在ArcGIS里用Mosaic To New Raster可以直接合并多个分幅;在QGIS里用Raster菜单下的Merge工具也可以。合并时有一个细节:不同来源的DEM重叠区高程可能不一致,比如SRTM1和ALOS DEM重叠带存在系统偏差。此时优先使用单一数据源拼全疆,避免交叉混用。
裁剪操作很常用。比如你研究区是天山某山区,只需要把整个新疆DEM裁剪出子集。用矢量边界做裁剪时,注意坐标系必须一致,否则裁剪结果会偏移。QGIS的Clip Raster by Mask Layer、ArcGIS的Extract by Mask都很好用。这里我还想提一下“DEM分割”这个操作,它一般不是指裁剪,而是指按地形特征把DEM分成不同区域,比如按坡度分级或按流域划分。可以使用栅格重分类或水文分析模块实现,和单纯的矩形裁剪是两码事。
6.3 空洞、负值和异常凸起修复
任何公开DEM在特定的地形和地表条件下都会出现空洞。SRTM在高山区会有“洞”,ALOS在峡谷会有异常尖峰。对于空洞,最常规的做法是用邻域均值填充滤波。ArcGIS的Focal Statistics配合Con函数能自动用周围有效像元填补空洞。如果空洞面积较大,建议用低分辨率DEM作为背景补充,比如用ASTER GDEM补SRTM的空洞,再用平滑滤波处理过渡。
负值问题在吐鲁番盆地非常典型,把海平面以下的负高程当成NoData是新手常犯的错。在填写NoData值或设置色带范围时,要注意保留负值区域。异常凸起往往表现为局部像元高程值特别大,比如ALOS DEM在雪线附近的错误匹配。可以用低通滤波消除,但滤波强度不要过大,否则真实峡谷和陡坎会被抹平。
7. 高频问题排查与避坑清单
| 问题现象 | 可能原因 | 解决办法 |
|---|---|---|
| 下载时提示登录超时或页面无响应 | 国外平台网络波动或账号过期 | 清理浏览器缓存后重新登录;使用批量下载清单分时段下载 |
| 下载下来打开是黑白色且有拉伸 | 高程值域未正确显示 | 在图层属性里使用“Stretch”自动拉伸,检查NoData值设置 |
| DEM拼接后出现明显接缝 | 相邻分幅坐标基准或精度不一致 | 使用同一数据源、同一坐标系拼接;必要时用羽化过渡消除接缝 |
| 山区DEM有洞或负异常 | 原始数据覆盖不足或传感器匹配错误 | 用多源数据补洞,使用低通滤波修复异常像元 |
| 12.5米ALOS数据范围与研究区不匹配 | 框选区域坐标偏移 | 检查输入的经纬度是否反了,新疆经度在73-96°,纬度在34-49° |
| 5米DEM网格感明显 | 点云密度不足或插值参数不当 | 增加点云密度,改用更平滑的插值算法,比如径向基函数插值 |
关于“autosar dem”这个搜索词我也多说一句。很多人在搜DEM时会看到汽车软件领域的“AUTOSAR DEM”,那是汽车电子里的诊断事件管理接口,跟GIS地形数据完全是两码事。搜索时建议加上“地形”“高程”或“GIS”等限定词,避免拿到完全无关的资料。
另一个常见误区是混淆DOM、DSM和DEM。DOM是正射影像,主要提供真实地表颜色纹理;DSM是数字表面模型,反映包括树木、建筑在内的地表最高面;DEM是数字高程模型,只反映裸地形。三者区别可以用一句话记:DEM是“去掉毛刺的地面”,DSM是“连房顶树尖一起算”,DOM则是“能看到东西长什么样的像片”。在项目中,用DOM做背景底图、用DEM做地形分析、用DSM做地表覆盖高度估计,才是标准组合。
最后再分享一个我自己的习惯
每次给新疆项目下载DEM,我都会先建一个目录模板,按“原始数据、中间成果、最终成果、元数据”四层放文件。原始数据永远单独存一份不修改,所有拼接、裁剪、补洞操作都输出到中间成果目录。这样即使中间某一步做坏了,也能随时回退到原始数据重来。地形数据处理最忌讳的就是在原文件上反复覆盖,等你意识到“这步滤波好像过头了”的时候,原始数据已经找不回来了。养成版本管理的习惯,比会多少个GIS命令都重要。