ArcGIS 制作中国区的数字高程DEM地图（附中国区STRM 90m DEM百度云免费下载链接）[通俗易懂]

全栈程序员-用户IM • 2022年10月22日下午10:00 • 未分类

大家好，又见面了，我是你们的朋友全栈君。如果您正在找激活码,请点击查看最新教程,关注关注公众号 “全栈程序员社区” 获取激活教程,可能之前旧版本教程已经失效.最新Idea2022.1教程亲测有效,一键激活。

Jetbrains全家桶1年46，售后保障稳定

数字高程模型（Digital Elevation Model)，简称DEM
网上有很多中国DEM的下载链接，要么收费，要么地图不对（缺藏南、湾湾等等），要么版本太老
所以自己做了一个，流程如下：

第一步：NASA官网下载STRM 90m文件包（http://srtm.csi.cgiar.org）

在这里插入图片描述 (i) 第二张图的Tile size就是在选区域下载，所以也可以选5x5degree，但就要选很多区域，而且下载的时候要一个个点，太麻烦了，30x30degree的话，只需要选六块就能包住中国；

(ii) 下载速度可能有点儿慢，值得等待；

第二步：将六块地图拼接

在这里插入图片描述
(i) 一张一张读到ArcGIS中（上图之所以看着怪是因为每个子图的colorbar还没统一），可能会卡，别着急

(ii) 像素类型选16_BIT_SIGNED （因为原始图像的类型就是16位有符号整型，保持一致；可以通过属性查看原图像的像素类型）为arcgis中DEM图进行拼接时像素类型怎么填？

(iii)波段选1（因为不是真彩色，全部合到一幅上）；在ArcGIS下如何合并dem?

这时候咱们的DEM图还是栅格数据集，需要转换成tif格式才能进行掩模裁剪

第三步：转换格式栅格数据集→TIF

网上有很多DEM合并办法（镶嵌）和裁减办法，但是我发现用镶嵌办法得到的是栅格数据集而不是TIF文件，如果用栅格数据集去裁减，那么就必须保存在默认路径（Default.gdb）里，否则图像就怪怪的完全不对（丢失空间参考信息解决ArcGIS栅格数据空间参考丢失问题）但这样的话，我又没办法导出，就很鸡肋，所以就要转成TIF来裁减
在这里插入图片描述

第四步：根据边界线掩模裁剪出DEM

如果你的待裁剪图像文件没有空间参考信息丢失的问题，那就可以直接掩模裁剪，请看：
只有边界线的情况下如何在Arcgis上截取目标区域（掩膜裁剪）

如果你的待裁剪图像文件有空间参考信息丢失的问题，要注意：
掩模裁剪，因为空间信息的丢失，裁出来的图会怪怪的也不对，调试了很久，最后发现用代码反而能裁出来，代码很简单如下（也是自己写的，改四个路径就行，如果需要批量裁减，可以看用Python对ASCII文件进行批量掩模裁剪）：

# -*- coding: utf-8 -*-
#此代码是用以掩模裁剪，从全球DEM中裁剪出中国区

import arcpy

arcpy.CheckOutExtension("spatial")

arcpy.gp.overwriteOutput=1

arcpy.env.workspace = "F:\\DEM_SRTM\\China_include_tif"               #工作空间

raster = "F:\\DEM_SRTM\\China_include_tif\\china_include.tif"         #需要裁减的图像

mask= "F:\\DEM_SRTM\\China_include_tif\\China_boundary.shp"           #掩膜文件

out= "F:\\DEM_SRTM\\China_include_tif\\china_only.tif"                #输出文件

arcpy.gp.ExtractByMask_sa(raster, mask, out)


print("All done")