- 积分
- 141
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2012-3-1
- 最后登录
- 1970-1-1
|
发表于 2013-9-20 19:35:23
|
显示全部楼层
方游 发表于 2012-4-16 13:06
很给力,多谢楼主分享!!
我使用了楼主说明的“'set mpdset cnworld'
'set mpdset cnscs'”后的图是这样子的,底图好怪哦,gs文件如下:
'reinit'
'set parea 1.0 5.0 1.2 5.0'
'open E:/2A12/2A12_data/CenterLonLat/clwmrprs.ctl'
'enable print E:/2A12/2A12_data/CenterLonLat/CenterLonLat_002.gmf'
'set mpdset cnworld'
'set mpdset cnscs'
'set grads off'
'set lon 100 180'
'set lat 0 40'
*'draw map'
'set cmax -100'
'd clwmrprs'
f1='E:/2A12/2A12_data/CenterLonLat/CenterLon.txt'
f2='E:/2A12/2A12_data/CenterLonLat/CenterLat.txt'
i=1
while(i<=1776)
aa=read(f1)
bb=read(f2)
aa1=sublin(aa,2)
bb1=sublin(bb,2)
lon=subwrd(aa1,1)
lat=subwrd(bb1,1)
'q w2xy 'lon' 'lat''
x1=subwrd(result,3)
y1=subwrd(result,6)
'set line 4'
'draw mark 3 'x1' 'y1' 0.03'
i=i+1
endwhile
'print'
'disable print'
'reinit'
结果如下:
底图的陆地和位置都非常有问题,请问该怎么样修改过来呢?
|
|