GrADS绘制站点分布图-升级版

zhaolanxiai

mofangbao 发表于 2015-9-25 08:25
这个要自己检查啦

谢谢了，我检查了，当在while(1)前面的i=2改成i=3就显示另外一部分图。。。。

∴☆°★°∴

好！顶一个！！！

侠女黒丫

'reinit'
'open d:\xiaoyu\sichuan\xiaoyubiao.ctl'
'set grid off'
'set grads off'
'set map 1 1 5'
'set mpdset sic'
'set lon 97 110.5'
'set lat 25.8 34.5'
'draw map'
'set gxout print'
'd b'
res = result
siz.1=0.1
siz.2=0.15
siz.3=0.2
siz.4=0.25
i=2
while(i)
    1n=sublin(res,i)
    if(1n='');break;endif
    i=i+1
    lat1=subwrd(1n,3)
    lon1=subwrd(1n,2)
    data=sublin(res,i)
    i=i+1
'q w2xy 'lon1' 'lat1''
x1=subwrd(result,3)
y1=subwrd(result,6)
data=subwrd(data,1)*1
'set line 5'
if(data<-1.5)
    'draw mark 2 'x1' 'y1' 'siz.4''
endif

if(data>-1 & data<-1.5)
    'draw mark 2 'x1' 'y1' 'siz.3''
endif
   
if(data>-1 & data<-0.5)
    'draw mark 2 'x1' 'y1' 'siz.2''
endif

if(data>-0.5 & data<0)
    'draw mark 2 'x1' 'y1' 'siz.1''
endif

if(data>0 & data<0.5)
    'draw mark 3 'x1' 'y1' 'siz.1''
endif

if(data>0 & data<0.5)
    'draw mark 3 'x1' 'y1' 'siz.1''
endif

if(data>0.5 & data<1)
    'draw mark 3 'x1' 'y1' 'siz.2''
endif

if(data>1 & data<1.5)
    'draw mark 3 'x1' 'y1' 'siz.3''
endif

if(data>1.5)
    'draw mark 3 'x1' 'y1' 'siz.4''
endif
endwhile

'enable print d:\xiaoyu\sichuan\gs\xyb.gif'
'print'
'printim d:\xiaoyu\sichuan\gs\xyb.gif white'
'disable print'

我按照楼主的方法写了类似gs
可是run之后出错。

错误是这样的

侠女黒丫

侠女黒丫 发表于 2015-12-18 15:29
'reinit'
'open d:\xiaoyu\sichuan\xiaoyubiao.ctl'
'set grid off'

请各位大神帮忙看一下，错误在哪里

gaowl

学习了，这个方法非常有意思。谢谢！

蓝冰草

实用，谢谢楼主～～～

要吃早饭

{:eb502:}{:eb502:}

aas314

真厉害  谢谢分享

微微一笑

楼主呀！怎样能使右下角的南海诸岛上也有大图的站点呀？参考家园里的帖子，可以画出填色图，但是这种站点图怎么弄呢？没有琢磨出来…

hugeal10000

这个不错啊!