Matlab: 众数插值

Introduction

土地利用覆盖分类矩阵也就是土地利用图，年代之间可能对应的空间范围有点点差异，尤其是边界地带可能不一致，这种不一致是指个别年份在这些地带无数据，像元数值是背景值，这给下面的多年分析就带来不便。如Fig. 1，左边是土地利用类型图，有限区域像元数值表示一定的土地类型，研究区之外表示无效区域以数值-99表达，右边是有效研究区范围的无数据地带，在左边以0数值表达。

Fig. 1

这里尝试以这类无数据像元毗邻的众数填充，规定毗邻的范围是八邻域（Fig. 2），八邻域之中频率高者即是该无数据像元的填充数值。个别极端情况下，某无数据像元毗邻八邻域均是无效区域或无数据像元，此时的邻域须向外围扩展一圈（如Fig. 2外围白色的像元），或再扩展一圈，直至可以得到有效的邻域数值为止。

Fig. 2

经此处理过程，土地利用图件的有效区域范围就不存在无数据像元了。

Example.rar.