在grdas图的指定经纬度编注实心点

筱筱筱筱筱雅痴

这是脚本：
'open J:\lw-data\grads\jiaoben\grid.ctl'
'open J:\lw-data\grads\jiaoben\rh.ctl'
'open J:\lw-data\grads\jiaoben\wind.ctl'
i=1
while(i<=560)
'set grads off'
'set grid off'
'set lon 112 120'
'set lat 36 43'
'set t 'i''
'set mpdset hires'
'set map 1 1 9'
'set xlopts 1 6 0.2'
'set ylopts 1 6 0.2'

'set gxout shaded'
'cnbasemap oacres(g,rh.2)'
'cbarn.gs'
'set gxout barb'
'set ccolor 1'
'set cthick 8'
'set digsiz 0.1'
'd -sin(winddir.3*3.14159/180)*wind.3*2.5;-cos(winddir.3*3.14159/180)*wind.3*2.5'
*'set gxout value'
*'set stid off'
*'set ccolor 1'
*'set digsiz 0.1'
*'set dignum 2'
*'d rh.3'
'q time'
m=subwrd(result,3)
'draw title 'm''
'set vpage 0 11 0 8.5'
*'set parea 1 10 1 8'
'q w2xy 117.12 39.04'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3 'x' 'y' 0.2'
'printim J:\lw-data\grads\gs\deal\micaps_rh_'m'.png white'
'c'
i=i+1
endwhile
'reinit'
;

筱筱筱筱筱雅痴

之前画图是可以标注的，用了cnbasemap指令之后就标注不了了，但是不知道怎么改，求解

river

你这个把cnbasemap放在最前面使用了之后，然后又来一个display，达到你只要中国境内数据的目的了吗？感觉这个顺序不对啊。你可以在这个顺序上试试

筱筱筱筱筱雅痴

river 发表于 2018-12-12 22:04
你这个把cnbasemap放在最前面使用了之后，然后又来一个display，达到你只要中国境内数据的目的了吗？感觉这 ...

你的意思是让我先display在cnbasemap吗

我如果不加这一段是可以运行出图的
'q w2xy 117.12 39.04'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3 'x' 'y' 0.2'
可是标注站点的这段指令加上去，就报错，标注点显示不出来

river

筱筱筱筱筱雅痴 发表于 2018-12-13 19:03
你的意思是让我先display在cnbasemap吗

我如果不加这一段是可以运行出图的

你仔细看下cnbasemap这个文件，本身里面自带了display命令，目的是屏蔽中国国境外的数据，所以其他的都先画，最后再用cnbasemap画等值线。试试行不行