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

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 7754|回复: 10

[图形美化] GrADS 读 FY2G PRE AWX 出图时出现问题

[复制链接]

新浪微博达人勋

发表于 2017-7-7 16:46:41 | 显示全部楼层 |阅读模式

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

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

x
本帖最后由 听海落雪 于 2017-7-7 16:46 编辑

我想用 FY2G有一小时降水估计资料画图,参考了论坛里的贴子,编写出了ctl和gs。在GrADS打开ctl,直接 'd pre',能够出正常的图,但是运行我编写的gs文件出图明显有问题。希望大神们不吝赐教。

ctl 如下:

dset F:\Data\GrADS\2017\20170624\PRE\FY2G_PRE_001_OTG_2017%m2%d2_%h200.AWX
title 2017 FY2G_PRE
options template
options yrev
* yrev表示y轴反向
undef 0.0
xdef 1001 linear 55 0.1
ydef 1001 linear -50 0.1
zdef 1 levels 850
tdef 48 linear 00z23JUN2017 1hr
vars 1
pre 0 -1,40,1 FY2G_PRE(m)
endvars



gs如下:

'reinit'

'open F:\Data\GrADS\2017\20170624\PRE\FY2G_PRE_001.ctl'

'set grads off'
'set grid off'

'set lon 96 106'
'set lat 25 33'

****** 1小时降水 ******
'set rgb 54 255 255 255'  
'set rgb 55 220 250 211'  
'set rgb 56 166 242 143'  
'set rgb 57 61 186 61'   
'set rgb 58 97 184 255'   
'set rgb 59 0 0 225'      
'set rgb 60 250 0 250'   
'set rgb 61 128 0 64'

********* 一页四图 *********

p=1
while(p<=4)

'run page q'p''
'set t 'p''

'set parea 0.5 9.0    1.0 8.0'
'set xlopts 1 5 0.26'
'set ylopts 1 5 0.26'
'set xlab on'
'set ylab on'
'set xlabs 96`3.|98`3.|100`3.|102`3.|104`3.|106`3.`1E'
'set ylabs 25`3.|27`3.|29`3.|31`3.|33`3.`1N'

'set gxout shaded'
'set csmooth on'
'set clevs 0 0.1 1.5 7 15 40 50'
'set ccols 54 55 56 57 58 59 60 61'
'd smth9(pre)'

'set mpdset country1'
'set map 15 1 2'
'draw map'
'set mpdset xnqy_sj'
'set map 15 1 2'
'draw map'
'set mpdset pxdq_dj'
'set map 8 1 5'
'draw map'

'query w2xy 102.2147369385 27.9398202896'
x1=subwrd(result,3)
y1=subwrd(result,6)
'set string 1 c 5 0'
'set strsiz 0.16'
'draw mark 1  'x1' 'y1'  0.25'

'query w2xy 103.4 32.7'
x11=subwrd(result,3)
y11=subwrd(result,6)
'set string 1 c 5 0'
'set strsiz 0.28'
'draw string 'x11' 'y11'  `2  'p+13'BJZ23JUN'

p=p+1
endwhile

************************************************

'printim F:\Data\GrADS\2017\20170624\pic\PRE_2314-2317.png white'





PRE_2314-2317.png

FY2G_PRE_001_OTG_20170623_0000.AWX

983.4 KB, 下载次数: 1, 下载积分: 金钱 -5

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

新浪微博达人勋

发表于 2017-7-7 21:13:40 | 显示全部楼层
你把这句'set clevs 0 0.1 1.5 7 15 40 50'先注释掉看看
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2017-7-7 22:14:25 | 显示全部楼层
river 发表于 2017-7-7 21:13
你把这句'set clevs 0 0.1 1.5 7 15 40 50'先注释掉看看

还是不行,颜色变了。
之前看到一篇网页,是这样写的“……再说CTT,PRE等一些因存储格式导致grads不能顺利读出的产品,原因有老师在之前帖子中已分析。由于强迫症同类产品想出一致的图,我的做法比较简单粗暴,该写AWX成grads可读即可。给出matlab处理的部分脚本……”


http://bbs.06climate.com/forum.php?mod=viewthread&tid=46113
PRE_2314-2317.png
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2017-8-19 15:07:20 | 显示全部楼层
AWX格式前面有说明信息,数据只能重新读写一下,过滤掉说明信息,否则grads读取不了
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2017-8-19 19:12:01 | 显示全部楼层
rabin_xu 发表于 2017-8-19 15:07
AWX格式前面有说明信息,数据只能重新读写一下,过滤掉说明信息,否则grads读取不了

您说是 CTL 文件?
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2017-8-21 08:30:29 | 显示全部楼层
你可以看看3.1,那里比较清楚

气象卫星分发产品及其格式规范AWX2.1.pdf

2.06 MB, 下载次数: 82, 下载积分: 金钱 -5

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

新浪微博达人勋

发表于 2017-8-21 08:31:15 | 显示全部楼层
还要扣积分啊,那自己去这个连接下载吧
http://satellite.nsmc.org.cn/PortalSite/StaticContent/DocumentDownload.aspx?TypeID=10
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2017-8-22 09:02:55 | 显示全部楼层
简单粗暴的那个帖子是我写的。。。脚本不知道放哪儿了,现在FY2G PRE AWX读取也出岔子了{:eb303:}{:eb303:}{:eb303:}
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2017-8-22 14:53:32 | 显示全部楼层
应该是整出来了,matlab脚本留在这里备份一下

list=dir(['./','FY2*.AWX']);   %必须与CTA对应
k=length(list);
for i=1:k
     str= strcat ('./', list(i).name);
end

file=list(i).name;
fid = fopen(file,'rb');
test = fread(fid,[1001 503],'ubit16');% 读取TBB数据
data=test(:,3:503);
for m = 1:1001
    for n = 1:501
        if(data(m,n)==64546)
        data(m,n)=0;
        end        
    end
end
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2017-8-22 16:24:37 | 显示全部楼层
rabin_xu 发表于 2017-8-21 08:30
你可以看看3.1,那里比较清楚

好的,谢谢
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

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

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

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