爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

搜索
查看: 3502|回复: 4

[脚本编辑] 怎么用中国自动站与CMORPH融合的逐时降水量0.1°网格数据集做插值?

[复制链接]
发表于 2016-12-18 17:07:03 | 显示全部楼层 |阅读模式

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

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

x
现有数据集SEVP_CLI_CHN_MERGE_CMP_PRE_HOUR_GRID_0.10-*.grd,怎么用这个来做空间插值啊,要先转化成nc格式么?有哪位高手处理过这个数据集么
密码修改失败请联系微信:mofangbao
发表于 2016-12-19 10:54:40 | 显示全部楼层
这个数据本来就是0.1*0.1分辨率的网格数据,不需要做插值吧?
除非你要改变网格分辨率,那就需要重新插值。
密码修改失败请联系微信:mofangbao
 楼主| 发表于 2016-12-20 14:09:04 | 显示全部楼层
栤唥de嶶笑 发表于 2016-12-19 10:54
这个数据本来就是0.1*0.1分辨率的网格数据,不需要做插值吧?
除非你要改变网格分辨率,那就需要重新插值 ...

最终我们想插值成0.01*0.01度的资料
密码修改失败请联系微信:mofangbao
发表于 2016-12-20 16:07:53 | 显示全部楼层
有相应的ctl吗?如果有,可以用grads中rc函数插值为更高经度的
密码修改失败请联系微信:mofangbao
发表于 2016-12-20 16:22:10 | 显示全部楼层
open语句打开ctl
*pr的经纬度范围和时间、层次设定
'set lon 70 140'   
'set lat 17 55'
'set z 1'
'set t 1 684'
*7100是x方向格点数,70是插值后的起始经度,0.01是经向分辨率,3900是y方向格点数,17是插值后的起始纬度,0.01是纬向分辨率,bl是从粗网格插值向细网格
'define pr=re(precip,7100,'linear',70,0.01,3900,'linear',17,0.01,bl)'
然后用fwriite写数据即可
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

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

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

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