一、多张dem影像拼接
(我好愚蠢哦,今天才知道arcgis怎么一次导入多个数据,以前都是一个一个点,急死人,我还一直觉得这点arcgis不如envi,原来是我自己的原因。。。)
1.打开要拼接的dem数据
2.拼接
步骤:ArcToolBox——>Data Management Tools——>Raster——>Raster Dataset——>Mosaic To New Raster
(注意:此处不要选择Mosaic,一定要选择镶嵌至新栅格)
(1)Input Rasters:自己要拼接的所有影像。可以直接点右侧文件夹全选数据。
(2)Output Location:输出位置。这个非常关键!!!需要提前建好一个文件夹,然后输出位置选择这个文件夹就行,不然你会无法进行下一步。
(3)Raster Dataset Name with Extension:输出的栅格数据集的名称,需要加扩展名。扩展名一般是**.tif**。
(4)Pixel Type:像素类型。这个需要与原始dem数据类型一致。查看方式如下:
右击任意一张dem图层,选择属性即可查看。
(5)Number of Bands:波段数。输入1即可。如果波段不一致会提示错误。
这是几个必填选项,其他的默认或不填即可。
拼接完成后的图像:
二:裁剪
拼接后往往需要裁剪以获得研究区的dem分布图像。
我的研究区是河南省,使用河南省矢量边界数据裁剪dem栅格数据。前提是二者投影一致。
**裁剪步骤:**ArcToolBox——>Spatial Analyst Tools——>Extraction——>Extract by Mask
(空间分析工具——>提取——>按掩膜提取)
三、DEM负值处理
这个其实我也很懵~
可以看到,我拼接后的图像dem范围是-64~2398的,但是拼接前单张影像最低值是-55,裁剪后得到的河南省区域内dem的范围也是-64-2398,上限没问题,下限这个负值有问题,可能是拼接过程中因为某种原因导致最后出现的异常值。但是这并不影响整个区域的高程分布,从后面可以看出。
我百度查看了一下,河南境内没有海拔低于0的值,所以这个负值必须得处理掉,处理方法:
ArcToolBox——>Spatial Analyst Tools——>Extraction——>Extract by Attributes
where clause:点击右侧SQL,输入:“Value”>0
这里Value就是指dem值。
结果如下:
可以看到最低值变为了5。
前面说到,负值不影响最后的高程分布是因为我调整色带后,发现高程分布符合河南省的地势地貌特征。所以负值应该就是拼接时某一像元的异常。(仅为个人猜测)
上图即为调整色带后的高程分布情况,根据我查看的多篇论文可以发现,这个结果完全正确。
但是那个负值的真正原因我不懂,希望有大佬可以告知~