surfer插值方法及提取插值結果 轉載

 

 

 

 

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有點不同,不夠人性化啊......)


看看 結果就出來了,上面一張圖是我已經算出來結果,然後在把變換對話框調出來的

然後,就沒有然後了

其實吧,還是有一點,這個插出來的結果,我其實是從一個數據組裏面抽出來的,原值當然是有的,閒來無事,跟原來的值對比了一下,還是有一點點差異的(降水資料,貌似這樣做有點沒意思)

剛剛是閒來無事,現在更閒,就把上面的東西,用當時自己寫的一個反距離的代碼算了一下,不同的結果對比一下


  1. 31.485714    28.82906455    29.8017
  2. 31.257143    29.40234631    30.01523
  3. 26.615385    29.11981737    29.01336
  4. 30.1875    29.97364826    30.54164
  5. 35.415385    32.25222953    33.02903
  6. 35.021053    31.64563599    32.40569
  7. 27.864706    31.09236399    30.74112
  8. 29    30.54069527    30.97822
  9. 23.3    31.68735399    32.36668
  10. 34.275    29.93878857    29.98241
  11. 29.916667    29.22721679    29.61226
  12. 26.088889    28.44353863    28.72874
  13. 29.23    28.62757898    28.55389
  14. 33.566667    28.54992658    28.48255
  15. 32.961538    32.12949803    32.5848
  16. 28.433333    28.70826568    28.04766
  17. 55.24    33.22774589    35.11396
  18. 31.66    31.22388293    30.44315
  19. 26.95    28.88664391    28.44494
  20. 30.15    28.23905972    27.22583
  21. 32.288889    28.30212281    28.1416
複製代碼


左邊是原始值,中間是surfer差值提取結果,右邊是自己寫的那個反距離的代碼算出來的,  懶得分析了....有空再說吧

 

發表評論
所有評論
還沒有人評論,想成為第一個評論的人麼? 請在上方評論欄輸入並且點擊發布.
相關文章