高原生物量图改按海拔连边:局部插值R²提高0.048,30公里外反而更差
核心概要
青藏高原东南缘的地上生物量估计中,按地理最近邻连边的空间图会把谷底森林与相邻高山草甸连在一起,因为水平距离近并不等于生态相似。研究先量化这种错连,再用海拔差阈值筛边重建图,在6638个样本、10个随机种子上把局部插值的R²从0.362提高到0.409,提升0.048。但在30公里缓冲留出的严格外推下,同一做法比未筛边的图低0.054,说明收益依赖测试点附近仍有训练邻居。
深度剖析
空间图不再默认“离得近就相似”,而是先检查每条边两端的海拔与生物量差,再决定是否保留。 此前的标准做法是按经纬度取k个最近邻连边,隐含假设水平距离能代表目标相似性;在海拔146至5383米的复杂地形里,这一假设会失效。 在6638个干净样本、k=10的66380条有向边上,21.41%的边两端海拔差超过800米,12.24%的边两端生物量差超过70兆克每公顷,而这些高生物量差边中有60.6%连接不同的土地覆盖类型。
筛边规则只用一条:两端海拔差不超过阈值τ的边才保留,其余删除。 这属于在构图阶段修正拓扑,而不是在模型结构上补救;海拔作为约束直接作用于边,而不是作为额外输入特征交给模型自行学习。 在10个随机种子的随机划分下,τ=250米时保留60698条边、仅6个孤立节点,R²为0.409±0.023,而未筛边的地理k最近邻图为0.362±0.019,配对差0.048,p=6.0×10⁻⁷,95%置信区间[0.039,0.056],十个种子方向一致。
同一套筛边在严格空间外推下方向反转,收益变成损失。 随机划分下测试点附近仍有训练邻居,图能把局部知识传过去;缓冲划分把测试点周围30或50公里内的训练样本全部移除,这种传递被切断。 30公里缓冲下筛边图比未筛边图低0.054(p=0.027),50公里缓冲下低0.083(p=0.007),十个种子方向一致;同一协议下局部支持指数在30公里处为0。
启示与展望
结论建立在青藏高原东南缘一个区域、2021年单年快照和最小间距1.19公里的固定采样密度上,阈值250米是该区域与采样密度下的经验取值。适用场景是采样域内的局部插值,例如区域内加密与空缺填补;严格外推场景下筛边图并不优于未筛边图。
作者指出阈值τ=250米是区域与采样密度特定的,跨区域可迁移性尚待多区域验证;GraphSAGE骨干下提升不显著(−0.001,p=0.642),说明收益依赖消息传递机制。生物量标签并非野外观测,而是由雷达与激光雷达参考数据经随机森林升尺度并空间聚合得到的派生产品,部分增益可能来自与标签平滑模式的对齐。
