爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 3799|回复: 5

[图形美化] 如何将以下这两个gs中的图叠加在一起,大神求指教~

[复制链接]

新浪微博达人勋

发表于 2018-1-11 16:28:44 | 显示全部楼层 |阅读模式

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

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

x
'reinit'
'reinit'
'sdfopen e:/lcx/ncdata/air.2011.nc'
'sdfopen e:/lcx/ncdata/air.2012.nc'
'sdfopen e:/lcx/ncdata/air.2013.nc'
'sdfopen e:/lcx/ncdata/air.2014.nc'
'sdfopen e:/lcx/ncdata/air.2015.nc'


'enable print e:/lcx/ncdata/wenzhou/air.gmf'

alt=471.8
yy=2015
mm=nov
dd=4

'set lon 73 136';'set lat 3 54'
'set lev 'alt''
'set time 06z'dd''mm''yy''
'set cint 5'
'd air.5'
'print'


'disable print'
'reinit'
'reinit'


这个是输出温度的

'reinit'
'reinit'
'open E:/lcx/frq.locate.ctl'

'set grads off'
'set parea 0.8 7.5 5.9 10.9'
'set vpage 0 8.5 0 11'
'set mproj scaled';
'set grads off';'set grid off'
'set lon 70 140';'set lat 15 55'
*'set xlpos -10';'set ylpos -10'



'set cmin 100000';'d loc'
'southsea_value loc'

'set lon 70 140';'set lat 15 55'
'set cmin 100000';'d loc'
id = 1
while(id<=7000)
    lonlat=read('E:/lcx/wenzhou.txt')
    lonlat=sublin(lonlat,2)
    ln=subwrd(lonlat,1);lt=subwrd(lonlat,2)

    'q w2xy 'ln' 'lt''
    x0=subwrd(result,3);y0=subwrd(result,6)
    'set line 2'
    'draw mark 3 'x0' 'y0' 0.03'
    id=id+1
endwhile

'printim E:/lcx/wenzhou.png white'

'reinit'
'reinit'

这个是输出坐标点的

现在我想把坐标点叠加在温度图上面,然后进行分析,如何把这两个图叠加在一块儿啊

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

新浪微博达人勋

发表于 2018-1-11 18:14:20 | 显示全部楼层
'reinit'
'sdfopen e:/lcx/ncdata/air.2011.nc'
'sdfopen e:/lcx/ncdata/air.2012.nc'
'sdfopen e:/lcx/ncdata/air.2013.nc'
'sdfopen e:/lcx/ncdata/air.2014.nc'
'sdfopen e:/lcx/ncdata/air.2015.nc'

'set grads off'
'set parea 0.8 7.5 5.9 10.9'
'set vpage 0 8.5 0 11'
'set mproj scaled';
'set grads off';'set grid off'

alt=471.8
yy=2015
mm=nov
dd=4

'set lon 73 136';'set lat 3 54'
'set lev 'alt''
'set time 06z'dd''mm''yy''
'set cint 5'
'd air.5'

'open E:/lcx/frq.locate.ctl'

'set lon 73 136';'set lat 3 54'
'set cmin 100000';'d loc.6'
'southsea_value loc'

'set lon 73 136';'set lat 3 54'
'set cmin 100000';'d loc.6'
id = 1
while(id<=7000)
    lonlat=read('E:/lcx/wenzhou.txt')
    lonlat=sublin(lonlat,2)
    ln=subwrd(lonlat,1);lt=subwrd(lonlat,2)

    'q w2xy 'ln' 'lt''
    x0=subwrd(result,3);y0=subwrd(result,6)
    'set line 2'
    'draw mark 3 'x0' 'y0' 0.03'
    id=id+1
endwhile


'printim E:/lcx/wenzhou.png white'
;


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

新浪微博达人勋

 楼主| 发表于 2018-1-11 19:46:29 | 显示全部楼层
river 发表于 2018-1-11 18:14
'reinit'
'sdfopen e:/lcx/ncdata/air.2011.nc'
'sdfopen e:/lcx/ncdata/air.2012.nc'

