气象要素空间插值方法优化.pdfVIP

第 23 卷 　第 3 期 地 　　理 　　研 　　究 Vol 23 , No 3 2004 年 5 月 GEO GRAPHICAL 　RESEARCH May , 2004 气象要素空间插值方法优化 封志明 , 杨艳昭 , 丁晓强 , 林忠辉 ( 中国科学院地理科学与资源研究所, 北京 100101) 摘要 : 在区域水土平衡模型的研究中 , 空间插值可提供每个计算栅格的气象要素资料。本文 运用反距离加权法 ( IDW) 和梯度距离反比法 ( GIDW) , 对 1961～2000 年甘肃省及其周围 85 个气象站点的多年平均温度与降雨量进行了内插。交叉验证结果表明 : 对于 IDW 和 GIDW , 二者温度插值的平均绝对误差 ( MAE) 分别为 2 28 ℃和 0 73 ℃, 平均相对误差 (MRE) 分别为 29 02 %和 941 % , 降雨插值的 MAE 值依次为 552mm 和 4 90mm , MRE 值 分别为 1943 %和 1780 % , GIDW 明显优于 IDW 。需要指出的是: 对于降雨 , 当其经纬度和 海拔高程的复相关系数大于 0 80 时 , GIDW 插值结果优于 IDW ; 否则相反。 关 　键 　词 : 距离反比法 ; 梯度距离反比法 ; 降雨; 温度 中图分类号: P4680 　　文章编号 : (2004) 根据已知气象数据对未知点进行插值预测是上世纪以来地理科学研究的一个热点。在 许多重要的国际会议中 , 气象要素插值已被作为专门的议程来进行探讨和研究[1 ] 。空间 插值方法是多种多样的 , 不同的方法插值结果　然不同。因此 , 如何根据数据与区域的特 征选择最优的方法一直是地学研究的一个难题。 对于温度和降雨 , 不同的学者用不同的方法进行了空间内插。主要有泰森多边形法、 距离反比法、克里格法、梯度距离平方反比法等。泰森多边形法算法最为简单 , 但容易生 成专题与专题之间明显的边界 ; 克里格法由南非地质学家 Krige 发明 , 其实质是最佳线性 无偏估计 , 应用广泛[2～5 ] 。但该方法无论是变程的确定、半变异模型的拟合 , 还是邻域 的选取 , 均耗时耗力 , 且有较强的主观性[6～8 ] ; 距离反比法由美国国家气象局于 1972 年 提出 , 其算法较为简单 , 应用也很广泛[9 ] 。当数据不存在各向异性时 , 其插值结果优于 克里格法[10 ] ; 梯度距离平方反比法物理意义明确 , 算法简单 , 应用广泛 , 当气象要素与 海拔高程和经纬度相关性较强时 , 其插值精度较高。 在作者进行区域水土平衡模型的研究中 , 其第一步就是对温度、降雨等气象数据进行 合理的插值。考虑到时间与经费的限制 , 力图寻求一种简单、灵活、实用的方法。为此 , 我们采用距离平方反比法和梯度距离平方反比法对甘肃省的降雨和温度数据进行了内插 , 并且提出了两种方法的优化方案 , 最后对各种插值方法进行了验证和选择。 1 　资料与方法 11 　数据分析与处理方法 本文以甘肃省为研究区域 , 选取甘肃省及其周围 5 个省区 129 个气象台站作为样本 　　收稿日期 : ; 修订日期 : 　　基金项目 : 中国科学院知识创新工程重要方向项目 ( KZCX3 - SW - 333) (