提取surfer插值结果方法简介

topmad

        经常用到插值，编写程序是很方便，但是，有些插值方法编起来确实很麻烦，ArcGis  ，SURFER 这些软件都可以对空间数据进行多种方式的插值，尤其是surfer，妹啊，为啥每次我一说srufer，就想起来言深深 ---------

---------------------------------------------------------------
********************************************
我是纯洁的分割线，下面抛开言深深，言归正传
********************************************
----------------------------------------------------------------
surfer中不同插值方法中还可以选择一些高级属性，比较方便

不多讲了，做一下利用surfer 提取出插值结果的方法吧

1 数据准备
   将已知数据按经度，纬度，要素值进行排列，保存为txt或excel等surfer可识别格式的文件，如图
  

  将待插值的地理坐标也按照该方式进行排列，待求要素用0补齐
  

2 插值开始鸟....

  打开surfer  
  网格---数据---选择已知数据文件，老五习惯用抽出值，就以它为例，然后插值方法，高级选项，自己选
  

  然后点确认，提示grd文件保存成功



接下来...

点击  数据--残差


然后选择刚刚生成的grd文件



然后会叫你再输入一个文件，这个时候，把刚刚准备好的那个仅有站点信息，没有值的文件选中




然后会出来这个


这个基本可以无视的....

然后点击确认
就有这个了


然后 在菜单栏里面选择 数据--变换  把插出来的残差值乘以负1，如图所示

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



大致做了个东西，希望对大家有用吧

------------------------------------------------------------------------------
***********言深深要出来了，分割线不纯洁了************
------------------------------------------------------------------------------


@言深深 你没事给编程提取一个呗....上次白跟你讨论半天了....哈


写的有错的，不严谨的地方，各位批评指正哈...

陌小默的

哈哈，原来还可以单独提取出来啊~~学习了，深深肯定会给力的！

mofangbao

你什么时候开始弄surfer了啊。。。

topmad

mofangbao 发表于 2012-5-7 21:22

                               
登录/注册后可看大图

你什么时候开始弄surfer了啊。。。

被逼的哇....

陌小默的

topmad 发表于 2012-5-7 21:26

                               
登录/注册后可看大图

被逼的哇....

逼上梁山~

言深深

哈哈，有点儿意思，原来残差就是算差的啊！不错啊，有这个提取程序就不用写了吧

兰溪之水

手机党也支持下…

topmad

mofangbao 发表于 2012-5-7 21:22

                               
登录/注册后可看大图

你什么时候开始弄surfer了啊。。。

哪里加错了

mofangbao

topmad 发表于 2012-5-8 07:56

                               
登录/注册后可看大图

哪里加错了

啥意思？

topmad

mofangbao 发表于 2012-5-8 08:12

                               
登录/注册后可看大图

啥意思？

你说的是体力跟认证币？