- 积分
- 15959
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2011-8-10
- 最后登录
- 1970-1-1
|
发表于 2015-9-10 09:45:21
|
显示全部楼层
你好,按照你basemap贴子里的方法,我做了一下我们这里的站点降水量,其它的都还好,就是两个问题,第一,我的图画出来之后,要么我要的里面是白的,要么里面是黑的。第二,我最后用printim存图的时候,相应的目录下面没有图出来。
我把我的GS发给你,麻烦你帮我看看。谢谢!!
'reinit'
'open D:\0\1\grid.ctl'
'open D:\0\1\sta.ctl'
'set grads off'
'set grid off'
'set lon 111.3 114.3'
'set lat 36.6 38.2'
'define rgrid=oacres(gst,rain.2,1.3,1.3)'
'set mpdset hires'
'set csmooth on'
'set gxout shaded'
'set rgb 26 166 242 143'
'set rgb 27 61 186 61'
'set rgb 28 97 184 255'
'set rgb 29 0 0 255'
'set rgb 30 250 0 250'
'd rgrid'
'set clevs 10 25 50 100'
'set ccols 26 27 28 29 30'
'province-basemap map rgrid 1'
'set mpdset map'
'draw map'
'cbarn 1 1 0.7 2.5'
cnt=23
j=1
while(j<cnt)
ret=read('D:\0\1\dksts1.txt')
aa1=sublin(ret,2)
name=subwrd(aa1,1)
lat=subwrd(aa1,2)
lon=subwrd(aa1,3)
'q w2xy 'lon' 'lat''
x = subwrd(result,3)
y = subwrd(result,6)
*写汉字
'writehz 'x+0.05' 'y-0.05' 'name' 1 2'
*画点
'draw mark 3 'x' 'y' 0.05'
j=j+1
endwhile
ret=close('D:\0\1\dksts1.txt')
'printim D:\0\1\00.png white'
*'reinit'
; |
|