在我们从事的项目工作中,通常会根据研究需要对研究区域的栅格数据(如遥感影像)进行不规则裁剪。在ArcGIS中利用不规则矢量多边形对栅格数据进行裁切主要有两种方法:
空间分析中的掩膜提取和栅格处理中的裁剪。
1空间分析中的掩膜提取
在ArcGIS中,空间分析中的掩膜提取是通过ArcToolbox->Spatial Analyst Tools->Extraction->Extract byMask实现的。在Extract by Mask对话框中设置好输入输出参数以后,直接点击OK即可得到掩膜后的栅格数据。然而看似简单的过程实则包含了几个比较棘手的问题。
1.1错位重叠
裁切前后两个栅格的像元并不一定是完全重叠的,有可能错位,也就是说对应像元的中心点不是一个点,如下图1所示。这时裁切后像元的值可能是裁切前相邻像元的值,所以对应值就可能不一致,尤其在高程变化比较大的地方表现得更为明显。
图1 栅格数据裁切像元不完全重叠
若想让裁切前后的每个像元都完全重叠,可在裁切界面的Environment Settings里设置snap raster成裁切前的数据即可。
1.2 Pixel Depth变化
用不规则矢量多边形裁切栅格数据时,Extract by Mask对话框中没有提供设置Nodata的选项,这导致有时会出现裁剪后的pixel