PostGIS实现DEM分析之地形崎岖度指数计算【ST_TRI】
一、函数概述
ST_TRI（Terrain Ruggedness Index），用于计算地形崎岖度指数。TRI 通过量化中心像素与其 8 邻域像素的高程差异，反映地形表面的粗糙程度和破碎度。该指标在生态栖息地评估、军事行动规划和灾害风险分析中具有重要应用。
二、参数解析
raster ST_TRI(raster rast,-- 输入DEM栅格integernband,-- 指定栅格波段（通常为高程波段）raster customextent,-- 自定义输出范围（可选，默认与输入一致）textpixeltype='32BF',-- 输出像素类型（默认32位浮点型）booleaninterpolate_nodata=FALSE-- 是否插值处理NoData);参数说明：
1. 计算逻辑：
- 对每个像素 P，计算其与 8 邻域像素的高程差绝对值之和，再取平均：
TRI(P) = [Σ| P.value - 邻域像素值 |] / 8- 结果值越大，地形越崎岖（如山区、破碎地貌）；值越小，地形越平坦（如平原、冰盖）。
2. 窗口限制：
- 函数固定使用3×3 邻域窗口（半径 1 像素），不可调整。
3. NoData 处理：
- 若窗口内存在 NoData，默认输出 NoData；启用interpolate_nodata可插值填充邻域值。
三、算法原理
1. 高程差计算：
- 对中心像素 P 及其 8 邻域像素 N₁,N₂,…,N₈：
差异值 = |P - N₁| + |P - N₂| + ... + |P - N₈|2. 指数计算：
TRI = 差异值 / 8- 本质为局部高程变异性度量，突出地形的微观起伏特征。
四、典型示例
示例 1：基础 TRI 计算
WITHdem_dataAS(SELECTST_SetValues(ST_MakeEmptyRaster(5,5,0,0,1,-1),1,1,1,ARRAY[[10,10,10,10,10],-- 平坦区域[10,15,15,15,10],-- 缓坡[10,15,20,15,10],-- 山顶（中心像素20）[10,15,15,15,10],-- 缓坡[10,10,10,10,10]-- 平坦区域]::doubleprecision[][])ASrast)SELECTST_TRI(rast,1,NULL,'32BF')AStri_rasterFROMdem_data;输出解析：
- 中心像素（20）的 TRI 值为：
[|20-15|+|20-15|+|20-15|+|20-15|+|20-10|+|20-10|+|20-10|+|20-10|]/8 = 5- 平坦区域 TRI 值为 0，缓坡区域 TRI 约为 2.5。
示例 2：插值处理 NoData 区域
SELECTST_TRI(rast,nband :=1,customextent :=NULL,pixeltype :='32BF',interpolate_nodata :=TRUE-- 插值填充邻域NoData)AStri_with_interpolationFROMdem_with_nodata;场景：处理含湖泊、建筑物等空洞的 DEM，避免空洞区域影响 TRI 计算。
示例 3：崎岖度分级统计
WITHtri_resultsAS(SELECTST_TRI(dem.rast,1)AStriFROMdem_table)-- 分级统计各地形崎岖度占比SELECTCASEWHENtri<1THEN'极平坦'WHENtriBETWEEN1AND3THEN'平缓'WHENtriBETWEEN3AND5THEN'中等崎岖'ELSE'高度崎岖'ENDASruggedness_class,COUNT(*)ASpixel_countFROMST_DumpValues(tri_results.tri)GROUPBYruggedness_class;五、输出值与地形类型关系
| TRI 值范围 | 地形特征描述 | 典型区域示例 |
|---|---|---|
| 0-1 | 极平坦 | 城市广场、湖面 |
| 1-3 | 平缓起伏 | 农田、低丘 |
| 3-5 | 中等崎岖 | 丘陵、梯田 |
| 5-10 | 高度崎岖 | 山区、峡谷 |
| >10 | 极其崎岖（破碎地形） | 喀斯特地貌、火山口 |
六、性能优化
1. 数据类型选择：
- 输出推荐使用32BF（浮点型）保留精度，若仅需分级可转为8BUI（整数型）：
ST_TRI(rast,pixeltype :='8BUI')2. 并行计算：
- 结合 PostgreSQL 并行查询加速大规模栅格处理：
SETmax_parallel_workers_per_gather=4;SELECTST_TRI(rast,1)FROMlarge_dem;3. 边界扩展：
- 对边缘像素，可先通过ST_Expand扩展栅格，避免边界效应：
WITHexpanded_demAS(SELECTST_Expand(rast,1)ASexpanded_rast-- 扩展1个像素FROMdem_table)SELECTST_TRI(expanded_rast,1)FROMexpanded_dem;4. 与其他指标结合：
- TRI 常与坡度（ST_Slope）、粗糙度（ST_Roughness）联合分析，全面描述地形特征。
七、应用场景
1. 生态栖息地评估：
- 野生动物（如山地山羊）偏好特定崎岖度的地形，TRI 用于栖息地适宜性建模。
2. 军事与机器人路径规划：
- 高 TRI 区域通行难度大，需避开或规划特殊通行策略。
3. 自然灾害风险分析：
- 高 TRI 区域易发生滑坡、泥石流，结合坡度和土壤类型可构建风险模型。
4. 户外休闲规划：
- 徒步路线设计中，TRI 值可量化路线难度，辅助路线规划。
八、总结
通过ST_TRI，可快速量化地形的局部变异性，为地形分析、生态研究和工程应用提供关键参数。合理处理边界和 NoData 区域，结合其他地形指标，能更全面地刻画地表形态。