当前位置：首页 > news >正文

告别ArcGIS！用GEE+QGIS搞定流域DEM下载与地形分析（附完整代码）

news 2026/6/6 10:39:24

告别ArcGIS！用GEE+QGIS搞定流域DEM下载与地形分析（附完整代码）

在GIS领域，数字高程模型（DEM）是地形分析的基础数据。传统上，ArcGIS凭借其完善的功能和稳定的性能，成为DEM处理的首选工具。但对于预算有限的个人用户或小型团队来说，ArcGIS高昂的授权费用往往成为一道难以逾越的门槛。今天，我将分享一套完全免费的替代方案——Google Earth Engine（GEE）与QGIS的组合工作流，让你无需ArcGIS也能高效完成从DEM获取到地形分析的全流程操作。

这套方案特别适合以下场景：

你需要快速获取特定流域或研究区域的DEM数据
你的预算有限，无法承担商业GIS软件的费用
你希望利用云端计算能力处理大规模地理数据
你需要进行基础的坡度、坡向等地形分析

1. 准备工作与环境配置

1.1 注册与配置GEE账号

Google Earth Engine是一个强大的地理空间分析平台，提供海量的遥感数据和处理能力。要使用GEE，首先需要：

访问Google Earth Engine官网并注册开发者账号
等待Google审核通过（通常需要1-2个工作日）
登录后进入代码编辑器界面（https://code.earthengine.google.com/）

1.2 安装QGIS

QGIS是一款开源的地理信息系统软件，功能强大且完全免费：

前往QGIS官网下载最新稳定版
按照向导完成安装（支持Windows、macOS和Linux）
首次启动后，建议安装以下插件：
- SCP（Semi-Automatic Classification Plugin）
- QuickMapServices（用于加载底图）

1.3 准备研究区域边界数据

你需要准备研究区域的矢量边界文件（如Shapefile），这是后续DEM裁剪的基础。确保你的边界文件：

采用WGS84坐标系（EPSG:4326）
文件大小不超过GEE的上传限制（通常10MB以内）
包含明确的几何信息

2. 在GEE中获取DEM数据

2.1 上传边界数据到GEE

在GEE代码编辑器中，点击"Assets"选项卡，然后选择"New"→"Shape files"上传你的边界文件。上传完成后，记下生成的asset路径（格式为projects/your-project-name/assets/your-file-name）。

2.2 编写DEM提取代码

以下是一个完整的DEM提取脚本，可以直接在GEE代码编辑器中使用：

// 1. 加载研究区域边界 var roi = ee.FeatureCollection('projects/YOUR_PROJECT/assets/YOUR_ASSET_NAME'); // 2. 获取SRTM DEM数据（全球30米分辨率） var dem = ee.Image('USGS/SRTMGL1_003').select('elevation'); // 3. 裁剪DEM到研究区域 var clippedDem = dem.clip(roi); // 4. 可视化设置 Map.centerObject(roi, 8); Map.addLayer(roi, {color: 'red'}, 'Study Area'); Map.addLayer(clippedDem, {min: 0, max: 5000, palette: ['blue', 'green', 'yellow', 'red']}, 'DEM'); // 5. 导出DEM到Google Drive Export.image.toDrive({ image: clippedDem, description: 'DEM_Export', scale: 30, region: roi.geometry(), maxPixels: 1e13, folder: 'GEE_Exports', crs: 'EPSG:4326' });

2.3 执行代码与数据导出

将上述代码粘贴到GEE代码编辑器中
替换YOUR_PROJECT和YOUR_ASSET_NAME为你的实际项目名称和资产名称
点击"Run"执行代码
在"Tasks"选项卡中点击"Run"开始导出任务
导出完成后，DEM文件将保存到你的Google Drive指定文件夹中

3. 在QGIS中进行地形分析

3.1 导入DEM数据

打开QGIS，创建一个新项目
通过"Layer"→"Add Layer"→"Add Raster Layer"导入从GEE下载的DEM文件
右键点击图层，选择"Properties"→"Symbology"调整显示样式

3.2 坡度分析

坡度是地表倾斜程度的度量，在QGIS中可以轻松计算：

打开"Processing Toolbox"（如果未显示，可通过"View"→"Panels"启用）
搜索并选择"Slope"工具
在参数设置中：
- 输入DEM图层
- 选择输出文件位置
- 设置输出格式（建议使用GeoTIFF）
点击"Run"执行计算

3.3 坡向分析

坡向表示坡面的朝向，对太阳辐射和植被分布等研究很重要：

# 在QGIS Python控制台中也可以使用以下代码进行坡向分析 processing.run("gdal:aspect", {'INPUT':'/path/to/your/dem.tif', 'ZEVENBERGEN':False, 'TRIG_ANGLE':False, 'ZERO_FLAT':False, 'COMPUTE_EDGES':False, 'OUTPUT':'/path/to/output/aspect.tif'})

3.4 地形阴影与可视化增强

为了更直观地展示地形特征，可以创建地形阴影：

使用"Raster"→"Terrain Analysis"→"Hillshade"工具
设置适当的方位角（Azimuth）和高度角（Altitude）
将生成的hillshade图层设置为DEM图层的底图，并调整透明度

4. 进阶技巧与性能优化

4.1 处理大区域DEM数据

当研究区域较大时，可能会遇到以下问题：

问题	解决方案
GEE导出限制	将区域分块处理，然后合并
内存不足	降低处理分辨率或使用QGIS的"Raster Calculator"分步计算
处理速度慢	启用QGIS的多线程处理（设置→Options→Processing）