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

最新推荐文章于 2024-12-18 17:18:24 发布

空中旋转篮球

最新推荐文章于 2024-12-18 17:18:24 发布

阅读量3.2k

点赞数 2

分类专栏： Python遥感图像处理文章标签： DEM 地形指数

本文链接：https://blog.csdn.net/soderayer/article/details/124676651

版权

Python遥感图像处理专栏收录该内容

52 篇文章 ¥29.90 ¥99.00

订阅专栏

超级会员免费看

本文介绍了使用Python处理数字高程模型（DEM）数据来计算地形指数，特别是3*3窗口下的地形高差指数。通过定义计算函数并进行卷积操作，实现了地形指数的计算。虽然文中未详细展开其他指数如TPI、TRI和粗糙度，但提到了它们的计算。最终展示了计算结果的全局和局部效果图，显示平原地区的地形指数较小，山区则较大。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

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-

了解本专栏

超级会员免费看