基于gis的浙江省年平均气温空间插值评估.docxVIP

基于gis的浙江省年平均气温空间插值评估 在空间插值方面的应用 全球气候变化引起了人们的广泛关注。气温是地球表层系统的综合环境指数，也是全球气候变化研究的重要参数。在全球变化研究的各个领域中得到了广泛应用。然而,由于人力、财力等因素的限制,现阶段各地建立的气象观测系统非常有限,且空间位置分布不均,很难满足各领域所要求的时间和空间分辨率的气温数据产品。因此,采用有限的气象站点对气象要素进行空间推算研究显得尤其重要。 1980s初,常采用数理公式来推算和模拟气温的空间分布,这些方法虽能够反映气温要素的整体变化趋势,但由于各地地形复杂,造成局部误差较大,尤其是在冬季,最大误差可达5 ℃以上,且计算繁琐,工作量大,制图困难。1990s后期,GIS技术的应用日趋成熟,与此同时逆距离权重法、普通克里金、样条插值法和趋势面等空间插值技术开始被用于气象要素的推算。林忠辉等、于贵瑞等、李军等、庄立伟等、于飞等、蔡福等等[12,13,14,15,16,17]学者运用上述插值方法分别对全国、东北和西南等不同区域的气象要素空间分布特征进行了研究。然而,这些研究所选用的空间插值方法,只适合于其所对应的研究区,并不能广泛的应用于其他地形地貌复杂的地区。 浙江省位于中国东部沿海,地貌类型多样,西南部地势高山脉较多,中部为丘陵地带,大小盆地分布于丘陵之间,东北部为冲积平原地势低平。多样化的地形条件,造成了区域内热量资源的空间分布差异很大。而且,至今未见有关基于空间插值方法的浙江省平均气温分析的文献。因此,本文利用研究区域内72个气象站点近38 a的年平均气温资料,采用逆距离权重法、普通克里金、样条法和多元回归方法进行空间插值,通过误差分析来对插值结果进行精度评价,旨在从中选出适宜于浙江省年平均气温空间插值的方法,制作浙江省年平均气温的精细化区划图,并分析了浙江省年平均气温的空间分布特征。 1 数据与方法的研究 1.1 数据来源和说明 研究区共有72个气象观测站,基本上是均匀分布于各地(见图1)。本文所用数据为浙江省72个气象观测站1971—2008年的逐日平均气温、72个气象观测站点的地理属性(经度、纬度和海拔高度)、浙江省1∶25万的行政区划图和1∶25万的数字化地形图。 1.2 网格化dem 对上述数据按以下步骤进行预处理:在1∶25万数字化地形图的基础上,运用ARCGIS9.2中的CREATETIN命令生成矢量化的地形不规则三角网(TIN),选用TINLATTICE命令进行空间内插生成浙江省的DEM,栅格大小为100 m×100 m。然后从DEM数据中提取出经度、纬度栅格图。本文旨在研究年平均气温的空间分布特征,故用1971—2008年的逐日平均气温资料计算出38 a的年平均气温数据。根据均匀分布的原则,从72个气象站点中选择62个气象站点,建立年平均气温的空间分布模型,再用另10个气象站点的年平均气温数据作为验证数据,用以评估模型的精度。 1.3 年平均气温的空间模拟 影响气温的因素是多方面的,其空间插值方法也有很多。本文在前人研究基础上,结合浙江本省的地域特性,采用了逆距离权重法、普通克里金、样条法和多元回归法共4种方法来对年平均气温进行空间模拟。 1.3.1 未开明点预测算法 简称IDW) IDW是基于“地理第一定律”的基本假设,即两个物体的相似性随它们间的距离增大而减小。假设已知样点对预测点值的预测都有局部性影响,其影响随距离增加而减小。离预测点近的已知样点在预测过程中所占权重大于离预测点远的已知样点。计算公式如下 Ζ(so)=Ν∑i=1λiΖ(si)Z(so)=∑i=1NλiZ(si), (1) 式中Z(so)为so处的预测值,Z(si)为在si处测量值;N为预测过程中要使用的预测点周围样点的数量;λi为预测过程中已知样点的权重。权重计算公式 λi=d-pio/Ν∑i=1d-pioλi=d?pio/∑i=1Nd?pio, (2) 式中dio为预测点so与各已知样点si之间的距离。这里p的最佳值通过求均方根预测误差的最小值求得,一般情况下该值取2。 1.3.2 确定权重系数 普通克立格法是以空间自相关性为基础,利用原始数据和半方差函数的结构性,对区域化变量的未知采样点进行线性无偏最优估计的一种方法。计算公式为 Ζ=n∑i=1λiΖ(xi)Z=∑i=1nλiZ(xi), (3) 其中,Z为待估点的预测值,λi为权重系数,n为已知的观测点数目,Z(xi)为观测点的观测值。 为满足无偏性和最优性两个条件,通过建立如下克立格方程组来确定权重系数。 {n∑j=1λjC(vi,vj)-μ=C(vi,V)n∑i=1λi=1?????????????∑j=1nλjC(vi,vj)?μ=C(vi,V)∑i=1nλi=1, (4) 其中,C(vi,vj)为观测点