爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

搜索

【已解决】grd用cdo转成nc格式,但是grd有pdef参数

[复制链接]
 楼主| 发表于 2018-3-16 09:47:12 | 显示全部楼层
下是一种 发表于 2018-3-15 15:05
最近也在用CDO将grb转nc,按照楼主的方法试了,但是一直有错误。
我的ctl文件是图1。显示错误是图2,是idx ...

没用过idx文件呢,不过看报错好像是cdo不支持对idx文件的转换,如果实在不行就只能用其他方法了
密码修改失败请联系微信:mofangbao
发表于 2018-3-18 09:40:19 | 显示全部楼层
楼主 您好  我用你说的办法 删除ctl文件中 pdef(pdef  147   79 lcc   32.00  110.00   74.00   40.00   30.00   60.00  110.00  60000.  60000.) 然后修改xdef ydef为模式格点数,用cdo转化模式的lambert 格点数据 为 nc数据 但最后画出的图形不对 其中
res@mpProjection           = "LambertConformal"
  res@mpLambertParallel1F    = 31.
  res@mpLambertParallel2F    = 32.
  res@mpLambertMeridianF     = 110.
上面三个数据我猜的
我看你的例子数据是 有经纬度数据,而我的原始二进制数据没有经纬度,是不是只要定义LambertParallel1 FLambertParallel2 FLambertMeridianF 这三个就可以了 不要定义全部的经纬度信息
密码修改失败请联系微信:mofangbao
发表于 2018-3-18 11:19:16 | 显示全部楼层
我的原始经纬度数据是二维的 且是不均匀的经纬度数据
密码修改失败请联系微信:mofangbao
 楼主| 发表于 2018-3-18 19:01:14 | 显示全部楼层
云大小子 发表于 2018-3-18 09:40
楼主 您好  我用你说的办法 删除ctl文件中 pdef(pdef  147   79 lcc   32.00  110.00   74.00   40.00   3 ...

pdef  147   79 lcc   32.00  110.00   74.00   40.00   30.00   60.00  110.00  60000.  60000.这些就包含了res@mpLambertParallel1F  
  res@mpLambertParallel2F  
  res@mpLambertMeridianF  这三个参数的设置,建议你可以百度一下grads的pdef的用法、
密码修改失败请联系微信:mofangbao
 楼主| 发表于 2018-3-18 19:02:31 | 显示全部楼层
云大小子 发表于 2018-3-18 11:19
我的原始经纬度数据是二维的 且是不均匀的经纬度数据

没有试过没有经纬度信息画图,但是你可以根据起始经纬度,和经纬度间隔自己创建一个经纬度数组
密码修改失败请联系微信:mofangbao
发表于 2018-3-18 22:36:40 | 显示全部楼层
心如止水的饭团 发表于 2018-3-18 19:02
没有试过没有经纬度信息画图,但是你可以根据起始经纬度,和经纬度间隔自己创建一个经纬度数组

我看了您的例子 pdef这一行我理解。你的原始的经纬度信息lat和lon都是二维的,由于是lambert投影等距离数据,所以经纬度信息不是等间距。当删除pdef,然后修改xdef和ydef为网格数,用cdo转化为nc数据后,这个nc数据的经纬度应该经过插值后就变成等间距的吧,原来的经纬度数据用不了,你是如何定义经纬度数据然后画图,能帮我解释一下吗。
密码修改失败请联系微信:mofangbao
发表于 2018-3-19 08:55:23 | 显示全部楼层
心如止水的饭团 发表于 2018-3-18 19:01
pdef  147   79 lcc   32.00  110.00   74.00   40.00   30.00   60.00  110.00  60000.  60000.这些就包 ...

我看到手册上的说明了 确实 pdef直接有这三个参数
密码修改失败请联系微信:mofangbao
 楼主| 发表于 2018-3-19 17:02:29 | 显示全部楼层
云大小子 发表于 2018-3-18 22:36
我看了您的例子 pdef这一行我理解。你的原始的经纬度信息lat和lon都是二维的,由于是lambert投影等距离数 ...

不知道你的ctl在未修改前,里面是否有xdef 和ydef的值,这些值我想就是他们的起始经纬度,经纬度格点数,和格点间隔大小。我想按照这些数据来作为你画图的经纬度应该是没问题的
密码修改失败请联系微信:mofangbao
发表于 2018-3-23 11:04:58 | 显示全部楼层
你好 我按照你的方法能转出nc数据,但是我用ncl画图画不出lambert投影图,我在学习阶段,你能共享一下你的绘图脚本吗,谢谢。我用ncl画图的时候总是提示res@lambert相关的属性不可用。D:\a.png
密码修改失败请联系微信:mofangbao
发表于 2018-3-23 11:08:59 | 显示全部楼层
我能画出图,但是 图形不是lambert的扇形,且有警告:warning:mpProjection is not a valid resource in lcnative_contour at this time
warning:mpLambertMeridianF is not a valid resource in lcnative_contour at this time
warning:mpLambertParallel1F is not a valid resource in lcnative_contour at this time
warning:mpLambertParallel2F is not a valid resource in lcnative_contour at this time
warning:mpDataBaseVersion is not a valid resource in lcnative_contour at this time

密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

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

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

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