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