Python遥感图像处理应用篇:DEM地形指数计算

Posted 2022-06-30 空中旋转篮球

1.DEM地形指数

DEM(数字高程模型)数据可以计算不同的地形指数，地形指数也有多种，可以表达不同的物理意义。本文主要计算指定窗口大小下，地形高差指数，比如3*3窗口下，最大值和最小值之差，计算方法也比较简单。其他地形指数还有很多，比如地形耐用指数（Terrain Ruggedness Index，TPI）、地形位置指数（Topographic Position Index，TRI）和粗糙度（roughness).GDAL使用DEM数据计算地形指数这篇文章做了详细介绍。这里不重复了。

2.Python实现过程和代码

定义计算函数，这里使用3*3窗格大小进行卷积计算，即计算得到每一个中心像元的地形指数值，这里filter直接完整写出来了，也可以使用双循环写出。如需计算其他地形指数，可以修改函数中的公式即可完成计算。

def GetMaxMinus(i,j,demarray):
    #计算3*3窗格最大值最小值之差
    filter33=[demarray[i-1][j-1],demarray[i-1][j],demarray[i-1][j+1],demarray[i][j-1],demarray[i][j],demarray[i][j+1],demarray[i+1][j-1],demarray[i+1][j],demarray[i+1][j+1]]
    filter33_max=max(filter33)
    filter33_min=min(fil