台风路径

酒国

之前在气象论坛里找到了画台风路径的方法，比较不错，这是其链接。http://bbs.06climate.com/forum.php?mod=viewthread&tid=18391&highlight=%CC%A8%B7%E7%C2%B7%BE%B6
在画的过程中，有个问题不明白，在楼主原帖留言没有答复，所以想请大家帮忙解决解决。

txt数据跟gs文件怎么联系呀？随便找一个ctl是什么意思啊？

无敌皓皓

grads通过a=read("XXX.txt")读取txt数据。
随便找一个ctl是为了画一个变量来召唤地图，否则没有地图就没办法画路径。但为了只显示地图不显示用来召唤地图随便画的那个变量，于是设置cmin特别大。

酒国

无敌皓皓 发表于 2015-3-9 22:07
grads通过a=read("XXX.txt")读取txt数据。
随便找一个ctl是为了画一个变量来召唤地图，否则没有地图就没办 ...

非常感谢你。。。还有个问题，他的gs文件里并没有a=read("XXX.txt"),要放在哪里呢？
麻烦你了。。。。。非常感谢
**使用方法：在ga>路径\path.gs 台风文件路径 输出台风图片路径
**把初始点的经纬度写改写在lat0,lon0处
function path(args)
'reinit'
infile=subwrd(args,1)
outpng=subwrd(args,2)
say infile
say outpng
'enable print 'outpng
'open e:\20110702\070218.ctl'
'set parea 1 11 1 8'
'set grads off'
'set lat 5 55'
'set lon 90 160'
'set cmin 10000'
'set xlopts 1 0.1 0.16'
'set ylopts 1 0.1 0.16'
'd TMP2m'
i=1
strong1=1
c1=1
c2=1
c3=1
c4=1
c5=1
c6=1
lat0=25.7
lon0=147.4
'q w2xy 'lon0' 'lat0
x0=subwrd(result,3)
y0=subwrd(result,6)
while(1)
info=read(infile)
code=sublin(info,1)
if(code!=0)
break
endif
info=sublin(info,2)
latty=subwrd(info,3)*0.1
lonty=subwrd(info,4)*0.1
tim=subwrd(info,1)
wind=subwrd(info,6)
strong0=subwrd(info,2)
'q w2xy 'lonty' 'latty
x=subwrd(result,3)
y=subwrd(result,6)
'set string 2 l 2 0'
'set strsiz 0.15 0.15'
'set line 1 1 5'
'draw line 'x0' 'y0' 'x' 'y
x0=x
y0=y
st=substr(tim,5,10)
if(wind>10&wind<=17.1&c1=1)
'draw wxsym 41 'x' 'y' 0.2 4 1'
c1=c1+1
endif
if(wind>17.1&wind<=24.4&c2=1)
'draw wxsym 41 'x' 'y' 0.2 9 1'
c2=c2+1
endif
if(wind>24.4&wind<=32.6&c3=1)
'draw wxsym 41 'x' 'y' 0.2 7 1'
c3=c3+1
endif
if(wind>32.6&wind<=41.4&c4=1)
'draw wxsym 41 'x' 'y' 0.2 8 1'
c4=c4+1
endif
if(wind>41.4&wind<=50.9&c5=1)
'draw wxsym 41 'x' 'y' 0.2 6 1'
c5=c5+1
endif
if(wind>50.9&c6=1)
'draw wxsym 41 'x' 'y' 0.2 2 1'
c6=c6+1
endif
*if(math_fmod(i,4)=0)
if(strong0!=strong1|i=1)
strong1=strong0
if(strong0=1)
'set string 1 l 1 1'
'draw string 'x' 'y' 'st'--TD'
'draw wxsym 41 'x' 'y' 0.2 4 1'
endif
if(strong0=2)
'set string 1 l 1 1'
'draw string 'x' 'y' 'st'--TS'
'draw wxsym 41 'x' 'y' 0.2 9 1'
endif
if(strong0=3)
'set string 1 l 1 1'
'draw string 'x' 'y' 'st'--STS'
'draw wxsym 41 'x' 'y' 0.2 7 1'
endif
if(strong0=4)
'set string 1 l 1 1'
'draw string 'x' 'y' 'st'--TY'
'draw wxsym 41 'x' 'y' 0.2 8 1'
endif
if(strong0=5)
'set string 1 l 1 1'
'draw string 'x' 'y' 'st'--STY'
'draw wxsym 41 'x' 'y' 0.2 6 1'
endif
if(strong0=6)
'set string 1 l 1 1'
'draw string 'x' 'y' 'st'--SuperTY'
'draw wxsym 41 'x' 'y' 0.2 2 1'
endif
endif
i=i+1
endwhile
'print'
'disable print'
;

无敌皓皓

info=read(infile)应该就是读取txt的，infile在开始的args里面：infile=subwrd(args,1)
我没有用过这个脚本不好意思。不过原理就是这样啦你可以自己修改啦～

酒国

无敌皓皓 发表于 2015-3-10 20:52
info=read(infile)应该就是读取txt的，infile在开始的args里面：infile=subwrd(args,1)
我没有用过这个脚 ...

再次感谢。现在可以出图了。。。但是怎么会这样。。多了两条长线。。。

酒国

酒国 发表于 2015-3-11 09:20
再次感谢。现在可以出图了。。。但是怎么会这样。。多了两条长线。。。

txt里面的数据没有超出画图的范围、、、、、、

nicocle

怎么把TXT跟GS联系起来？

IAMQWW1

楼主呢，能不能解释一下？

酒国

IAMQWW1 发表于 2015-9-6 15:32
楼主呢，能不能解释一下？

程序就要意思是先调去去地图，当然了，随便一个ctl一般都符合，然后把台风路径的txt通过函数来调取出来，对里面的内容进行判断，在那个地方，是什么强度

mm_cat

酒国 发表于 2015-3-11 09:20
再次感谢。现在可以出图了。。。但是怎么会这样。。多了两条长线。。。

你好，我也遇到了多出线的问题，请问你最后解决了吗？是怎么弄的？
这里删不了... 已经解决了，不麻烦了，谢谢！