GDAL使用DEM数据计算地形指数

零、        前言

本文是接上文《GDAL使用DEM数据计算坡度坡向》，还是一篇关于DEM计算地形指数的一篇文章。这里所要计算的地形指数主要包括以下三个指数：地形耐用指数（Terrain Ruggedness Index，TPI）、地形位置指数（Topographic Position Index，TRI）和粗糙度（roughness）。

上面三个地形指数都是在一个3×3的窗口中计算的，3×3的窗口定义如图1所示，下文中的公式中的符号与图1中的符号一致。下文中的算法及其含义具体可以参考Wilson等人2007年的著作《Marine Geodesy》。

图1 3×3的地形窗口及其编号

一、        地形指数简介与计算公式