GIS数据获取：全国34省DEM高程数据免费下载（ASTER GDEM V003，30m精度，百度网盘+夸克网盘）
摘要：本文提供基于 ASTER GDEM V003 的全国34省DEM高程数据完整下载指南。包括34省成果裁切数据（10.61GB，GeoTIFF格式）和1123幅原始分幅数据（49GB），两种获取方式——公众号关键词自动回复和YouGIS数据助手平台。附完整34省数据清单（编码/大小/图幅数/高程范围）、Python下载验证脚本、坐标系检查代码和常见问题解答。
关键词：DEM数据下载、ASTER GDEM、高程数据、GeoTIFF、30m DEM、免费GIS数据、数字高程模型、百度网盘
标签：
GIS遥感空间数据库数据处理地理信息系统
数据简介
获取方式一：公众号关键词（推荐）
关注公众号，在对话框发送对应关键词，系统自动回复网盘链接（百度网盘 + 夸克网盘）。
成果数据（已裁切）
发送 DEM-省份名
DEM-北京 → 110000.tif (19MB)
DEM-河南 → 410000.tif (156MB)
DEM-四川 → 510000.tif (658MB)
DEM-新疆 → 650000.tif (1833MB)
原始分幅数据
发送 DEM-N{纬度}E{经度}
DEM-N30E114 → ASTGTMV003_N30E114_dem.tif（武汉附近）
DEM-N39E116 → ASTGTMV003_N39E116_dem.tif（北京）
关键词不区分大小写，省份支持简称（如 DEM-内蒙、DEM-黑龙）。
获取方式二：YouGIS数据助手
支持按省份或编号检索，查看数据元信息（分辨率/坐标系/范围），适合批量下载多个省份。
成果数据完整清单
34省数据（按大小排序）
总计：10,868 MB (10.61 GB) | 124.8亿有效像素
Python下载验证脚本
批量验证下载的DEM文件
import rasterio
import os
from pathlib import Path
def verify_dem(tif_path, expected_crs='EPSG:4326'):
"""验证下载的DEM文件完整性和坐标系"""
with rasterio.open(str(tif_path)) as src:
arr = src.read(1)
valid = arr[arr != src.nodata]
# 检查坐标系
actual_crs = str(src.crs)
crs_ok = expected_crs in actual_crs
return {
'file': Path(tif_path).name,
'valid': True,
'crs': actual_crs,
'crs_ok': crs_ok,
'dtype': str(src.dtypes[0]),
'pixels': valid.size,
'elev_min': int(valid.min()),
'elev_max': int(valid.max()),
'size_mb': os.path.getsize(str(tif_path)) // (1024 * 1024)
}
# 批量验证
DOWNLOAD_DIR = Path('a_output')
for tif in sorted(DOWNLOAD_DIR.glob('*.tif')):
info = verify_dem(tif)
status = '✅' if info['crs_ok'] else '⚠️ 坐标系异常'
print(f"{status} {info['file']}: "
f"{info['pixels']:>12,} px | "
f"{info['elev_min']:>5}~{info['elev_max']:<5}m | "
f"{info['size_mb']:>4}MB | "
f"CRS: {info['crs']}")
GDAL命令行快速检查坐标系
# 检查GeoTIFF的坐标系和元信息
python -c "import rasterio; src=rasterio.open('410000.tif'); print(f'CRS: {src.crs}\nSize: {src.width}x{src.height}\nBBox: {src.bounds}')"
# 或使用gdalinfo（如已安装GDAL命令行工具）
# gdalinfo 410000.tif | head -20
读取并可视化
import rasterio
import numpy as np
import matplotlib.pyplot as plt
with rasterio.open('410000.tif') as src:
dem = src.read(1)
dem = np.where(dem == -9999, np.nan, dem)
plt.figure(figsize=(10, 8))
plt.imshow(dem, cmap='terrain')
plt.colorbar(label='高程 (m)')
plt.title('河南省 DEM高程图')
plt.savefig('henan_dem.png', dpi=150, bbox_inches='tight')
工程注意事项
1. 噪声像素处理
650000.tif（新疆）存在104个极端噪声像素（-32287~32008m），使用前需过滤：
arr[(arr < -1000) | (arr > 9000)] = -9999
2. 分辨率匹配
与其他栅格数据叠加分析时，需确保分辨率一致：
# GDAL重采样到不同分辨率（命令行）
# gdalwarp -tr 0.001 0.001 -r bilinear input.tif output_resampled.tif
# 0.001° ≈ 100m，-r bilinear 适合连续高程数据
重采样方法选择：
bilinear：连续表面（高程、温度），推荐
nearest：分类数据（土地利用），保持类别不混
cubic：更平滑但计算量大
3. 引用规范
学术使用时建议引用：
NASA/METI (2019). ASTER Global Digital Elevation Model Version 3.
Retrieved from https://asterweb.jpl.nasa.gov/gdem.asp
网盘说明
两种网盘链接均通过公众号关键词或YouGIS获取。
常见问题
Q1: 坐标系是什么？需要投影转换吗？
A: WGS84（EPSG:4326），经纬度坐标。如果需要计算面积/体积/坡度，需先投影到Albers等积投影或UTM投影。
Q2: ASTER GDEM和SRTM哪个好？
A: 各有优势。ASTER覆盖到N83°（SRTM只到N60°），对中国高纬度地区覆盖更好。SRTM在中低纬度精度略优。建议高纬度用ASTER，低纬度可交叉验证。
Q3: 和SRTM对比技术参数？
Q4: 商业使用是否合规？
A: ASTER GDEM由NASA和METI发布，无需特别许可即可使用。详见 https://asterweb.jpl.nasa.gov/gdem.asp
参考文献
NASA/METI. ASTER GDEM Version 3. (2019). https://asterweb.jpl.nasa.gov/gdem.asp
NASA EarthData Search. https://earthdata.nasa.gov
地理空间数据云. https://www.gscloud.cn
YouGIS数据助手. https://yougis.com.cn/res/home
觉得有用就点赞收藏，有问题评论区留言，看到就回 👊
评论区