grads 经纬度坐标转换输出的问题

占戈示申

通过学习台风路径的gs,目前提取出了几个时次的低压中心的经纬度，每次只能用 q w2xy 语句去转换。一个一个输入真的麻烦，有没有人知道如何批量的转换呢？或者有人知道如何把w2xy语句下面输出的x=**,y=**，如何输出一个到txt中吗？

下面的是自己写的一个gs ，但是每次运行都是错误，精通的大虾，可以帮忙修改修改吗？
'reinit'
'open G:\07\2005.ctl'
'set lat 40 60'
'set lon 90 140'
'set lev 500'
'set xlopts 1 4 0.15'
'set ylopts 1 4 0.15'
'set mpdset cnworld'
'set grads off'
'set grid off'
'set cmax 0'
'd hgtprs'
aa='G:\07\xy.txt'
a1=read(aa)
bi=sublin(a1,1)
i=1
while(i<=2)
la=subwrd(b1,i)
lo=subwrd(b1,i+1)
'q w2xy 'lo' 'la''
x1=subwrd(result,3)
y1=subwrd(result,6)
i=i+1
'grads2ascii x1 G:\grads\gs\1x1\ing\x1.txt a %.2f'
'grads2ascii y1 G:\grads\gs\1x1\ing\y1.txt a %.2f'
endwhile
ff=close(aa)
;

占戈示申

{:5_275:}{:5_275:}

侠女黒丫

我也正在纠结这个问题。
之前画27个云南站点图就是一个一个转换的。
现在又是28个四川站点图，要疯了~

占戈示申

侠女黒丫 发表于 2015-12-17 15:32
我也正在纠结这个问题。
之前画27个云南站点图就是一个一个转换的。
现在又是28个四川站点图，要疯了~

这个问题我解决了，通过grads读取txt文件，然后循环，用sublin和subwrd语句，可以读取转换出来的

侠女黒丫

占戈示申 发表于 2015-12-17 20:46
这个问题我解决了，通过grads读取txt文件，然后循环，用sublin和subwrd语句，可以读取转换出来的

可以告知一下解决办法吗，谢谢了~

占戈示申

侠女黒丫 发表于 2015-12-18 14:20
可以告知一下解决办法吗，谢谢了~

刚找到这个gs，希望还是可以帮到你。前一阵子弄的，应该可以用。有什么不对的可以联系我。

侠女黒丫

占戈示申 发表于 2015-12-19 21:21
刚找到这个gs，希望还是可以帮到你。前一阵子弄的，应该可以用。有什么不对的可以联系我。

非常感谢！！

侠女黒丫

占戈示申 发表于 2015-12-19 21:21
刚找到这个gs，希望还是可以帮到你。前一阵子弄的，应该可以用。有什么不对的可以联系我。

呃，，，运行出错了。。
我已经按照你的方式分别放好经纬度、

出错提示是这样的：
error in sublin: 2nd argumenr invalid
    error occurred on line 22
    in file d:\xiaoyu\sichuan\jwdzh.gs

好像是第二参数无效？！
我的gs是这样的：
'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'
**以上只是为了可以单独显示地图

i=2
while(i<=28)
aa='d:\xiaoyu\sichuan\lon.txt'
bb='d:\xiaoyu\sichuan\lat.txt'
****当初是经纬度分别存放在两个不同的txt文件里面了
ln1=read(aa)
lt1=read(bb)
**分别打开两个文件
ln=sublin(ln1,'i')
lo=subwrd(ln,1)
lt=sublin(lt1,'i')
la=subwrd(lt,1)
**一次读取一行，然后循环。

'q w2xy 'lo' 'la''
**转换经纬度
x1=subwrd(result,3)
y1=subwrd(result,6)
say x1
say y1
i=i+1
endwhile
ff=close(aa)
gg=close(bb)
;

其中第22行是ln=sublin(ln1,'i')

并没有改太多，不知为何出错
能否请同学帮忙看一下

侠女黒丫

侠女黒丫 发表于 2015-12-25 20:23
呃，，，运行出错了。。
我已经按照你的方式分别放好经纬度、

还有一句错误信息是：
no hardcopy metafile open

占戈示申

侠女黒丫 发表于 2015-12-25 20:33
还有一句错误信息是：
no hardcopy metafile open

'reinit'
i=1
while(i<=5)
a='C:\OpenGrADS\ts\123.txt'
b=read(a)
c=sublin(b,2)
say c
i=i+1
endwhile
;
你再试验一下这个呢？那个是我挺长时间一场写的，可能有些错误吧，这个我刚才试了一下可以读取txt内的文本的。