用grads画散点图超出地图边界

姗姗来迟

用grads画散点图时，散点超过地图边界，我觉得可能是对set parea的理解不清楚，但是具体不知道怎么改？
此外，想问一下，画散点图时必须要先画其它图吗？不画的话地图就显示不出来，请各位指导一下，非常感谢。

lqouc

这个不是散点，是站点填色，看起来是投影的问题。
提问请给出相关代码。

姗姗来迟

lqouc 发表于 2019-1-11 10:56
这个不是散点，是站点填色，看起来是投影的问题。
提问请给出相关代码。

非常感谢您的回复，下面是我的gs代码
'c'
'reinit'
'open E:\sjj\39mon01.ctl'
'enable print E:\example\scatter\yu99scatter.gmf'
lon.1=70
lon.2=140
lat.1=15
lat.2=55
'set lon 'lon.1' 'lon.2''
'set lat 'lat.1' 'lat.2''
'set parea 2 7 3 8'
'set xlopts 1 5 0.12'
'set ylopts 1 5 0.12'
'set clopts 1 5 0.08'
'set mpdset cnworld'
'set map 15 1 4'
'run E:\example\scatter\sevencolor.gs'
'set grid off'
'set grads off'
'set ylint 10'   
'set xlint 20'   
'set gxout shaded'
'set cmin 40000'
'd z'
pull down
'close 1'
'open E:\example\scatter\st693y99.ctl'
'set parea 2 7 3.85 7.2'
'set gxout scatter'
'set vrange 'lon.1' 'lon.2' '
'set vrange2 'lat.1' 'lat.2''
*'set digsize 0.13'
'set xlab off'
'set ylab off'
'set frame off'
***************
'set t 1 74'
'set ccolor 21'
'set cmark 3'
'set digsize 0.1'
'd  jd;wd'

ctl文件代码如下：
dset E:\example\scatter\st693y99.grd
TITLE position of the yuzhi99%
undef 32766
XDEF 1 levels 1
YDEF 1 levels 1
ZDEF 1 levels 0 1
TDEF 693 linear 1Jan2011 1dy
VARS 3
jd 0 99  
wd 0 99  
yu 0 99  
ENDVARS

lqouc

你还真是用散点画的啊。。。。。

lqouc

那需要改投影，set mpdset lonlat 我记得大概是这个

姗姗来迟

lqouc 发表于 2019-1-12 10:06
那需要改投影，set mpdset lonlat 我记得大概是这个

您好，您说的是set mproj latlon这个命令吗？ 查过说其表示用固定的投影角进行lat/lon投影，不理解它的具体用法，另外，代码中问题在哪里，还希望您给指出

lqouc

先把两个set parea 改的一致了。如果不行再在gs里面加上set mproj scaled。

姗姗来迟

lqouc 发表于 2019-1-12 14:31
先把两个set parea 改的一致了。如果不行再在gs里面加上set mproj scaled。

非常感谢，加了set mproj scaled果然就好了，我自己弄了一上午都不行，您一下就把问题解决了。

lqouc

姗姗来迟 发表于 2019-1-12 15:05
非常感谢，加了set mproj scaled果然就好了，我自己弄了一上午都不行，您一下就把问题解决了。

不用客气，欢迎常来论坛

yuyuxuan

lqouc 发表于 2019-1-12 10:05
你还真是用散点画的啊。。。。。

请问一下  在画scatter图的时候，所设置的vars1和vars2是什么呀，谢谢了