批量对逐日天气进行反距离加权

  背景介绍  

       在使用作物模型模拟时，需要输入逐日天气数据，但是由于数据获取有限及气象站点分布不均匀，常常会出现农业站点附近10km（或20km）内并不存在气象站点，此时就需要搜寻该农业站点附近最近的若干个气象站点进行插值，从而得到该站点的天气数据。本文已经根据距离远近，获取需要进行插值的20个农业站点最近的10个气象站点信息，同时获取了对应的全部气象站点2000-2010年的每日天气资料，然后利用R语言对逐日天气数据进行反距离加权插值，最后生成DSSAT格式的WTH文件。

  数据介绍  

1)农业站点及最近10个气象站点信息

2)气象站点的天气资料

 代码实现  

  最终结果  

      得到了20个农业站点对应的WTH文件，一个农业站点一年一个文件（前四个字母为该农业站点对应的英文名称，如50888为SRBQ；5-6字符是年份后两位，04表示2004年）。