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差值提取结果，右边是自己写的那个反距离的代码算出来的,  懒得分析了....有空再说吧