- 积分
- 1485
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2014-3-14
- 最后登录
- 1970-1-1
|
发表于 2014-6-15 20:05:09
|
显示全部楼层
谢啦!!我终于发现我的问题出在哪儿了,我没把站点的经纬度数值计算到度,而是小数部分直接用分了,所以出错了,但现在使用了basemap命令,仍无法把超出边界的部分去掉。可以帮我看看我的gs文件哪儿出错了吗?'reinit'
'open d:\diqingjiangshui\10\diqing.ctl'
'open d:\diqingjiangshui\4\grid.ctl'
'open d:\diqingjiangshui\5\eof1.ctl'
'enable print d:\diqingjiangshui\11\eof1.gmf'
'set parea 1 10 1.5 7.75'
'set grads off'
'set grid off'
'set frame on'
'set lon 98 101'
'set lat 26.5 29.5'
'set ylopts 1 5 0.22'
'set xlopts 1 5 0.13'
'define a=oacres(g.2,eof1.3,30)'
'define a1=maskout(a,g.2)'
'define aa=smth9(a1)'
'd aa'
'set gxout contour'
'set clab forced'
'set cint 0.05'
'basemap aa'
'cbarn 1.5 0'
'draw shp d:/diqingjiangshui/10/diqing2.shp'
'print'
'disable print'
; |
|