谢谢大佬,很厉害,但我进行批量处理的时候结果只有第一条的坐标点点上去了,后面的一直搞不上去,您能帮忙看一下吗?
'reinit'
'reinit'
'sdfopen e:/lcx/ncdata/air.2011.nc'
'sdfopen e:/lcx/ncdata/air.2012.nc'
'sdfopen e:/lcx/ncdata/air.2013.nc'
'sdfopen e:/lcx/ncdata/air.2014.nc'
'sdfopen e:/lcx/ncdata/air.2015.nc'

'enable print e:/lcx/ncdata/wenzhou/hh.gmf'

'set grads off'
'set parea 0.8 7.5 5.9 10.9'
'set vpage 0 8.5 0 11'
'set mproj scaled';
'set grads off';'set grid off'

alt=471.8
yy=2015
mm=nov
dd=4

'set lon 73 136';'set lat 3 54'
'set lev 'alt''
'set time 06z'dd''mm''yy''
'set cint 5'
'd air.5'

'open E:/lcx/frq.locate.ctl'

'set lon 73 136';'set lat 3 54'
'set cmin 100000';'d loc.6'
'southsea_value loc'

'set lon 73 136';'set lat 3 54'
'set cmin 100000';'d loc.6'
id = 1
while(id<=7000)
    lonlat=read('E:/lcx/wenzhou.txt')
    lonlat=sublin(lonlat,2)
    ln=subwrd(lonlat,1);lt=subwrd(lonlat,2)

    'q w2xy 'ln' 'lt''
    x0=subwrd(result,3);y0=subwrd(result,6)
    'set line 2'
    'draw mark 3 'x0' 'y0' 0.03'
    id=id+1
endwhile
'print'
'c'

'set grads off'
'set parea 0.8 7.5 5.9 10.9'
'set vpage 0 8.5 0 11'
'set mproj scaled';
'set grads off';'set grid off'


alt=511.9
yy=2015
mm=sep
dd=29

'set lon 73 136';'set lat 3 54'
'set lev 'alt''
'set time 12z'dd''mm''yy''
'set cint 5'
'd air.5'

'open E:/lcx/frq.locate.ctl'

'set lon 73 136';'set lat 3 54'
'set cmin 100000';'d loc.8'
'southsea_value loc'

'set lon 73 136';'set lat 3 54'
'set cmin 100000';'d loc.8'
id = 1
while(id<=7000)
    lonlat=read('E:/lcx/wenzhou.txt')
    lonlat=sublin(lonlat,2)
    ln=subwrd(lonlat,1);lt=subwrd(lonlat,2)

    'q w2xy 'ln' 'lt''
    x0=subwrd(result,3);y0=subwrd(result,6)
    'set line 2'
    'draw mark 3 'x0' 'y0' 0.03'
    id=id+1
endwhile
'print'
'c'




'reinit'
'reinit'

这是两个,结果只有第一个可以显示坐标,第二个就不行了,因为我要处理的比较多,所以一个一个的就太慢了,但是不知道问题出现在哪儿~
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2018-1-11 22:02:32 | 显示全部楼层
小霸王和小霸王 发表于 2018-1-11 19:46
谢谢大佬,很厉害,但我进行批量处理的时候结果只有第一条的坐标点点上去了,后面的一直搞不上去,您能帮 ...

没有报错信息?你出来的图呢?
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2018-1-12 17:55:21 | 显示全部楼层
river 发表于 2018-1-11 22:02
没有报错信息?你出来的图呢?

这个好像没办法附图。显示的是query error:syntax is query w2xy lon lat x y size
  hard output file is closed marktype x y size
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2018-1-12 17:55:25 | 显示全部楼层
river 发表于 2018-1-11 22:02
没有报错信息?你出来的图呢?

这个好像没办法附图。显示的是query error:syntax is query w2xy lon lat x y size
  hard output file is closed marktype x y size
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

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

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

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