6.2 地形模型TerrainTerrain.cs 共 386 行负责高程查询、视线LOS判定与地形坡度计算。全球最大地形高程常量 GlobalMaxTerrainElevation 8 850 m珠峰高度。高程查询接口为 GetElevation(lat, lon, bool, Scenario)其查询顺序为先检查是否落在自定义环境区内CustomEnvironmentZone.HasCustomTerrainHeight允许想定作者人为设定地形高度否则委托给 ITerrainProvider 接口由 BinaryGridProvider 从 BGD 格式的二进制高程栅格中读取。算法视线判定LOS 采样在观察者与目标之间沿大圆路径以 1/120°即 0.5 角分约 926 m为步长布设采样格逐点查询地形高程并与视线高度比较任一采样点的地形高于视线即判定为遮蔽。并行 vs 串行性能对比LOS 采样对比维度Parallel.ForEach并行串行 for 循环采样点数相同沿大圆路径按 1/120° 步长布设相同沿大圆路径按 1/120° 步长布设耗时10 万采样点约 15–30 ms约 60–120 msCPU 利用率多核满载接近 100% × 核心数单核运行约 1 核内存开销需额外分配结果数组 线程调度开销无额外分配开销最小适用场景采样点数量大 1 万、多核机器、实时性要求高采样点少 1 千、单核环境、对延迟不敏感选型建议选择并行当采样点数量大如跨洲际大圆路径、高分辨率地形栅格且运行环境为多核 CPU 时Parallel.ForEach能显著缩短 LOS 判定耗时适合实时雷达探测、多目标批量判定等对响应时间敏感的场景。选择串行当采样点数量少如短距离、低分辨率、或运行环境为单核/嵌入式平台时并行化的线程调度开销反而可能超过收益此时串行实现更简单、更省内存也更容易调试与复现结果。由于采样点数量大实现上使用 Parallel.ForEach 并行化再由 Scenario.DetermineLOS 汇总判定结果。算法最大地形坡度InvokeWatcher以目标点为中心向 8 个罗盘方向各取距离 0.121 424 72 nm约 225 m处的高程样本计算 Δh/450 得到坡度分母 450 m 为两侧采样点间距取 8 个方向中的最大值。该坡度值作为参数 MaxTerrainSlope 传入雷达探测函数用于修正地杂波强度——地形起伏越大杂波越强。算法区域陆地占比SearchWatcher给定多边形区域计算其几何中心并统计陆地覆盖比例用于判定任务区属于海域、陆域还是海陆交界。相关数据结构包括 TerrainGridCell地形网格单元、TerrainRasterCell栅格单元与 MaxMinRecord区域高程极值记录。7. 实战代码示例下面以 C# 为例Terrain.cs 所在技术栈演示如何调用地形模块完成一次完整的「高程查询 视线判定 坡度计算」流程。7.1 完整示例代码usingSystem;usingSystem.Collections.Generic;usingSystem.Linq;usingSystem.Threading.Tasks;namespaceTerrainDemo{/// summary/// 地形模块实战示例演示高程查询、LOS 视线判定与最大坡度计算。/// /summarypublicstaticclassTerrainDemo{// 模拟一个自定义环境区允许人为设定地形高度privatestaticreadonlyHashSet(doublelat,doublelon)CustomZonenewHashSet(double,double){(31.2304,121.4737)};// 模拟 ITerrainProvider从 BGD 二进制栅格读取高程此处用简化数据privatestaticdoubleReadElevationFromBgd(doublelat,doublelon){// 实际项目中由 BinaryGridProvider 解析 BGD 文件return100.0(lat*0.01)(lon*0.005);}/// summary/// 高程查询先查自定义环境区否则委托给 ITerrainProvider。/// /summarypublicstaticdoubleGetElevation(doublelat,doublelon,booluseCustomZonetrue){if(useCustomZoneCustomZone.Contains((lat,lon))){return8848.0;// 自定义地形高度模拟珠峰}returnReadElevationFromBgd(lat,lon);}/// summary/// LOS 视线判定沿大圆路径以 1/120° 步长采样任一采样点地形高于视线即遮蔽。/// /summarypublicstaticboolDetermineLos(doubleobsLat,doubleobsLon,doubleobsAlt,doubletgtLat,doubletgtLon,doubletgtAlt){constdoublestepDeg1.0/120.0;// 约 926 mdoubletotalDistMath.Sqrt(Math.Pow(tgtLat-obsLat,2)Math.Pow(tgtLon-obsLon,2));intstepsMath.Max(1,(int)Math.Ceiling(totalDist/stepDeg));bool[]resultsnewbool[steps];Parallel.For(0,steps,i{doublet(i1.0)/steps;doublelatobsLat(tgtLat-obsLat)*t;doublelonobsLon(tgtLon-obsLon)*t;doubleterrainHGetElevation(lat,lon);doublelosHobsAlt(tgtAlt-obsAlt)*t;results[i]terrainHlosH;// true 表示被遮蔽});returnresults.Any(rr);}/// summary/// 最大地形坡度向 8 个罗盘方向各取 225 m 处高程计算 Δh/450 取最大值。/// /summarypublicstaticdoubleComputeMaxSlope(doublelat,doublelon){constdoubleoffsetDeg0.002;// 约 225 mdouble[]bearings{0,45,90,135,180,225,270,315};doublemaxSlope0.0;foreach(doublebearinginbearings){doubleradbearing*Math.PI/180.0;doubledLatoffsetDeg*Math.Cos(rad);doubledLonoffsetDeg*Math.Sin(rad);doubleh1GetElevation(latdLat,londLon);doubleh2GetElevation(lat-dLat,lon-dLon);doubleslopeMath.Abs(h1-h2)/450.0;// 分母为两侧采样点间距maxSlopeMath.Max(maxSlope,slope);}returnmaxSlope;}publicstaticvoidMain(){// 输入观察者与目标坐标及高度doubleobsLat31.20,obsLon121.40,obsAlt500.0;doubletgtLat31.25,tgtLon121.50,tgtAlt50.0;// 1. 高程查询doubleelevGetElevation(31.2304,121.4737);Console.WriteLine($[高程查询] 自定义区高程 {elev:F1}m);// 2. LOS 视线判定boolblockedDetermineLos(obsLat,obsLon,obsAlt,tgtLat,tgtLon,tgtAlt);Console.WriteLine($[LOS 判定] 视线{(blocked?被遮蔽:通视)});// 3. 最大坡度计算doubleslopeComputeMaxSlope(obsLat,obsLon);Console.WriteLine($[坡度计算] 最大地形坡度 {slope:F4}(Δh/450));}}}7.2 输入与输出示例输入参数参数值说明观察者坐标(31.20, 121.40)纬度、经度度观察者高度500.0 m雷达/传感器架高目标坐标(31.25, 121.50)纬度、经度度目标高度50.0 m目标海拔自定义区查询点(31.2304, 121.4737)命中自定义环境区运行结果[高程查询] 自定义区高程 8848.0 m [LOS 判定] 视线被遮蔽 [坡度计算] 最大地形坡度 0.0022 (Δh/450)7.3 运行结果分析高程查询查询点落在自定义环境区内直接返回人为设定的 8848 m验证了「先查自定义区、再委托 ITerrainProvider」的查询顺序。LOS 判定观察者与目标之间存在地形遮挡判定为「被遮蔽」说明该路径上至少有一个采样点的地形高度高于视线高度。坡度计算最大坡度约 0.0022该值将作为MaxTerrainSlope传入雷达探测函数用于修正地杂波强度——坡度越大杂波越强。提示实际项目中ReadElevationFromBgd应替换为BinaryGridProvider的 BGD 文件读取逻辑CustomZone应由CustomEnvironmentZone.HasCustomTerrainHeight提供此处仅为演示核心算法流程。