请选择 进入手机版 | 继续访问电脑版
爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 54887|回复: 66

[源代码] 上传一个自己编写的利用反距离权重插值函数

  [复制链接]

新浪微博达人勋

发表于 2012-5-7 11:47:54 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

x
本帖最后由 cnphoenix 于 2012-5-7 11:56 编辑

        这些天正准备对缺测数据进行插值,经过清风的点拨,决定还是进行空间插值所得到的结果较好,于是编写了一个利用反距离权重法计算的点插值函数。
        欢迎讨论!

sub_idw_int.f90

0 Bytes, 下载次数: 235, 下载积分: 金钱 -5

售价: 1 贡献  [记录]

点评

值得买一下  发表于 2012-5-7 14:20

评分

参与人数 5金钱 +38 贡献 +10 体力 +20 收起 理由
wellionx + 1
yingyefang + 2
topmad + 10 + 2 + 20
mofangbao + 15 + 6
言深深 + 10 + 2

查看全部评分

密码修改失败请联系微信:mofangbao

新浪微博达人勋

0
早起挑战累计收入
发表于 2012-5-7 14:20:02 | 显示全部楼层
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2012-5-7 11:56:16 | 显示全部楼层
支持,我童鞋也在 做,估计看到这个要乐了~
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2012-5-7 12:51:10 | 显示全部楼层
个人觉得用BARNES的插值比较方便一点。
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 成长值: 0
发表于 2012-5-7 14:04:37 | 显示全部楼层
反距离权重是不是此前@topmad也做过?
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2012-5-7 14:21:13 | 显示全部楼层
言深深 发表于 2012-5-7 14:04
反距离权重是不是此前@topmad也做过?

恩,原来做过一次,再来学一下这位兄长的代码...
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2012-5-7 16:30:33 | 显示全部楼层
刚刚看了一下楼主的代码,应该是格点数据的插值方法吧  很清晰
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2012-5-7 17:12:59 | 显示全部楼层
topmad 发表于 2012-5-7 16:30
刚刚看了一下楼主的代码,应该是格点数据的插值方法吧  很清晰

我准备用来对日值数据进行插值
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2012-5-7 18:48:27 | 显示全部楼层
用BARNES方法将格点资料插值到站点资料上

  1. DO it=1,NT
  2. DO k=1,LEV
  3.     DO index=1,64
  4.         w=0.0
  5.         DO i=1,NX1
  6.         DO j=1,NY1
  7.             r=(2.5*i+97.5-station(index).lon)**2+(2.5*j+12.5-station(index).lat)**2
  8.             IF (r>6.25) THEN
  9.                 CYCLE
  10.             ELSE
  11.                 station(index).ncep(k,it)=station(index).ncep(k,it)+exp(-r/6.25)*ncep(i,j,k,it)
  12.                 w=w+exp(-r/6.25)
  13.             ENDIF            
  14.         ENDDO
  15.         ENDDO
  16.         station(index).ncep(k,it)=station(index).ncep(k,it)/w
  17.     ENDDO
  18. ENDDO
  19. ENDDO
  20. WRITE(*,*) 'NCEP INTERPOLATION OK!'





因为格点分辨率是2.5*2.5,故选6.25为阈值,对应的数据分辨率改一下阈值就好。
阈值大了,会滤去一些信息;阈值小了,可能有效半径内没有格点值

评分

参与人数 1金钱 +10 贡献 +1 收起 理由
topmad + 10 + 1 给你弄到代码里面

查看全部评分

密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2012-5-7 19:49:09 | 显示全部楼层
易小凯 发表于 2012-5-7 18:48
用BARNES方法将格点资料插值到站点资料上

对于栅格点上的插值不懂,学习了
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表