http://bbs.06climate.com/forum.php?mod=viewthread&tid=7345
1 數據準備
將已知數據按經度,緯度,要素值進行排列,保存爲txt或excel等surfer可識別格式的文件,如圖
<ignore_js_op>
將待插值的地理座標也按照該方式進行排列,待求要素用0補齊
<ignore_js_op>
2 插值開始鳥....
打開surfer
網格---數據---選擇已知數據文件,老五習慣用抽出值,就以它爲例,然後插值方法,高級選項,自己選
<ignore_js_op>
然後點確認,提示grd文件保存成功
<ignore_js_op>
接下來...
點擊 數據--殘差
<ignore_js_op>
然後選擇剛剛生成的grd文件
<ignore_js_op>
然後會叫你再輸入一個文件,這個時候,把剛剛準備好的那個僅有站點信息,沒有值的文件選中
<ignore_js_op>
然後會出來這個
<ignore_js_op>
這個基本可以無視的....
然後點擊確認
就有這個了
<ignore_js_op>
然後 在菜單欄裏面選擇 數據--變換 把插出來的殘差值乘以負1,如圖所示
<ignore_js_op>
(這個表格跟excel有點不同,不夠人性化啊......)
看看 結果就出來了,上面一張圖是我已經算出來結果,然後在把變換對話框調出來的
然後,就沒有然後了
其實吧,還是有一點,這個插出來的結果,我其實是從一個數據組裏面抽出來的,原值當然是有的,閒來無事,跟原來的值對比了一下,還是有一點點差異的(降水資料,貌似這樣做有點沒意思)
剛剛是閒來無事,現在更閒,就把上面的東西,用當時自己寫的一個反距離的代碼算了一下,不同的結果對比一下
- 31.485714 28.82906455 29.8017
- 31.257143 29.40234631 30.01523
- 26.615385 29.11981737 29.01336
- 30.1875 29.97364826 30.54164
- 35.415385 32.25222953 33.02903
- 35.021053 31.64563599 32.40569
- 27.864706 31.09236399 30.74112
- 29 30.54069527 30.97822
- 23.3 31.68735399 32.36668
- 34.275 29.93878857 29.98241
- 29.916667 29.22721679 29.61226
- 26.088889 28.44353863 28.72874
- 29.23 28.62757898 28.55389
- 33.566667 28.54992658 28.48255
- 32.961538 32.12949803 32.5848
- 28.433333 28.70826568 28.04766
- 55.24 33.22774589 35.11396
- 31.66 31.22388293 30.44315
- 26.95 28.88664391 28.44494
- 30.15 28.23905972 27.22583
- 32.288889 28.30212281 28.1416
左邊是原始值,中間是surfer差值提取結果,右邊是自己寫的那個反距離的代碼算出來的, 懶得分析了....有空再說吧