问题:在使用GlobeLand30数据的过程中,经常需要通过镶嵌的方法来获得完整的研究区范围的土地覆被分类数据,但是直接进行镶嵌后发现,无论选择ArcGIS中哪一种镶嵌方法,都存在黑边(即拼接接缝处存在缝隙)。很多人认为这是数据本身存在的问题,但是其实在绝大部分情况下并不是数据的问题,可以通过增加一些简单的步骤来消除这些裂缝。
产生原因
仔细观察拼接后土地利用数据可以发现,在“裂缝”位置上的栅格值为0或255;回顾原始土地利用数据可以发现在图像的边缘位置上(镶嵌时会发生重合的位置),存在的无效值往往同时包括0和255,这导致无论使用Maximum或是Minimum方法都无法使得镶嵌后的结果不含无效值。所以我们需要将这两个值统一至一个无效值,才能合理地完成镶嵌操作。
具体操作
-
首先在ArcGIS中加载需要镶嵌的数据集;
这里以n50_30_2000lc030.tif
、n51_30_2000lc030.tif
为例 -
打开【Spatial Analyst工具】->【地图代数】->【栅格计算器】
分别对两个图层进行计算将值为255的栅格转为0:
n50_30_2000lc030.tif
:
Con("n50_30_2000lc030.tif" == 255,0,"n50_30_2000lc030.tif")
n51_30_2000lc030.tif
:
Con("n51_30_2000lc030.tif" == 255,0,"n51_30_2000lc030.tif")
-
打开【数据管理工具】->【栅格】->【栅格数据集】->【镶嵌至新栅格】
将刚刚使用【栅格计算器】得到的两个新栅格添加进来,按需要设定输出位置、具有扩展名的栅格数据集名称、栅格数据的空间参考、像素类型、像元大小以及波段数(本例中为1);最重要的是将镶嵌运算符设置为MAXIMUM
。
这样就可以得到镶嵌好的,没有“裂缝”的GlobeLand30土地覆被数据了。