Chapter03第三章 空间平滑和空间插值.pdfVIP

第三章 空间平滑和空间插值 35 本章介绍基于GIS 的空间分析中两个常用操作：空间平滑和空间插值。空间平滑和空间插 值关系密切，它们都可以用于显示空间分布态式及空间分布趋势，二者还共享某些算法（如核 密度估计法Find/Replace All ）。空间平滑和空间插值的方法有很多种，本章只介绍其中最常用 的几种。 空间平滑与移动平均在概念上类似（移动平均是求一个时间段内的均值），而空间平滑术 是一个空间窗口内计算平均值。第 3.1 节介绍空间平滑的概念和方法，第 3.2 节是案例分析 3A ，用空间平滑法研究中国南方/泰语地名(Find/Replace all)分布。空间插值是用某些点的已知 数值来估算其他点的未知数值。第3.3 节介绍了基于点的空间插值，第3.4 节为案例3B，演示 了一些常用的点插值法。案例3B 所用数据与3A 相同，是案例3A 工作的延伸。第3.5 节介绍 基于面的空间插值，用一套面域数值（一般面单元较小）来估算另一个面域的数值（范围较 大）。面插值可用于数据融合以及不同面域单元的数据整合。第 3.6 节为案例 3C ，介绍两种 简单的面插值法。第3.7 节为小结。 3.1 空间平滑 与移动平均法计算一个时间段的平均值 （例如：五日平均温度）相似，空间平滑是将某点 周围地区（定义为一个空间窗口）的平均值作为该点的平滑值，以此减少空间变异。空间平滑 适用面很广。其中一种应用是处理小样本问题，我们在第八章会详细讨论。对于那些人口较少 的地区，由于小样本事件中随机误差的影响，癌症或谋杀等稀有事件发生率的估算不够可靠。 对于某些地区，这样的事情发生一次就可导致一个高发生率，而对于另外许多地区，没有发生 这种事情的结果是零发生率。另外一种应用是将离散的点数据转化为连续的密度图，从而考察 点数据的空间分布模式，可参见下面的第 3.2 节。本节介绍两种空间平滑方法（移动搜索法及 核密度估计法），附录3 介绍经验贝叶斯估计。 1 3.1.1 移动搜索法 36 移动搜索法 （FCA ）是以某点为中心画一个圆或正方形作为滤波窗口，用窗口内的平均值 （或数值密度）作为该点的值。将窗口在研究区内移动，直到得到所有位置的平均值。平均值 的变动性较小，从而实现空间上的平滑效果。FCA 也可以用于可达性测量（见第五章第 5.2 节）等其他研究中。 图3.1 是由 72 个网格单元组成的研究区。以网格 53 为中心的圆定义了一个包含 33 个网 格的窗口（如果网格的中心在圆内，则它属于这个圆），从而这 33 个网格的平均值为网格 33 的空间平滑值。将圆心在研究区内不同网格中心之间移动，就得到所有网格的空间平滑值。例 如，以网格 56 为中心的等半径的圆包含的 33 个网格定义了网格 56 的滤波窗口。值得注意的 是，研究区边界附近的滤波窗口包含网格较少，从而平滑度较低。这种效果称为边缘效应。 窗口的大小很重要，需要仔细确定。较大的窗口得到较强的空间平滑效果，从而更好地反 映了区域全局分布态势，而不是局部差异；较小窗口则得到相反的结果。我们可以通过试验不 同大小的窗口来寻求一个合适的窗口。 案例3A 详细介绍了FCA 法在ArcGIS 中的应用。我们先计算所有点之间的距离 （例如欧 式距离），然后提取那些小于或等于阈值的距离。在计算距离时，我们可以选用阈值距离作为 搜索半径，从而直接得到阈值范围内的距离。这里我们先计算所有点之间的距离，然后通过属 37 性值(距离)大小提取不同窗口内的观察值, 比较灵活些。在ArcGIS 中，我们可以基于提取的距 离表通过汇总起始点计算属性值的均值来计算在起始总周边范围内观察点的平均值。因为距离 表只包含那些阈值范围内的距离，只有在窗口之内的观察值才参与了汇总操作，实现了搜索的 效果。这样，我们可以省掉逐个画圆并查询圆内点的反复计算过程。 3.1.2 核估计 核估计与 FCA 的方法类似。两种方法都要用一个滤波窗口来定义近邻对象。所不同的 是，在 FCA 法中，所有对象的权重相同，而在核估计法中，