Arcpy实现DEM数据镶嵌与裁切——以陕西省ASTER GDEM30数据为例

一、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)