【新手求助】地形灰度图叠加降水等值线图

mlittlej

试了好多办法想在地形底图上叠加降水，这个还好看点。不过现在这个gs有两个问题1.grads可以画出来，但是在我保存的路径里图片打不开无法预览。2.等值线是灰色的，我想让它显示彩色的。ps降水是站点插值的。拜托大家帮我看看娜丽不行，或者还有什么别的办法地形叠加降水好看点。谢谢啦！！！！！
'reinit'
'open d:/opengrads/topo10.ctl'
'enable print e:/lmj/diejia.png white'
'set parea 1 10 1 8'
'set mpdset hire'
'set map 1 1 5'
'set xlopts 1 5.5 0.13'
'set ylopts 1 5.5 0.13'
'set grads off'
'set csmooth on'
'set lon 73 97'
'set lat 34 50'
'set mpdset hires'
'd:/opengrads/grey.gs'
'set gxout shaded'
'd tp'
'cnbasemap tp'
'cbarn'
'print'
'disable print'
'close 1'
'open d:/opengrads/grid.ctl'
'open d:/opengrads/sta.ctl'
'enable print e:/lmj/diejia.png white'
'set lon 73 97'
'set lat 34 50'
'set gxout contour'
'set ccolor rainbow'
'set cmin 0'
'set grads off'
'set grid off'
'd oacres(g,rain.2)'
'province-basemap xinj'
'print'
'disable print'
'close 1'
;

river

首先既然是叠加，那就是画在一张图片上，为什么要分开使用两个enable print 命令。其次你要画的只是中国的一个省份，覆盖这个省份以外的地方使用province-basemap 足够了，不需要先使用basemap命令再使用province-basemap命令。最后就是楼主先理解一下grads的绘图的基础知识，尤其是enable print 命令

mofangbao

楼上说的是对的。
还有你黑白的是不是由于这个gs引起（'d:/opengrads/grey.gs'）的，感觉楼主还不太懂GrADS的基础，多看看基础知识吧

mlittlej

river 发表于 2013-7-17 23:34

                               
登录/注册后可看大图

首先既然是叠加，那就是画在一张图片上，为什么要分开使用两个enable print 命令。其次你要画的只是中国的一 ...

'reinit'
'open d:/opengrads/topo10.ctl'
'set parea 1 10 1 8'
'set mpdset hires'
'set map 1 1 5'
'set xlopts 1 5.5 0.13'
'set ylopts 1 5.5 0.13'
'set grads off'
'set csmooth off'
'set lon 73 97'
'set lat 34 50'
'set mpdset hires'
'd:/opengrads/grey.gs'
'set gxout shaded'
'd tp'
'cnbasemap tp'
'cbarn'
'printim e:/lmj/diejia.gif gif white'
'close 1'
'open d:/opengrads/grid.ctl'
'open d:/opengrads/sta.ctl'
'set lon 73 97'
'set lat 34 50'
'set gxout contour'
'set ccolor rainbow'
'set cmin 0'
'set grads off'
'set grid off'
'd oacres(g,rain.2)'
'province-basemap xinj'
'printim e:/lmj/diejia.gif gif white'
;
现在能画出来图，但是等值线还是灰度的
*脚本由GrADS调色盘自动生成-气象家园论坛
'set rgb 21 255 255 255'
'set rgb 22 234 234 234'
'set rgb 23 223 223 223'
'set rgb 24 193 193 193'
'set rgb 25 172 172 172'
'set rgb 26 150 150 150'
'set rgb 27 129 129 129'
'set rgb 28 108 108 108'
'set rgb 29 88 88 88'
'set rgb 30 66 66 66'
'set rgb 31 45 45 45'
'set rgb 32 25 25 25'
'set rgb 33 170 170 170'
'set rbcols 21 22 23 24 25 26 27 28 29 30 31 32 33'
return
这个是grey.gs
是不是还要set一下color
真的不好意思，确实基础比较差。。。麻烦再帮我看看。。。谢谢了~~~

mlittlej

mofangbao 发表于 2013-7-18 08:12

                               
登录/注册后可看大图

楼上说的是对的。
还有你黑白的是不是由于这个gs引起（'d:/opengrads/grey.gs'）的，感觉楼主还不太懂GrAD ...

确实是临时急着用才看的。。。这几天也努力的翻论坛了。。。发帖救助也觉得挺不好意思的。。。抱歉麻烦了啊。。。

river

mlittlej 发表于 2013-7-18 08:36

                               
登录/注册后可看大图

'reinit'
'open d:/opengrads/topo10.ctl'
'set parea 1 10 1 8'

cnbasemap命令本身就在带了display了，所以把前面一句'd tp'去掉试一下

mlittlej

river 发表于 2013-7-19 16:08

                               
登录/注册后可看大图

cnbasemap命令本身就在带了display了，所以把前面一句'd tp'去掉试一下

乱试竟然被我试出来了！！！谢谢你啦！！！

schliezer

mlittlej 发表于 2013-7-18 08:36

                               
登录/注册后可看大图

'reinit'
'open d:/opengrads/topo10.ctl'
'set parea 1 10 1 8'

我的问题恰好相反，灰度图总无法输出。在图形窗口是灰度图，但输出到文件夹里以后就成了彩色的，你知道怎么弄吗？附gs：
'reinit'
'open D:\gradsdata\ncepfnl_1_6h\fnl_20120719_00_00_c.ctl'
'set grid off '
'set lat 20 70'
'set lon 50 140'
'set display greyscale white'
'set csmooth on '
'set mpdset cnworld'
t=1
while(t<=13)
   'c'
   'set t 't''
   'set gxout shaded'   
   'set lev 200'   
   'd mag(UGRDprs,VGRDprs)'   'cbarn 1 0 5.5 8'
   'set gxout contour'
   'set cint 4'
   'set ccolor 1'
   'set cthick 4'
   'set lev 500'
   'd hgtprs/10'
   'set gxout barb'  
   'set ccolor 1'
   'set lev 700'
   'd skip(UGRDprs,5);skip(VGRDprs,5)'
   'printim  D:\gradsdata\ncepfnl_1_6h\0719图\叠加图\两层叠加'%t%'.gif  white'
   t=t+1
endwhile
;

mlittlej

schliezer 发表于 2014-2-11 11:14

                               
登录/注册后可看大图

我的问题恰好相反，灰度图总无法输出。在图形窗口是灰度图，但输出到文件夹里以后就成了彩色的，你知道怎 ...

才看到，不好意思啊。。。你试试那个grads调色板，我就是用的那个。要不我把gs给你？

schliezer

mlittlej 发表于 2014-3-6 14:05

                               
登录/注册后可看大图

才看到，不好意思啊。。。你试试那个grads调色板，我就是用的那个。要不我把gs给你？

那我看看你的gs吧，谢谢！