- 积分
- 304
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2012-12-27
- 最后登录
- 1970-1-1
|
发表于 2017-5-24 15:19:14
|
显示全部楼层
很早之前的问题,当时哪里出问题也忘记了。
刚才找到gs文件:
'reinit'
'open d:\lunwen\qushi\rqs.ctl'
'set grid off'
'set grads off'
'set map 1 1 5'
'set mpdset cnworld'
'set lon 85 110.5'
'set lat 20 35'
'd r'
**以上只是为了可以单独显示地图
i=1
while(i<=90)
aa='d:\xiaoyu\xinan\lon1.txt'
bb='d:\xiaoyu\xinan\lat1.txt'
****当初是经纬度分别存放在两个不同的txt文件里面了
ln1=read(aa)
lt1=read(bb)
**分别打开两个文件
ln=sublin(ln1,2)
lo=subwrd(ln,1)
lt=sublin(lt1,2)
la=subwrd(lt,1)
'q w2xy 'lo' 'la''
**转换经纬度
x1=subwrd(result,3)
y1=subwrd(result,6)
say x1
say y1
file1="d:\xiaoyu\xinan\zbx1.txt"
file2="d:\xiaoyu\xinan\zby1.txt"
rc=write(file1,x1,append)
rc=write(file1,a,append)
rc=write(file2,y1,append)
rc=write(file2,b,append)
say m
i=i+1
endwhile
ff=close(aa)
gg=close(bb)
;
|
|