ArcGIS中栅格空值计算与替换方法

在研究或工作中，经常会遇到栅格数据中某些像元数据缺失或者不能被使用，也会因为异常值而将部分像元定义为空值(NoData)。再此情况下，我们需要重新计算空值像元中正常值，从某种意义上也是对空值像元的内向插值。

整体思路是根据空值像元周围值来计算空值像元，类似于线性内插。具体方法如下：
1、首先确定空值栅格数据的空值范围，根据空值范围、类型、数据结构选择不同邻域统计工具。需要注意栅格数据的背景值与NoData的区分，最好去掉背景值。

2、在大多数情况下可以选择Block Statistics或者Focal Statistics或者Filter等，邻域类型可以选择矩形或者圆形、环形，需要根据自己数据以及研究目的来确定。邻域半径大小越大，空值越平滑，但统计结果栅格的分辨率越小，同原始栅格范围出入较大。

需要注意的是，由于分辨率调整使得栅格数据行列数变化，这会同原始栅格范围有出入，此种情况可以选择通过重采样之后，求取此栅格同原始栅格的交集，保证范围大小一致。此过程会对空值处像元值产生变化，但可忽略不计。
3、根据邻域统计之后的栅格数据，替换原始栅格空值像元。