fy-2e tbb(相当黑体温度)绘制方法（初步）--期待完善

喜上眉梢

想用GRADS处理云图的想法有好几年了，但始终未得法，看了http://bbs.06climate.com/forum.php?mod=viewthread&tid=2528贴后，深受启示（在此首先向nx花儿-LJQ深表感谢），经过多次周折得到了数据资料源http://satellite.cma.gov.cn/PortalSite/Ord/Satellite.aspx?seriesCode=FY2X&satellitecode=FY2E，不过要先注册，邮箱激活后用户才可以下载到所需数据，不过要等24小时才能到邮箱中（急切中只下了一个时次且仅是TBB的，其余的数据未做试验），nx花儿-LJQ 的CTL改写如下（仅一个时次） ,如下：
dset J:\tbb\FY2E_TBB_IR1_OTG_20100530_0700.AWX
title 2010 twpice fy2e_tbb
options template
options yrev
* yrev表示y轴反向
undef 0.0
xdef 1201 linear  45 0.1
ydef 1201 linear -60 0.1
zdef 1 linear 1000 1
tdef 1 linear 23z30may2010 1hr
vars  1
tbb   1    -1,40,1    Geopotential Heights (m)
endvars
编写了绘图GS ,如下：
*绘制fy2e tbb 云图
'reinit'
'open J:\tbb\tbb2010053007.ctl'
'set mpdset lowres xinjiang_bj1_grads hetian_diqu_border'
'set csmooth on'
'set poli on'
'set mpdraw on'
'set map 13 6 7'
'set xlint 5'
'set ylint 5'
'set xlopts 1 7 0.14'
'set ylopts 1 7 0.14'


'set grads off'
'set grid off'
'set map 13 6 7'
'set xlint 10'
'set ylint 5'
'set xlopts 1 7 0.14'
'set ylopts 1 7 0.14'
'set lon 50 120'
'set lat 20 60'
'set gxout shaded'
****   清风---GSEDIT 降水量颜色 -------------------------------
'set rgb 21 249 252 215'
'set rgb 22 209 249 205'
'set rgb 24 167 250 200'
'set rgb 25 18 215 239'
'set rgb 26 24 160 224'
'set rgb 27 10 134 229'
'set rgb 28 3 121 242'
'set rgb 29 62 44 236'
'set rgb 30 86 31 226'
'set rgb 31 59 25 187'
'set rgb 32 62 3 128'
'set rbcols 21 22 24 25 26 27 28 29 30 31 32'
*****------------------------------------------
**  新颜色(MICAPS 调色板)
**本想偷个懒，使用micaps下面的调色版PAL颜色，确提示出错，未搞懂，期待高手解决！！！
*
*来自   D:\MICAPS3\modual1\cloudgpf\pal
*'set rgb 1     91     10     10'
*'set rgb 2     91     10     10'
*'set rgb 3     91     10     10'
*'set rgb 4     91     10     10'
*'set rgb 5     91     10     10'
*'set rgb 6     91     10     10'
*'set rgb 7     91     10     10'
*'set rgb 8     91     10     10'
*'set rgb 9     92     11     11'
*     。。。。。有很多。。。。。
*'set rgb 255   10     10     10'
*'set rgb 256   0      0      0 '
****
'set cmin 160'
'd tbb'
'cbarn 0.8 1 10.6'
*皮山
'q w2xy 78.283 37.617'
x=subwrd(result,3)
y=subwrd(result,6)
'set strsiz 0.2'   
'wrtchs 'x+0.05' 'y-0.1' 皮山 1.3 2 c'
'printim J:\tbb\tbb2010053007h.emf x1024 y768 white'
'enable print  J:\tbb\tbb2010053007h.gmf'
'print'
'disable print'
'clear'

以上GS出图效果如下：
上图速度好慢。。。。。。。
不足之处颜色设置还是没有达到MICAPS下的效果，有空再整吧，先到这，肚子抗议了，请老大见谅啊。

mofangbao

替你上图了，等过两天回来了我也试试

繁星一点

楼主这个是画单个时次的吧，能不能提示一下画多个时次时gs文件要怎么改一下呢，我自己编的gs文件批处理时老是出错，又查不出错在哪里

传说中的谁

亲，grads里面0-16的颜色值已经被系统占用了，所以你要定义颜色必须得是17开始，到99，也就是set rgb x中的x取值在17~99之间。

传说中的谁

呃~~~忘记了，感谢楼主分享

lacie翼

繁星一点 发表于 2011-12-10 21:27

                               
登录/注册后可看大图

楼主这个是画单个时次的吧，能不能提示一下画多个时次时gs文件要怎么改一下呢，我自己编的gs文件批处理时老 ...

ctl用批处理
grads
i=1
while(i<=n)
'set t 'i
……
'd tbb-173.15'
'print'
'c'
i=i+1
endwhile
'disable print'

飞雪残剑

好东西，收了！我也在做TBB,一起交流下？我的QQ:405305221！

喜上眉梢

自发完后到今天上午网速都不畅，看到大家的指导，我抽空再完善下。

冬日情愫+

请教一下：9210格式和标称格式有什么区别呢？

noisy

尝试出来啦