? 背景介紹??

? ? ? ?在使用作物模型模擬時(shí)，需要輸入逐日天氣數(shù)據(jù)，但是由于數(shù)據(jù)獲取有限及氣象站點(diǎn)分布不均勻，常常會(huì)出現(xiàn)農(nóng)業(yè)站點(diǎn)附近10km（或20km）內(nèi)并不存在氣象站點(diǎn)，此時(shí)就需要搜尋該農(nóng)業(yè)站點(diǎn)附近最近的若干個(gè)氣象站點(diǎn)進(jìn)行插值，從而得到該站點(diǎn)的天氣數(shù)據(jù)。本文已經(jīng)根據(jù)距離遠(yuǎn)近，獲取需要進(jìn)行插值的20個(gè)農(nóng)業(yè)站點(diǎn)最近的10個(gè)氣象站點(diǎn)信息，同時(shí)獲取了對(duì)應(yīng)的全部氣象站點(diǎn)2000-2010年的每日天氣資料，然后利用R語(yǔ)言對(duì)逐日天氣數(shù)據(jù)進(jìn)行反距離加權(quán)插值，最后生成DSSAT格式的WTH文件。

? 數(shù)據(jù)介紹??

1)農(nóng)業(yè)站點(diǎn)及最近10個(gè)氣象站點(diǎn)信息

2)氣象站點(diǎn)的天氣資料

?代碼實(shí)現(xiàn)??

? 最終結(jié)果??

? ? ? 得到了20個(gè)農(nóng)業(yè)站點(diǎn)對(duì)應(yīng)的WTH文件，一個(gè)農(nóng)業(yè)站點(diǎn)一年一個(gè)文件（前四個(gè)字母為該農(nóng)業(yè)站點(diǎn)對(duì)應(yīng)的英文名稱，如50888為SRBQ；5-6字符是年份后兩位，04表示2004年）。

?