- 积分
- 1527
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2012-8-21
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
其实就是一个台风中心位置标注的事,以前曾写个几个都成功了,这次是放在一个循环里面,老是出不来结果,提示经纬度变化直角坐标有问题,捯饬了一晚上,实在是找不出哪里有错,故请教论坛里的大侠们,赐教。先把gs及报错信息放上:
*—————————-------------
********draw the symbol of typhoon*******
lp=1
while(lp<=35)
* '文件共有35行,即35个经纬度点
parap=read('F:/besttrack/lj/1973_0303.txt')
parap=sublin(parap,2)
*get day time latitude lontitude**
ymd1=subwrd(parap,1)
day1=substr(ymd1,7,2)
time1=substr(ymd1,9,2)
if(d1=day1&t2=time1)
lat1=subwrd(parap,3)
lon1=subwrd(parap,4)
say 'lon1 ='lon1
say 'lat1 ='lat1
*注这一步已经出来了,显示lon1=125.6,lat1=19.7
*转成直角坐标
'q w2xy 'lon1' 'lat1''
x1=subwrd(result,3)
y1=subwrd(result,6)
say 'laty=' y1
say 'lonx='x1
注:这一步就不行了,x1,y1出不来
'draw mark 8 'x1' 'y1' 0.12'
*'draw wxsym 3 'x1' 'y1' 0.08 2 5'
endif
lp=lp+1
endwhile
*----------------------
实在想不通,为什么出不来?求助
|
|