一、DEM数据获取
1.1 数据源介绍
ASTER GDEM,全称为先进星载热发射和反射辐射仪全球数字高程模型,与SRTM一样为数字高程DEM。该数据是根据 NASA的新一代对地观测卫星Terra的详尽观测结果制作完成的,数据覆盖范围为北纬83°到南纬83°之间的所有陆地区域,达到了地球陆地表面的99%。
投影方式:UTM/WGS84
覆盖范围:全球
空间分辨率:1弧度秒 (约30 米)
精度:垂直精度20米,水平精度30米
1.2 数据下载
打开地理空间数据云网站,点击进入官网,点击菜单栏数据资源选项,公开可下载的数据列表中有ASTER GDEM 30M 分辨率数字高程数据。
点击菜单栏高级检索,我们首先选择数据集:ASTER GDEM,然后采取按行政区范围方式选中陕西省,然后点击检索进行数据查找。如图所示,总共搜索到40幅DEM数据,依次点击右边的下载符号下载即可。
下载好的DEM数据为压缩包形式:
二、Arcmap处理
将下载好的DEM数据解压后加载到Arcmap中,如下图所示:
2.1 DEM镶嵌
通过【镶嵌至新栅格】 工具进行DEM拼接。工具路径:ArcToolbox——>数据管理工具——>栅格——>栅格数据集——>镶嵌至新栅格。
镶嵌工具参数设置:
DEM镶嵌结果:
2.2 DEM裁切
通过【按掩膜提取】 工具进行DEM裁切。工具路径:ArcToolbox——>Spatial Analyst Tools工具——>提取工具——>按掩膜提取。
利用陕西省行政区划面状矢量文件作为掩膜范围进行裁切。掩膜工具参数设置:
DEM裁剪结果:
三、Arcpy处理
3.1 函数介绍
【镶嵌至新栅格】 工具对应函数定义:
【按掩膜提取】 工具对应函数定义:
3.2 核心代码
#DEM批量镶嵌至新栅格
env.workspace = src_path
rasterlist = arcpy.ListRasters("*", "IMG")
data = []
for raster in rasterlist:
data.append(raster)
arcpy.MosaicToNewRaster_management(data, src_path, "DEM1.img","","16_BIT_SIGNED", "30", "1", "LAST","FIRST")
#按掩膜提取
input_raster = os.path.join(src_path, "DEM-315 .img")
clip_shp = os.path.join(src_path, "sx.shp")
out_raster = os.path.join(src_path, "sx.img")
#outExtractByMask = arcpy.sa.ExtractByMask(input_raster,clip_shp)
#outExtractByMask.save(out_raster)